X-ray scattering measurements on imploding CH spheres at the National Ignition Facility

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}.

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.

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⁶.

Nonideal effect of free electrons on ionization equilibrium and radiative property in dense plasmas

The thermodynamic as well as optical properties of strongly coupled plasmas depend crucially on the average degree of ionization and the ionic state composition, which, however, cannot be determined by using the normal Saha equation usually used for the ideal plasmas. Hence, an adequate treatment of the ionization balance and the charge state distribution of strongly coupled plasmas is still a challenge for theory due to the interactions between the electrons and ions and among the electrons themselves. Based on a local density temperature-dependent ion-sphere model, the Saha equation approach is extended to the regime of strongly coupled plasmas by taking into account the free-electron-ion interaction, the free-free-electron interaction, the nonuniform free-electron space distribution, and the free-electron quantum partial degeneracy. All the quantities, including the bound orbitals with ionization potential depression, free-electron distribution, and bound and free-electron partition function contributions, are calculated self-consistently in the theoretical formalism. This study shows that the ionization equilibrium is evidently modified by considering the above nonideal characteristics of the free electrons. Our theoretical formalism is validated by the explanation of a recent experimental measurement of the opacity of dense hydrocarbon.

Accurate temperature diagnostics for matter under extreme conditions

The experimental investigation of matter under extreme densities and temperatures, as in astrophysical objects and nuclear fusion applications, constitutes one of the most active frontiers at the interface of material science, plasma physics, and engineering. The central obstacle is given by the rigorous interpretation of the experimental results, as even the diagnosis of basic parameters like the temperature T is rendered difficult at these extreme conditions. Here, we present a simple, approximation-free method to extract the temperature of arbitrarily complex materials in thermal equilibrium from X-ray Thomson scattering experiments, without the need for any simulations or an explicit deconvolution. Our paradigm can be readily implemented at modern facilities and corresponding experiments will have a profound impact on our understanding of warm dense matter and beyond, and open up a variety of appealing possibilities in the context of thermonuclear fusion, laboratory astrophysics, and related disciplines.

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.

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.

High Pressure Hydrocarbons Revisited: From van der Waals Compounds to Diamond

Methane and other hydrocarbons are major components of the mantle regions of icy planets. Several recent computational studies have investigated the high-pressure behaviour of specific hydrocarbons. To develop a global picture of hydrocarbon stability, to identify relevant decomposition reactions, and probe eventual formation of diamond, a complete study of all hydrocarbons is needed. Using density functional theory calculations we survey here all known C-H crystal structures augmented by targeted crystal structure searches to build hydrocarbon phase diagrams in the ground state and at elevated temperatures. We find that an updated pressure-temperature phase diagram for methane is dominated at intermediate pressures by CH 4 :H 2 van der Waals inclusion compounds. We discuss the P-T phase diagram for CH and CH 2 (i.e., polystyrene and polyethylene) to illustrate that diamond formation conditions are strongly composition dependent. Finally, crystal structure searches uncover a new CH 4 (H 2 ) 2 van der Waals compound, the most hydrogen-rich hydrocarbon, stable between 170 and 220 GPa.

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.

Using time-resolved penumbral imaging to measure low hot spot x-ray emission signals from capsule implosions at the National Ignition Facility

We have developed and fielded a new x-ray pinhole-imaging snout that exploits time-resolved penumbral imaging of low-emission hot spots in capsule implosion experiments at the National Ignition Facility. We report results for a series of indirectly driven Be capsule implosions that aim at measuring x-ray Thomson scattering (XRTS) spectra at extreme density conditions near stagnation. In these implosions, x-ray emission at stagnation is reduced by 100-1000× compared to standard inertial confinement fusion (ICF) implosions to mitigate undesired continuum background in the XRTS spectra. Our snout design not only enables measurements of peak x-ray emission times, t _o , where standard ICF diagnostics would not record any signal, but also allows for inference of hot spot shapes. Measurement of t _o is crucial to account for shot-to-shot variations in implosion velocity and therefore to benchmark the achieved plasma conditions between shots and against radiation hydrodynamic simulations. Additionally, we used differential filtering to infer a hot spot temperature of 520 ± 80 eV, which is in good agreement with predictions from radiation hydrodynamic simulations. We find that, despite fluctuations of the x-ray flash intensity of up to 5×, the emission time history is similar from shot to shot and slightly asymmetric with respect to peak x-ray emission.

