Crystal structures of ethylene glycol and ethylene glycol monohydrate

Unlocking a new hydrogen-bonding marker: C-O bond shortening in vicinal diols revealed by rotational spectroscopy

The conformational space of cis-1,2-cyclohexanediol, a model molecule for cyclic vicinal diols, was investigated using rotational spectroscopy and density functional theory calculations. Four low energy conformers within an energy window of 5 kJ mol-1 were identified computationally. A rotational spectrum of jet-cooled cis-1,2-cyclohexanediol was recorded with a chirped pulse Fourier transform microwave spectrometer. Two sets of rotational transitions were observed and could be assigned to conformers of cis-1,2-cyclohexanediol. The non-observation of other low energy conformers was explained by conformational conversion barrier height calculations and results from experimental spectra recorded with different carrier gases. Eight isotopologues, including those with 13C and 18O, of the lowest energy conformer were observed, allowing the determination of the semi-experimental equilibrium structure, reSE. Interestingly, the structural analysis revealed that the C-O bond length of the intramolecular hydrogen-bond donor is shorter than that of the acceptor. This appears to be a general characteristic of vicinal diols and can be used as a novel hydrogen-bond marker in such compounds.

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.

Spiral Square Nanosheets Assembled from Ru Clusters

Spiral two-dimensional (2D) nanosheets exhibit unique physical and chemical phenomena due to their twisted structures. While self-assembly of clusters is an ideal strategy to form hierarchical 2D structures, it is challenging to form spiral nanosheets. Herein, we first report a screw dislocation involved assembled method to obtain 2D spiral cluster assembled nanosheets (CANs) with uniform square morphology. The 2D spiral Ru CANs with a length of approximately 4 μm and thickness of 20.7 ± 3.0 nm per layer were prepared via the assembly of 1-2 nm Ru clusters in the presence of molten block copolymer Pluronic F127. Cryo-electron microscopy (cryo-EM) and high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) demonstrate the existence of screw dislocation in the spiral assembled structure. The X-ray absorption fine structure spectrum indicates that Ru clusters are Ru³⁺ species, and Ru atoms are mainly coordinated with Cl with a coordination number of 6.5. Fourier-transform infrared (FT-IR) spectra and solid-state nuclear magnetic resonance hydrogen spectra (¹H NMR) indicate that the assembly process of Ru clusters is formed by noncovalent interactions, including hydrogen bonding and hydrophilic interactions. Additionally, the Ru-F127 CANs exhibit excellent photothermal conversion performance in the near-infrared (NIR) region.

Hydrogen-Bond Dynamics and Water Structure in Aqueous Ethylene Glycol Solution via Two-Dimensional Raman Correlation Spectroscopy

The hydrogen-bond (H-bond) dynamics and water structural transitions in aqueous ethylene glycol (EG) solution were investigated on the basis of concentration- and temperature-dependent two-dimensional Raman correlation spectroscopy (2D Raman-COS). At room temperature, EG-induced enhancement of the water structure when the EG/water molar ratio is less than 1:28 resulted from the hydrophobic effect around the methylene groups of EG. The decrease in the temperature caused an enhancement of the Raman peak at about 3200 cm^-1, representing an increase in the orderliness of water molecules. Further analysis of the water-specific structures by 2D Raman-COS reveals that the strong H-bond structure preferentially responds to external perturbations and induces a weak H-bond structural transition in water. Finally, EG-induced water structural transitions were calculated by the density functional theory (DFT). Hopefully, 2D Raman-COS combined with DFT calculations would advance the study of solute-induced water structural transitions in water-organic chemistry.

Raman structural study of ethylene glycol and 1,3-propylene glycol aqueous solutions

Raman spectra of ethylene glycol (EG) and 1,3-propylene glycol (1,3-PG) aqueous solutions with the diol content from 10 to 90 mol% were measured. The diol content weakly influences the EG and 1,3-PG Raman bands in the spectra of the solutions in the region 250-1800 cm^-1. This fact means that the conformational compositions of both the diols do not change significantly with dissolving in water. The intensity of the OH stretching band with respect to the diol bands intensities is the linear function of the ratio of the mole contents of water and the diol in the solutions. The spectral region 2800-3800 cm^-1 can be used to evaluate the chemical composition of these binary solutions. DFT modeling of the Raman spectra of EG molecule in water shell confirms the prevalence of the gauche-conformation of EG in the aqueous solutions.

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%,...

Role of the OH-group in the adsorption properties of methanol, ethanol, and ethylene glycol on 15-atom 3d, 4d, and 5d transition-metal clusters

