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

Citations Over TimeTop 10% of 2013 papers

Abstract

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.

Related Papers

• → Secondary structure dependence of amyloid‐β(1–40) on simulation techniques and force field parameters(2021)24 cited
• → Optimized all-atom force field for alkynes within the OPLS-AA framework(2021)15 cited
• → Optimisation of OPLS–UA force-field parameters for protein systems using protein data bank(2010)14 cited
• → Condensed phase properties of n-pentadecane emerging from application of biomolecular force fields(2020)4 cited
• → B211 Molecular dynamics simulation of Ethylene Glycol and Water mixtures(2014)