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)能量分解结果表明, 在分子间相互作用中, 静电作用与诱导作用占主导地位, 而诱导作用与复合物的电荷转移之间具有良好的正相关性.     

Symmetry-adapted perturbation theory analysis of the N...HX hydrogen bonds

The main aim of the study was the detailed investigation of the interaction energy decomposition in dimers and trimers containing N...HX bonds of different types. The study of angular dependence of interaction energy terms partitioned according to the symmetry-adapted perturbation theory (SAPT) was performed for the dimers containing N...HX bonds as mentioned above: ammonia-HX (X = F, Cl, Br) and pyridine-HF complexes. It was found that the electrostatic and induction terms exhibit strong angular dependence, while the exchange contributions are less affected. The dispersion terms are virtually nondirectional. In addition, the three-body SAPT interaction energy analysis for the mixed acid-base NH3...(HF)2 and (NH3)2...HF trimers revealed strong differences between interactions of similar strength but different types (i.e., hydrogen bond and general electrostatic interaction). The importance of the induction terms for the nonadditivity of the interaction energy in strongly polar systems was confirmed.     

The boron-boron single bond in diborane(4) as a non-classical electron donor for hydrogen bonding

An ab initio study of an isomer of diborane(4) [B(2)H(4)] has been carried out at MP2/aug-cc-pVTZ to investigate the ground-state properties of this unusual molecule, a derivative of which has been described in the recent literature. The geometric, electronic and orbital characteristics of B(2)H(4)(4) have been analyzed using AIM, NBO, and ELF methodologies. A region with a high concentration of electron density is located near and along the B-B bond, on the opposite side of this bond relative to the bridging H atoms. This site serves as an electron-donor site to electrophiles, resulting in hydrogen-bonded complexes of B(2)H(4) with proton donors HF, HNC, HCl, HCN, and HCCH, and a van der Waals complex with H(2). These complexes have C(2v) symmetry and stabilization energies that vary from 2 to 27 kJ mol(-1). The SAPT2 energy decomposition analysis shows that the relative importance of the various terms that contribute to the interaction energy depends on the strength of the interaction.     

Cooperative effects, strengths of hydrogen bonds, and intermolecular interactions in circular cis, trans-cyclotriazane clusters (n = 3-8)

Based on Becke's three parameter functional [J. Chem. Phys. 98, 5648 (1993)] of density functional theory (DFT) with the correlation of Lee-Yang-Parr [Phys. Rev. B 37, 785 (1988)] (DFT/B3LYP), the natural bond orbital (NBO) analysis, the Bader's theory of atoms in molecule (AIM), our calculations indicate that as cluster size (n) increases, the n-dependent cooperative changes in the lengths of the N...H H bonds (HBs) and N-H bonds, the N-H stretching frequencies and intensities, and the n(N)-->sigma*(N-H) charge transfers are observed to be pervasive in the circular cis, trans-cyclotriazane clusters (n = 3-8), which is very different from the linear cis, trans-cyclotriazane clusters reported in previous work. According to the NBO and AIM theories, the cooperativity of the intermolecular n(N)-->sigma*(N-H) interaction leads to the n-dependent N...H contractions. In this way, the stronger N...H bond is formed, as reflected in the increase in their rho(r(cp)) values. This increased electron density is translated into the improved capacity to concentrate electrons at the HB bond critical point (BCP), i.e., a higher potential energy V(r(cp)). On the other hand, stronger repulsion is also activated to counteract the contraction, which is reflected in the increased G(r(cp)) value that gives the tendency of the system to dilute electrons at the HB BCP. In terms of the three-body symmetry-adapted perturbation theory (three-body SAPT), the induction nonadditivity accounts for up to 97% of the nonadditive energy in the circular trimer. It can believed that the marked cooperativity of the n(N)-->sigma*(N-H) interactions is of nonadditive induction in nature. The N...H formation and nature of cooperativity in the circular clusters differ from those in the linear clusters that have been reported previously. According to the SAPT(DFT) method which is a combination of SAPT with the asymptotically corrected DFT, the cis, trans-cyclotriazane systems should contain remarkable dispersion interactions. However, the short-range dispersion cannot be reproduced thoroughly by DFT/B3LYP. A quantum cluster equilibrium model illustrates the neglected dispersion energies and the nonadditive energies can affect markedly the properties of the liquid consisting of the circular clusters.     

