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     

Optical rotation and two chiral C atoms of aliphatic cyclic dipeptides

Accurate geometries structures and total energies have been determined for the conformers of cyclo(L-Pro-Gly), cyclo(L-Ala-L-Ala), and cyclo(L-Pro-Ala) in the gaseous phase, using HF and B3LYP correlation methods at 6−31++G(d), 6−311++G(d, p), 6−311++G(2d, 2p) and aug-cc-pvdz basis sets. High level computations MP2 with 6−311++G(2d, 2p) basis set indicate that the relative stabilities of the available conformers can be determined correctly at the B3LYP/6−311++g(2d, 2p) level of theory. We have also described the implementation of DFT and HF theory for calculations of the optical rotation at 589.3 nm. In L-Ala-L-Ala, and L-Pro-Ala molecules, they have two chiral C (C*), so we discuss the different effect of two chiral C to optical activity of cydo(L-Pro-Gly), cyclo(L-Ala-L-Ala), and cyclo(L-Pro-Ala).     

HNO(H<Subscript>2</Subscript>O)<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript> (<Emphasis Type="Italic">n</Emphasis> = 1–4) clusters: a theoretical study

The geometrical structure, binding energy, and vibrational spectra of small clusters of nitrosyl hydride (HNO) and water molecules, HNO(H₂O) _n , where (n = 1–4), have been investigated at the MP2 level of theory, using 6-311++G(2d,2p) basis set. We located three dimers, six trimers, nine tetramers, and three pentamers at the MP2/6-311++G(2d,2p) computational level. Particular attention is given to existence and magnitude of NH···O blue-shifting hydrogen bonds. Blue shifts of the NH stretching frequency upon complex formation in the ranges between 28 and 151 cm⁻¹ is predicted. Cooperative effect in terms of stabilization energy along with the many-body interaction energies analysis was performed for the studied clusters. The Atoms in Molecules (AIM) theory was also applied to explain the nature of the complexes.     

Complexes of ammonia with propane and cyclopropane: electrostatic guidelines for ab initio treatment

A model based on the molecular electrostatic potential (MESP) is employed for the investigation of structures and energies of complexes of ammonia with propane and cyclopropane. The electrostatic model geometries are employed as starting points for an ab initio investigation at the self-consistent field and second-order M?ller-Plesset (MP2) levels. The most stable structures of C₃H₆..NH₃ and C₃H₈..NH₃ complexes have the interaction energies of 10.07 kJ/mol and 8.15 kJ/mol, respectively, at the MP2/6-31G(d,p) level. The energy rank order of the structures is not altered with the use of the 6-31++G(d,p) basis set, and the basis␣set superposition error has little effect. The interaction energy decomposition analysis shows that the electrostatic component is dominant over the other ones. MESP topography thus seems to offer valuable hints for predicting the structures of weakly bonded complexes. Received: 8 July 1998 / Accepted: 4 August 1998 / Published online: 2 November 1998     

A theoretical study of the interactions between N,N-dimethylformamide and xylenes

The ab initio and density functional (DFT) methods were performed on binary systems of N,N-dimethylformamide (DMF) with xylenes (o-, or m-, or p-xylene), and seven stable configurations were obtained with no imaginary frequencies. To obtain the interaction energies of these complexes, single-point energy calculations with basis set superposition error (BSSE) correction were carried out at B3LYP/6-31G* and MP2/6-31G* levels. The structures, Chelpg (charges from electrostatic potentials using a grid-based method) charge distribution and bond characteristics of the mentioned complexes were calculated. The results indicated the presence of double C–H···O hydrogen bonds between DMF and xylenes in these complexes and the interaction energies of hydrogen bonding between DMF and xylene systems decreased in the following sequence: DMF–o-xylene: a1 > DMF–m-xylene: b1 > DMF–p-xylene: c1.     

Self-assembly of DMF with chloromethane and their structures: a theoretical study

