Density maximum and polarizable models of water

Melting points of water models: Current situation

By using the direct coexistence method, we have calculated the melting points of ice I h at normal pressure for three recently proposed water models, namely, TIP3P-FB, TIP4P-FB, and TIP4P-D. We obtained T m = 216 K for TIP3P-FB, T m = 242 K for TIP4P-FB, and T m = 247 K for TIP4P-D. We revisited the melting point of TIP4P/2005 and TIP5P obtaining T m = 250 and 274 K, respectively. We summarize the current situation of the melting point of ice I h for a number of water models and conclude that no model is yet able to simultaneously reproduce the melting temperature of ice I h and the temperature of the maximum in density at room pressure. This probably points toward our both still incomplete knowledge of the potential energy surface of water and the necessity of incorporating nuclear quantum effects to describe both properties simultaneously.

Maximum in density of electrolyte solutions: Learning about ion-water interactions and testing the Madrid-2019 force field

In this work, we studied the effect of Li⁺, Na⁺, K⁺, Mg²⁺, and Ca²⁺ chlorides and sulfates on the temperature of maximum density (TMD) of aqueous solutions at room pressure. Experiments at 1 molal salt concentration were carried out to determine the TMD of these solutions. We also performed molecular dynamics simulations to estimate the TMD at 1 and 2 m with the Madrid-2019 force field, which uses the TIP4P/2005 water model and scaled charges for the ions, finding an excellent agreement between experiment and simulation. All the salts studied in this work shift the TMD of the solution to lower temperatures and flatten the density vs temperature curves (when compared to pure water) with increasing salt concentration. The shift in the TMD depends strongly on the nature of the electrolyte. In order to explore this dependence, we have evaluated the contribution of each ion to the shift in the TMD concluding that Na⁺, Ca²⁺, and SO₄ ^2- seem to induce the largest changes among the studied ions. The volume of the system has been analyzed for salts with the same anion and different cations. These curves provide insight into the effect of different ions upon the structure of water. We claim that the TMD of electrolyte solutions entails interesting physics regarding ion-water and water-water interactions and should, therefore, be considered as a test property when developing force fields for electrolytes. This matter has been rather unnoticed for almost a century now and we believe it is time to revisit it.

Polarizable TIP7P Water Model with Perturbation Charges Evaluated from ABEEM

A polarizable version of the rigid seven-site (TIP7P) water model with the atom-bond electronegativity equalization method (ABEEM) is proposed. The model uses direct polarization, where an isolated water monomer in the equilibrium geometry is assumed as a reference state and the polarization of the monomer arises from interacting with other molecules as a perturbation of the reference state. The charge on each site of the monomer splits into reference charge and perturbation charge. The perturbation charge arises only because of other reference charges. The interaction of the perturbation charge with other perturbation charges is replaced using polarization scaling to enhance the interaction of perturbation charge with the reference charges of the sites from other molecules. The perturbation charges are updated by evaluating explicit expressions once. This direct polarization is time-reversible because the charge update is independent of the charges in previous simulation steps. A Slater-type damping function moderates the short-range electrostatics to treat charge diffusion. The Ewald method corrects the long-range electrostatics both in the nuclei movement and in electronegativity equalization to diminish the size effect. The water model is parameterized by fitting the ab initio results of water clusters and the experimental results of water monomers and thermodynamic properties for liquid water. Owing to polarizability, the model performs better than the TIP7P model in terms of vaporization enthalpy, isothermal compressibility, and shear viscosity of the liquid phase. It performs better at the melting point of ice but slightly worse under critical conditions than the TIP7P model. Direct polarization has a low time complexity of O(N) and is a good choice for ABEEM to improve its computational efficiency.

Hydration and Dispersion Forces in Hydroxypropylcellulose Phase Behavior

Many-body polarization and hydration forces can strongly affect the equilibrium structure and energetics of mixed phases. Accurately reproducing both forces presents a challenge to force field models because it requires balancing hydrogen bonding at short range with many-body orientational order and dispersive attraction at long range. This work reports the first comparison of experimental measurements of the pressure-area isotherm for hydroxypropylcellulose (HPC) against molecular dynamics results with four different force field models-united-atom, all-atom (OPLS and CHARMM), and Drude oscillator models. All force fields exhibit the experimentally determined, exponentially shaped repulsive force at short range. Above a critical temperature of about 40 °C and a lattice spacing of around 14 Å, HPC experiments show a reversible, heat-induced polymer aggregation into an ordered phase driven by loss of water. The nonpolarizable force fields do not display the critical point and instead show biphasic behavior at all temperatures tested. This indicates net attractive forces at intermediate lattice spacings. In contrast, the Drude polarizable force field shows positive osmotic pressure and a single, homogeneous phase over all temperatures and spacings tested. Analysis of structural data from our simulations provides several clues to help interpret these findings. Although all force fields show similar water-water hydrogen bond numbers in the mixed phase, the polarizable model predicts that water-HPC hydrogen bonds are much more favorable than HPC-HPC hydrogen bonds when polymers are dispersed. At high density, water is driven out and replaced by HPC-HPC hydrogen bonds. The polarizable force field shows that both effects have a stronger dependence on polymer density than any of the nonpolarizable models. Our observations support the conclusion that hydration forces are coupled to the polymer coordination number by local, structural waters and that long-range dispersive attraction is overestimated by pairwise additive models.

