Molecular dynamics simulation of liquid ethylene glycol and its aqueous solution

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.

Raman analysis of aqueous solutions of ethylene glycol and 1,3-propylene glycol: Fundamental and applied aspects

Aqueous solutions of glycols, on the one hand, are widely used in many applications; on the other hand, they can serve as simple and representative models for studying intra- and intermolecular hydrogen bonds. In this work, we analyze the possibilities and limitations of Raman spectroscopy for fundamental and applied researches of such solutions on the examples of ethylene glycol (EG) and 1,3-propylene glycol (1,3-PG). It is shown that Raman spectroscopy is an effective tool for monitoring temporal changes in the structure of glycol solutions deposited on substrates. This study demonstrates that the water content in the solutions on the substrates decreases rapidly with time, and the rate of this decrease depends on the chemical structure of both glycol and substrate. It was found that the reduction in the water content leads to slight decrease in the contents of gauche-conformers in the backbones of EG and 1,3-PG molecules. It is shown that use of the 1064 nm excitation ensures a reliable Raman analysis of automotive antifreezes containing various dyes, in particular determination of the relative contents of water and glycol.

Response of Ionic Hydration Structure and Selective Transport Behavior to Aqueous Solution Chemistry during Nanofiltration

The effect of aqueous solution chemistry on the ionic hydration structure and its corresponding nanofiltration (NF) selectivity is a research gap concerning ion-selective transport. In this study, the hydration distribution of two typical monovalent anions (Cl- and NO3-) under different aqueous solution chemical conditions and the corresponding transmembrane selectivity during NF were investigated by using in situ liquid time-of-flight secondary ion mass spectrometry in combination with molecular dynamics simulations. We demonstrate the inextricable link between the ion hydration structure and the pore steric effect and further find that ionic transmembrane transport can be regulated by breaking the balance between the hydrogen bond network (i.e., water-water) and ion hydration (i.e., ion-water) interactions of hydrated ion. For strongly hydrated (H2O)nCl- with more intense ion-water interactions, a higher salt concentration and coexisting ion competition led to a larger hydrated size and, thus, a higher ion rejection by the NF membrane, whereas weakly hydrated (H2O)nNO3- takes the reverse under the same conditions. Stronger OH--anion hydration competition resulted in a smaller hydrated size of (H2O)nCl- and (H2O)nNO3-, showing a lower observed average hydration number at pH 10.5. This study deepens the long-overlooked understanding of NF separation mechanisms, concerning the hydration structure.

Photoisomerization of Spiropyran in Solution: A Surface Hopping Investigation

We performed a computational investigation of the photoisomerization of an unsubstituted spiropyran in solution using surface hopping molecular dynamics simulations in a semiempirical framework. To bring out the solvent effects on the excited state dynamics, we have run simulations in three different environments: chloroform, methanol, and ethylene glycol. We show that the interaction with a moderately polar solvent such as chloroform has little effect on the dynamics when compared to vacuum results previously obtained by our group. At variance, the interaction with the protic solvents considered considerably affects the reaction mechanism, the quantum yield, and the excited state lifetimes.

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.

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.

Structural and Dynamical Properties of Liquids in Confinements: A Review of Molecular Dynamics Simulation Studies

Molecular dynamics (MD) simulations are a powerful tool for detailed studies of altered properties of liquids in confinement, in particular, of changed structures and dynamics. They allow, on one hand, for perfect control and systematic variation of the geometries and interactions inherent in confinement situations and, on the other hand, for type-selective and position-resolved analyses of a huge variety of structural and dynamical parameters. Here, we review MD simulation studies on various types of liquids and confinements. The main focus is confined aqueous systems, but also ionic liquids and polymer and silica melts are discussed. Results for confinements featuring different interactions, sizes, shapes, and rigidity will be presented. Special attention will be given to situations in which the confined liquid and the confining matrix consist of the same type of particles and, hence, disparate liquid-matrix interactions are absent. Findings for the magnitude and the range of wall effects on molecular positions and orientations and on molecular dynamics, including vibrational motion and structural relaxation, are reviewed. Moreover, their dependence on the parameters of the confinement and their relevance to theoretical approaches to the glass transition are addressed.

A Relation between the Formation of a Hydrogen-Bond Network and a Time-Scale Separation of Translation and Rotation in Molecular Liquids

