Structures and energies of D-galactose and galabiose conformers as calculated by ab initio and semiempirical methods

Optimized geometries and total energies of some conformers of alpha- and beta-D-galactose have been calculated using the RHF/6-31G* ab initio method. Vibrational frequencies were computed at the 6-31G* level for the conformers that favor internal hydrogen bonding, in order to evaluate their enthalpies, entropies, Gibbs free energies, and then their structural stabilities. The semiempirical AM1, PM3, MNDO methods have also been performed on the conformers GG, GT, and TG of alpha- and beta-D-galactose. In order to test the reliability of each semiempirical method, the obtained structures and energies from the AM1, PM3, and MNDO methods have been compared to those achieved using the RHF/6-31G* ab initio method. The MNDO method has not been investigated further, because of the large deviation in the structural parameters compared with those obtained by the ab initio method for the galactose. The semiempirical method that has yielded the best results is AM1, and it has been chosen to perform structural and energy calculations on the galabiose molecule (the disaccharides constituted by two galactose units alpha 1,4 linked). The goal of such calculations is to draw the energy surface maps for this disaccharide. To realize each map, 144 different possible conformations resulting from the rotations of the two torsional angles psi and phi of the glycosidic linkage are considered. In each calculation, at each increment of psi and phi, using a step of 30 degrees from 0 to 330 degrees, the energy optimization is employed. In this article, we report also calculations concerning the galabiose molecule using different ab initio levels such as RHF/6-31G*, RHF/6-31G**, and B3Lyp/6-31G*.     

Ab initio Study on Ionization Energies of 3-Amino-1-propanol

Fourteen conformers of 3-amino-1-propanol as the minima on the potential energy surface are examined at the MP2/6-311++G** level. Their relative energies calculated at B3LYP, MP3 and MP4 levels of theory indicated that two most stable conformers display the in-tramolecular OH…N hydrogen bonds. The vertical ionization energies of these conformers calculated with ab initio electron propagator theory in the P3/aug-cc-pVTZ approximation are in agreement with experimental data from photoelectron spectroscopy. Natural bond orbital analyses were used to explain the differences of IEs of the highest occupied molec-ular ortibal of conformers. Combined with statistical mechanics principles, conformational distributions at various temperatures are obtained and the temperature dependence of pho-toelectron spectra is interpreted.     

Structure and dynamics of the aniline-argon complex as derived from its potential energy surface

The structure and dynamics of the van der Waals (vdW) complex of aniline (An) with argon (Ar) are studied using ab initio methods. The inversion potential of the aniline-argon (AnAr) complex perturbed by the weak vdW interaction is calculated taking into account subtle corrections from the zero-point energy of the vdW modes and from the frequency shifts of the An normal modes modified by the complexation. The intermolecular potential energy surface (PES) of the AnAr complex is determined by performing a large-scale computation of the interaction energy and the fitting of the analytical many-body expansion to the set of single-point interaction energies. The PES determined shows two deep local minima corresponding to the anti and syn AnAr conformers. The difference in the energies of these two minima is only 15 cm-1, but it is sufficient to localize the inversion wave functions and to form the two conformers. In the conformer anti (syn) of lower (higher) energy, Ar is bound to the An ring opposite (adjacent) the amino-hydrogens. In the additional local minima higher in energy, Ar approaches the aniline ring between the C-H bonds near its plane. An additional local minimum is located opposite of nitrogen between the two N-H bonds. The high-energy minima are, however, too flat to form stable conformers. The perturbation of the interaction of Ar with the phenyl ring by the NH2 group is described by the vdW hole, which is responsible for unusually strong intermode mixing in the excited intermolecular vibrational states. The analysis of these states calculated for the ground (S0) as well as the first excited electronic state (S1) resolves difficulties faced earlier with the assignment of the observed vibronic bands of AnAr.     

On the conformation of the propranolol molecule

The structure of the propranolol molecule has been optimized within the AM1 and PM3 semiempirical framework followed by ab initio HF/6-31G^* refinement. On each calculation level the conformational space was sampled to search for the lowest-energy conformer(s) from among a few hundreds of conformers at the semiempirical step and next from among a few dozens of conformers at the ab initio level. Finally, five stable conformers were found; each stabilized by one or two of the three possible hydrogen bonds. The geometrical and electronic parameters were established and found to differ only slightly in the structures with the hydrogen bond either present or not.     

