Scaling Law for the Onset of Solidification at Extreme Undercooling.

Quasi-isentropic compression enables one to study the solidification of metastable liquid states that are inaccessible through other experimental means. The onset of this nonequilibrium solidification is known to depend on the compression rate and material-specific factors, but this complex interdependence has not been well characterized. In this study, we use a combination of experiments, theory, and computational simulations to derive a general scaling law that quantifies this dependence. One of its applications is a novel means to elucidate melt temperatures at high pressures.

Model for the Solid-Liquid Interfacial Free Energy at High Pressures.

The free energy involved in the formation of an interface between two phases (e.g., a solid-liquid interface) is referred to as the interfacial free energy. For the case of solidification, the interfacial free energy dictates the height of the energy barrier required to nucleate stable clusters of the newly forming solid phase and is essential for producing an accurate solidification kinetics model using classical nucleation theory (CNT)-based methods. While various methods have been proposed for modeling the interfacial free energy for solid-liquid interfaces in prior literature, many of these formulations involve making restrictive assumptions or approximations, such as the system being at or near equilibrium (i.e., the system temperature is approximately equal to the melt temperature) or that the system is at pressures close to atmospheric. However, these approximations and assumptions may break down in highly non-equilibrium situations, such as in dynamic-compression experiments where metastable liquids that are undercooled by hundreds of kelvin or overpressurized by several gigapascals or more are formed before eventually solidifying. We derive a solid-liquid interfacial free-energy model for such high-pressure conditions by considering the enthalpies of interactions between pairs of atoms or molecules. We also consider the contribution of interface roughness (disordering) by incorporating a multilayer interface model known as the Temkin n-layer model. Our formulation is applicable to a diverse variety of materials, and we demonstrate it by developing models specifically for two different materials: water and gallium. We apply our interfacial free-energy formulation to CNT-based kinetics simulations of several suites of dynamic-compression experiments that cause liquid water to solidify to the high-pressure solid polymorph ice VII and have found good agreement to the observed kinetics with only minor empirical fitting.

Measuring the melting curve of iron at super-Earth core conditions.

The discovery of more than 4500 extrasolar planets has created a need for modeling their interior structure and dynamics. Given the prominence of iron in planetary interiors, we require accurate and precise physical properties at extreme pressure and temperature. A first-order property of iron is its melting point, which is still debated for the conditions of Earth's interior. We used high-energy lasers at the National Ignition Facility and in situ x-ray diffraction to determine the melting point of iron up to 1000 gigapascals, three times the pressure of Earth's inner core. We used this melting curve to determine the length of dynamo action during core solidification to the hexagonal close-packed (hcp) structure. We find that terrestrial exoplanets with four to six times Earth's mass have the longest dynamos, which provide important shielding against cosmic radiation.

Numerical modeling of solid-cluster evolution applied to the nanosecond solidification of water near the metastable limit.

Classical nucleation theory (CNT) is a promising way to predictively model the submicrosecond kinetics of phase transitions that occur under dynamic compression, such as the suite of experiments performed over the past two decades on the solidification of liquid water to the high-pressure ice VII phase. Myint et al. [Phys. Rev. Lett. 121, 155701 (2018)] presented the first CNT-based model for these types of rapid phase transitions, but relied on an empirical scaling parameter in their transient induction model to simulate the lag time that occurs prior to the onset of significant formation of ice VII clusters in the system. To build on that study, we model the liquid water-ice VII phase transformation using a numerical discretization scheme to solve the Zel'dovich-Frenkel partial differential equation, which is a fundamental CNT-based kinetic equation that describes the statistical time-dependent behavior of solid cluster formation. The Zel'dovich-Frenkel equation inherently accounts for transience in the nucleation kinetics and eliminates the need for the empirical scaling factor used by Myint et al. One major result of this research is that transience is found to play a relatively small role in the nucleation process for the dynamic-compression time scales of the liquid water-ice VII experiments being simulated. Instead, we show that it is possible to properly model the lag time using steady-state CNT by making small refinements to the interfacial free energy value. We have also developed a new dimensionless parameter that may be applied a priori to predict whether or not transient nucleation will be important in a given dynamic-compression experiment.