Hydrogen-bonding interactions in acetonitrile oligomers using density functional theory method

Hydrogen-bonding interaction in acetonitrile oligomers is studied using density functional theory method. Two types of hydrogen-bonded oligomers are considered viz. cyclic and ladder. Different levels are used to optimize the geometry of acetonitrile monomer and found that at B3LYP/aug-cc-pvtz level the geometrical parameters and vibrational frequencies are in agreement with the experimental determinations. The BSSE corrected total energies of acetonitrile oligomers show that the cyclic structures are more stable than the ladder and the hydrogen bonds in former are stronger than those in the latter. Many-body analysis approach was used to study the nature of interactions between different molecules in these oligomers. It is found that the contribution from many-body energies to the binding energy of a complex is different in cyclic and ladder structures. An increase and decrease in the energy per hydrogen bond with cluster size for the cyclic and ladder structures, respectively, indicates the positive and negative hydrogen-bond cooperativity, respectively.     

Density functional theory study of formaldehyde oligomers

Hydrogen‐bonded formaldehyde oligomers (dimer to pentamer) are studied using density functional theory (DFT), the B3LYP method, and the 6‐311+G* basis set. Many‐body interaction energies are obtained to study the contribution of many‐body terms to binding energy. The basis set superposition error (BSSE)‐corrected total energies are ?229.08170, ?343.61410, ?458.16660, and ?572.70901 hartrees for dimer, trimer, tetramer, and pentamer, respectively, with corresponding binding energies ?2.55, ?4.86, ?6.99, and ?9.49 kcal/mol. Two‐body energies have been found to contribute significantly to the total binding energy in dimer to pentamer, whereas higher‐order interaction energies are negligible. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.     

DFT study on cooperativity in the interactions of hydrazoic acid clusters

Density functional theory at the B3LYP level with the 6‐311G** basis set is performed to calculate the systems consisting of up to four hydrazoic acid molecules. The dimers are found to exhibit cyclic and chain structures with N … H contacts at ca. 2.1–2.7 Å. However, there are only cyclic structures with N … H contacts at ca. 2.0–2.3 Å and 2.0–2.1 Å in the trimer and tetramer, respectively. Hydrogen bond distances in the trimer and tetramer are shorter than those in the cyclic dimer as a result of the stronger interaction between molecules. The contribution of cooperative effect to the interaction energy is significant. After the correction of the basis set superposition error and zero‐point energy, the binding energies are ?10.69, ?29.34, and ?54.26 kJ·mol^?1 for the most stable dimer, trimer, and tetramer, respectively. The calculated IR shifts for N? H stretching mode increase with the size of the cluster growths, reaching more than 200 cm^?1 in the tetramer. For the most stable clusters, the transition from the monomer to dimer, dimer to trimer, and trimer to tetramer involve changes of ?14.40, ?25.68, and ?31.88 kJ·mol^?1 for the enthalpies at 298.15 K and 1atm, respectively. We also perform Mulliken populations analysis and find the Mulliken populations on intermolecular N … H increasing in the sequence of the dimer, trimer, and tetramer. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 279–286, 2003     

Thermodynamics of forming water clusters at various temperatures and pressures by Gaussian-2, Gaussian-3, complete basis set-QB3, and complete basis set-APNO model chemistries; implications for atmospheric chemistry

The Gaussian-2, Gaussian-3, complete basis set- (CBS-) QB3, and CBS-APNO methods have been used to calculate Delta H degrees and Delta G degrees values for neutral clusters of water, (H(2)O)(n), where n = 2-6. The structures are similar to those determined from experiment and from previous high-level calculations. The thermodynamic calculations by the G2, G3, and CBS-APNO methods compare well against the estimated MP2(CBS) limit. The cyclic pentamer and hexamer structures release the most heat per hydrogen bond formed of any of the clusters. While the cage and prism forms of the hexamer are the lowest energy structures at very low temperatures, as temperature is increased the cyclic structure is favored. The free energies of cluster formation at different temperatures reveal interesting insights, the most striking being that the cyclic trimer, cyclic tetramer, and cyclic pentamer, like the dimer, should be detectable in the lower troposphere. We predict water dimer concentrations of 9 x 10(14) molecules/cm(3), water trimer concentrations of 2.6 x 10(12) molecules/cm(3), tetramer concentrations of approximately 5.8 x 10(11) molecules/cm(3), and pentamer concentrations of approximately 3.5 x 10(10) molecules/cm(3) in saturated air at 298 K. These results have important implications for understanding the gas-phase chemistry of the lower troposphere.     

