Predictions of x-ray scattering spectra for warm dense matter

X-ray Thomson scattering absolute intensity from the f-sum rule in the imaginary-time domain

We present a formally exact and simulation-free approach for the normalization of X-ray Thomson scattering (XRTS) spectra based on the f-sum rule of the imaginary-time correlation function (ITCF). Our method works for any degree of collectivity, over a broad range of temperatures, and is applicable even in nonequilibrium situations. In addition to giving us model-free access to electronic correlations, this new approach opens up the intriguing possibility to extract a plethora of physical properties from the ITCF based on XRTS experiments.

Plasma ionization balance in chemical-picture and average-atom models

We propose an approximate method to calculate ion partition functions in the context of the chemical-picture representation of plasmas as an interacting mixture of various ions and free electrons under the local-thermodynamic-equilibrium conditions. The method uses the superconfiguration approach and implies that the first-order corrections to the energies of excited electron configurations due to the electron-electron interaction may be replaced by a similar first-order correction to the energy of the basic configuration of an ion with the same number of bound electrons. The method enables one to significantly speed up the calculations and generally provides quite accurate results. Using the method proposed, plasma ionization balance and average ion charges calculated on the base of the chemical-picture representation show a good agreement with the relevant average-atom data. For the case of weak electron-ion nonideality, we provide approximate relations between the chemical-picture and average-atom values of the average ion charge, chemical potential, and plasma-density depression of ionization potential.

Observing the onset of pressure-driven K-shell delocalization

The gravitational pressure in many astrophysical objects exceeds one gigabar (one billion atmospheres)^1-3, creating extreme conditions where the distance between nuclei approaches the size of the K shell. This close proximity modifies these tightly bound states and, above a certain pressure, drives them into a delocalized state⁴. Both processes substantially affect the equation of state and radiation transport and, therefore, the structure and evolution of these objects. Still, our understanding of this transition is far from satisfactory and experimental data are sparse. Here we report on experiments that create and diagnose matter at pressures exceeding three gigabars at the National Ignition Facility⁵ where 184 laser beams imploded a beryllium shell. Bright X-ray flashes enable precision radiography and X-ray Thomson scattering that reveal both the macroscopic conditions and the microscopic states. The data show clear signs of quantum-degenerate electrons in states reaching 30 times compression, and a temperature of around two million kelvins. At the most extreme conditions, we observe strongly reduced elastic scattering, which mainly originates from K-shell electrons. We attribute this reduction to the onset of delocalization of the remaining K-shell electron. With this interpretation, the ion charge inferred from the scattering data agrees well with ab initio simulations, but it is significantly higher than widely used analytical models predict⁶.

Average-atom model with Siegert states

In plasmas, electronic states can be well-localized bound states or itinerant free states, or something in between. In self-consistent treatments of plasma electronic structure such as the average-atom model, all states must be accurately resolved in order to achieve a converged numerical solution. This is a challenging numerical and algorithmic problem in large part due to the continuum of free states which is relatively expensive and difficult to resolve accurately. Siegert states are an appealing alternative. They form a complete eigenbasis with a purely discrete spectrum while still being equivalent to a representation in terms of the usual bound states and free states. However, many of their properties are unintuitive, and it is not obvious that they are suitable for self-consistent plasma electronic structure calculations. Here it is demonstrated that Siegert states can be used to accurately solve an average-atom model and offer advantages over the traditional finite-difference approach, including a concrete physical picture of pressure ionization and continuum resonances.

Effective Static Approximation: A Fast and Reliable Tool for Warm-Dense Matter Theory

