Density Functional Theory Study on Hydrogen Bonding Interaction of Catechin-(H2O)n

Density functional theory B3LYP method with 6‐31G* basis set has been used to optimize the geometries of the catechin, water and catechin‐(H₂O)_n complexes. The vibrational frequencies have been studied at the same level to analyze these complexes. Six and eleven stable structures for the catechin‐H₂O and catechin‐(H₂O)₂ have been found, respectively. Theories of atoms in molecules (AIM) and natural bond orbital (NBO) have been utilized to investigate the hydrogen bonds involved in all the systems. The interaction energies of all the complexes corrected by basis set superposition error, are from ?13.27 to ?83.56 kJ/mol. All calculations also indicate that there are strong hydrogen‐bonding interactions in catechin‐water complexes. The strong hydrogen‐bonding contributes to the interaction energies dominantly. The O–H stretching motions in all the complexes are red‐shifted relative to that of the monomer.     

Density Functional Theory Study on Interaction between Catechin and Thymine

The interacting patterns and mechanism of the catechin and thymine have been investigated with the density functional theory Becke's three-parameter nonlocal exchange functional and the Lee, Yang, and Parr nonlocal correlation functional (B3LYP) method by 6-31+G*basis set. Thirteen stable structures for the catechin-thymine complexes have been found which form two hydrogen bonds at least. The vibrational frequencies are also studied at the same level to analyze these complexes. The results indicated that catechin interactedwith thymine by three different hydrogen bonds as N-H…O、C-H…O、O-H…O and the complexes are mainly stabilized by the hydrogen bonding interactions. Theories of atoms in molecules and natural bond orbital have been adopted to investigate the hydrogen bondsinvolved in all systems. The interaction energies of all complexes have been corrected for basis set superposition error, which are from -18.15 kJ/mol to -32.99 kJ/mol. The results showed that the hydrogen bonding contribute 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.     

Density Functional Theory Study of the Interaction between Thymine and Luteolin

The density function B3LYP method has been used to optimize the geometries of the luteolin, thymine and luteolin‐thymine complexes at 6‐31+G?? basis. The vibrational frequencies have been studied at the same level to analyze these seventeen complexes, respectively. Theories of atoms in molecules (AIM) and natural bond orbital (NBO) have been utilized to investigate the hydrogen bonds involved in all the systems. The interaction energies of the complexes corrected by basis set superposition error are between ?93.00–?76.69 kJ/mol. The calculating results indicate that strong hydrogen bonding interactions have been found in the luteolin‐thymine complexes.     

儿茶素-胞嘧啶分子间相互作用的密度泛函研究

采用密度泛函理论的B3LYP方法,在6-31+G*基组水平上研究了儿茶素-胞嘧啶分子间相互作用机制,得到稳定的儿茶素-胞嘧啶复合物11个.计算结果表明氢键对于复合物的稳定性起着重要的作用,并且当复合物形成2个或更多的氢键时,氢键的类型及强度共同决定着复合物的稳定性.我们还应用了分子中的原子(AIM)理论和自然键轨道(NBO)理论对这11种复合物中氢键的性质和特征进行了分析.通过研究发现,所有的氢键复合物进行基组重叠误差(BSSE)校正后的相互作用能为-17.35～-43.27kJ/mol,相互作用能主要由氢键所贡献.振动分析显示,氢键的形成使相对应键的对称伸缩振动频率减小,说明这些复合物中形成的氢键都是正常的红移型氢键,与实验结果相一致.     

Density functional theory study on the interaction of catechin and cytosine

The interacting patterns and mechanism of the catechin and cytosine have been investigated using the density functional theory B3LYP method with 6-31+G* basis set.Eleven stable structures of the catechin-cytosine complexes have been found respectively.The results indicate that the complexes are mainly stabilized by the hydrogen bonding interactions.Theories of atoms in molecules(AIM) and natural bond orbital(NBO) have been utilized to investigate the hydrogen bonds involved in all the systems.The interactio...     

A polarizable dipole–dipole interaction model for evaluation of the interaction energies for NH···OC and CH···OC hydrogen‐bonded complexes

