Ab initio studies of structural features not easily amenable to experiment

Nuclear quantum effects in gas-phase ethylene glycol

Path integral molecular simulations are used to explore the nuclear quantum effects (NQEs) on the structure, dihedral landscape and infrared spectrum of ethylene glycol. The simulations are carried out on a new reaction surface Hamiltonian-based model potential energy surface, with special focus on the role of the OCCO and HOCC dihedrals. In contrast with classical simulations, we analyse how the intramolecular interactions between the OH groups change due to zero-point effects as well as temperature. These are found to be weak. The NQEs on the free energy profile along the OCCO dihedral are analysed, where notable effects are seen at low temperatures and found to be correlated with the radii of gyration of the atoms. Finally, the power spectrum of the molecule from path integral simulations is compared with the experimental infrared spectrum, yielding a good agreement of band positions.

Computational analysis of the far infrared spectral region of various deuterated varieties of ethylene glycol

The far infrared spectra of three deuterated isotopologues of ethylene glycol, CH2OD-CH2OD, CH2OD-CH2OH and CH2OH-CH2OD, where the latter two species differ in their intramolecular hydrogen-bonding arrangement, are studied using highly correlated ab initio methods, vibrational second order perturbation theory and a variational procedure of reduced dimensionality. New subroutines suitable for the study of large systems with more than two interacting large amplitude motions were implemented and applied. The molecular symmetry of ethylene glycol decays by the formation of weak intramolecular bonds producing very asymmetrical stable structures. Three internal rotations contribute to the formation of a very anisotropic potential energy surface and to the puzzling distribution of the rovibrational energy levels. The ground vibrational state rotational constants and the centrifugal distortion constants (S-reduction, Ir representation) corresponding to the aGg' (G1) and gGg' (G2) forms are provided for the studied isotopologues. The low-lying vibrational levels up to 550 cm-1 are obtained variationally for J = 0. Two series of sublevels of the ground vibrational state are obtained: eight components localized in G1 lying between 0.0 and 0.3 cm-1 and eight sublevels localized in G2 lying between 138.1 and 138.4 cm-1. The gap between both sets is lower in CH2OD-CH2OD and more dispersed in the monodeuterated variety.

Computational Evidence Suggests That 1-Chloroethanol May Be an Intermediate in the Thermal Decomposition of 2-Chloroethanol into Acetaldehyde and HCl

The dehalogenation of 2-chloroethanol (2ClEtOH) in the gas phase with and without the participation of catalytic water molecules has been investigated using methods rooted into the density functional theory. The well-known HCl elimination leading to vinyl alcohol (VA) was compared to the alternative elimination route toward oxirane and shown to be kinetically and thermodynamically more favorable. However, the isomerization of VA to acetaldehyde in the gas phase, in the absence of water, was shown to be kinetically and thermodynamically less favorable than the recombination of VA and HCl to form the isomeric 1-chloroethanol (1ClEtOH) species. At the ωB97X-D/cc-pVTZ level of calculation, this species is more stable than 2ClEtOH by about 6 kcal mol^-1 at 298 K, and the reaction barrier for VA to 1ClEtOH is 23 kcal mol^-1 versus 55 kcal mol^-1 for the direct transformation of VA to acetaldehyde. In a successive step, 1ClEtOH can decompose directly to acetaldehyde and HCl with a lower barrier (29 kcal mol^-1) than that of VA to the same products (55 kcal mol^-1). The calculations were repeated using a single ancillary water molecule (W) in the complexes 2ClEtOH_W and 1ClEtOH_W. The latter adduct is now more stable than 2ClEtOH_W by about 8 kcal mol^-1 at 298 K, implying that the water molecule increased the already higher stability of 1ClEtOH in the gas phase. However, this catalytic water molecule lowers dramatically the barrier for the interconversion of VA to acetaldehyde (from 55 to 7 kcal mol^-1). This barrier is now smaller than the one for the conversion to 1ClEtOH (which also decreases, but not so much, from 23 to 13 kcal mol^-1). Thus, it is concluded that while 1ClEtOH may be a plausible intermediate in the gas phase dehalogenation of 2ClEtOH, it is unlikely that it plays a major role in water complexes (or, by inference, aqueous solution). It is also shown that neither in the gas phase nor in the cluster with one water molecule, the oxirane path is more favorable than the VA alcohol path. Additionally, a direct conversion of 2ClEtOH to 1ClEtOH through a transition state which resembles a VA molecule in a complex with a chlorine atom and a hydrogen atom on both sides of this planar species was found. This reaction path has also lower activation energy than the conversion to oxirane but not as low as the conversion to VA.

