气相中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)分析了三种复合物中氢键的电子密度拓扑性质.     

气相中开壳型(CH3)2S(O)…HOO红移氢键复合物的结构与性质

在DFT-B3LYP/6-311++G**水平上分别求得(CH3)2S…HOO和(CH3)2O…HOO开壳型氢键复合物势能面上的稳定构型. 频率分析表明, 与单体HOO自由基相比, 复合物中H10-O11键伸缩振动频率发生显著的红移, 红移值分别为424.21和374.22 cm-1. 在MP2/6-311++G**水平计算得到, 含基组重叠误差(BSSE)校正和零点振动能(ZPVE)校正的相互作用能分别为-24.68和-31.01 kJ·mol-1. 自然键轨道(NBO)理论分析表明, 在(CH3)2S…HOO复合物中, 引起H10-O11键变长的因素包括两种电荷转移: (1) LP(S1)1→σ*(H10-O11); (2) LP(S1)2→σ*(H10-O11), 其中LP(S1)2→σ*(H10-O11)占主要作用, 总的结果是使σ*(H10-O11)的自然布居数增加了37.27 me; 在(CH3)2O…HOO中也有相似的电荷转移的超共轭作用. AIM理论分析表明, S1…H10间和O1…H10间都存在键鞍点, ▽2ρ(r)分别为0.06196和0.03745, 说明这种相互作用介于共价键和离子键之间, 偏于静电作用.     

CH3SH与HOCl分子间氢键和卤键的结构与性质

在DFT-B3LYP及MP2/6-311＋＋G**水平上分别求得CH3SH…HOCl氢键复合物和CH3SH…ClOH卤键复合物势能面上的稳定构型. 频率分析表明, 与单体HOCl相比, 在两种复合物中, Cl(9)—O(7)和H(8)—O(7)键伸缩振动频率发生显著的红移. 经MP2/6-311＋＋G**水平计算的含基组重叠误差(BSSE)校正的气相中相互作用能分别为－19.23和－6.85 kJ•mol－1. 自然键轨道理论(NBO)分析表明, 在CH3SH…ClOH卤键复合物中, 引起Cl(9)—O(7)键变长的因素包括2种电荷转移: (i)孤对电子LP[S(1)]1→σ*[Cl(9)—O(7)]; (ii)孤对电子LP[S(1)]2→σ*[Cl(9)—O(7)], 其中孤对电子LP[S(1)]2→σ*[Cl(9)—O(7)]转移占主要作用, 总的结果是使σ*[Cl(10)—O(11)]的自然布居数增加, 同时O(7)和Cl(9)原子s成分均增加的杂化重优具有与电荷转移作用相同的“拉长效应”; 在CH3SH…HOCl氢键复合物中也存在类似的电荷转移, 但是O(7)原子的再杂化效应不同于前者. 自然键共振理论(NRT)进行键序分析表明, 在氢键复合物和卤键复合物中, H(8)—O(7)和Cl(9)—O(7)键的键序都减小. 通过分子中原子理论(AIM)分析了复合物中氢键和卤键的电子密度拓扑性质.     

CH_3SH和CH_3SLi分子间的锂键和氢键共存型相互作用

在CH3SLi+CH3SH势能面上求得锂键和氢键共存型复合物的两种稳定构型.频率分析表明,与单体相比复合物中S(5)—Li(6)键伸缩振动频率发生红移,而C(8)—H(10)键伸缩振动频率发生蓝移.经B3LYP/6-311++G**,MP2/6-311++G**及MP2/AUG-CC-PVDZ水平计算的含基组重叠误差(BSSE)校正的复合物Ⅰ中相互作用能分别为-58.99,-57.87和-62.89kJ·mol-1.采用自然键轨道(NBO)理论,分析了复合物中单体轨道间的电荷转移,电子密度重排及其与相关键键长变化的本质等.采用分子中的原子(AIM)理论分析了复合物中氢键和锂键的电子密度拓扑性质.在极化连续模型(PCM)下,考察了溶剂化效应.结果表明,所考察的水、二甲亚砜、乙醇和乙醚等四种溶剂均使单体间的相互作用能增大,且溶剂对复合物中的锂键结构及其振动频率具有显著的影响,而对复合物中的氢键的振动频率影响不大.     