We present an effective static approximation (ESA) to the local field correction (LFC) of the electron gas that enables highly accurate calculations of electronic properties like the dynamic structure factor S(q,ω), the static structure factor S(q), and the interaction energy v. The ESA combines the recent neural-net representation by T. Dornheim et al., [J. Chem. Phys. 151, 194104 (2019)JCPSA60021-960610.1063/1.5123013] of the temperature-dependent LFC in the exact static limit with a consistent large wave-number limit obtained from quantum Monte Carlo data of the on-top pair distribution function g(0). It is suited for a straightforward integration into existing codes. We demonstrate the importance of the LFC for practical applications by reevaluating the results of the recent x-ray Thomson scattering experiment on aluminum by Sperling et al. [Phys. Rev. Lett. 115, 115001 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.115001]. We find that an accurate incorporation of electronic correlations in terms of the ESA leads to a different prediction of the inelastic scattering spectrum than obtained from state-of-the-art models like the Mermin approach or linear-response time-dependent density functional theory. Furthermore, the ESA scheme is particularly relevant for the development of advanced exchange-correlation functionals in density functional theory.

Structure and dynamics of warm dense aluminum: a molecular dynamics study with density functional theory and deep potential

We perform a systematic study on the structure and dynamics of warm dense aluminum (Al) at temperatures ranging from 0.5 to 5.0 eV with molecular dynamics utilizing both density functional theory (DFT) and the deep potential (DP) method. On one hand, unlike the Thomas-Fermi kinetic energy density functional (KEDF), we find that the orbital-free DFT method with the Wang-Teter non-local KEDF yields properties of warm dense Al that agree well with the Kohn-Sham DFT method, enabling accurate orbital-free DFT simulations of warm dense Al at relatively low temperatures. On the other hand, the DP method constructs a deep neural network that has a high accuracy in reproducing short- and long-ranged properties of warm dense Al when compared to the DFT methods. The DP method is orders of magnitudes faster than DFT and is well-suited for simulating large systems and long trajectories to yield accurate properties of warm dense Al. Our results suggest that the combination of DFT methods and the DP model is a powerful tool for accurately and efficiently simulating warm dense matter.

Permutation-blocking path-integral Monte Carlo approach to the static density response of the warm dense electron gas

The static density response of the uniform electron gas is of fundamental importance for numerous applications. Here we employ the recently developed ab initio permutation blocking path integral Monte Carlo (PB-PIMC) technique [T. Dornheim et al., New J. Phys. 17, 073017 (2015)10.1088/1367-2630/17/7/073017] to carry out extensive simulations of the harmonically perturbed electron gas at warm dense matter conditions. In particular, we investigate in detail the validity of linear response theory and demonstrate that PB-PIMC allows us to obtain highly accurate results for the static density response function and, thus, the static local field correction. A comparison with dielectric approximations to our new ab initio data reveals the need for an exact treatment of correlations. Finally, we consider a superposition of multiple perturbations and discuss the implications for the calculation of the static response function.

Pair potentials for warm dense matter and their application to x-ray Thomson scattering in aluminum and beryllium

Ultrafast laser experiments yield increasingly reliable data on warm dense matter, but their interpretation requires theoretical models. We employ an efficient density functional neutral-pseudoatom hypernetted-chain (NPA-HNC) model with accuracy comparable to ab initio simulations and which provides first-principles pseudopotentials and pair potentials for warm-dense matter. It avoids the use of (i) ad hoc core-repulsion models and (ii) "Yukawa screening" and (iii) need not assume ion-electron thermal equilibrium. Computations of the x-ray Thomson scattering (XRTS) spectra of aluminum and beryllium are compared with recent experiments and with density-functional-theory molecular-dynamics (DFT-MD) simulations. The NPA-HNC structure factors, compressibilities, phonons, and conductivities agree closely with DFT-MD results, while Yukawa screening gives misleading results. The analysis of the XRTS data for two of the experiments, using two-temperature quasi-equilibrium models, is supported by calculations of their temperature relaxation times.

Average-atom treatment of relaxation time in x-ray Thomson scattering from warm dense matter

