DFT study on the intermolecular interactions between Aun (n = 2–4) and thymine

Density functional B3LYP method with 6-31++G** basis set is applied to optimize the geometries of the luteolin, water and luteolin–(H₂O)_n complexes. The vibrational frequencies are also studied at the same level to analyze these complexes. We obtained four steady luteolin–H₂O, nine steady luteolin–(H₂O)₂ and ten steady luteolin–(H₂O)₃, respectively. Theories of atoms in molecules (AIM) and natural bond orbital (NBO) are used to investigate the hydrogen bonds involved in all the systems. The interaction energies of all the complexes corrected by basis set superposition error, are within −13.7 to −82.5 kJ/mol. The strong hydrogen bonding mainly contribute to the interaction energies, Natural bond orbital analysis is performed to reveal the origin of the interaction. All calculations also indicate that there are strong hydrogen bonding interactions in luteolin–(H₂O)_n complexes. The OH stretching modes of complexes are red-shifted relative to those of the monomer.     

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.     

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.     

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.     

BrO-H2O 及 HOBr-H2O复合物偶合方式的理论研究

在6-311++G(d,p)水平上采用四种方法讨论了两种BrO-H₂O和三种HOBr-H₂O复合物的构型性质。在两种BrO-H₂O复合物中，结合能为11.37–13.92 J/mol的复合物2 （电子态为²A′）最稳定，该复合物是通过BrO中的Br原子和水中的O原子结合的。三种HOBr-H₂O复合物中，复合物3和4的结合能约为16.30–21.32 J/mol，三种复合物的稳定顺序为：复合物3 ≈ 复合物4 > 复合物5。     

O?O bond homolysis in hydrogen peroxide

O O bond homolysis in hydrogen peroxide (H₂O₂) has been studied using theoretical methods of four conceptually different types: hybrid DFT (B3LYP, M06‐2X), double‐hybrid DFT (B2‐PLYP), coupled‐cluster (CCSD(T)), and multiconfigurational (CASPT2). In addition, the effects of basis set size have also been analyzed. For all of these methods, the O O bond homolysis in hydrogen peroxide has been found to proceed through hydrogen bonded radical pair complexes. Reaction barriers for collapse of the radical pairs to hydrogen peroxide are minute, leading to an overall very flat potential energy surface. However, hydrogen bonding energies in the radical pair complex expressed as the energy difference to two separate hydroxyl radicals are sizeable and exceed 10 kJ/mol for all theoretical methods considered in this study. © 2017 Wiley Periodicals, Inc.     

A comparative study of some lithium and hydrogen‐bonded complexes: Ab initio and QTAIM studies

The nature of the interactions of cyanide with lithium and hydrogen halides was investigated using ab initio calculations and topological analysis of electron density. The computed properties of the lithium‐bonded complexes RCN···LiX (R = H, F, Cl, Br, C?CH, CH?CH₂, CH₃, C₂H₅; X = Cl, Br) were compared with those of corresponding hydrogen‐bonded complexes RCN···HX. The results show that both types of intermolecular interactions are “closed‐shell” noncovalent interactions. The effect of substitution on the interaction energy and electron density at the bond critical points of the lithium and hydrogen bonding interactions is similar. In comparison, the interaction energies of lithium‐bonded complexes are more negative than those of hydrogen‐bonded counterparts. The electrostatic interaction plays a more important role in the lithium bond than in the hydrogen bond. On complex formation, the net charge and energy of the Li atom decrease and the atomic volume increases, while the net charge and energy of the H atom increase and the atomic volume decreases. © 2013 Wiley Periodicals, Inc.     

1,2,4-三氮杂苯-(H2O)n复合物氢键相互作用的密度泛函理论研究

用密度泛函理论方法在B3LYP/6-31++G**水平上对1,2,4-三氮杂苯-(H₂O)_n (n＝1, 2, 3)氢键复合物的基态进行了结构优化和能量计算, 结果表明复合物之间存在较强的氢键作用, 所有稳定复合物结构中形成一个N…H—O氢键并终止于弱O…H—C氢键的氢键水链的构型最稳定. 同时, 用含时密度泛函理论方法(TD-DFT)在TD-B3LYP/6-31++G**水平上计算了1,2,4-三氮杂苯单体及其氢键复合物的单重态第一¹(n, π^*)垂直激发能.     

Nonadditivity of methyl group in single‐electron hydrogen bond of methyl radical‐water complex