GeH4与HX(X＝F,Cl,Br)间二氢键复合物的理论研究

对GeH4与HX形成的二氢键复合物的结构特征及本质进行了探讨.在MP2/6-311 ++G(3 d,3p)水平优化、频率验证得到复合物的分子结构,通过分子的几何参数及电子密度拓扑分析,确认GeH4与卤化氢已形成了二氢键复合物.MP2/6-311 ++ G(3d,3p)水平下进行BSSE校正后的结合能为3.281到4.5...     

(CH3)2S与HOCI分子间的卤键和氢键相互作用

在DFT-B3LYP/6-311++G**水平上分别求得(CH3)2S…ClOH卤键复合物和(CH3)2S…HOCl氢键复合物势能面上的稳定构型.频率分析表明,与单体HOC1相比,在两种复合物中,10C1-11O和12H-11O键伸缩振动频率发生显著的红移.经MP2/6-311++G**水平计算的含基组重叠误差(BSSE)校正的气相中相互作用能分别为-11.69和-24.16 kJ·mol-1.自然键轨道理论(NBO)分析表明,在(CH3)2s…ClOH卤键复合物中,引起10Cl-11O键变长的因素包括两种电荷转移:(i)孤对电子LP(1S)1→σ*(10C1-11O);(ⅱ)孤对电子LP(1s)2→σ*(10Cl-11O),其中孤对电子LP(lS)2→σ*(10Cl-11O)转移占主要作用,总的结果是使σ(10Cl-11O)的自然布居数增加0.14035e,同时11O原子的再杂化使其与10Cl成键时s成分增加,即具有与电荷转移作用同样的"拉长效应";在(CH3)2s…HOCl氢键复合物中也存在类似的电荷转移,但是11O原子的再杂化不同于前者.自然键共振理论(NRT)进行键序分析表明,在卤键复合物和氢键复合物中,10Cl-11O和12H-11O键的键序都减小.通过分子中原子理论(AIM)分析了复合物中卤键和氢键的电子密度拓扑性质.     

（CH3）2S与HOCl分子间的卤键和氢键相互作用

在DFT-B3LYP/6-311＋＋G＊＊水平上分别求得（CH3）2S…ClOH卤键复合物和（CH3）2S…HOCl氢键复合物势能面上的稳定构型.频率分析表明,与单体HOCl相比,在两种复合物中,10Cl—11O和12H—11O键伸缩振动频率发生显著的红移.经MP2/6-311＋＋G＊＊水平计算的含基组重叠误差（BSSE）校正的气相中相互作用能分别为-11.69和-24.16kJ·mol^-1.自然键轨道理论（NBO）分析表明,在（CH3）2S…ClOH卤键复合物中,引起10Cl—11O键变长的因素包括两种电荷转移：（i）孤对电子LP（1S）1→σ＊（10Cl—11O）;（ii）孤对电子LP（1S）2→σ＊（10Cl—11O）,其中孤对电子LP（1S）2→σ＊（10Cl—11O）转移占主要作用,总的结果是使σ＊（10Cl—11O）的自然布居数增加0.14035e,同时11O原子的再杂化使其与10Cl成键时s成分增加,即具有与电荷转移作用同样的“拉长效应”;在（CH3）2S…HOCl氢键复合物中也存在类似的电荷转移,但是11O原子的再杂化不同于前者.自然键共振理论（NRT）进行键序分析表明,在卤键复合物和氢键复合物中,10Cl—11O和12H—11O键的键序都减小.通过分子中原子理论（AIM）分析了复合物中卤键和氢键的电子密度拓扑性质.     

(CH3)2S与HOCl分子间的卤键和氢键相互作用

