In this work, we present a general form for the electron-ion diffractive potential derived from the quantum pair density matrix and fit to the improved Kelbg potential for atomic numbers up to Z = 54. We apply classical molecular dynamics using the improved Kelbg potential for carbon with various forms of the Pauli potential to compute internal energies and pressures for hot, dense plasma conditions. Our results are compared to an equation of state model based on path integral Monte Carlo and density functional theory simulations to examine the extent to which the improved Kelbg potential reproduces the internal energy and pressure of carbon plasmas. The regions of validity for carbon agree generally with those derived previously for hydrogen once pressure ionization effects are incorporated. Based on our carbon results and previously published hydrogen studies, we discuss the general applicability and limitations of these potentials for equation of state studies in warm dense matter and high energy density plasmas.

We present anab initiofree energy model derived from a fully relativistic density functional theory (DFT) electronic structure with dynamic magnetism forδ-plutonium (face-centered cubic, fcc). The DFT model is extended with orbital-orbital interaction in a parameter free orbital polarization (OP) mechanism consistent with previous modeling of plutonium. Gibbs free energy is built from components associated with the temperature dependence of the electronic structure and the corresponding electronic entropy, lattice vibrations within an anharmonic lattice dynamics model, and dynamical fluctuations of the magnetization density, i.e. magnetic fluctuations. The fluctuation model consists of transverse and longitudinal modes driven by temperature induced excitations of the DFT + OP electronic structure. Theab initiomodel thus incorporates fluctuating states beyond the electronic ground state. Thanks to the dynamic magnetism, the theory predicts excellent thermodynamic properties and a Gibbs free energy in accord with CALPHAD and semi-empirical modeling developed from the thermodynamic observables. The magnetic fluctuations further explain anomalous behaviors of the thermal expansion in plutonium. Specifically, a thermal expansion for theδ-plutonium system turning from positive to negative at temperatures above room temperature, a tendency for gallium to reduce and remove the negative thermal expansion depending on composition, and a positive thermal expansion for the high temperatureϵphase.

Using pulsed-power magnetic field sources to compress gallium to gigapascal pressures on nanosecond timescales, we report here experiments on shockless dynamic compression of a liquid metal. Time-resolved velocimetry data reveal signatures of rapid freezing from a metastable liquid state, and we demonstrate that the kinetics of this nonequilibrium solidification can be accurately simulated with a computational modeling framework we have developed in previous studies, where classical nucleation theory is coupled with hydrodynamics. Notably, velocity traces in some of our experiments show evidence of a phase transition, while others do not, even though other types of evidence suggest that solidification may be occurring in all of them. We explain how predictions made by our models regarding the presence or absence of these phase-transition signatures motivated additional experiments that later confirmed the theoretical predictions. Our analysis shows that due to the rapid, quasi-isentropic nature of the loading path, our experiments were able to compress liquid gallium to metastable states that are undercooled below the equilibrium melt temperature by more than 300 K and exhibit pressures that approach five times the equilibrium melt pressure. The understanding gained in this study should form the basis for future dynamic-compression experiments aimed at interrogating melt curves at high pressures.

We present a multiscale simulation framework that couples the Finite Element Method with molecular dynamics. Bypassing traditional equations of state (EOS) by using in-line atomistic simulations, the method offers the advantage of incorporating detailed microscale physics not easily represented with coarse-grained models. Coupling consistency with the continuum code is ensured through the use of lifting and restriction operators, in line with heterogeneous multiscale methods. The concurrent continuum-atomistic framework is validated through comparison with experimental results and conventional EOS models, and demonstrated in a shock-driven hydrodynamic flow simulation under extreme conditions. We further evaluate the framework's usability by comparing it to state-of-the-art EOS models of deuterium. A computational performance study reveals that the atomistic EOS evaluation is a feasible alternative to conventional approaches, and demonstrates a weak scaling of 99