Theoretical Study on Hydrogen Interaction between 5-Fluorouracil and Glycine

The geometries,electronic structure,IR spectrum and other properties of hydrogen interaction between 5-fluorouracil and glycine were studied at the B3LYP/6-31+G* level.Single point energy calculations were executed at the B3LYP/6-311++G** and B3LYP/aug-cc-pvdz levels,and natural bond orbital (NBO) analysis was carried out at the B3LYP/6-31+G* level.Finally,the hydrogen bonds were discussed via AIM electronic density topology analysis.     

Molecular structure and infrared spectra of (HXeCN)n (n=2, 3 or 4)

The structure, energetics, and vibrational spectra of the (HXeCN)2 dimer were investigated at the CCSD(T), MP2 and B3LYP levels. Such properties of the (HXeCN)3 trimer and (HXeCN)4 tetramer were investigated at the B3LYP level. The dimer, trimer, and tetramer were predicted to have a C2h, C2v, and D2d structure, respectively. In all of these oligomers, the N?Xe intermonomeric interaction is the most important one for holding the monomers together. Included with the ZPVE and BSSE, the stabilization energy of the dimer is 12.36 kcal/mol at the CCSD(T) level, while those of the dimer, trimer, and tetramer are 10.42, 18.23, and 31.34 kcal/mol, respectively, at the B3LYP level. At the B3LYP level, with respect to those of the isolated monomer, the C-Xe and Xe-H asymmetric stretching frequencies are shifted by -11.2 and +128.0 cm(-1) for the dimer, -51.6, +220.7 and -11.5, +96.6 cm(-1) for the trimer, and -14.1 and +201.8 cm(-1) for the tetramer.     

NTO二聚体分子间相互作用的理论研究

在DFT-B3LYP/6-311＋＋G**水平上求得NTO二聚体势能面上六种优化构型和电子结构. 经基组叠加误差(BSSE)和零点能(ZPE)校正, 求得分子间最大相互作用能为－53.66 kJ/mol. 二子体系间的电荷转移很少. 由自然键轨道分析揭示了相互作用的本质. 对优化构型进行振动分析, 并基于统计热力学求得200.0～800.0 K温度范围从单体形成二聚体的热力学性质变化. 发现二聚主要由强氢键所贡献, 但结合能大小并不为氢键所完全决定. 二聚过程在较低温度或常温下能自发进行.     

Interaction with glycine increases stability of a mutagenic tautomer of uracil. A density functional theory study

The most stable structures for the gas-phase complexes of minor tautomers of uracil (U) with glycine (G) were characterized at the density functional B3LYP/6-31++G level of theory. These are cyclic structures stabilized by two hydrogen bonds. The relative stability of isolated tautomers of uracil was rationalized by using thermodynamic and structural arguments. The stabilization energies for complexes between the tautomers of U and G result from interplay between the stabilizing two-body interaction energies and destabilizing one-body terms. The latter are related to the energies of (i) tautomerization of the unperturbed moieties and (ii) distortions of the resulting rare tautomers in the complex. The two-body term describes the interaction energy between distorted tautomers. The two-body interaction energy term correlates with perturbations of length of the proton-donor bonds as well as with deprotonation enthalpies and proton affinities of the appropriate monomer sites. It was demonstrated that the relative instability of rare tautomers of uracil is diminished due to their interactions with glycine. In particular, the instability of the third most stable tautomer (U(III)) is decreased from 11.9 kcal/mol for non-interacting uracil to 6.7 kcal/mol for uracil in a complex with the zwitterionic tautomer of glycine. A decrease of instability by 5.2 kcal/mol could result in an increase of concentration of U(III) by almost 5 orders of magnitude. This is the tautomer with proton donor and acceptor sites matching guanine rather than adenine. Moreover, kinetic characteristics obtained for the glycine-assisted conversion of the most stable tautomer of uracil (U(I)) to U(III) indicate that the U(I)<-->U(III) thermodynamic equilibrium could be easily attained at room temperature. The resulting concentration of this tautomer falls in a mutationally significant range.     