在DFT-B3LYP/6-311++G**水平上分别求得(CH3)2S…ClOH卤键复合物和(CH3)2S…HOCl氢键复合物势能面上的稳定构型. 频率分析表明, 与单体HOCl相比, 在两种复合物中, 10Cl—11O和12H—11O键伸缩振动频率发生显著的红移. 经MP2/6-311++G**水平计算的含基组重叠误差(BSSE)校正的气相中相互作用能分别为-11.69和-24.16 kJ·mol-1. 自然键轨道理论(NBO)分析表明, 在(CH3)2S…ClOH卤键复合物中, 引起10Cl—11O键变长的因素包括两种电荷转移: (i) 孤对电子LP(1S)1→σ*(10Cl—11O); (ii) 孤对电子LP(1S)2→σ*(10Cl—11O), 其中孤对电子LP(1S)2→σ*(10Cl—11O)转移占主要作用, 总的结果是使σ*(10Cl—11O)的自然布居数增加0.14035e, 同时11O原子的再杂化使其与10Cl成键时s成分增加, 即具有与电荷转移作用同样的“拉长效应”; 在(CH3)2S…HOCl氢键复合物中也存在类似的电荷转移, 但是11O原子的再杂化不同于前者. 自然键共振理论(NRT)进行键序分析表明, 在卤键复合物和氢键复合物中, 10Cl—11O和12H—11O键的键序都减小. 通过分子中原子理论(AIM)分析了复合物中卤键和氢键的电子密度拓扑性质.     

CH3SH和CH3SLi分子间的锂键和氢键共存型相互作用

在CH₃SLi＋CH₃SH势能面上求得锂键和氢键共存型复合物的两种稳定构型. 频率分析表明, 与单体相比复合物中S(5)—Li(6)键伸缩振动频率发生红移, 而C(8)—H(10)键伸缩振动频率发生蓝移. 经B3LYP/6-311＋＋G**, MP2/6-311＋＋G**及MP2/AUG-CC-PVDZ水平计算的含基组重叠误差(BSSE)校正的复合物?中相互作用能分别为－58.99, －57.87和－62.89 kJ•mol^－1. 采用自然键轨道(NBO)理论, 分析了复合物中单体轨道间的电荷转移, 电子密度重排及其与相关键键长变化的本质等. 采用分子中的原子(AIM)理论分析了复合物中氢键和锂键的电子密度拓扑性质．在极化连续模型(PCM)下, 考察了溶剂化效应. 结果表明, 所考察的水、二甲亚砜、乙醇和乙醚等四种溶剂均使单体间的相互作用能增大, 且溶剂对复合物中的锂键结构及其振动频率具有显著的影响, 而对复合物中的氢键的振动频率影响不大.     

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)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.     

Does the molecular structure of CaH2 affect the dihydrogen bonding in CaH2 … HY (Y = CH3, C2H3, C2H, CN, and NC) complexes? A quantum chemistry study using MP2 and B3LYP methods

Second-order M ller-Plesset(MP2) and density functional theory(DFT) calculations have been carried out in order to investigate the structures and properties of dihydrogen-bonded CaH 2 HY(Y = CH 3,C 2 H 3,C 2 H,CN,and NC) complexes.Our calculations revealed two possible structures for CaH 2 in CaH 2 HY complexes:linear(I) and bent(II).The bond lengths,interaction energies,and strengths for H H interactions obtained by both MP2 and B3LYP methods are quite close to each other.It was found that the interaction energy decreases with increasing electron density at the Ca-H bond critical point.Atom-in-molecule(AIM) results show that for all of Ca-H H-Y interactions considered here,the Laplacian of the electron density at the H H bond critical point is positive,indicating the electrostatic nature of these Ca-H H-Y dihydrogen bonded systems.     

Theoretical investigation of C–HM interactions in organometallic complexes: A natural bond orbital (NBO) study

