Premium
A polarizable empirical force field for molecular dynamics simulation of liquid hydrocarbons
Author(s) -
Szklarczyk Oliwia M.,
Bachmann Stephan J.,
van Gunsteren Wilfred F.
Publication year - 2014
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/jcc.23551
Subject(s) - molecular dynamics , force field (fiction) , polarizability , field (mathematics) , statistical physics , chemical physics , dynamics (music) , computational chemistry , physics , chemistry , molecule , quantum mechanics , mathematics , acoustics , pure mathematics
Electronic polarizability is usually treated implicitly in molecular simulations, which may lead to imprecise or even erroneous molecular behavior in spatially electronically inhomogeneous regions of systems such as proteins, membranes, interfaces between compounds, or mixtures of solvents. The majority of available molecular force fields and molecular dynamics simulation software packages does not account explicitly for electronic polarization. Even the simplest charge‐on‐spring (COS) models have only been developed for few types of molecules. In this work, we report a polarizable COS model for cyclohexane, as this molecule is a widely used solvent, and for linear alkanes, which are also used as solvents, and are the precursors of lipids, amino acid side chains, carbohydrates, or nucleic acid backbones. The model is an extension of a nonpolarizable united‐atom model for alkanes that had been calibrated against experimental values of the density, the heat of vaporization and the Gibbs free energy of hydration for each alkane. The latter quantity was used to calibrate the parameters governing the interaction of the polarizable alkanes with water. Subsequently, the model was tested for other structural, thermodynamic, dielectric, and dynamic properties such as trans / gauche ratios, excess free energy, static dielectric permittivity, and self‐diffusion. A good agreement with the experimental data for a large set of properties for each considered system was obtained, resulting in a transferable set of polarizable force‐field parameters for CH 2 , CH 3 , and CH 4 moieties. © 2014 Wiley Periodicals, Inc.