Conformational stability and intramolecular hydrogen bonding in 1,2-ethanediol and 1,4-butanediol

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.

Transferable Anisotropic Mie Potential Force Field for Alkanediols

The development of force fields for polyfunctional molecules, such as alkanediols, requires a careful account of different average intramolecular conformations for gas states compared to dense liquid states, where intra- and intermolecular hydrogen bonds compete. In the present work, the transferable anisotropic Mie (TAMie) potential is extended to 1,n-alkanediols. Using the convention that intramolecular nonbonded interactions up to and including the third neighbor are excluded, all force field parameters developed previously for 1-alcohols were transferred to 1,5-pentanediol and beyond, with good agreement with experimental phase equilibrium data. To obtain trans-gauche ratios of 1,2-ethanediol and 1,3-propanediol that are consistent with experimental results, the propensities for intra- and intermolecular hydrogen bonds had to be balanced. This was achieved by parameterizing the intramolecular dihedral energy functions governing the O-C-C-O and O-C-C-C angles while intramolecular charge-charge interactions were active. All partial charges belonging to a functional group are collected in a charge group and all interactions among two charge groups are evaluated even if they are separated by less than three bonds. With this approach, it is possible to apply the nonbonded parameters from 1-alcohols to alkanediols without further refinement. The agreement with experimental phase equilibrium and shear viscosity data is of similar quality as for the 1-alcohols and the trans-gauche ratio agrees with literature results from spectroscopic measurements and ab initio calculations.

Nuclear quantum effects in gas-phase 2-fluoroethanol

Torsional motions along the FCCO and HOCC dihedrals lead to the five unique conformations of 2-fluoroethanol, of which the conformer that is gauche along both dihedrals has the lowest energy. In this work, we explore how nuclear quantum effects (NQEs) manifest in the structural parameters of the lowest energy conformer, in the intramolecular free energy landscape along the FCCO and HOCC dihedrals, and also in the infrared spectrum of the title molecule, through the use of path integral simulations. We have first developed a full dimensional potential energy surface using the reaction surface Hamiltonian framework. On this potential, we have carried out path integral molecular dynamics simulations at several temperatures starting from the minimum energy well to explore structural influences of NQEs including geometrical markers of the interaction between the OH and F groups. From the computed free energy landscapes, significant reduction of the torsional barrier is found at low temperature near the cis region of the dihedrals, which can be understood through the trends in the radii of gyration of the atomic ring polymers. We find that the inclusion of NQEs in the computation of infrared spectrum is important to obtain good agreement with the experimental band positions.

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

Conformer Selection Upon Dilution with Water: The Fascinating Case of Liquid Ethylene Glycol Studied via Molecular Dynamics Simulations

The aqueous solution of ethylene glycol (EG) is a binary liquid mixture that displays rich conformational and structural behaviour, which has not yet been adequately explored through atomistic molecular dynamics simulations. Herein, employing an accurate force field for EG, several physical properties of this solution are calculated to be in quantitative agreement with experimental data. While 79 % of molecules in neat liquid EG exist with their central OCCO dihedral in the gauche state, this fraction increases to 89 % in the dilute aqueous solution, largely in response to the increase in the static dielectric constant of the solution from that of neat liquid EG. The increase in gauche conformers increases the mean dipole moment of EG molecules in the solution which is additionally contributed by specific conformational states of the two terminal HOCC dihedral angles.

Liquid Ethylene Glycol: Prediction of Physical Properties, Conformer Population and Interfacial Enrichment with a Refined Non-Polarizable Force Field

Periodic density functional theory based molecular dynamics simulations confirm the fraction of molecules in neat liquid ethylene glycol with their central OCCO dihedral in the trans conformation to be 21%,...

Subtle hydrogen bonds: benchmarking with OH stretching fundamentals of vicinal diols in the gas phase