The chemical nature of the (plus sign in circle)C-H...X- (X=Cl or Br) interaction in imidazolium halide ionic liquids

The C-H...X (X=Cl or Br) interaction is traditionally characterized as a relatively weak interaction. However, this interaction becomes very strong in the imidazolium-based halide ionic liquids [J. Phys. Chem. 123, 174501 (2005)]. This strong interaction had been attributed to the electrostatic interaction between the imidazolium cation and the halide anion. In this paper, the chemical nature of the (plus sign in circle)C-H...Cl(-) and ( plus sign in circle)C-H...Br(-) interactions is investigated by atoms in molecules (AIM) and natural bond orbital (NBO) analyses. The AIM calculations indicate that in the EmimX complexes, the (plus sign in circle)C-H...Cl(-) and (plus sign in circle)C-H...Br(-) interactions have some covalent character, especially the (plus sign in circle)C-H...Cl(-) interaction. Mulliken, ChelpG charge, and natural bond orbital population analyses for these two kinds of interactions indicate that the charge transfer is important in the interaction of the cation with the anion. In addition, the NBO analysis demonstrated that the stabilization energy is due to an n-->sigma(C-H) (*) orbital interaction. However, in the Emim2X and Emim3X complexes, the calculated results suggested a dominant electrostatic character for the (plus sign in circle)C-H...Cl(-) and (plus sign in circle)C-H...Br(-) interactions.     

Unusual halogen-bonded complex FBrdelta+...delta+BrF and hydrogen-bonded complex FBrdelta+...delta+HF formed by interactions between two positively charged atoms of different polar molecules

Using ab initio calculations, the authors' predicted for the first time that the halogen-bonded complex FBrdelta+...delta+BrF and hydrogen-bonded complex FBrdelta+...delta+HF formed by the interactions between two positively charged atoms of different polar molecules can be stable in gas phase. It shows that halogen bond or hydrogen bond not only exists between oppositely charged atoms but also between like-charged atoms. That the attraction arising from the special halogen bond or hydrogen bond can exceed the electrostatic repulsion between two contact positively charged atoms stabilizes the complex. Of course, from the point of view of physics they can consider the interactions in FBrdelta+...delta+BrF and FBrdelta+...delta+HF as mainly the sum of the long range molecular interactions, namely, electrostatic, induction, and dispersion with some short-range repulsion. They found that the intermolecular electron correlation contribution representing dispersion interaction plays a crucial role in the stabilities of seemingly repulsive complexes FBrdelta+...delta+BrF and FBrdelta+...delta+HF.     

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.     

Mutual interplay between pnicogen–π and tetrel bond in PF<Subscript>3</Subscript>⊥X–Pyr…SiH<Subscript>3</Subscript>CN complexes: NMR,SAPT, AIM,NBO, and MEP analysis

The mutual interplay between pnicogen–π and tetrel bond in the formation of PF₃⊥X–Pyr…SiH₃CN ternary complexes has been investigated via a computational chemistry at MP2/aug-cc-pVDZ level of theory. We proved by computational NMR data the effect of electron-withdrawing and electron-donating substituents on ^1tJ(N-Si) across ¹⁵N...³⁵Si tetrel bonds was investigated at M06-2X/aug-cc-pVDZ levels of theory in PF₃⊥CN–Pyr…SiH₃CN complex. The nature of the interactions has been studied by means of symmetry-adapted perturbation theory (SAPT) and molecular electrostatic potentials (MEP). The electrostatic interaction played a major role in the change of tetrel bond interaction strength in the X–Pyr…SiH₃CN binary systems, whereas the change of pnicogen–π strength in the PF₃⊥X–Pyr complexes was caused jointly by the dispersion interactions. Energy decomposition indicates that the percentage of the electrostatic term in the tetrel bond system constitutes in the total attractive binding energies, while the percentage of the dispersion term in the pnicogen bonding constitutes in the attractive binding energies. In addition, atoms in molecules (AIM) and natural bond orbital (NBO) analyses were also performed to unveil the mechanism of these interactions in the title complexes.     

Study on the nature of interaction of BrCl with HF,H2O,and NH3

