In dynamic-compression experiments, the line-imaging Velocity Interferometer System for Any Reflector (VISAR) is a well-established diagnostic used to probe the velocity history, including wave profiles derived from dynamically compressed interfaces and wavefronts, depending on material optical properties. Knowledge of the velocity history allows for the determination of the pressure achieved during compression. Such a VISAR analysis is often based on Fourier transform techniques and assumes that the recorded interferograms are free from image distortions. In this paper, we describe the VISAR diagnostic installed at the HED-HIBEF instrument located at the European XFEL along with its calibration and characterization. It comprises a two-color (532, 1064 nm), three-arm (with three velocity sensitivities) line imaging system. We provide a procedure to correct VISAR images for geometric distortions and evaluate the performance of the system using Fourier analysis. We finally discuss the spatial and temporal calibrations of the diagnostic. As an example, we compare the pressure extracted from the VISAR analysis of shock-compressed polyimide and silicon.

Iron is a major component of the cores of the Earth and inhabited exoplanets. Its phase diagram at extreme pressures (P) and temperatures (T) is the subject of extensive debate. While recent experiments provide evidence for the stability of the body-centered cubic (bcc) phase, several theoretical studies point to the stability (even though marginal) of the hexagonal close-packed phase (hcp). None of those studies considered the itinerant magnetism of iron at extreme conditions. We compute the high-pressure phase diagram of Fe using density functional theory-based molecular dynamics (DFT MD) in which the paramagnetic nature of Fe is treated within the model of thermally induced longitudinal spin fluctuations (LSF). The LSF DFT MD with 16 valence electrons favors bcc phase stability. Two-phase large-scale simulations with quantum accurate machine learning potentials provide us with both melting and hcp-bcc phase boundaries. The computed phase diagram agrees with most of the experimental data and solves most of the numerous controversies. We conclude that the account for magnetism results in the new physics of iron under extreme conditions and brings the theory in agreement with experiment and seismic data. We expect that the approach we use can be applied to other metals where itinerant magnetism is important.

Diamond possesses exceptional physical properties due to its remarkably strong carbon-carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum-accurate multimillion atom molecular dynamics (MD) simulations, we have uncovered the extreme metastability of diamond at very high pressures, significantly exceeding its range of thermodynamic stability. We predict the post-diamond BC8 phase to be experimentally accessible only within a narrow high pressure-temperature region of the carbon phase diagram. The diamond to BC8 transformation proceeds through premelting followed by BC8 nucleation and growth in the metastable carbon liquid. We propose a double-shock compression pathway for BC8 synthesis, which is currently being explored in experiments at the National Ignition Facility.

Iron is the major element of the Earth's core and the cores of Earth-like exoplanets. The crystal structure of iron, the major component of the Earth's solid inner core (IC), is unknown under the high pressures (P) (3.3-3.6 Mbar) and temperatures (T) (5000-7000 K) and conditions of the IC and exoplanetary cores. Experimental and theoretical data on the phase diagram of iron at these extreme PT conditions are contradictory. Though some of the large-scale ab initio molecular dynamics (AIMD) simulations point to the stability of the body-centered cubic (bcc) phase, the latest experimental data are often interpreted as evidence for the stability of the hexagonal close-packed (hcp) phase. Applying large-scale AIMD, we computed the properties of iron phases at the experimental pressures and temperatures reported in the experimental papers. The use of large-scale AIMD is critical since the use of small bcc computational cells (less than approximately 1000 atoms) leads to the collapse of the bcc structure. Large-scale AIMD allowed us to compare the measured and computed coordination numbers as well as the measured and computed structural factors. This comparison, in turn, allowed us to suggest that the computed density, coordination number, and structural factors of the bcc phase are in agreement with those observed in experiments, which were previously assigned either to the liquid or hcp phase.

Knowledge of high-pressure melting curves of silicate minerals is critical for modeling the thermal-chemical evolution of rocky planets. However, the melting temperature of davemaoite, the third most abundant mineral in Earth's lower mantle, is still controversial. Here, we investigate the melting curves of two minerals, MgSiO3 bridgmanite and CaSiO3 davemaoite, under their stability field in the mantle by performing first-principles molecular dynamics simulations based on the density functional theory. The melting curve of bridgmanite is in excellent agreement with previous studies, confirming a general consensus on its melting temperature. However, we predict a much higher melting curve of davemaoite than almost all previous estimates. Melting temperature of davemaoite at the pressure of core-mantle boundary (similar to 136 gigapascals) is about 7700(150) K, which is approximately 2000 K higher than that of bridgmanite. The ultrarefractory nature of davemaoite is critical to reconsider many models in the deep planetary interior, for instance, solidification of early magma ocean and geodynamical behavior of mantle rocks.

