Why many semiempirical molecular orbital theories fail for liquid water and how to fix them

Water is an extremely important liquid for chemistry and the search for more accurate force fields for liquid water continues unabated. Neglect of diatomic differential overlap (NDDO) molecular orbital methods provide and intriguing generalization of classical force fields in this regard because they can account both for bond breaking and electronic polarization of molecules. However, we show that most standard NDDO methods fail for water because they give an incorrect description of hydrogen bonding, water's key structural feature. Using force matching, we design a reparameterized NDDO model and find that it qualitatively reproduces the experimental radial distribution function of water, as well as various monomer, dimer, and bulk properties that PM6 does not. This suggests that the apparent limitations of NDDO models are primarily due to poor parameterization and not to the NDDO approximations themselves. Finally, we identify the physical parameters that most influence the condensed phase properties. These results help to elucidate the chemistry that a semiempirical molecular orbital picture of water must capture. We conclude that properly parameterized NDDO models could be useful for simulations that require electronically detailed explicit solvent, including the calculation of redox potentials and simulation of charge transfer and photochemistry. © 2015 Wiley Periodicals, Inc.

[1]  Antoine Marion,et al.  Water interactions with hydrophobic groups: assessment and recalibration of semiempirical molecular orbital methods. , 2014, The Journal of chemical physics.

[2]  G. A. Jeffrey,et al.  An Introduction to Hydrogen Bonding , 1997 .

[3]  A. D. Buckingham,et al.  A theory of the dielectric polarization of polar substances , 1956, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[4]  Feng Wang,et al.  Application of a Drude model to the binding of excess electrons to water clusters , 2002 .

[5]  F. Weigend,et al.  Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculations , 2002 .

[6]  Haibo Yu,et al.  Accounting for polarization in molecular simulation , 2005, Comput. Phys. Commun..

[7]  Arieh Warshel,et al.  Polarizable Force Fields:  History, Test Cases, and Prospects. , 2007, Journal of chemical theory and computation.

[8]  Steven J. Stuart,et al.  Dynamical fluctuating charge force fields: Application to liquid water , 1994 .

[9]  J. Stewart Optimization of parameters for semiempirical methods I. Method , 1989 .

[10]  Gregory A Voth,et al.  Effective force fields for condensed phase systems from ab initio molecular dynamics simulation: a new method for force-matching. , 2004, The Journal of chemical physics.

[11]  Jeff R. Hammond,et al.  Accurate dipole polarizabilities for water clusters n=2-12 at the coupled-cluster level of theory and benchmarking of various density functionals. , 2009, The Journal of chemical physics.

[12]  Michiel Sprik,et al.  Ab initio molecular dynamics simulation of liquid water: Comparison of three gradient‐corrected density functionals , 1996 .

[13]  J. Pople,et al.  Approximate Self-Consistent Molecular Orbital Theory. I. Invariant Procedures , 1965 .

[14]  S. Papson “Model” , 1981 .

[15]  Peter A. Kollman,et al.  Hydrogen bonded dimers and polymers involving hydrogen fluoride, water, and ammonia , 1970 .

[16]  S. Nosé A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[17]  Benoît Roux,et al.  Modeling induced polarization with classical Drude oscillators: Theory and molecular dynamics simulation algorithm , 2003 .

[18]  D. Truhlar,et al.  Direct dynamics calculations with NDDO (neglect of diatomic differential overlap) molecular orbital theory with specific reaction parameters , 1991 .

[19]  David E. Woon,et al.  Gaussian basis sets for use in correlated molecular calculations. IV. Calculation of static electrical response properties , 1994 .

[20]  Kari Laasonen,et al.  ‘‘Ab initio’’ liquid water , 1993 .

[21]  John S. Muenter,et al.  The dipole moment of water. I. Dipole moments and hyperfine properties of H2O and HDO in the ground and excited vibrational states , 1991 .

[22]  M. Holz,et al.  Biological applications of scanning tunnelling microscopy , 1993 .

