Convergence of the Many-Body Expansion for Energy and Forces for Classical Polarizable Models in the Condensed Phase.

We analyze convergence of energies and forces for the AMOEBA classical polarizable model when evaluated as a many-body expansion (MBE) against the corresponding N-body parent potential in the context of a condensed-phase water simulation. This is in contrast to most MBE formulations based on quantum mechanics, which focus only on convergence of energies for gas-phase clusters. Using a single water molecule as a definition of a body, we find that truncation of the MBE at third order, 3-AMOEBA, captures direct polarization exactly and yields apparent good convergence of the mutual polarization energy. However, it renders large errors in the magnitude of polarization forces and requires at least fourth-order terms in the MBE to converge toward the parent potential gradient values. We can improve the convergence of polarization forces for 3-AMOEBA by embedding the polarization response of dimers and trimers within a complete representation of the fixed electrostatics of the entire system. We show that the electrostatic embedding formalism helps identify the specific configurations involving linear hydrogen-bonding arrangements that are poorly convergent at the 3-body level. By extending the definition of a body to be a large water cluster, we can reduce errors in forces to yield an approximate polarization model that is up to 10 times faster than the parent potential. The 3-AMOEBA model offers new ways to investigate how the properties of bulk water depend on the degree of connectivity in the liquid.

[1]  Jules W. Moskowitz,et al.  Water Molecule Interactions , 1970 .

[2]  Ryan M. Richard,et al.  Understanding the many-body expansion for large systems. II. Accuracy considerations. , 2016, The Journal of chemical physics.

[3]  P. J. Bygrave,et al.  The embedded many-body expansion for energetics of molecular crystals. , 2012, The Journal of chemical physics.

[4]  Louis A. Hageman,et al.  Iterative Solution of Large Linear Systems. , 1971 .

[5]  Andreas W Götz,et al.  On the representation of many-body interactions in water. , 2015, The Journal of chemical physics.

[6]  Donald G Truhlar,et al.  Electrostatically Embedded Many-Body Expansion for Large Systems, with Applications to Water Clusters. , 2007, Journal of chemical theory and computation.

[7]  Pengyu Y. Ren,et al.  Systematic improvement of a classical molecular model of water. , 2013, The journal of physical chemistry. B.

[8]  H. Eugene Stanley,et al.  Low-Density "Patches" in the Hydrogen-Bond Network of Liquid Water: Evidence from Molecular-Dynamics Computer Simulations , 1982 .

[9]  N. Govind,et al.  Electronic-structure calculations by first-principles density-based embedding of explicitly correlated systems , 1999 .

[10]  Frederick R. Manby,et al.  A Simple, Exact Density-Functional-Theory Embedding Scheme , 2012, Journal of chemical theory and computation.

[11]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[12]  John M Herbert,et al.  A generalized many-body expansion and a unified view of fragment-based methods in electronic structure theory. , 2012, The Journal of chemical physics.

[13]  Kazuo Kitaura,et al.  The Fragment Molecular Orbital Method: Practical Applications to Large Molecular Systems , 2009 .

[14]  John M Herbert,et al.  Aiming for benchmark accuracy with the many-body expansion. , 2014, Accounts of chemical research.

[15]  Wei Wang,et al.  Fast evaluation of polarizable forces. , 2005, The Journal of chemical physics.

[16]  Kazuo Kitaura,et al.  Extending the power of quantum chemistry to large systems with the fragment molecular orbital method. , 2007, The journal of physical chemistry. A.

[17]  Qing Peng,et al.  Time-dependent potential-functional embedding theory. , 2014, The Journal of chemical physics.

[18]  Alán Aspuru-Guzik,et al.  Prediction and Calculation of Crystal Structures: Methods and Applications , 2014 .

[19]  Frederick R Manby,et al.  Exact nonadditive kinetic potentials for embedded density functional theory. , 2010, The Journal of chemical physics.

[20]  W. C. Swope,et al.  A computer simulation method for the calculation of equilibrium constants for the formation of physi , 1981 .

[21]  Alfons Geiger,et al.  Aspects of the percolation process for hydrogen‐bond networks in water , 1979 .

[22]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole-based Molecular Mechanics for Organic Molecules. , 2011, Journal of chemical theory and computation.

[23]  Chen Huang,et al.  Quantum mechanical embedding theory based on a unique embedding potential. , 2011, The Journal of chemical physics.

[24]  John Z H Zhang,et al.  New Advance in Computational Chemistry:  Full Quantum Mechanical ab Initio Computation of Streptavidin-Biotin Interaction Energy. , 2003, The journal of physical chemistry. B.

[25]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

[26]  Frederick R Manby,et al.  Density functional theory embedding for correlated wavefunctions: improved methods for open-shell systems and transition metal complexes. , 2012, The Journal of chemical physics.

[27]  Donald G Truhlar,et al.  Electrostatically embedded many-body method for dipole moments, partial atomic charges, and charge transfer. , 2012, Physical chemistry chemical physics : PCCP.

[28]  Jiali Gao,et al.  A molecular-orbital derived polarization potential for liquid water , 1998 .

[29]  K. Kitaura,et al.  Fragment molecular orbital method: an approximate computational method for large molecules , 1999 .

[30]  Clifford E. Dykstra,et al.  Electrostatic interaction potentials in molecular force fields , 1993 .

[31]  Teresa Head-Gordon,et al.  Advanced potential energy surfaces for condensed phase simulation. , 2014, Annual review of physical chemistry.

[32]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation , 2003 .

[33]  G. Beran,et al.  Accurate Molecular Crystal Lattice Energies from a Fragment QM/MM Approach with On-the-Fly Ab Initio Force Field Parametrization. , 2011, Journal of chemical theory and computation.

[34]  Spencer R Pruitt,et al.  Fragmentation methods: a route to accurate calculations on large systems. , 2012, Chemical reviews.

[35]  John M Herbert,et al.  Understanding the many-body expansion for large systems. I. Precision considerations. , 2014, The Journal of chemical physics.

[36]  Chen Huang,et al.  Potential-functional embedding theory for molecules and materials. , 2011, The Journal of chemical physics.

[37]  H. Ågren,et al.  The principles of infrared-x-ray pump-probe spectroscopy. Applications on proton transfer in core-ionized water dimers. , 2005, The Journal of chemical physics.

[38]  Jiali Gao,et al.  Toward a Molecular Orbital Derived Empirical Potential for Liquid Simulations , 1997 .

[39]  P. J. Bygrave,et al.  Energy benchmarks for water clusters and ice structures from an embedded many-body expansion. , 2013, The Journal of chemical physics.

[40]  Mark S. Gordon,et al.  The Effective Fragment Potential Method: A QM-Based MM Approach to Modeling Environmental Effects in Chemistry , 2001 .

[41]  Mark S. Gordon,et al.  An effective fragment method for modeling solvent effects in quantum mechanical calculations , 1996 .

[42]  Emily A Carter,et al.  Self-consistent embedding theory for locally correlated configuration interaction wave functions in condensed matter. , 2006, The Journal of chemical physics.

[43]  Teresa Head-Gordon,et al.  An efficient and stable hybrid extended Lagrangian/self-consistent field scheme for solving classical mutual induction. , 2015, The Journal of chemical physics.