BeH2与HX(X=F,Cl, Br,I)形成的二氢键复合物的结构特征与本质

对BeH2与HX(X=F, Cl, Br, I)形成的二氢键复合物的结构特征及本质进行了探讨. 在MP2/6-311++G(3d,3p)水平优化、频率验证, 得到复合物的分子结构, 用分子间距离及电子密度拓扑理论确认BeH2与卤化氢已形成了二氢键型复合物. 在MP2/6-311++G(3d, 3p)水平下进行基函数重叠误差(BSSE)校正后的结合能在-14.468 kJ·mol-1到-5.464 kJ·mol-1之间.用对称匹配微扰理论(SAPT)对复合物的结合能进行分解, 结果表明, BeH2…HX二氢键复合物中静电能对总吸引能的贡献都是最主要的, 但交换排斥能、诱导能、色散能对总结合能的贡献也很重要. 从BeH2…HF到BeH2…HI, 诱导能对总吸引能的贡献从37.8%逐渐减小到24.0%. 而色散能对总吸引能的贡献从BeH2…HF体系中的16.0%逐渐增加到BeH2…HI体系中的33.8%.     

Proton-transfer and H2-elimination reactions of main-group hydrides EH4- (E = B, Al, Ga) with alcohols

The reaction of the isostructural anions of group 13 hydrides EH4- (E = B, Al, Ga) with proton donors of different strength (CH3OH, CF3CH2OH, and CF3OH) was studied with different theoretical methods [DFT/B3LYP and second-order M?ller-Plesset (MP2) using the 6-311++G(d,p) basis set]. The results show the general mechanism of the reaction: the dihydrogen-bonded (DHB) adduct (EH...HO) formation leads through the activation barrier to the next concerted step of H2 elimination and alkoxo product formation. The structures, interaction energies (calculated by different approaches including the energy decomposition analysis), vibrational E-H modes, and electron-density distributions were analyzed for all of the DHB adducts. The transition state (TS) is the dihydrogen complex stabilized by a hydrogen bond with the anion [EH3(eta2-H2)...OR-]. The single exception is the reaction of BH4- with CF3OH exhibiting two TSs separated by a shallow minimum of the BH3(eta2-H2)...OR- intermediate. The structures and energies of all of the species were calculated, leading to the establishment of the potential energy profiles for the reaction. A comparison is made with the mechanism of the proton-transfer reaction to transition-metal hydrides. The solvent influence on the stability of all of the species along the reaction pathway was accounted for by means of polarizable conductor calculation model calculations in tetrahydrofuran (THF). Although in THF the DHB intermediates, the TSs, and the products are destabilized with respect to the separated reactants, the energy barriers for the proton transfer are only slightly affected by the solvent. The dependence of the energies of the DHB complexes, TSs, and products as well as the energy barriers for the H2 release on the central atom and the proton donor strength is also discussed.     

Intermolecular HH vibrations of dihydrogen bonded complexes H3EH(-)...HOR in the low-frequency region: theory and IR spectra

The results of DFT calculations of harmonic and anharmonic frequencies of the dihydrogen bonded (DHB) complexes H 3EH (-)...HOR (E = B, Al, Ga and HOR = CH 3OH, CF 3CH 2OH) in gas phase and in low polar medium (by CPCM model) in comparison with the partners are presented. Normal coordinate analysis of the low-frequency modes was carried out to assign the new vibrations induced by DHB formation by the potential energy distribution values. Among them, the intermolecular H...H stretching vibrations only have individual modes. The influence of central atom mass and isotope and the strength of the proton donor effects were determined. The systems convenient for IR studies were chosen from the calculation predictions. The spectral investigation was made on the BH 4 (-)/ROH complexes (ROH = CH 2FCH 2OH (MFE), CF 3CH 2OH (TFE), (CF 3) 2CHOH (HFIP)). The results of temperature dependence, isotope substitution, and influence of the proton-donor strength studies agree with the theoretical conclusions. Combination of experimental and theoretical approaches allowed determining for the first time the intermolecular stretching mode characterizing intrinsic DHB vibrations.     

Theoretical study on stabilities of multiple hydrogen bonded dimers

