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

Development of a new quantum trajectory molecular dynamics framework

An extension to the wave packet description of quantum plasmas is presented, where the wave packet can be elongated in arbitrary directions. A generalized Ewald summation is constructed for the wave packet models accounting for long-range Coulomb interactions and fermionic effects are approximated by purpose-built Pauli potentials, self-consistent with the wave packets used. We demonstrate its numerical implementation with good parallel support and close to linear scaling in particle number, used for comparisons with the more common wave packet employing isotropic states. Ground state and thermal properties are compared between the models with differences occurring primarily in the electronic subsystem. Especially, the electrical conductivity of dense hydrogen is investigated where a 15% increase in DC conductivity can be seen in our wave packet model compared with other models. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Equivalence between pressure- and structure-defined ionization in hot dense carbon

The determination of the ionization of a system in the hot dense regime is a long standing issue. Recent studies have shown inconsistencies between standard predictions using average atom models and evaluations deduced from electronic transport properties computed with quantum molecular dynamics simulations [Bethkenhagen et al., Phys. Rev. Res. 2, 023260 (2020)]2643-156410.1103/PhysRevResearch.2.023260. Here, we propose a definition of the ionization based on its effect on the plasma structure as given by the pair distribution function (PDF), and on the concept of effective one-component plasma (eOCP). We also introduce a definition based on the total pressure and on a modelization of the electronic pressure. We show the equivalence of these definitions on two studies of carbon along the 100 eV isotherm and the 10 g/cm^{3} isochor. Simulations along the 100 eV isotherm are obtained with the newly implemented Ext. First principles molecular dynamics (Fpmd) method in Abinit for densities ranging from 1 to 500 g/cm^{3}and along the 10 g/cm^{3} isochor with the recently published Spectral quadrature DFT (Sqdft) simulations, between 8 and 860 eV. The resulting ionizations are compared to the predictions of the average-atom code Qaam which is based on the muffin-tin approximation. A disagreement between the eOCP and the actual PDFs (non-OCP behavior) is interpreted as the onset of bonding in the system.

An Investigation into the Approximations Used in Wave Packet Molecular Dynamics for the Study of Warm Dense Matter

Wave packet molecular dynamics (WPMD) has recently received a lot of attention as a computationally fast tool with which to study dynamical processes in warm dense matter beyond the Born–Oppenheimer approximation. These techniques, typically, employ many approximations to achieve computational efficiency while implementing semi-empirical scaling parameters to retain accuracy. We investigated three of the main approximations ubiquitous to WPMD: a restricted basis set, approximations to exchange, and the lack of correlation. We examined each of these approximations in regard to atomic and molecular hydrogen in addition to a dense hydrogen plasma. We found that the biggest improvement to WPMD comes from combining a two-Gaussian basis with a semi-empirical correction based on the valence-bond wave function. A single parameter scales this correction to match experimental pressures of dense hydrogen. Ultimately, we found that semi-empirical scaling parameters are necessary to correct for the main approximations in WPMD. However, reducing the scaling parameters for more ab-initio terms gives more accurate results and displays the underlying physics more readily.

First-Principles Determination of Electron-Ion Couplings in the Warm Dense Matter Regime

We present first-principles calculations of the rate of energy exchanges between electrons and ions in nonequilibrium warm dense plasmas, liquid metals, and hot solids, a fundamental property for which various models offer diverging predictions. To this end, a Kubo relation for the electron-ion coupling parameter is introduced, which includes self-consistently the quantum, thermal, nonlinear, and strong coupling effects that coexist in materials at the confluence of solids and plasmas. Most importantly, like other Kubo relations widely used for calculating electronic conductivities, the expression can be evaluated using quantum molecular dynamics simulations. Results are presented and compared to experimental and theoretical predictions for representative materials of various electronic complexity, including aluminum, copper, iron, and nickel.

Dynamical structure factor of strongly coupled ions in a dense quantum plasma

