Continuum Lowering and Fermi-Surface Rising in Strongly Coupled and Degenerate Plasmas

Nonideal effects on ionization potential depression and ionization balance in dense Al and Au plasmas

For low-density plasmas, the ionization balance can be properly described by the normal Saha equation in the chemical picture. For dense plasmas, however, nonideal effects due to the interactions between the electrons and ions and among the electrons themselves affect the ionization potential depression and the ionization balance. With the increasing of plasma density, the pressure ionization starts to play a more obvious role and competes with the thermal ionization. Based on a local-density temperature-dependent ion-sphere model, we develop a unified and self-consistent theoretical formalism to simultaneously investigate the ionization potential depression, the ionization balance, and the charge states distributions of the dense plasmas. In this work, we choose Al and Au plasmas as examples as Al is a prototype light element and Au is an important heavy element in many research fields such as in the inertial confinement fusion. The nonideal effect of the free electrons in the plasmas is considered by the single-electron effective potential contributed by both the bound electrons of different charge states and the free electrons in the plasmas. For the Al plasmas, we can reconcile the results of two experiments on measuring the ionization potential depression, in which one experiment can be better explained by the Stewart-Pyatt model while the other fits better with the Ecker-Kröll model. For dense Au plasmas, the results show that the double peak structure of the charge state distribution appears to be a common phenomenon. In particular, the calculated ionization balance shows that the two- and three-peak structures can appear simultaneously for denser Au plasmas above ∼30g/cm^{3}.

Bound-State Breaking and the Importance of Thermal Exchange-Correlation Effects in Warm Dense Hydrogen

Hydrogen at extreme temperatures and pressures is of key relevance for cutting-edge technological applications, with inertial confinement fusion research being a prime example. In addition, it is ubiquitous throughout our universe and naturally occurs in a variety of astrophysical objects. In the present work, we present exact ab initio path integral Monte Carlo (PIMC) results for the electronic density of warm dense hydrogen along a line of constant degeneracy across a broad range of densities. Using the well-known concept of reduced density gradients, we develop a new framework to identify the breaking of bound states due to pressure ionization in bulk hydrogen. Moreover, we use our PIMC results as a reference to rigorously assess the accuracy of a variety of exchange-correlation (XC) functionals in density functional theory calculations for different density regions. Here, a key finding is the importance of thermal XC effects for the accurate description of density gradients in high-energy-density systems. Our exact PIMC test set is freely available online and can be used to guide the development of new methodologies for the simulation of warm dense matter and beyond.

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.

First-principles study of L-shell iron and chromium opacity at stellar interior temperatures

Recently developed free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter has been applied for studying L-shell opacity of iron and chromium at T=182 eV. We use Mermin-Kohn-Sham density functional theory with a ground-state and a fully-temperature-dependent generalized gradient approximation exchange-correlation (XC) functionals. It is demonstrated that the role of XC at such a high-T is negligible due to the total free energy of interacting systems being dominated by the noninteracting free-energy term, in agreement with estimations for the homogeneous electron gas. Our DFT predictions are compared with the radiative emissivity and opacity of the dense plasma model, with the real-space Green's function method, and with experimental measurements. Good agreement is found between all three theoretical methods, and in the bound-continuum region for Cr when compared with the experiment, while the discrepancy between direct DFT calculations and the experiment for Fe remains essentially the same as for plasma-physics models.

Probing atomic physics at ultrahigh pressure using laser-driven implosions

Spectroscopic measurements of dense plasmas at billions of atmospheres provide tests to our fundamental understanding of how matter behaves at extreme conditions. Developing reliable atomic physics models at these conditions, benchmarked by experimental data, is crucial to an improved understanding of radiation transport in both stars and inertial fusion targets. However, detailed spectroscopic measurements at these conditions are rare, and traditional collisional-radiative equilibrium models, based on isolated-atom calculations and ad hoc continuum lowering models, have proved questionable at and beyond solid density. Here we report time-integrated and time-resolved x-ray spectroscopy measurements at several billion atmospheres using laser-driven implosions of Cu-doped targets. We use the imploding shell and its hot core at stagnation to probe the spectral changes of Cu-doped witness layer. These measurements indicate the necessity and viability of modeling dense plasmas with self-consistent methods like density-functional theory, which impact the accuracy of radiation transport simulations used to describe stellar evolution and the design of inertial fusion targets.