Molecular structure and conformational stability of allylisocyanate—post-Hartree–Fock and density functional theory studies

The molecular structure and conformational stability of allylisocyanate (CH₂CHCH₂NCO) molecule was studied using the ab initio and DFT methods. The geometries of possible conformers, C-gauche (δ=120°, θ=0°) (δ=C=C–C–N and θ=C–C–N=C) and C-cis N-trans (δ=0° and θ=180°) were optimized employing HF/6-31G^*, MP2/6-31G^* levels of theory of ab initio and BLYP, B3LYP, BPW91 and B3PW91 methods of DFT implementing the atomic basis set 6-311+G(d,p). The structural and physical parameters of the above conformers were discussed with the experimental and theoretical values of the related molecules, methylisocyanate and 3-fluoropropene. It has been found that the N=C=O bond angle is not linear as the experimental result for both the conformers and the theoretical bond angle is 173°. The rotational potential energy surfaces have been performed at the HF/6-31G^*, and MP2/6-31G^* levels of theory. The Fourier decomposition potentials were analysed at the HF/6-31G^*, and MP2/6-31G^* levels of theory. The HF/6-31G^* level of theory predicted that the C-gauche conformer is more stable than the C-cis N-trans conformer by 0.41 kJ/mol, but the MP2 and DFT methods predicted the C-cis N-trans conformer is found to be more stable than the C-gauche conformer. The calculated chemical hardness value at the HF/6-31G^* level of theory predicted the C-cis N-trans form is more stable than C-gauche form, whereas the chemical hardness value at the MP2/6-31G^* level of theory favours the slight preference towards the C-gauge conformer.     

Ab initio conformational maps in the gas phase and aqueous solution for a prototype of the glycosidic linkage

Ab initio conformational maps for methoxyethoxymethane (MEM) in both the gas phase and aqueous solution have been constructed using two different approaches. The results obtained allow us to conclude that a rigid conformational map is able to predict the regions of the minima, in the potential energy surface of MEM, in full agreement with those found in the relaxed conformational map, in both phases studied. This is a good indication that ab initio rigid conformational maps may be reliably used to sort the stablest conformers of disaccharides in aqueous solution. Besides that, in the MEM case, the solvation effects do not give rise to any new local minimum in its potential energy surface, but just change the relative energies of the stablest conformers found in the gas phase. This may be an indication that even in aqueous solution the anomeric effect is still the determinant effect defining the conformation of the molecule.Proceedings of the 11th International Congress of Quantum Chemistry satellite meeting in honor of Jean-Louis Rivail     

Molecular structures of the two most stable conformers of free glycine

The equilibrium molecular structures of the two lowest-energy conformers of glycine, Gly-Ip and Gly-IIn, have been characterized by high-level ab initio electronic structure computations, including all-electron cc-pVTZ CCSD(T) geometry optimizations and 6-31G* MP2 quartic force fields, the latter to account for anharmonic zero-point vibrational effects to isotopologic rotational constants. Based on experimentally measured vibrationally averaged effective rotational constant sets of several isotopologues and our ab initio data for structural constraints and zero-point vibrational shifts, least-squares structural refinements were performed to determine improved Born-Oppenheimer equilibrium (r(e)) structures of Gly-Ip and Gly-IIn. Without the ab initio constraints even the extensive set of empirical rotational constants available for 5 and 10 isotopologues of Gly-Ip and Gly-IIn, respectively, cannot satisfactorily fix their molecular structure. Excellent agreement between theory and experiment is found for the rotational constants of both conformers, the rms residual of the final fits being 7.8 and 51.6 kHz for Gly-Ip and Gly-IIn, respectively. High-level ab initio computations with focal point extrapolations determine the barrier to planarity separating Gly-IIp and Gly-IIn to be 20.5 +/- 5.0 cm(-1). The equilibrium torsion angle tau(NCCO) of Gly-IIn, characterizing the deviation of its heavy-atom framework from planarity, is (11 +/- 2) degrees. Nevertheless, in the ground vibrational state the effective structure of Gly-IIn has a plane of symmetry.     

Chiral self-recognition: direct spectroscopic detection of the homochiral and heterochiral dimers of propylene oxide in the gas phase