Developing a long-duration Zn K- α source for x-ray scattering experiments

We are developing a long-duration K-α x-ray source at the Omega laser facility. Such sources are important for x-ray scattering measurements at small scattering angles where high spectral resolution is required. To date, He-α x-ray sources are the most common probes in scattering experiments, using ns-class lasers to heat foils to keV temperatures, resulting in K-shell emission from He-like charge states. The He-α spectrum can be broadened by emission from multiple charge states and lines (e.g., He-like, Li-like, Be-like). Here, we combine the long duration of He-α sources with the narrow spectral bandwidth of cold K-α emission. A Ge foil is irradiated by the Omega laser, producing principally Ge He-α emission, which pumps Zn K-α emission at 8.6 keV from a nearby Zn layer. Using this technique, we demonstrate a long-duration Zn K-α source suitable for scattering measurements. Our experimental results show a 60% reduction in spectral bandwidth compared to a standard Zn He-α source, significantly improving the measurement precision of scattering experiments with small inelastic shifts.

A platform for x-ray Thomson scattering measurements of radiation hydrodynamics experiments on the NIF

We present an experimental design for a radiation hydrodynamics experiment at the National Ignition Facility that measures the electron temperature of a shocked region using the x-ray Thomson scattering technique. Previous National Ignition Facility experiments indicate a reduction in Rayleigh-Taylor instability growth due to high energy fluxes, compared to the shocked energy flux, from radiation and electron heat conduction. In order to better quantify the effects of these energy fluxes, we modified the previous experiment to allow for non-collective x-ray Thomson scattering to measure the electron temperature. Photometric calculations combined with synthetic scattering spectra demonstrate an estimated noise.

Absolute Equation-of-State Measurement for Polystyrene from 25 to 60 Mbar Using a Spherically Converging Shock Wave

We have developed an experimental platform for the National Ignition Facility that uses spherically converging shock waves for absolute equation-of-state (EOS) measurements along the principal Hugoniot. In this Letter, we present one indirect-drive implosion experiment with a polystyrene sample that employs radiographic compression measurements over a range of shock pressures reaching up to 60 Mbar (6 TPa). This significantly exceeds previously published results obtained on the Nova laser [R. Cauble et al., Phys. Rev. Lett. 80, 1248 (1998)PRLTAO0031-900710.1103/PhysRevLett.80.1248] at a strongly improved precision, allowing us to discriminate between different EOS models. We find excellent agreement with Kohn-Sham density-functional-theory-based molecular dynamics simulations.

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.

Ionization-potential depression and dynamical structure factor in dense plasmas

The properties of a bound electron system immersed in a plasma environment are strongly modified by the surrounding plasma. The modification of an essential quantity, the ionization energy, is described by the electronic and ionic self-energies, including dynamical screening within the framework of the quantum statistical theory. Introducing the ionic dynamical structure factor as the indicator for the ionic microfield, we demonstrate that ionic correlations and fluctuations play a critical role in determining the ionization potential depression. This is, in particular, true for mixtures of different ions with large mass and charge asymmetry. The ionization potential depression is calculated for dense aluminum plasmas as well as for a CH plasma and compared to the experimental data and more phenomenological approaches used so far.

First-principles equation of state and shock compression predictions of warm dense hydrocarbons

We use path integral Monte Carlo and density functional molecular dynamics to construct a coherent set of equations of state (EOS) for a series of hydrocarbon materials with various C:H ratios (2:1, 1:1, 2:3, 1:2, and 1:4) over the range of 0.07-22.4gcm^{-3} and 6.7×10^{3}-1.29×10^{8}K. The shock Hugoniot curve derived for each material displays a single compression maximum corresponding to K-shell ionization. For C:H = 1:1, the compression maximum occurs at 4.7-fold of the initial density and we show radiation effects significantly increase the shock compression ratio above 2 Gbar, surpassing relativistic effects. The single-peaked structure of the Hugoniot curves contrasts with previous work on higher-Z plasmas, which exhibit a two-peak structure corresponding to both K- and L-shell ionization. Analysis of the electronic density of states reveals that the change in Hugoniot structure is due to merging of the L-shell eigenstates in carbon, while they remain distinct for higher-Z elements. Finally, we show that the isobaric-isothermal linear mixing rule for carbon and hydrogen EOS is a reasonable approximation with errors better than 1% for stellar-core conditions.