We present a multiscale simulation framework that couples the finite-element method with molecular dynamics. Bypassing traditional equations of state (EOS) by using in-line atomistic simulations, the method offers the advantage of incorporating detailed microscale physics not easily represented with coarse-grained models. Coupling consistency with the continuum code is ensured through the use of lifting and restriction operators, in line with heterogeneous multiscale methods. The concurrent continuum-atomistic framework is validated through comparison with experimental results and conventional EOS models, and demonstrated in a shock-driven hydrodynamic flow simulation under extreme conditions. We further evaluate the framework's usability by comparing it to state-of-the-art EOS models of deuterium. A computational performance study reveals that the atomistic EOS evaluation is a feasible alternative to conventional approaches, and demonstrates a weak scaling of 99% efficiency. These results highlight the framework's potential for large-scale multiscale modeling across a broad range of materials and conditions.

We have developed a first-passage kinetic Monte Carlo approach for materials aging to investigate the sensitivity of void swelling to model parameters, including helium bubble density and size distribution. In addition to explicitly accounting for the spatial distribution of individual point defects, bubbles, and voids, our approach can simulate total doses equivalent to 100 years of natural aging on statistically representative volumes of materials. This technique enables us to study the effects on swelling and radiation damage evolution due to temperature and dose rate (as altered in artificially aged experiments), differences in effective interaction radii between vacancies and interstitials, and varying defect diffusion activation energies, while providing more detailed information than previous rate-equation based approaches. Our results indicate that spatial effects that are not modeled in mean-field rate theories could play a significant role in void swelling initiation and growth for certain regimes of model parameters.

On December 5, 2022, an indirect drive fusion implosion on the National Ignition Facility (NIF) achieved a target gain G_{target} of 1.5. This is the first laboratory demonstration of exceeding "scientific breakeven" (or G_{target}>1) where 2.05 MJ of 351 nm laser light produced 3.1 MJ of total fusion yield, a result which significantly exceeds the Lawson criterion for fusion ignition as reported in a previous NIF implosion [H. Abu-Shawareb et al. (Indirect Drive ICF Collaboration), Phys. Rev. Lett. 129, 075001 (2022)PRLTAO0031-900710.1103/PhysRevLett.129.075001]. This achievement is the culmination of more than five decades of research and gives proof that laboratory fusion, based on fundamental physics principles, is possible. This Letter reports on the target, laser, design, and experimental advancements that led to this result.

We review and discuss Lawrence Livermore's multiphase, multicomponent mixture EOS generation code, MM Mix, highlighting its evolving capabilities and improved computational efficiency. The new capabilities include the possibility to treat line compounds (i.e., multicomponent systems with fixed stoichiometry), and a streamlined and practical method for constructing distinct ion-thermal and electron-thermal contributions to the total mixture EOS needed for radiationhydrodynamic simulations (of inertial-confinement fusion, for example), which require separate ion and electron temperatures. Dramatic speedups in the running of MM Mix are reported that benefit from parallelization on multicore processors and conversion from Python to modern C++, and preliminary predictions of the composition-dependence of the high-pressure melt temperature in iron nickel (Fe Ni) mixtures are presented as an illustrative example.

The classical Richtmyer-Meshkov instability is a hydrodynamic instability characterizing the evolution of an interface following shock loading. In contrast to other hydrodynamic instabilities such as Rayleigh-Taylor, it is known for being unconditionally unstable: regardless of the direction of shock passage, any deviations from a flat interface will be amplified. In this article, we show that for negative Atwood numbers, there exist special sequences of shocks which result in a nearly perfectly suppressed instability growth. We demonstrate this principle computationally and experimentally with stepped fliers and phase transition materials. A fascinating immediate corollary is that in specific instances a phase transitioning material may self-suppress RMI.

We report measurements of the compressibility of ramp compressed tantalum to a final stress of 2.3 TPa corresponding to threefold volumetric compression. Using these data, we extended the experimental constraint on the Ta cold compression curve by an order of magnitude in pressure. By combining the resulting data with previous measurements of shock compression and ambient pressure heating, we construct an experimentally bounded and thermodynamically consistent equation of state model for Ta which has 2% uncertainty in pressure at 1 TPa. We therefore propose Ta as an in situ pressure scale for laser-heated static compression experiments which were recently able to reach terapascal pressures and thousands of degrees Kelvin. Our new equation of state of Ta is experimentally constrained at extreme pressures and temperatures relevant to a wide range of planetary interiors and will allow for more accurate comparison between experimental measurements and theory at extreme conditions.

