Incorporating Phase-Dependent Polarizability in Non-Additive Electrostatic Models for Molecular Dynamics Simulations of the Aqueous Liquid-Vapor Interface.

We discuss a new classical water force field that explicitly accounts for differences in polarizability between liquid and vapor phases. The TIP4P-QDP (4-point transferable intermolecular potential with charge dependent-polarizability) force field is a modification of the original TIP4P-FQ fluctuating charge water force field of Rick et al.(1) that self-consistently adjusts its atomic hardness parameters via a scaling function dependent on the M-site charge. The electronegativity (χ) parameters are also scaled in order to reproduce condensed-phase dipole moments of comparable magnitude to TIP4P-FQ. TIP4P-QDP is parameterized to reproduce experimental gas-phase and select condensed-phase properties. The TIP4P-QDP water model possesses a gas phase polarizability of 1.40 Å(3) and gas-phase dipole moment of 1.85 Debye, in excellent agreement with experiment and high-level ab initio predictions. The liquid density of TIP4P-QDP is 0.9954(±0.0002) g/cm(3) at 298 K and 1 atmosphere, and the enthalpy of vaporization is 10.55(±0.12) kcal/mol. Other condensed-phase properties such as the isobaric heat capacity, isothermal compressibility, and diffusion constant are also calculated within reasonable accuracy of experiment and consistent with predictions of other current state-of-the-art water force fields. The average molecular dipole moment of TIP4P-QDP in the condensed phase is 2.641(±0.001) Debye, approximately 0.02 Debye higher than TIP4P-FQ and within the range of values currently surmised for the bulk liquid. The dielectric constant, ε = 85.8 ± 1.0, is 10% higher than experiment. This is reasoned to be due to the increase in the condensed phase dipole moment over TIP4P-FQ, which estimates ε remarkably well. Radial distribution functions for TIP4P-QDP and TIP4P-FQ show similar features, with TIP4P-QDP showing slightly reduced peak heights and subtle shifts towards larger distance interactions. Since the greatest effects of the phase-dependent polarizability are anticipated in regions with both liquid and vapor character, interfacial simulations of TIP4P-QDP were performed and compared to TIP4P-FQ, a static polarizability analog. Despite similar features in density profiles such as the position of the GDS and interfacial width, enhanced dipole moments are observed for the TIP4P-QDP interface and onset of the vapor phase. Water orientational profiles show an increased preference (over TIP4P-FQ) in the orientation of the permanent dipole vector of the molecule within the interface; an enhanced z-induced dipole moment directly results from this preference. Hydrogen bond formation is lower, on average, in the bulk for TIP4P-QDP than TIP4P-FQ. However, the average number of hydrogen bonds formed by TIP4P-QDP in the interface exceeds that of TIP4P-FQ, and observed hydrogen bond networks extend further into the gaseous region. The TIP4P-QDP interfacial potential, calculated to be -11.98(±0.08) kcal/mol, is less favorable than that for TIP4P-FQ by approximately 2% as a result of a diminished quadrupole contribution. Surface tension is calculated within a 1.3% reduction from the experimental value. Results reported demonstrate TIP4P-QDP as a model comparable to the popular TIP4P-FQ while accounting for a physical effect previously neglected by other water models. Further refinements to this model, as well as future applications are discussed.

[1]  Jacek Korchowiec,et al.  Molecular hardness and softness parameters and their use in chemistry , 1988 .

[2]  Harry A. Stern,et al.  Development of a polarizable force field for proteins via ab initio quantum chemistry: First generation model and gas phase tests , 2002, J. Comput. Chem..

[3]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[4]  L. E. Chirlian,et al.  Atomic charges derived from electrostatic potentials: A detailed study , 1987 .

[5]  T. Halgren,et al.  Polarizable force fields. , 2001, Current opinion in structural biology.

[6]  Charles L. Brooks,et al.  Fluctuating charge force fields: recent developments and applications from small molecules to macromolecular biological systems , 2006 .

[7]  S. Nosé.  A molecular dynamics method for simulations in the canonical ensemble , 1984 .

[8]  G. Richmond,et al.  Understanding the Effects of Hydrogen Bonding at the Vapor−Water Interface: Vibrational Sum Frequency Spectroscopy of H2O/HOD/D2O Mixtures Studied Using Molecular Dynamics Simulations , 2007 .

[9]  L. Dang,et al.  Molecular Mechanism of Ion Binding to the Liquid/Vapor Interface of Water , 2002 .

[10]  L. Dang,et al.  Many-body interactions in liquid methanol and its liquid/vapor interface: A molecular dynamics study , 2003 .

[11]  M. Berkowitz,et al.  Structure and dynamics of Cl−(H2O)20 clusters: The effect of the polarizability and the charge of the ion , 1992 .

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

