Density-functional calculations of transport properties in the nondegenerate limit and the role of electron-electron scattering

Average-atom approach for transport properties of shocked argon in the presence of a magnetic field

We present electron transport calculations of shocked argon based on an average-atom modeling of the plasma and compare them with measurements, involving both incident and reflected shock waves. Since the corresponding experiments are subject to a 5 T magnetic field, the impact of the latter on the Rankine-Hugoniot equations is taken into account, starting from the magnetoresistive hydrodynamics, and the resistivity tensor is deduced from the Boltzmann equation. The resistivity tensor yields the electrical and Hall resistivities. Our average-atom code Paradisio provides the quantities required for the calculation of electrical resistivity within the Ziman-Evans formalism, as well as for the Hall resistivity. We obtain good agreement between calculated conductivities and experimental values, both for the incident and reflected shocks. Our values of the Hall constant are compared to experimental values derived from Hall voltage measurements, as well as to theoretical ones based on the quantum statistical linear-relaxation-time approach.

Design of the first fusion experiment to achieve target energy gain G>1

In this work we present the design of the first controlled fusion laboratory experiment to reach target gain G>1 N221204 (5 December 2022) [Phys. Rev. Lett. 132, 065102 (2024)10.1103/PhysRevLett.132.065102], performed at the National Ignition Facility, where the fusion energy produced (3.15 MJ) exceeded the amount of laser energy required to drive the target (2.05 MJ). Following the demonstration of ignition according to the Lawson criterion N210808, experiments were impacted by nonideal experimental fielding conditions, such as increased (known) target defects that seeded hydrodynamic instabilities or unintentional low-mode asymmetries from nonuniformities in the target or laser delivery, which led to reduced fusion yields less than 1 MJ. This Letter details design changes, including using an extended higher-energy laser pulse to drive a thicker high-density carbon (also known as diamond) capsule, that led to increased fusion energy output compared to N210808 as well as improved robustness for achieving high fusion energies (greater than 1 MJ) in the presence of significant low-mode asymmetries. For this design, the burnup fraction of the deuterium and tritium (DT) fuel was increased (approximately 4% fuel burnup and a target gain of approximately 1.5 compared to approximately 2% fuel burnup and target gain approximately 0.7 for N210808) as a result of increased total (DT plus capsule) areal density at maximum compression compared to N210808. Radiation-hydrodynamic simulations of this design predicted achieving target gain greater than 1 and also the magnitude of increase in fusion energy produced compared to N210808. The plasma conditions and hotspot power balance (fusion power produced vs input power and power losses) using these simulations are presented. Since the drafting of this manuscript, the results of this paper have been replicated and exceeded (N230729) in this design, together with a higher-quality diamond capsule, setting a new record of approximately 3.88MJ of fusion energy and fusion energy target gain of approximately 1.9.

Analysing the dynamic structure of warm dense matter in the imaginary-time domain: theoretical models and simulations

Rigorous diagnostics of experiments with warm dense matter are notoriously difficult. A key method is X-ray Thomson scattering (XRTS), but the interpretation of XRTS measurements is usually based on theoretical models that entail various approximations. Recently, Dornheim et al. [Nat. Commun. 13, 7911 (2022)] introduced a new framework for temperature diagnostics of XRTS experiments that is based on imaginary-time correlation functions. On the one hand, switching from the frequency to the imaginary-time domain gives one direct access to a number of physical properties, which facilitates the extraction of the temperature of arbitrarily complex materials without relying on any models or approximations. On the other hand, the bulk of theoretical work in dynamic quantum many-body theory is devoted to the frequency domain, and, to the best of our knowledge, the manifestation of physics properties within the imaginary-time density-density correlation function (ITCF) remains poorly understood. In the present work, we aim to fill this gap by introducing a simple, semi-analytical model for the imaginary-time dependence of two-body correlations within the framework of imaginary-time path integrals. As a practical example, we compare our new model to extensive ab initio path integral Monte Carlo results for the ITCF of a uniform electron gas, and find excellent agreement over a broad range of wavenumbers, densities and temperatures. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Lawson Criterion for Ignition Exceeded in an Inertial Fusion Experiment

For more than half a century, researchers around the world have been engaged in attempts to achieve fusion ignition as a proof of principle of various fusion concepts. Following the Lawson criterion, an ignited plasma is one where the fusion heating power is high enough to overcome all the physical processes that cool the fusion plasma, creating a positive thermodynamic feedback loop with rapidly increasing temperature. In inertially confined fusion, ignition is a state where the fusion plasma can begin "burn propagation" into surrounding cold fuel, enabling the possibility of high energy gain. While "scientific breakeven" (i.e., unity target gain) has not yet been achieved (here target gain is 0.72, 1.37 MJ of fusion for 1.92 MJ of laser energy), this Letter reports the first controlled fusion experiment, using laser indirect drive, on the National Ignition Facility to produce capsule gain (here 5.8) and reach ignition by nine different formulations of the Lawson criterion.