The theoretical description of spectral signatures for weakly bound hydrogen contacts between alcohol groups is challenging and remains poorly characterised. By combining Raman jet spectroscopy with appropriately scaled harmonic DFT predictions and relaxation path analyses for 16 vicinal diols (ethylene glycol (ethane-1,2-diol), propane-1,2-diol, 3,3,3-trifluoro-propane-1,2-diol, rac-butane-2,3-diol, 2-methyl-propane-1,2-diol, 2-methyl-butane-2,3-diol, pinacol (2,3-dimethyl-butane-2,3-diol), 3-butene-1,2-diol, 1-phenyl-ethane-1,2-diol, trans-cyclobutane-1,2-diol, trans-cyclopentane-1,2-diol, trans-cyclohexane-1,2-diol, trans-cycloheptane-1,2-diol, cis-cyclohexane-1,2-diol, 1-(1-hydroxy-1-methylethyl)-cyclopentanol and [1,1'-bicyclopentyl]-1,1'-diol), 69 conformational assignments become possible in a two-tier approach with a 5 diol training and an 11 diol test set. The latter reveals systematic deviations for ring strain and secondary π interactions, but otherwise a remarkably robust correction and correlation model based on hybrid DFT with a minimally augmented triple-zeta basis set is obtained, whereas GGA functionals perform significantly worse. Raw experimental data in the 3560-3700 cm^-1 wavenumber range as well as computed geometries of all conformations invite further vibrational and structural benchmarking at the onset of hydrogen bonding. Beyond this diol-probed threshold, the accurate prediction of hydrogen bond induced shifts of different magnitudes remains one of the challenges for DFT functionals.

Intramolecular H-Bond Dynamics of Catechol Investigated by THz High-Resolution Spectroscopy of Its Low-Frequency Modes

Catechol is an oxygenated aromatic volatile organic compound and a biogenic precursor of secondary organic aerosols. Monitoring this compound in the gas phase is desirable due to its appreciable reactivity with tropospheric ozone. From a molecular point of view, this molecule is attractive since the two adjacent hydroxy groups can interchangeably act as donor and acceptor in an intramolecular hydrogen bonding due to the tunnelling between two symmetrically equivalent structures. Using synchrotron radiation, we recorded a rotationally-resolved Fourier Transform far-infrared (IR) spectrum of the torsional modes of the free and bonded -OH groups forming the intramolecular hydrogen bond. Additionally, the room temperature, pure rotational spectrum was measured in the 70–220 GHz frequency range using a millimeter-wave spectrometer. The assignment of these molecular transitions was assisted by anharmonic high-level quantum-chemical calculations. In particular, pure rotational lines belonging to the ground and the four lowest energy, vibrationally excited states were assigned. Splitting due to the tunnelling was resolved for the free -OH torsional state. A global fit combining the far-IR and millimeter-wave data provided the spectroscopic parameters of the low-energy far-IR modes, in particular those characterizing the intramolecular hydrogen bond dynamics.

Molecular Conformation and Hydrogen Bond Formation in Liquid Ethylene Glycol

The ethylene glycol (EG) molecule, HOCH₂CH₂OH, adopts a conformation where the central OCCO dihedral is exclusively gauche in the gaseous and crystalline states, but in the liquid state, for close to 20% of the molecules, the central OCCO adopts the energetically unfavorable trans conformation. Here we report calculations, based on ab initio molecular dynamics simulations, on the thermodynamics associated with hydrogen bond formation in the liquid state of EG between donor-acceptor pairs with different molecular conformations. We establish an operational, geometric definition of hydrogen bonds in liquid EG from an analysis of the proton NMR data and show that the key feature, irrespective of the conformation, is marked directionality with almost linear ∠HO···O angles. The free energy for hydrogen bond formation estimated as the potential of mean force for the reversible work associated with the passage from a hypothetical state where hydrogen bonding is absent and donor-acceptor pairs are randomly oriented to the hydrogen-bonded state where the pairs are oriented showed comparable magnitudes irrespective of the molecular conformation of either the donor or acceptor. The results suggest that the presence of the trans conformer in liquid EG would require an understanding of its role in the extended hydrogen-bonded network of the liquid.

Probing the Endo/Exo Isomeric Variants of a Binary Complex of Cyclopropylamine with p-Fluorophenol by LIF Spectroscopy: A Comparative Study with Ammonia Complex

Conformational isomers of an O-H···N hydrogen-bonded binary complex between para-fluorophenol (pFP) and a nonrigid primary amine base, cyclopropylamine (CPA), have been probed by means of laser-induced fluorescence (LIF) spectroscopy in a supersonic jet expansion. Two closely spaced electronic origin bands have been identified in the measured LIF excitation spectrum, and their assignments have been corroborated by making comparisons with the spectra of the parent pFP-NH₃ complex recorded under the same expansion condition. The observation is consistent with the presence of endo and exo isomeric variants of the complex predicted by electronic structure theory methods, and the endo isomer is stabilized by ∼2 kcal/mol additionally owing to the formation of a C-H···O and a C-H···π type of weak hydrogen bonds between the two moieties. In the fluorescence excitation (FE) spectrum, the low-frequency bands for different intermolecular modes gain substantial intensity, and this spectral feature is in contrast to that of the pFP-NH₃ complex. The Franck-Condon intensity of the bands has been simulated invoking Duschinsky rotation scheme, taking into consideration the ground- and excited-state geometries.

The reduced cohesion of homoconfigurational 1,2-diols

Homochiral encounters of vicinal diols are blocked from relaxing to the heterochiral global minimum dimer structure in supersonic jet expansions.

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.

How many solvent molecules are required to solvate chiral 1,2-diols with hydrogen bonding solvents? A VCD spectroscopic study

A VCD spectroscopic analysis of selected model systems for solute–solvent interactions of chiral diols with hydrogen bonding solvents DMSO and ACN.

Dealing with Hydrogen Bonding on the Conformational Preference of 1,3-Aminopropanols: Experimental and Molecular Dynamics Approaches

This study expands the knowledge on the conformational preference of 1,3-amino alcohols in the gas phase and in solution. By employing Fourier transform infrared spectroscopy, nuclear magnetic resonance (NMR) spectroscopy, density functional theory (DFT) calculations, quantum theory of atoms in molecules (QTAIM), natural bond orbital (NBO) analysis, and molecular dynamics (MD), the compounds 3-aminopropan-1-ol (1), 3-methylaminopropan-1-ol (2), and 3-dimethylaminopropan-1-ol (3) are evaluated. The results show that the most stable conformation of each compound in the gas phase and in nonpolar solvents exhibited an O-H···N intramolecular hydrogen bond (IHB). Based on the experimental and theoretical OH-stretching frequencies, the IHB becomes stronger from 1 to 3. In addition, from the experimental NMR J-couplings, the IHB conformers are predominant in nonbasic solvents, representing 70-80% of the conformational equilibrium, while in basic solvents, such conformers only represent 10%. DFT calculations and QTAIM analysis in the gas phase support the occurrence of IHBs in these compounds. The MD simulation indicates that the non-hydrogen-bonded conformers are the lowest energy conformations in the solution because of molecular interactions with the solvent, while they are absent in the implicit solvation model based on density. NBO analysis suggests that methyl groups attached on the nitrogen atom affect the charge transfer energy involved in the IHB. This effect occurs mostly because of a decrease in the s-character of the LPN orbital along with weakening of the charge transfer from LPN to σ*OH, which is caused by an increase in the C-C-N bond angle.

Cis/Trans Energetics in Epoxide, Thiirane, Aziridine and Phosphirane Containing Cyclopentanols: Effects of Intramolecular OH⋯O, S, N and P Contacts

A recent computational analysis of the stabilizing intramolecular OH⋯O contact in 1,2-dialkyl-2,3-epoxycyclopentanol diastereomers has been extended to thiiriane, aziridine and phosphirane analogues. Density functional theory (DFT), second-order Møller-Plesset perturbation theory (MP2) and CCSD(T) coupled-cluster computations with simple methyl and ethyl substituents indicate that electronic energies of the cis isomers are lowered by roughly 3 to 4 kcal mol⁻¹ when the OH group of these cyclopentanol systems forms an intramolecular contact with the O, S, N or P atom on the adjacent carbon. The results also suggest that S and P can participate in these stabilizing intramolecular interactions as effectively as O and N in constrained molecular environments. The stabilizing intramolecular OH⋯O, OH⋯S, OH⋯N and OH⋯P contacts also increase the covalent OH bond length and significantly decrease the OH stretching vibrational frequency in every system with shifts typically on the order of −41 cm⁻¹.

Anharmonic calculations of frequencies and intensities of OH stretching vibrations of (R)-1,3-butanediol conformers in the fundamentals and first overtones by density functional theory

The frequencies and absorption intensities of the five kinds of conformers of 1,3-butanediol with the same carbon skeleton (GG') were calculated by anharmonic calculation for the fundamentals and first overtones of OH stretching vibrations. The four kinds of conformers form intramolecular hydrogen bonds and one conformer did not. Intramolecular hydrogen bond formation shifted the frequency of fundamental and first overtone of H-bonding OH stretching vibration to the lower frequency. The absorption intensities of the fundamentals as well as the vibrational anharmonicities increased upon hydrogen bond formation, while the intensities of first overtones decreased. The differences of conformers were clearly seen in the frequencies of the first overtones of free OH.

Blue- and Red-Shifting Hydrogen Bonding: A Gas Phase FTIR and Ab Initio Study of RR'CO···DCCl3 and RR'S···DCCl3 Complexes

Blue-shifting H-bonded (C-D···O) complexes between CDCl₃ and CH₃HCO, (CH₃)₂CO, and C₂H₅(CH₃)CO, and red-shifting H-bonded (C-D···S) complexes between CDCl₃ with (CH₃)₂S and (C₂H₅)₂S have been identified by Fourier transform infrared spectroscopy in the gas phase at room temperature. With increasing partial pressure of the components, a new band appears in the C-D stretching region of the vibrational spectra. The intensity of this band decreases with an increase in temperature at constant pressure, which provides the basis for identification of the H-bonded bands in the spectrum. The C-D stretching frequency of CDCl₃ is blue-shifted by +7.1, +4, and +3.2 cm^-1 upon complexation with CH₃HCO, (CH₃)₂CO, and C₂H₅(CH₃)CO, respectively, and red-shifted by -14 and -19.2 cm^-1 upon complexation with (CH₃)₂S and (C₂H₅)₂S, respectively. By using quantum chemical calculations at the MP2/6-311++G** level, we predict the geometry, electronic structural parameters, binding energy, and spectral shift of H-bonded complexes between CDCl₃ and two series of compounds named RCOR' (H₂CO, CH₃HCO, (CH₃)₂CO, and C₂H₅(CH₃)CO) and RSR' (H₂S, CH₃HS, (CH₃)₂S, and (C₂H₅)₂S) series. The calculated and observed spectral shifts follow the same trends. With an increase in basicity of the H-bond acceptor, the C-D bond length increases, force constant decreases, and the frequency shifts to the red from the blue. The potential energy scans of the above complexes are done, which show that electrostatic attraction between electropositive D and electron-rich O/S causes bond elongation and red shift, and the electronic and nuclear repulsions lead to bond contraction and blue shifts. The dominance of the two opposing forces at the equilibrium geometry of the complex determines the nature of the shift, which changes both in magnitude and in direction with the basicity of the hydrogen-bond acceptor.

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.

Trends of intramolecular hydrogen bonding in substituted alcohols: a deeper investigation

NBO analyses were used to perform a deeper investigation on the effects that drive IAHB strength in substituted acyclic alcohols.