卤化氢HX（X=F,Cl,Br,I）与C2H2相互作用本质的量子化学研究

在MP2水平上,采用全电子基组,对C2H2与HX(X=F,Cl,Br,I)相互作用进行了研究.构型优化同时进行频率验证,得到4个T型结构的稳定复合物,相互作用能在-12.761～-7.086kJ/mol之间.自然键轨道(NBO)与分子中的原子(AIM)理论分析表明,形成复合物分子间的电荷转移量都很少,最大仅为0.009a.u.,作用强度与氢键类似.对称性匹配微拢理论(SAPT)能量分解数据表明,对于C2H2…HX(X=F,Cl,Br,I)体系,从F到I,静电作用逐渐减弱,色散作用逐渐增强;相互作用能中对吸引能的贡献主要为静电能和色散能,二者之和占到80%以上,诱导能所占的比例很小,卤化氢与乙炔分子间相互作用的本质为静电作用和色散作用.     

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)都是很小的正值, 说明这种相互作用介于共价键和离子键之间, 偏静电作用为主.     

苯-卤素(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, 色散作用逐渐减弱, 诱导作用逐渐增强, 表明在电子相关水平上将苯-卤素体系称为"电荷转移复合物"的说法并不确切.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要.在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型,B3LYP方法得到了三种构型(复合物Ⅰ,Ⅱ和Ⅲ),而MP2方法只能得到一种构犁(复合物Ⅱ).在复合物Ⅰ和Ⅲ中,HSO单元中的1H原子作为质子供体.与O3分子中的端基O原子作为质子受体相互作用,形成红移氢键复合物;而在复合物Ⅱ中,虽与复合物Ⅰ和Ⅲ中具有相间的质子供体和质子受体,却形成了蓝移氢键复合物.B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重甍误差(BSSE)和零点振动能(ZPVE)校正,其值在-3.37到-4.55 kJ·mol-1之间.采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

N2O…HOCl复合物结构和性质的理论研究

分别在DFT-B3LYP和MP2/6-311++G**水平上求得HOCl + N2O复合物势能面上的六种(S1, S2, S3, S4, S5和S6)和四种(S1, S2, S4和S5)构型. 频率分析表明,其中的S1和S3为过渡态,其它为稳定构型. 在复合物S3, S5 和S6中,HOCl 单体的σ*(5O-6H)作为质子供体,与N2O单体中作为质子受体的3O原子相互作用,形成氢键结构,而在氢键复合物S2中, 质子受体为N2O单体中的端1N原子;复合物S1中,HOCl分子的σ*(5O-4Cl)作为质子供体与N2O分子中的端1N原子(质子受体)相互作用,形成卤键结构,而复合物S4中的卤键结构的质子受体为N2O分子中的端3O原子. 经B3LYP/6-311++G**水平上的计算,考虑了基组重叠误差(BSSE)校正的单体间相互作用能在-1.56 ~ -8.73 kJ·mol-1之间. 采用自然键轨道理论(NBO)对两种单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要. 在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型, B3LYP方法得到了三种构型(复合物I, II和III), 而MP2方法只能得到一种构型(复合物II). 在复合物I和III中, HSO单元中的1H原子作为质子供体, 与O3分子中的端基O原子作为质子受体相互作用, 形成红移氢键复合物; 而在复合物II中, 虽与复合物I和III中具有相同的质子供体和质子受体, 却形成了蓝移氢键复合物. B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-3.37到-4.55 kJ·mol-1之间. 采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

HOCl…HCOCl复合物的结构和电子性质

在DFT-B3LYP/6-311++G**水平上求得HOCl+HCOCl复合物势能面上的四种稳定构型(S1, S2, S3和S4). 其中, 在复合物S1和S3中, HOCl单体的5H原子作为质子供体, 与HCOCl单体中作为质子受体的1O原子相互作用, 形成红移氢键复合物; 在复合物S4中, HOCl单体的7Cl原子作为质子供体, 与HCOCl单体中作为质子受体的1O原子相互作用, 形成红移卤键复合物; 而在复合物S2中, 同时存在2C—3H…6O蓝移氢键和4Cl…5O相互作用. 在MP2/6-311++G**水平上计算的单体间的相互作用能考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-5.05与-14.76 kJ·mol-1之间. 采用自然键轨道理论(NBO)对两种单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.     

HOCl…HCOCl复合物的结构和电子性质

在DFT-B3LYP/6-311++G**水平上求得HOCl+HCOCl复合物势能面上的四种稳定构型(S1,S2,S3和S4).其中,在复合物S1和S3中,HOCl单体的5H原子作为质子供体,与HCOCl单体中作为质子受体的10原子相互作用,形成红移氢键复合物;在复合物S4中,HOCl单体的7Cl原子作为质子供体,与HCOCl单体中作为质子受体的IO原子相互作用,形成红移卤键复合物;而在复合物S2中,同时存在2C-3H…6O蓝移氢键和4Cl…5O相互作用.在MP2/6-311++G**水平上计算的单体间的相互作用能考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正,其值在-5.05与-14.76 kJ·mol-1之间.采用自然键轨道理论(NBO)对两种单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.     

CaH_2…HY（Y=CH_3,C_2H_3,C_2H,CN,NC）体系双氢键性质的理论研究

在B3LYP/6-311＋＋G（3df,2p）及MP2/6-311＋＋G（3df,2p）水平上探讨了以CaH2为质子受体及多种常见质子供体构成的CaH2…HY（Y=CH3,C2H3,C2H,CN,NC）弱相互作用体系的几何结构以及所形成的双氢键的性质.MP2方法计算得到线型的CaH2分子结构,在B3LYP方法计算下优化得到的CaH2部分则是一个非常明显的弯曲结构.NBO分析计算表明,线型结构（MP2）中Ca原子杂化轨道主要源自于sp杂化而弯曲结构（B3LYP）则主要为sd杂化,这与之前报道的相关研究结论相一致.基于不同CaH2分子结构计算得到的CaH2…HY体系的H…H间距离、相互作用能以及电荷转移性质等非常接近.MP2和B3LYP方法计算得到的CaH2…HY体系H…H间距离均处于0.15～0.24nm,相互作用能处于1～44kJmol-1范围.CaH2…HY体系的相互作用能变化趋势为CaH2…CH4-CaH2…C2H4-CaH2…C2H2-CaH2…HCN-CaH2…HNC,H…H距离的变化规律与相互作用能变化趋势相反,即键能越大H…H距离越小.相较于HY单体,形成的CaH2…HY复合物体系的质子供体键H-Y的键长增大,质子供体的键长增加量△r（H-Y）与H…H间距离间有这样的关系：H…H间距离越小,△r（H-Y）越大.双氢键AIM拓扑数据表明形成CaH2…HY复合物后,Ca-H键临界点的电子密度（ρc）较单体时减小且与双氢键相互作用能成线性关系,即复合物中Ca-H键临界点的ρc数值越小同时双氢键键能越大.NPA结果表明,CaH2…HY体系中分子间的电荷转移总量随键能增大而增大.通过H…H间距离、相互作用能以及AIM分析,我们认为CaH2…CH4和CaH2…C2H4体系可划分为弱的范德华相互作用,而CaH2…C2H2,CaH2…HCN和CaH2…HNC体系形成的双氢键强度上相当于常规氢键,且最值得注意的是CaH2分子结构的差异性对CaH2…HY体系性质的影响微乎其微.     

杯[4]吡咯-阴离子体系相互作用的理论研究——杯[4]吡咯-卤素阴离子体系的密度泛函研究

用密度泛函B3LYP/LANL2DZ方法对自由杯[4]吡咯的最低能量构象和卤素阴离子-杯[4]吡咯复合物进行了计算研究.结果表明,杯[4]吡咯与卤素离子能通过彼此间的相互作用形成复合物,并且这种相互作用在本质上应为分子间的氢键相互作用;杯[4]吡咯与卤素阴离子形成的复合物在构型变化、电荷转移、前线轨道及其作用方式、成键布居以及能量和热力学参数等方面均按元素周期律有规律地变化,杯[4]吡咯与卤素阴离子间的相互作用沿元素周期依次减少.     

Ab initio study of ternary radical–molecule complexes between HCN(HNC) and HO(HS) species

MP2 calculations with the cc-pVTZ basis set were used to analyze the intermolecular interactions in ternary radical–molecule complexes between HCN(HNC) and HO(HS) species, in gas phase and in water media. Particular attention was given to parameters such as the cooperative energies and many-body interaction energies. The results indicate that hydrogen bonding between two HCN(HNC) molecules gives more stability to triads than hydrogen bonding between HCN(HNC) and OH(SH) species. The electronic properties of the complexes were analyzed using the parameters derived from the atoms in molecules methodology. The water media has an enhancing influence on the stabilities of studied clusters versus the ones obtained in gas phase.     

Comparative analysis of blue‐shifted hydrogen bond versus conventional hydrogen bond in methyl radical complexes

Methyl radical complexes H₃C…HCN and H₃C…HNC have been investigated at the UMP2(full)/aug‐cc‐pVTZ level to elucidate the nature of hydrogen bonds. To better understand the intermolecular H‐bond interactions, topological analysis of electron density at bond critical points (BCP) is executed using Bader's atoms‐in‐molecules (AIM) theory. Natural bond orbital (NBO) analysis has also been performed to study the orbital interactions and change of hybridization. Theoretical calculations show that there is no essential difference between the blue‐shift H‐bond and the conventional one. In H₃C…HNC complex, rehybridization is responsible for shortening of the N? H bond. The hyperconjugative interaction between the single electron of the methyl radical and N? H antibonding orbital is up to 7.0 kcal/mol, exceeding 3.0 kcal/mol, the upper limit of hyperconjugative n(Y)→σ*(X–H) interaction to form the blue‐shifted H‐bond according to Alabugin's theory. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Structures,vibrational frequencies,topologies, and energies of hydrogen bonds in cysteine‐formaldehyde complexes

The hydrogen bonding interactions between cysteine (Cys) and formaldehyde (FA) were studied with density functional theory regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules and natural bond orbital analyses were employed to elucidate the interaction characteristics in the Cys‐FA complexes. The intramolecular hydrogen bonds (H‐bonds) formed between the hydroxyl and the N atom of cysteine moiety in some Cys‐FA complexes were strengthened because of the cooperativity. Most of intermolecular H‐bonds involve the O atom of cysteine/FA moiety as proton acceptors, while the strongest H‐bond involves the O atom of FA moiety as proton acceptor, which indicates that FA would rather accept proton than providing one. The H‐bonds formed between the CH group of FA and the S atom of cysteine in some complexes are so weak that no hydrogen bonding interactions exist among them. In most of complexes, the orbital interaction of H‐bond is predominant during the formation of complex. The electron density (ρ_b) and its Laplace (?²ρ_b) at the bond critical point significantly correlate with the H‐bond parameter δR, while a linearly relationship between the second‐perturbation energy E(2) and ρ_b has been found as well. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

The properties of weak and strong dihydrogen-bonded D-H...H-A complexes.

The properties of six dihydrogen-bonded (DHB) dimers with the BeH2 molecule as a proton acceptor were calculated by MP2, CCSD(T) and B3LYP methods. The structural, energetic and spectroscopic parameters are presented and analyzed in terms of their possible correlation with the interaction energy and the intermolecular H...H separation. The symmetry-adapted perturbation theory (SAPT) calculations were performed to gain more insight into the nature of the H...H interactions. The studied complexes are divided into three groups based on the calculated intermolecular distances and the interaction energies which range from approximately -1 to -42 kJ mol(-1). The analysis of the interaction energy components indicates that, in contrast to conventional hydrogen bonds, the induction energy is the most important term in the BeH2NH4+ complex. On the other hand, there is no sharp boundary between the DHB complexes classified as hydrogen bonded and van der Waals systems. The complexation-induced changes in vibrational frequencies and in proton shielding constants show a relationship with the interaction energy. The values of the 2hJXH and 3hJBeX coupling constants correlate well with the interaction energy and with the intermolecular distance.     

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.     

Halogen bond versus hydrogen bond: The many‐body interactions approach

The analysis of interrelation between halogen bond and hydrogen bond in the (RX)(HNC)(HCN) complexes (R = CH₃, CF₃ and X = Cl, Br, I) was performed on the basis of DFT calculations. Both two‐body additive contributions and three‐body nonadditive contributions to the total interaction energy were discussed. QTAIM was used for topological analysis of electron density. Additionally, QTAIM analysis of electron density was performed for both two‐ and three‐body complexes. The electron charge transfer in trimers showed the dual character of the fragment with halogen atom involved into the investigated interactions—it acts as Lewis acid and Lewis base, depending on the type of interaction considered. The effect of cooperativity of X‐ and H‐bonding was assessed on the basis of many‐body interaction energy and electron density analysis. Additionally, an alternative two‐body model with the same situation (in the context of intermolecular interactions) is investigated. The anti‐cooperative effect was found also for this model.     

Theoretical Study on the Hydrogen Bonding Interactions in Complexes of 5‐Hydroxytryptamine with Water

The energies, geometries and harmonic vibrational frequencies of 1:1 5‐hydroxytryptamine‐water (5‐HT‐H₂O) complexes are studied at the MP2/6‐311++G(d,p) level. Natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM) analyses and the localized molecular orbital energy decomposition analysis (LMO‐EDA) were performed to explore the nature of the hydrogen‐bonding interactions in these complexes. Various types of hydrogen bonds (H‐bonds) are formed in these 5‐HT‐H₂O complexes. The intermolecular C4H5^5‐HT···O^w H‐bond in HTW3 is strengthened due to the cooperativity, whereas no such cooperativity is found in the other 5‐HT‐H₂O complexes. H‐bond in which nitrogen atom of amino in 5‐HT acted as proton donors was stronger than other H‐bonds. Our researches show that the hydrogen bonding interaction plays a vital role on the relative stabilities of 5‐HT‐H₂O complexes.     

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.     

Hydrogen-bond interactions in hydrated 6-selenoguanine tautomers: a theoretical study

A comprehensive theoretical investigation has been performed to study the six most stable complexes of isolated, mono, and hexahydrated 6-selenoguanine tautomers. The ground state geometries are studied at the density-functional theory and Møller–Plesset Perturbation theory implementing the 6-311++G (2d, 2p) basis set. The intermolecular distances between the water molecule and the acceptor atom of 6-selenoguanine is about 0.6 Å longer for hydrogen bonds involving selenium atom. The relative Gibbs free energy of the 6-selenoguanine tautomers favors the selenone tautomer. The majority of the stable monohydrated complexes are the one in which the oxygen atom of water accepts the acidic N7-H proton while donating a proton to the carbonyl selenium atom of 6-selenoguanine; the interaction toward N7-H being stronger than that with the selenium site. The amino group planarity has been found to be increased in the hydrated complexes. The examination of molecular orbital reveals a moderate band gap between the donor and acceptor atoms of isolated and hydrated complexes. An excellent linear correlation is found to exist between electron density and laplacian of electron density with hydrogen-bond length through atoms in molecule analysis. The natural bond orbital analysis shows a maximum charge transfer of 0.060e for selenium acceptors and around 0.025e for selenium donors.