J. L. Aragones

Anomalies in water as obtained from computer simulations of the TIP4P/2005 model: density maxima, and density, isothermal compressibility and heat capacity minima

The so-called thermodynamic anomalies of water form an integral part of the peculiar behaviour of this both important and ubiquitous molecule. In this paper our aim is to establish whether the recently proposed TIP4P/2005 model is capable of reproducing a number of these anomalies. Using molecular dynamics simulations we investigate both the maximum in density and the minimum in the isothermal compressibility along a number of isobars. It is shown that the model correctly describes the decrease in the temperature of the density maximum with increasing pressure. At atmospheric pressure the model exhibits an additional minimum in density at a temperature of about 200K, in good agreement with recent experimental work on super-cooled confined water. The model also presents a minimum in the isothermal compressibility close to 310K. We have also investigated the atmospheric pressure isobar for three other water models; the SPC/E and TIP4P models also present a minimum in the isothermal compressibility, although at a considerably lower temperature than the experimental one. For the temperature range considered no such minimum is found for the TIP5P model.

Solubility of NaCl in water by molecular simulation revisited

In this paper, the solubility of NaCl in water is evaluated by using computer simulations for three different force fields. The condition of chemical equilibrium (i.e., equal chemical potential of the salt in the solid and in the solution) is obtained at room temperature and pressure to determine the solubility of the salt. We used the same methodology that was described in our previous work [E. Sanz and C. Vega, J. Chem. Phys. 126, 014507 (2007)] although several modifications were introduced to improve the accuracy of the calculations. It is found that the predictions of the solubility are quite sensitive to the details of the force field used. Certain force fields underestimate the experimental solubility of NaCl in water by a factor of four, whereas the predictions of other force fields are within 20% of the experimental value. Direct coexistence molecular dynamic simulations were also performed to determine the solubility of the salt. Reasonable agreement was found between the solubility obtained from free energy calculations and that obtained from direct coexistence simulations. This work shows that the evaluation of the solubility of salts in water can now be performed in computer simulations. The solubility depends on the ion-ion, ion-water, and water-water interactions. For this reason, the prediction of the solubility can be quite useful in future work to develop force fields for ions in water.

Anomalies in bulk supercooled water at negative pressure

Significance Water is the most familiar liquid, and arguably the most complex. Anomalies of supercooled water have been measured during decades, and competing interpretations proposed. Yet, a decisive experiment remains elusive, because of unavoidable crystallization into ice. We investigate the state of water that is both supercooled and under mechanical tension, or negative pressure. Liquids under negative pressure can be found in plants or fluid inclusions in minerals. Using such water inclusions in quartz, we report, to our knowledge, the first measurements on doubly metastable water down to −15°C and around −100 MPa. We observe sound velocity anomalies that can be reproduced quantitatively with molecular dynamics simulations. These results suggest the possibility to rule out two proposed scenarios for water anomalies, and put further constraints on the remaining ones. Water anomalies still defy explanation. In the supercooled liquid, many quantities, for example heat capacity and isothermal compressibility κT, show a large increase. The question arises if these quantities diverge, or if they go through a maximum. The answer is key to our understanding of water anomalies. However, it has remained elusive in experiments because crystallization always occurred before any extremum is reached. Here we report measurements of the sound velocity of water in a scarcely explored region of the phase diagram, where water is both supercooled and at negative pressure. We find several anomalies: maxima in the adiabatic compressibility and nonmonotonic density dependence of the sound velocity, in contrast with a standard extrapolation of the equation of state. This is reminiscent of the behavior of supercritical fluids. To support this interpretation, we have performed simulations with the 2005 revision of the transferable interaction potential with four points. Simulations and experiments are in near-quantitative agreement, suggesting the existence of a line of maxima in κT (LMκT). This LMκT could either be the thermodynamic consequence of the line of density maxima of water [Sastry S, Debenedetti PG, Sciortino F, Stanley HE (1996) Phys Rev E 53:6144–6154], or emanate from a critical point terminating a liquid–liquid transition [Sciortino F, Poole PH, Essmann U, Stanley HE (1997) Phys Rev E 55:727–737]. At positive pressure, the LMκT has escaped observation because it lies in the “no man’s land” beyond the homogeneous crystallization line. We propose that the LMκT emerges from the no man’s land at negative pressure.

Properties of ices at 0 K: A test of water models

The properties of ices Ih, II, III, V, and VI at zero temperature and pressure are determined by computer simulation for several rigid water models (SPC/E, TIP5P, TIP4P/Ice, and TIP4P/2005). The energies of the different ices at zero temperature and pressure (relative to the ice II energy) are compared to the experimental results of Whalley [J. Chem. Phys. 81, 4087 (1984)]. TIP4P/Ice and TIP4P/2005 provide a qualitatively correct description of the relative energies of the ices at these conditions. In fact, only these two models provide the correct ordering in energies. For the SPC/E and TIP5P models, ice II is the most stable phase at zero temperature and pressure whereas for TIP4P/Ice and TIP4P/2005 ice Ih is the most stable polymorph. These results are in agreement with the relative stabilities found at higher temperatures. The solid-solid phase transitions at 0 K are determined. The predicted pressures are in good agreement with those obtained from free energy calculations.

Plastic crystal phases of simple water models.