[13]  L. Renee Olano,et al.  Fluctuating charge normal modes: An algorithm for implementing molecular dynamics simulations with polarizable potentials , 2005, J. Comput. Chem..

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

[15]  Alexander D. MacKerell.  Empirical force fields for biological macromolecules: Overview and issues , 2004, J. Comput. Chem..

[16]  R. T. Sanderson,et al.  An Interpretation of Bond Lengths and a Classification of Bonds. , 1951, Science.

[17]  Wilfried J. Mortier,et al.  Electronegativity-equalization method for the calculation of atomic charges in molecules , 1986 .

[18]  Brad A. Bauer,et al.  Condensed-phase properties of n-alkyl-amines from molecular dynamics simulations using charge equilibration force fields , 2008 .

[19]  J. Tse,et al.  Transformation of ice VIII to amorphous ice by ‘‘melting’’ at low temperature , 1989 .

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

[21]  Ilan Benjamin,et al.  Chemical Reactions and Solvation at Liquid Interfaces: A Microscopic Perspective. , 1996, Chemical reviews.

[22]  Harry A. Stern,et al.  Fluctuating Charge, Polarizable Dipole, and Combined Models: Parameterization from ab Initio Quantum Chemistry , 1999 .

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

[24]  Michael L. Klein,et al.  Effective pair potentials and the properties of water , 1989 .

[25]  G. Richmond,et al.  Depth profiling of water molecules at the liquid-liquid interface using a combined surface vibrational spectroscopy and molecular dynamics approach. , 2007, Journal of the American Chemical Society.

[26]  Alan K. Soper,et al.  A new determination of the structure of water at 25°C , 1986 .

[27]  K. Eisenthal,et al.  Protein Adsorption at Interfaces Detected by Second Harmonic Generation , 2000 .

[28]  Ronald M. Levy,et al.  Computer simulations of the dielectric properties of water: Studies of the simple point charge and transferrable intermolecular potential models , 1989 .

[29]  L. Dang,et al.  A Mechanism for Ion Transport across the Water/Dichloromethane Interface: A Molecular Dynamics Study Using Polarizable Potential Models , 2001 .

[30]  B. Berne,et al.  Combined fluctuating charge and polarizable dipole models: Application to a five-site water potential function , 2001 .

[31]  Sandeep Patel,et al.  Origin and control of superlinear polarizability scaling in chemical potential equalization methods. , 2008, The Journal of chemical physics.

[32]  Greg L. Hura,et al.  Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. , 2004, The Journal of chemical physics.

[33]  Alexander D. MacKerell,et al.  Polarizable empirical force field for alkanes based on the classical Drude oscillator model. , 2005, The journal of physical chemistry. B.

[34]  Gilles H. Peslherbe,et al.  Importance of polarization in quantum mechanics/molecular mechanics descriptions of electronic excited states: NaI(H2O)n photodissociation dynamics as a case study. , 2008, The journal of physical chemistry. B.

[35]  I-Feng W. Kuo,et al.  An ab Initio Molecular Dynamics Study of the Aqueous Liquid-Vapor Interface , 2004, Science.

[36]  D. J. Price,et al.  A modified TIP3P water potential for simulation with Ewald summation. , 2004, The Journal of chemical physics.

[37]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[38]  A. Pohorille,et al.  Comment on "Study on the liquid-vapor interface of water. I. Simulation results of thermodynamic properties and orientational structure". , 1989, The Journal of chemical physics.

[39]  George A. Kaminski,et al.  Development of an Accurate and Robust Polarizable Molecular Mechanics Force Field from ab Initio Quantum Chemistry , 2004 .

[40]  Nohad Gresh,et al.  Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy. , 2007, Journal of chemical theory and computation.

[41]  Karl K. Irikura,et al.  Predicting Unexpected Chemical Reactions by Isopotential Searching , 2000 .

[42]  Charles L Brooks,et al.  Revisiting the hexane-water interface via molecular dynamics simulations using nonadditive alkane-water potentials. , 2006, The Journal of chemical physics.

[43]  Mauro C. C. Ribeiro,et al.  Fluctuating charge model for polyatomic ionic systems: A test case with diatomic anions , 1999 .

[44]  J. Kirkwood,et al.  The Statistical Mechanical Theory of Surface Tension , 1949 .

[45]  David E. Smith,et al.  Computer simulations of NaCl association in polarizable water , 1994 .

[46]  Sandeep Patel,et al.  Electrostatic properties of aqueous salt solution interfaces: a comparison of polarizable and nonpolarizable ion models. , 2008, The journal of physical chemistry. B.