The stability of hydrogen-bonded complexes, DMF–H _n CCl_4−n (n = 1–3), has been investigated by several theoretical methods including the MP2 level of ab initio theory at various basis sets from 6-31+G* to 6-311++G**. Two stable configurations (respectively a and b) were obtained for each complex with no imaginary frequencies. The minimum energy structure of these complexes has also been analyzed by means of the atoms in molecule theory at MP2/6-311++G** level. It is found that C–H···O hydrogen bonding exists in these systems and that the intensity of HB interaction gradually increases with successive chlorination. Computed results indicate that these complexes automatically assemble into different stable configurations. For the complexes under consideration, their stabilities can be mainly ascribed to the intermolecular HB interaction. The present work is helpful to clearly understand the interaction mechanism of these complexes in theory.     

Intermolecular potentials and force constants from ab initio energies – Application to the N-H…O=C hydrogen bonds in formamide dimers

The ab initio energies and force constants of 38 geometrically optimized formamide dimers which differ in the lengths of the hydrogen bonds, are evaluated using the program GAUSSIAN 90 with the 6–31G** basis set. A potential energy function was fitted simultaneously to the dimerization energies (including vibrational energy contributions to association energies) and the force constants of the N-H…O=C bridge. As an application, the broadening of the signals in vibrational spectra of liquid formamide was simulated by a superposition of spectra of different formamide oligomers. Received: 12 November 1996 / Revised: 27 May 1997 / Accepted: 27 May 1997     

Searching blue-shifted O–H···Y hydrogen bond in CH<Subscript>3</Subscript>OH complexes with CF<Subscript>4</Subscript>, C<Subscript>2</Subscript>F<Subscript>2</Subscript>, OC,Ne, and He: a theoretical study

The hydrogen-bonded structures of the CH₃OH complexes with CF₄, C₂F₂, OC, Ne, and He are designated as the starting points for geometry optimizations without and with counterpoise (CP) correction at MP2 level of theory with the basis sets 6-31+G*, 6-31++G**, and 6-311++G**, respectively. Tight convergence criteria are applied throughout all geometry optimizations in order to reduce the computational errors. According to the optimizations without CP correction, a blue-shifted O–H···Y (where Y = F, O, Ne, or He) hydrogen bond exists in all these five complexes. The magnitudes of blue shifts of ν(O–H) of the former four complexes with respect to that of CH₃OH are reduced greatly when the polarization and diffuse functions of the hydrogen atoms are added (results from 6-31+G* versus those from 6-31++G**). However, for the complexes CH₃OH–CF₄ and CH₃OH–C₂F₂, our optimizations using the CP corrections did not find the hydrogen-bonded structure to be a stationary point. The energy minimum of both the complexes corresponds to a non-hydrogen-bonded structure.     

An ab initio molecular orbital study of the structures and energies of neutral and charged bimolecular complexes of NH3 with the hydrides AHn (A = N,O, F,P, S,and Cl)

Hartree-Fock 6-31G(d) structures for the neutral, positive ion, and negative ion bimolecular complexes of NH₃ with the first- and second-row hydrides AH_n (AH_n = NH₃, OH₂, FH, PH₃, SH₂, and ClH) have been determined. All of the stable neutral complexes except (NH₃)₂, the positive ion complexes with NH₃ as the proton acceptor, and the negative ion complexes containing first-row anions exhibit conventional hydrogen bonded structures with essentially linear hydrogen bonds and directed lone pairs of electrons. The positive ion complex NH₄⁺ …? OH₂ has the dipole moment vector of H₂O instead of a lone pair directed along the intermolecular line, while the complexes of NH₄⁺ with SH₂, FH, and ClH have structures intermediate between the lone-pair directed and dipole directed forms. The negative ion complexes containing second-row anions have nonlinear hydrogen bonds. The addition of diffuse functions on nonhydrogen atoms to the valence double-split plus polarization 6-31G(d,p) basis set usually decreases the computed stabilization energies of these complexes. Splitting d polarization functions usually destabilizes these complexes, whereas splitting p polarization functions either has no effect or leads to stabilization. The overall effect of augmenting the 6-31G(d,p) basis set with diffuse functions on nonhydrogen atoms and two sets of polarization functions is to lower computed stabilization energies. Electron correlation stabilizes all of these complexes. The second-order Møller–Plesset correlation term is the largest term and always has a stabilizing effect, whereas the third and fourth-order terms are smaller and often of opposite sign. The recommended level of theory for computing the stabilization energies of these complexes is MP2/6-31+G(2d,2p), although MP2/6-31+G(d,p) is appropriate for the negative ion complexes.     

