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.

Close-Packed Ices in Nanopores

Water molecules in any of the ice polymorphs organize themselves into a perfect four-coordinated hydrogen-bond network at the expense of dense packing. Even at high pressures, there seems to be no way to reconcile the ice rules with the close packing. Here, we report several close-packed ice phases in carbon nanotubes obtained from molecular dynamics simulations of two different water models. Typically they are in plastic states at high temperatures and are transformed into the hydrogen-ordered ice, keeping their close-packed structures at lower temperatures. The close-packed structures of water molecules in carbon nanotubes are identified with those of spheres in a cylinder. We present design principles of hydrogen-ordered, close-packed structures of ice in nanotubes, which suggest many possible dense ice forms with or without nonzero polarization. In fact, some of the simulated ices are found to exhibit ferroelectric ordering upon cooling.

Electrolyte Permeability in Plastic Ice VII

Deep brines in water-rich planets form when electrolytes diffuse from the rocky interior through layers of thick dense ice such as ice VII and the hypothesized plastic ice VII. We perform classical molecular dynamics simulations of Li⁺, Na⁺, and K⁺ alkali ions and F^- and Cl^- halide ions in plastic ice VII at conditions similar to water-rich super-Earths, icy moons, and ocean worlds. We find that plastic ice VII is permeable to electrolytes on geological timescales. Diffusion occurs via jumps between adjacent voids in the bcc crystal structure and is governed by molecular rotations. An exception to this mechanism is Na⁺ which, at variance with other ions, can substitute water molecules on lattice positions. The bulk modulus of pristine plastic ice VII is dependent on the pace of molecular rotations: when the rotations are slow, the bulk modulus is 1 order of magnitude lower compared to the bulk modulus at conditions of fast rotations, hence providing direct evidence of the role of molecular rotations in determining elastic properties. Electrolytes affect the bulk modulus only at high-concentration conditions and slow molecular rotations. Our results show that plastic ice VII may facilitate the development of brines in water-rich planets and ocean worlds, with a clear significance for their potential to support exobiology and for the chemical evolution of their aqueous reservoirs.

Molecular rotations trigger a glass-to-plastic fcc heterogeneous crystallization in high-pressure water

We report a molecular dynamics study of the heterogeneous crystallization of high-pressure glassy water using (plastic) ice VII as a substrate. We focus on the thermodynamic conditions P ∈ [6-8] GPa and T ∈ [100-500] K, at which (plastic) ice VII and glassy water are supposed to coexist in several (exo)planets and icy moons. We find that (plastic) ice VII undergoes a martensitic phase transition to a (plastic) fcc crystal. Depending on the molecular rotational lifetime τ, we identify three rotational regimes: for τ > 20 ps, crystallization does not occur; for τ ∼ 15 ps, we observe a very sluggish crystallization and the formation of a considerable amount of icosahedral environments trapped in a highly defective crystal or in the residual glassy matrix; and for τ < 10 ps, crystallization takes place smoothly, resulting in an almost defect-free plastic fcc solid. The presence of icosahedral environments at intermediate τ is of particular interest as it shows that such a geometry, otherwise ephemeral at lower pressures, is, indeed, present in water. We justify the presence of icosahedral structures based on geometrical arguments. Our results represent the first study of heterogeneous crystallization occurring at thermodynamic conditions of relevance for planetary science and unveil the role of molecular rotations in achieving it. Our findings (i) show that the stability of plastic ice VII, widely reported in the literature, should be reconsidered in favor of plastic fcc, (ii) provide a rationale for the role of molecular rotations in achieving heterogeneous crystallization, and (iii) represent the first evidence of long-living icosahedral structures in water. Therefore, our work pushes forward our understanding of the properties of water.

Engineering the Thermodynamic Stability and Metastability of Mesophases of Colloidal Bipyramids through Shape Entropy

We report several types of entropy-driven phase transition behaviors in hard bipyramid systems using Monte Carlo simulations. Bipyramidal nanoparticle shapes are synthesizable from gold and silver, with sizes ranging from tens to hundreds of nanometers. We report numerous colloidal crystalline phases with varying symmetries and complexities as the bipyramid aspect ratio and base polygon are varied. Some bipyramids are mesogenic and undergo either monotropic or enantiotropic phase transitions. We show that such mesophase behavior can be modulated by tuning the bipyramid aspect ratio. In addition, we report stepwise kinetic crystallization and melting pathways that occur via an intermediate mesophase as the system gains or loses order in successive stages. Our results demonstrate that complex phase transition behavior involving mesophases can be driven by entropy alone. Importantly, our results can guide the synthesis of bipyramid shapes able to assemble target structures and can be used to engineer the kinetic pathways to and from those structures to involve or avoid mesophases.

Temperature- and pressure-dependence of the hydrogen bond network in plastic ice VII

We model, via classical molecular dynamics simulations, the plastic phase of ice VII across a wide range of the phase diagram of interest for planetary investigations. Although structural and dynamical properties of plastic ice VII are mostly independent on the thermodynamic conditions, the hydrogen bond network (HBN) acquires a diverse spectrum of topologies distinctly different from that of liquid water and of ice VII simulated at the same pressure. We observe that the HBN topology of plastic ice carries some degree of similarity with the crystal phase, stronger at thermodynamic conditions proximal to ice VII, and gradually lessening upon approaching the liquid state. Our results enrich our understanding of the properties of water at high pressure and high temperature, and may help in rationalizing the geology of