By counterpoise‐corrected optimization method, the interactions of BrCl with the first‐row hydrides (HF, H₂O, NH₃) have been investigated at the MP2/6–311++G(3d,3p) level. To understand that the X? Br‐type (X = F, O, N) structure is more stable than the corresponding hydrogen‐bonded structure in these complexes, the electronic properties were also investigated. Symmetry‐adapted perturbation theory (SAPT) analysis has been carried out to understand the nature of the weak hydrogen bond and X? Br‐type interactions. On the other hand, for the weak hydrogen‐bonded complexes and the X? Br‐type complexes charges transfer is well correlated with the total induction energies. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

On the nature of stabilization in weak, medium, and strong charge-transfer complexes: CCSD(T)/CBS and SAPT calculations

Weak, medium, and strong charge-transfer (CT) complexes containing various electron donors (C(2)H(4), C(2)H(2), NH(3), NMe(3), HCN, H(2)O) and acceptors (F(2), Cl(2), BH(3), SO(2)) were investigated at the CCSD(T)/complete basis set (CBS) limit. The nature of the stabilization for these CT complexes was evaluated on the basis of perturbative NBO calculations and DFT-SAPT/CBS calculations. The structure of all of the complexes was determined by the counterpoise-corrected gradient optimization performed at the MP2/cc-pVTZ level, and most of complexes possess a linear-like contact structure. The total stabilization energies lie between 1 and 55 kcal/mol and the strongest complexes contain BH(3) as an electron acceptor. When ordering the electron donors and electron acceptors on the basis of these energies, we obtain the same order as that based on the perturbative E2 charge-transfer energies, which provides evidence that the charge-transfer term is the dominant energy contribution. The CCSD(T) correction term, defined as the difference between the CCSD(T) and MP2 interaction energies, is mostly small, which allows the investigation of the CT complexes of this type at the "cheap" MP2/CBS level. In the case of weak and medium CT complexes (with stabilization energy smaller than about 15 kcal/mol), the dominant stabilization originates in the electrostatic term; the dispersion as well as induction and δ(HF) terms covering the CT energy contribution are, however, important as well. For strong CT complexes, induction energy is the second (after electrostatic) most important energy term. The role of the induction and δ(HF) terms is unique and characteristic for CT complexes. For all CT complexes, the CCSD(T)/CBS and DFT-SAPT/CBS stabilization energies are comparable, and surprisingly, it is true even for very strong CT complexes with stabilization energy close to 50 kcal/mol characteristic by substantial charge transfer (more than 0.3 e). It is thus possible to conclude that perturbative DFT-SAPT analysis is robust enough to be applied even for dative-like complexes with substantial charge transfer.     

SiH_4与AB型卤素互化物之间的反转氢键

采用从头算方法对SiH4与AB型卤素互化物(ClF、BrF、IF、ICl、IBr、BrCl)形成的复合物的结构特征及本质进行了探讨.在MP2/6-311++G(3d,3p)水平上优化复合物的分子结构,并进行频率验证.通过分子间距离、自然键轨道(NBO)净电荷迁移数及分子图,确认SiH4与卤素互化物形成反转氢键型复合物.在MP2/6-311++G(3d,3p)水平上进行基组重叠误差(BSSE)校正后的作用能为-5.113--9.468kJ·mol-1.用对称匹配微扰理论(SAPT)对作用能进行分解,结果显示,诱导能对总吸引能的贡献在55.0%到72.2%之间,是最主要的贡献部分,静电能和色散能对总吸引能的贡献都小于25.0%.     

Nature of a hydride-halogen bond. A SAPT-, QTAIM-, and NBO-based study

The nature of a hydride-halogen bond is investigated by means of high-level quantum mechanical calculations expended with symmetry-adapted perturbation theory (SAPT), quantum theory of atoms in molecules (QTAIM), and natural bond orbital (NBO) methods. As model hydride-halogen bonded systems complexes between either LiH or HBeH and either XCF(3) or XCCH (X = F, Cl, Br, I) are used. It is shown that the formation of a hydride-halogen bond leads to the elongation of the R(δ+)-H(δ-) hydride bond, which is accompanied by the blue shift of the ν(R-H) stretching vibration frequency and the increase of the IR intensity of this mode. All these effects, although untypical in the case of, e.g., hydrogen bonds, can be considered as rather typical for hydride-halogen bonded systems. The decomposition of the interaction energy based on the SAPT method clearly indicates the dominant role of the induction term, thus the inductive nature of a hydride-halogen bond in opposition to previous findings. NBO-based analysis indicates the charge transfer from the hydride molecule to the more remote parts of the halogen donor and that the elongation of the R-H bond is caused by the charge outflow from the σ(RH) bonding orbital.     

