Phase Diagram of TIP4P/2005 Water at High Pressure

Hydrogen-bond linking is crucial for growing ice VII embryos

We use molecular dynamics simulations to examine the homogeneous nucleation of ice VII from metastable liquid water. An unsupervised machine learning classification identifies two distinct local structures composing Ice VII nuclei. The seeding method, combined with the classical nucleation theory (CNT), predicts the solid-liquid interfacial free energy, consistent with the value from the mold integration method. Meanwhile, the nucleation rates estimated from the CNT framework and brute force spontaneous nucleations are inconsistent, and we discuss the reasons for this discrepancy. Structural and dynamical heterogeneities suggest that the potential birthplace for an ice VII embryo is relatively ordered, although not necessarily relatively immobile. Moreover, we demonstrate that without the formation of hydrogen-bond links, ice VII embryos do not grow.

Observation of a Plastic Crystal in Water-Ammonia Mixtures under High Pressure and Temperature

Solid mixtures of ammonia and water, the so-called ammonia hydrates, are thought to be major components of solar and extra-solar icy planets. We present here a thorough characterization of the recently reported high pressure (P)-temperature (T) phase VII of ammonia monohydrate (AMH) using Raman spectroscopy, X-ray diffraction, and quasi-elastic neutron scattering (QENS) experiments in the ranges 4-10 GPa, 450-600 K. Our results show that AMH-VII exhibits common structural features with the disordered ionico-molecular alloy (DIMA) phase, stable above 7.5 GPa at 300 K: both present a substitutional disorder of water and ammonia over the sites of a body-centered cubic lattice and are partially ionic. The two phases however markedly differ in their hydrogen dynamics, and QENS measurements show that AMH-VII is characterized by free molecular rotations around the lattice positions which are quenched in the DIMA phase. AMH-VII is thus a peculiar crystalline solid in that it combines three types of disorder: substitutional, compositional, and rotational.

High-temperature and high-pressure plastic phase of ice at the boundary of liquid water and ice VII

Simultaneous high-temperature and high-pressure studies reveal phase transformation of bulk liquid water to an ice-VII-like structure having an eight coordination. It was demonstrated through this numerical study that the observed high-temperature and high-pressure phase of water obtained upon shock compression and equilibration has high rotational diffusion and thereby the hydrogen dynamics of these crystal structures are significantly complex compared with ice VII. The current work provides new characterization methods for the numerically observed plastic crystal phase of ice at the boundary of the liquid water and ice VII phases in which the molecules have a defined lattice position but rotate freely. It is anticipated that the present work will provide important data and guide new theoretical and experimental investigations in the search for plastic crystal phases of water. The power spectra plots of bulk liquid water subjected to different temperature and pressure conditions have also been presented in this numerical study, demonstrating significant differences between these high-temperature and high-pressure shock-equilibrated phases and those of pure ice VII at 10 GPa and liquid water at ambient temperature and pressure, as well as at elevated pressures and temperatures.

Time-resolved optical probing of the non-equilibrium supercritical state in molecular media under ns laser-plasma impact

We proposed a complex method based on a combination of shadow photography and time-resolved Raman spectroscopy to observe the non-stationary laser-induced supercritical state in molecular media. Shadow photography is applied for retrieving pressure values, while Raman spectroscopy with molecular dynamics for temperature estimation. Time resolution of 0.25 ns is achieved by varying the delay between the pump (creating an extreme energy delivery) and the probe laser pulses by the self-made digital delay electronic circuit . The proposed method was employed in liquid carbon dioxide and water. Under nanosecond laser pulse impact, the estimated temperatures and pressures (∼700 K and ∼0.5 GPa) achieved in media are higher than the critical parameters of the samples.

Insight into Liquid Polymorphism from the Complex Phase Behavior of a Simple Model