DFT computational studies on rotation barriers,tautomerism, intramolecular hydrogen bond,and solvent effects in 8‐hydroxyquinoline

Intramolecular hydrogen binding interactions in 8‐hydroxyquinoline, both in its zwitterionic tautomer and in the rotamer without the intramolecular hydrogen bond (IHB), have been computed using the B3LYP and MPW1K density functionals. The rotation of the O? H bond and intramolecular proton transfer reactions were studied theoretically. The following theory levels have been applied: B3LYP/6‐31G(d,p), B3LYP/6‐311++G(d,p), MPW1K/6‐311++G(d,p), and MPW1K/6‐311++G(2d,3p)//MPW1K/6‐311++G(d,p). Natural bond orbital (NBO) analysis has also been carried out. The effect of medium (benzene, chloroform, tetrahydrofuran, 1,2‐dichloroethane, acetone, water) was simulated using the self‐consistent reaction field (SCRF) method within the framework of the polarizable continuum model (PCM), at the MPW1K/6‐311++G(d,p) level. The evolution of geometry, relative energies, heights of rotation (around the O? H bond) and tautomerization barriers, IHB energies, and ΔG(solv) have been systematically investigated. The results obtained have shown the failure to neglect some changes of the above characteristics in polar media with respect to the gaseous phase. The series of stability of the forms under study in the gaseous phase remains the same in solution. Thus, in spite of the important role of the solvent electrostatic effects, the intrinsic stability of those species overcomes the solvent effects. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

The influence of formamide (model of protein unit) on the intramolecular proton transfer in the DNA simple base guanine: A density functional theory study

For the purpose of investigating the influence of protein unit on the intramolecular proton transfer (IPT) reactions in the simple base guanine, a simple model (formamide) of peptides is designed to biological system investigations, and five complexes of formamide–guanine (FG1, FG2, FG3, FG4, and FG5) are determined at the B3LYP/6‐311++G(d,p) level of theory. For comparison, HF and MP2 methods are also used in this paper. The proton transfer (PT) reaction processes of guanine and FGs have been investigated employing the B3LYP/6‐311++G(d,p) level of theory. The selected thermodynamic and kinetic parameters, such as the activation energies (E_a), changes of enthalpy (ΔH) and changes of free energies (ΔG), as well as the equilibrium constants (K_p) for those reaction processes, have also been obtained by calculational means. The calculated results indicate that the assisted and protected effects of formamide on IPT in guanine are site‐dependent. CH1 is the lowest activation energy needed PT process no matter where the formamide molecule is located in. The activation energy of CH1 with formamide in S2 is the lowest one (153.3 KJ/mol), whereas the one of CH5 with formamide in S5 is the highest (318.3 KJ/mol). © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Ab initio study of substituent effect on the addition of hydrogen fluoride to fluoroethylenes

An ab initio 3-21G study of the direct addition of HF to C₂H_nF_(4–n), with n = 0 to 4, has been performed to investigate the effect of the substituent on the reaction. Geometry optimization of all charge-transfer complexes and transition states has been done. Standard analysis of activation energies of addition reactions, vibrational and thermodynamical analysis, as well as Morokuma energy decomposition, BSSE correction, PMO analysis, and Pauling bond orders were used to explain the results. A subset of the reactions, including that of C₂H₄ as reference one and the two most favorable cases, was also studied at the MP2/6–31G(d,p)//HF/6–31G(d,p) level. The barriers so obtained are in agreement with the indirectly found from experimental data. It was found that the effect of the substituent is not monotonic for the additions. Decomposition of the interaction energy is shown to be adequate to explain this nonmonotonic behavior. The implications for laser chemistry of the addition of hydrogen halides to fluorosubstituted olefins is briefly discussed.     

