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

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.

Bound electron screening effect on ion-ion potential of warm and hot dense matter

The effects of bound electron screening in warm and hot dense matter are investigated analytically and a theoretical description of screened short-range repulsion is given meanwhile. An empirical ion-ion potential including the classic charge screening and chemical bond attraction at various temperatures and densities is proposed. By solving hypernetted chain equations and comparing the obtained radial distribution function (RDF) with ab initio simulations, the proposed ion-ion potential is found to be promising over a wide range of temperatures and densities for warm dense aluminum and iron. The elastic scattering amplitude and the x-ray absorption near the edge structure of warm dense aluminum calculated from the obtained RDF are in good agreement with experiment results.

Equation of state modeling with pseudoatom molecular dynamics

Using a modified version of the pseudoatom molecular-dynamics approach, the silicon and oxygen equations of state were generated and then employed to construct the equation of state of silicon dioxide. The results are supported by the close agreement with ab initio simulations of the silicon pressure and experimental shock Hugoniot of silicon dioxide. Ion thermal contributions to thermodynamic functions provided by the PAMD simulations are compared to their counterparts obtained with the one-component plasma and charged-hard-sphere approximations.

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.

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.

Equation of state of dense plasmas with pseudoatom molecular dynamics

We present an approximation for calculating the equation of state (EOS) of warm and hot dense matter that is built on the previously published pseudoatom molecular dynamics (PAMD) model of dense plasmas [Starrett et al., Phys. Rev. E 91, 013104 (2015)PLEEE81539-375510.1103/PhysRevE.91.013104]. While the EOS calculation with PAMD was previously limited to orbital-free density functional theory (DFT), the new approximation presented here allows a Kohn-Sham DFT treatment of the electrons. The resulting EOS thus includes a quantum mechanical treatment of the electrons with a self-consistent model of the ionic structure, while remaining tractable at high temperatures. The method is validated by comparisons with pressures from ab initio simulations of Be, Al, Si, and Fe. The EOS in the Thomas-Fermi approximation shows remarkable thermodynamic consistency over a wide range of temperatures for aluminum. We calculate the principal Hugoniots of aluminum and silicon up to 500 eV. We find that the ionic structure of the plasma has a modest effect that peaks at temperatures of a few eV and that the features arising from the electronic structure agree well with ab initio simulations.

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.

Ionic Transport Coefficients of Dense Plasmas without Molecular Dynamics

We present a theoretical model that allows a fast and accurate evaluation of ionic transport properties of realistic plasmas spanning from warm and dense to hot and dilute conditions, including mixtures. This is achieved by combining a recent kinetic theory based on effective interaction potentials with a model for the equilibrium radial density distribution based on an average atom model and the integral equations theory of fluids. The model should find broad use in applications where nonideal plasma conditions are traversed, including inertial confinement fusion, compact astrophysical objects, solar and extrasolar planets, and numerous present-day high energy density laboratory experiments.