Search

X-ray Thomson scattering (XRTS) has emerged as a powerful tool for the diagnostics of matter under extreme conditions. In principle, it gives one access to important system parameters such as the temperature, density, and ionization state, but the interpretation of the measured XRTS intensity usually relies on theoretical models and approximations. In this work, we show that it is possible to extract the Rayleigh weight -- a key property that describes the electronic localization around the ions -- directly from the experimental data without the need for any model calculations or simulations. As a practical application, we consider an experimental measurement of strongly compressed Be at the National Ignition Facility (NIF) [D\"oppner \emph{et al.}, \textit{Nature} \textbf{618}, 270-275 (2023)]. In addition to being interesting in their own right, our results will open up new avenues for diagnostics from \emph{ab initio} simulations, help to further constrain existing chemical models, and constitute a rigorous benchmark for theory and simulations.

Accurate identification of ice phases is essential for understanding various physicochemical phenomena. However, such classification for structures simulated with molecular dynamics is complicated by the complex symmetries of ice polymorphs and thermal fluctuations. For this purpose, both traditional order parameters and data-driven machine learning approaches have been employed, but they often rely on expert intuition, specific geometric information, or large training datasets. In this work, we present an unsupervised phase classification framework that combines a score-based denoiser model with a subsequent model-free classification method to accurately identify ice phases. The denoiser model is trained on perturbed synthetic data of ideal reference structures, eliminating the need for large datasets and labeling efforts. The classification step utilizes the Smooth Overlap of Atomic Positions (SOAP) descriptors as the atomic fingerprint, ensuring Euclidean symmetries and transferability to various structural systems. Our approach achieves a remarkable 100\% accuracy in distinguishing ice phases of test trajectories using only seven ideal reference structures of ice phases as model inputs. This demonstrates the generalizability of the score-based denoiser model in facilitating phase identification for complex molecular systems. The proposed classification strategy can be broadly applied to investigate structural evolution and phase identification for a wide range of materials, offering new insights into the fundamental understanding of water and other complex systems., Comment: 20 pages, 3 figures

An accurate description of information is relevant for a range of problems in atomistic machine learning (ML), such as crafting training sets, performing uncertainty quantification (UQ), or extracting physical insights from large datasets. However, atomistic ML often relies on unsupervised learning or model predictions to analyze information contents from simulation or training data. Here, we introduce a theoretical framework that provides a rigorous, model-free tool to quantify information contents in atomistic simulations. We demonstrate that the information entropy of a distribution of atom-centered environments explains known heuristics in ML potential developments, from training set sizes to dataset optimality. Using this tool, we propose a model-free UQ method that reliably predicts epistemic uncertainty and detects out-of-distribution samples, including rare events in systems such as nucleation. This method provides a general tool for data-driven atomistic modeling and combines efforts in ML, simulations, and physical explainability., Comment: v2.0

We combine nanosecond laser shock compression with \emph{in-situ} picosecond X-ray diffraction to provide structural data on iron up to 275 GPa. We constrain the extent of hcp-liquid coexistence, the onset of total melt, and the structure within the liquid phase. Our results indicate that iron, under shock compression, melts completely by 258(8) GPa. A coordination number analysis indicates that iron is a simple liquid at these pressure-temperature conditions. We also perform texture analysis between the ambient body-centered-cubic (bcc) $\alpha$, and the hexagonal-closed-packed (hcp) high-pressure $\epsilon-$phase. We rule out the Rong-Dunlop orientation relationship (OR) between the $\alpha$ and $\epsilon-$phases. However, we cannot distinguish between three other closely related ORs: Burger's, Mao-Bassett-Takahashi, and Potter's OR. The solid-liquid coexistence region is constrained from a melt onset pressure of 225(3) GPa from previously published sound speed measurements and full melt (246.5(1.8)-258(8) GPa) from X-ray diffraction measurements, with an associated maximum latent heat of melting of 623 J/g. This value is lower than recently reported theoretical estimates and suggests that the contribution to the earth's geodynamo energy budget from heat release due to freezing of the inner core is smaller than previously thought. Melt pressures for these nanosecond shock experiments are consistent with gas gun shock experiments that last for microseconds, indicating that the melt transition occurs rapidly.

Hydrocarbon mixtures are extremely abundant in the Universe, and diamond formation from them can play a crucial role in shaping the interior structure and evolution of planets. With first-principles accuracy, we first estimate the diamond nucleation rate in pure liquid carbon, and then reveal the nature of chemical bonding in hydrocarbons at extreme conditions. We finally establish the pressure-temperature phase boundary where diamond can form from hydrocarbon mixtures with different atomic fractions of carbon. Notably, we find a depletion zone at pressures above 200~GPa and temperatures below 3000~K-3500~K where diamond formation is thermodynamically favorable regardless of the carbon atomic fraction, due to a phase separation mechanism. The cooler condition of the interior of Neptune compared to Uranus means that the former is much more likely to contain the depletion zone. Our findings can help explain the dichotomy of the two ice giants manifested by the low luminosity of Uranus, and lead to a better understanding of (exo-)planetary formation and evolution.

Most experimentally known high-pressure ice phases have a body-centred cubic (bcc) oxygen lattice. Our atomistic simulations show that, amongst these bcc ice phases, ices VII, VII' and X are the same thermodynamic phase under different conditions, whereas superionic ice VII'' has a first-order phase boundary with ice VII'. Moreover, at about 300 GPa, ice X transforms into the Pbcm phase with a sharp structural change but no apparent activation barrier, whilst at higher pressures the barrier gradually increases. Our study thus clarifies the phase behaviour of the high-pressure insulating ices and reveals peculiar solid-solid transition mechanisms not known in other systems.

We derive a continuum-level plasticity model for polycrystalline materials in the high energy density regime, based on a single dislocation density and single mobility mechanism, with an evolution model for the dislocation density. The model is formulated explicitly in terms of quantities connected closely with equation of state (EOS) theory, in particular the shear modulus and Einstein temperature, which reduces the number of unconstrained parameters while increasing the range of applicability. The least constrained component is the Peierls barrier $E_P$, which is however accessible by atomistic simulations. We demonstrate an efficient method to estimate the variation of $E_P$ with compression, constrained to fit a single flow stress datum. The formulation for dislocation mobility accounts for some or possibly all of the stiffening at high strain rates usually attributed to phonon drag. The configurational energy of the dislocations is accounted for explicitly, giving a self-consistent calculation of the conversion of plastic work to heat. The configurational energy is predicted to contribute to the mean pressure, and may reach several percent in the terapascal range, which may be significant when inferring scalar EOS data from dynamic loading experiments. The bulk elastic strain energy also contributes to the pressure, but appears to be much smaller. Although inherently describing the plastic relaxation of elastic strain, the model can be manipulated to estimate the flow stress as a function of mass density, temperature, and strain rate, which is convenient to compare with other models and inferences from experiment. The deduced flow stress reproduces systematic trends observed in elastic waves and instability growth experiments, and makes testable predictions of trends versus material and crystal type over a wide range of pressure and strain rate., Comment: Various corrections; updates and additions to references; response to referee comments

Accurately modeling dense plasmas over wide ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/cm$^3$ to 16 g/cm$^3$ and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature (SQ) method. We find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the SESAME database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.

Atom-in-jellium calculations of the Einstein frequency in condensed matter and of the equation of state were used to predict the variation of shear modulus from zero pressure to ~$10^7$ g/cm$^3$, for several elements relevant to white dwarf (WD) stars and other self-gravitating systems. This is by far the widest range reported electronic structure calculation of shear modulus, spanning from ambient through the one-component plasma to extreme relativistic conditions. The predictions were based on a relationship between Debye temperature and shear modulus which we assess to be accurate at the o(10%) level, and is the first known use of atom-in-jellium theory to calculate a shear modulus. We assessed the overall accuracy of the method by comparing with experimental measurements and more detailed electronic structure calculations at lower pressures.

Most water in the universe may be superionic, and its thermodynamic and transport properties are crucial for planetary science but difficult to probe experimentally or theoretically. We use machine learning and free energy methods to overcome the limitations of quantum mechanical simulations, and characterize hydrogen diffusion, superionic transitions, and phase behaviors of water at extreme conditions. We predict that a close-packed superionic phase with mixed stacking is stable over a wide temperature and pressure range, while a body-centered cubic phase is only thermodynamically stable in a small window but is kinetically favored. Our phase boundaries, which are consistent with the existing-albeit scarce-experimental observations, help resolve the fractions of insulating ice, different superionic phases, and liquid water inside of ice giants.

Atom-in-jellium calculations of the electron states, and perturbative calculations of the Einstein frequency, were used to construct equations of state (EOS) from around $10^{-5}$ to $10^7$g/cm$^3$ and $10^{-4}$ to $10^{6}$eV for elements relevant to white dwarf (WD) stars. This is the widest range reported for self-consistent electronic shell structure calculations. Elements of the same ratio of atomic weight to atomic number were predicted to asymptote to the same $T=0$ isotherm, suggesting that, contrary to recent studies of the crystallization of WDs, the amount of gravitational energy that could be released by separation of oxygen and carbon is small. A generalized Lindemann criterion based on the amplitude of the ion-thermal oscillations calculated using atom-in-jellium theory, previously used to extrapolate melt curves for metals, was found to reproduce previous thermodynamic studies of the melt curve of the one component plasma with a choice of vibration amplitude consistent with low pressure results. For elements for which low pressure melting satisfies the same amplitude criterion, such as Al, this melt model thus gives a likely estimate of the melt curve over the full range of normal electronic matter; for the other elements, it provides a useful constraint on the melt locus.

We present an accurate and efficient real-space formulation of the Hellmann-Feynman stress tensor for $\mathcal{O}(N)$ Kohn-Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi-Dirac distribution becomes smoother and density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to $\mathcal{O}(N)$ methods. We then describe its evaluation within the $\mathcal{O}(N)$ infinite-cell Clenshaw-Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic as well as insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results, and show convergence with respect to mesh size to established planewave results. We employ the new formulation to compute the viscosity of hydrogen at a million kelvin from Kohn-Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods., Comment: 13 pages, 7 figures

Diamond is used extensively as a component in high energy density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first principles theoretical EOS models, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Greuneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. The resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.

Ru and Rh are interesting cases for comparing equations of state (EOS), because most general purpose EOS are semi-empirical, relying heavily on shock data, and none has been reported for Ru. EOS were calculated for both elements using all-electron atom-in-jellium theory, and cold compression curves were calculated for the common crystal types using the multi-ion pseudopotential approach. Previous EOS constructed for these elements used Thomas-Fermi (TF) theory for the electronic behavior at high temperatures, which neglects electronic shell structure; the atom-in-jellium EOS exhibited pronounced features from the excitation of successive electron shells. Otherwise, the EOS matched surprisingly well, especially considering the lack of experimental data for Ru. The TF-based EOS for Ru may however be inaccurate in the multi-terapascal range needed for some high energy density experiments. The multi-ion calculations predicted that the hexagonal close-packed phase of Ru remains stable to at least 2.5 TPa and possibly 10 TPa, and that its c/a should gradually increase to the ideal value. A method was devised to estimate the variation in Debye temperature from the cold curve, and thus estimate the ion-thermal EOS without requiring relatively expensive dynamical force calculations, in a form convenient for adjusting EOS or phase boundaries. The Debye temperature estimated in this way was similar to the result from atom-in-jellium calculations. We also predict the high-pressure melt loci of both elements.

Equations of state (EOS) calculated from a computationally efficient atom-in-jellium treatment of the electronic structure have recently been shown to be consistent with more rigorous path integral Monte Carlo (PIMC) and quantum molecular dynamics (QMD) simulations of metals in the warm dense matter regime. Here we apply the atom-in-jellium model to predict wide-ranging EOS for the cryogenic liquid elements nitrogen, oxygen, and fluorine. The principal Hugoniots for these substances were surprisingly consistent with available shock data and Thomas-Fermi (TF) EOS for very high pressures, and exhibited systematic variations from TF associated with shell ionization effects, in good agreement with PIMC, though deviating from QMD and experiment in the molecular regime. The new EOS are accurate much higher in pressure than previous widely-used models for nitrogen and oxygen in particular, and should allow much more accurate predictions for oxides and nitrides in the liquid, vapor, and plasma regime, where these have previously been constructed as mixtures containing the older EOS.

Atom-in-jellium calculations of the Einstein frequency were used to calculate the mean displacement of an ion over a wide range of compression and temperature. Expressed as a fraction of the Wigner-Seitz radius, the displacement is a measure of the asymptotic freedom of the ion at high temperature, and thus of the change in heat capacity from 6 to 3 quadratic degrees of freedom per atom. A functional form for free energy was proposed based on the Maxwell-Boltzmann distribution as a correction to the Debye free energy, with a single free parameter representing the effective density of potential modes to be saturated. This parameter was investigated using molecular dynamics simulations, and found to be ~0.2 per atom. In this way, the ion-thermal contribution can be calculated for a wide-range equation of state (EOS) without requiring a large number of molecular dynamics simulations. Example calculations were performed for carbon, including the sensitivity of key EOS loci to ionic freedom.

Recent path-integral Monte Carlo and quantum molecular dynamics simulations have shown that computationally efficient average-atom models can predict thermodynamic states in warm dense matter to within a few percent. One such atom-in-jellium model has typically been used to predict the electron-thermal behavior only, although it was previously developed to predict the entire equation of state (EOS). We report completely atom-in-jellium EOS calculations for Be, Al, Si, Fe, and Mo, as elements representative of a range of atomic number and low-pressure electronic structure. Comparing the more recent method of pseudo-atom molecular dynamics, atom-in-jellium results were similar: sometimes less accurate, sometimes more. All these techniques exhibited pronounced effects of electronic shell structure in the shock Hugoniot which are not captured by Thomas-Fermi based EOS. These results demonstrate the value of a hierarchical approach to EOS construction, using average-atom techniques with shell structure to populate a wide-range EOS surface efficiently, complemented by more rigorous 3D multi-atom calculations to validate and adjust the EOS.

The validity of the widely used linear mixing approximation for the equations of state (EOS) of planetary ices is investigated at pressure-temperature conditions typical for the interior of Uranus and Neptune. The basis of this study are ab initio data ranging up to 1000 GPa and 20 000 K calculated via density functional theory molecular dynamics simulations. In particular, we calculate a new EOS for methane and EOS data for the 1:1 binary mixtures of methane, ammonia, and water, as well as their 2:1:4 ternary mixture. Additionally, the self-diffusion coefficients in the ternary mixture are calculated along three different Uranus interior profiles and compared to the values of the pure compounds. We find that deviations of the linear mixing approximation from the results of the real mixture are generally small; for the thermal EOS they amount to 4% or less. The diffusion coefficients in the mixture agree with those of the pure compounds within 20% or better. Finally, a new adiabatic model of Uranus with an inner layer of almost pure ices is developed. The model is consistent with the gravity field data and results in a rather cold interior ($\mathrm{T_{core}} \mathtt{\sim}$ 4000 K).

Density Functional Theory calculations using the quasi-harmonic approximation have been used to calculate the solid Hugoniot of two polytypes of boron carbide up to 100 GPa. Under the assumption that segregation into the elemental phases occurs around the pressure that the B11Cp(CBC) polytype becomes thermodynamically unstable with respect to boron and carbon, two discontinuities in the Hugoniot, one at 50 GPa and one at 90-100 GPa, are predicted. The former is a result of phase segregation, and the latter a phase transition within boron. First principles molecular dynamics (FPMD) simulations were employed to calculate the liquid Hugoniot of B4C up to 1.5 TPa, and the results are compared to recent experiments carried out at the Omega Laser Facility up to 700 GPa [Phys. Rev. B 94, 184107 (2016)]. A generally good agreement between theory and experiment was observed. Analysis of the FPMD simulations provides evidence for an amorphous, but covalently bound, fluid below 438 GPa, and an atomistic fluid at higher pressures and temperatures., Comment: 10 pages, 6 figures

We simulate bulk materials under strong currents by following in real-time the dynamics of the electrons under an electric field. By changing the intensity of the electric field, our method can model, for the first time, non-linear effects in the conductivity from first principles. To illustrate our approach, we show calculations that predict that liquid aluminum exhibits negative-differential conductivity for current densities of the order of $10^{12}-10^{13}~\mathrm{A/cm^2}$. We find that the change in the non-linear conductivity emerges from a competition between the accumulation of charge around the nuclei that increases the scattering of the conduction electrons, and a decreasing scattering cross-section at high currents., Comment: 8 pages, 8 figures

Megabar (1 Mbar = 100 GPa) laser shocks on precompressed samples allow reaching unprecedented high densities and moderately high 10000-100000K temperatures. We describe here a complete analysis framework for the velocimetry (VISAR) and pyrometry (SOP) data produced in these experiments. Since the precompression increases the initial density of both the sample of interest and the quartz reference for pressure-density, reflectivity and temperature measurements, we describe analytical corrections based on available experimental data on warm dense silica and density-functional-theory based molecular dynamics computer simulations. Using our improved analysis framework we report a re-analysis of previously published data on warm dense hydrogen and helium, compare the newly inferred pressure, density and temperature data with most advanced equation of state models and provide updated reflectivity values., Comment: 13 pages, 16 figures, 4 tables

We present a 5-phase equation of state (EOS) for elemental carbon. The phases considered are: diamond, BC8, simple-cubic, simple-hexagonal, and the liquid/plasma state. The solid phase free energies are constrained by density functional theory (DFT) calculations. Vibrational contributions to the free energy of each solid phase are treated within the quasiharmonic framework. The liquid free energy model is constrained by fitting to a combination of DFT molecular dynamics performed over the range 10,000 K < T < 100,000 K, and path integral quantum Monte Carlo calculations for T > 100,000 K (both for {\rho} between 3 and 12 g/cc, with select higher-{\rho} DFT calculations as well). The liquid free energy model includes an atom-in-jellium approach to account for the effects of ionization due to temperature and pressure in the plasma state, and an ion-thermal model which includes the approach to the ideal gas limit. The precise manner in which the ideal gas limit is reached is greatly constrained by both the highest temperature DFT data and the path integral data, forcing us to discard an ion-thermal model we had used previously in favor of a new one. Predictions are made for the principal Hugoniot and the room-temperature isotherm, and comparisons are made to recent experimental results., Comment: 41 pages, 14 figures, Supplementary Material

For planets other than Earth, interpretation of the composition and structure depends largely on comparing the mass and radius with the composition expected given their distance from the parent star. The composition implies a mass-radius relation which relies heavily on equations of state calculated from electronic structure theory and measured experimentally on Earth. We lay out a method for deriving and testing equations of state, and deduce mass-radius and mass-pressure relations for key materials whose equation of state is reasonably well established, and for differentiated Fe/rock. We find that variations in the equation of state, such as may arise when extrapolating from low pressure data, can have significant effects on predicted mass- radius relations, and on planetary pressure profiles. The relations are compared with the observed masses and radii of planets and exoplanets. Kepler-10b is apparently 'Earth- like,' likely with a proportionately larger core than Earth's, nominally 2/3 of the mass of the planet. CoRoT-7b is consistent with a rocky mantle over an Fe-based core which is likely to be proportionately smaller than Earth's. GJ 1214b lies between the mass-radius curves for H2O and CH4, suggesting an 'icy' composition with a relatively large core or a relatively large proportion of H2O. CoRoT-2b is less dense than the hydrogen relation, which could be explained by an anomalously high degree of heating or by higher than assumed atmospheric opacity. HAT-P-2b is slightly denser than the mass-radius relation for hydrogen, suggesting the presence of a significant amount of matter of higher atomic number. CoRoT-3b lies close to the hydrogen relation. The pressure at the center of Kepler-10b is 1.5+1.2-1.0 TPa. The central pressure in CoRoT-7b is probably close to 0.8TPa, though may be up to 2TPa., Comment: Added more recent exoplanets. Tidied text and references. Added extra "rock" compositions. Responded to referee comments