X-ray Thomson scattering measurements from hohlraum-driven spheres on the OMEGA laser

X-ray Thomson scattering (XRTS) is a powerful diagnostic for probing warm and hot dense matter. We present the design and results of the first XRTS experiments with hohlraum-driven CH₂ targets on the OMEGA laser facility at the Laboratory for Laser Energetics in Rochester, NY. X-rays seen directly from the XRTS x-ray source overshadow the elastic scattering signal from the target capsule but can be controlled in future experiments. From the inelastic scattering signal, an average plasma temperature is inferred that is in reasonable agreement with the temperatures predicted by simulations. Knowledge gained in this experiment shows a promising future for further XRTS measurements on indirectly driven OMEGA targets.

Improving a high-efficiency, gated spectrometer for x-ray Thomson scattering experiments at the National Ignition Facility

We are developing x-ray Thomson scattering for applications in implosion experiments at the National Ignition Facility. In particular we have designed and fielded MACS, a high-efficiency, gated x-ray spectrometer at 7.5-10 keV [T. Döppner et al., Rev. Sci. Instrum. 85, 11D617 (2014)]. Here we report on two new Bragg crystals based on Highly Oriented Pyrolytic Graphite (HOPG), a flat crystal and a dual-section cylindrically curved crystal. We have performed in situ calibration measurements using a brass foil target, and we used the flat HOPG crystal to measure Mo K-shell emission at 18 keV in 2nd order diffraction. Such high photon energy line emission will be required to penetrate and probe ultra-high-density plasmas or plasmas of mid-Z elements.

Atomic bound state and scattering properties of effective momentum-dependent potentials

Effective classical dynamics provide a potentially powerful avenue for modeling large-scale dynamical quantum systems. We have examined the accuracy of a Hamiltonian-based approach that employs effective momentum-dependent potentials (MDPs) within a molecular-dynamics framework through studies of atomic ground states, excited states, ionization energies, and scattering properties of continuum states. Working exclusively with the Kirschbaum-Wilets (KW) formulation with empirical MDPs [C. L. Kirschbaum and L. Wilets, Phys. Rev. A 21, 834 (1980)0556-279110.1103/PhysRevA.21.834], optimization leads to very accurate ground-state energies for several elements (e.g., N, F, Ne, Al, S, Ar, and Ca) relative to Hartree-Fock values. The KW MDP parameters obtained are found to be correlated, thereby revealing some degree of transferability in the empirically determined parameters. We have studied excited-state orbits of electron-ion pair to analyze the consequences of the MDP on the classical Coulomb catastrophe. From the optimized ground-state energies, we find that the experimental first- and second-ionization energies are fairly well predicted. Finally, electron-ion scattering was examined by comparing the predicted momentum transfer cross section to a semiclassical phase-shift calculation; optimizing the MDP parameters for the scattering process yielded rather poor results, suggesting a limitation of the use of the KW MDPs for plasmas.

Nuclear-plus-interference-scattering effect on the energy deposition of multi-MeV protons in a dense Be plasma

The nuclear plus interference scattering (NIS) effect on the stopping power of hot dense beryllium (Be) plasma for multi-MeV protons is theoretically investigated by using the generalized Brown-Preston-Singleton (BPS) model, in which a NIS term is taken into account. The analytical formula of the NIS term is detailedly derived. By using this formula, the density and temperature dependence of the NIS effect is numerically studied, and the results show that the NIS effect becomes more and more important with increasing the plasma temperature or density. Different from the cases of protons traveling through the deuterium-tritium plasmas, for a Be plasma, a prominent oscillation valley structure is observed in the NIS term when the proton's energy is close to E_{p}=7MeV. Furthermore, the penetration distance is remarkably reduced when the NIS term is considered.