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

Ionization calculations using classical molecular dynamics

By performing an ensemble of molecular dynamics simulations, the model-dependent ionization state is computed for strongly interacting systems self-consistently. This is accomplished through a free energy minimization framework based on the technique of thermodynamic integration. To illustrate the method, two simple models applicable to partially ionized hydrogen plasma are presented in which pair potentials are employed between ions and neutral particles. Within the models, electrons are either bound in the hydrogen ground state or distributed in a uniform charge-neutralizing background. Particular attention is given to the transition between atomic gas and ionized plasma, where the effect of neutral interactions is explored beyond commonly used models in the chemical picture. Furthermore, pressure ionization is observed when short-range repulsion effects are included between neutrals. The developed technique is general, and we discuss the applicability to a variety of molecular dynamics models for partially ionized warm dense matter.

Quantifying ionization in hot dense plasmas

Ionization is a problematic quantity in that it does not have a well-defined thermodynamic definition and yet it is a key parameter within plasma modeling. One still therefore aims to find a consistent and unambiguous definition for the ionization state. Within this context we present finite-temperature density functional theory calculations of the ionization state of carbon in CH plasmas using two potential definitions: one based on counting the number of continuum electrons, and another based on the optical conductivity. Differences of up to 10% are observed between the two methods. However, including "Pauli forbidden" transitions in the conductivity reproduces the counting definition, suggesting such transitions are important to evaluate the ionization state.

Investigating mechanisms of state localization in highly ionized dense plasmas

We present experimental observations of K_{β} emission from highly charged Mg ions at solid density, driven by intense x rays from a free electron laser. The presence of K_{β} emission indicates the n=3 atomic shell is relocalized for high charge states, providing an upper constraint on the depression of the ionization potential. We explore the process of state relocalization in dense plasmas from first principles using finite-temperature density functional theory alongside a wave-function localization metric, and find excellent agreement with experimental results.

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.

Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics

Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/cm^{3} to 16 g/cm^{3} and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method. We find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the sesame database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT-MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.