In this report, we describe rotational spectroscopic and high-level ab initio studies of the 1:1 chiral molecular adduct of propylene oxide dimer. The complexes are bound by weak secondary hydrogen bonds, that is, the O(epoxy)...H-C noncovalent interactions. Six homochiral and six heterochiral conformers were predicted to be the most stable configurations where each monomer acts as a proton acceptor and a donor simultaneously, forming two six- or five-membered intermolecular hydrogen-bonded rings. Rotational spectra of six, that is, three homochiral and heterochiral conformer pairs, out of the eight conformers that were predicted to have sufficiently large permanent electric dipole moments were measured and analyzed. The relative conformational stability order and the signs of the chiral recognition energies of the six conformers were determined experimentally and were compared to the ab initio computational results. The experimental observations and the ab initio calculations suggest that the concerted effort of these weak secondary hydrogen bonds can successfully lock the subunits in a particular orientation and that the overall binding strength is comparable to a classic hydrogen bond.     

Gaseous arginine conformers and their unique intramolecular interactions

Extensive ab initio calculations were employed to characterize stable conformers of gaseous arginine, both the canonical and zwitterionic tautomers. Step-by-step geometry optimizations of possible single-bond rotamers at the B3LYP/6-31G(d), B3LYP/6-31++G(d,p), and MP2/6-31++G(d,p) levels yield numerous structures that are more stable than any known ones. The final electronic energies of the conformers were determined at the CCSD/6-31++G(d,p) level. The lowest energies of the canonical and zwitterionic structures are lower than the existing values by 2.0 and 2.3 kcal/mol, respectively. The relative energies, rotational constants, dipole moments, and harmonic frequencies of the stable conformers remain for future experimental verification. The conformational distributions at various temperatures, estimated according to thermodynamic principles, consist almost exclusively of the newly found structures. One striking feature is the occurrence of blue-shifting hydrogen bonds in all six of the most stable conformers. A unique feature of important conformations is the coexistence of dihydrogen and blue- and red-shifting hydrogen bonds. In addition to the hydrogen bonds, the stereoelectronic effects were also found to be important stabilization factors. The calculated and measured proton affinities agree within the theoretical and experimental uncertainties, affirming the high quality of our conformational search. The theoretical gas-phase basicity of 245.9 kcal/mol is also in good agreement with the experimental value of 240.6 kcal/mol. The extensive searches establish firmly that gaseous arginine exists primarily in the canonical and not the zwitterionic form.     

Conformational analysis of primary ethylene ozonide by gradient and multiconfigurational scf calculations

Four conformers of ethylene primary ozonide have been studied by ab initio gradient and MC SCF calculations, using gaussian-type basis functions. The MC SCF results indicate that the conformers are not as close in energy as suggested from single-determinant SCF calculations. The oxygen-oxygen and carbon-oxygen half-chair structures are much lower in energy than the carbon-carbon half-chair.     

Structure,conformation and NMR studies on 1,2-dioxane and halogen substituted 1,2-dioxane molecules

Molecular structure and conformational stability of chair and twist conformers of 1,2-dioxane and halogen substituted compounds of the 1,2-dioxane have been studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 1,2-dioxane, 3,6-difluoro, 3,6-dichloro, 3,3,6,6-tetrafluoro and 3,3,6,6-tetrachloro 1,2-dioxane compounds were optimized at HF, MP2, B3LYP and B3PW91 levels of theory by implementing 6-31G* basis set. To study the effect of polar medium, self-consistent reaction field theory is used to optimize the conformers at B3LYP/6-31G* level of theory. The geometrical parameters of chair and twist conformers have been discussed in the light of interaction between lone pair electrons present in the oxygen and substituted halogen atoms. The relative stability of the conformers have been studied using relative energy, maximum hardness principle and thermodynamical quantities. The 13C-NMR chemical shift study for carbon atoms in the title compounds are calculated and the results have been discussed.     

Quantum chemical studies on structure, conformation and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene

The molecular structure, conformational stability and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene in gas phase have been investigated using ab initio and density functional theory methods. The molecular geometries and energetics of possible conformers were obtained by employing MP2, B3LYP and B3PW91 levels of theory implementing 6-31G* basis set. The relative stabilities of the conformations were evaluated from the energy differences of the structure. Chemical hardness (η) and chemical potential (μ) were calculated at HF/6-31G* level of theory for all the positional and geometrical isomers to study the maximum hardness principle. Each optimized structure has been tested against the imaginary frequencies at MP2/6-31G* level of theory in order to be sure they are located at energy minimum.     

