In this work, we studied the effect of Li+, Na+, K+, Mg2+, and Ca2+ 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+, Ca2+, and [Formula: see text] 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.

One of the most accepted hypothesis to explain the anomalous behavior of water is the presence of a critical point between two liquids, the liquid–liquid critical point (LLCP), buried within the deep supercooled regime. Unfortunately, such hypothesis is hard to be experimentally confirmed due to fast freezing. Here, we show that the TIP4P/Ice water potential shifted by 400 bar can reproduce with unprecedented accuracy the experimental isothermal compressibility of water and its liquid equation of state for a wide pressure and temperature range. We find, both by extrapolation of response function maxima and by a Maxwell construction, that the location of the model LLCP is consistent with previous calculations. According to the pressure shift needed to recover the experimental behavior of supercooled water, we estimate the experimental LLCP to be located around 1250 bar and 195 K. We use the model to estimate the ice nucleation rate (J) in the vicinity of the hypothesized LLCP experimental location and obtain J = 1024 m−3 s−1. Thereby, experiments where the ratio between the cooling rate and the sample volume is equal or larger than the estimated nucleation rate could probe liquid–liquid equilibrium before freezing. Such conditions are not accessible in common experiments with microdroplets cooled at a few kelvin per second, but they could be, for instance, using nanodroplets of around 50 nm radius observed in a millisecond timescale.

No abstract