[47]  Nohad Gresh,et al.  Toward a Separate Reproduction of the Contributions to the Hartree-Fock and DFT Intermolecular Interaction Energies by Polarizable Molecular Mechanics with the SIBFA Potential. , 2007, Journal of chemical theory and computation.

[48]  Ursula Rothlisberger,et al.  Polarization effects and charge transfer in the KcsA potassium channel. , 2006, Biophysical chemistry.

[49]  Alexander D. MacKerell,et al.  CHARMM fluctuating charge force field for proteins: II Protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model , 2004, J. Comput. Chem..

[50]  M. McGrath,et al.  Spatial correlation of dipole fluctuations in liquid water , 2007 .

[51]  R. Saykally,et al.  Enhanced Concentration of Polarizable Anions at the Liquid Water Surface: SHG Spectroscopy and MD Simulations of Sodium Thiocyanide , 2005 .

[52]  A. Morita,et al.  Molecular Dynamics Study of Gas−Liquid Aqueous Sodium Halide Interfaces. I. Flexible and Polarizable Molecular Modeling and Interfacial Properties , 2007 .

[53]  Peter Itskowitz,et al.  Chemical Potential Equalization Principle: Direct Approach from Density Functional Theory , 1997 .

[54]  Haiyan Liu,et al.  Molecular dynamics simulations of liquid methanol and methanol–water mixtures with polarizable models , 2006, J. Comput. Chem..

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

[56]  Alexander D. MacKerell,et al.  Determination of Electrostatic Parameters for a Polarizable Force Field Based on the Classical Drude Oscillator. , 2005, Journal of chemical theory and computation.

[57]  Mark S. Gordon,et al.  An interpretation of the enhancement of the water dipole moment due to the presence of other water molecules. , 2008, The journal of physical chemistry. A.

[58]  Sandeep Patel,et al.  A nonadditive methanol force field: bulk liquid and liquid-vapor interfacial properties via molecular dynamics simulations using a fluctuating charge model. , 2005, The Journal of chemical physics.

[59]  L. Dang,et al.  The nonadditive intermolecular potential for water revised , 1992 .

[60]  Charles L. Brooks,et al.  CHARMM fluctuating charge force field for proteins: I parameterization and application to bulk organic liquid simulations , 2004, J. Comput. Chem..

[61]  Hydrogen-bond dynamics in the air-water interface. , 2005, The journal of physical chemistry. B.

[62]  M. Berkowitz,et al.  Ion solvation in water clusters , 1993 .

[63]  Ruhong Zhou,et al.  Parametrizing a polarizable force field from ab initio data. I. The fluctuating point charge model , 1999 .

[64]  A. Pohorille,et al.  Surface potential of the water liquid-vapor interface. , 1988, The Journal of chemical physics.

[65]  Nohad Gresh,et al.  Energetics of Zn2+ binding to a series of biologically relevant ligands: A molecular mechanics investigation grounded on ab initio self‐consistent field supermolecular computations , 1995, J. Comput. Chem..

[66]  Alexander D. MacKerell,et al.  A simple polarizable model of water based on classical Drude oscillators , 2003 .

[67]  Sagar A. Pandit,et al.  Mixed bilayer containing dipalmitoylphosphatidylcholine and dipalmitoylphosphatidylserine: lipid complexation, ion binding, and electrostatics. , 2003, Biophysical journal.

[68]  Nohad Gresh,et al.  Modeling of Copper(II) Complexes with the SIBFA Polarizable Molecular Mechanics Procedure. Application to a New Class of HIV-1 Protease Inhibitors , 2003 .

[69]  P. Procacci,et al.  A transferable polarizable electrostatic force field for molecular mechanics based on the chemical potential equalization principle , 2002 .

[70]  T. R. Dyke,et al.  Partially deuterated water dimers: Microwave spectra and structure , 1980 .

[71]  Pengyu Y. Ren,et al.  Towards accurate solvation dynamics of divalent cations in water using the polarizable amoeba force field: From energetics to structure. , 2006, The Journal of chemical physics.

[72]  Bin Chen,et al.  Structure and dynamics of the aqueous liquid-vapor interface: a comprehensive particle-based simulation study. , 2006, The journal of physical chemistry. B.

[73]  Michiel Sprik,et al.  Hydrogen bonding and the static dielectric constant in liquid water , 1991 .

[74]  Thomas F. Miller,et al.  Quantum diffusion in liquid water from ring polymer molecular dynamics. , 2005, The Journal of chemical physics.

[75]  G. Loglio,et al.  Temperature dependence of the surface tension of water by the equilibrium ring method , 1972 .

[76]  Alexander D. MacKerell,et al.  Additive and Classical Drude Polarizable Force Fields for Linear and Cyclic Ethers. , 2007, Journal of chemical theory and computation.

