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.

Thomson scattering spectra measure the response of plasma particles to incident radiation. In warm dense matter, which is opaque to visible light, x-ray Thomson scattering (XRTS) enables a detailed probe of the electron distribution and has been used as a diagnostic for electron temperature, density, and plasma ionization. In this work, we examine the sensitivities of inelastic XRTS signatures to modeling details including the dynamic collision frequency and the electronic density of states. Applying verified Monte Carlo inversion methods to dynamic structure factors obtained from time-dependent density functional theory, we assess the utility of XRTS signals as a way to inform the dynamic collision frequency, especially its direct-current (DC) limit, which is directly related to the electrical conductivity.

Magnetic reconnection is an important process in astrophysical environments, as it re-configures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the layer. In this paper, we perform 2D & 3D simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively-cooled regime. These simulations are performed in GORGON, an Eulerian resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and P$_{1/3}$ multi-group radiation transport. The load produces highly collisional, super-Alfv\'enic $(M_{A} \approx 1.5)$, supersonic $(M_S \approx 4-5)$ plasma flows which generate a reconnection layer ($L/{\delta} \approx 100, S_L \approx 400$). The reconnection layer undergoes radiative collapse when the radiative losses exceed Ohmic and compressional heating ${\tau}_{cool}^{-1}/{\tau}_A^{-1} \approx 100$; this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Finally, the current sheet is unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

Real-time time-dependent density functional theory (TDDFT) is presently the most accurate available method for computing electronic stopping powers from first principles. However, obtaining application-relevant results often involves either costly averages over multiple calculations or ad hoc selection of a representative ion trajectory. We consider a broadly applicable, quantitative metric for evaluating and optimizing trajectories in this context. This methodology enables rigorous analysis of the failure modes of various common trajectory choices in crystalline materials. Although randomly selecting trajectories is common practice in stopping power calculations in solids, we show that nearly 30% of random trajectories in an FCC aluminium crystal will not representatively sample the material over the time and length scales feasibly simulated with TDDFT, and unrepresentative choices incur errors of up to 60%. We also show that finite-size effects depend on ion trajectory via "ouroboros" effects beyond the prevailing plasmon-based interpretation, and we propose a cost-reducing scheme to obtain converged results even when expensive core-electron contributions preclude large supercells. This work helps to mitigate poorly controlled approximations in first-principles stopping power calculations, allowing 1-2 order of magnitude cost reductions for obtaining representatively averaged and converged results.

Simulations and diagnostics of high-energy-density plasmas and warm dense matter rely on models of material response properties, both static and dynamic (frequency-dependent). Here, we systematically investigate variations in dynamic electron-ion collision frequencies $\nu(\omega)$ in warm dense matter using data from a self-consistent-field average-atom model. We show that including the full quantum density of states, strong collisions, and inelastic collisions lead to significant changes in $\nu(\omega)$. These changes result in red shifts and broadening of the plasmon peak in the dynamic structure factor, an effect observable in x-ray Thomson scattering spectra, and modify stopping powers around the Bragg peak. These changes improve the agreement of computationally efficient average-atom models with first-principles time-dependent density functional theory in warm dense aluminum, carbon, and deuterium.

Bound-bound transitions can occur when localized atomic orbitals are thermally depleted, allowing excitations that would otherwise be forbidden at zero temperature. We predict signatures of bound-bound transitions in x-ray Thomson scattering measurements of laboratory-accessible warm dense conditions. Efficient average-atom models amended to include quasibound states achieve continuity of observables under ionization and agree with time-dependent density functional theory in their prediction of these scattering signatures, which hold compelling diagnostic potential for high-energy-density experiments., Comment: 6 pages, 3 figures. Preprint for PNP17

Finite-temperature Kohn--Sham density-functional theory (KS-DFT) is a widely-used method in warm dense matter (WDM) simulations and diagnostics. Unfortunately, full KS-DFT-molecular dynamics models scale unfavourably with temperature and there remains uncertainty regarding the performance of existing approximate exchange-correlation (XC) functionals under WDM conditions. Of particular concern is the expected explicit dependence of the XC functional on temperature, which is absent from most approximations. Average-atom (AA) models, which significantly reduce the computational cost of KS-DFT calculations, have therefore become an integral part of WDM modelling. In this paper, we present a derivation of a first-principles AA model from the fully-interacting many-body Hamiltonian, carefully analysing the assumptions made and terms neglected in this reduction. We explore the impact of different choices within this model -- such as boundary conditions and XC functionals -- on common properties in WDM, for example equation-of-state data, ionization degree and the behaviour of the frontier energy levels. Furthermore, drawing upon insights from ground-state KS-DFT, we discuss the likely sources of error in KS-AA models and possible strategies for mitigating such errors., Comment: To be published in Physical Review Research. 29 pages, 15 figures. Derivation of Hamiltonian clarified since v3 and added appendix B on computation of free energy

X-Ray Thomson Scattering (XRTS) is an important experimental technique used to measure the temperature, ionization state, structure, and density of warm dense matter (WDM). The fundamental property probed in these experiments is the electronic dynamic structure factor (DSF). In most models, this is decomposed into three terms [Chihara, J. Phys. F: Metal Phys. {\bf 17}, 295 (1987)] representing the response of tightly bound, loosely bound, and free electrons. Accompanying this decomposition is the classification of electrons as either bound or free, which is useful for gapped and cold systems but becomes increasingly questionable as temperatures and pressures increase into the WDM regime. In this work we provide unambiguous first principles calculations of the dynamic structure factor of warm dense beryllium, independent of the Chihara form, by treating bound and free states under a single formalism. The computational approach is real-time finite-temperature time-dependent density functional theory (TDDFT) being applied here for the first time to WDM. We compare results from TDDFT to Chihara-based calculations for experimentally relevant conditions in shock-compressed beryllium., Comment: 12 pages, 8 figures, and 1 table. 6 pages main manuscript