Weak intramolecular interaction effects on the torsional spectra of ethylene glycol, an astrophysical species

An elaborate variational procedure of reduced dimensionality based on explicitly correlated coupled clusters calculations is applied to understand the far infrared spectrum of ethylene-glycol, an astrophysical species. This molecule can be classified in the double molecular symmetry group G8 and displays nine stable conformers, gauche and trans. In the gauche region, the effect of the potential energy surface anisotropy due to the formation of intramolecular hydrogen bonds is relevant. For the primary conformer, stabilized by a hydrogen bond, the ground vibrational state rotational constants are computed to be A0 = 15 369.57 MHz, B0 = 5579.87 MHz, and C0 = 4610.02 MHz corresponding to differences of 6.3 MHz, 7.2 MHz, and 3.5 MHz from the experimental parameters. Ethylene glycol displays very low torsional energy levels whose classification is not straightforward and requires a detailed analysis of the torsional wavefunctions. Tunneling splittings are significant and unpredictable due to the anisotropy of the potential energy surface PES. The ground vibrational state splits into 16 sublevels separated ∼142 cm(-1). The splitting of the "G1 sublevels" was calculated to be ∼0.26 cm(-1) in very good agreement with the experimental data (0.2 cm(-1) = 6.95 MHz). Transitions corresponding to the three internal rotation modes allow assignment of previously observed Q branches. Band patterns, calculated between 362.3 cm(-1) and 375.2 cm(-1), 504 cm(-1) and 517 cm(-1), and 223.3 cm(-1) and 224.1 cm(-1), that correspond to the tunnelling components of the v21 fundamental (v21 = OH-torsional mode), are assigned to the prominent experimental Q branches.

1H NMR spectra of ethane-1,2-diol and other vicinal diols in benzene: GIAO/DFT shift calculations

The proton NMR spectra of several 1,2-diols in benzene have been analysed so as to associate each magnetically nonequivalent proton with its chemical shift. The shifts and coupling constants of the OH and methylene protons of ethane-1,2-diol have been determined in a wide range of solvents. The conformer distribution and the proton NMR shifts of these 1,2-diols in benzene have been computed on the basis of density functional theory. The solvent is included using the integral-equation-formalism polarizable continuum model implemented in Gaussian 09. Relative Gibbs energies for all stable conformers are calculated at the Perdew, Burke and Enzerhof (PBE)0/6-311 + G(d,p) level, and shifts are calculated using the gauge-including atomic orbital method with the PBE0/6-311 + G(d,p) geometry but using the cc-pVTZ basis set. Previous calculations on ethane-1,2-diol and propane-1,2-diol have been corrected and extended. New calculations on tert-butylethane-1,2-diol, phenylethane-1,2-diol, butane-2,3-diols (dl and meso) and cyclohexane-1,2-diols (cis and trans) are presented. Overall, the computed NMR shifts are in good agreement with experimental values for the OH protons but remain systematically high for CH protons. Some results based on the Gaussian 03 solvation model are included for comparison.

Additive empirical force field for hexopyranose monosaccharides

We present an all-atom additive empirical force field for the hexopyranose monosaccharide form of glucose and its diastereomers allose, altrose, galactose, gulose, idose, mannose, and talose. The model is developed to be consistent with the CHARMM all-atom biomolecular force fields, and the same parameters are used for all diastereomers, including both the alpha- and beta-anomers of each monosaccharide. The force field is developed in a hierarchical manner and reproduces the gas-phase and condensed-phase properties of small-molecule model compounds corresponding to fragments of pyranose monosaccharides. The resultant parameters are transferred to the full pyranose monosaccharides, and additional parameter development is done to achieve a complete hexopyranose monosaccharide force field. Parametrization target data include vibrational frequencies, crystal geometries, solute-water interaction energies, molecular volumes, heats of vaporization, and conformational energies, including those for over 1800 monosaccharide conformations at the MP2/cc-pVTZ//MP2/6-31G(d) level of theory. Although not targeted during parametrization, free energies of aqueous solvation for the model compounds compare favorably with experimental values. Also well-reproduced are monosaccharide crystal unit cell dimensions and ring pucker, densities of concentrated aqueous glucose systems, and the thermodynamic and dynamic properties of the exocyclic torsion in dilute aqueous systems. The new parameter set expands the CHARMM additive force field to allow for simulation of heterogeneous systems that include hexopyranose monosaccharides in addition to proteins, nucleic acids, and lipids.