Polarizable Drude Model with s-Type Gaussian or Slater Charge Density for General Molecular Mechanics Force Fields

Gas-phase electric properties of molecules can be computed routinely using wave function methods or density functional theory (DFT). However, these methods remain computationally expensive for high-throughput screening of the vast chemical space of virtual compounds. Therefore, empirical force fields are a more practical choice in many cases, particularly since force field methods allow one to routinely predict the physicochemical properties in the condensed phases. This work presents Drude polarizable models, to increase the physical realism in empirical force fields, where the core particle is treated as a point charge and the Drude particle is treated either as a 1 s-Gaussian or a ns-Slater ( n = 1, 2, 3) charge density. Systematic parametrization to large high-quality quantum chemistry data obtained from the open access Alexandria Library ( https://doi.org/10.5281/zenodo.1004711 ) ensures the transferability of these parameters. The dipole moments and isotropic polarizabilities of the isolated molecules predicted by the proposed Drude models are in agreement with experiment with accuracy similar to DFT calculations at the B3LYP/aug-cc-pVTZ level of theory. The results show that the inclusion of explicit polarization into the models reduces the root-mean-square deviation with respect to DFT calculations of the predicted dipole moments of 152 dimers and clusters by more than 50%. Finally, we show that the accuracy of the electrostatic interaction energy of the water dimers can be improved systematically by the introduction of polarizable smeared charges as a model for charge penetration.

Hydrothermal properties of the COS/D2 water model: a polarizable charge-on-spring water model, at elevated temperatures and pressures

Molecular simulations have been conducted to assess the pVT properties and static permittivity of the charge-on-spring polarizable water model COS/D2 at hydrothermal conditions from 300 to 450 °C and bulk densities of 0.001 to 1.0 g cm⁻³.

Polarizable six-point water models from computational and empirical optimization

Tröster et al. (J. Phys. Chem B 2013, 117, 9486-9500) recently suggested a mixed computational and empirical approach to the optimization of polarizable molecular mechanics (PMM) water models. In the empirical part the parameters of Buckingham potentials are optimized by PMM molecular dynamics (MD) simulations. The computational part applies hybrid calculations, which combine the quantum mechanical description of a H2O molecule by density functional theory (DFT) with a PMM model of its liquid phase environment generated by MD. While the static dipole moments and polarizabilities of the PMM water models are fixed at the experimental gas phase values, the DFT/PMM calculations are employed to optimize the remaining electrostatic properties. These properties cover the width of a Gaussian inducible dipole positioned at the oxygen and the locations of massless negative charge points within the molecule (the positive charges are attached to the hydrogens). The authors considered the cases of one and two negative charges rendering the PMM four- and five-point models TL4P and TL5P. Here we extend their approach to three negative charges, thus suggesting the PMM six-point model TL6P. As compared to the predecessors and to other PMM models, which also exhibit partial charges at fixed positions, TL6P turned out to predict all studied properties of liquid water at p0 = 1 bar and T0 = 300 K with a remarkable accuracy. These properties cover, for instance, the diffusion constant, viscosity, isobaric heat capacity, isothermal compressibility, dielectric constant, density, and the isobaric thermal expansion coefficient. This success concurrently provides a microscopic physical explanation of corresponding shortcomings of previous models. It uniquely assigns the failures of previous models to substantial inaccuracies in the description of the higher electrostatic multipole moments of liquid phase water molecules. Resulting favorable properties concerning the transferability to other temperatures and conditions like the melting of ice are also discussed.

A systematic development of a polarizable potential of water

Based on extensive studies of existing potentials we propose a new molecular model for water. The new model is rigid and contains three Gaussian charges. Contrary to other models, all charges take part in the polarization of the molecule. They are connected by harmonic springs to their gas-phase positions: the negative charge to a prescribed point on the main axis of the molecule; the positive charges to the hydrogens. The mechanical equilibrium between the electrostatic forces and the spring forces determines the polarization of the molecule which is established by iteration at every timestep. The model gives excellent estimates for ambient liquid properties and reasonably good results from high-pressure solids to gas-phase clusters. We present a detailed description of the development of this model and a large number of calculated properties compared to the estimates of the nonpolarizable TIP4P∕2005 [J. L. F. Abascal and C. Vega, J. Chem. Phys. 123, 234505 (2005)], the polarizable GCPM [P. Paricaud, M. Predota, A. A. Chialvo, and P. T. Cummings, J. Chem. Phys. 122, 244511 (2005)], and our earlier BKd3 model [P. T. Kiss and A. Baranyai, J. Chem. Phys. 137, 084506 (2012)]. The best overall performance is shown by the new model.

The Microscopic Physical Cause for the Density Maximum of Liquid Water

The existence of a density maximum at 277 K is probably the most prominent anomaly among the many very special thermodynamic properties of liquid water. While usually attributed to so-called hydrogen bonding, the microscopic physical cause of this prominent anomaly is still elusive. Here we show that the density anomaly is caused by those short-range electrostatic forces, which are generated by the quadrupole and higher moments of the charge distributions present in liquid-phase water molecules. This conclusion derives from 20 ns replica exchange molecular-dynamics simulations with closely related polarizable four-, five-, and six-point water models. As soon as the model complexity suffices to represent the higher electrostatic moments with sufficient accuracy, the density temperature profile n(T) calculated for T ∈ [250,320] K at the standard pressure 1 bar locks in to the experimental observation. The corresponding six-point model is, therefore, the most simple available cartoon for liquid-phase water molecules.

Systematic improvement of a classical molecular model of water

We report the iAMOEBA ("inexpensive AMOEBA") classical polarizable water model. The iAMOEBA model uses a direct approximation to describe electronic polarizability, in which the induced dipoles are determined directly from the permanent multipole electric fields and do not interact with one another. The direct approximation reduces the computational cost relative to a fully self-consistent polarizable model such as AMOEBA. The model is parameterized using ForceBalance, a systematic optimization method that simultaneously utilizes training data from experimental measurements and high-level ab initio calculations. We show that iAMOEBA is a highly accurate model for water in the solid, liquid, and gas phases, with the ability to fully capture the effects of electronic polarization and predict a comprehensive set of water properties beyond the training data set including the phase diagram. The increased accuracy of iAMOEBA over the fully polarizable AMOEBA model demonstrates ForceBalance as a method that allows the researcher to systematically improve empirical models by efficiently utilizing the available data.

Polarizable water models from mixed computational and empirical optimization

Here we suggest a mixed computational and empirical approach serving to optimize the parameters of complex and polarizable molecular mechanics (PMM) models for complicated liquids. The computational part of the parameter optimization relies on hybrid calculations combining density functional theory (DFT) for a solute molecule with a PMM treatment of its solvent environment at well-defined thermodynamic conditions. As an application we have developed PMM models for water featuring ν = 3, 4, and 5 points of force action, a Gaussian inducible dipole and a Buckingham potential at the oxygen, the experimental liquid phase geometry, the experimental gas phase polarizability α(exp)(g) = 1.47 ų, and, for ν = 4 and 5, the gas phase value μ(exp)(g) = 1.855 D for the static dipole moment. The widths of the Gaussian dipoles and, for ν = 4 and 5, also the electrostatic geometries of these so-called TLνP models are derived from self-consistent DFT/PMM calculations, and the parameters of the Buckingham potentials (and the static TL3P dipole moment) are estimated from molecular dynamics (MD) simulations. The high quality of the resulting models is demonstrated for the observables targeted during optimization (potential energy per molecule, pressure, radial distribution functions) and a series of predicted properties (quadrupole moments, density at constant pressure, dielectric constant, diffusivity, viscosity, compressibility, heat capacity) at certain standard conditions. Remaining deficiencies and possible ways for their removal are discussed.

Electronically coarse-grained model for water

We introduce an electronically coarse-grained description of water representing all long range, many-body electronic responses via an embedded quantum oscillator. Leading-order response coefficients and gas phase electrostatic moments are exactly reproduced. Molecular dynamics, using electronic path integral sampling, shows that this framework is sufficient for a realistic liquid to emerge naturally with transferability extending further to nonambient state points and to the free water surface. The model allows the strength of many-body dispersion and polarization to be adjusted independently and these are found to have significant effects on the condensed phase.

Six-site polarizable model of water based on the classical Drude oscillator

A polarizable water model, SWM6, was developed and optimized for liquid phase simulations under ambient conditions. Building upon the previously developed SWM4-NDP model, additional sites representing oxygen lone-pairs were introduced. The geometry of the sites is assumed to be rigid. Considering the large number of adjustable parameters, simulated annealing together with polynomial fitting was used to facilitate model optimization. The new water model was shown to yield the correct self-diffusion coefficient after taking the system size effect into account, and the dimer geometry is better reproduced than in the SWM4 models. Moreover, the experimental oxygen-oxygen radial distribution is better reproduced, indicating that the new model more accurately describes the local hydrogen bonding structure of bulk phase water. This was further validated by its ability to reproduce the experimental nuclear magnetic shielding and related chemical shift of the water hydrogen in the bulk phase, a property sensitive to the local hydrogen bonding structure. In addition, comparison of the liquid properties of the SWM6 model is made with those of a number of widely used additive and polarizable models. Overall, improved balance between the description of monomer, dimer, clustered, and bulk phase water is obtained with the new model compared to its SWM4-NDP polarizable predecessor, though application of the model requires an approximately twofold increase on computational resources.