Chiral discrimination in hydrogen‐bonded complexes of hydrogen peroxide with methyl hydroperoxide: Theoretical study

For the first time, the discrimination of different chiral forms of 1:1 complexes with hydrogen peroxide and methyl hydroperoxide have been investigated using density functional theory (DFT) and Møller–Plesset type 2 (MP2) methods at varied basis set levels from 6‐31+G(d,p) to 6‐31++G(2d,2p). Three pairs of chiral enantiomers were considered. The optimized geometric parameters, interaction energies, and chirodiastatic energies for various isomers at different levels are estimated. To take into account the water solvation effect, the polarized continuum model (PCM) method has been used to evaluate the ΔG_solv. The gas phase results show that the heterochiral six‐membered ring complex (structure I) and homochiral five‐membered ring complexes (structures IV and V) are preferred configurations for the three pairs of chiral enantiomers. The solvation effect on six‐membered ring complexes (structures I and II) shows nonsignificant changes in the configurations preferred, but on five‐membered ring complexes, the homo‐/heterochiral preference is found to be inverse in the polar solvent. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Hydrogen-bonding interaction of urea with DNA bases: A density functional theory study

This work deals with the interaction between urea and DNA bases (adenine, thymine, guanine, and cytosine). The optimized geometries, binding energies, and harmonic vibrational frequencies are calculated using the DFT/B3LYP functional combined with the 6–31+G(d,p) basis set. Their interactions are studied aiming to understand more about the nature of the intercalation binding forces between urea and DNA. Fourteen stable complexes are found on the potential energy surface. The structures are cyclic; they are stabilized by NH...O/N and CH...O interactions. The binding energies range from −19.9 kJ·mol⁻¹ to −74.0 kJ·mol⁻¹. The obtained formation energies indicate that Urea:G and Urea:C are more favorable than Urea:T and Urea:A. In addition, the Atoms in Molecules theory is performed to study the hydrogen bonds in the complexes.     

Study of the formamide–methanol dimer with ab initio and density functional theory methods

For the first time, the structures and energies for the hydrogen bonding of a 1:1 complex formed between formamide and methanol molecules have been computed with various pure and hybrid density functional theory (DFT) and ab initio methods at varied basis set levels from 6‐31g to 6‐31+g(d,p). Five reasonable geometries on the potential energy surface of methanol and formamide system are considered and their relative stability is discussed. The infrared (IR) spectrum frequencies, IR intensities, and vibrational frequency shifts are reported. From the systematic studies, it is found that all the DFT methods selected here correctly compute the dimerization energies and geometries, with the B3P86 method predicting the hydrogen bond lengths relatively shorter and BPW91 yielding the interaction energies relatively lower. Finally, the solvent effects on the geometries of the formamide–methanol complexes have also been investigated using self‐consistent reaction field (SCRF) calculations with five different DFT methods at the 6‐31+g(d,p) basis set level. The results indicate that the polarity of the solvent has played an important role on the structures and relative stabilities of different isomers. Moreover, the basis set superposition error correction is critical to the interaction energies in the polar solvents. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Long‐range π‐type hydrogen bond in dimers CH2?HF,CH2O?H2O,and CH2O?NH3

