Abstract

Machine learning algorithms have been demonstrated to predict atomistic properties approaching the accuracy of quantum chemical calculations at significantly less computational cost. Difficulties arise, however, when attempting to apply these techniques to large systems, or systems possessing excessive conformational freedom. In this article, the machine learning method kriging is applied to predict both the intra‐atomic and interatomic energies, as well as the electrostatic multipole moments, of the atoms of a water molecule at the center of a 10 water molecule (decamer) cluster. Unlike previous work, where the properties of small water clusters were predicted using a molecular local frame, and where training set inputs (features) were based on atomic index, a variety of feature definitions and coordinate frames are considered here to increase prediction accuracy. It is shown that, for a water molecule at the center of a decamer, no single method of defining features or coordinate schemes is optimal for every property. However, explicitly accounting for the structure of the first solvation shell in the definition of the features of the kriging training set, and centring the coordinate frame on the atom‐of‐interest will, in general, return better predictions than models that apply the standard methods of feature definition, or a molecular coordinate frame. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Highlights

  • Molecular dynamics (MD) simulations are an important tool in understanding the dynamical evolution of condensed matter systems

  • By exploiting the physical indistinguishability of the water molecules within the cluster to redefine the features of the training set, we show that training density can be increased, at no additional cost

  • All training sets contained the same conformational information, careful definition of training set features to account for local structure led to improvements in mean absolute prediction error of up to $75% for certain properties

Read more

Summary

Introduction

Molecular dynamics (MD) simulations are an important tool in understanding the dynamical evolution of condensed matter systems. In contrast to computationally expensive ab initio MD techniques,[1,2] most condensed phase simulations currently rely on one of a number of parameterizable force fields, including CHARMM[3] and AMBER,[4] among others Such force fields commonly treat electrostatics as pairwise interactions between point charges, and may describe the energy variation of a molecular bond under compression or elongation through a simple Hooke potential. The necessity of including the underlying quantum mechanical effects in the description of simulated water molecules has been acknowledged by several authors.[10,11]

Methods
Results
Conclusion

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.