Structure and phase behavior of high-density ice from molecular-dynamics simulations with the ReaxFF potential

We report a molecular dynamics simulation study of dense ice modeled by the reactive force field (ReaxFF) potential, focusing on the possibility of phase changes between crystalline and plastic phases as observed in earlier simulation studies with rigid water models. It is demonstrated that the present model system exhibits phase transitions, or crossovers, among ice VII and two plastic ices with face-centered cubic (fcc) and body-centered cubic (bcc) lattice structures. The phase diagram derived from the ReaxFF potential is different from those of the rigid water models in that the bcc plastic phase lies on the high-pressure side of ice VII and does the fcc plastic phase on the low-pressure side of ice VII. The phase boundary between the fcc and bcc plastic phases on the pressure, temperature plane extends to the high-temperature region from the triple point of ice VII, fcc plastic, and bcc plastic phases. Proton hopping, i.e., delocalization of a proton, along between two neighboring oxygen atoms in dense ice is observed for the ReaxFF potential but only at pressures and temperatures both much higher than those at which ice VII-plastic ice transitions are observed.

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.

Computational Study on Homogeneous Melting of Benzene Phase I

Molecular-dynamics simulations are used for examining the microscopic details of the homogeneous melting of benzene phase I. The equilibrium melting temperatures of our model were initially determined using the direct-coexistence method. Homogeneous melting at a higher temperature is achieved by heating a defect- and surfacefree crystal. The temperature-dependent potential energy and lattice parameters do not indicate a premelting phase even under superheated conditions. Further, statistical analyses using induction times computed from 200 melting trajectories were conducted, denoting that the homogeneous melting of benzene occurs stochastically, and that there is no intermediate transient state between the crystal and liquid phases. Additionally, the critical nucleus size is estimated using the seeding approach, along with the local bond order parameter. We found that the large diffusive motion arising from defect migration or neighbor-molecule swapping is of little importance during nucleation. Instead, the orientational disorder activated using the flipping motion of the benzene plane results in the melting nucleus.

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.

Phase Diagram of TIP4P/2005 Water at High Pressure

We report a new ice phase that forms spontaneously at the interface between ice VII and liquid water in molecular dynamics simulations of TIP4P/2005 water. The new phase is structurally quite similar to an ice phase originally found to be a precursor in the course of the homogeneous nucleation of ice VII in liquid water. A part of the water molecules in these ice phases can rotate easily because the number of hydrogen bonds in them is less than four, and thus they can be regarded as partial plastic phases. A rough estimate suggests that these phases are thermodynamically more stable than either ice VI or ice VII for 3 GPa < P < 18 GPa at T = 300 K. Although the partial plastic phases would be metastable states at any point in the phase diagram of real water, they might be realized experimentally with the aid of dopants and/or solid surfaces.

Liquid-liquid phase separation of N-isopropylpropionamide aqueous solutions above the lower critical solution temperature

We investigate driving forces of the liquid–liquid phase separation of N-isopropylpropionamide (NiPPA) aqueous solutions above the lower critical solution temperature using molecular dynamics simulations. Spontaneous phase separations of the model aqueous solution with a modified OPLS-AA force field are observed above the experimentally determined cloud point. The destabilization toward the phase separation is confirmed by temperature dependence of the long-wavelength limit of the concentration-concentration structure factor, the dominant component of which is found to be an increasing effective attraction between NiPPA molecules. At varying temperatures, the potentials of mean force (PMFs) between a pair of NiPPA molecules at infinite dilution are obtained and decomposed into the nonpolar and Coulombic contributions. The nonpolar contribution, arising essentially from molecular volume, promotes association of NiPPA molecules with increasing temperature while the Coulombic one antagonizes the association. Thus, our analysis leads to a conclusion that the driving force of thermally induced aggregation of NiPPA molecules is the temperature dependence of the nonpolar contribution in PMF between NiPPA molecules, not the temperature dependence of the number or strength of hydrogen bonds between NiPPA and water molecules.

Solid-liquid critical behavior of a cylindrically confined Lennard-Jones fluid

Extensive molecular dynamics simulations have been performed to study the phase behavior of Lennard-Jones particles confined in a quasi-one-dimensional hydrophobic nanopore. We provide unambiguous evidence for a solid-liquid critical point by investigating (i) isotherms in the pressure-volume plane, (ii) the spontaneous solid-liquid phase separation below a certain temperature, (iii) diverging heat capacity and isothermal compressibility as a certain point is approached, (iv) continuous change of dynamical and structural properties above the point, (v) the finite-size scaling analysis of the density distribution below and above the point. The result combined with earlier studies of confined water suggests that the solid-liquid critical point is not uncommon in quasi-one- and quasi-two-dimensional fluids.

The atomistic mechanism of hcp-to-bcc martensitic transformation in the Ti–Nb system revealed by molecular dynamics simulations

The atomic mechanism of the hcp-to-bcc martensitic transformation and the corresponding crystallographic correlation has been elucidated.