Structure and conformation of cyclobutylsilane: an electron diffraction and ab initio study

A gas electron diffraction study of cyclobutylsilane results in a mixture of equatorial and axial conformers, with the equatorial confomer slightly more stable (Δ G = 0.8 ± 0.4 kJ mol⁻¹). The cyclobutyl ring is distorted with the adjacent bonds longer (C₁---C₂ = 1.573 (4) Å) than the opposite bonds (C₂---C₃ = 1.557 (4) Å). The experimental values for the energy difference between the two conformers and for the geometric parameters are reproduced very well by ab initio calculations. The importance of silicon 3d orbitals in the interpretation of ring distortion is ambiguous, but on the basis of the ab initio calculations the participation of silicon 3d functions is negligible.     

Ab initio rate constants from hyperspherical quantum scattering: application to H + CH4 --> H2 + CH3

A general and practical procedure is described for calculating rate constants for chemical reactions using a minimal number of ab initio calculations and quantum-dynamical computations. The method exploits a smooth interpolating functional developed in the hyperspherical representation. This functional is built from two Morse functions and depends on a relatively small number of parameters with respect to conventional functionals developed to date. Thus only a small number of ab initio points needs to be computed. The method is applied to the H + CH4 --> H2 + CH3 reaction. The quantum scattering calculations are performed treating explicitly the bonds being broken and formed. All the degrees of freedom except the breaking and forming bonds are optimized ab initio and harmonic vibrational frequencies and zero-point energies for them are calculated at the MP2(full) level with a cc-pVTZ basis set. Single point energies are calculated at a higher level of theory with the same basis set, namely CCSD(T, full). We report state-to-state cross sections and thermal rate constants for the title reaction and make comparisons with previous results. The calculated rate constants are in good agreement with experiments.     

Studies of intramolecular hydrogen bonding in protonated and non-protonated hexahydropyridinobenzodioxins

The conformational stability of hexahydropyridobenzodioxin and related derivatives in both protonated and non-protonated forms have been investigated by means of ab initio molecular orbital methods as well as semi-empirical AM1 and PM3 methods. One of the cis conformers (cis2e) has been found to be most stable due to the formation of an intramolecular hydrogen bond, other conformers including the trans isomer cannot form this interaction but are of different stability because of the orientation of the polar oxygens and the nitrogen. The effect of the intramolecular hydrogen bonding on the stability of hexahydropyridobenzodioxin and its methylated derivatives has been examined using various basis sets levels. In protonated form, both the semi-empirical and ab initio calculations give excellent agreement in energetic order; however, different orderings of conformer stabilities are observed by different computational methods in non-protonated form. The results provide insight into the intramolecular hydrogen bonding in computational studies of biologically important molecules.     

Ab initio and DFT studies of hydrogen bond interactions in difluoroacetic acid dimer

Density functional theory (DFT) and ab initio calculations were performed for difluoroacetic acid (DFA). Eight theoretically possible conformers were considered and by using conformational analysis only three stable conformers were found. The hydrogen bonding interaction of DFA complex has been investigated using DFT and ab initio methods for cis conformers. Stabilization energies of dimers including basis set superposition error and ZPE were found in the range 8.89–13.08 kcal mol⁻¹. It was found that EFC dimer is slightly more stable. Red shift of O–H bond in the range −226.3 to 505.7 cm⁻¹ predicted for dimers. The natural bond orbital analysis was applied to characterize nature of the interaction.     

Conformations of allyl amine: theory vs experiment

The relative stabilities of the five conformers of allyl amine, a medium-size aliphatic molecule, were estimated by applying ab initio quantum mechanical methods at several levels of theory. The second-order M?ller-Plesset perturbation method (MP2), quadratic configuration interaction including single and double excitations (QCISD), coupled-cluster with single and double excitations (CCSD) and CCSD plus perturbative triple excitations [CCSD(T)] were applied. The Dunning correlation consistent basis sets (through aug-cc-pVQZ and cc-pV5Z) were employed. The MP2 energies relative to the energy of the cis-trans conformer reported here appear to approach the basis set limit. The predicted allyl amine conformer energies approaching the Hartree-Fock basis set limit are 158 cm-1 (cis-gauche), -5 cm-1 (gauche-trans), and -146 cm-1 (gauche-gauche). The same three relative energies near the MP2 basis set limit are 135, 103, and 50 cm-1, respectively. The analogous energies deduced from experiment are 173 +/- 12, 92 +/- 8, and 122 +/- 5 cm-1. The theoretical results obtained in the present study suggest that satisfactory predictions of the conformer energetics of allyl amine may be achieved only by theoretical methods that incorporate consideration of correlation effects in conjunction with large basis sets. Evaluation of the zero-point vibrational energy corrections is critical, due to the very small classical energy differences between the five conformers of allyl amine. Agreement between theory and experiment for the gauche-gauche conformational energy remains problematical.     