Using four basis bets, (6‐311G(d,p), 6‐31+G(d,p), 6‐31++G(2d,2p), and 6‐311++G(3df,3pd), the optimized structures with all real frequencies were obtained at the MP2 level for the dimers CH₂O? HF, CH₂O? H₂O, CH₂O? NH₃, and CH₂O? CH₄. The structures of CH₂O? HF, CH₂O? H₂O, and CH₂O? NH₃ are cycle‐shaped, which result from the larger bend of σ‐type hydrogen bonds. The bend of σ‐type H‐bond O…H? Y (Y?F, O, N) is illustrated and interpreted by an attractive interaction of a chemically intuitive π‐type hydrogen bond. The π‐type hydrogen bond is the interaction between one of the H atoms of CH₂O and lone pair(s) on the F atom in HF, the O atom in H₂O, or the N atom in NH₃. In contrast with the above three dimers, for CH₂O? CH₄, because there is not a π‐type hydrogen bond to bend its linear hydrogen bond, the structure of CH₂O? CH₄ is noncyclic shaped. The interaction energy of hydrogen bonds and the π‐type H‐bond are calculated and discussed at the CCSD (T)/6‐311++G(3df,3pd) level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theoretical binding affinities and spectroscopy of complexes formed by cyclobis(paraquat-<Emphasis Type="Italic">p</Emphasis>-anthrancene) with some pharmaceutical molecules

Theoretical investigation on the stabilities and spectroscopic properties of the complexes formed by cyciobis(paraquat-p-anthracene) with pharmaceutical molecules were performed using the semi-empirical PM3 and B3LYP/3-21G methods. Based on the B3LYP/3-21G optimized geometries, the energies of the complexes were calculated at B3LYP/6-31G(d) level. The binding energies of the complexes were computed after the correction of basis set superposition error (BSSE). The energy gaps of the complexes are decreased due to the formation of the hydrogen bonds. The stretching vibrations of the C-H bonds adjacent to the hydrogen bonds in the IR spectra of the complexes calculated with PM3 method are red-shifted compared with those of the host. The chemical shifts of α-C and β-C atoms in the complexes calculated at B3LYP/3-21G level are shifted downfield due to the formation of the hydrogen bonds and the electron-withdrawing effect of the nitrogen atoms. The aromaticities of the complexes are improved because of the enlargement of the conjugation system and the overlap of electron cloud based on the nuclear independent chemical shifts (NICS) calculated at B3LYP/3-21G level.     

The thermodynamic properties of delta-lactones

The enthalpies of formation of δ-hexanolactone and δ-nonanolactone were determined by combustion calorimetry. Conformational analysis and quantum-chemical calculations of equilibrium structures, fundamental vibrations, moments of inertia, and total energies were performed for δ-pentanolactone (I), δ-hexanolactone (II), and δ-nonanolactone (III) by the B3LYP/6-311G(d,p), B3LYP/6-311++G(d,p), and G3MP2 methods. The experimental IR spectra and calculated vibrational frequencies were used to suggest the assignment of vibrational frequencies of stable conformations. The thermodynamic properties of I–III in the ideal gas state were determined over the temperature range 0–1500 K. A thermodynamic analysis of mutual isomerization in the gas and liquid phases over the temperature range 298.15–900 K and liquid-phase polymerization of γ- and δ-pentanolactones and 4-pentenoic acid over the temperature range 298.15–500 K was performed.     

Theoretical investigation on the structural and electronic properties of complexes formed by thymine and 2′-deoxythymidine with different anions

Hydrogen bonding interactions between thymine nucleobase and 2′-deoxythymidine nucleoside (dT) with some biological anions such as F⁻ (fluoride), Cl⁻ (chloride), OH⁻ (hydroxide), and NO₃ ⁻ (nitrate) have been explored theoretically. In this study, complexes have been studied by density functional theory (B3LYP method and 6-311++G (d,p) basis set). The relevant geometries, energies, and characteristics of hydrogen bonds (H-bonds) have been systematically investigated. There is a correlation between interaction energy and proton affinity for complexes of thymine nucleobase. The nature of all the interactions has been analyzed by means of the natural bonding orbital (NBO) and quantum theory atoms in molecules (QTAIM) approaches. Donors, acceptors, and orbital interaction energies were also calculated for the hydrogen bonds. Excellent correlations between structural parameter (δR) and electron density topological parameter (ρ _b) as well as between E(2) and ρ _b have been found. It is interesting that hydrogen bonds with anions can affect the geometry of thymine and 2′-deoxythymidine molecules. For example, these interactions can change the bond lengths in thymine nucleobase, the orientation of base unit with respect to sugar ring, the furanose ring puckering, and the C1′–N1 glycosidic linkage in dT nucleoside. Thus, it is necessary to obtain a fundamental understanding of chemical behavior of nucleobases and nucleosides in presence of anions.     

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