Atoms and molecules under extreme temperature and pressure can be investigated using dense plasmas achieved by laser-driven implosion. Here the authors report spectral change of copper in billions atmosphere pressure that can only be explained by a self-consistent approach.

First-principles equation of state of CHON resin for inertial confinement fusion applications

A wide-range (0 to 1044.0 g/cm^{3} and 0 to 10^{9} K) equation-of-state (EOS) table for a CH_{1.72}O_{0.37}N_{0.086} quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. The simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

Multi-Configuration Calculation of Ionization Potential Depression

The modelling of ionization potential depression in warm and hot dense plasmas constitutes a real theoretical challenge due to ionic coupling and electron degeneracy effects. In this work, we present a quantum statistical model based on a multi-configuration description of the electronic structure in the framework of Density Functional Theory. We discuss different conceptual issues inherent to the definition of ionization potential depression and compare our results with the famous and widely-used Ecker-Kröll and Stewart-Pyatt models.

Plasma environmental effects in the atomic structure for simulating x-ray free-electron-laser-heated solid-density matter

High energy density (HED) matter exists extensively in the Universe, and it can be created with extreme conditions in laboratory facilities such as x-ray free-electron lasers (XFEL). In HED matter, the electronic structure of individual atomic ions is influenced by a dense plasma environment, and one of the most significant phenomena is the ionization potential depression (IPD). Incorporation of the IPD effects is of great importance in accurate modeling of dense plasmas. All theoretical treatments of IPD so far have been based on the assumption of local thermodynamic equilibrium, but its validity is questionable in ultrafast formation dynamics of dense plasmas, particularly when interacting with intense XFEL pulses. A treatment of transient IPD, based on an electronic-structure calculation of an atom in the presence of a plasma environment described by classical particles, has recently been proposed [Phys. Rev. E 103, 023203 (2021)2470-004510.1103/PhysRevE.103.023203], but its application to and impact on plasma dynamics simulations have not been investigated yet. In this work, we extend XMDYN, a hybrid quantum-classical approach combining Monte Carlo and molecular dynamics, by incorporating the proposed IPD treatment into plasma dynamics simulations. We demonstrate the importance of the IPD effects in theoretical modeling of aluminum dense plasmas by comparing two XMDYN simulations: one with electronic-structure calculations of isolated atoms (without IPD) and the other with those of atoms embedded in a plasma (with IPD). At equilibrium, the mean charge obtained in the plasma simulation with IPD is in good agreement with the full quantum-mechanical average-atom model. The present approach promises to be a reliable tool to simulate the creation and nonequilibrium evolution of dense plasmas induced by ultraintense and ultrashort XFEL pulses.

The Strong Enhancement of Electron-Impact Ionization Processes in Dense Plasma by Transient Spatial Localization

Recent experiments have observed much higher electron–ion collisional ionization cross sections and rates in dense plasmas than predicted by the current standard atomic collision theory, including the plasma screening effect. We suggest that the use of (distorted) plane waves for incident and scattered electrons is not adequate to describe the dissipation that occurs during the ionization event. Random collisions with free electrons and ions in plasma cause electron matter waves to lose their phase, which results in the partial decoherence of incident and scattered electrons. Such a plasma-induced transient spatial localization of the continuum electron states significantly modifies the wave functions of continuum electrons, resulting in a strong enhancement of the electron–ion collisional ionization of ions in plasma compared to isolated ions. Here, we develop a theoretical formulation to calculate the differential and integral cross sections by incorporating the effects of plasma screening and transient spatial localization. The approach is then used to investigate the electron-impact ionization of ions in solid-density magnesium plasma, yielding results that are consistent with experiments. In dense plasma, the correlation of continuum electron energies is modified, and the integral cross sections and rates increase considerably. For the ionization of Mg9+e+1s22s2S→1s21S+2e, the ionization cross sections increase several-fold, and the rates increase by one order of magnitude. Our findings provide new insight into collisional ionization and three-body recombination and may aid investigations of the transport properties and nonequilibrium evolution of dense plasma.