Charged versus Neutral Hydrogen‐Bonded Complexes: Is There a Difference in the Nature of the Hydrogen Bonds?

A theoretical study on some carboxylic acid dimers formed by positively or negatively charged molecules has been carried out by using DFT methods. The resulting dimers possess either a charge of +2 or ?2. In addition, the corresponding neutral complexes have also been considered. The electron density distribution described by the atoms in molecules and the natural bond orbital methods, as well as the electric field maps of the systems, have been analyzed and compared without finding significant differences between the neutral and ionic complexes. The interaction energy along the dissociation path of the charged dimers shows both a local minimum and a local maximum, defining a stability region between them. When this energetic profile is recalculated by removing the repulsion between the charged groups, it resembles to those of the neutral molecules. Hence, the characteristics of the charged dimers are similar to those of the neutral ones: the addition of a repulsion term for the charged groups permits to retrieve the energetic profiles dependence with the distance in the charged system. The interacting quantum atom (IQA) method has been used to calculate the interaction energy terms, including the classic Coulombic term between the whole molecules and the corresponding of the carboxylic acid groups. The IQA results show repulsive electrostatic interactions when the whole molecules are considered in the ionic complexes, but attractive ones between the carboxylic groups in both neutral and ionic complexes.     

苯-卤素(X2, X=F, Cl, Br, I)相互作用本质的对称性匹配微扰理论研究

在MP2水平下对被定义为"电荷转移复合物(CTC)"的苯(C6H6)-卤素分子X2(X＝F, Cl, Br, I)相互作用体系进行了量子化学研究. 在优化所得C6H6-X2(X＝F, Cl, Br, I)复合物的平衡几何结构中, 卤素分子X2接近垂直指向苯环上碳-碳双键的中心. 自然键轨道(NBO)分析结果表明, 苯-卤素体系中电荷转移的数量很少. 对称性匹配微扰理论(Symmetry-adapted perturbation theory, SAPT) 能量分解结果显示, 在4个复合物体系中, 静电作用的贡献相对较小(只占总吸引作用的20%左右), 对于C6H6-F2体系, 色散作用是其主要吸引作用, 对于C6H6-Cl2, C6H6-Br2和C6H6-I2 体系, 诱导作用则是其主要的吸引作用, 从F到I, 色散作用逐渐减弱, 诱导作用逐渐增强, 表明在电子相关水平上将苯-卤素体系称为"电荷转移复合物"的说法并不确切.     

Theoretical investigation on charge‐assisted halogen bonding interactions in the complexes of bromocarbons with some anions

A series of dimeric complexes formed between bromocarbon molecules and two anions (Br^? and CN^?) have been investigated by using MP2 method. The quantum theory of atoms in molecules (QTAIM) and the second‐order perturbation natural bond orbital (NBO) approaches were applied to analyze the electron density distributions of these complexes and to explore the nature of charge‐assisted halogen bonding interactions. As anticipated, these interactions are significantly stronger relative to the corresponding neutral ones. The results derived from ab initio calculations described herein reveal a major contribution from the electrostatic interaction on the stability of the systems considered. Beside the electrostatic interaction, the charge‐transfer force and the second‐order orbital interaction also play an important role in the formation of the complexes, as a NBO analysis suggested. The presence of halogen bonds in the complexes has been identified in terms of the QTAIM methodology, and several linear relationships have been established to provide more insight into charge‐assisted halogen bonding interactions. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Prediction and characterization of a chalcogen-hydride interaction with metal hybrids as an electron donor in F2CS-HM and F2CSe-HM (M = Li, Na, BeH, MgH, MgCH3) complexes