We report the appearance of two plastic crystal phases of water at high pressure and temperature using computer simulations. In one of them the oxygen atoms form a body centered cubic structure (bcc) and in the other they form a face centered cubic structure (fcc). In both cases the water molecules were able to rotate almost freely. We have found that the bcc plastic crystal transformed into a fcc plastic crystal via a Martensitic phase transition when heated at constant pressure. We have performed the characterization and localization in the phase diagram of these plastic crystal phases for the SPC/E, TIP4P, and TIP4P/2005 water potential models. For TIP4P/2005 model free energy calculations were carried out for the bcc plastic crystal and fcc plastic crystal using a new method (which is a slight variation of the Einstein crystal method) proposed for these types of solid. The initial coexistence points for the SPC/E and TIP4P models were obtained using Hamiltonian Gibbs-Duhem integration. For all of these models these two plastic crystal phases appear in the high pressure and temperature region of the phase diagram. It would be of interest to study if such plastic crystal phases do indeed exist for real water. This would shed some light on the question of whether these models can describe satisfactorily the high pressure part of the phase diagram of water, and if not, where and why they fail.

Computer simulation study of the structure of LiCl aqueous solutions: test of non-standard mixing rules in the ion interaction.

Aqueous solutions of LiCl have recently received much attention in connection with the study of the anomalies of supercooled water and its polyamorphism. From the point of view of computer simulation, there is need for a force field that can reproduce the structural and dynamical properties of this solution, and more importantly it is also simple enough to use in large scale simulations of supercooled states. We study by molecular dynamics the structure of the LiCl-water solutions with the force field proposed by Joung and Cheatham (J. Phys. Chem. B 2008, 112, 9020) appropriate for the water TIP4P-Ew model potential. We found that this force field does not reproduce the experimental ion pairing when the Lorentz-Berthelot (LB) rules are used. By incorporating deviations to the LB rules to obtain the crossed interactions between the ions, it is possible to get agreement with experiment. We have studied how the modification of the LB rule affects the structural and thermodynamic properties of the solution at increasing concentration of the solution from the low (around 2%) to medium (around 14%) concentration regimes. We also tested the transferability of the Joung and Cheatham force field to the water TIP4P/2005 model that works very well for supercooled water.

Free energy calculations for molecular solids using GROMACS

In this work, we describe a procedure to evaluate the free energy of molecular solids with the GROMACS molecular dynamics package. The free energy is calculated using the Einstein molecule method that can be regarded as a small modification of the Einstein crystal method. Here, the position and orientation of the molecules is fixed by using an Einstein field that binds with harmonic springs at least three non-collinear atoms (or points of the molecule) to their reference positions. The validity of the Einstein field is tested by performing free-energy calculations of methanol, water (ice), and patchy colloids molecular solids. The free energies calculated with GROMACS show a very good agreement with those obtained using Monte Carlo and with previously published results.

The phase diagram of water from quantum simulations.

The phase diagram of water has been calculated from the TIP4PQ/2005 model, an empirical rigid non-polarisable model. The path integral Monte Carlo technique was used, permitting the incorporation of nuclear quantum effects. The coexistence lines were traced out using the Gibbs-Duhem integration method, once having calculated the free energies of the liquid and solid phases in the quantum limit, which were obtained via thermodynamic integration from the classical value by scaling the mass of the water molecule. The resulting phase diagram is qualitatively correct, being displaced to lower temperatures by 15-20 K. It is found that the influence of nuclear quantum effects is correlated to the tetrahedral order parameter.

Note: Free energy calculations for atomic solids through the Einstein crystal/molecule methodology using GROMACS and LAMMPS

In this work the free energy of solid phases is computed for the Lennard-Jones potential and for a model of NaCl. The free energy is evaluated through the Einstein crystal/molecule methodologies using the molecular dynamics programs: GROMACS and LAMMPS. The obtained results are compared with the results obtained from Monte Carlo. Good agreement between the different programs and methodologies was found. The procedure to perform the free energy calculations for the solid phase in the molecular dynamic programs is described. Since these programs allow to study any continuous intermolecular potential (when given in a tabulated form) this work shows that for isotropic potentials (describing for instance atomic solids or colloidal particles) free energy calculations can be performed on a routinely basis using GROMACS and/or LAMMPS.

Calculation of the melting point of alkali halides by means of computer simulations

In this paper, we study the liquid-solid coexistence of NaCl-type alkali halides, described by interaction potentials such as Tosi-Fumi (TF), Smith-Dang (SD), and Joung-Cheatham (JC), and compute their melting temperature (T(m)) at 1 bar via three independent routes: (1) liquid/solid direct coexistence, (2) free-energy calculations, and (3) Hamiltonian Gibbs-Duhem integration. The melting points obtained by the three routes are consistent with each other. The calculated T(m) of the Tosi-Fumi model of NaCl is in good agreement with the experimental value as well as with other numerical calculations. However, the other two models considered for NaCl, SD and JC, overestimate the melting temperature of NaCl by more than 200 K. We have also computed the melting temperature of other alkali halides using the Tosi-Fumi interaction potential and observed that the predictions are not always as close to the experimental values as they are for NaCl. It seems that there is still room for improvement in the area of force-fields for alkaline halides, given that so far most models are still unable to describe a simple yet important property such as the melting point.