Dense plasma opacity via the multiple-scattering method

The calculation of the optical properties of hot dense plasmas with a model that has self-consistent plasma physics is a grand challenge for high energy density science. Here we exploit a recently developed electronic structure model that uses multiple scattering theory to solve the Kohn-Sham density functional theory equations for dense plasmas. We calculate opacities in this regime, validate the method, and apply it to recent experimental measurements of opacity for Cr, Ni, and Fe. Good agreement is found in the quasicontinuum region for Cr and Ni, while the self-consistent plasma physics of the approach cannot explain the observed difference between models and the experiment for Fe.

Local field correction to ionization potential depression of ions in warm or hot dense matter

An analytical self-consistent approach was recently established to predict the ionization potential depression (IPD) in multicomponent dense plasmas, which is achieved by considering the self-energy of ions and electrons within the quantum statistical theory. In order to explicitly account for the exchange-correlation effect of electrons, we incorporate the effective static approximation of local field correction (LFC) within our IPD framework through the connection of dynamical structure factor. The effective static approximation poses an accurate description for the asymptotic large wave number behavior with the recently developed machine learning representation of static LFC induced from the path-integral Monte Carlo data. Our calculation shows that the introduction of static LFC through dynamical structure factor brings a nontrivial influence on IPD at warm/hot dense matter conditions. The correlation effect within static LFC could provide up to 20% correction to free-electron contribution of IPD in the strong coupling and degeneracy regime. Furthermore, a new screening factor is obtained from the density distribution of free electrons calculated within the average-atom model, with which excellent agreements are observed with other methods and experiments at warm/hot dense matter conditions.

Enhancement of electron-impact ionization induced by warm dense environments

Studies have shown significant discrepancies between the recent experiment [Berg et al., Phys. Rev. Lett. 120, 055002 (2018)PRLTAO10.1103/PhysRevLett.120.055002] and current theoretical calculations on the electron-impact ionization cross section of ions in warm dense magnesium. Here, we present a systematic study the effects of the ionic correlations and free-electron screening on the electron-impact ionization of ions in warm dense matter. The ionic correlation and the free-electron screening effects yield additional Hermitian terms to the calculation of the ionic central-force-field potential, which significantly change the electronic structure compared with that of the isolated ion. In calculating the electron-impact ionization, we describe the impact and ionized electrons using a damped-distorted wave function, which considers the momentum relaxation of free electrons due to collisions with other free electrons and ions. We reproduce the electron-impact ionization process for Mg^{7+} in the solid-density plasma and increase the ionization cross section by one order of magnitude compared with that of the isolated ion, which excellently agrees with the experimental result of Berg et al.

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