Quantum Mechanical Analysis of 1,2-Ethanediol Conformational Energetics and Hydrogen Bonding

A proper understanding of the conformational energetics of 1,2-ethanediol (ethylene glycol) is important to the construction of molecular mechanics force fields for the treatment of carbohydrates since these biologically important molecules have a prevalence of vicinal hydroxyl groups. In the present study, quantum mechanical analysis of the 10 unique minimum-energy conformations of ethylene glycol is performed by using 10 model chemistries ranging from HF/6-311++G(d,p) up to a hybrid method that approximates CCSD(T)/cc-pVQZ. In addition, natural bond orbital (NBO) analysis of these conformations with deletion of pairings of CO bond/antibonding and lone pair/antibonding orbitals is used to investigate contributions from the "gauche" effect to ethylene glycol conformational energetics. MP2 with the "correlation consistent" basis sets and DFT/6-311++G(d,p) do the best job of matching the approximate CCSD(T)/cc-pVQZ energies while MP2/6-31G(d) and Hartree-Fock both fare poorly. NBO analysis shows the conformational energies to be independent of the deletion of matrix elements associated with (i) CO bonding and antibonding orbital interactions and (ii) lone pair and antibonding orbital interactions, whereas the energetic ordering correlates with geometric parameters consistent with internal hydrogen bonds. Thus, the present results suggest that standard molecular mechanics potential energy functional forms, which lack explicit terms to account for stereoelectronic effects, are appropriate for carbohydrates.

Weak intramolecular interactions in ethylene glycol identified by vapor phase OH-stretching overtone spectroscopy

Vapor phase OH-stretching overtone spectra of ethylene glycol were recorded to investigate weak intramolecular hydrogen bonding. The spectra were recorded with conventional absorption spectroscopy and laser photoacoustic spectroscopy in the first to fourth OH-stretching overtone regions. The room-temperature spectra are dominated by two conformers that show weak intramolecular hydrogen bonding. A less abundant third conformer, with no sign of hydrogen bonding, is also observed. Vapor phase spectra of the ethylene-d(4) glycol isotopomer were also recorded and used to identify an interfering resonance between CH-stretching and OH-stretching states in the fourth overtone. Anharmonic oscillator local mode calculations of the OH-stretching transitions have provided an accurate simulation of the observed spectra. The local mode parameters were calculated with coupled cluster ab initio methods. The calculations facilitate assignment of the different conformers in the spectra and illustrate the effect of the intramolecular hydrogen bonding.

On the Extent of Intramolecular Hydrogen Bonding in Gas-Phase and Hydrated 1,2-Ethanediol

We investigate the quantum dynamical nature of hydrogen bonding in 1,2-ethanediol and monohydrated 1,2-ethanediol using different levels of ab initio theory. Global full-dimensional potential energy surfaces were constructed from PW91/cc-pVDZ, B3LYP/cc-pVDZ, and MP2/cc-pVDZ ab initio data for gas-phase and monohydrated 1,2-ethanediol, using a modified Shepard interpolation scheme. Zero-point energies and nuclear vibrational wave functions were calculated on these surfaces using the quantum diffusion Monte Carlo algorithm. The nature of intra- and intermolecular hydrogen bonding in these molecules was investigated by considering a ground-state nuclear vibrational wavefunction with reduced complete nuclear permutation and inversion (CNPI) symmetry. Separate wavefunction histograms were determined from the ground-state nuclear vibrational wavefunction by projection into bondlength coordinates. The O-H and O-O wavefunction histograms and vibrationally averaged distances were then used to probe the extent of intra- and intermolecular hydrogen bonding. From these data, we conclude that gas-phase ethanediol may possess a weak hydrogen bond, with a relatively short O-O distance but no detectable proton delocalization. Monohydrated ethanediol was found to exhibit no intramolecular hydrogen bonding but instead possessed two intermolecular hydrogen bonds, indicated by both shortening of the O-O distance and significant proton delocalization. The degree of proton delocalization and shortening of the vibrationally averaged O-O distance was found to be dependent on the ab initio method used to generate the potential energy surface (PES) data set.