The adsorption of alcohols on transition-metal (TM) substrates has received the attention of many researchers due to the applications of alcohols in several technological fields. However, our atomic-level understanding is still far from satisfactory, in particular for the interaction of alcohols with finite-size TM clusters, where new effects can arise due to the presence of quantum-size effects. In this work, we report a theoretical investigation of the adsorption properties of methanol, ethanol, and ethylene glycol on 12 different 3d, 4d, and 5d TM₁₅ clusters based on density functional theory calculations within the semi-empirical D3 van der Waals corrections. From the correlation analysis of all the lowest- and high-energy configurations, we identified the adsorption modes of methanol, ethanol, and ethylene glycol on the TM₁₅ clusters, in which the OH group binds to the cationic TM sites via the O-TM and H-TM interactions. Due to the relatively weak alcohol-TM₁₅ interaction, the changes induced on the TM₁₅ clusters are small, except for Au₁₅ and Ru₁₅, where the bare cluster changes its structure to a nearby minimum in the potential energy surface. The adsorption energy for the alcohol/TM₁₅ systems is correlated to the combination of several parameters, in which the main contribution is connected with the O-TM interaction and the HOTM angles. Furthermore, the TM electronegativity is an important descriptor for the methanol and ethanol adsorption energies, while charge transfer is important for ethylene glycol.

New Insights into the Structure of Glycols and Derivatives: A Comparative X-Ray Diffraction, Raman and Molecular Dynamics Study of Ethane-1,2-Diol, 2-Methoxyethan-1-ol and 1,2-Dimethoxy Ethane

In this study, we report a detailed experimental and theoretical investigation of three glycol derivatives, namely ethane-1,2-diol, 2-methoxyethan-1-ol and 1,2-dimethoxy ethane. For the first time, the X-ray spectra of the latter two liquids was measured at room temperature, and they were compared with the newly measured spectrum of ethane-1,2-diol. The experimental diffraction patterns were interpreted very satisfactorily with molecular dynamics calculations, and suggest that in liquid ethane-1,2-diol most molecules are found in gauche conformation, with intramolecular hydrogen bonds between the two hydroxyl groups. Intramolecular H-bonds are established in the mono-alkylated diol, but the interaction is weaker. The EDXD study also evidences strong intermolecular hydrogen-bond interactions, with short O···O correlations in both systems, while longer methyl-methyl interactions are found in 1,2-dimethoxy ethane. X-ray studies are complemented by micro Raman investigations at room temperature and at 80 °C, that confirm the conformational analysis predicted by X-ray experiments and simulations.

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.

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 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.

Effects of monohydric alcohols and polyols on the thermal stability of a protein

The thermal stability of a protein is lowered by the addition of a monohydric alcohol, and this effect becomes larger as the size of hydrophobic group in an alcohol molecule increases. By contrast, it is enhanced by the addition of a polyol possessing two or more hydroxyl groups per molecule, and this effect becomes larger as the number of hydroxyl groups increases. Here, we show that all of these experimental observations can be reproduced even in a quantitative sense by rigid-body models focused on the entropic effect originating from the translational displacement of solvent molecules. The solvent is either pure water or water-cosolvent solution. Three monohydric alcohols and five polyols are considered as cosolvents. In the rigid-body models, a protein is a fused hard spheres accounting for the polyatomic structure in the atomic detail, and the solvent is formed by hard spheres or a binary mixture of hard spheres with different diameters. The effective diameter of cosolvent molecules and the packing fractions of water and cosolvent, which are crucially important parameters, are carefully estimated using the experimental data of properties such as the density of solid crystal of cosolvent, parameters in the pertinent cosolvent-cosolvent interaction potential, and density of water-cosolvent solution. We employ the morphometric approach combined with the integral equation theory, which is best suited to the physical interpretation of the calculation result. It is argued that the degree of solvent crowding in the bulk is the key factor. When it is made more serious by the cosolvent addition, the solvent-entropy gain upon protein folding is magnified, leading to the enhanced thermal stability. When it is made less serious, the opposite is true. The mechanism of the effects of monohydric alcohols and polyols is physically the same as that of sugars. However, when the rigid-body models are employed for the effect of urea, its addition is predicted to enhance the thermal stability, which conflicts with the experimental fact. We then propose, as two essential factors, not only the solvent-entropy gain but also the loss of protein-solvent interaction energy upon protein folding. The competition of changes in these two factors induced by the cosolvent addition determines the thermal-stability change.

On the perturbation of the H-bonding interaction in ethylene glycol clusters upon hydration

Ab initio and density functional methods have been employed to study the structure, stability, and spectral properties of various ethylene glycol (EG(m)) and ethylene glycol-water (EG(m)W(n)) (m = 1-3, n = 1-4) clusters. The effective fragment potential (EFP) approach was used to explore various possible EG(m)W(n) clusters. Calculated interaction energies of EG(m)W(n) clusters confirm that the hydrogen-bonding interaction between EG molecules is perturbed by the presence of water molecules and vice versa. Further, energy decomposition analysis shows that both electrostatic and polarization interactions predominantly contribute to the stability of these clusters. It was found from the same analysis that ethylene glycol-water interaction is predominant over the ethylene glycol-ethylene glycol and water-water interactions. Overall, the results clearly illustrate that the presence of water disrupts the ethylene glycol-ethylene glycol hydrogen bonds.