The Stability Limit and other Open Questions on Water at Negative Pressure

Stability and Metastability of Liquid Water in a Machine-Learned Coarse-Grained Model with Short-Range Interactions

Coarse-grained water models are ∼100 times more efficient than all-atom models, enabling simulations of supercooled water and crystallization. The machine-learned monatomic model ML-BOP reproduces the experimental equation of state (EOS) and ice-liquid thermodynamics at 0.1 MPa on par with the all-atom TIP4P/2005 and TIP4P/Ice models. These all-atom models were parametrized using high-pressure experimental data and are either accurate for water's EOS (TIP4P/2005) or ice-liquid equilibrium (TIP4P/Ice). ML-BOP was parametrized from temperature-dependent ice and liquid experimental densities and melting data at 0.1 MPa; its only pressure training is from compression of TIP4P/2005 ice at 0 K. Here we investigate whether ML-BOP replicates the experimental EOS and ice-water thermodynamics along all pressures of ice I. We find that ML-BOP reproduces the temperature, enthalpy, entropy, and volume of melting of hexagonal ice up to 400 MPa and the EOS of water along the melting line with an accuracy that rivals that of both TIP4P/2005 and TIP4P/Ice. We interpret that the accuracy of ML-BOP originates from its ability to capture the shift between compact and open local structures to changes in pressure and temperature. ML-BOP reproduces the sharpening of the tetrahedral peak of the pair distribution function of water upon supercooling, and its pressure dependence. We characterize the region of metastability of liquid ML-BOP with respect to crystallization and cavitation. The accessibility of ice crystallization to simulations of ML-BOP, together with its accurate representation of the thermodynamics of water, makes it promising for investigating the interplay between anomalies, glass transition, and crystallization under conditions challenging to access through experiments.

Polymer-Water Interaction Enabled Intelligent Moisture Regulation in Hydrogels

The water-vapor transition is critical for hydrogels in a collection of applications. However, how the polymer-water interaction along with the nature of the structure affect the macroscopic water-vapor transition remains a challenging question to answer. In this work, we tested the moisture transfer behaviors of a series of hydrogels at different humidities and found some hydrogels capable of lowering their surface vapor pressure to stop dehydration at low humidity and absorbing water from ambient air to recover toward initial states at high humidity. Through molecular dynamic simulations, we demonstrate that water inside these hydrogels undergoes increasing intensive intermolecular bonding during evaporation. The increased intermolecular bonding reduces the vapor pressure of the hydrogels and leads to the self-regulation. More interestingly, we demonstrate the self-regulation is closely related to the Young's modulus of hydrogels. These results provide further insight into the mechanism of the water-vapor transition in hydrogels and show potential in a broad range of future applications.

The effects of ultrasound pressure and temperature fields in millisecond bubble nucleation

A phenomenological implementation of Classical Nucleation Theory (CNT) is employed to investigate the connection between high intensity focused ultrasound (HIFU) pressure and temperature fields with the energetic requirements of bubble nucleation. As a case study, boiling histotripsy in tissue-mimicking phantoms is modelled. The physics of key components in the implementation of CNT in HIFU conditions such as the derivation of nucleation pressure thresholds and approximations regarding the surface tension of the liquid are reviewed and discussed. Simulations show that the acoustic pressure is the ultimate trigger for millisecond bubble nucleation in boiling histotripsy, however, HIFU heat deposition facilitates nucleation by lowering nucleation pressure thresholds. Nucleation thus occurs preferentially at the regions of highest heat deposition within the HIFU field. This implies that bubble nucleation subsequent to millisecond HIFU heat deposition can take place at temperatures below 100 °C as long as the focal HIFU peak negative pressure exceeds the temperature-dependent nucleation threshold. It is also found that the magnitude of nucleation pressure thresholds decreases with decreasing frequencies. Overall, results indicate that it is not possible to separate thermal and mechanical effects of HIFU in the nucleation of bubbles for timescales of a few milliseconds. This methodology provides a promising framework for studying time and space dependencies of the energetics of bubble nucleation within a HIFU field.

How Solutes Modify the Thermodynamics and Dynamics of Filling and Emptying in Extreme Ink-Bottle Pores

We investigate the filling and emptying of extreme ink-bottle porous media-micrometer-scale pores connected by nanometer-scale pores-when changing the pressure of the external vapor, in a case where the pore liquid contains solutes. These phenomena are relevant in diverse contexts, such as the weathering of building materials and artwork, aerosol formation in the atmosphere, and the hydration of soils and plants. Using model systems made of vein-shaped microcavities interconnected by a mesoporous matrix, we show experimentally that the presence of a nonvolatile solute shifts the condensation and evaporation transitions and in a way that is consistent with a modified Kelvin-Laplace equation that takes into account the osmotic pressure of the solution. Emptying occurs far below saturation, when the Kelvin stress, mediated by the large curvature of the liquid-vapor interfaces in the nanopores, is negative enough to induce spontaneous bubble nucleation in the microveins. Filling, on the other hand, occurs close to equilibrium (i.e., at saturation, p_sat for pure water and p_s < p_sat for a solution), driven by the weak capillary pressure of the liquid-vapor interface in the microveins. Interestingly, solutes allow the system to reach situations where the vapor is supersaturated with respect to the solution ( p_s < p < p_sat). We show that in that latter situation, a condensation layer covers the outer surface of the porous system, preventing the generation of Kelvin stresses but inducing osmotic stresses and flows that are vapor pressure-dependent. The timescales and dynamics reflect these different driving forces: emptying proceeds through discrete, stochastic nucleation events with very fast, unsteady bubble growth associated with a poroelastic relaxation process, while filling occurs collectively in all veins of the sample through a slower steady-state process driven by a combination of osmosis and capillarity. The dynamics can however be rendered symmetrical between filling and emptying if bubbles pre-exist during emptying, a case that we explore using cycling of the vapor pressure around equilibrium.