A series of self-constituted multiple hydrogen bonded (MHB) complexes has been investigated systematically by density functional theory (PBE1PBE /6-31G**), the Morokuma energy decomposition method (HF/6-31G**) and MP2 (6-31G** and 6-311++G**) calculation. We have discovered that (i) for doubly hydrogen bonded (DHB) complexes, both the interaction energy and stability increase with the charge transfer energy; (ii) for quadruple hydrogen bonded (QHB) complexes, cooperativity is the most important factor determining stability of the complex: stronger cooperative energy correlates well with larger interaction energy and thus more stable complex and vice versa; (iii) correlation energy plays an important role in intermolecular interactions. The correlation energy, mainly consisting of dispersive energy, also exhibits cooperativity in MHB dimers: positive for M-aadd and generally negative for other complexes.     

Theoretical study of dihydrogen bonds between (XH)2, X = Li, Na, BeH, and MgH, and weak hydrogen bond donors (HCN, HNC, and HCCH)

The dihydrogen-bonded (DHB) complexes formed by (XH)2, with X = Li, Na, BeH, and MgH, with one, two, and four protonic molecules (HCN, HNC, and HCCH) have been studied. These complexes have been compared to those of the XH monomers with the same hydrogen bond donor molecules. The energetic results have been rationalized based on the electrostatic potential of the isolated hydridic systems. The electron density properties have been analyzed within the AIM methodology, both at the bond critical points and the integrated values at the atomic basins. Exponential relationships between several properties calculated at the bond critical points (rho,nabla2rho, lambdai, G, and V) and variation of integrated properties (energy, charge, and volume) vs the DHB distance have been obtained.     

Theoretical study of (CH...C)- hydrogen bonds in CH(4-n)X(n) (X = F, Cl; n = 0, 1, 2) systems complexed with their homoconjugate and heteroconjugate carbanions

