Hydrogen bonds with pi and sigma electrons as the multicenter proton acceptors: high level ab initio calculations

Pi-electrons of acetylene and sigma-electrons of molecular hydrogen were investigated as Lewis bases in different complexes. Hence high level ab initio calculations were performed up to the MP2/6-311++G(3df,3pd) level of approximation. It was found that species analyzed possess characteristics typical for H-bonded systems. The Bader theory was additionally applied; bond paths between proton and pi-electrons of acetylene or sigma-electrons of molecular hydrogen were detected with the corresponding bond critical points attributed to the proton-acceptor interactions. Numerous correlations between topological, geometrical and energetic parameters were also found. For example, the H...pi or H...sigma interaction is stronger for the shorter corresponding distance between the proton and the middle of C[triple bond]C or H-H bond. It is connected with the greater elongation of C[triple bond]C or H-H bonds and the greater transfer of electron charge from the Lewis base to the Lewis acid.     

How to determine whether intramolecular H...H interactions can be classified as dihydrogen bonds

Different types of intramolecular H...H interactions have been analyzed using the MP2/6-311++G(d,p) level of approximation. These are C-H...H-B, C-H...H-Al, C-H...H-C, C-H...H-O, O-H...H-Al and O-H...H-B contacts. Quantum theory of atoms in molecules and natural bond orbitals methods were applied to better understand the nature of these interactions. It was found that some of the species analyzed possess the characteristics of typical hydrogen bonds, such as the O-H...O ones. The electron charge transfer from the Lewis base to the antibonding X-H (for example O-H) orbital of the Lewis acid is one such characteristic. The NBO method may be considered decisive to classify any system as dihydrogen bonded.     

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.     

HArF和HNO二聚体相互作用氢键的理论研究

The hydrogen bonding interaction of 1:1 dimer formed between HNO and HArF molecule has been completely investigated in the present study using Second-order M?ller-Plesset Perturbation (MP2) method in conjunction with 6-311+G**, 6-311++G** and 6-311++G(2d,2p) basis sets. The standard and CP-corrected calculations have been employed to determine the equilibrium structures, the vibrational frequencies and interaction energies. The interaction energies of the dimers were also calculated at G2MP2 level. Two stable structures are found as the minima. Dimer I(H···F)is a five-membered cyclic hydrogen bonded structure and is more stable than the Dimer II(H···O). The blue-shifted N-H···F hydrogen bond is confirmed with standard and CP-corrected calculations by the MP2 and DFT methods in conjunction with different basis sets. The results obtained at MP2 in conjunction with different basis sets show there is a red-shifted hydrogen bond (Ar-H···O) in the Dimer II(H···O). The topological and electronic properties, the origin of red- and blue-shifted hydrogen bonds were investigated at MP2/6-311++G(2d,2p) with CP corrected calculations. From the NBO analysis, the reasonable explanations for the red- and blue-shifted hydrogen bonds were proposed.     

Red- and blue-shifted hydrogen bonds: the Bent rule from quantum theory of atoms in molecules perspective

MP2/6-311++G(3pd,3df) calculations were performed on complexes of acetylene and fluoroform acting as the proton donating systems and different Lewis bases being the proton acceptors since these complexes are linked through C-H···Y hydrogen bonds. Quantum Theory of Atoms in Molecules (QTAIM) is applied to explain the nature of these interactions. The characteristics of bond critical points are presented for these complexes. The inter-relations between energetic and geometrical parameters as well as the parameters derived from the Natural Bond Orbital (NBO) theory are analyzed here. Red- and blue-shifted hydrogen bonds are detected for the complexes investigated and the differences between those interactions are analyzed from the QTAIM perspective. It is shown that such differences are in agreement with the Bent rule. The position of the bond critical point of the proton donating bond is connected with the nature of hydrogen bonding, that is, if it is blue- or red-shifted.     

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.     

Properties of the halogen-hydride interaction: an ab initio and "atoms in molecules" analysis

X-Cl...H-Y interactions are analyzed by applying ab initio methods as well as the Bader theory. All calculations were performed using Pople's basis sets (6-311++G(2df,2pd) and 6-311++G(3df,3pd)) as well as the Dunning-type bases (aug-cc-pVDZ and aug-cc-pVTZ) within the MP2 method. For the complexes analyzed here, X-Cl and H-Y may be treated as a Lewis acid and a Lewis base, respectively. The Cl...H interactions are rather weak or at most moderate since, for the strongest interaction of the F3...HLi complex, the binding energy calculated at the MP2/6-311++G(3df,3pd) level of approximation amounts to -3.4 kcal/mol, and the H...Cl distance is equal to 2.65 A, less than the corresponding sum of van der Waals radii. These interactions may be classified as halogen-hydride interactions. However, some of the complexes analyzed, especially F3SiCl...HBeF and F3SiCl...HBeF, are very weakly bound, probably by typical van der Waals interactions.     

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

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