We present a free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin-Kohn-Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo-Greenwood optic data with a single atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo-Greenwood calculations with AIMD snapshots. Kubo-Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomenon of redshift-to-blueshift in K-L (1s→2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s→2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

Transient ionization potential depression in nonthermal dense plasmas at high x-ray intensity

The advent of x-ray free-electron lasers (XFELs), which provide intense ultrashort x-ray pulses, has brought a new way of creating and analyzing hot and warm dense plasmas in the laboratory. Because of the ultrashort pulse duration, the XFEL-produced plasma will be out of equilibrium at the beginning, and even the electronic subsystem may not reach thermal equilibrium while interacting with a femtosecond timescale pulse. In the dense plasma, the ionization potential depression (IPD) induced by the plasma environment plays a crucial role for understanding and modeling microscopic dynamical processes. However, all theoretical approaches for IPD have been based on local thermal equilibrium (LTE), and it has been controversial to use LTE IPD models for the nonthermal situation. In this work, we propose a non-LTE (NLTE) approach to calculate the IPD effect by combining a quantum-mechanical electronic-structure calculation and a classical molecular dynamics simulation. This hybrid approach enables us to investigate the time evolution of ionization potentials and IPDs during and after the interaction with XFEL pulses, without the limitation of the LTE assumption. In our NLTE approach, the transient IPD values are presented as distributions evolving with time, which cannot be captured by conventional LTE-based models. The time-integrated ionization potential values are in good agreement with benchmark experimental data on solid-density aluminum plasma and other theoretical predictions based on LTE. The present work is promising to provide critical insights into nonequilibrium dynamics of dense plasma formation and thermalization induced by XFEL pulses.

Energy Flow in Thin Shell Implosions and Explosions

Energy flow and balance in convergent systems beyond petapascal energy densities controls the fate of late-stage stars and the potential for controlling thermonuclear inertial fusion ignition. Time-resolved x-ray self-emission imaging combined with a Bayesian inference analysis is used to describe the energy flow and the potential information stored in the rebounding spherical shock at 0.22 PPa (2.2 Gbar or billions of atmospheres pressure). This analysis, together with a simple mechanical model, describes the trajectory of the shell and the time history of the pressure at the fuel-shell interface, ablation pressure, and energy partitioning including kinetic energy of the shell and internal energy of the fuel. The techniques used here provide a fully self-consistent uncertainty analysis of integrated implosion data, a thermodynamic-path independent measurement of pressure in the petapascal range, and can be used to deduce the energy flow in a wide variety of implosion systems to petapascal energy densities.

Ion-acoustic dipolar vortex in degenerate magnetoplasma with ions/electrons thermal corrections

Nonlinear evolution and propagation characteristics of an ion-acoustic (IA) dipolar vortex are examined in a magnetoplasma that comprises partially degenerate electrons and dynamic ions. The adiabatic ions reduce the quantum magnetohydrodynamic equations to a modified momentum equation. The latter admits a new solution and leads to an evolution equation for the description of a coherent IA dipolar vortex. Numerical analysis reveals that variations in the relevant plasma parameters alter the evolution condition for a stable vortex. The electron exchange-correlation enhances excitation of the vortex as it extends the parametric regime associated with the stable vortex solution. On the contrary, the ion temperature lessens the nonlinear evolution of the vortex. Importantly, a degree of enhancement in the exchange-correlation potential leads to the wave dispersion and-in return-widens the vortex potential spatially. The novel nature of the derived results stresses upon the comprehension of the ordered structures in the compact stars, the magnetic mirrors, the ionosphere, the pinch devices, etc., where thermal corrections significantly impact the waves dynamics.

Interspecies radiative transition in warm and superdense plasma mixtures

Superdense plasmas widely exist in planetary interiors and astrophysical objects such as brown-dwarf cores and white dwarfs. How atoms behave under such extreme-density conditions is not yet well understood, even in single-species plasmas. Here, we apply thermal density functional theory to investigate the radiation spectra of superdense iron–zinc plasma mixtures at mass densities of ρ = 250 to 2000 g cm⁻³ and temperatures of kT = 50 to 100 eV, accessible by double-shell–target implosions. Our ab initio calculations reveal two extreme atomic-physics phenomena—firstly, an interspecies radiative transition; and, secondly, the breaking down of the dipole-selection rule for radiative transitions in isolated atoms. Our first-principles calculations predict that for superdense plasma mixtures, both interatomic radiative transitions and dipole-forbidden transitions can become comparable to the normal intra-atomic Kα-emission signal. These physics phenomena were not previously considered in detail for extreme high-density plasma mixtures at super-high energy densities.

Matter at extremely high density and pressure behaves differently than at ambient conditions. Here the authors use first-principles calculations to show the existence of interspecies radiative and dipole-forbidden transitions in warm and superdense plasma mixture of iron and zinc.

Ionization potential depression and Pauli blocking in degenerate plasmas at extreme densities

New facilities explore warm dense matter (WDM) at conditions with extreme densities (exceeding ten times condensed matter densities) so that electrons are degenerate even at temperatures of 10-100 eV. Whereas in the nondegenerate region correlation effects such as Debye screening are relevant for the ionization potential depression (IPD), new effects have to be considered in degenerate plasmas. In addition to the Fock shift of the self-energies, the bound-state Pauli blocking becomes important with increasing density. Standard approaches to IPD such as Stewart-Pyatt and widely used opacity tables (e.g., OPAL) do not contain Pauli blocking effects for bound states. The consideration of degeneracy effects leads to a reduction of the ionization potential and to a higher degree of ionization. As an example, we present calculations for the ionization degree of carbon plasmas at T = 100 eV and extreme densities up to 40 g/cm^{3}, which are relevant to experiments that are currently scheduled at the National Ignition Facility.

Nuclear excitation by electron capture in optical-laser-generated plasmas

The process of nuclear excitation by electron capture in plasma environments generated by the interaction of ultrastrong optical lasers with solid-state samples is investigated theoretically. With the help of a plasma model, we perform a comprehensive study of the optimal parameters for the most efficient nuclear excitation and determine the corresponding laser setup requirements. We discern between the low-density plasma regime, modeled by scaling laws, and the high-density regime, for which we perform particle-in-cell calculations. As a nuclear transition case study we consider the 4.85-keV nuclear excitation starting from the long-lived ^{93m}Mo isomer. Our results show that the optimal plasma and laser parameters are sensitive to the chosen observable and that measurable rates of nuclear excitation and isomer depletion of ^{93m}Mo should be already achievable at laser facilities existing today.

Path integral Monte Carlo simulations of warm dense aluminum

We perform first-principles path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of aluminum. Our equation of state (EOS) simulations cover a wide density-temperature range of 0.1-32.4gcm^{-3} and 10^{4}-10^{8} K. Since PIMC and DFT-MD accurately treat effects of the atomic shell structure, we find two compression maxima along the principal Hugoniot curve attributed to K-shell and L-shell ionization. The results provide a benchmark for widely used EOS tables, such as SESAME, QEOS, and models based on Thomas-Fermi and average-atom techniques. A subsequent multishock analysis provides a quantitative assessment for how much heating occurs relative to an isentrope in multishock experiments. Finally, we compute heat capacity, pair-correlation functions, the electronic density of states, and 〈Z〉 to reveal the evolution of the plasma structure and ionization behavior.

Validating Continuum Lowering Models via Multi-Wavelength Measurements of Integrated X-ray Emission

X-ray emission spectroscopy is a well-established technique used to study continuum lowering in dense plasmas. It relies on accurate atomic physics models to robustly reproduce high-resolution emission spectra, and depends on our ability to identify spectroscopic signatures such as emission lines or ionization edges of individual charge states within the plasma. Here we describe a method that forgoes these requirements, enabling the validation of different continuum lowering models based solely on the total intensity of plasma emission in systems driven by narrow-bandwidth x-ray pulses across a range of wavelengths. The method is tested on published Al spectroscopy data and applied to the new case of solid-density partially-ionized Fe plasmas, where extracting ionization edges directly is precluded by the significant overlap of emission from a wide range of charge states.

Tailoring Laser-Generated Plasmas for Efficient Nuclear Excitation by Electron Capture

The optimal parameters for nuclear excitation by electron capture in plasma environments generated by the interaction of ultrastrong optical lasers with solid matter are investigated theoretically. As a case study we consider a 4.85 keV nuclear transition starting from the long-lived ^{93m}Mo isomer that can lead to the release of the stored 2.4 MeV excitation energy. We find that due to the complex plasma dynamics, the nuclear excitation rate and the actual number of excited nuclei do not reach their maximum at the same laser parameters. The nuclear excitation achievable with a high-power optical laser is up to twelve and up to six orders of magnitude larger than the values predicted for direct resonant and secondary plasma-mediated excitation at the x-ray free electron laser, respectively. Our results show that the experimental observation of the nuclear excitation of ^{93m}Mo and the subsequent release of stored energy should be possible at laser facilities available today.