环双(对-蒽基-对草快)的分子识别与谱学性质

环双(对-蒽基-对草快)是一种新型的缺电子大环仿生主体, 分子识别是其最重要的应用之一. 考察主体对一系列客体分子如水、氨、醇及杂环等的识别能力, 用密度泛函理论(DFT)中的B3LYP/3-21G基组对主客体复合物的结构进行优化. 在B3LYP/6-31G(d)水平上进行单点能计算, 校正后得到复合物的结合能. 用B3LYP/3-21G方法计算13C和3He化学位移. 结果表明, 主体对客体分子的识别主要靠客体上的杂原子与主体上的氢原子之间的氢键进行. 复合物的稳定化能受氢键的数目和距离影响. 氢键的形成导致部分复合物LUMO与HOMO能隙增大, 同时导致与氢键相连的C—H键上C原子的化学位移向低场移动. 复合物的芳香性与其结合能的大小及结合方式有关. 主体的芳香性因其与客体之间的弱相互作用而提高, 但太强的相互作用及客体在主体空腔内都将影响主体的环电流, 从而削弱其芳香性.     

An exploration of conformational search of leucine molecule and their vibrational spectra in gas phase using ab initio methods

An extensive exploration of the conformational space has been carried out to characterize all possible gas phase structures of leucine. A total of 324 unique trial structures for canonical leucine were generated by considering all possible combinations of single bond rotamers. All trial structures were optimized at the B3LYP/6-311G* level of the DFT method. A total of 77 unique and stationary canonical conformers were found. Further, 15 most stable conformers were reoptimized at B3LYP/6-311++G** level and their respective relative energies, vertical ionization energies, hydrogen bonding patterns, rotational constants and dipole moments were calculated. A single point energy calculations for leucine conformers have also been done at both B3LYP/6-311++G(2df, p) and MP2/6-311++G(2df, p) levels. The good agreement between our estimates of rotational constants for two most stable conformers and available experimental measurements supports the reliability of the B3LYP/6-311++G** level of theory for describing the conformational behavior of leucine molecule. The proton affinity and gas phase basicity were also determined. Using the statistical approach, conformational distributions at various temperatures have also been performed and analyzed. Vibrational spectra were also calculated. It is also observed that zwitterions of leucine are not stable in gas phase.     

H‐bonding and stacking interactions between chloroquine and temozolomide

The interactions between temozolomide and chloroquine were examined via Dispersion‐Corrected Density Functional Theory and MP2 methods. Chloroquine was considered in both its lowest energy structure and in a local minimum where its aromatic system and secondary amine group are free to interact directly with temozolomide. The accessibility of these two components to intermolecular interaction makes the lowest energy dimer of this local monomer minimum competitive in total energy with that involving chloroquine's most stable monomer geometry. In either case, the most stable heterodimer places the aromatic ring systems of the two molecules parallel and directly above one another in a stacked geometry. Most of the local minima are also characterized by a stacked geometry as well. Comparison between B3LYP and B3LYP‐D binding energies confirms dispersion is a primary factor in stabilizing these structures. © 2016 Wiley Periodicals, Inc.     

Interaction between RNA segment (adenine‐uracil) and model of protein unit (formamide): A density‐functional theory study