We study the relation between the translational and rotational motions of liquids, which is anticipated in the framework of the Stokes-Einstein-Debye (SED) treatment. For this purpose, we exploit the fact that ¹H field-cycling nuclear magnetic resonance relaxometry and molecular dynamics simulations provide access to both modes of motion. The experimental and computational findings are fully consistent and show that the time-scale separation between translation and rotation increases from the van der Waals liquid o-terphenyl over ethylene glycol to the hydrogen-bonded liquid glycerol, indicating an increasing degree of breakdown of the SED relation. Furthermore, the simulation results for two ethylene glycol models with different molecular conformations indicate that the translation is more retarded than the rotation when the density of intermolecular hydrogen bonds increases. We conclude that an increasing connectivity of a hydrogen-bond network leads to an increasing time-scale separation and, thus, to a stronger SED violation.

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

Molecular Dynamics Simulations of Water, Silica, and Aqueous Mixtures in Bulk and Confinement

Aqueous systems are omnipresent in nature and technology. They show complex behaviors, which often originate in the existence of hydrogen-bond networks. Prominent examples are the anomalies of water and the non-ideal behaviors of aqueous solutions. The phenomenology becomes even richer when aqueous liquids are subject to confinement. To this day, many properties of water and its mixtures, in particular, under confinement, are not understood. In recent years, molecular dynamics simulations developed into a powerful tool to improve our knowledge in this field. Here, our simulation results for water and aqueous mixtures in the bulk and in various confinements are reviewed and some new simulation data are added to improve our knowledge about the role of interfaces. Moreover, findings for water are compared with results for silica, exploiting that both systems form tetrahedral networks.

New Force Field Model for Propylene Glycol: Insight to Local Structure and Dynamics

In this work we developed a new force field model (FFM) for propylene glycol (PG) based on the OPLS all-atom potential. The OPLS potential was refined using quantum chemical calculations, taking into account the densities and self-diffusion coefficients. The validation of this new FFM was carried out based on a wide range of physicochemical properties, such as density, enthalpy of vaporization, self-diffusion coefficients, isothermal compressibility, surface tension, and shear viscosity. The molecular dynamics (MD) simulations were performed over a large range of temperatures (293.15-373.15 K). The comparison with other force field models, such as OPLS, CHARMM27, and GAFF, revealed a large improvement of the results, allowing a better agreement with experimental data. Specific structural properties (radial distribution functions, hydrogen bonding and spatial distribution functions) were then analyzed in order to support the adequacy of the proposed FFM. Pure propylene glycol forms a continuous phase, displaying no microstructures. It is shown that the developed FFM gives rise to suitable results not only for pure propylene glycol but also for mixtures by testing its behavior for a 50 mol % aqueous propylene glycol solution. Furthermore, it is demonstrated that the addition of water to the PG phase produces a homogeneous solution and that the hydration interactions prevail over the propylene glycol self-association interactions.

Lifshitz phase: the microscopic structure of aqueous and ethanol mixtures of 1,n-diols

We study binary mixtures of ethylene glycol and 1,3-propandiol with water or ethanol using computer simulations.

Hydrogen Bonding and Dielectric Spectra of Ethylene Glycol-Water Mixtures from Molecular Dynamics Simulations

Mixtures of ethylene glycol with water are a prominent example of media with variable viscosity. Classical molecular dynamics simulations at room temperature were performed for mixtures of ethylene glycol (EG) and water with EG mole fractions of xE = 0.0, 0.1, 0.2, 0.4, 0.6, 0.9, 1.0. The calculated dielectric loss spectra were in qualitative agreement with experiment. We found a slightly overestimated slowdown of the dynamics with increasing EG concentration compared to experimental data. Statistics of the hydrogen bond network and hydrogen bond lifetimes were derived from suitable time correlation functions and also show a slowdown in the dynamics with increasing xE. A similar picture is predicted for the time scales of EG conformer changes and for molecular reorientation. A slight blue shift was obtained for the power spectra of the molecular center of mass motion. The results were used to give a qualitative interpretation of the origin of three different relaxation times that appear in experimental complex dielectric spectra and of their change with xE.

Association of symmetrical alkane diols with pyridine: DFT/GIAO calculation of 1 H NMR chemical shifts