Ice nucleation triggered by negative pressure

Homogeneous ice nucleation needs supercooling of more than 35 K to become effective. When pressure is applied to water, the melting and the freezing points both decrease. Conversely, melting and freezing temperatures increase under negative pressure, i.e. when water is stretched. This study presents an extrapolation of homogeneous ice nucleation temperatures from positive to negative pressures as a basis for further exploration of ice nucleation under negative pressure. It predicts that increasing negative pressure at temperatures below about 262 K eventually results in homogeneous ice nucleation while at warmer temperature homogeneous cavitation, i. e. bubble nucleation, dominates. Negative pressure occurs locally and briefly when water is stretched due to mechanical shock, sonic waves, or fragmentation. The occurrence of such transient negative pressure should suffice to trigger homogeneous ice nucleation at large supercooling in the absence of ice-nucleating surfaces. In addition, negative pressure can act together with ice-inducing surfaces to enhance their intrinsic ice nucleation efficiency. Dynamic ice nucleation can be used to improve properties and uniformity of frozen products by applying ultrasonic fields and might also be relevant for the freezing of large drops in rainclouds.

A comprehensive scenario of the thermodynamic anomalies of water using the TIP4P/2005 model

The striking behavior of water has deserved it to be referred to as an "anomalous" liquid. The water anomalies are greatly amplified in metastable (supercooled and/or stretched) regions. This makes difficult a complete experimental description since, beyond certain limits, the metastable phase necessarily transforms into the stable one. Theoretical interpretation of the water anomalies could then be based on simulation results of well validated water models. But the analysis of the simulations has not yet reached a consensus. In particular, one of the most popular theoretical scenarios-involving the existence of a liquid-liquid critical point (LLCP)-is disputed by several authors. In this work, we propose to use a number of exact thermodynamic relations which may shed light on this issue. Interestingly, these relations may be tested in a region of the phase diagram which is outside the LLCP thus avoiding the problems associated to the coexistence region. The central property connected to other water anomalies is the locus of temperatures at which the density along isobars attain a maximum (TMD line) or a minimum (TmD). We have performed computer simulations to evaluate the TMD and TmD for a successful water model, namely, TIP4P/2005. We have also evaluated the vapor-liquid (VL) spinodal in the region of large negative pressures. The shape of these curves and their connection to the extrema of some response functions, in particular the isothermal compressibility and heat capacity at constant pressure, provides very useful information which may help to elucidate the validity of the theoretical proposals. In this way, we are able to present for the first time a comprehensive scenario of the thermodynamic water anomalies for TIP4P/2005 and their relation to the vapor-liquid spinodal. The overall picture shows a remarkable similarity with the corresponding one for the ST2 water model, for which the existence of a LLCP has been demonstrated in recent years. It also provides a hint as to where the long-sought for extrema in response functions might become accessible to experiments.

Stability Limit of Water by Metastable Vapor-Liquid Equilibrium with Nanoporous Silicon Membranes

Liquid can sustain mechanical tension as its pressure drops below the vapor-liquid coexistence line and becomes less than zero, until it reaches the stability limit-the pressure at which cavitation inevitably occurs. For liquid water, its stability limit is still a subject of debate: the results obtained by researchers using a variety of techniques show discrepancies between the values of the stability limit and its temperature dependence as temperature approaches 0 °C. In this work, we present a study of the stability limit of water by the metastable vapor-liquid equilibrium (MVLE) method with nanoporous silicon membranes. We also report on an experimental system which enables tests of the temperature dependence of the stability limit with MVLE. The stability limit we found increases monotonically (larger tension) as temperature approaches 0 °C; this trend contradicts the centrifugal result of Briggs but agrees with the experiments by acoustic cavitation. This result confirms that a quasi-static method can reach stability values similar to that from the dynamic stretching technique, even close to 0 °C. Nevertheless, our results fall in the range of ∼ -20 to -30 MPa, a range that is consistent with the majority of experiments but is far less negative than the limit obtained in experiments involving quartz inclusions and that predicted for homogeneous nucleation.

Drying by cavitation and poroelastic relaxations in porous media with macroscopic pores connected by nanoscale throats

We investigate the drying dynamics of porous media with two pore diameters separated by several orders of magnitude. Nanometer-sized pores at the edge of our samples prevent air entry, while drying proceeds by heterogeneous nucleation of vapor bubbles--cavitation--in the liquid in micrometer-sized voids within the sample. We show that the dynamics of cavitation and drying are set by the interplay of the deterministic poroelastic mass transport in the porous medium and the stochastic nucleation process. Spatiotemporal patterns emerge in this unusual reaction-diffusion system, with temporal oscillations in the drying rate and variable roughness of the drying front.

Anomalies in bulk supercooled water at negative pressure

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.