Various possible structures of adenine‐uracil‐formamide hydrogen‐bond complexes were optimized at 6‐311++G(d,p) level, and the binding energies of these complexes were also calculated at DFT B3LYP/6‐311++G(d,p) level. Eight stable cyclic structures being involved in the interaction are found on the potential energy surface. By analyzing the structure, NPA charge and interaction energy of complexes, we obtain the most stable geometry structure. The results show that the interactions between formamide and adenine‐uracil (A‐U) base pair affect the stabilities of the base pairs. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

平面五元水分子簇的量子化学研究

选用Gaussian03的B3LYP/6-31G(d,p)、DMol3的BLYP/DNP和deMon的BLYP/TZVP等方法计算了甲烷水合物(结构-1)中平面五元水分子簇的结合能和氢键能,作了基组重叠误差(BSSE)和色散能(dispersion)的修正,估算了次级相互作用的贡献.在DMol3程序中使用了大型数值基组DNP,将基组重叠误差降至最低.在Gaussi-an03的B3LYP/6-31G(d,p)计算中,采用平衡法(Counterpoise)校正基组重叠误差.两种计算方法给出了一致的结果,证实了在使用6-31G(d,p)基组时,一对水分子在平衡距离的基组重叠误差高达8 kJ/mol.为估算色散能的贡献,使用了新近发展的包含色散能的密度泛函的DFT程序deMon计算了五元水分子簇.用多种方法计算出了经基组重叠误差和色散能修正的五元水分子簇的分子间结合能和氢键能的较为精确的势能超曲面,为甲烷和其他气体水合物的分子动力学模拟提供了依据.     

Electric field effects on water clusters (n = 3-5): systematic ab initio study of structures, energetics, and transition states

The structures, energetics, and transition states of water clusters (trimer to pentamer, n = 3-5) are investigated as a function of electric field by using ab initio calculations. With an increasing strength of the field, the most stable cyclic structures of trimer, tetramer, and pentamer open up to align their dipole moments along the direction of the field. For the lower strength (below 0.3 V/angstroms) of the electric field, the dipole moment of each water monomer is along the same direction with the field, while it retains the cyclic structure. For the higher strength of the field, to have a higher dipole moment for the cluster along the field direction, each cyclic structure opens up to form a linear chain or "water wire." We have investigated the transition state structures between the cyclic and linear forms for the field strengths of 0.3-0.4 V/angstroms where both cyclic and linear forms are energetically comparable.     

A density-functional theory investigation of 3-nitro-1,2,4-triazole-5-one dimers and crystal

Density-functional method with different basis sets was applied to the study of the highly efficient and low sensitive explosive 3-nitro-1,2,4-triazole-5-one (NTO) in both gaseous dimer and its bulk state. The binding energies have been corrected for the basis set superposition errors. Six stable dimers (II-VII) were located. The corrected binding energy of the most stable dimer VII is predicted to be -53.66 kJ/mol at the B3LYP/6-311++G(**) level. It was found that the structures of the more stable dimers (V-VII) are through the hydrogen bonding interaction between the carbonyl oxygen and the azole hydrogen of 3-nitro-1,2,4-triazole-5-one. The changes of Gibbs free energies (DeltaG) in the processes from the monomer to the dimers at 298.15 K are 8.51, 0.90, 0.35, -8.74, -10.67, and -11.06 kJ/mol for dimers from II to VII, respectively. Dimers V-VII, possessing cyclic structures, can be spontaneously produced from the isolated monomer at room temperature. The lattice energy is -156.14 kJ/mol, and this value becomes to -150.43 kJ/mol when a 50% correction of the basis set superposition error was adopted. The frontier bands are quite flat. Judged from the value of band gap of 4.0 eV, it may be predicted that 3-nitro-1,2,4-triazole-5-one is an insulator. Most atoms in NTO, with the exception of C(5) atom and the nitro atoms, make up the upper valence bands. In contrast, the lower conduction bands mainly consist of the nitro N and O atoms. The population of the C-NO(2) bond is much less than those of the other bonds and the detonation may be initiated by the breakdown of this bond.     

Density functional theory study of self-association of N-methylformamide and its effect on intramolecular and intermolecular geometrical parameters and the cis/trans population