The solid Earth's inner core (IC) is a sphere with a radius of about 1300 km in the center of the Earth. The information about the IC comes mainly from seismic studies. The composition of the IC is obtained by matching the seismic data and properties of candidate phases subjected to high pressure (P) and temperature (T). The close match between the density of the IC and iron suggests that the main constituent of the IC is iron. However, the stable phase of iron is still a subject of debate. One such iron phase, the body-centered cubic phase (bcc), is dynamically unstable at pressures of the IC (330-364 GPa) and low T but gets stabilized at high T characteristic of the IC (5000-7000 K). So far, ab initio molecular dynamics (AIMD) studies attempted to compute the bcc elastic properties for a small (order of 102) number of atoms. The mechanism of the bcc stabilization cannot be enabled in such cells and that has led to erroneous results. Here we apply AIMD to compute elastic moduli and sound velocities of the Fe bcc phase for a 2000 Fe atom computational cell, which is a cell of unprecedented size for ab initio calculations of iron. Unlike in previous ab initio calculations, both the longitudinal and the shear sound velocities of the Fe bcc phase closely match the properties of the IC material at P = 360 GPa and T = 6600 K, likely the PT conditions in the IC. The calculated density of the bcc iron at these PT conditions is just 3% higher than the density of the IC material according to the Preliminary Earth Model. This suggests that the widely assumed amount of light elements in the IC may need a reconsideration. The anisotropy of the bcc phase is an exact match to the most recent seismic studies. 

A spectral neighbor analysis (SNAP) machine learning interatomic potential (MLIP) has been developed for simulations of carbon at extreme pressures (up to 5 TPa) and temperatures (up to 20 000 K). This was achieved using a large database of experimentally relevant quantum molecular dynamics (QMD) data, training the SNAP potential using a robust machine learning methodology, and performing extensive validation against QMD and experimental data. The resultant carbon MLIP demonstrates unprecedented accuracy and transferability in predicting the carbon phase diagram, melting curves of crystalline phases, and the shock Hugoniot, all within 3% of QMD. By achieving quantum accuracy and efficient implementation on leadership-class high-performance computing systems, SNAP advances frontiers of classical MD simulations by enabling atomic-scale insights at experimental time and length scales.

Billion atom molecular dynamics (MD) using quantum-Accurate machine-learning Spectral Neighbor Analysis Potential (SNAP) observed long-sought high pressure BC8 phase of carbon at extreme pressure (12 Mbar) and temperature (5,000 K). 24-hour, 4650 node production simulation on OLCF Summit demonstrated an unprecedented scaling and unmatched real-world performance of SNAP MD while sampling 1 nanosecond of physical time. Efficient implementation of SNAP force kernel in LAMMPS using the Kokkos CUDA backend on NVIDIA GPUs combined with excellent strong scaling (better than 97% parallel efficiency) enabled a peak computing rate of 50.0 PFLOPs (24.9% of theoretical peak) for a 20 billion atom MD simulation on the full Summit machine (27,900 GPUs). The peak MD performance of 6.21 Matom-steps/node-s is 22.9 times greater than a previous record for quantum-Accurate MD. Near perfect weak scaling of SNAP MD highlights its excellent potential to advance the frontier of quantum-Accurate MD to trillion atom simulations on upcoming exascale platforms.

The crystal structure of iron, the major component of the Earth's inner core (IC), is unknown under the IC high pressure (P) (3.3-3.6 Mbar) and temperature (T) (5000-7000 K). Experimental and theoretical data on the phase diagram of iron at these extreme PT conditions are contradictory. Applying quasi-ab initio and ab initio molecular dynamics we computed free energies of the body-centered cubic (bcc), hexagonal close-packed (hcp), and liquid phases. The ionic free energies, computed for the embedded-atom model, were corrected for electronic entropy. Such correction brings the melting temperatures of the hcp iron in very good agreement with previous ab initio data. This validates the calculation of the bcc phase, where fully ab initio treatment is not technically possible due to large sizes required for convergence. The resulting phase diagram shows stabilization of the bcc phase prior to melting in the pressure range of the IC. The melting temperature of the bcc phase is equal to 7190 K at the pressure 360 GPa.

Using molecular dynamics simulation and evolutionary metadynamic calculations, a series of structures were revealed that possessed enthalpies and Gibbs energies lower than those of aragonite but higher than those of calcite. The structures are polytypes of calcite, differing in the stacking sequence of close packed (cp) Ca layers. The two- and six-layered polytypes have hexagonal symmetry P6(3)22 and were named hexarag and hexite, respectively. Hexarag is similar to aragonite, but with all the triangles placed on the middle distance between the cp layers. On the basis of the structures found, a two-step mechanism for the transformation of aragonite to calcite is suggested. In the first step, CO3 triangles migrate to halfway between the Ca layers with the formation of hexarag. In the second step, the two-layered cp (hcp) hexarag structure transforms into three-layered cp (fcc) calcite through a series of many-layered polytypes. The topotactic character of the transformation of aragonite to calcite, with [001] of aragonite being parallel to [0001] of calcite, is consistent with the suggested mechanism. High-temperature X-ray powder diffraction experiments did not reveal hexarag reflections. To assess the possibility of the formation of the polytypes found in nature or experiments, a TEM analysis of ground aragonite was performed. A grain was found that had six superstructure reflections in a direction perpendicular to the plane of the cp layer. This grain is believed to correspond to one of the predicted polytypes, with the diffuse character of the diffraction spots indicating a partial disordering of the cp layer stacking. A topological analysis was also performed, along with energy calculations, of the metastable high-pressure polymorphs CaCO3-II, -III, -Mb, and -VI. The similarity of CaCO3-II, -II, and -Mb to the calcite structure and the small energy difference explain the metastable formation of these polymorphs during the cold compression of calcite. On the basis of the performed analysis, the evolution of the CaCO3 cation array at calcite to a post-aragonite transformation is described.