Proton nuclear magnetic resonance (NMR) shifts of the free diol and of its 1 : 1 and 1 : 2 hydrogen-bonded complexes with pyridine have been computed for five symmetrical alkane diols on the basis of density functional theory, by applying the gauge-including atomic orbital method to geometry-optimized conformers. For certain conformers, intramolecular OH···OH interactions, evidenced by high NMR OH proton shifts, are further enhanced on going from the free diol to the corresponding 1 : 1 diol/pyridine complex. This is confirmed by atoms-in-molecules and non-covalent interaction plots. The computed OH and CH proton shifts for the diol and the two complexes correlate well with values obtained by analysing data from the NMR titration of the diols in benzene against pyridine. Shift values for the diols in neat pyridine are calculated by weighting the shifts of the various protons in the three forms (free diol, 1 : 1 and 1 : 2 diol/pyridine complexes) according to the experimentally determined association constants. The results are in good agreement with those observed, and after empirical scaling, the root mean square difference is 0.18 ppm. Copyright © 2016 John Wiley & Sons, Ltd.

Ethylene glycol revisited: Molecular dynamics simulations and visualization of the liquid and its hydrogen-bond network

Molecular dynamics simulations of liquid ethylene glycol described by the OPLS-AA force field were performed to gain insight into its hydrogen-bond structure. We use the population correlation function as a statistical measure for the hydrogen-bond lifetime. In an attempt to understand the complicated hydrogen-bonding, we developed new molecular visualization tools within the Vish Visualization shell and used it to visualize the life of each individual hydrogen-bond. With this tool hydrogen-bond formation and breaking as well as clustering and chain formation in hydrogen-bonded liquids can be observed directly. Liquid ethylene glycol at room temperature does not show significant clustering or chain building. The hydrogen-bonds break often due to the rotational and vibrational motions of the molecules leading to an H-bond half-life time of approximately 1.5 ps. However, most of the H-bonds are reformed again so that after 50 ps only 40% of these H-bonds are irreversibly broken due to diffusional motion. This hydrogen-bond half-life time due to diffusional motion is 80.3 ps. The work was preceded by a careful check of various OPLS-based force fields used in the literature. It was found that they lead to quite different angular and H-bond distributions.

Polarizable empirical force field for acyclic polyalcohols based on the classical Drude oscillator

A polarizable empirical force field for acyclic polyalcohols based on the classical Drude oscillator is presented. The model is optimized with an emphasis on the transferability of the developed parameters among molecules of different sizes in this series and on the condensed-phase properties validated against experimental data. The importance of the explicit treatment of electronic polarizability in empirical force fields is demonstrated in the cases of this series of molecules with vicinal hydroxyl groups that can form cooperative intra- and intermolecular hydrogen bonds. Compared to the CHARMM additive force field, improved treatment of the electrostatic interactions avoids overestimation of the gas-phase dipole moments resulting in significant improvement in the treatment of the conformational energies and leads to the correct balance of intra- and intermolecular hydrogen bonding of glycerol as evidenced by calculated heat of vaporization being in excellent agreement with experiment. Computed condensed phase data, including crystal lattice parameters and volumes and densities of aqueous solutions are in better agreement with experimental data as compared to the corresponding additive model. Such improvements are anticipated to significantly improve the treatment of polymers in general, including biological macromolecules.

Molecular dynamics simulation of interactions between a sodium dodecyl sulfate micelle and a poly(ethylene oxide) polymer

We have performed atomistic molecular dynamics simulations of an anionic sodium dodecyl sulfate (SDS) micelle and a nonionic poly(ethylene oxide) (PEO) polymer in aqueous solution. The micelle consisted of 60 surfactant molecules, and the polymer chain lengths varied from 20 to 40 monomers. The force field parameters for PEO were adjusted by using 1,2-dimethoxymethane (DME) as a model compound and matching its hydration enthalpy and conformational behavior to experiment. Excellent agreement with previous experimental and simulation work was obtained through these modifications. The simulated scaling behavior of the PEO radius of gyration was also in close agreement with experimental results. The SDS-PEO simulations show that the polymer resides on the micelle surface and at the hydrocarbon-water interface, leading to a selective reduction in the hydrophobic contribution to the solvent-accessible surface area of the micelle. The association is mainly driven by hydrophobic interactions between the polymer and surfactant tails, while the interaction between the polymer and sulfate headgroups on the micelle surface is weak. The 40-monomer chain is mostly wrapped around the micelle, and nearly 90% of the monomers are adsorbed at low PEO concentration. Simulations were also performed with multiple 20-monomer chains, and gradual addition of polymer indicates that about 120 monomers are required to saturate the micelle surface. The stoichiometry of the resulting complex is in close agreement with experimental results, and the commonly accepted "beaded necklace" structure of the SDS-PEO complex is recovered by our simulations.