The dynamical structure factor (DSF) of strongly coupled ions in dense plasmas with partially and strongly degenerate electrons is investigated. The main focus is on the impact of electronic correlations (nonideality) on the ionic DSF. The latter is computed by carrying out molecular dynamics (MD) simulations with a screened ion-ion interaction potential. The electronic screening is taken into account by invoking the Singwi-Tosi-Land-Sjölander approximation, and it is compared to the MD simulation data obtained considering the electronic screening in the random phase approximation and using the Yukawa potential. We find that electronic correlations lead to lower values of the ion-acoustic mode frequencies and to an extension of the applicability limit with respect to the wave-number of a hydrodynamic description. Moreover, we show that even in the limit of weak electronic coupling, electronic correlations have a nonnegligible impact on the ionic longitudinal sound speed. Additionally, the applicability of the Yukawa potential with an adjustable screening parameter is discussed, which will be of interest, e.g., for the interpretation of experimental results for the ionic DSF of dense plasmas.

Structural characteristics of strongly coupled ions in a dense quantum plasma

The structural properties of strongly coupled ions in dense plasmas with moderately to strongly degenerate electrons are investigated in the framework of the one-component plasma model of ions interacting through a screened pair interaction potential. Special focus is put on the description of the electronic screening in the Singwi-Tosi-Land-Sjölander (STLS) approximation. Different cross-checks and analyses using ion potentials obtained from ground-state quantum Monte Carlo data, the random phase approximation (RPA), and existing analytical models are presented for the computation of the structural properties, such as the pair distribution and the static structure factor, of strongly coupled ions. The results are highly sensitive to the features of the screened pair interaction potential. This effect is particularly visible in the static structure factor. The applicability range of the screened potential computed from STLS is identified in terms of density and temperature of the electrons. It is demonstrated that at r_{s}>1, where r_{s} is the ratio of the mean interelectronic distance to the Bohr radius, electronic correlations beyond RPA have a nonnegligible effect on the structural properties. Additionally, the applicability of the hypernetted chain approximation for the calculation of the structural properties using the screened pair interaction potential is analyzed employing the effective coupling parameter approach.

Theory of Thomson scattering in inhomogeneous plasmas

A self-consistent kinetic theory of Thomson scattering of an electromagnetic field by a nonuniform plasma is derived. We show that not only the imaginary part, but also the time and space derivatives of the real part of the dielectric susceptibility determine the amplitude and the width of the Thomson scattering spectral lines. As a result of inhomogeneity, these properties become asymmetric with respect to inversion of the sign of the frequency. Our theory provides a method of a remote probing and measurement of electron density gradients in plasma; this is based on the demonstrated asymmetry of the Thomson scattering lines.

Ion-ion dynamic structure factor, acoustic modes, and equation of state of two-temperature warm dense aluminum

The ion-ion dynamical structure factor and the equation of state of warm dense aluminum in a two-temperature quasiequilibrium state, with the electron temperature higher than the ion temperature, are investigated using molecular-dynamics simulations based on ion-ion pair potentials constructed from a neutral pseudoatom model. Such pair potentials based on density functional theory are parameter-free and depend directly on the electron temperature and indirectly on the ion temperature, enabling efficient computation of two-temperature properties. Comparison with ab initio simulations and with other average-atom calculations for equilibrium aluminum shows good agreement, justifying a study of quasiequilibrium situations. Analyzing the van Hove function, we find that ion-ion correlations vanish in a time significantly smaller than the electron-ion relaxation time so that dynamical properties have a physical meaning for the quasiequilibrium state. A significant increase in the speed of sound is predicted from the modification of the dispersion relation of the ion acoustic mode as the electron temperature is increased. The two-temperature equation of state including the free energy, internal energy, and pressure is also presented.

Crossover from Classical to Fermi Liquid Behavior in Dense Plasmas

We explore the crossover from classical plasma to quantum Fermi liquid behavior of electrons in dense plasmas. To this end, we analyze the evolution with density and temperature of the momentum lifetime of a test electron introduced in a dense electron gas. This allows us (1) to determine the boundaries of the crossover region in the temperature-density plane and to shed light on the evolution of scattering properties across it, (2) to quantify the role of the fermionic nature of electrons on electronic collisions across the crossover region, and (3) to explain how the concept of the Coulomb logarithm emerges at a high enough temperature but disappears at a low enough temperature.

Bayesian inference of x-ray diffraction spectra from warm dense matter with the one-component-plasma model