Design of an inertial fusion experiment exceeding the Lawson criterion for ignition

We present the design of the first igniting fusion plasma in the laboratory by Lawson's criterion that produced 1.37 MJ of fusion energy, Hybrid-E experiment N210808 (August 8, 2021) [Phys. Rev. Lett. 129, 075001 (2022)10.1103/PhysRevLett.129.075001]. This design uses the indirect drive inertial confinement fusion approach to heat and compress a central "hot spot" of deuterium-tritium (DT) fuel using a surrounding dense DT fuel piston. Ignition occurs when the heating from absorption of α particles created in the fusion process overcomes the loss mechanisms in the system for a duration of time. This letter describes key design changes which enabled a ∼3-6× increase in an ignition figure of merit (generalized Lawson criterion) [Phys. Plasmas 28, 022704 (2021)1070-664X10.1063/5.0035583, Phys. Plasmas 25, 122704 (2018)1070-664X10.1063/1.5049595]) and an eightfold increase in fusion energy output compared to predecessor experiments. We present simulations of the hot-spot conditions for experiment N210808 that show fundamentally different behavior compared to predecessor experiments and simulated metrics that are consistent with N210808 reaching for the first time in the laboratory "ignition."

Electronic transport coefficients from density functional theory across the plasma plane

We investigate the thermopower and Lorenz number of hydrogen with Kohn-Sham density functional theory (DFT) across the plasma plane toward the near-classical limit, i.e., weakly degenerate and weakly coupled states. Our results are in concordance with certain limiting values for the Lorentz plasma, a model system which only considers electron-ion scattering. Thereby, we clearly show that the widely used method of calculating transport properties via the Kubo-Greenwood (KG) formalism does not capture electron-electron scattering processes. Our discussion also addresses the inadequateness of assuming a Drude-like frequency behavior for the conductivity of nondegenerate plasmas by revisiting the relaxation time approximation within kinetic theory.

Virial expansion of the electrical conductivity of hydrogen plasmas

The low-density limit of the electrical conductivity σ(n,T) of hydrogen as the simplest ionic plasma is presented as a function of the temperature T and mass density n in the form of a virial expansion of the resistivity. Quantum statistical methods yield exact values for the lowest virial coefficients which serve as a benchmark for analytical approaches to the electrical conductivity as well as for numerical results obtained from density functional theory-based molecular dynamics simulations (DFT-MD) or path-integral Monte Carlo simulations. While these simulations are well suited to calculate σ(n,T) in a wide range of density and temperature, in particular, for the warm dense matter region, they become computationally expensive in the low-density limit, and virial expansions can be utilized to balance this drawback. We present new results of DFT-MD simulations in that regime and discuss the account of electron-electron collisions by comparison with the virial expansion.

Kinetic model for electron-ion transport in warm dense matter

We present a model for electron-ion transport in warm dense matter that incorporates Coulomb coupling effects into the quantum Boltzmann equation of Uehling and Uhlenbeck through the use of a statistical potential of mean force. Although the model presented here can be derived rigorously in the classical limit [S. D. Baalrud and J. Daligault, Phys. Plasmas 26, 082106 (2019)PHPAEN1070-664X10.1063/1.5095655], its quantum generalization is complicated by the uncertainty principle. Here we apply an existing model for the potential of mean force based on the quantum Ornstein-Zernike equation coupled with an average-atom model [C. E. Starrett, High Energy Density Phys. 25, 8 (2017)1574-181810.1016/j.hedp.2017.09.003]. This potential contains correlations due to both Coulomb coupling and exchange, and the collision kernel of the kinetic theory enforces Pauli blocking while allowing for electron diffraction and large-angle collisions. We use the Uehling-Uhlenbeck equation to predict the momentum and temperature relaxation times and electrical conductivity of solid density aluminum plasma based on electron-ion collisions. We present results for density and temperature conditions that span the transition from classical weakly-coupled plasma to degenerate moderately-coupled plasma. Our findings agree well with recent quantum molecular dynamics simulations.

Ionization and transport in partially ionized multicomponent plasmas: Application to atmospheres of hot Jupiters