The nature of C–HM agostic interactions in model metal complexes [M²⁺(CH₂CH₃)(PH₃)_nCl] (where M = Sc, Ti, V, Ti, Cr, Mn, Fe, Co, Ni, Cu, Zn; n = 1, 2, 3, 4) was studied with the natural bond orbital analysis (NBO) approach using density functional theory (DFT) optimized geometries at the B3LYP/6-31G(d,p) level of theory. The effect of nature of metal, coordination number, oxidation state and ligand field effects on the agostic interaction is examined. A set of 20 crystal structures of organometallic complexes taken from the Cambridge Structural Database (CSD) was studied computationally employing AIM theory and NBO analysis, and the applicability of these methods was critically accessed in demarcating the two types of interaction.     

Density Functional Theory Study of the Interaction between Guanine and Catechin

The interacting patterns and mechanism of the catechin and guanine have been investigated with the density functional theory B3LYP method by 6‐31G* basis set. Fourteen stable structures for the catechin‐guanine complexes have been found which form two hydrogen bonds at least. The results indicate that the complexes are mainly stabilized by the hydrogen bonding interactions. At the same time, the number and strength of hydrogen bond play a co‐determinant parts in the stability of the complexes which can form two or more hydrogen bonds. Theories of atoms in molecules (AIM) and natural bond orbital (NBO) have been adopted to investigate the hydrogen bonds involved in all systems. The interaction energies of all complexes have been corrected for basis set superposition error (BSSE), ranging from ?38.86 to ?14.56 kJ/mol. The results showed that the hydrogen bonding contributes to the interaction energies dominantly. The corresponding bonds stretching motions in all complexes are red‐shifted relative to that of the monomer, which is in agreement with experimental results.     

Intramolecular hydrogen bond in the hydroxycyclohexadienyl peroxy radicals

The hydroxycyclohexadienyl peroxy radicals (HO? C₆H₆? O₂) produced from the reaction of OH‐benzene adduct with O₂ were studied with density functional theory (DFT) calculations to determine their characteristics. The optimized geometries, vibrational frequencies, and total energies of 2‐hydroxycyclohexadienyl peroxy radical II_s and 4‐hydroxycyclohexadienyl peroxy radical III_s were calculated at the following theoretical levels, B3LYP/6‐31G(d), B3LYP/6‐311G(d,p), and B3LYP/6‐311+G(d,p). Both were shown to contain a red‐shifted intramolecular hydrogen bond (O? H … O? H bond). According to atoms‐in‐molecules (AIM) analysis, the intramolecular hydrogen bond in the 2‐hydroxycyclohexadienyl peroxy radical II_s is stronger than that one in 4‐hydroxycyclohexadienyl peroxy radical III_s, and the former is the most stable conformation among its isomers. Generally speaking, hydrogen bonding in these radicals plays an important role to make them more stable. Based on natural bond orbital (NBO) analysis, the stabilization energy between orbitals is the main factor to produce red‐shifted intramolecular hydrogen bond within these peroxy radicals. The hyperconjugative interactions can promote the transfer of some electron density to the O? H antibonding orbital, while the increased electron density in the O? H antibonding orbital leads to the elongation of the O? H bond and the red shift of the O? H stretching frequency. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

The competition between the intramolecular hydrogen bond and π‐electron delocalization in trifluoroacetylacetone—A theoretical study

Conformational study of trifluoroacetylacetone was carried out using the HF, B3LYP, and MP2 methods with the 6‐31G(d, p) and 6‐311++G(d, p) basis sets. All of the results show that the chelated enol structures (E11 and E31) have extra stability with respect to the other forms and one of them (E11) is global minimum. The energy gap between the chelated forms is in the range 0.7–5.9 kJ mol^?1. Theoretical calculations show that this compound has an asymmetric double minimum potential energy surface which is in contrast with the electron diffraction result. Moreover, the computational results predict that due to the withdrawing effect of CF₃ group, hydrogen bond in trifluoroacetylacetone is weaker than the acetylacetone. Because of the more stability of E11, it is expected that the hydrogen bond energy in E11 is greater than the E31, but at all of the computational levels with most extended basis set the converse results were observed. These results clearly show that the hydrogen bond is not a superior parameter in conformational preference and the contribution of resonance is probably greater than the hydrogen bond. Finally, the analysis of this system by quantum theory of atoms in molecules and natural bond orbital methods fairly support the ab initio results. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Structures and Electronic Properties of HOCl···HCOCl Complexes