We show that the Bayesian inference of recently measured x-ray diffraction spectra from laser-shocked aluminum [L. B. Fletcher et al., Nat. Photon. 9, 274 (2015)10.1038/nphoton.2015.41] with the one-component-plasma (OCP) model performs remarkably well at estimating the ionic density and temperature. This statistical approach requires many evaluations of the OCP static structure factor, which were done using a recently derived analytic fit. The atomic form factor is approximated by an exponential function in the diffraction window of the first peak. The electronic temperature is then estimated from a comparison of this approximated form factor with the electronic structure of an average atom model. Out-of-equilibrium states, with electrons hotter than ions, are diagnosed for the spectra obtained early after the pump, whereas at a late time delay the plasma is at thermal equilibrium. Apart from the present findings, this OCP-based modeling of warm dense matter has an important role to play in the interpretation of x-ray Thomson scattering measurements currently performed at large laser facilities.

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.

Dynamic conductivity and plasmon profile of aluminum in the ultra-fast-matter regime

We use an explicitly isochoric two-temperature theory to analyze recent x-ray laser scattering data for aluminum in the ultra-fast-matter (UFM) regime up to 6 eV. The observed surprisingly low conductivities are explained by including strong electron-ion scattering effects using the phase shifts calculated via the neutral-pseudo-atom model. The difference between the static conductivity for UFM-Al and equilibrium aluminum in the warm-dense matter state is clearly brought out by comparisons with available density-fucntional+molecular-dynamics simulations. Thus the applicability of the Mermin model to UFM is questioned. The static and dynamic conductivity, collision frequency, and the plasmon line shape, evaluated within the simplest Born approximation for UFM aluminum, are in good agreement with experiment.

Unified Concept of Effective One Component Plasma for Hot Dense Plasmas

Orbital-free molecular dynamics simulations are used to benchmark two popular models for hot dense plasmas: the one component plasma (OCP) and the Yukawa model. A unified concept emerges where an effective OCP (EOCP) is constructed from the short-range structure of the plasma. An unambiguous ionization and the screening length can be defined and used for a Yukawa system, which reproduces the long-range structure with finite compressibility. Similarly, the dispersion relation of longitudinal waves is consistent with the screened model at vanishing wave number but merges with the OCP at high wave number. Additionally, the EOCP reproduces the overall relaxation time scales of the correlation functions associated with ionic motion. In the hot dense regime, this unified concept of EOCP can be fruitfully applied to deduce properties such as the equation of state, ionic transport coefficients, and the ion feature in x-ray Thomson scattering experiments.

Strongly-coupled plasmas formed from laser-heated solids

We present an analysis of ion temperatures in laser-produced plasmas formed from solids with different initial lattice structures. We show that the equilibrium ion temperature is limited by a mismatch between the initial crystallographic configuration and the close-packed configuration of a strongly-coupled plasma, similar to experiments in ultracold neutral plasmas. We propose experiments to demonstrate and exploit this crystallographic heating in order to produce a strongly coupled plasma with a coupling parameter of several hundred.

Microscopic nonlinear relativistic quantum theory of absorption of powerful x-ray radiation in plasma

The microscopic quantum theory of plasma nonlinear interaction with the coherent shortwave electromagnetic radiation of arbitrary intensity is developed. The Liouville-von Neumann equation for the density matrix is solved analytically considering a wave field exactly and a scattering potential of plasma ions as a perturbation. With the help of this solution we calculate the nonlinear inverse-bremsstrahlung absorption rate for a grand canonical ensemble of electrons. The latter is studied in Maxwellian, as well as in degenerate quantum plasma for x-ray lasers at superhigh intensities and it is shown that one can achieve the efficient absorption coefficient in these cases.

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.

Ab initio calculation of the ion feature in x-ray Thomson scattering

The spectrum of x-ray Thomson scattering is proportional to the dynamic structure factor. An important contribution is the ion feature which describes elastic scattering of x rays off electrons. We apply an ab initio method for the calculation of the form factor of bound electrons, the slope of the screening cloud of free electrons, and the ion-ion structure factor in warm dense beryllium. With the presented method we can calculate the ion feature from first principles. These results will facilitate a better understanding of x-ray scattering in warm dense matter and an accurate measurement of ion temperatures which would allow determining nonequilibrium conditions, e.g., along shock propagation.

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.