[77]  Steven J. Stuart,et al.  Fluctuating charge force fields for aqueous solutions , 1995 .

[78]  B. Roux,et al.  Absolute hydration free energy scale for alkali and halide ions established from simulations with a polarizable force field. , 2006, The journal of physical chemistry. B.

[79]  L. Dang,et al.  MOLECULAR DYNAMICS STUDY OF WATER CLUSTERS, LIQUID, AND LIQUID-VAPOR INTERFACE OF WATER WITH MANY-BODY POTENTIALS , 1997 .

[80]  Michele Parrinello,et al.  Water Molecule Dipole in the Gas and in the Liquid Phase , 1999 .

[81]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[82]  Dean R. Haeffner,et al.  Electron distribution in water , 2000 .

[83]  Nohad Gresh,et al.  Comparative binding energetics of Mg2+, Ca2+, Zn2+, and Cd2+ to biologically relevant ligands: Combined ab initio SCF supermolecule and molecular mechanics investigation , 1996, J. Comput. Chem..

[84]  Kazimierz Krynicki,et al.  Pressure and temperature dependence of self-diffusion in water , 1978 .

[85]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[86]  Y. Shen,et al.  Surface properties probed by second-harmonic and sum-frequency generation , 1989, Nature.

[87]  D. Tildesley,et al.  The free surface of water: molecular orientation, surface potential and nonlinear susceptibility , 1997 .

[88]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .

[89]  Michael W. Mahoney,et al.  A five-site model for liquid water and the reproduction of the density anomaly by rigid, nonpolarizable potential functions , 2000 .

[90]  A. Panagiotopoulos,et al.  Gibbs ensemble Monte Carlo simulations of coexistence properties of a polarizable potential model of water , 2002 .

[91]  Y. Fukuda,et al.  Observation of ultrafast spin-lattice relaxation in Tm2+-doped CaF2 and SrF2 crystals by optical means , 2000 .

[92]  Akihiro Morita.  Water polarizability in condensed phase: Ab initio evaluation by cluster approach , 2002, J. Comput. Chem..

[93]  Pengyu Y. Ren,et al.  Ion solvation thermodynamics from simulation with a polarizable force field. , 2003, Journal of the American Chemical Society.

[94]  Yang Zhong,et al.  Thermodynamic and structural properties of methanol–water solutions using nonadditive interaction models , 2008, J. Comput. Chem..

[95]  P. Shepson,et al.  The role of Br2 and BrCl in surface ozone destruction at polar sunrise. , 2001, Science.

[96]  R. Saykally,et al.  On the nature of ions at the liquid water surface. , 2006, Annual review of physical chemistry.

[97]  Sandeep Patel,et al.  Structure, thermodynamics, and liquid-vapor equilibrium of ethanol from molecular-dynamics simulations using nonadditive interactions. , 2005, The Journal of chemical physics.

[98]  Peter A. Kollman,et al.  Ion solvation in polarizable water: molecular dynamics simulations , 1991 .

[99]  D. Tobias,et al.  Propensity of soft ions for the air/water interface , 2004 .

[100]  R. Saykally,et al.  Confirmation of enhanced anion concentration at the liquid water surface , 2004 .

[101]  P. Senet,et al.  A Hirshfeld partitioning of polarizabilities of water clusters. , 2006, The Journal of chemical physics.

[102]  A. Soper,et al.  Joint structure refinement of x-ray and neutron diffraction data on disordered materials: application to liquid water , 2007, Journal of physics. Condensed matter : an Institute of Physics journal.

[103]  M. Matsumoto,et al.  Study on liquid–vapor interface of water. I. Simulational results of thermodynamic properties and orientational structure , 1988 .

[104]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[105]  Laurence S. Rothman,et al.  Dipole moment of water from Stark measurements of H2O, HDO, and D2O , 1973 .

[106]  S. Rick.  Simulations of ice and liquid water over a range of temperatures using the fluctuating charge model , 2001 .

[107]  Johann Gasteiger,et al.  Electronegativity equalization: application and parametrization , 1985 .

[108]  Bruce J. Berne,et al.  Dynamical Fluctuating Charge Force Fields: The Aqueous Solvation of Amides , 1996 .

[109]  M. Paluch.  Electrical properties of free surface of water and aqueous solutions. , 2000, Advances in colloid and interface science.

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

[111]  P. Kollman,et al.  Structure and Properties of Neat Liquids Using Nonadditive Molecular Dynamics: Water, Methanol, and N-Methylacetamide , 1995 .

[112]  Paul Tavan,et al.  The polarizability of point-polarizable water models: density functional theory/molecular mechanics results. , 2008, The journal of physical chemistry. B.