In this article, a polarizable dipole–dipole interaction model is established to estimate the equilibrium hydrogen bond distances and the interaction energies for hydrogen‐bonded complexes containing peptide amides and nucleic acid bases. We regard the chemical bonds N? H, C?O, and C? H as bond dipoles. The magnitude of the bond dipole moment varies according to its environment. We apply this polarizable dipole–dipole interaction model to a series of hydrogen‐bonded complexes containing the N? H···O?C and C? H···O?C hydrogen bonds, such as simple amide‐amide dimers, base‐base dimers, peptide‐base dimers, and β‐sheet models. We find that a simple two‐term function, only containing the permanent dipole–dipole interactions and the van der Waals interactions, can produce the equilibrium hydrogen bond distances compared favorably with those produced by the MP2/6‐31G(d) method, whereas the high‐quality counterpoise‐corrected (CP‐corrected) MP2/aug‐cc‐pVTZ interaction energies for the hydrogen‐bonded complexes can be well‐reproduced by a four‐term function which involves the permanent dipole–dipole interactions, the van der Waals interactions, the polarization contributions, and a corrected term. Based on the calculation results obtained from this polarizable dipole–dipole interaction model, the natures of the hydrogen bonding interactions in these hydrogen‐bonded complexes are further discussed. © 2013 Wiley Periodicals, Inc.     

甲醛与甲酰胺相互作用的从头算研究

The optimizations geometries and vibrational frequencies of H2CO,HCONH2 and acquired 3 complexes between H2CO?HCONH2 have been calculated by using the ab initio method at the MP2/6-31G（ d）and MP2 （FC）/6-311++G（d,p）level. The non-minimum structures with negative vibrational frequencies are excluded. The lowest energy conformer of these complexes is a cyclic structure with N - H?O and C - H?O hydrogen bonds in a common plane. No significant changes are observed in the geometries of the monomers in their complexed state. The most characteristic geometrical properties of the complex are the lengthening of the contacting N-H bonds by 0.4-1.1 pm,and the general shortening of the contacting C-H bonds by 0.3-0.4 pm with respect to the monomers. The interaction energies of complexes have been corrected by the basis set superposition error （BSSE）using the full Boys-Bernardi counterpoise correction scheme. The corrected complex interaction energies of 3 structures at MP2/6-311++G（2df,3p）/ / MP2（FC）/6-311++G（d,p）level are -29.94, -16.10 and -18.45 kJ/mol,respectively. The interaction energies indicate that C - H?O is a weak hydrogen bond. The results of natural bond orbital population analysis reveals that there is only a small charge-transfer in the process of forming the complexes. The results of natural bond orbital analysis and atom in the molecule scheme appear quite significant in view of their importance for understanding the mechanisms of intermolecular interaction leading to hydrogen bonding. The results of molecular interaction energy decomposition analysis show that the electrostatic interaction plays an essential role in stabilizing the H2CO?HCONH2 complexes.     

Red and blue shifted hydridic bonds

By performing MP2/aug‐cc‐pVTZ ab initio calculations for a large set of dimer systems possessing a R? H hydridic bond involved in diverse types of intermolecular interactions (dihydrogen bonds, hydride halogen bonds, hydride hydrogen bonds, and charge‐assisted hydride hydrogen bonds), we show that this is rather an elongation than a shortening that a hydride bond undergoes on interaction. Contrary to what might have been expected on the basis of studies in uniform electric field, this elongation is accompanied by a blue instead of red shift of the R? H stretching vibration frequency. We propose that the “additional” elongation of the R? H hydridic bond results from the significant charge outflow from the sigma bonding orbital of R? H that weakens this bond. The more standard red shift obtained for stronger complexes is explained by means of the Hermansson's formula and the particularly strong electric field produced by the H‐acceptor molecule. © 2014 Wiley Periodicals, Inc.     

A Comparative Study of Hydrogen Bonding Using Density Functional Theory

The hydrogen‐bond energies of water dimer and water‐formaldehyde complexes have been studied using density functional theory (DFT). Basis sets up to aug‐cc‐pVXZ (X=D, T, Q) were used. It was found that counterpoise corrected binding energies using the aug‐cc‐pVDZ basis set are very close to those predicted with the aug‐cc‐pVQZ set. Comparative studies using various DFT functionals on these two systems show that results from B3LYP, mPW1PW91 and PW91PW91 functionals are in better agreements with those predicted using high‐level ab initio methods. These functionals were applied to the study of hydrogen bonding between guanine (G) and cytosine (C), and between adenine (A) and thy mine (T) base pairs. With the aug‐cc‐pVDZ basis set, the predicted binding energies of base pairs are in good agreement with the most elaborate ab initio results.     