We study ionization and transport processes in partially ionized multicomponent plasmas. The plasma composition is calculated via a system of coupled mass-action laws. The electronic transport properties are determined by the electron-ion and electron-neutral transport cross sections. The influence of electron-electron scattering is considered via a correction factor to the electron-ion contribution. Based on these data, the electrical and thermal conductivities as well as the Lorenz number are calculated. For the thermal conductivity, we consider also the contributions of the translational motion of neutral particles and of the dissociation, ionization, and recombination reactions. We apply our approach to a partially ionized plasma composed of hydrogen, helium, and a small fraction of metals (Li, Na, Ca, Fe, K, Rb, and Cs) as typical for atmospheres of hot Jupiters. We present results for the plasma composition and the transport properties as a function of density and temperature and then along typical P-T profiles for the outer part of the hot Jupiter HD 209458b. The electrical conductivity profile allows revising the Ohmic heating power related to the fierce winds in the planet's atmosphere. We show that the higher temperatures suggested by recent interior models could boost the conductivity and thus the Ohmic heating power to values large enough to explain the observed inflation of HD 209458b.

Ultrafast multi-cycle terahertz measurements of the electrical conductivity in strongly excited solids

Key insights in materials at extreme temperatures and pressures can be gained by accurate measurements that determine the electrical conductivity. Free-electron laser pulses can ionize and excite matter out of equilibrium on femtosecond time scales, modifying the electronic and ionic structures and enhancing electronic scattering properties. The transient evolution of the conductivity manifests the energy coupling from high temperature electrons to low temperature ions. Here we combine accelerator-based, high-brightness multi-cycle terahertz radiation with a single-shot electro-optic sampling technique to probe the evolution of DC electrical conductivity using terahertz transmission measurements on sub-picosecond time scales with a multi-undulator free electron laser. Our results allow the direct determination of the electron-electron and electron-ion scattering frequencies that are the major contributors of the electrical resistivity.

The electrical conductivity is critical to understand warm dense matter, but the accurate measurement is extremely challenging. Here the authors use multi-cycle THz pulses to measure the conductivity of gold foils strongly heated by free-electron laser, determining the individual contributions of electron-electron and electron-ion scattering.

Finite-size effects in the reconstruction of dynamic properties from ab initio path integral Monte Carlo simulations

We systematically investigate finite-size effects in the dynamic structure factor S(q,ω) of the uniform electron gas obtained via the analytic continuation of ab initio path integral Monte Carlo data for the imaginary-time density-density correlation function F(q,τ). Using the recent scheme by Dornheim et al. [Phys. Rev. Lett. 121, 255001 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.255001], we find that the reconstructed spectra are not afflicted with any finite-size effects for as few as N=14 electrons both at warm dense matter (WDM) conditions and at the margins of the strongly correlated electron liquid regime. Our results further corroborate the high quality of our current description of the dynamic density response of correlated electrons, which is of high importance for many applications in WDM theory and beyond.

Transport properties of warm dense neon and krypton at high pressures

The transport properties of warm dense neon (Ne) and krypton (Kr) are studied by combining self-consistent fluid variational theory (SFVT) with linear response theory (LRT). The components are determined using the SFVT, and the transport parameters, including the electrical conductivity, thermal conductivity, and thermopower, are calculated with the LRT. The relevant scattering mechanisms, including electron-ion, electron-electron, and electron-atom scatterings, are taken into account. An effective potential model in combination with the Muffin-tin model is introduced to further improve the description for electron-atom scattering, which not only includes static, exchange, and polarization interactions but also considers the plasma environmental effects. It is found that for electron-atom scattering, the influence of the plasma density is significant at lower scattering energies but the effects are different for electron-Ne and electron-Kr scattering. For electron-Kr scattering, a plasma density-dependent Ramsauer-Townsend minimum is observed. The obtained transport parameters are compared with the available experiments and other simulations. The plasma phase transition of warm dense Kr is revisited from multiple perspectives based on the numerical simulation results for the electrical conductivity and thermopower. These observations may help one to better understand the transport properties of warm dense noble gases and are an important guide for future experimental designs and theoretical developments.

Probing the Electronic Structure of Warm Dense Nickel via Resonant Inelastic X-Ray Scattering

The development of bright free-electron lasers (FEL) has revolutionized our ability to create and study matter in the high-energy-density (HED) regime. Current diagnostic techniques have been successful in yielding information on fundamental thermodynamic plasma properties, but provide only limited or indirect information on the detailed quantum structure of these systems, and on how it is affected by ionization dynamics. Here we show how the valence electronic structure of solid-density nickel, heated to temperatures of around 10 of eV on femtosecond timescales, can be probed by single-shot resonant inelastic x-ray scattering (RIXS) at the Linac Coherent Light Source FEL. The RIXS spectrum provides a wealth of information on the HED system that goes well beyond what can be extracted from x-ray absorption or emission spectroscopy alone, and is particularly well suited to time-resolved studies of electronic-structure dynamics.