B3LYP/6-311++G** and MP2/6-311++G** calculations were used to analyze the interaction between hypochlorous acid (HOCl) and formyl chloride (HCOCl). The results showed that there were four equilibrium geometries (S1, S2, S3, and S4) optimized at B3LYP/6-311++G** level, and all the equilibrium geometries were confirmed to be in stable states by analytical frequency calculations. Complexes S1 and S3 use the 5H atom of HOCl as proton donor and the terminal 1O atom of HCOCl as acceptor to form red shift hydrogen bond systems. However, the blue-shifted hydrogen bond (2C-3H···6O) coexists with 4Cl···5O interaction in structures S2. As for S4, it uses the 7Cl atom of HOCl as proton donor and the terminal 1O atom of HCOCl as acceptor to form red shift halogen bond system. Interaction energies between monomers in the four complexes corrected with basis set superposition error (BSSE) and zero-point vibrational energy (ZPVE) lie in the range from −5.05 to −14.76 kJ·mol⁻¹ at MP2/6-311++G** level. The natural bond orbital (NBO) and atoms in molecules (AIM) theories have also been applied to explain the structures and the properties of the complexes.     

Hydrogen bonds and ionic interactions in Guanidine/Guanidinium complexes: a computational case study

It is frequently said that hydrogen bonds (HBs) are enhanced by ionic interactions and in this article we intend to determine the degree at which this reinforcement happens. Considering our interest in the Guanidine(neutral)/Guanidinium(cation) system and its particular nature, all the possible 1:1 complexes with the Chloride(anion)/Hydrochloric acid(neutral) system have been studied at different levels of computation (B3LYP with 6-31+G* and TZVP basis sets; MP2 with 6-31+G*, 6-311++G** and aug-cc-pVDZ basis sets; CBS-QB3 and G3MP2). The nature of these interactions established in all the systems and, when possible, at all the levels of computation used in this study, has been analyzed using Atoms in Molecules and Natural Bond Orbital methodologies. By examining the interaction energy, the electron density at the bond critical bonds, the atomic energy, the charge transfer, the orbital energy, and the deformation energy we can conclude that HBs are stronger when the ionic interaction is stronger. Thus, both interactions do not work in an independent manner but one reinforces the other to different degrees depending on the nature of the charges present. Several correlations with the interaction energy have been found and a partition of the contributions of both the HB and ionic forces to the total interactions is proposed.     

Comparison of triel bonds with different chalcogen electron donors: Its dependence on triel donor and methyl substitution

The complexes between R₃Tr (Tr = B, Al, and Ga; R = H, F, Cl, and Br) and H₂X (X = O, S, and Se) were theoretically studied. The interaction energies of R₃Al⋯H₂X and R₃Ga⋯H₂X are consistent with the electronegativity of the halogen atom R (R ≠ H), but an opposite dependence is found for R₃B⋯H₂X. The triel bond of R₃Tr⋯H₂X is weaker for the heavier chalcogen donor. The dependence of triel bonding strength on the triel atom is complicated, depending on the nature of R and X. The methyl substitution of H₂X causes a substantial increase in the interaction energy from −5.74 kcal/mol to −22.88 kcal/mol, and its effect is relevant to the nature of Tr, X, and R groups. For the S and Se donors, the increased percentage of interaction energy is almost the same due to the methyl substitution, which is larger than that of the O analogue. In most triel-bonded complexes, electrostatic dominates and polarization has comparable contribution. However, polarization plays a dominant role in R₃B⋯ and R₃B⋯ (R = Cl and Br; R′ = H and Me).