z-logo
open-access-imgOpen Access
How Short Is Too Short for the Interactions of a Water Potential? Exploring the Parameter Space of a Coarse-Grained Water Model Using Uncertainty Quantification
Author(s) -
Liam C. Jacobson,
Robert M. Kirby,
Valeria Molinero
Publication year - 2014
Publication title -
the journal of physical chemistry b
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.864
H-Index - 392
eISSN - 1520-6106
pISSN - 1520-5207
DOI - 10.1021/jp5012928
Subject(s) - parameter space , degrees of freedom (physics and chemistry) , monatomic ion , statistical physics , water model , sensitivity (control systems) , thermal diffusivity , computer science , biological system , mathematics , physics , molecular dynamics , thermodynamics , statistics , chemistry , computational chemistry , quantum mechanics , electronic engineering , biology , engineering
Coarse-grained models are becoming increasingly popular due to their ability to access time and length scales that are prohibitively expensive with atomistic models. However, as a result of decreasing the degrees of freedom, coarse-grained models often have diminished accuracy, representability, and transferability compared with their finer grained counterparts. Uncertainty quantification (UQ) can help alleviate this challenge by providing an efficient and accurate method to evaluate the effect of model parameters on the properties of the system. This method is useful in finding parameter sets that fit the model to several experimental properties simultaneously. In this work we use UQ as a tool for the evaluation and optimization of a coarse-grained model. We efficiently sample the five-dimensional parameter space of the coarse-grained monatomic water (mW) model to determine what parameter sets best reproduce experimental thermodynamic, structural and dynamical properties of water. Generalized polynomial chaos (gPC) was used to reconstruct the analytical surfaces of density, enthalpy of vaporization, radial and angular distribution functions, and diffusivity of liquid water as a function of the input parameters. With these surfaces, we evaluated the sensitivity of these properties to perturbations of the model input parameters and the accuracy and representability of the coarse-grained models. In particular, we investigated what is the optimum length scale of the water-water interactions needed to reproduce the properties of liquid water with a monatomic model with two- and three-body interactions. We found that there is an optimum cutoff length of 4.3 Å, barely longer than the size of the first neighbor shell in water. As cutoffs deviate from this optimum value, the ability of the model to simultaneously reproduce the structure and thermodynamics is severely diminished.

The content you want is available to Zendy users.

Already have an account? Click here to sign in.
Having issues? You can contact us here
Accelerating Research

Address

John Eccles House
Robert Robinson Avenue,
Oxford Science Park, Oxford
OX4 4GP, United Kingdom