We systematically explored the phase behavior of the hard-core two-scale ramp model suggested by Jagla [Phys. Rev. E 63, 061501 (2001)PRESCM1539-375510.1103/PhysRevE.63.061501] using a combination of the nested sampling and free energy methods. The sampling revealed that the phase diagram of the Jagla potential is significantly richer than previously anticipated, and we identified a family of new crystalline structures, which is stable over vast regions in the phase diagram. We showed that the new melting line is located at considerably higher temperature than the boundary between the low- and high-density liquid phases, which was previously suggested to lie in a thermodynamically stable region. The newly identified crystalline phases show unexpectedly complex structural features, some of which are shared with the high-pressure ice VI phase.

Structure and dynamics of water plastic crystals from computer simulations

Water has a rich phase diagram with several crystals, as confirmed by experiments. High-pressure and high-temperature water is of interest for Earth's mantle and exoplanetary investigations. It is in this region of the phase diagram of water that new plastic crystal phases of water have been revealed via computer simulations by both classical forcefields and ab initio calculations. However, these plastic phases still remain elusive in experiments. Here, we present a complete characterization of the structure, dynamics, and thermodynamics of the computational plastic crystal phases of water using molecular dynamics and the two-phase thermodynamic method and uncover the interplay between them. The relaxation times of different reorientational correlation functions are obtained for the hypothetical body-centered-cubic and face-centered-cubic plastic crystal phases of water at T = 440 K and P = 8 GPa. Results are compared to a high pressure liquid and ice VII phases to improve the understanding of the plastic crystal phases. Entropy results indicate that the fcc crystal is more stable compared to the bcc structure under the studied conditions.

Unraveling nucleation pathway in methane clathrate formation

Methane clathrates are widespread on the ocean floor of the Earth. A better understanding of methane clathrate formation has important implications for natural-gas exploitation, storage, and transportation. A key step toward understanding clathrate formation is hydrate nucleation, which has been suggested to involve multiple evolution pathways. Herein, a unique nucleation/growth pathway for methane clathrate formation has been identified by analyzing the trajectories of large-scale molecular dynamics (MD) simulations. In particular, ternary water-ring aggregations (TWRAs) have been identified as fundamental structures for characterizing the nucleation pathway. Based on this nucleation pathway, the critical nucleus size and nucleation timescale can be quantitatively determined. Specifically, a methane hydration layer compression/shedding process is observed to be the critical step in (and driving) the nucleation/growth pathway, which is manifested through overlapping/compression of the surrounding hydration layers of the methane molecules, followed by detachment (shedding) of the hydration layer. As such, an effective way to control methane hydrate nucleation is to alter the hydration layer compression/shedding process during the course of nucleation.

Translational and Rotational Diffusion in Liquid Water at Very High Pressure: A Simulation Study

Translational and rotational diffusion coefficients of liquid water have been computed from molecular dynamics simulation with a recent polarizable potential at 298, 400, and 550 K at very high pressure. At 298 K, the model reproduces the initial increase and the occurrence of a maximum for the translational and rotational diffusion coefficients when the pressure increases. At 400 and 550 K, translational and rotational diffusion coefficients are found to monotonically decrease when pressure increases in the gigapascal range, with the translational coefficient decreasing faster than the rotational one. At 400 K, such an evolution of the rotational diffusion coefficient contrasts with quasielastic neutron scattering results predicting a near independence of the rotational diffusion with a pressure increase above ≃0.5 GPa. An interpretation is proposed to explain this discrepancy. The pressure dependence of the translation-rotation coupling is analyzed. The anisotropy of rotational diffusion is investigated by computing the rotational diffusion tensor in a molecular system of axes and the reorientational correlation times of rank 1 and rank 2 of the inertia axes and of the OH bond vector. Deviation of the simulation data with respect to the predictions of the isotropic Debye model of rotational diffusion are quantified and can be used to estimate experimental rotational diffusion coefficients from experimental reorientational correlation times.

Molecular dynamics simulation of proteins under high pressure: Structure, function and thermodynamics

BACKGROUND

Molecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.

SCOPE OF REVIEW

First, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.

MAJOR CONCLUSIONS

Recent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.

GENERAL SIGNIFICANCE

MD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.

