Effects of Temperature on the Properties of Glycerol: A Computer Simulation Study of Five Different Force Fields

We perform molecular dynamics simulations of glycerol (propane-1,2,3-triol) at normal pressure and a wide range of temperatures (300–460 K) and study the sensitivity of simulation results to the force field (FF) considered. We employ five commonly used FFs: (i) AMBER, (ii) CHARMM22, and (iii) three...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:The journal of physical chemistry. B 2014-09, Vol.118 (38), p.11284-11294
Hauptverfasser: Jahn, David A, Akinkunmi, Frederick O, Giovambattista, Nicolas
Format: Artikel
Sprache:eng
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:We perform molecular dynamics simulations of glycerol (propane-1,2,3-triol) at normal pressure and a wide range of temperatures (300–460 K) and study the sensitivity of simulation results to the force field (FF) considered. We employ five commonly used FFs: (i) AMBER, (ii) CHARMM22, and (iii) three versions of the OPLS-AA FF (OPLS1, OPLS2, and OPLS3). We study thermodynamic (density ρ(T), thermal expansion coefficient α P (T), isobaric specific heat c P (T)), dynamic (diffusion coefficient D(T)), as well as structural properties (molecular conformations and hydrogen-bond statistics). In comparison with experiments, FFs i and iii provide reasonable estimations of ρ(T) with deviations of ≤4.5%; for FF ii, deviations in density are more pronounced, ≤9%. Values of α P (T) vary considerably among the FFs; e.g., deviations are ≤9% for OPLS1-FF and ≤60% for FF ii. For all models studied, values of c P (T) are approximately twice the corresponding experimental values. Diffusion coefficients are very sensitive to the FFs considered. Specifically, for FFs i and ii and OPLS3, the values of D(T) are remarkably close to the experimental values over the whole range of temperatures studied. Instead, in the cases of OPLS1 and OPLS2-FFs, D(T) is underestimated by approximately 2 orders of magnitude. Interestingly, in all cases, D(T) can be well described by a Vogel–Tamman–Fulcher equation, as observed in experiments. We present a detailed characterization of glycerol backbone conformation based on the traditional classification introduced by Bastiansen, defined in terms of glycerol’s OCCC dihedral angles. All FFs indicate that the conformer population varies smoothly with temperature. However, the FFs provide very different conformer distributions. This implies that, from the microscopic point of view, these glycerol models may provide very different liquid environments for, for example, guest biomolecules and hence may play a relevant role in interpreting simulation results involving glycerol-based solutions. We also discuss the statistics of inter- and intramolecular hydrogen bonds (HBs). The FFs are qualitatively comparable regarding HB statistics; however, quantitative differences remain. For example, molecules form a total of 5.5–7 HBs at T = 350 K, depending on the FF considered, including at least one intramolecular HB.
ISSN:1520-6106
1520-5207
DOI:10.1021/jp5059098