Ab initio molecular-dynamics simulation of aqueous proton solvation and transport revisited.

The solvation and transport of the hydrated excess proton is studied using the Car-Parrinello molecular-dynamics (CPMD) simulation method. The simulations were performed using BLYP and HCTH gradient-corrected exchange-correlation energy functionals. The fictitious electronic mass was chosen to be small enough so that the underlying water structural and dynamical properties were converged with respect to this important CPMD simulation parameter. An unphysical overstructuring of liquid water in the CPMD simulations using the BLYP functional resulted in the formation of long-lived hydrogen-bonding structures involving the excess proton and a particular (special) water oxygen. The excess proton was observed to be attracted to the special oxygen through the entire length of the BLYP CPMD simulations. Consequently, the excess proton diffusion was limited by the mobility of the special oxygen in the slowly diffusing water network and, in turn, the excess proton self-diffusion coefficient was found to be significantly below the experimental value. On the other hand, the structural properties of liquid water in the HCTH CPMD simulation were seen to be in better agreement with experiment, although the water and excess proton diffusions were still well below the experimental value.

[1]  G. Voth,et al.  Car–Parrinello molecular dynamics simulation of liquid water: New results , 2002 .

[2]  Fred A. Hamprecht,et al.  Development and assessment of new exchange-correlation functionals , 1998 .

[3]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[4]  Gregory A. Voth,et al.  The computer simulation of proton transport in water , 1999 .

[5]  Aaron Lefohn,et al.  A Multistate Empirical Valence Bond Approach to a Polarizable and Flexible Water Model , 2001 .

[6]  L. Dang.  Solvation of the hydronium ion at the water liquid/vapor interface , 2003 .

[7]  Carlo Adamo,et al.  Predicting proton transfer barriers with density functional methods , 1999 .

[8]  Hadas Lapid,et al.  A bond-order analysis of the mechanism for hydrated proton mobility in liquid water. , 2005, The Journal of chemical physics.

[9]  Bin Chen,et al.  Liquid Water from First Principles: Investigation of Different Sampling Approaches , 2004 .

[10]  Michele Parrinello,et al.  On the Quantum Nature of the Shared Proton in Hydrogen Bonds , 1997, Science.

[11]  Martins,et al.  Efficient pseudopotentials for plane-wave calculations. , 1991, Physical review. B, Condensed matter.

[12]  Kari Laasonen,et al.  Ab initio molecular dynamics simulation of the solvation and transport of hydronium and hydroxyl ions in water , 1995 .

[13]  M. Klein,et al.  Ab initio molecular dynamics study of proton transfer in a polyglycine analog of the ion channel gramicidin A. , 1996, Biophysical journal.

[14]  Joseph A Morrone,et al.  Ab initio molecular dynamics study of proton mobility in liquid methanol , 2002 .

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

[16]  Free energy of liquid water on the basis of quasichemical theory and ab initio molecular dynamics. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  Nicholas C. Handy,et al.  Development of New Exchange-Correlation Functionals. 2 , 1998 .

[18]  J. Kress,et al.  Ab initio molecular dynamics and quasichemical study of H+(aq). , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[19]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[20]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[21]  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 .

[22]  M. Holz,et al.  Isotope effect on the translational and rotational motion in liquid water and ammonia , 2001 .

[23]  Michiel Sprik,et al.  New generalized gradient approximation functionals , 2000 .

[24]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[25]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[26]  A. Soper,et al.  Ions in water: the microscopic structure of a concentrated HCl solution. , 2004, The Journal of chemical physics.

[27]  Daniel Borgis,et al.  Transport and spectroscopy of the hydrated proton: A molecular dynamics study , 1999 .

[28]  Sandro Scandolo,et al.  How well do Car-Parrinello simulations reproduce the Born-Oppenheimer surface ? , 2001, cond-mat/0107415.

[29]  J. Grotendorst,et al.  Modern methods and algorithms of quantum chemistry : winterschool 21. - 25. February 2000 Forschungszentrum Jülich : proceedings / org. by John von Neumann Institute for Computing , 2000 .

[30]  Michele Parrinello,et al.  Structural, electronic, and bonding properties of liquid water from first principles , 1999 .

[31]  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.

[32]  N. Agmon,et al.  The Grotthuss mechanism , 1995 .

[33]  Johannes Grotendorst,et al.  Modern methods and algorithms of quantum chemistry , 2000 .

[34]  M. Parrinello,et al.  The nature of the hydrated excess proton in water , 1999, Nature.

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

[36]  Jens Ulstrup,et al.  Kinetics of Proton Transport in Water , 2003 .

[37]  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.