Particle-in-cell simulation method for macroscopic degenerate plasmas

Nowadays hydrodynamic equations coupled with external equation of states provided by quantum mechanical calculations is a widely used approach for simulations of macroscopic degenerate plasmas. Although such an approach is proven to be efficient and shows many good features, especially for large scale simulations, it encounters intrinsic challenges when involving kinetic effects. As a complement, here we have invented a fully kinetic numerical approach for macroscopic degenerate plasmas. This approach is based on first principle Boltzmann-Uhling-Uhlenbeck equations coupled with Maxwell's equation, and is eventually achieved via an existing particle-in-cell simulation code named LAPINS. In this approach, degenerate particles obey Fermi-Dirac statistics and nondegenerate particles follow the typical Maxwell-Boltzmann statistics. The equation of motion of both degenerate and nondegenerate particles are governed by long range collective electromagnetic fields and close particle-particle collisions. Especially, Boltzmann-Uhling-Uhlenbeck collisions ensure that evolution of degenerate particles is enforced by the Pauli exclusion principle. The code is applied to several benchmark simulations, including electronic conductivity for aluminium with varying temperatures from 2 eV to 50 eV, thermalization of alpha particles in a cold fuel shell in inertial confinement fusion, and rapid heating of solid sample by short and intense laser pulses.

Model of electron transport in dense plasmas spanning temperature regimes

We present a new model of electron transport in warm and hot dense plasmas which combines the quantum Landau-Fokker-Planck equation with the concept of mean-force scattering. We obtain electrical and thermal conductivities across several orders of magnitude in temperature, from warm dense matter conditions to hot, nondegenerate plasma conditions, including the challenging crossover regime between the two. The small-angle approximation characteristic of Fokker-Planck collision theories is mitigated to good effect by the construction of accurate effective Coulomb logarithms based on mean-force scattering, which allows the theory to remain accurate even at low temperatures, as compared with high-fidelity quantum simulation results. Electron-electron collisions are treated on equal footing as electron-ion collisions. Their accurate treatment is found to be essential for hydrogen, and is expected to be important to other low-Z elements. We find that electron-electron scattering remains influential to the value of the thermal conductivity down to temperatures somewhat below the Fermi energy. The accuracy of the theory seems to falter only for the behavior of the thermal conductivity at very low temperatures due to a subtle interplay between the Pauli exclusion principle and the small-angle approximation as they pertain to electron-electron scattering. Even there, the model is in fair agreement with ab initio simulations.

Correlations between conduction electrons in dense plasmas

Most treatments of electron-electron correlations in dense plasmas either ignore them entirely (random phase approximation) or neglect the role of ions (jellium approximation). In this work, we go beyond both these approximations to derive a formula for the electron-electron static structure factor which properly accounts for the contributions of both ionic structure and quantum-mechanical dynamic response in the electrons. The result can be viewed as a natural extension of the quantum Ornstein-Zernike theory of ionic and electronic correlations, and it is suitable for dense plasmas in which the ions are classical and the conduction electrons are quantum-mechanical. The corresponding electron-electron pair distribution functions are compared with the results of path integral Monte Carlo simulations, showing good agreement whenever no strong electron resonance states are present. We construct approximate potentials of mean force which describe the effective screened interaction between electrons. Significant deviations from Debye-Hückel screening are present at temperatures and densities relevant to high-energy density experiments involving warm and hot dense plasmas. The presence of correlations between conduction electrons is likely to influence the electron-electron contribution to the electrical and thermal conductivity. It is expected that excitation processes involving the conduction electrons (e.g., free-free absorption) will also be affected.

Heat-release equation of state and thermal conductivity of warm dense carbon by proton differential heating

Warm dense carbon is generated at 0.3-2.0 g/cc and 1-7 eV by proton heating. The release equation of state (EOS) after heating and thermal conductivity of warm dense carbon are studied experimentally in this regime using a Au/C dual-layer target to initiate a temperature gradient and two picosecond time-resolved diagnostics to probe the surface expansion and heat flow. Comparison between the data and simulations using various EOSs and thermal conductivity models is quantified with a statistical χ^{2} analysis. Out of seven EOS tables and five thermal conductivity models, only L9061 with the Lee-More model provides a probability above 50% to match all data.

Comment on "Isochoric, isobaric, and ultrafast conductivities of aluminum, lithium, and carbon in the warm dense matter regime"

