Structure and dynamics properties of liquid ethylene glycol from molecular dynamics simulations

Exploring the structural feature of water, alcohols, and their binary mixtures with concrete atomic charge assignments in Dreiding forcefield

The water, alcohols, and their binary mixtures are widely used in molecular simulations. However, the Dreiding force field lacks a generally accepted method for assigning atomic charges to these solvents during simulations. In this study, we propose a universal charge assignment for water and eight water-miscible alcohols in Dreiding. Through extensive molecular simulations, we demonstrate the good accuracy of our charge assignments in displaying characteristic of these solvents and their mixtures, including liquid density and structure. Moreover, we investigate equilibrium snapshot, radial distribution function, coordination number and hydrogen bonding, all of which confirm the miscibility of alcohols with water or ethanol. Notably, we reveal that the structure diversity among different mixtures can be attributed to distinctive characteristic of alcohols, including molecular volume, as well as the number and position of hydroxyls.

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.

Molecular Dynamics Study of the Effect of Charge and Glycosyl on Superoxide Anion Distribution near Lipid Membrane

To examine the effects of membrane charge, the electrolyte species and glycosyl on the distribution of negatively charged radical of superoxide anion (·O₂^-) around the cell membrane, different phospholipid bilayer systems containing ·O₂^- radicals, different electrolytes and phospholipid bilayers were constructed through Charmm-GUI and Amber16. These systems were equilibrated with molecular dynamics by using Gromacs 5.0.2 to analyze the statistical behaviors of ·O₂^- near the lipid membrane under different conditions. It was found that in the presence of potassium rather than sodium, the negative charge of the phospholipid membrane is more likely to rarefy the superoxide anion distribution near the membrane surface. Further, the presence of glycosyl significantly reduced the density of ·O₂^- near the phospholipid bilayer by 78.3% compared with that of the neutral lipid membrane, which may have a significant contribution to reducing the lipid peroxidation from decreasing the ·O₂^- density near the membrane.

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.

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.

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

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.