The single-point total energy (E) of several acyclic and cyclic oligomers of N-methylformamide (NMF) was computed by the first time without any geometrical restriction, using the B3LYP6-31G* method of the density functional theory in order to determine the effect of self-association on intramolecular geometrical parameters of cis- and trans-NMF, the intermolecular distances of the hydrogen-bonding chains formed by NMF as well as intermolecular association energies including counterpoise corrections. It is concluded that liquid NMF exists mainly as polymers formed by self-association of trans-NMF units, whereas the cis-NMF isomer occurs as isolated units inserted along the chains. These computational results are in accordance with the experimentally determined predominance (ca. 90%) of trans-NMF population by means of (1)H- NMR and other spectroscopic techniques, but in severe contradiction with a recent interpretation of x-ray diffraction data on liquid NMF, postulating a cyclic trimer of cis-NMF (c-C(3)) as the predominating species. The counterpoise-corrected values of the association energy, DeltaE(CP), calculated for cyclic oligomers, increase with the polymerization degree (n) revealing a high grade of cooperative effect for amidic hydrogen-bonded chains. Noteworthy, the difference between the DeltaE(CP) values of the cyclic cis- and trans-homooligomers of NMF is positive for n=2 and 3 but negative for n > or =4.     

One‐Pot Multimolecular Macrocyclization for the Expedient Synthesis of Macrocyclic Aromatic Pentamers by a Chain Growth Mechanism

POCl₃‐mediated one‐pot macrocyclization allows the highly selective formation of five‐residue macrocycles that are rigidified by internally placed intramolecular hydrogen bonds. Mechanistic investigation by using tailored competition experiments and kinetic simulation provides a comprehensive model, supporting a chain‐growth mechanism underlying the one‐pot formation of aromatic pentamers, whereby the successive addition of a bifunctional monomer unit onto either another monomer or the growing oligomeric backbone is faster than other types of bimolecular condensations involving oligomers longer than monomers. DFT calculations at the B3LYP/6‐31G* level reveal the five‐residue pentamer to be the most stable with respect to alternative four‐, six‐, and seven‐residue macrocycles. These novel mechanistic insights may become useful in analyzing other macrocyclization, oligomerization, and ploymerization reactions.     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Difluoroamine Complex

Introduction　　ApplicationsofabinitiocalculationstointermolecularinteractionsincludingeitherweakvanderWaalsorstrongerhydrogenbondinghavedrawnmuchattentioninthepastdecadesbecausetheyareimportantinawiderangeofphysical,chemicalandbiologicalfields .1 5Inre centyears ,wehaveappliedtheintermolecularinteractionstoenergeticsystemsandobtainedsomemeaningfulinfor mationthatisvaluableforthestudyofenergeticmateri als .6 14 Thebehaviorofmolecularcomplexesisusuallybe tweentwoextremes :thegasphaseandthecrys…     

Mechanisms of hydrolysis-oligomerization of aluminum alkoxide Al(OC3H7)3

As one of the representative superinsulating materials, the aluminum trioxypropyl Al(OC(3)H(7))(3) aerogel may be applied in launch vehicles and manned spacecrafts. In this study, the structures and hydrolysis mechanisms of the monomer, dimers, and trimers of Al(OC(3)H(7))(3) in neutral and alkaline environments were studied at the B3LYP/6-31G(d,p) level by using the CPCM solvation model to understand the fundamental chemistry of Al(OC(3)H(7))(3) hydrolysis and oligomerization. Our calculation shows that the first-order hydrolyses of the monomer and oligomers are energetically favorable in both alkaline and neutral solutions. In alkaline solutions, they are more apt to oligomerize than to hydrolyze due to high energy barriers and large binding energies in the formation of anionic species. For the oligomers under neutral condition (1) Al(OC(3)H(7))(3) is linked by four-membered Al-O rings with pentacoordinated bridging and tetracoordinated Al atoms, (2) the hydrolyzed propoxy groups will be expelled by solvent molecules, and (3) partly hydrolyzed species can condense to oligomers with bridging OH groups or O atoms.