[23]  Martin Korth,et al.  Third-Generation Hydrogen-Bonding Corrections for Semiempirical QM Methods and Force Fields , 2010 .

[24]  M. Tuckerman,et al.  Dynamical properties of liquid water from ab initio molecular dynamics performed in the complete basis set limit. , 2007, The Journal of chemical physics.

[25]  Yan Zhou,et al.  Reaction Path Force Matching: A New Strategy of Fitting Specific Reaction Parameters for Semiempirical Methods in Combined QM/MM Simulations. , 2014, Journal of chemical theory and computation.

[26]  Larry A. Curtiss,et al.  Studies of molecular association in H2O and D2O vapors by measurement of thermal conductivity , 1979 .

[27]  Darrin M York,et al.  Quantum mechanical study of solvent effects in a prototype SN2 reaction in solution: Cl- attack on CH3Cl. , 2014, The Journal of chemical physics.

[28]  John S. Muenter,et al.  THE STRUCTURE OF WATER DIMER FROM MOLECULAR BEAM ELECTRIC RESONANCE SPECTROSCOPY: PARTIALLY DEUTERATED DIMERS , 1977 .

[29]  Joost VandeVondele,et al.  The influence of temperature and density functional models in ab initio molecular dynamics simulation of liquid water. , 2005, The Journal of chemical physics.

[30]  Shawn T. Brown,et al.  Advances in methods and algorithms in a modern quantum chemistry program package. , 2006, Physical chemistry chemical physics : PCCP.

[31]  T. H. Dunning Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

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

[33]  Michael C. Zerner,et al.  Semiempirical Molecular Orbital Methods , 2007 .

[34]  Philip Ball,et al.  Water: Water — an enduring mystery , 2008, Nature.

[35]  Karl N. Kirschner,et al.  Modeling of magic water clusters (H2O)20 and (H2O)21H+ with the PM3 quantum‐mechanical method , 1993, J. Comput. Chem..

[36]  Alan K. Soper,et al.  The radial distribution functions of water and ice from 220 to 673 K and at pressures up to 400 MPa , 2000 .

[37]  J. Stewart Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements , 2007, Journal of molecular modeling.

[38]  C. Vega,et al.  A general purpose model for the condensed phases of water: TIP4P/2005. , 2005, The Journal of chemical physics.

[39]  M. Holz,et al.  Temperature-dependent self-diffusion coefficients of water and six selected molecular liquids for calibration in accurate 1H NMR PFG measurements , 2000 .

[40]  Shoshana J. Wodak,et al.  Molecular dynamics simulation of polarizable water by an extended Lagrangian method , 1992 .

[41]  Bertrand Guillot,et al.  A reappraisal of what we have learnt during three decades of computer simulations on water , 2002 .

[42]  A. R. H. Goodwin,et al.  A Database for the Static Dielectric Constant of Water and Steam , 1995 .

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

[44]  James B. Adams,et al.  Interatomic Potentials from First-Principles Calculations: The Force-Matching Method , 1993, cond-mat/9306054.

[45]  Peter Pulay,et al.  Second-order Møller–Plesset calculations with dual basis sets , 2003 .

[46]  C. Cramer,et al.  PM3‐SM3: A general parameterization for including aqueous solvation effects in the PM3 molecular orbital model , 1992 .

[47]  William F. Murphy,et al.  The Rayleigh depolarization ratio and rotational Raman spectrum of water vapor and the polarizability components for the water molecule , 1977 .

[48]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[49]  F. Weigend,et al.  RI-MP2: first derivatives and global consistency , 1997 .

[50]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[51]  James J. P. Stewart,et al.  Optimization of parameters for semiempirical methods VI: more modifications to the NDDO approximations and re-optimization of parameters , 2012, Journal of Molecular Modeling.

[52]  Julian Tirado-Rives,et al.  Potential energy functions for atomic-level simulations of water and organic and biomolecular systems. , 2005, Proceedings of the National Academy of Sciences of the United States of America.