Dharma-wardana et al. [M. W. C. Dharma-wardana et al., Phys. Rev. E 96, 053206 (2017)2470-004510.1103/PhysRevE.96.053206] recently calculated dynamic electrical conductivities for warm dense matter as well as for nonequilibrium two-temperature states termed "ultrafast matter" (UFM) [M. W. C. Dharma-wardana, Phys. Rev. E 93, 063205 (2016)2470-004510.1103/PhysRevE.93.063205]. In this Comment we present two evident reasons why these UFM calculations are neither suited to calculate dynamic conductivities nor x-ray Thomson scattering spectra in isochorically heated warm dense aluminum. First, the ion-ion structure factor, a major input into the conductivity and scattering spectra calculations, deviates strongly from that of isochorically heated aluminum. Second, the dynamic conductivity does not show a non-Drude behavior which is an essential prerequisite for a correct description of the absorption behavior in aluminum. Additionally, we clarify misinterpretations by Dharma-wardana et al. concerning the conductivity measurements of Gathers [G. R. Gathers, Int. J. Thermophys. 4, 209 (1983)IJTHDY0195-928X10.1007/BF00502353].

Ab initio Path Integral Monte Carlo Results for the Dynamic Structure Factor of Correlated Electrons: From the Electron Liquid to Warm Dense Matter

The accurate description of electrons at extreme density and temperature is of paramount importance for, e.g., the understanding of astrophysical objects and inertial confinement fusion. In this context, the dynamic structure factor S(q,ω) constitutes a key quantity as it is directly measured in x-ray Thomson scattering experiments and governs transport properties like the dynamic conductivity. In this work, we present the first ab initio results for S(q,ω) by carrying out extensive path integral Monte Carlo simulations and developing a new method for the required analytic continuation, which is based on the stochastic sampling of the dynamic local field correction G(q,ω). In addition, we find that the so-called static approximation constitutes a promising opportunity to obtain high-quality data for S(q,ω) over substantial parts of the warm dense matter regime.

Thermal Conductivity and Electrical Resistivity of Solid Iron at Earth's Core Conditions from First Principles

We compute the thermal conductivity and electrical resistivity of solid hcp Fe to pressures and temperatures of Earth's core. We find significant contributions from electron-electron scattering, usually neglected at high temperatures in transition metals. Our calculations show a quasilinear relation between the electrical resistivity and temperature for hcp Fe at extreme high pressures. We obtain thermal and electrical conductivities that are consistent with experiments considering reasonable error. The predicted thermal conductivity is reduced from previous estimates that neglect electron-electron scattering. Our estimated thermal conductivity for the outer core is 77±10  W m^{-1} K^{-1} and is consistent with a geodynamo driven by thermal convection.

Dynamic electron-ion collisions and nuclear quantum effects in quantum simulation of warm dense matter

The structural, thermodynamic and transport properties of warm dense matter (WDM) are crucial to the fields of astrophysics and planet science, as well as inertial confinement fusion. WDM refers to the states of matter in a regime of temperature and density between cold condensed matter and hot ideal plasmas, where the density is from near-solid up to ten times solid density, and the temperature between 0.1 and 100 eV. In the WDM regime, matter exhibits moderately or strongly coupled, partially degenerate properties. Therefore, the methods used to deal with condensed matter and isolated atoms need to be properly validated for WDM. It is therefore a big challenge to understand WDM within a unified theoretical description with reliable accuracy. Here, we review the progress in the theoretical study of WDM with state-of-the-art simulations, i.e. quantum Langevin molecular dynamics and first principles path integral molecular dynamics. The related applications for WDM are also included.

Thermal conductivity measurements of proton-heated warm dense aluminum

Thermal conductivity is one of the most crucial physical properties of matter when it comes to understanding heat transport, hydrodynamic evolution, and energy balance in systems ranging from astrophysical objects to fusion plasmas. In the warm dense matter regime, experimental data are very scarce so that many theoretical models remain untested. Here we present the first thermal conductivity measurements of aluminum at 0.5–2.7 g/cc and 2–10 eV, using a recently developed platform of differential heating. A temperature gradient is induced in a Au/Al dual-layer target by proton heating, and subsequent heat flow from the hotter Au to the Al rear surface is detected by two simultaneous time-resolved diagnostics. A systematic data set allows for constraining both thermal conductivity and equation-of-state models. Simulations using Purgatorio model or Sesame S27314 for Al thermal conductivity and LEOS for Au/Al release equation-of-state show good agreement with data after 15 ps. Discrepancy still exists at early time 0–15 ps, likely due to non-equilibrium conditions.