Density functional theory study of the interaction between formamide and uracil

The hydrogen bonding of complexes formed between the formamide and uracil molecule has been fully investigated in the present study using the density functional theory (DFT) method at varied basis set levels from 6‐31G to 6‐311++G(d,p). Eight stable cyclic structures with two hydrogen bonds involved in the interaction are found on the potential energy surface (PES). The four structures are seven‐membered rings; the others are eight‐membered rings. The eight‐membered ring is preferred over the seven‐membered one by analyzing the hydrogen bond lengths and the interaction energies. The infrared (IR) spectrum frequencies, IR intensities, and the vibrational frequency shifts are reported. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

木犀草素与胞嘧啶相互作用的理论研究

采用密度泛函理论中的B3LYP方法,在6-31＋G＊基组水平上对木犀草素、胞嘧啶、木犀草素-胞嘧啶复合物进行结构优化和振动频率分析,得到了12种稳定复合物.并应用分子中的原子理论（AIM）分析、自然键轨道（NBO）理论分析得到复合物氢键性质和特征.通过基组重叠误差（BSSE）校正后的相互作用能、成键临界点电荷密度、二阶...     

The effect of methylation on the hydrogen-bonding and stacking interaction of nucleic acid bases

Methylated nucleosides play an important role in DNA/RNA function, and may affect republication by interrupting the base-pairing and base-stacking. In order to investigate the effect of methylation on the interaction between nucleic acid bases, this work presents the hydrogen-bonding and stacking interactions between 5-methylcytosine and guanine (G), cytosine (C) and G, 1-methyladenine and thymine (T), as well as adenine and T. Geometry optimization and potential energy surface scan have been performed for the involved complexes by MP2 calculations. The interaction energies, which were corrected for the basis-set superposition error by the full Boys–Bernardi counterpoise correction scheme, were used to evaluate the interaction intensity of these nucleic acid bases. The atoms in molecules theory and natural bond orbital analysis have been performed to study the hydrogen bonds in these complexes. The result shows that the methyl substitute contributes the stability to these complexes because it enhances either the hydrogen bonding or the staking interaction between nucleic acid bases studied.     

带电组氨酸侧链与DNA碱基间非键作用强度的理论研究

采用MP2方法和6-31+G(d,p)基组优化得到了带有一个正电荷的组氨酸侧链与4个DNA碱基间形成的18个氢键复合物的气相稳定结构, 从文献中获取了组氨酸侧链与DNA碱基间形成的12个堆积和T型复合物的气相稳定结构, 使用包含基组重叠误差(BSSE)校正的MP2方法和aug-cc-pVTZ基组及密度泛函理论M06-2X-D3方法和aug-cc-pVDZ基组计算了这些复合物的结合能. 研究结果表明, 包含BSSE校正的M06-2X-D3方法和aug-cc-pVDZ基组能够给出较准确的结合能; 气相条件下, 组氨酸侧链与同种DNA碱基间的离子氢键作用明显强于堆积作用和T型作用, 组氨酸侧链最易通过离子氢键与胞嘧啶C和鸟嘌呤G作用形成氢键复合物, 组氨酸与胞嘧啶C和鸟嘌呤G间的T型作用强于与腺嘌呤A和胸腺嘧啶T间的离子氢键作用; 水相条件下, 组氨酸侧链与同种DNA碱基间的离子氢键作用仍明显强于堆积作用和T型作用, 组氨酸侧链更易与胞嘧啶C和鸟嘌呤G相互作用形成氢键复合物, 但是最强的组氨酸侧链与胞嘧啶C间的T型作用明显弱于与腺嘌呤A和胸腺嘧啶T间的离子氢键作用, 说明水相条件下组氨酸侧链与DNA碱基间主要通过离子氢键作用形成氢键复合物.     

Hydrogen‐bonded complexes of nicotine with simple alcohols