乙醛二聚体内蓝移型氢键的理论研究

利用理论方法研究了乙醛二聚体内的氢键. 在MP2/6-31＋G(d), B3LYP/6-31＋G(d), B3LYP/6-311＋＋G(d,p)和B3LYP/6-311＋＋G(3df,2p)水平上, 利用常规方法和均衡校正方法对3种稳定的乙醛二聚体进行了几何优化和振动频率计算. 计算结果表明: 在二聚体A和C中乙醛中C—H键强烈收缩, 存在显著的C—H…O蓝移型氢键. 自然键轨道(NBO)分析表明, 电子供体轨道和电子受体轨道之间相互作用的稳定化能、分子内电子密度重排、轨道再杂化和结构重组是决定氢键红移和蓝移的主要因素. 其中, 轨道间稳定化能属于键伸长效应, 分子内电子密度重排、轨道再杂化和电子受体内部结构重组属于键收缩效应. 在二聚体A和C中, 由于键收缩效应处于优势地位导致C—H…O蓝移氢键存在.     

乙醇-水分子团簇C2H5OH(H2O)n (n=1-9)稳定结构的量子化学研究

采用密度泛函理论B3LYP方法, 在B3LYP/6-311++G(2d,2p)//B3LYP/6-311++G(d,p)基组水平上对乙醇-水分子团簇(C₂H₅OH(H₂O)_n (n=1-9))的各种性质进行研究, 如: 优化的几何构型、结构参数、氢键、结合能、平均氢键强度、自然键轨道(NBO)电荷分布、团簇的生长规律等. 结果表明, 从二维(2-D)环状结构到三维(3-D)笼状结构的过渡出现在n=5的乙醇-水分子团簇中. 此外, 利用团簇结合能的二阶差分、形成能、能隙等性质, 发现在n=6时乙醇-水分子团簇的最低能量结构稳定性较好, 可能为幻数结构. 最后, 为了进一步探讨氢键本质, 将C₂H₅OH(H₂O)_n (n=2-9)最低能量结构的各种性质与纯水分子团簇(H₂O)_n (n=3-10)比较, 结果表明前者与后者中的水分子之间氢键相似.     

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.     

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.     

Weakly bound complexes of N2O: an ab initio theoretical analysis toward the design of N2O receptors

Ab initio calculations at MP2/6-311++G(2d,2p) and MP2/6-311++G(3df,3pd) computational levels have been used to analyze the interactions between nitrous oxide and a series of small and large molecules that act simultaneously as hydrogen bond donors and electron donors. The basis set superposition error (BSSE) and zero point energy (ZPE) corrected binding energies of small N2O complexes (H2O, NH3, HOOH, HOO*, HONH2, HCO2H, H2CO, HCONH2, H2CNH, HC(NH)NH2, SH2, H2CS, HCSOH, HCSNH2) vary between -0.93 and -2.90 kcal/mol at MP2/6-311++G(3df,3pd) level, and for eight large complexes of N2O they vary between -2.98 and -3.37 kcal/mol at the MP2/6-311++G(2d,2p) level. The most strongly bound among small N2O complexes (HCSNH2-N2O) contains a NH..N bond, along with S-->N interactions, and the most unstable (H2S-N2O) contains just S-->N interactions. The electron density properties have been analyzed within the atoms in molecules (AIM) methodology. Results of the present study open a window into the nature of the interactions between N2O with other molecular moieties and open the possibility to design N2O abiotic receptors.     

The role of chalcogen-chalcogen interactions in the intrinsic basicity acidity of beta-chalcogenovinyl(thio)aldehydes HC([double bond]X)[bond]CH[double bond]CH[bond]CYH

The intrinsic acidity and basicity of a series of beta-chalcogenovinyl(thio)aldehydes HC([double bond]X)[bond]CH[double bond]CH[bond]CYH (X=O, S; Y=Se, Te) were investigated by B3LYP/6-311+G(3df,2p) density functional and G2(MP2) calculations on geometries optimized at the B3LYP/6-31G(d) level for neutral molecules and at the B3LYP/6-31+G(d) level for anions. The results showed that selenovinylaldehyde and selenovinylthioaldehyde should behave as Se bases in the gas phase, because the most stable neutral conformer is stabilized by an X[bond]H...Se (X=O, S) intramolecular hydrogen bond (IHB). In contrast the Te-containing analogues behave as oxygen or sulfur bases, because the most stable conformer is stabilized by typical X...Y[bond]H chalcogen-chalcogen interactions. These compounds have a lower basicity than expected because either chalcogen-chalcogen interactions or IHBs become weaker upon protonation. Similarly, they are also weaker acids than expected because deprotonation results in a significantly destabilized anion. Loss of the proton from the X[bond]H or Y[bond]H groups is a much more favorable than from the C[bond]H groups. Therefore, for Se compounds the deprotonation process results in loss of the X[bond]H...Se (X=O, S) IHBs present in the most stable neutral conformer, while for Te-containing compounds the stabilizing X...Y[bond]H chalcogen-chalcogen interaction present in the most stable neutral conformer becomes repulsive in the corresponding anion.     