A novel type of σ-hole bonding has been predicted and characterized in F(2)CS-HM and F(2)CSe-HM (M = Li, Na, BeH, MgH) complexes at the MP2/aug-cc-pVTZ level. This interaction, termed a chalcogen-hydride interaction, was analyzed in terms of geometric, energetic and spectroscopic features of the complexes. It exhibits similar properties to hydrogen bonding and halogen bonding. The methyl group in metal hydrides makes a positive contribution to the formation of chalcogen-hydride bonded complexes. In the F(2)CSe-HLi-OH(2) complex, the chalcogen-hydride bonding shows synergetic effects with lithium bonding. These complexes have been analyzed with the atoms in molecules (AIM) theory and symmetry adapted perturbation theory (SAPT) method. The results show that the chalcogen-hydride bonding is dominated with an electrostatic interaction.     

传统氢键(X—H…Y)还是X—H…π相互作用？ ——无机苯与卤化氢相互作用本质的量子化学研究

在MP2/6-311＋＋G**水平对无机苯(B3N3H6)与卤化氢HX (X＝F, Cl, Br, I)相互作用体系进行了系统研究. 结果表明在B3N3H6-HX (X＝F, Cl, Br, I)体系的平衡几何结构中, HX的H原子倾向于指向B3N3H6环上的N原子, 且从HF到HI相互作用强度依次减弱. 与苯-卤化氢体系比较, 除与HF相互作用B3N3H6较C6H6强外, 其余体系B3N3H6均较C6H6弱(结合能数值相差4 kJ/mol左右). 对称匹配微扰理论(SAPT)能量分解结果说明静电、诱导和色散力对描述B3N3H6-卤化氢体系的相互作用都很重要, 从HF到HI静电能占总吸引作用能的百分比逐渐减少, 色散能占总吸引作用能的百分比逐渐增加, 这种变化趋势与苯-卤化氢体系比较类似, 表明B3N3H6与卤化氢的相互作用随着卤素原子序数的递增, 传统氢键作用趋势减弱, X—H…π相互作用趋势增强.     

Halogen as halogen-bonding donor and hydrogen-bonding acceptor simultaneously in ring-shaped H3N·X(Y)·HF (X = Cl, Br and Y = F, Cl, Br) complexes

A series of ring-shaped molecular complexes formed by H(3)N, HF and XY (X = Cl, Br and Y = F, Cl, Br) have been investigated at the MP2/aug-cc-pVTZ level of theory. Their optimized geometry, stretching mode, and interaction energy have been obtained. We found that each complex possesses two red-shifted hydrogen bonds and one red-shifted halogen bond, and the two hydrogen bonds exhibit strong cooperative effects on the halogen bond. The cooperativity among the NH(3)···FH, FH···XY and H(3)N···XY interactions leads to the formations of these complexes. The AIM analysis has been performed at the CCSD(T)/aug-cc-pVQZ level of theory to examine the topological characteristics at the bond critical point and at the ring critical point, confirming the coexistence of the two hydrogen bonds and one halogen bond for each complex. The NBO analysis carried out at the B3LYP/aug-cc-pVTZ level of theory demonstrates the effects of hyperconjugation, hybridization, and polarization coming into play during the hydrogen and halogen bonding formations processes, based on which a clockwise loop of charge transfer was discovered. The molecular electrostatic potential has been employed to explore the formation mechanisms of these molecular complexes.     

Intermolecular Weak Interactions in HTeXH Dimers (X=O,S, Se,Te): Hydrogen Bonds,Chalcogen–Chalcogen Contacts and Chiral Discrimination

A theoretical study of the HTeXH (X=O, S, Se and Te) monomers and homodimers was carried out by means of second‐order Møller‐Plesset perturbation theory (MP2) computational methods. In the case of monomers, the isomerization energy from HTeXH to H₂Te=X and H₂X=Te (X=O, S, Se, and Te) and the rotational transition‐state barriers were obtained. Due to the chiral nature of these compounds, homo and heterochiral dimers were found. The electron density of the complexes was characterized with the atoms‐in‐molecules (AIM) methodology, finding a large variety of interactions. The charge transfer within the dimers was analyzed by means of natural bond orbitals (NBO). The density functional theory‐symmetry adapted perturbation theory (DFT‐SAPT) method was used to compute the components of the interaction energies. Hydrogen bonds and chalcogen–chalcogen interactions were characterized and their influence analyzed concerning the stability and chiral discrimination of the dimers.