Ab initio and density functional theory studies have been performed on the hydrogen‐bonded complexes of neutral and protonated nicotine with ethanol, methanol, and trifluromethanol to explore their relative stability in a systematic way. Among all the hydrogen‐bonded nicotine complexes considered here, protonated forms in nicotine–ethanol and nicotine–methanol, and neutral form in nicotine–trifluromethanol complexes have been found to be the most stable. In the former two complexes, the proton attached to the pyrrolidine nitrogen acts as a strong hydrogen bond donor, whereas the pyrrolidine nitrogen atom acts as a hydrogen bond acceptor in the latter case. Neutral complex of nicotine with trifluromethanol has been found to possess a very short hydrogen bond (1.57 Å) and basis set superposition error corrected hydrogen bond energy value of 19 kcal/mol. The nature of the various hydrogen bonds formed has been investigated through topological aspects using Bader's atoms in molecules theory. From the calculated topological results, excellent linear correlation is shown to exist among the hydrogen bond length, electron density, and its Laplacian at the bond critical points for all the complexes considered. The natural bond orbital analysis has been carried out to investigate the charge transfer in the nicotine alcohol complexes. In contrast to the blue shifting behavior that is generally exhibited by other C? H···O hydrogen bonds involving sp³ carbon atom, the C? H···O hydrogen bond in the protonated nicotine–ethanol and methanol complexes has been found to be proper with red shifting in nature. © 2011 Wiley Periodicals, Inc.     

Theoretical study on the interactions of halogen‐bonds and pnicogen‐bonds in phosphine derivatives with Br2, BrCl,and BrF

The MP2 ab initio quantum chemistry methods were utilized to study the halogen‐bond and pnicogen‐bond system formed between PH₂X (X = Br, CH₃, OH, CN, NO₂, CF₃) and BrY (Y = Br, Cl, F). Calculated results show that all substituent can form halogen‐bond complexes while part substituent can form pnicogen‐bond complexes. Traditional, chlorine‐shared and ion‐pair halogen‐bonds complexes have been found with the different substituent X and Y. The halogen‐bonds are stronger than the related pnicogen‐bonds. For halogen‐bonds, strongly electronegative substituents which are connected to the Lewis acid can strengthen the bonds and significantly influenced the structures and properties of the compounds. In contrast, the substituents which connected to the Lewis bases can produce opposite effects. The interaction energies of halogen‐bonds are 2.56 to 32.06 kcal·mol^?1; The strongest halogen‐bond was found in the complex of PH₂OH???BrF. The interaction energies of pnicogen‐bonds are in the range 1.20 to 2.28 kcal·mol^?1; the strongest pnicogen‐bond was found in PH₂Br???Br₂ complex. The charge transfer of lp(P) ? σ*(Br? Y), lp(F) ? σ*(Br? P), and lp(Br) ? σ*(X? P) play important roles in the formation of the halogen‐bonds and pnicogen‐bonds, which lead to polarization of the monomers. The polarization caused by the halogen‐bond is more obvious than that by the pnicogen‐bond, resulting in that some halogen‐bonds having little covalent character. The symmetry adapted perturbation theory (SAPT) energy decomposition analysis showes that the halogen‐bond and pnicogen‐bond interactions are predominantly electrostatic and dispersion, respectively.     

Quantum Study on the Hydrogen Bonding Interaction of Luteolin-（CH_3OH）_n

The optimized geometries and vibration frequencies of luteolin,methanol and luteolin-(CH3OH)n complexes have been investigated by density functional theory using B3LYP method.Four stable luteolin-CH3OH complexes,six stable luteolin-(CH3OH)2 complexes and four stable luteolin-(CH3OH)3 complexes have been obtained.The theories of atoms in molecules(AIM) and natural bond orbital(NBO) have been used to analyze the hydrogen bonds of these compounds,and their interaction energies corrected by basis set superposition error are between-8.046 and-76.124 kJ/mol.The calculation results indicate strong hydrogen bonding interactions in the luteolin-(CH3OH)n complexes.Then the nuclear magnetic resonance(NMR) and electronic absorption spectrum of luteolin have been calculated,and the results are in agreement with the experimental data.     

Theoretical Study of the Interplay between Lithium Bond and Hydrogen Bond in Complexes Involved with HLi and HCN