This work deals with a theoretical study of the (CH...C)- hydrogen bonds in CH4, CH3X, and CH2X2 (X = F, Cl) complexed with their homoconjugate and heteroconjugate carbanions. The properties of the complexes are calculated with the B3LYP method using the 6-311++G(d,p) or 6-311++G(2df,2p) basis sets. The deprotonation enthalpies (DPE) of the CH bond or the proton affinities of the carbanions (PA(C-) are calculated as well. All the systems with the exception of the CH4...CHCl2(-) one are characterized by a double minimum potential. In some of the complexes, the (CH(b)...C)- hydrogen bond is linear. In other systems, such as CH3F...CH2F- and CH3F...CHF2(-), there is a large departure from linearity, the systems being stabilized by electrostatic interactions between the nonbonded H of the neutral molecule and the F atom of the carbanion. In the transition state, the (CH(b)...C)- bond is linear, and there is a large contraction of the intermolecular C...C distance. The binding energies vary within a large range, from -1.4 to -11.1 kcal mol(-1) for the stable complexes and -8.6 to -44.1 kcal mol(-1) for the metastable complexes. The energy barriers to proton transfer are between 5 and 20 kcal mol(-1) for the heteroconjugate systems and between 3.8 and 8.3 kcal mol(-1) for the homoconjugate systems. The binding energies of the linear complexes depend exponentially on 1.5DPE - PA(C-), showing that the proton donor is more important than the proton acceptor in determining hydrogen bond strength. The NBO analysis indicates an important electronic reorganization in the two partners. The elongations of the CH bond resulting from the interaction with the carbanion depend on the occupation of the sigma*(CH(b)) antibonding orbitals and on the hybridization of the C bonded to H(b). The frequency shifts of the nu(CH)(A1) stretching vibration range between 15 and 1150 cm(-1). They are linearly correlated to the elongation of the CH(b) bond.     

HCN(HNC)与NH3, H2O和HF分子间相互作用的理论研究

在MP2/aug-cc-pVTZ水平上, 对HCN(HNC)与NH3, H2O和HF分子间可能存在的氢键型复合物进行了全自由度能量梯度优化, 通过在相同水平上的频率验证分析发现了稳定的分子间相互作用形式是HCN(HNC)作为质子供体或作为质子受体形成的复合物. 基组重叠误差对总相互作用能的影响均小于3.34 kJ/mol. 通过自然键轨道(NBO)分析, 研究了单体和复合物中的原子电荷和电荷转移对分子间相互作用的影响. 对称性匹配微扰理论(SAPT, Symmetry Adapted Perturbation Theory)能量分解结果表明, 在分子间相互作用中, 静电作用与诱导作用占主导地位, 而诱导作用与复合物的电荷转移之间具有良好的正相关性.     

4(3H)-嘧啶酮及其类似物互变异构的理论研究

采用密度泛函理论,在B3LYP/6-311G**基组水平上,计算并考察了4(3H)-嘧啶酮及其类似物(5-氟-4(3H)-嘧啶酮、4-巯基嘧啶和5-氟-4-巯基嘧啶)醇式结构和酮式结构进行结构互变质子迁移过程中的2种可能途径:(a)分子内质子迁移;(b)水助质子迁移.计算结果表明,途经b所需要的活化能较小.研究还表明,氢键在降低反应活化能方面起着重要的作用.     

Dihydrogen bonding vs metal-σ interaction in complexes between H2 and metal hydride

The complexes formed by hydrogen with metal hydrides (LiH, NaH, BeH(2), MgH(2), BH(3), AlH(3), Li(2)H(2), Na(2)H(2), Be(2)H(4), and Mg(2)H(4)) have been theoretically studied at the MP2/aug-cc-pVTZ, MP2/aug-cc-pVQZ and CCSD(T)/aug-cc-pVTZ//CCSD/aug-cc-pVTZ levels of theory. The hydrogen molecule can act as a Lewis acid or base. In the first case, a dihydrogen bonded complex is obtained and in the second an interaction between the σ-bond of the hydrogen molecule and an empty orbital of the metal atoms is found. Quantum theory of atoms in molecules and natural bond orbitals methods have been applied to analyze the intermolecular interactions. Additionally, the cooperativity effects are analyzed for selected complexes with two H(2) molecules where both kinds of interactions exist simultaneously.     

CH3SH…HOO开壳型氢键复合物的结构及其电子密度拓扑性质

在DFT-B3LYP/6-311++G**水平下求得CH3SH…HOO复合物势能面上的稳定构型. 计算结果表明, 在HOO以其O8—H7作为质子供体与CH3SH分子中的S5原子为质子受体形成的氢键复合物1和2中, O8—H7明显被“拉长”, 且其伸缩振动频率发生显著的红移, 红移值分别为330.1和320.4 cm-1; 在CH3SH分子以其S5—H6作为质子供体与HOO的端基O9原子为质子受体形成的氢键复合物3和4中, 也存在类似的情况, 但S5—H6伸缩振动频率红移不大. 经MP2/6-311++G**水平计算的4种复合物含BSSE校正的相互作用能分别为-20.81, -20.10, -4.46和-4.52 kJ/mol. 自然键轨道理论(NBO)分析表明, 在CH3SH…HOO复合物1和2中, 引起H7—O8键长增加的因素包括两种电荷转移, 即孤对电子n1(S5)→σ*(H7—O8)和孤对电子n2(S5)→σ*(H7—O8), 其中后者为主要作用. 在复合物3和4中也有相似的电荷转移情况, 但轨道间的相互作用要弱一些. AIM理论分析结果表明, 4个复合物中的S5…H7间和O9…H6间都存在键鞍点, 且其Laplacian量▽2ρ(r)都是很小的正值, 说明这种相互作用介于共价键和离子键之间, 偏静电作用为主.     

How short can the H...H intermolecular contact be? New findings that reveal the covalent nature of extremely strong interactions

Ab initio calculations at the MP2/6-311++G(d,p) and MP2/aug-cc-pVDZ//MP2/aug-cc-pVTZ levels have been performed for the following complexes: H2OH+...HBeH, H2OH+...HBeBeH, H2OH+...HBeF, HClOH+...HBeH, Cl2OH+...HBeH, and Cl2OH+...HBeF. For all dimers considered, extremely short H...H intermolecular contacts (1.0-1.3 A) were obtained. These are the shortest intermolecular distances which have ever been reported, with binding energies within the range of 13.7-24.3 kcal/mol (MP2/aug-cc-pVDZ//MP2/aug-cc-pVTZ level). The interaction energies of the complexes analyzed were also extrapolated to the complete basis set (CBS) limit. To explain the nature of such strong interactions, the Bader theory was applied, and the characteristics of the bond critical points (BCPs) were analyzed. It was pointed out that for the major part of the H...H contacts considered here the Laplacian of the electron density at H...H BCP is negative indicating the partly covalent nature of such a connection. The term "covalent character of the hydrogen bond" used sometimes in recent studies is discussed. An analysis of the interaction energy components for dihydrogen bonded systems considered indicates that in contrast to conventional hydrogen bonded systems the attractive electrostatic term is outweighed by the repulsive exchange energy term and that the higher order delocalization energy term is the most important attractive term.     

Cooperativity in hydrogen-bonded interactions: ab initio and "atoms in molecules" analyses

The H(2)CO...(HF)(n) (n = 1, ..., 9) complexes were investigated using the MP2 method and the following basis sets: 6-311++G(d,p), aug-cc-pVDZ and aug-cc-pVTZ. It was found that the cooperativity effect enhances significantly the F-H...O hydrogen bond; in some of cases one can detect the covalent nature of hydrogen bonding. To deepen the nature of the interactions investigated, the scheme of decomposition of the interaction energy was applied; for stronger H-bonds where the cooperativity is more important, the delocalization energy term increases. The ratio of delocalization energy to electrostatic energy increases for stronger hydrogen bonds where the proton...acceptor distance is shorter. The Bader theory was also applied, and it was found that for stronger H-bonds the electronic energy density at the proton...acceptor bond critical point is negative and may be attributed to the partly covalent interaction.     

Theoretical analysis of the hydrogen bond of imidazolium C(2)-H with anions

The intermolecular interaction energies of ion pairs of imidazolium-based ionic liquids were studied by MP2/6-311G level ab initio calculations. Although the hydrogen bond between the C(2) hydrogen atom of an imidazolium cation and anion has been regarded as an important interaction in controlling the structures and physical properties of ionic liquids as in the cases of conventional hydrogen bonds, the calculations show that the nature of the C(2)-H...X interaction is considerably different from that of conventional hydrogen bonds. The interaction energies of the imidazolium cation with neighboring anions in the four crystals of ionic liquids were calculated. The size of the interaction is determined mainly by the distance between the imidazolium ring and anion. The calculated interaction energy is nearly inversely proportional to the distance, which shows that the charge-charge interaction is the dominant interaction in the attraction. The orientation of the anion relative to the C(2)-H bond does not greatly affect the size of the interaction energy. Calculated interaction energy potentials of 1,3-dimethylimidazolium tetrafluoroborate ([dmim][BF(4)]) complexes show that the C(2)-H bond does not prefer to point toward a fluorine atom of the BF(4). This shows that the C(2)-H...X hydrogen bond is not essential for the attraction.     

Electronic properties of multifurcated bent hydrogen bonds CH3...Y and CH2...Y

H-bonding angle angleYHX has an important effect on the electronic properties of the H-bond Y...HX, such as intra- and intermolecular hyperconjugations and rehybridization, and topological properties of electron density. We studied the multifurcated bent H-bonds of the proton donors H3CZ (Z = F, Cl, Br), H2CO and H2CF2 with the proton acceptors Cl(-) and Br(-) at the four high levels of theory: MP2/6-311++G(d,p), MP2/6-311++G(2df,2p), MP2/6-311++G(3df,3pd) and QCISD/6-311++G(d,p), and found that they are all blue-shifted. These complexes have large interaction energies, 7-12 kcal mol(-1), and large blue shifts, delta r(HC) = -0.0025 --0.006 A and delta v(HC) = 30-90 cm(-1). The natural bond orbital analysis shows that the blue shifts of these H-bonds Y...HnCZ are mainly caused by three factors: rehybridization; indirect intermolecular hyperconjugation n(Y) -->sigma*(CZ), in that the electron density from n(Y) of the proton acceptor is transferred not to sigma*(CH), but to sigma*(CZ) of the donor; intramolecular hyperconjugation n(Z) -->sigma*(CH), in that the electron density in sigma*(CH) comes back to n(Z) of the donor such that the occupancy in sigma*(CH) decreases. The topological properties of the electron density of the bifurcated H-bonds Y...H2CZ are similar to those of the usual linear H-bonds, there is a bond critical point between Y and each hydrogen, and a ring critical point inside the tetragon YHCH. However, the topological properties of electron density of the trifurcated H-bonds Y...H3CZ are essentially different from those of linear H-bonds, in that the intermolecular bond critical point, which represents a closed-shell interaction, is not between Y and hydrogen, but between Y and carbon.     

Infrared studies on hydrogen bonding interaction between acrylic esters with 1-dodeconal and phenol in carbon tetrachloride

The influence of the hydroxyl group in 1-dodeconal and phenol on the carbonyl vibration of representative compounds belonging to methyl methacrylate, ethyl methacrylate and butyl methacrylate in carbon tetrachloride has been studied by Fourier transform infrared spectroscopic method. The integrated intensities and change in dipole moments for O–H and C=O bonds were calculated. The formation constants of the 1 : 1 complexes have been calculated using Nash's method. The values of the formation constant and Gibbs energy vary with the ester chain length, which suggests that the strengths of the intermolecular hydrogen bonds (O–H ··· O=C) are dependent on the alkyl group of the acrylic ester and the results show that the proton accepting ability of acrylic esters is in the order methyl methacrylate < ethyl methacrylate < butyl methacrylate. The strength of the intermolecular C=O: HO bonds is also shown to be dependent on the basicity of the C=O group of acrylic esters and the acidity of the proton donor.     

pi-H...O hydrogen bonds: multicenter covalent pi-H interaction acts as the proton-donating system

The MP2 method and the Pople-style basis sets 6-311++G(d,p), 6-311++G(2df,2pd), and 6-311++G(3df,3pd) were used to perform calculations on H3O+...C2H2 and C2H3+...C2H2 complexes and related species. Hydrogen bonds existing for the analyzed complexes were investigated as well as related pi-H...O --> pi...H-O and pi-H...pi --> pi...H-pi proton-transfer processes. For some of the complexes analyzed the multicenter pi-H interaction possessing the properties of a covalent bond acts as a proton donor; more generally it is classified as the Lewis acid. The quantum theory of "atoms in molecules" (QTAIM) was also applied to deepen the nature of these interactions in terms of characteristics of bond critical points. The pi-H...O, O-H...pi, and pi-H...pi interactions analyzed here may be classified as hydrogen bonds since their characteristics are the same as or at least similar to those of typical hydrogen bonds. H...pi interactions are common in crystal structures of organic and organometallic compounds. The analyses performed here show a continuum of such interactions since there are H...pi contacts possessing the characteristics of weak intermolecular interactions on the one hand and pi-H multicenter covalent bonds on the other. Ab initio and QTAIM results support the latter statements.     

Unveiling the Noncovalent Interaction of Thiazol-2-ylidene and Its Derivatives as N-heterocyclic Carbene with Different Proton Donor Molecules

The importance of noncovalent interaction has gained attention in various domains covering drug and novel catalyst design. The present study mainly characterizes the role of hydrogen bond (H-bond) and other intermolecular interactions in different (1 : 1) complex analogues formed between the N-aryl-thiazol-2-ylidene (YR) and five proton donor (HX) molecules. The analysis of the singlet-triplet energy gap ( ) confirmed the stability of the singlet state for this class of N-aryl-thiazol-2-ylidenes than the triplet state. The interaction energy values of the YR-HX complexes follow the order: YR-NH₃<YR-HCN<YR-H₂O<YR-MeOH<YR-HF. In addition, substituting the H-atom of the N−H bond with bulky groups (−R) leads to an increase in the interaction energy of the YR-HX complexes. Hence, it was found that the replacement of N-atom in N-heterocyclic carbene (NHC) by S-atom forming N-aryl-thiazol-2-ylidene results in comparable intermolecular interactions with proton donor molecules similar to imidazole-2-ylidene (NHC). The current study enlightened the role of noncovalent interactions in carbene complexes with proton donor molecules. We hope that our work on carbene chemistry will pave the way for its application in the designing and synthesis of efficient catalysts.     

One-dimensional BeH2 polymers: infrared spectra and theoretical calculations

Laser-ablated beryllium atoms react with H2 upon co-condensation in excess hydrogen and neon to form BeH2 and (BeH2)2, which are identified through isotopic substitution and DFT calculations. Unreacted Be atoms isolated in solid neon or hydrogen are excited to the 1P0 state and react further with H2 to enhance the BeH2 and (BeH2)2 concentrations and produce (BeH2)n polymers. The series of strong infrared-active parallel Be-H-Be bridge-bond stretching modes observed for (BeH2)n polymers suggests one-dimensional structures, and this conclusion is supported by DFT calculations. The computed polymerization energy per BeH2 unit is about 33 kcal/mol.     

Cooperative effects of noncovalent bonds to the Br atom of halogen-bonded H(3)N[ellipsis (horizontal)]BrZ and HCN[ellipsis (horizontal)]BrZ (Z = F, Br) complexes

A series of complexes formed between halogen-bonded H(3)N∕HCN[ellipsis (horizontal)]BrZ (Z = Br, F) dimers and H(3)N∕HCN[ellipsis (horizontal)]BrZ[ellipsis (horizontal)]XY (XY = HF, ClF, BeH(2), LiF) trimers were investigated at the MP2 and B3LYP levels of theory using a 6-31++G(d,p) basis set. Optimized structures, interaction energies, and other properties of interest were obtained. The addition of XY to the H(3)N∕HCN[ellipsis (horizontal)]BrZ dyad leads to enhanced intermolecular binding with respect to the isolated monomers. This enhanced binding receives contributions from the electrostatic and inductive forces between the constituent pairs, with, in some instances, substantial three-body non-additive contributions to the binding energy. It was found that the XY = LiF interaction causes the greatest distortion of the H(3)N∕HCN[ellipsis (horizontal)]BrZ halogen bond from the preferred linear orientation and also provides the strongest binding energy via the nonadditive energy.     

Intermolecular interactions in uracil–nitrous acid complexes: structures,binding energy,topological properties,and nuclear magnetic resonance study

Molecular interactions between uracil and nitrous acid (U–NA) [C₄N₂O₂H₄? NO₂H] have been studied using B3LYP, B3PW91, and MP2 methods with different basis sets. The optimized geometries, harmonic vibrational frequencies, charge transfer, topological properties of electron density, nucleus‐independent chemical shift (NICS), and nuclear magnetic resonance one‐ and two‐bonds spin–spin coupling constants were calculated for U–NA complexes. In interaction between U and NA, eight cyclic complexes were obtained with two intermolecular hydrogen bonds N(C)HU…N(O) and OH_NA…O_U. In these complexes, uracil (U) simultaneously acts as proton acceptor and proton donor. The most stable complexes labeled, UNA1 and UNA2, are formed via NH bond of U with highest acidity and CO group of U with lowest proton affinity. There is a relationship between hydrogen bond distances and the corresponding frequency shifts. The solvent effect on complexes stability was examined using B3LYP method with the aug‐cc‐pVDZ basis set by applying the polarizable continuum model (PCM). The binding energies in the gas phase have also been compared with solvation energies computed using the PCM. Natural bond orbital analysis shows that in all complexes, the charge transfer takes place from U to NA. The results predict that the Lone Pair (LP)(O)_U → σ*(O? H) and LP(N(O)_NA → σ*(N(C)? H)_U donor–acceptor interactions are most important interactions in these complexes. Atom in molecule analysis confirms that hydrogen bond contacts are electrostatic in nature and covalent nature of proton donor groups decreases upon complexation. The relationship between spin–spin coupling constant (^1hJ_H…_Y and ^2hJ_H…_Y) with interaction energy and electronic density at corresponding hydrogen bond critical points and H‐bonds distances are investigated. NICS used for indicating of aromaticity of U ring upon complexation. © 2013 Wiley Periodicals, Inc.