Alcohols, ethers, carbohydrates, and related compounds. III. The 1,2-dimethoxyethane system

Ethylene glycol, its dimethyl ether, and some related compounds have been studied using the MM4 molecular mechanics force field. The MM4 calculated structural and energetic results have been brought into satisfactory agreement with a considerable number of experimental data and MP2/6-311++G(2d,2p) ab initio calculations. The heats of formation of these compounds are also well calculated. The MM4 ethylene glycol conformations in particular are in good agreement, both geometrically and in terms of energy, with those from the ab initio calculations. The corresponding dimethyl ether is of special interest, because it has been suggested that the trans-gauche conformation is unusually stable due to the hydrogen bonding of a hydrogen on a methyl group with the more distant oxygen. It is shown in the present work that while this conformation is more stable than might have been expected, the energy is adequately calculated by MM4 without using any hydrogen bonding between the Cbond;H bond and the oxygen. If such hydrogen bonding occurs, it amounts to no more than about 0.5 kcal/mol in energy, and is too small to detect with certainty. Additionally, energetic relationships in trans-1,2-dimethoxycyclohexane, 1,3,5,7-tetraoxadecalin, and 3-methoxytetrahydropyran have been studied, and the calculated results are compared with experimental information, which is adequately reproduced.

Ab initio conformational studies on diols and binary diol-water systems using DFT methods. Intramolecular hydrogen bonding and 1:1 complex formation with water

Studies on the conformational equilibrium for the following diols, ethane-1,2-diol (12EG, CAS 107-21-1), 2R-D-(-)-propane-1,2-diol (12PG, CAS 4254-14-2), (2S,3S)-L-(+)-butane-2,3-diol (L23BD, CAS 19132-06-0), and (2S,3R)-meso-butane-2,3-diol (m23BD, CAS 5341-95-7), are described using Gaussian ab initio calculations involving density functional theory (DFT) methods. We also report in this article results on the stability and conformation for the 1:1 water-diol complex formed by ethane-1,2-diol, propane-1,2-diol, and L- and meso-butane-2,3-diol. The relative stability of the intramolecular (internal) hydrogen bond in a range of diols (n = 2 to 6), based on ab initio geometry optimization and determination of the -O...H- distance, dOH, and -O-H...O- angle, theta, increases through the sequence 1,2 approximately equals 2,3 < 1,3 < 1,4 approximately equals 1,5 approximately equals 1,6, as judged from the bond linearity and -O...H- separation. Quantum mechanical and topological analysis of possible intramolecular hydrogen bonding in this complete series of diols provides convincing evidence for this in diols in which the hydroxyl groups are separated by three or more carbon atoms, that is, in (n, n+m) diols for m > or = 2, but not for ethane-1,2-diol or other vicinal diols, which do not satisfy Popelier's topological and electron density criteria based on the AIM theory of Bader. Based on these criteria it is unlikely that vicinal diols are in fact capable of forming an intramolecular hydrogen bond, in spite of geometric and spectroscopic data in the literature suggesting otherwise.

The Rotational-Torsional Spectrum of the g'Gg Conformer of Ethylene Glycol: Elucidation of an Unusual Tunneling Path

The microwave spectrum of the energetically unfavored g'Gg conformer of ethylene glycol (CH(2)OH&bond;CH(2)OH) is reported. This spectrum is dominated by an interconversion geared-type large-amplitude motion during which each OH group in turn forms the intramolecular hydrogen bond. The microwave spectrum has been analyzed with the help of a Watson-type Hamiltonian plus a 1.4-GHz tunneling splitting. The rotational dependence of this tunneling splitting has been examined using an IAM approach and this yielded qualitative information on the tunneling path the molecule uses to interconvert between its two most stable conformers. Unexpectedly, but in agreement with ab initio calculations, when tunneling occurs between the energetically equivalent g'Gg and gGg' conformers, the OH groups are rotated stepwise through 240 degrees in the sense of a flip-flop rather than a concerted rotation and the molecule goes through the more stable g'Ga and aGg' forms. The electronic reasons for preferring a long rather than a short rotational path via a gGg form are discussed using calculated adiabatic vibrational modes. Copyright 2001 Academic Press.