- Key Laboratory of Molecular and Nano Probes, College of Chemistry, Chemical Engineering and Materials Science, Collaborative Innovation Center of Functionalized Probes for Chemical Imaging, Ministry of Education, Shandong Normal University, Jinan, China
Noble-gas bonding represents curiosity. Some xenon hydrides, such as HXeY (Y = Cl, Br, I) and their hydrogen-bonded complexes HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CN, CCH), have been identified in matrixes by observing H-Xe frequencies or its monomer-to-complex blue shifts. However, the H-Xe bonding in HXeY is not yet completely understood. Previous theoretical studies provide two answers. The first one holds that it is a classical covalent bond, based on a single ionic structure H-Xe+ Y−. The second one holds that it is resonance bonding between H-Xe+ Y− and H− Xe+-Y. This study investigates the H-Xe bonding, via unusual blue-shifted phenomena, combined with some NBO/NRT calculations for chosen hydrogen-bonded complexes HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CN, CCH). This study provides new insights into the H-Xe bonding in HXeY. The H-Xe bond in HXeY is not a classical covalent bond. It is a charge-shift (CS) bond, a new class of electron-pair bonds, which is proposed by Shaik and Hiberty et al. The unusual blue shift in studied hydrogen-bonded complexes is its H-Xe CS bonding character in IR spectroscopy. It is expected that these studies on the H-Xe bonding and its IR spectroscopic property might assist the chemical community in accepting this new-class electron-pair bond concept.
Introduction
The chemical bond is the most central concept in chemistry. The model of chemical bonding can help chemists to understand and design matter. Despite the apparent utility of the model of chemical bonding, it is incomplete. Developing bonding models is undoubtedly important to our understanding of novel molecules, such as the molecules of noble gas chemistry (Grandinetti, 2018).
The challenge to noble-gas bonding comes mainly from the inertness of noble-gas atoms. However, significant progress has been made in noble-gas chemistry during the past 20 years. Besides noble-gas hydrides (Khriachtchev et al., 2009), noble gas–noble metal complexes have been studied for their thermodynamic and kinetic stabilities in theory (Jana et al., 2017, 2018a,b; Pan et al., 2019; Saha et al., 2019). Experimentally, about 30 noble-gas hydrides had been identified by the beginning of 1995 (Pettersson et al., 1995; Räsänen et al., 2000; Lundell et al., 2002; Feldman et al., 2003; Khriachtchev et al., 2003, 2009; Duarte and Khriachtchev, 2017). About 10 hydrogen-bonded complexes between HNgY (Ng = Xe, Kr) and H2O/HCl/HBr/HI/HCCH have been well-characterized via infrared (IR) spectroscopy. Interestingly, these complexes show unusual shifts of the H-Xe stretching vibration. For example, the H-Xe stretching mode of the complex HXeI···HCCH exhibits a blue shift of 49 cm−1, in comparison with the monomer HXeI (Zhu et al., 2015). Some other complexes such as HXeBr···H2O and HXeI···H2O are characterized by much larger experimental blue shifts of the H-Xe stretching frequency (>100 cm−1) (Tsuge et al., 2014). The largest blue shift (300 cm−1) has been reported for the H-Kr stretching mode of the complex HKrCl···HCl (Corani et al., 2009). These experimental findings provide theoretical chemistry researchers with excellent opportunities to develop a bonding model for noble-gas hydrides.
In pioneering theoretical work, Last and George put forward one simple ionic structure model H-Ng+Y− (Last and George, 1988), where H-Ng+ belongs to a classical covalent bond, to be exact, an electron-sharing bond, while the interaction between H-Ng+ and Y− comes from electrostatic attraction. Energy decomposition analyses (EDAs) carried out by Frenking for HArF provided a consistent bonding picture (Lein et al., 2004). Besides this ionic structure model, empirical resonance bonding models were proposed by Räsänen's group (Pettersson et al., 1999) and Alabugin's group (Alabugin et al., 2004). Recently, our group has carried out the resonance bonding analyses for noble-gas hydrides (Zhang et al., 2016) by using Natural Bond Orbital (NBO) and Natural Resonance Theory (NRT) methods (Glendening and Weinhold, 1998a,b; Glendening et al., 1998, 2001, 2013a,b; Weinhold and Klein, 2014; Weinhold et al., 2016). We found that each molecule HNgY could be best described as three structures, H-Xe+Y−, H−Xe+-Y, and H∧Y, where the first two resonance structures mix to form hyperbonding of H-Xe+Y−↔ H−Xe+-Y, as proposed by Weinhold and Landis (2005b). Such bonding provides a picture of resonance covalency for noble-gas hydrides (Weinhold and Klein, 2012, 2014; Landis and Weinhold, 2013; Zhang et al., 2015, 2017a; Jiao and Weinhold, 2019).
It is important to note that they are two variants of covalent bonding, which are related to two kinds of electron-pair bonds. They are our familiar electron-sharing bonds and dative bonds. The difference between these two types of bonds is due to the origin of the bonding electrons. In electron-sharing bonds, each fragment provides one electron. In dative bonds, both electrons come from one fragment, which donates two electrons to the vacant orbital of the other. Also note that besides our familiar electron-sharing bonds and dative bonds there is one new kind of electron-pair bonds, charge-shift (CS) bonds. The concept of charge-shift bonds was first proposed by Shaik and Hiberty et al. in 1992, to describe a new-class electron-pair bonds, such as F2 (Shaik et al., 1992). Ten years later, they further presented experimental manifestations of CS bonding (Shaik et al., 2005, 2009). They have applied the concept of CS bonds to the understanding of bonding and stabilities of several hypervalent molecules (Braïda and Hiberty, 2013; Braïda et al., 2014), such as XeF2. In 2018, Grandinetti pointed out that a stabilization like HNgY is, generally, known as CS bonding (Grandinetti, 2018). Indeed, the decades following their original work saw more CS bonding molecules. Very recently, they published their latest review on CS bonds (Shaik et al., 2020).
Given that XeF2 and HXeY are similar in geometrical and electronic structures, two obvious questions arise: (1) Is the H-Xe bond in HXeY a charge-shift bond? (2) What is the H-Xe CS bonding character in IR spectroscopy? This study explores these two questions through blue-shifted phenomena, with the help of NBO/NRT analyses. We choose HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CN, CCH) as study systems. Some of them are identified in the matrix experiment. We analyze the relationship between the H-Xe bonding and H-Xe blue shifts. We aim to understand the H-Xe bonding in HXeY and to extend CS bonding concepts to noble-gas hydrides.
This paper will be organized as follows. Firstly, we summarize the computational details and discuss the geometrical structural details and H-Xe IR spectroscopic properties for our studied monomers and complexes. Secondly, we analyze the resonance bonding of HXeY, especially for the bonding between H and Xe. Thirdly, we analyze the H-Xe bond order. And we determine that the H-Xe bond order must include the contributions of two resonance structures. Fourthly, we confirm that the H-Xe in HXeY is a charge-shift bond, and further analyze its CS bonding character. Finally, we present the concluding summary, with emphasis on the H-Xe CS bonding and its bonding character.
Computational Details
The geometry optimization and vibrational frequency calculations were carried out with the Gaussian 09 program at the level of the second order MØller-Plesset perturbation theory (MP2) (Head-Gordon et al., 1988; Frisch et al., 2009). The def2-TZVPPD basis set was used. This basis set, taken from the EMSL basis set library (Feller, 1996; Schuchardt et al., 2007), is the triple-zeta-valence basis set augmented with two sets of polarization and diffuse basis functions. No imaginary frequencies were found in any case, which confirms that our obtained structures are true local minima on the potential energy surface. The NBO and NRT were employed to analyze the bonding of our studied systems with the NBOPro 6.0 program (Glendening and Weinhold, 1998a,b; Weinhold, 2012, 2013). Directed NBO analyses could provide the second order perturbation energy of a donor-acceptor interaction in the best Lewis structure. For any other resonance structure, we use $CHOOSE keylist to calculate the second order perturbation energy of a donor-acceptor interaction. It is important to attach the $NRTSTR keylist in NRT analyses, for it insures a consistent set of reference structures for NRT comparisons of studied complexes. The NBO-based natural resonance theory complements and extends NBO analyses to other resonance structures. By using $NRTSTR keylist to specify key structures as reference structures, we can obtain accurate weightings (ωI, ωII, ωIII…) of resonance structures, and NRT bond orders (bAB) that express the strength of resonance-weighted chemical bonds between any atom pair. More importantly, the NBO/NRT-based models provide a framework for analyzing chemical bonding in terms of familiar concepts, such as Lewis structures, resonance, and donor-acceptor interactions. Herein, we use NBO/NRT methods to analyze H-Xe bonding. Besides, the NBOview 2.0 module was acquired to obtain the orbital overlap graphics.
Results and Discussions
Geometry and Blue Shifts
In general, four structures need to be considered for each of the complexes HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CN, CCH). Take HXeY···HCl, as shown in Figure 1, as one representative example. In the first two structures (A and B), the HCl moiety is closed to the halogen atom of the monomer HXeY to form a bent structure. Structure A is stabilized with the Cl-H···Y hydrogen bond, whereas Structure B is dependent on the Xe-Y···Cl halogen bond. Structure C, whose halogen atom in the moiety is collinear with the monomer HXeY, is formed by the Cl atom interacting with the H atom. Structure D is stabilized with the dihydrogen Cl-H···H-Xe bond with all atoms in line. Experimentally, it has been observed that several infrared absorption bands originate from the H-Xe stretching mode for our chosen complexes. With the aid of quantum chemical calculations, they are assigned to Structure A (Lignell et al., 2008; Tsuge et al., 2013, 2014; Zhu et al., 2015). Therefore, the following analyses are restricted to Structure A.
As shown in Structure A in Figure 1, HXeY maintains its original linear structure for all the HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CN, CCH) optimized geometry. At the MP2/def2-TZVPPD level, the calculated H-Xe and Xe-Y bond lengths (RH−Xe, RXe−Y) and vibrational frequency shifts of H-Xe stretching mode in complexes HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CN, CCH) as well as in monomers HXeY are collected in Table S1 and compared with the available experimental data. Note that Räsänen's group has computed RH−Xe and RXe−Y as well as a variety of H-Xe frequency blue shifts for some of the complexes HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CN, CCH) at the CCSD(T) level. The difference in calculated H-Xe bond lengths at both MP2 and CCSD(T) levels is <0.060 Å, whereas for Xe-Y bond lengths the largest difference is 0.048 Å in the monomer HXeI. These comparisons show that the MP2/def2-TZVPPD is an appropriate model chemistry for this study. Thus, the following discussion will be based on the calculated data at the MP2/def2-TZVPPD level.
Data in Table S1 show that the complexation has non-negligible influences on monomers HXeY. On the one hand, the bond lengths of H-Xe in complexes become shorter as compared with the value in the corresponding HXeY monomer. Taking HXeCl species as one example, H-Xe bond lengths RH−Xe in complexes are 1.646, 1.645, 1.644, 1.646, 1.648, and 1.655 Å for HX = H2O/HCl/HBr/HI/HCN/HCCH, respectively, while the H-Xe bond length of the monomer HXeCl is 1.666 Å. This indicates that the interaction between H and Xe atoms in monomers HXeY is weaker than that in corresponding complexes HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CN, CCH). On the contrary, the bond lengths of Xe-Y within complexes are slightly larger than that in the corresponding monomer HXeY, which means that the interaction between Xe and Y atoms is stronger in monomers than in corresponding complexes. It is worthwhile noting that the calculated RH−Xe ranges from 1.644 to 1.708 Å, slightly larger than the covalent limits Rcov[rcov(H)+rcov(Xe)] for 1.63 Å and significantly shorter than VdW limits Rvdw[rvdw(H)+rvdw(Xe)] for 3.86 Å (Pyykkö, 2015; Rahm et al., 2016). This indicates strong covalency of the H-Xe bond.
On the other hand, the complexation leads to experimentally observable blue shifts of the H-Xe stretching frequency, which seems to be normal phenomena according to Khriachtchev's group studies. The calculated and experimental vibrational frequency shifts of the H-Xe stretching mode are also collected in Table S1. Taking hydrogen-bonded complexes HXeCl···H2O, HXeBr···H2O and HXeI···H2O as one example group, calculations at the MP2 level predict blue shifts of 101, 116, and 139 cm−1, respectively. The experimental data are 82, 101, and 138 cm−1, respectively. Obviously, the calculated complexation-induced spectroscopic shift of the H-Xe stretching mode at the MP2 level is in qualitative agreement with the experimental values. These preliminary structural and blue-shifted analyses provide the backdrop for the following exploration of H-Xe bonding for Xe hydrides.
H-Xe Resonance Bonding in HXeY
Earlier studies have found that the NBO/NRT method is a helpful tool to explore the resonance bonding because NBO analyses can provide the best Natural Lewis Structure (NLS), identify donor-acceptor orbital interactions, and estimate the second-order perturbation energy [E(2)] of each donor-acceptor orbital interaction (Glendening et al., 2013a,b; Weinhold, 2013; Weinhold et al., 2016).
Our studies begin with NBO/NRT analyses for HXeY species. The results show that each of the studied HXeY could be better described as a hybrid of the three structures (I,II,and III) as shown in Figure 2, where the H-Xe+:Y− (I) is the best NLS. Table S2 lists three types of donor-acceptor interactions and the value of the second-order perturbation energies [E(2)] of the studied HXeY. For the best NLS H-Xe+:Y− (I), the donor-acceptor interaction () takes place between the lone pair orbital of Y (nY) and the antibonding orbital of H-Xe moiety (). Pay particular attention to such a delocalization interaction. It represents resonance mixing between H-Xe+:Y− (I) and H:− Xe+-Y (II), where the latter corresponds to the lone pair of H atom (nH) delocalizing to the antibonding orbital of Xe-Y moiety (). As proposed by Weinhold et al., these two structures make up resonance bonding. Additionally, there is a non-negligible long-bonding structure H∧Y (III), shown in the final entry in Figure 2. It arises from the delocalization of Xe atom lone pair (nXe) to the antibonding orbital of H-Y moiety (). Special attention is paid to the unusual values of E(2) for and interactions. Such results show that they are no longer suitable for the description of low-order perturbative NLS limit. Even in the unavailable value of some second-order perturbation energies, the importance of these three donor-acceptor interactions can be exhibited by the orbital overlap contour diagrams or 3-D surface views. Figure 3 presents one illustrative example. These results are consistent with our previous studies for HXeY noble-gas hydrides (Y = Cl, Br, I) at the B3LYP level of theory (Zhang et al., 2016).
Figure 3. NBO orbital contour diagrams and 3-D surface views of donor-acceptor interactions for HXeCl.
Once again, the importance of resonance bonding H-Xe+ Y− ↔ H− Xe+-Y is emphasized. On the one hand, resonance bonding is an essential feature of H-Xe bonding in HXeY. On the other hand, we note that earlier studies on the bonding of noble-gas hydrides have already analyzed the leading resonance structure H-Xe+ Y− (Pérez-Peralta et al., 2009; Juarez et al., 2011). In the next section, we will carry out detailed analyses on H− Xe+-Y, as well as H-Xe+ Y−.
H-Xe Bond in HXeY Is Not a Classical Covalent Bond
HXeY is a particularly simple and interesting molecule. Its H-Xe stretching modes provide experimental probes to learn more about the bonding of noble-gas hydrides. It was reported that the H-Xe stretching vibration frequencies of HXeY species are blue shifts upon complexation. In this section, we will explore the H-Xe bonding via this complexation effect, combining with quantitative NRT analyses on chosen hydrogen-bonded complexes. Table 1 displays the weighting of three resonance structures upon complexation for Xe cases. Obviously, the moiety H2O/HCl/HBr/HI/HCN/HCCH complexing with HXeY has a significant influence on the weighting of three resonance structures for HXeY. To be specific, for the long-bonding structure (Weinhold et al., 2016; Zhang et al., 2017b, 2018), its weighting always decreases upon complexation. This indicates that the complexation of the moiety H2O/HCl/HBr/HI/HCN/HCCH studied here is not beneficial to the stability of long-bonding in HXeY. For two other resonance structures, the ωI decreases while ωII increases relative to the corresponding monomer in studied hydrogen-bonded complexes except for the complex HXeCl···HCCH. For example, the weightings of H-Xe+:I− and H:− Xe+-I in the hydrogen-bonded complex HXeI···HI are 62.8 and 19.2%, compared with 66.4 and 9.5% in the monomer HXeI. In contrast, the complex HXeCl···HCCH shows a different trend. As it is easily seen, the weightings of H-Xe+:Cl− and H:− Xe+-Cl structures in the monomer HXeCl are 74.3 and 10.9% respectively; in the complex HXeCl···HCCH they are 76.1 and 9.7%, respectively. In short, for all of our studied complexes, the HXeY complexation with small molecules always leads to a decrease of the weighting about the long-bonding structure. For weightings of these two other resonance structures, one structure always shows a decreasing trend while the other exhibits an increasing trend for all of the studied complexes. In most of the Xe cases studied here, the complexation results in an increasing weighting of the resonance structure H-Xe+:Cl−, with a decreasing weighting of the resonance structure H:− Xe+-Cl. In HXeCl···HCCH the situation is different. An opposite trend is seen for the complex HXeCl···HCCH. From preliminary NBO analyses on our studied complexes, peculiarities in HCCH complexes shown here, in Tables 1–3, and Figure 4, may be due to nXe(d) → donor-acceptor interaction.
Table 1. The weighting and its change of three resonance structures (H-Xe+ :Y−, H:− Xe+-Y, H∧Y) upon HXeY complexation with H2O/HCl/HBr/HI/HCN/HCCH, compared with the monomer HXeY.
Figure 4. Correlation plot for the sum of electron-sharing and dative H-Xe/Kr bond orders [bH−Xe/Kr (electron-sharing) + bH−Xe/Kr (dative)]–H-Xe/Kr bond lengths (RH−Xe/Kr).
Table 2 lists natural bond orders of the H-Xe and H−Y bonds in the monomer HXeY and its complex. It is worthwhile noting that the weighting and the corresponding bond order are equivalent for our studied cases in NBO/NRT framework if the former is expressed as a fraction rather than a percentage. For instance, the weightings of the H-Xe+:Cl− in the monomer HXeCl and in the complex HXeCl···HCl are 74.3% and 71.8%, with corresponding bond order bH−Xe 0.743 and 0.718, respectively. Generally, the larger the bond order is, the stronger the bond is. The data obtained from current NBO programs show the decrease of H-Xe bond orders, reflecting the weakening H-Xe bonds upon complexation. Obviously, these calculated results do not reflect the experimental fact: the strengthened H-Xe bond. This disagreement confirms that the H-Xe bond in HXeY is not a classical covalent bond.
Table 2. The bond orders (b) of H-Xe and H∧Y bonds for the monomer and their hydrogen-bonded complex, with the changes shown in parentheses.
This conclusion deserves some illustration. According to the natural bond orders' definition (Weinhold and Landis, 2005a), the H-Xe bond in HXeY is regarded as a classical covalent bond. The resonance structure H:− Xe+-Y does not contribute to the H-Xe bond order at all. The H-Xe bond order is only contributed to by the resonance structure H-Xe+:Y−. Thus, the H-Xe bond order in Table 2 reflects merely the electron-sharing contribution in the H-Xe bond. It is insufficient to reflect the real strength between H and Xe in HXeY.
H-Xe Bond in HXeY Is a CS Bond
The above studies on the H-Xe bond in HXeY confirm that it is not a classical covalent bond. Then, is it a CS bond? To address this question, we first need to solve the problem of the H-Xe bond order. Let us return to the original NBO/NRT theory. In the framework of NRT theory (Weinhold and Landis, 2005a), the H-Xe total bond strength in HXeY could be written in resonance-averaged form D(H-Xe) = ωIx DI + ωIIx DII, where “ωI” and “ωII” respectively correspond to the weighting of H-Xe+ Y− and H:− Xe+-Y; “DI” and “DII” represent the H-Xe bond strength in H-Xe+ Y− and in H:− Xe+-Y respectively. As we know, a bond order is roughly proportional to the bond strength or the bond length. If introducing two arbitrary constants is presumed to be expressed in terms of the same factor k, DI, and DII can be written in the form, DI= k x bI and DII= k x bII. Importantly, we obtain that b(H-Xe) = ωIx bI + ωIIx bII, where bI represents the H-Xe bond order in resonance structure I, as we know, bI= 1; bII refers to the H-Xe bond order in resonance structure II. Therefore, a calculation that includes two resonance structures is needed to deal with the H-Xe bond order. Note that the procedure employed to calculate the H-Xe CS bond order can be found in Supplementary Material “Explanation of the Procedure Employed to Calculate the BO of the H-Xe Bond.”
To obtain bII, we have to carry out some analyses on H:− Xe+-Y. For this structure, where does the H-Xe bonding originate from? Here we emphasize that the H:− Xe+-Y structure is a natural Lewis structure, in the NBO/NRT language. Consideration of the antibond of the Xe+-Y could lead to extension of the elementary Lewis structure concept to include hyperconjugative delocalization corrections in simple NBO perturbative estimates (Weinhold, 2012). Such hyperconjugative delocalization forms one starting point for our understanding of bonding about the H:− Xe+-Y structure. As shown in Figure 3, one donor-acceptor interaction exists in the H:− Xe+-Y structure. On the basis of this result, we propose that the H-Xe bonding about the structure H:− Xe+-Y is attributed to this donor-acceptor interaction. In our familiar language, it is dative bonding due to hyperconjunctive interaction.
However, the question was still left open: How do we estimate the degree of this H-Xe dative covalency? For this question, our research ideas are from the natural bond orders' definition proposed by Weinhold and Landis (2005b). If such a single dative bond order bII is defined as 1, the NRT bond order of the dative structure is equal to the corresponding fractional weighting. Taking HXeCl as one illustrative example, the dative weightings of monomer HXeCl and hydrogen-bonded complex HXeCl···HCl are 10.9 and 14.6%, respectively. This simple method lets us respectively estimate the dative bond orders: 0.109 and 0.146.
The above discussion is an effort to rationalize bII. We now begin a discussion with a sum of ωIx bI, and ωIIx bII. For convenience, we use bH−Xe (electron-sharing) = ωIx bI,, bH−Xe (dative) = ωIIx bII. As shown in NRT theory, it is a sum of bH−Xe (dative) and bH−Xe (electron-sharing) that can reflect the H-Xe total strength in HXeY. Figure 4 shows correlation plots between the H-Xe bond order [bH−Xe (dative) + bH−Xe (electron-sharing) and the H-Xe bond length RH−Xe for our studied Xe species. One additional example is also shown in Figure 4 for Kr analogs. Good correlation shown in Figure 4 for our studied species except Y = I, provides evidence to support our estimated method. Not enough good correlation in Y = I case may be due to a larger coupling effect between in HXeY and in its H-bonding complex.
On the basis of the data of the H-Xe bond order, we can analyze the H-Xe total bonding of HXeY. From preceding NRT analyses, we have shown that the electron-sharing bond orders [bH−Xe (electron-sharing)] in the monomer and in the complex are 0.743 and 0.718, respectively, and 0.109 and 0.146 for the dative bond orders [bH−Xe (dative)]. Thus, summing bH−Xe (dative) to bH−Xe (electron-sharing) yields 0.852 and 0.864 for the HXeCl monomer and the HXeCl···HCl complex, respectively. Obviously, the HXeCl complexing with the molecule HCl leads to an increase of the H-Xe bond order. The same is true for other complexes studied here. Table 3 lists the total bond order of the H-Xe bond in HXeY and in its hydrogen-bonded complexes. Note that it includes dative contribution and electron-sharing contribution of H-Xe bonding. As shown in Table 3, the HXeY complexing with the small molecule leads to an increase of the H-Xe bond order. In other words, it is an enhancement of the H-Xe bond. Such a result is consistent with experimental observations about our studied cases.
Table 3. The H-Xe bond orders: bH−Xe (electron-sharing), bH−Xe (dative), and the total [bH−Xe (electron-sharing) + bH−Xe (dative)], for the monomer and their hydrogen-bonded complex, with the changes shown in parentheses.
It becomes clear that the bonding between H and Xe in HXeY must meet two conditions. First, the resonance bonding between H− Xe+-Y and H-Xe+ Y− is a necessary condition. Second, for the H-Xe bond order, including contribution of these two resonance structures is essential. These two conditions are, in effect, consistent with emphases in original CS bonding concept paper (Shaik et al., 1992). The mixed covalent-ionic description, such as F· ·F↔F− F+, is an essential feature of CS bonding, wherein most, if not the entire, bond energy is provided by the covalent-ionic resonance energy. And both the covalent and ionic structures must be treated explicitly and on an equal footing. Thus, we conclude that the H-Xe bond in HXeY is a charge-shift bond.
It should be noted, however, that neither the H-Xe bond length nor its bond strength can be reliable probes for addressing the question on whether the bond between H and Xe in HXeY should be classified as a CS bond or a classical covalent bond.
H-Xe CS Bonding Character
Now that the H-Xe bond in HXeY is a charge-shift bond, a new and unique form bonding, one question which arises is whether the unusual H-Xe blue shift is its CS bonding character in IR spectroscopy.
Deep analyses on the data in Table 3 show that the complexation leads to an increase of bH−Xe (dative) while bH−Xe (electron-sharing) decreases, for most Xe cases. The competition of these two opposite factors results in the strengthened H-Xe bonds, corresponding to the blue shifts. Here, we want to point out that unlike other Xe complexes, the HXeCl···HCCH complex shows a difference in the dominant factor. It is the electron-sharing factor that dominates the H-Xe frequent shifts. The overall effect of two opposing factors is still an enhancement of the H-Xe bond and a blue shift for the H-Xe stretching frequency. These analyses reflect that the monomer-to-complex blue shifts of H-Xe stretching modes for HXeY species should be attributed to the balance of dative and electron-sharing covalency in H-Xe bonds, corresponding to H− Xe+-Y and H-Xe+ Y− resonance structures. In brief, the H-Xe frequent shifts in HXeY···H2O/HCl/HBr/HI/HCCH/HCN hydrogen-bonded complexes is controlled by a balance of two factors acting in opposite directions.
All in all, blue shifts of H-Xe vibrational frequencies are controlled by a balance of two opposing factors for dative and electron-sharing covalency. The blue shift in our studied complexes can be seen as a normal spectroscopic phenomenon. It is natural for us to conclude that the H-Xe blue shift is the H-Xe CS bonding character in IR spectroscopy.
Summary and Conclusions
The H-Xe bonding in HXeY has been a debated question in the chemical community. The usual answer is that it is classically covalent in character, or to be exact, it is electron-sharing. The unusual blue shifts of complexes HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CN, CCH) definitely reflect the unusual features of H-Xe bonding in noble gas hydrides. Via observed blue-shifted phenomena, we have computationally investigated the H-Xe bonding in HXeY from an NBO/NRT perspective.
We establish that the resonance bonding between H-Xe+ Y− and H− Xe+-Y is an essential feature of the H-Xe bonding, and that the H-Xe bonding in these two resonance structures must be considered explicitly. Specifically, its bonding includes the donor-acceptor interaction in H− Xe+-Y, as well as the electron-sharing interaction in H-Xe+ Y−; the H-Xe bond order is contributed to by these two resonance structures. We confirm that the H-Xe bond in HXeY is not a classical covalent bond but a charge-shift bond, and that H-Xe blue shifts is a normal spectroscopic phenomenon.
Our conclusions are (1) the H-Xe bond in HXeY is a charge-shift bond. (2) The H-Xe blue shift in its hydrogen-complexes is its CS bonding character in IR spectroscopy.
The first conclusion is consistent with ab initio VB methods' insight into the F-Xe CS bonding in XeF2 (Braïda and Hiberty, 2013). But we note a little difference in the understanding of CS bonding mechanism. We stress the point that the H-Xe CS bonding is due to the resonance between H-Xe+:I− and H:− Xe+-I, based on the natural Lewis structures' concept, whereas Shaik et al. (1992) think that CS bonding is due to strong mixing between the covalent structure and the ionic structure, based on Pauling-type Lewis structures' concepts.
Finally, we want to point out that H-Xe CS bonding is significantly different from some two-structure resonance bonding, such as hydrogen-bonding, although there is a formal resemblance in their resonance description. We have noticed that the bond order of hydrogen-bonding in literature (Jiao and Weinhold, 2019) is only considered a contribution from one resonance structure; the other does not contribute to bond orders of hydrogen-bonding at all. Why is CS bonding not important in hydrogen-bonding? More studies are on the way to answer such a question and to generalize CS bonding models. We believe that more surprise will be gained via NBO/NRT methods, in particular, new-type NBO-based NRT methods on analyses for larger CS bonding species (Glendening et al., 2019).
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.
Author Contributions
GZ, CS, and DC made contributions to the design of the study. YS, XZ, JS and LF performed the calculations. GZ and YS wrote the draft of the manuscript. All the authors contributed to the manuscript revision, read and approved the submitted version.
Funding
This work was supported by the Natural Science Foundation of Shandong Province (No. ZR2015BM025).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Professor Frank Weinhold of the University of Wisconsin-Madison for his continuous help. Especially his help and important suggestions during the last few months, which helped us make great improvements to our manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2020.00277/full#supplementary-material
References
Alabugin, I. V., Manoharan, M., and Weinhold, F. (2004). Blue-shifted and red-shifted hydrogen bonds in hypervalent rare-gas FRg–H···Y sandwiches. J. Phys. Chem. A 108, 4720–4730. doi: 10.1021/jp049723l
Braïda, B., and Hiberty, P. C. (2013). The essential role of charge-shift bonding in hypervalent prototype XeF2. Nat. Chem. 5, 417–422. doi: 10.1038/nchem.1619
Braïda, B., Ribeyre, T., and Hiberty, P. C. (2014). A valence bond model for electron-rich hypervalent species: application to SFn (n = 1, 2, 4), PF5, and ClF3. Chem. Eur. J. 20, 9643–9649. doi: 10.1002/chem.201402755
Corani, A., Domanskaya, A., Khriachtchev, L., Räsänen, M., and Lignell, A. (2009). Matrix-isolation and ab initio study of the HKrCl···HCl complex. J. Phys. Chem. A 113, 10687–10692. doi: 10.1021/jp9044622
Duarte, L., and Khriachtchev, L. (2017). An aromatic noble-gas hydride: C6H5CCXeH. Sci. Rep. 7:3130. doi: 10.1038/s41598-017-02869-9
Feldman, V. I., Sukhov, F. F., Orlov, A. Y., and Tyulpina, I. V. (2003). Experimental evidence for the formation of HXeCCH: the first hydrocarbon with an inserted rare-gas atom. J. Am. Chem. Soc. 125, 4698–4699. doi: 10.1021/ja034585j
Feller, D. (1996). The role of databases in support of computational chemistry calculations. J. Comput. Chem. 17, 1571–1586. doi: 10.1002/(SICI)1096-987X(199610)17:13<1571::AID-JCC9>3.0.CO;2-P
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., et al. (2009). Gaussian 09, Revision B.01. Wallingford CT: Gaussian Inc.
Glendening, E. D., Badenhoop, J. K., Reed, A. E., Carpenter, J. E., Bohmann, J. A., Morales, C. M., et al. (2001). NBO 5.0., Theoretical Chemistry Institute. Madison: University of Wisconsin.
Glendening, E. D., Badenhoop, J. K., and Weinhold, F. (1998). Natural resonance theory: III. chemical applications. J. Comput. Chem. 19, 628–646. doi: 10.1002/(SICI)1096-987X(19980430)19:6<628::AID-JCC5>3.0.CO;2-T
Glendening, E. D., Landis, C. R., and Weinhold, F. (2013b). NBO 6.0: natural bond orbital analysis program. J. Comput. Chem. 34, 1429–1437. doi: 10.1002/jcc.23266
Glendening, E. D., Landis, C. R., and Weinhold, F. (2019). Resonance theory reboot. J. Am. Chem. Soc. 141, 4156–4166. doi: 10.1021/jacs.8b12336
Glendening, E. D., and Weinhold, F. (1998a). Natural resonance theory: I. general formalism. J. Comput. Chem. 19, 593–609. doi: 10.1002/(SICI)1096-987X(19980430)19:6<593::AID-JCC3>3.0.CO;2-M
Glendening, E. D., and Weinhold, F. (1998b). Natural resonance theory: II. natural bond order and valency. J. Comput. Chem. 19, 610–627. doi: 10.1002/(SICI)1096-987X(19980430)19:6<610::AID-JCC4>3.0.CO;2-U
Glendening, E. D., Badenhoop, J. K., Reed, A. E., Carpenter, J. E., Bohmann, J. A., Morales, C. M., et al. (2013a). NBO 6.0. Madison, WI: Theoretical Chemistry Institute, University of Wisconsin. Available online at: http://nbo6.chem.wisc.edu/ (accessed March 24, 2018).
Grandinetti, F. (2018). Noble Gas Chemistry: Structure, Bonding, and Gas-phase Chemistry. Weinheim: Wiley-VCH. doi: 10.1002/9783527803552
Head-Gordon, M., Pople, J. A., and Frisch, M. J. (1988). MP2 energy evaluation by direct methods. J. Chem. Phys. Lett. 153, 503–506. doi: 10.1016/0009-2614(88)85250-3
Jana, G., Pan, S., Merino, G., and Chattaraj, P. K. (2017). MNgCCH (M = Cu, Ag, Au; Ng = Xe, Rn): the first set of compounds with M–Ng–C bonding motif. J. Phys. Chem. A 121, 6491–6499. doi: 10.1021/acs.jpca.7b04993
Jana, G., Pan, S., Merino, G., and Chattaraj, P. K. (2018a). Noble gas inserted metal acetylides (metal = Cu, Ag, Au). J. Phys. Chem. A 122, 7391–7401. doi: 10.1021/acs.jpca.8b05404
Jana, G., Pan, S., Osorio, E., Zhao, L., Merino, G., and Chattaraj, P. K. (2018b). Cyanide–isocyanide isomerization: stability and bonding in noble gas inserted metal cyanides (metal = Cu, Ag, Au). Phys. Chem. Chem. Phys. 20, 18491–18502. doi: 10.1039/C8CP02837K
Jiao, Y., and Weinhold, F. (2019). What is the nature of supramolecular bonding? comprehensive NBO/NRT picture of halogen and pnicogen bonding in RPH2···IF/FI complexes (R = CH3, OH, CF3, CN, NO2). Molecules. 24:2090. doi: 10.3390/molecules24112090
Juarez, R., Zavala-Oseguera, C., Jimenez-Halla, J. O. C., Bickelhaupt, F. M., and Merino, G. (2011). Radon hydrides: structure and bonding. Phys. Chem. Chem. Phys. 13, 2222–2227. doi: 10.1039/C0CP01488E
Khriachtchev, L., Räsänen, M., and Gerber, R. B. (2009). Noble-gas hydrides: new chemistry at low temperatures. Acc. Chem. Res. 42, 183–191. doi: 10.1021/ar800110q
Khriachtchev, L., Tanskanen, H., Lundell, J., Pettersson, M., Kiljunen, H., and Räsänen, M. (2003). Fluorine-free organoxenon chemistry: HXeCCH, HXeCC, and HXeCCXeH. J. Am. Chem. Soc. 125, 4696–4697. doi: 10.1021/ja034485d
Landis, C. R., and Weinhold, F. (2013). 3c/4e σ-type long-bonding: a novel transitional motif toward the metallic delocalization limit. Inorg. Chem. 52, 5154–5166. doi: 10.1021/ic4000395
Last, I., and George, T. F. (1988). Electronic states of the XenHCl systems in gas and condensed phases. J. Chem. Phys. 89, 3071–3078. doi: 10.1063/1.454963
Lein, M., Frunzke, J., and Frenking, G. (2004). Christian Klixbull Jørgensen and the nature of the chemical bond in HArF. Struct. Bond. 181–191. doi: 10.1007/b11312
Lignell, A., Lundell, J., Khriachtchev, L., and Räsänen, M. (2008). Experimental and computational study of HXeY···HX complexes (X, Y = Cl and Br) : an example of exceptionally large complexation effect. J. Phys. Chem. A 112, 5486–5494. doi: 10.1021/jp801363u
Lundell, J., Cohen, A., and Gerber, R. B. (2002). Quantum chemical calculations on novel molecules from xenon insertion into hydrocarbons. J. Phys. Chem. A 106, 11950–11955. doi: 10.1021/jp026777r
Pan, S., Jana, G., Merino, G., and Chattaraj, P. K. (2019). Noble-noble strong union: gold at its best to make a bond with a noble gas atom. ChemistryOpen 8. 173–187. doi: 10.1002/open.201800257
Pérez-Peralta, N., Juárez, R., Cerpa, E., Bickelhaupt, F. M., and Merino, G. (2009). Bonding of xenon hydrides. J. Phys. Chem. A 113, 9700–9706. doi: 10.1021/jp903266a
Pettersson, M., Lundell, J., and Räsänen, M. (1995). Neutral rare-gas containing charge-transfer molecules in solid matrices. I. HXeCl, HXeBr, HXeI, and HKrCl in Kr and Xe. J. Chem. Phys. 102, 6423–6431. doi: 10.1063/1.469357
Pettersson, M., Lundell, J., and Räsänen, M. (1999). New rare-gas-containing neutral molecules. Eur. J. Inorg. Chem. 1999, 729–737. doi: 10.1002/(SICI)1099-0682(199905)1999:5<729::AID-EJIC729>3.0.CO;2-M
Pyykkö, P. (2015). Additive covalent radii for single-, double-, and triple-bonded molecules and tetrahedrally bonded crystals: a summary. J. Phys. Chem. A 119, 2326–2337. doi: 10.1021/jp5065819
Rahm, M., Hoffmann, R., and Ashcroft, N. W. (2016). Atomic and ionic radii of elements 1-96. Chem. Eur. J. 22, 14625–14632. doi: 10.1002/chem.201602949
Räsänen, M., Khriachtchev, L., Pettersson, M., Runeberg, N., and Lundell, J. (2000). A stable argon compound. Nature. 406, 874–876. doi: 10.1038/35022551
Saha, R., Jana, G., Pan, S., Merino, G., and Chattaraj, P. K. (2019). How far can one push the noble gases towards bonding?: a personal account. Molecules. 24:2933. doi: 10.3390/molecules24162933
Schuchardt, K. L., Didier, B. T., Elsethagen, T., Sun, L., Gurumoorthi, V., Chase, J., et al. (2007). Basis set exchange: a community database for computational sciences. J. Chem. Inf. Model. 47, 1045–1052. doi: 10.1021/ci600510j
Shaik, S., Danovich, D., Galbraith, J. M., Braida, B., Wu, W., and Hiberty, P. C. (2020). Charge-shift bonding: a new and unique form of bonding. Angew. Chem. Int. Ed. 59, 984–1001. doi: 10.1002/anie.201910085
Shaik, S., Danovich, D., Silvi, B., Lauvergnat, D. L., and Hiberty, P. C. (2005). Charge-shift bonding—a class of electron-pair bonds that emerges from valence bond theory and is supported by the electron localization function approach. Chem. Eur. J. 11, 6358–6371. doi: 10.1002/chem.200500265
Shaik, S., Danovich, D., Wu, W., and Hiberty, P. C. (2009). Charge-shift bonding and its manifestations in chemistry. Nat. Chem. 1, 443–449. doi: 10.1038/nchem.327
Shaik, S., Maitre, P., Sini, G., and Hiberty, P. C. (1992). The charge-shift bonding concept. Electron-pair bonds with very large ionic-covalent resonance energies. J. Am. Chem. Soc. 114, 7861–7866. doi: 10.1021/ja00046a035
Tsuge, M., Berski, S., Räsänen, M., Latajka, Z., and Khriachtchev, L. (2013). Experimental and computational study of the HXeI···HY complexes (Y = Br and I). J. Chem. Phys. 138:104314. doi: 10.1063/1.4794309
Tsuge, M., Berski, S., Räsänen, M., Latajka, Z., and Khriachtchev, L. (2014). Matrix-isolation and computational study of the HXeY···H2O complexes (Y = Cl, Br, and I). J. Chem. Phys. 140:044323. doi: 10.1063/1.4862692
Weinhold, F. (2012). Natural bond orbital analysis: a critical overview of relationships to alternative bonding perspectives. J. Comput. Chem. 33, 2363–2379. doi: 10.1002/jcc.23060
Weinhold, F. (2013). NBOPro 6.0. Madison: Theoretical Chemistry Institute, U. Wisconsin. Available online at: http://nbo6.chem.wisc.edu/~nbopro_adv.htm (accessed March 24, 2018).
Weinhold, F., and Klein, R. A. (2012). What is a hydrogen bond? mutually consistent theoretical and experimental criteria for characterizing H-bonding interactions. Mol. Phys. 110, 565–579. doi: 10.1080/00268976.2012.661478
Weinhold, F., and Klein, R. A. (2014). What is a hydrogen bond? resonance covalency in the supramolecular domain. Chem. Educ. Res. Pract. 15, 276–285. doi: 10.1039/C4RP00030G
Weinhold, F., and Landis, C. R. (2005a). “Natural resonance structures and weightings,” in Valency and Bonding: A Natural Bond Orbital Donor-acceptor Perspective. (New York, NY: Cambridge University Press), 32–36.
Weinhold, F., and Landis, C. R. (2005b). “The NBO donor-acceptor picture of hypervalency: hyperbonds,” in Valency and Bonding: A Natural Bond Orbital Donor-acceptor Perspective. (New York, NY: Cambridge University Press), 281–285.
Weinhold, F., Landis, C. R., and Glendening, E. D. (2016). What is NBO analysis and how is it useful? Int. Rev. Phys. Chem. 35, 399–440. doi: 10.1080/0144235X.2016.1192262
Zhang, G., Fu, L., Li, H., Fan, X., and Chen, D. (2017a). Insight into the bonding mechanism and the bonding covalency in noble gas-noble metal halides: an NBO/NRT investigation. J. Phys. Chem. A 121, 5183–5189. doi: 10.1021/acs.jpca.7b02047
Zhang, G., Li, H., Weinhold, F., and Chen, D. (2016). 3c/4e -type long-bonding competes with ω-bonding in noble-gas hydrides HNgY (Ng = He, Ne, Ar, Kr, Xe, Rn; Y = F, Cl, Br, I): a NBO/NRT perspective. Phys. Chem. Chem. Phys. 18, 8015–8026. doi: 10.1039/C5CP07965A
Zhang, G., Song, J., Fu, L., Tang, K., Su, Y., and Chen, D. (2018). Understanding and modulating the high-energy properties of noble-gas hydrides from their long-bonding: an NBO/NRT investigation on HNgCO+/CS+/OSi+ and HNgCN/NC (Ng = He, Ar, Kr, Xe, Rn) molecules. Phys. Chem. Chem. Phys. 20, 10231–10239. doi: 10.1039/C8CP00306H
Zhang, G., Yue, H., Weinhold, F., Wang, H., Li, H., and Chen, D. (2015). Resonance character of Cu/Ag/Au bonding in small molecule···M-X (X=F, Cl, Br, CH3, CF3) complexes. Chem. Phys. Chem. 16, 2424–2431. doi: 10.1002/cphc.201500211
Zhang, G., Zhang, S., and Chen, D. (2017b). Long-bonding in HNgCN/NC (Ng = noble gas) molecules: an NBO/NRT investigation. J. Phys. Chem. A 121, 5524–5532. doi: 10.1021/acs.jpca.7b03177
Keywords: chemical bonding, charge-shift bonds, electron-pair bonds, resonance bonding, noble-gas hydride, hypervalent molecule, blue shifts, NBO/NRT methods
Citation: Zhang G, Su Y, Zou X, Fu L, Song J, Chen D and Sun C (2020) Charge-Shift Bonding in Xenon Hydrides: An NBO/NRT Investigation on HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CCH, CN) via H-Xe Blue-Shift Phenomena. Front. Chem. 8:277. doi: 10.3389/fchem.2020.00277
Received: 04 February 2020; Accepted: 23 March 2020;
Published: 23 April 2020.
Edited by:
Sudip Pan, University of Marburg, GermanyReviewed by:
Felice Grandinetti, University of Tuscia, ItalyGabriel Merino, Centro de Investigacion y de Estudios Avanzados–Unidad Mérida, Mexico
Lili Zhao, Nanjing Tech University, China
Copyright © 2020 Zhang, Su, Zou, Fu, Song, Chen and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guiqiu Zhang, Z3F6aGFuZyYjeDAwMDQwO3NkbnUuZWR1LmNu; Chuanzhi Sun, c3VuY3omI3gwMDA0MDtzZG51LmVkdS5jbg==