Conformational behavior of beta-proline oligomers

Conformations of the monomer, dimer, and hexamer of beta-proline ((S) pyrrolidine-3-carboxylic acid) were determined using ab initio molecular orbital calculations at the RHF/6-31G level of theory. The calculated minima are in good agreement with experimental data for the system and imply that the conformations could be controlled through chemical modification at Calpha, Cgamma, or Cdelta. The monomer and dimer are small and flexible with many low-energy minima. In the hexamer, two forms of regular secondary structure are preferred: left-handed helices with cis-peptide bonds and right-handed helices with trans-peptide bonds. This is similar to the behavior of alpha-proline helices, except that the relationship between the peptide rotamer and the handedness of the helix is reversed. Therefore, helices of the enantiomer of beta-proline ((R)-pyrrolidine-3-carboxylic acid) should exhibit the same behavior as alpha-proline helices. Through understanding the conformational behavior of beta-proline in various environments, it may be possible to use these protein mimics to inhibit various protein-protein recognition events. To estimate these effects, SCRF energies for the conformers were determined in dielectrics corresponding to water, methanol, and chloroform. It appears that the cis helices are more favorably solvated than the trans helices, but the cause is not clear.     

Atomic partitioning of the dissociation energy of the P-O(H) bond in hydrogen phosphate anion (HPO4(2-)): disentangling the effect of Mg2+

This paper has three goals: (1) to provide a first step in understanding the atomic basis of the role of magnesium in facilitating the dissociation of the P-O bond in phosphorylated biochemical fuel molecules (such as ATP or GTP), (2) to compare second-order M?ller-Plesset perturbation theory (MP2) results with those obtained at the more economical density functional theory (DFT) level for a future study of larger more realistic models of ATP/GTP, and (3) to examine the calculation of atomic total energies from atomic kinetic energies within a Kohn-Sham implemention of DFT, as compared to ab initio methods. A newly described method based on the quantum theory of atoms in molecules (QTAIM), which is termed the "atomic partitioning of the bond dissociation energy" (APBDE), is applied to a simple model of phosphorylated biological molecules (HPO42-). The APBDE approach is applied in the presence and in the absence of magnesium. It is found that the P-O(H) bond in the magnesium complex is shorter, exhibits a higher stretching frequency, and has a higher electron density at the bond critical point than in the magnesium-free hydrogen phosphate anion. Though these data would seem to suggest a stronger P-O(H) bond in the magnesium complex compared to the magnesium-free case, the homolytic breaking of the P-O(H) bond in the complex is found to be easier, i.e., has a lower BDE. This effect is the result of the balance of several atomic contributions to the BDE induced by the magnesium cation, which stabilizes the dissociation product more than it stabilizes the intact model molecule.     

Theoretical studies on the hydration of formic acid by ab initio and ABEEMσπ fluctuating charge model

The interaction between formic acid (FA) and water was systemically investigated by atom-bond electronegativity equalization method fused into molecular mechanics (ABEEMσπ/MM) and ab initio methods. The geometries of 20 formic acid–water complexes (FA–water) were obtained using B3LYP/aug-cc-pVTZ level optimizations, and the energies were determined at the MP2/aug-cc-pVTZ level with basis set superposition error (BSSE) and zero-point vibrational energy (ZPVE) corrections. The ABEEMσπ potential model gives reasonable properties of these clusters when compared with the present ab initio data. For interaction energies, the root mean square deviation is 0.74 kcal/mol, and the linear coefficient reaches 0.993. Next, FA in aqueous solution was also studied. The hydrogen-bonding pattern due to the interactions with water has been analyzed in detail. Furthermore, the ABEEMσπ charges changed when H₂O interacted with the FA molecule, especially at the sites where the hydrogen bonds form. These results show that the ABEEMσπ fluctuating charge model is fine giving the overall characteristic hydration properties of FA–water systems in good agreement with the high-level ab initio calculations.