The nonadditivity of methyl group in the single‐electron hydrogen bond of the methyl radical‐water complex has been studied with quantum chemical calculations at the UMP2/6‐311++G(2df,2p) level. The bond lengths and interaction energies have been calculated in the four complexes: CH₃? H₂O, CH₃CH₂? H₂O, (CH₃)₂CH? H₂O, and (CH₃)₃C? H₂O. With regard to the radicals, tert‐butyl radical forms the strongest hydrogen bond, followed by iso‐propyl radical and then ethyl radical; methyl radical forms the weakest hydrogen bond. These properties exhibit an indication of nonadditivity of the methyl group in the single‐electron hydrogen bond. The degree of nonadditivity of the methyl group is generally proportional to the number of methyl group in the radical. The shortening of the C···H distance and increase of the binding energy in the (CH₃)₂CH? H₂O and (CH₃)₃C? H₂O complexes are less two and three times as much as those in the CH₃CH₂? H₂O complex, respectively. The result suggests that the nonadditivity among methyl groups is negative. Natural bond orbital (NBO) and atom in molecules (AIM) analyses also support such conclusions. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Intermolecular interactions in uracil–nitrous acid complexes: structures,binding energy,topological properties,and nuclear magnetic resonance study

Molecular interactions between uracil and nitrous acid (U–NA) [C₄N₂O₂H₄? NO₂H] have been studied using B3LYP, B3PW91, and MP2 methods with different basis sets. The optimized geometries, harmonic vibrational frequencies, charge transfer, topological properties of electron density, nucleus‐independent chemical shift (NICS), and nuclear magnetic resonance one‐ and two‐bonds spin–spin coupling constants were calculated for U–NA complexes. In interaction between U and NA, eight cyclic complexes were obtained with two intermolecular hydrogen bonds N(C)HU…N(O) and OH_NA…O_U. In these complexes, uracil (U) simultaneously acts as proton acceptor and proton donor. The most stable complexes labeled, UNA1 and UNA2, are formed via NH bond of U with highest acidity and CO group of U with lowest proton affinity. There is a relationship between hydrogen bond distances and the corresponding frequency shifts. The solvent effect on complexes stability was examined using B3LYP method with the aug‐cc‐pVDZ basis set by applying the polarizable continuum model (PCM). The binding energies in the gas phase have also been compared with solvation energies computed using the PCM. Natural bond orbital analysis shows that in all complexes, the charge transfer takes place from U to NA. The results predict that the Lone Pair (LP)(O)_U → σ*(O? H) and LP(N(O)_NA → σ*(N(C)? H)_U donor–acceptor interactions are most important interactions in these complexes. Atom in molecule analysis confirms that hydrogen bond contacts are electrostatic in nature and covalent nature of proton donor groups decreases upon complexation. The relationship between spin–spin coupling constant (^1hJ_H…_Y and ^2hJ_H…_Y) with interaction energy and electronic density at corresponding hydrogen bond critical points and H‐bonds distances are investigated. NICS used for indicating of aromaticity of U ring upon complexation. © 2013 Wiley Periodicals, Inc.     

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.     

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     

Theoretical study on the interlay of hydrogen bonds in the trimers involving HCN and water

In this study, the hydration of hydrogen cyanide (HCN) has been investigated by means of quantum chemical ab initio calculations at the MP2/6‐311++G(3df,2p) level. Various HCN· · ·( H₂O)₂ and (HCN)₂· · ·H₂O complexes were optimized. Geometrics and energetics in these complexes have been analyzed. The hydration of the H atom leads to an elongation of the N?C and C? H bonds, whereas the hydration of the N atom results in a contraction of the N?C bond and a little elongation of the C? H bond. The interaction energy between each molecule pair in the trimers (except in HCN? H₂O? H₂O trimer) is increased relative to that in the respective dimer. The cooperativity of hydrogen bond in HCN? H₂O? H₂O trimer plays a negative contribution to the total interaction energy of the complex, whereas that in the other trimers is a positive contribution. Geometry and energy in H₂O? H₂OO? HCN? H₂O tetramer have also been analyzed. The binding energies in the trimers and tetramer have been studied by means of many‐body interaction analysis. The mechanism of HCN hydration was suggested. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

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     

Dual Role of the 1,2,3‐Triazolium Ring as a Hydrogen‐Bond Donor and Anion–π Receptor in Anion‐Recognition Processes

Several bis(triazolium)‐based receptors have been synthesized as chemosensors for anion recognition. The central naphthalene core features two aryltriazolium side‐arms. NMR experiments revealed differences between the binding modes of the two triazolium rings: one triazolium ring acts as a hydrogen‐bond donor, the other as an anion–π receptor. Receptors 9²⁺?2BF₄ ^? (C₆H₅), 11²⁺?2BF₄ ^? (4‐NO₂?C₆H₄), and 13²⁺?2BF^4? (ferrocenyl) bind HP₂O₇^3? anions in a mixed‐binding mode that features a combination of hydrogen‐bonding and anion–π interactions and results in strong binding. On the other hand, receptor 10²⁺?2 BF₄ ^? (4‐CH₃O?C₆H₄) only displays combined Csp₂?H/anion–π interactions between the two arms of the receptors and the bound anion rather than triazolium (CH)⁺???anion hydrogen bonding. All receptors undergo a downfield shift of the triazolium protons, as well as the inner naphthalene protons, in the presence of H₂PO₄^? anions. That suggests that only hydrogen‐bonding interactions exist between the binding site and the bound anion, and involve a combination of cationic (triazolium) and neutral (naphthalene) C?H donor interactions. Theoretical calculations relate the electronic structure of the substituent on the aromatic group with the interaction energies and provide a minimum‐energy conformation for all the complexes that explains their measured properties.     

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.     

Density functional study of HO(H2O)n (n = 1–3) clusters

The hydrogen bonding complexes HO(H₂O)_n (n = 1–3) were completely investigated in the present study using DFT and MP2 methods at varied basis set levels from 6‐31++G(d,p) to 6‐311++G(2d,2p). For n = 1 two, for n = 2 two, and for n = 3 five reasonable geometries are considered. The optimized geometric parameters and interaction energies for various complexes at different levels are estimated. The infrared spectrum frequencies and IR intensities of the most stable structures are reported. Finally, thermochemistry studies are also carried out. The results indicate that the formation and the number of hydrogen bonding have played an important role in the structures and relative stabilities of different complexes. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

4-羟甲基吡啶－水红移氢键的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-311++G**基函数对4-羟甲基吡啶与水形成的1:1和1:2（摩尔比）氢键复合物进行了理论计算研究,分别得到稳定的4-羟甲基吡啶－H₂O和4-羟甲基吡啶－(H₂O)₂氢键复合物3个和8个。经基组重叠误差和零点振动能校正后,最稳定的1:1和1:2氢键复合物的相互作用能分别为-20.536和-44.246 kJ/mol。振动分析显示O-H···N(O)氢键的形成使复合物中O-H键对称伸缩振动频率红移(减小)。自然键轨道分析表明,4-羟甲基吡啶与水形成最稳定的1:1和1:2氢键复合物时,分子间电荷转移分别为0.02642 e 和0.03813 e 。含时密度泛函理论TD-B3LYP/ 6-311++G**计算显示,相对于4-羟甲基吡啶单体分子,氢键H-OH···N和H-OH···OH的形成分别使最大吸收光谱波长兰移8~16纳米和红移4~11纳米。     

Characterization of two types of intermolecular interactions on halogen monoxide monohydrates

Monohydrates of halogen monoxides ClO·H₂O and BrO·H₂O have been studied by means of DFT (B3LYP) and ab initio (MP2) correlated calculations with aug‐cc‐pVnZ basis sets ranging from triple‐ up to quintuple‐ζ. These complexes might be formed in the troposphere and stratosphere and participate in chemical reactions involved in ozone depletion. Two stable structures are found that differ in the intermolecular interaction which takes place, namely: conventional XO…HOH hydrogen bond and OX…OH₂ halogen bond. We demonstrate that both types of interactions participate in the formation of these complexes yet all the computational methods tested predict a slightly greater stability for the latter OX…O link. Both intermolecular interactions are characterized upon analyzing electron density distribution, charge transfer effects, and electron localization domains. These analyses reveal the central role played by electron redistribution. Because of this, the greater spatial extent of the electron density in Cl or Br as compared to H could be the main cause to yield a slightly greater stability for the O? X…O halogen bond with respect to the O…H? O hydrogen bond. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Theoretical Study on the Structures and Properties of Hydrogen Bonding Complexes between Diazines and Water

Density functional theory B3LYP method and second-order Moller-Plesset perturbation theory MP2 method were employed to obtain the optimized geometries of the ground state and interaction energy for diazines and water complexes. The results show that the ground state complexes have strong hydrogen bonding interaction with -20.99, -16.73 and -15.31 kJ/mol after basis set superposition error and zero-point vibration energy correction for pyridazine-water, pyrimidine-water and pyrazine-water, respectively, and large red-shift for the symmetric H-O stretching vibration frequencies due to the formation of N…H-O hydrogen bond in the diazine-water complexes. The NBO analysis indicates that intermolecular charge transfer are 0.0316, 0.0255 and 0.0265 e respectively. In addition, the first singlet （n,n＊） vertical excitation energy of the monomer and the hydrogen bonding complexes between diazines and water was investigated by time-dependent density functional theory.