For the Ta equation of state (EOS) model introduced in our recent publication (Gorman et al., 2023), we have made a fairly accurate analytic fit to this EOS over a limited range. The detailed functional form and associated parameters of this fit are presented here.

We describe the construction of a wide-ranged multiphase equation of state for elemental iron, consisting of four solid phases and one fluid phase. The free-energy models for the phases are constrained by fitting to a broad swath of ambient, static, and dynamic high-pressure experimental data, including recent ramp-compression and shock-ramp recrystallization experiments. In order to better describe the conditions near Earth's inner-core boundary, ab initio electronic-structure calculations of various types have been carried out for the hcp and liquid phases and used as additional constraints. Given the variability of reported iron melting temperatures in the 1-3 Mbar range, we construct EOS model variations based on two separate Tmelt(P) curves: lower-Tmelt and higher-Tmelt. Our analysis shows that only the higher-Tmelt(P) data are consistent with reported shock-melting pressure measurements. Furthermore, our examination of the recrystallization data from a recent shock-ramp x-ray diffraction study motivates a family of EOS models that address experimental uncertainties and provide a range of possible Tmelt(P) curves that have a somewhat lower melt temperature above 5 Mbar than the one reported in that study. This family of multiphase EOS models allow us to suggest new experiments to further reduce the uncertainties on the high-pressure melt curve, and they should also enable more accurate predictions of complex high-P, T processes involving iron. We show their potential utility by applying these EOSs to the investigation of properties of interest to planetary science, such as the vaporization entropy, critical temperature and density, Earth inner-core boundary temperature plus inner- and outer-core density deficits.

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.

We report on measurements of the ion-electron energy-transfer cross section utilizing low-velocity ion stopping in high-energy-density plasmas at the OMEGA laser facility. These measurements utilize a technique that leverages the close relationship between low-velocity ion stopping and ion-electron equilibration. Shock-driven implosions of capsules filled with D^{3}He gas doped with a trace amount of argon are used to generate densities and temperatures in ranges from 1×10^{23} to 2×10^{24} cm^{-3} and from 1.4 to 2.5 keV, respectively. The energy loss of 1-MeV DD tritons and 3.7-MeV D^{3}He alphas that have velocities lower than the average velocity of the thermal electrons is measured. The energy loss of these ions is used to determine the ion-electron energy-transfer cross section, which is found to be in excellent agreement with quantum-mechanical calculations in the first Born approximation. This result provides an experimental constraint on ion-electron energy transfer in high-energy-density plasmas, which impacts the modeling of alpha heating in inertial confinement fusion implosions, magnetic-field advection in stellar atmospheres, and energy balance in supernova shocks.

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.

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 models of equilibrium multiphase EOS, 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 -10 TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Gruneisen 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.

with previously published lower strain data (101–107 s-1), we observe a change in rate dependence, suggesting a transition from thermally activated to defect-limited (phonon drag) dislocation motion occurring at a strain rate of about 105 s-1.

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 g/cm, 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.

all co-existing phases combined with deviatoric stress averaging through a volume fraction weighted ow stress and a volume fraction weighted shear modulus. Implementations of the framework based on equilibrium phase fractions (i.e. instantaneous kinetic rate) and for finite rate transition kinetics are proposed.

For more than half a century, researchers around the world have been engaged in attempts to achieve fusion ignition as a proof of principle of various fusion concepts. As recently reported, a burning plasma state, where the alpha-heating in the plasma is the primary source of heating, was achieved in laboratory experiments. Following the Lawson criterion, an ignited plasma is one where the fusion heating power is high enough to overcome all the physical processes that cool the fusion plasma, creating a positive thermodynamic feedback loop with rapidly increasing temperature. In inertially confined fusion, ignition is a state where the fusion plasma can begin ``burn propagation'' into surrounding cold fuel, enabling the possibility of high energy gain. While ``scientific breakeven'' (i.e. unity target gain) has not yet been achieved, this talk reports the first controlled fusion experiment on the National Ignition Facility to produce capsule gain greater than unity (here 5.8) and reach ignition by many different formulations of the Lawson criterion. In the talk, we will discuss some key basic physics inertial confinement fusion (ICF) principles behind the burning plasma and ignition results as well as discuss future challenges.