Stepwise Homogeneous Melting of Benzene Phase I at High Pressure

We investigate, using molecular dynamics simulations, the spontaneous homogeneous melting of benzene phase I under a high pressure of 1.0 GPa. We find an apparent stepwise transition via a metastable crystal phase, unlike the direct melting observed at ambient pressure. The transition to the metastable phase is achieved by rotational motions, without the diffusion of the center of mass of benzene. The metastable crystal completely occupies the whole space and maintains its structure for at least several picoseconds, so that the phase seems to have a local free energy minimum. The unit cell is found to be unique—no such crystalline structure has been reported so far. Furthermore, we discuss the influence of pressure control on the melting behavior.

Pressure response of the THz spectrum of bulk liquid water revealed by intermolecular instantaneous normal mode analysis

The radial distribution functions of liquid water are known to change significantly their shape upon hydrostatic compression from ambient conditions deep into the kbar pressure regime. It has been shown that despite their eye-catching changes, the fundamental locally tetrahedral fourfold H-bonding pattern that characterizes ambient water is preserved up to about 10 kbar (1 GPa), which is the stability limit of liquid water at 300 K. The observed increase in coordination number comes from pushing water molecules into the first coordination sphere without establishing an H-bond, resulting in roughly two such additional interstitial molecules at 10 kbar. THz spectroscopy has been firmly established as a powerful experimental technique to analyze H-bonding in aqueous solutions given that it directly probes the far-infrared lineshape and thus the prominent H-bond network mode around 180 cm^-1. We, therefore, set out to assess pressure effects on the THz response of liquid water at 10 kbar in comparison to the 1 bar (0.1 MPa) reference, both at 300 K, with the aim to trace back the related lineshape changes to the structural level. To this end, we employ the instantaneous normal mode approximation to rigorously separate the H-bonding peak from the large background arising from the pronounced librational tail. By exactly decomposing the total molecular dynamics into hindered translations, hindered rotations, and intramolecular vibrations, we find that the H-bonding peak arises from translation-translation and translation-rotation correlations, which are successively decomposed down to the level of distinct local H-bond environments. Our utmost detailed analysis based on molecular pair classifications unveils that H-bonded double-donor water pairs contribute most to the THz response around 180 cm^-1, whereas interstitial waters are negligible. Moreover, short double-donor H-bonds have their peak maximum significantly shifted toward higher frequencies with respect to such long H-bonds. In conjunction with an increasing relative population of these short H-bonds versus the long ones (while the population of other water pair classes is essentially pressure insensitive), this explains not only the blue-shift of the H-bonding peak by about 20-30 cm^-1 in total from 1 bar to 10 kbar but also the filling of the shallow local minimum of the THz lineshape located in between the network peak and the red-wing of the librational band at 1 bar. Based on the changing populations as a function of pressure, we are also able to roughly estimate the pressure-dependence of the H-bond network mode and find that its pressure response and thus the blue-shifting are most pronounced at low kbar pressures.

Rapid freezing of water under dynamic compression

Understanding the behavior of materials at extreme pressures is a central issue in fields like aerodynamics, astronomy, and geology, as well as for advancing technological grand challenges such as inertial confinement fusion. Dynamic compression experiments to probe high-pressure states often encounter rapid phase transitions that may cause the materials to behave in unexpected ways, and understanding the kinetics of these phase transitions remains an area of great interest. In this review, we examine experimental and theoretical/computational efforts to study the freezing kinetics of water to a high-pressure solid phase known as ice VII. We first present a detailed analysis of dynamic compression experiments in which water has been observed to freeze on sub-microsecond time scales to ice VII. This is followed by a discussion of the limitations of currently available molecular and continuum simulation methods in modeling these experiments. We then describe how our phase transition kinetics models, which are based on classical nucleation theory, provide a more physics-based framework that overcomes some of these limitations. Finally, we give suggestions on future experimental and modeling work on the liquid-ice VII transition, including an outline of the development of a predictive multiscale model in which molecular and continuum simulations are intimately coupled.