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.

Electrostatic properties of aqueous salt solution interfaces: a comparison of polarizable and nonpolarizable ion models.

The effects of ion force field polarizability on the interfacial electrostatic properties of approximately 1 M aqueous solutions of NaCl, CsCl, and NaI are investigated using molecular dynamics simulations employing both nonpolarizable and Drude-polarizable ion sets. Differences in computed depth-dependent orientational distributions, "permanent" and induced dipole and quadrupole moment profiles, and interfacial potentials are obtained for both ion sets to further elucidate how ion polarizability affects interfacial electrostatic properties among the various salts relative to pure water. We observe that the orientations and induced dipoles of water molecules are more strongly perturbed in the presence of polarizable ions via a stronger ionic double layer effect arising from greater charge separation. Both anions and cations exhibit enhanced induced dipole moments and strong z alignment in the vicinity of the Gibbs dividing surface (GDS) with the magnitude of the anion induced dipoles being nearly an order of magnitude larger than those of the cations and directed into the vapor phase. Depth-dependent profiles for the trace and z z components of the water molecular quadrupole moment tensors reveal 40% larger quadrupole moments in the bulk phase relative to the vapor which mimics a similar observed 40% increase in the average water dipole moment. Across the GDS, the water molecular quadrupole moments increase nonmonotonically (in contrast to the water dipoles) and exhibit a locally reduced contribution just below the surface due to both orientational and polarization effects. Computed interfacial potentials for the nonpolarizable salts yield values 20-60 mV more positive than pure water and increase by an additional 30-100 mV when ion polarizability is included. A rigorous decomposition of the total interfacial potential into ion monopole, water and ion dipole, and water quadrupole components reveals that a very strong, positive ion monopole contribution is offset by negative contributions from all other potential sources. Water quadrupole components modulated by the water density contribute significantly to the observed interfacial potential increments and almost entirely explain observed differences in the interfacial potentials for the two chloride salts. By lumping all remaining nonquadrupole interfacial potential contributions into a single "effective" dipole potential, we observe that the ratio of quadrupole to "effective" dipole contributions range from 2:1 in CsCl to 1:1.5 in NaI, suggesting that both contributions are comparably important in determining the interfacial potential increments. We also find that oscillations in the quadrupole potential in the double layer region are opposite in sign and partially cancel those of the "effective" dipole potential.

Revised charge equilibration potential for liquid alkanes.

We present a revised liquid alkane force field based on the charge equilibration formalism for incorporating electrostatic nonadditive effects arising from local polarization. The model is a revision of earlier work by Patel and Brooks, specifically addressing deficiencies in the dihedral potential, electrostatic, and Lennard-Jones (van der Waals) parameters of the force field. We discuss refinement of the alkane backbone torsion potential to match high-level ab initio relative conformational energetics for pentane, hexane, and heptane. We further discuss refinement of the electrostatic and Lennard-Jones (van der Waals) parameters to reproduce the experimental polarizability, liquid density, and vaporization enthalpy of hexane. Finally, we calculate bulk liquid properties including densities, vaporization enthalpies, self-diffusion constants, isothermal compressibilities, constant pressure heat capacities, and NMR T 1 relaxation times for a series of linear alkanes ranging from hexane to pentadecane based on the current revised model. We also compute free energies of hydration for pentane, hexane, and heptane. The revised force field offers a significantly improved overall description of these properties relative to the original parametrization. The current alkane force field represents a platform for ongoing development of a CHARMM (Chemistry at Harvard Molecular Mechanics) polarizable force field for lipids and integral membrane proteins.

Thermodynamic and structural properties of methanol-water solutions using nonadditive interaction models.

We study bulk structural and thermodynamic properties of methanol-water solutions via molecular dynamics simulations using novel interaction potentials based on the charge equilibration (fluctuating charge) formalism to explicitly account for molecular polarization at the atomic level. The study uses the TIP4P-FQ potential for water-water interactions, and the CHARMM-based (Chemistry at HARvard Molecular Mechanics) fluctuating charge potential for methanol-methanol and methanol-water interactions. In terms of bulk solution properties, we discuss liquid densities, enthalpies of mixing, dielectric constants, self-diffusion constants, as well as structural properties related to local hydrogen bonding structure as manifested in radial distribution functions and cluster analysis. We further explore the electronic response of water and methanol in the differing local environments established by the interaction of each species predominantly with molecules of the other species. The current force field for the alcohol-water interaction performs reasonably well for most properties, with the greatest deviation from experiment observed for the excess mixing enthalpies, which are predicted to be too favorable. This is qualitatively consistent with the overestimation of the methanol-water gas-phase interaction energy for the lowest-energy conformer (methanol as proton donor). Hydration free energies for methanol in TIP4P-FQ water are predicted to be -5.6 +/- 0.2 kcal/mol, in respectable agreement with the experimental value of -5.1 kcal/mol. With respect to solution microstructure, the present cluster analysis suggests that the microscale environment for concentrations where select thermodynamic quantities reach extremal values is described by a bipercolating network structure.

Hydration free energies of monovalent ions in transferable intermolecular potential four point fluctuating charge water: an assessment of simulation methodology and force field performance and transferability.

Hydration free energies of nonpolarizable monovalent atomic ions in transferable intermolecular potential four point fluctuating charge (TIP4P-FQ) are computed using several commonly employed ion-water force fields including two complete model sets recently developed for use with the simple water model with four sites and Drude polarizability and TIP4P water models. A simulation methodology is presented which incorporates a number of finite-system free energy corrections within the context of constant pressure molecular dynamics simulations employing the Ewald method and periodic boundary conditions. The agreement of the computed free energies and solvation structures with previously reported results for these models in finite droplet systems indicates good transferability of ion force fields from these water models to TIP4Q-FQ even when ion polarizability is neglected. To assess the performance of the ion models in TIP4P-FQ, we compare with consensus values for single-ion hydration free energies arising from recently improved cluster-pair estimates and a reevaluation of commonly cited, experimentally derived single-ion hydration free energies; we couple the observed consistency of these energies with a justification of the cluster-pair approximation in assigning single-ion hydration free energies to advocate the use of these consensus energies as a benchmark set in the parametrization of future ion force fields.

Population size bias in descendant-weighted diffusion quantum Monte Carlo simulations.

We consider the influence of population size on the accuracy of diffusion quantum Monte Carlo simulations that employ descendant weighting or forward walking techniques to compute expectation values of observables that do not commute with the Hamiltonian. We show that for a simple model system, the d-dimensional isotropic harmonic oscillator, the population size must increase rapidly with d in order to ensure that the simulations produce accurate results. When the population size is too small, expectation values computed using descendant-weighted diffusion quantum Monte Carlo simulations exhibit significant systematic biases.