The possible covalent nature of N-H...O hydrogen bonds in formamide dimer and related systems: an ab initio study

The N-H...O hydrogen bonds are analyzed for formamide dimer and its simple fluorine derivatives representing a wide spectrum of more or less covalent interactions. The calculations were performed at the MP2/6-311++G(d,p) level of approximation. To explain the nature of such interactions, the Bader theory was also applied, and the characteristics of the bond critical points (BCPs) were analyzed: the electron density at BCP and its Laplacian, the electron energy density at BCP and its components, the potential electron energy density, and the kinetic electron energy density. These parameters are used to justify the statement that some of the interactions analyzed are partly covalent in nature. An analysis of the interaction energy components for the systems considered indicates that the covalent character of the hydrogen bond is manifested by a markedly increased contribution of the delocalization term relative to the electrostatic interaction energy. Moreover, the ratio of stabilizing the delocalization/electrostatic contributions grows linearly with the decreasing lengths of the hydrogen bond.     

Theoretical investigation of the weakly dihydrogen bonded complexes FArCCH...HBeX (X = H, F, Cl, Br)

An ab initio computational study of the properties of four linear dihydrogen-bonded complexes formed between the first compound with an Ar-C chemical bond (FArCCH) and HBeX (X = H, F, Cl, and Br) molecules was undertaken at the MP2/6-311++G(2d,2p) level of theory. The calculated complexation energy at MP2 and G2(MP2) levels decreases in the order HBeH...HCCArF > BrBeH...HCCArF > ClBeH...HCCArF > FBeH...HCCArF. The intermolecular stretching frequency, and shifts within the monomers, are compared with the energetic strength of complexation.     

N, P, and As ylides and aza- and arsa-Wittig reactions from topological analyses of electron density

The nature of bonding between N, P, and As constituent atoms in ylide systems with the R(3)XYR' formula (X = N, P, As; Y = N, P, As; R = F, H; R' = H, CH(3)) has been characterized by ab initio (MP2/6-311++G**) and density functional theory (B3LYP/6-311++G**) calculations. Its electronic structure has been analyzed through electron density with the quantum theory of atoms in molecules and the electron localization function (ELF). The characteristics of the central bond are inspected with the calculated rotational barriers. The results show that N has a behavior different from that of the remaining pnicogen atoms (P, As), where the bond is much stronger. Fluorine substituents strengthen the X-Y bond, reduce the bond distance, and increase the electron density in the central bond so that the substituent pulls charge from the bond in the pnicogen X atom. For the N-pnicogen ylides, the results showed different bonding characters between F and X atoms; depending on the position of the F atom, the difference of the bond character is sensed by the basin synaptic order, as it is deduced from the analysis of the ELF basins. The energy profiles of the rotational barriers have been calculated at the MP2/6-311++G** level, indicating that the electronegativity of the substituents is a relevant factor that has consequences in the characteristics of the X-Y bond.     

Ab initio study of hydrogen-bond formation between cyclic ethers and selected amino acid side chains

Binding energies for hydrogen-bonded complexes of six cyclic ethers with five hydrogen-bond donor molecules that mimic selected amino acid side chains have been calculated at the MP2/6-31G*, MP2/6-31+G*, MP2/6-311++G**(single point), and MP2/aug-cc-pvtz levels, using geometries obtained with or without counterpoise corrections throughout the geometry optimization. The calculated basis set superposition error (BSSE) amounts to 10-20% and 5-10% of the uncorrected binding energies for the neutral and ionic species, respectively, at the MP2/aug-cc-pvtz level. The authors conclude that the O...H distances in the hydrogen bonds and binding energies for the studied systems may be determined with uncertainties of up to 0.08 A and 1-2 kcal/mol, respectively, in comparison with the MP2/aug-cc-pvtz values at a reasonable computational cost by performing standard geometry optimization at the MP2/6-31+G* level. Hydrogen-bond formation energies are more negative for cyclic ethers compared to their counterparts with a C=C double bond in the ring next to the oxygen atom. The less negative hydrogen-bonding energy and the increased O...H separation have been attributed to the reduced basicity of the ether oxygen when the lone pairs can enter conjugation with the pi-electrons of the Calpha=Cbeta double bond. The present study is the first step toward the development of an affordable computational level for estimating the binding energies of natural product, fused ring ether systems to the human estrogen receptor.