The lithium‐ and hydrogen‐bonded complex of HLi? NCH? NCH is studied with ab initio calculations. The optimized structure, vibrational frequencies, and binding energy are calculated at the MP2 level with 6‐311++G(2d,2p) basis set. The interplay between lithium bonding and hydrogen bonding in the complex is investigated with these properties. The effect of lithium bonding on the properties of hydrogen bonding is larger than that of hydrogen bonding on the properties of lithium bonding. In the trimer, the binding energies are increased by about 19 % and 61 % for the lithium and hydrogen bonds, respectively. A big cooperative energy (?5.50 kcal mol^?1) is observed in the complex. Both the charge transfer and induction effect due to the electrostatic interaction are responsible for the cooperativity in the trimer. The effect of HCN chain length on the lithium bonding has been considered. The natural bond orbital and atoms in molecules analyses indicate that the electrostatic force plays a main role in the lithium bonding. A many‐body interaction analysis has also been performed for HLi? (NCH)_N (N=2–5) systems.     

The influence of CH…π interaction on hydrogen bonding ability of –CONH<Subscript>2</Subscript> functional group of benzamide

Interplay between CH…π and hydrogen bond interactions of benzamide has been investigated by quantum mechanical calculations. The effect of the substituents on geometrical parameters has also been studied at the B3LYP level with 6-311++G(d,p) basis set. The electron-withdrawing substituents enhance the total interaction energy of the complexes. The results indicated that the cooperativity of interactions leads to extra stability of the ternary complexes. The CH…π interaction and the hydrogen bond energies have been estimated using the electron densities calculated by the atoms in molecules (AIM) method at hydrogen bond critical points. The strength of hydrogen bonding increases in the presence of CH…π interaction in the ternary complexes. The effect of CH…π interaction on the hydrogen bond interaction has also been studied by the natural bond orbital, AIM and the molecular electrostatic potential analyses.     

Density functional studies on hydrogen‐bonded clusters of hydrogen halides and the interaction on halide anions

Density functional theory (DFT) calculations have been performed to study the structures and stability of X^?·(HX)_n=2–5 clusters where X = F, Cl, Br at B3LYP/6‐311++G** level of theory. The presence of halide ions in these clusters disintegrates the hydrogen halide clusters. All the hydrogen halides are then hydrogen bonded to the centrally placed halide ions, thereby forming multiple hydrogen bonds. The interaction energies have been corrected for the basis set superposition error (BSSE) using Boy's counterpoise correction method. Evidence for the destruction of hydrogen bonds in hydrogen halide clusters due to the presence of halide ions is further obtained from topological analysis and natural bond orbital analysis. The chemical hardness and chemical potential have been calculated for all the anion clusters. The above analysis reveals that hydrogen bonding in these systems is not an essentially electrostatic interaction. The nature of the stabilization interactions operative in these multiple hydrogen‐bonded clusters has been explained in terms of many‐body contribution to interaction energies. From these studies, an attempt has been made to understand the nature of the molecular properties resulting from different electronegativities of the halogens. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Density functional theory and topological analysis on the hydrogen bonding interactions in cysteine‐thymine complexes

Hydrogen bonding interactions between amino acids and nucleic acid bases constitute the most important interactions responsible for the specificity of protein binding. In this study, complexes formed by hydrogen bonding interactions between cysteine and thymine have been studied by density functional theory. The relevant geometries, energies, and IR characteristics of hydrogen bonds (H‐bonds) have been systematically investigated. The quantum theory of atoms in molecule and natural bond orbital analysis have also been applied to understand the nature of the hydrogen bonding interactions in complexes. More than 10 kinds of H‐bonds including intra‐ and intermolecular H‐bonds have been found in complexes. Most of intermolecular H‐bonds involve O (or N) atom as H‐acceptor, whereas the H‐bonds involving C or S atom usually are weaker than other ones. Both the strength of H‐bonds and the structural deformation are responsible for the stability of complexes. Because of the serious deformation, the complex involving the strongest H‐bond is not the most stable structures. Relationships between H‐bond length (ΔR_X‐H), frequency shifts (Δv), and the electron density (ρ_b) and its Laplace (?²ρ_b) at bond critical points have also been investigated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011