The influence of finite relaxation times on Thomson scattering from warm dense plasmas is examined within the framework of the average-atom approximation. Presently most calculations use the collision-free Lindhard dielectric function to evaluate the free-electron contribution to the Thomson cross section. In this work, we use the Mermin dielectric function, which includes relaxation time explicitly. The relaxation time is evaluated by treating the average atom as an impurity in a uniform electron gas and depends critically on the transport cross section. The calculated relaxation rates agree well with values inferred from the Ziman formula for the static conductivity and also with rates inferred from a fit to the frequency-dependent conductivity. Transport cross sections determined by the phase-shift analysis in the average-atom potential are compared with those evaluated in the commonly used Born approximation. The Born approximation converges to the exact cross sections at high energies; however, differences that occur at low energies lead to corresponding differences in relaxation rates. The relative importance of including relaxation time when modeling x-ray Thomson scattering spectra is examined by comparing calculations of the free-electron dynamic structure function for Thomson scattering using Lindhard and Mermin dielectric functions. Applications are given to warm dense Be plasmas, with temperatures ranging from 2 to 32 eV and densities ranging from 2 to 64 g/cc.

X-ray Thomson Scattering in Warm Dense Matter without the Chihara Decomposition

X-ray Thomson scattering 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. In most models, this is decomposed into three terms [J. Chihara, J. Phys. F 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.

Temperature relaxation in dense plasmas

We present a model to calculate temperature-relaxation rates in dense plasmas. The electron-ion interaction potential and the thermodynamic data of interest are provided by an average-atom model. This approach allows the study of the temperature relaxation in a two-temperature electron-ion system.

Models of the elastic x-ray scattering feature for warm dense aluminum

The elastic feature of x-ray scattering from warm dense aluminum has recently been measured by Fletcher et al. [Nature Photonics 9, 274 (2015)]10.1038/nphoton.2015.41 with much higher accuracy than had hitherto been possible. This measurement is a direct test of the ionic structure predicted by models of warm dense matter. We use the method of pseudoatom molecular dynamics to predict this elastic feature for warm dense aluminum with temperatures of 1-100 eV and densities of 2.7-8.1g/cm^{3}. We compare these predictions to experiments, finding good agreement with Fletcher et al. and corroborating the discrepancy found in analyses of an earlier experiment of Ma et al. [Phys. Rev. Lett. 110, 065001 (2013)]PRLTAO0031-900710.1103/PhysRevLett.110.065001. We also evaluate the validity of the Thomas-Fermi model of the electrons and of the hypernetted chain approximation in computing the elastic feature and find them both wanting in the regime currently probed by experiments.

Atomic and optical properties of warm dense copper

The emission of x rays from warm dense matter is of great interest for both spectroscopic diagnostics and development of intense x-ray sources. We report the results from the collisional-radiative steady-state (CRSS) modeling of atomic and optical properties of copper plasmas at near-solid and solid-state density for a range of temperatures. The CRSS model is validated against the available data on the average charge state and shifts of energy levels in aluminum and the opacity and emissivity spectra of carbon and aluminum plasmas. The average charge states, number density of ion species, and free electrons as a function of temperature are investigated for the solid-density copper plasma. Due to the dense plasma environment the four outer electrons are found to be unbounded even in the low-temperature limit ∼1eV. As the temperature changes from 1 to 100 eV, the predominant species vary from fivefold- to twelvefold-ionized copper ions. The opacity and emissivity spectra of dense copper plasmas are studied using the local thermodynamic equilibrium (LTE) and non-LTE approaches. It is found that the non-LTE effects are important in the spectral region of soft x rays emitted from the K shell. The emissivity in spectral lines is completely suppressed, indicating the importance of the energy-dissipating radiative processes in this soft x-ray region. Line broadening and redshifts of the K- and L-shell spectral lines toward higher wavelengths are observed with the increase of plasma density. These results have important implications for understanding the radiative properties of warm dense copper and can be useful for future experimental studies.

Out-of-equilibrium conditions in x-ray Thomson scattering experiments

We study out-of-equilibrium conditions in recent x-ray Thomson scattering experiments performed in warm dense matter. We use an effective one-component plasma model to characterize the states in which electron and ion temperatures are different. An estimation of the ion temperature is obtained. This method is tested against two recent experiments. Strong out-of-equilibrium conditions are found.

Ion-ion dynamic structure factor of warm dense mixtures

The ion-ion dynamic structure factor of warm dense matter is determined using the recently developed pseudoatom molecular dynamics method [Starrett et al., Phys. Rev. E 91, 013104 (2015)]. The method uses density functional theory to determine ion-ion pair interaction potentials that have no free parameters. These potentials are used in classical molecular dynamics simulations. This constitutes a computationally efficient and realistic model of dense plasmas. Comparison with recently published simulations of the ion-ion dynamic structure factor and sound speed of warm dense aluminum finds good to reasonable agreement. Using this method, we make predictions of the ion-ion dynamical structure factor and sound speed of a warm dense mixture-equimolar carbon-hydrogen. This material is commonly used as an ablator in inertial confinement fusion capsules, and our results are amenable to direct experimental measurement.

Average-atom model combined with the hypernetted chain approximation applied to warm dense matter

We have combined the average-atom model with the hypernetted chain approximation (AAHNC) to describe the electronic and ionic structure in the warm dense matter regime. On the basis of the electronic and ionic structures, the x-ray Thomson scattering (XRTS) spectrum is calculated using the random-phase approximation. While the electronic structure is described within the average-atom model, the effects of other ions on the electronic structure are considered using an integral equation method of the theory of liquids, namely the hypernetted chain approximation. The ion-ion pair potential is calculated using the modified Gordon-Kim model based on the electronic density distribution. Finally, the electronic and ionic structures are determined self-consistently. The XRTS spectrum is calculated according to the Chihara formula, where the scattering contributions are divided into three components: elastic, bound-free, and free-free. Comparison of the present AAHNC results with other theoretical models and experimental data shows very good agreement. Thus the AAHNC model can give a reasonable description of the electronic and ionic structure in warm dense matter.

Evidence for out-of-equilibrium states in warm dense matter probed by x-ray Thomson scattering

A recent and unexpected discrepancy between ab initio simulations and the interpretation of a laser shock experiment on aluminum, probed by x-ray Thomson scattering (XRTS), is addressed. The ion-ion structure factor deduced from the XRTS elastic peak (ion feature) is only compatible with a strongly coupled out-of-equilibrium state. Orbital free molecular dynamics simulations with ions colder than the electrons are employed to interpret the experiment. The relevance of decoupled temperatures for ions and electrons is discussed. The possibility that it mimics a transient, or metastable, out-of-equilibrium state after melting is also suggested.

Resistivity saturation in warm dense matter

Electrical resistivity is shown to saturate in solid-density aluminum in the warm dense matter regime. Calculations are done using the average-atom model SCAALP and the finite-temperature Ziman-Evans formula for electrical resistivity. The mean free path is estimated using the Drude law. This mean free path is shown to present a minimum of the order of the interatomic spacing.

Integral equation model for warm and hot dense mixtures

In a previous work [C. E. Starrett and D. Saumon, Phys. Rev. E 87, 013104 (2013)] a model for the calculation of electronic and ionic structures of warm and hot dense matter was described and validated. In that model the electronic structure of one atom in a plasma is determined using a density-functional-theory-based average-atom (AA) model and the ionic structure is determined by coupling the AA model to integral equations governing the fluid structure. That model was for plasmas with one nuclear species only. Here we extend it to treat plasmas with many nuclear species, i.e., mixtures, and apply it to a carbon-hydrogen mixture relevant to inertial confinement fusion experiments. Comparison of the predicted electronic and ionic structures with orbital-free and Kohn-Sham molecular dynamics simulations reveals excellent agreement wherever chemical bonding is not significant.

Linear and quadratic static response functions and structure functions in Yukawa liquids

We compute linear and quadratic static density response functions of three-dimensional Yukawa liquids by applying an external perturbation potential in molecular dynamics simulations. The response functions are also obtained from the equilibrium fluctuations (static structure factors) in the system via the fluctuation-dissipation theorems. The good agreement of the quadratic response functions, obtained in the two different ways, confirms the quadratic fluctuation-dissipation theorem. We also find that the three-point structure function may be factorizable into two-point structure functions, leading to a cluster representation of the equilibrium triplet correlation function.