Search

We present a model of the electron thermal conductivity of a laser-produced plasma. The model, supported by Vlasov-Fokker-Planck simulations, predicts that laser absorption reduces conductivity by forcing electrons out of a Maxwell-Boltzmann equilibrium, which results in the depletion of both low-velocity bulk electrons and high-velocity tail electrons. We show that both the bulk and tail electrons approximately follow super-Gaussian distributions, but with distinct exponents that each depend on the laser intensity and wavelength through the parameter α=Zv_{E}^{2}/v_{T}^{2}. For a value of α=0.5, tail depletion reduces the thermal conductivity to half its zero-intensity value. We present our results as simple analytic fits that can be readily implemented in any radiation-hydrodynamics code or used to correct the local limit of nonlocal conduction models.

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. We utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code lilac to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. Results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

Laser-direct-drive fusion target designs with solid deuterium-tritium (DT) fuel, a high-Z gradient-density pusher shell (GDPS), and a Au-coated foam layer have been investigated through both 1D and 2D radiation-hydrodynamic simulations. Compared with conventional low-Z ablators and DT-push-on-DT targets, these GDPS targets possess certain advantages of being instability-resistant implosions that can be high adiabat (α≥8) and low hot-spot and pusher-shell convergence (CR&#95;{hs}≈22 and CR&#95;{PS}≈17), and have a low implosion velocity (v&#95;{imp}<3×10^{7}cm/s). Using symmetric drive with laser energies of 1.9 to 2.5MJ, 1D lilac simulations of these GDPS implosions can result in neutron yields corresponding to ≳50-MJ energy, even with reduced laser absorption due to the cross-beam energy transfer (CBET) effect. Two-dimensional draco simulations show that these GDPS targets can still ignite and deliver neutron yields from 4 to ∼10MJ even if CBET is present, while traditional DT-push-on-DT targets normally fail due to the CBET-induced reduction of ablation pressure. If CBET is mitigated, these GDPS targets are expected to produce neutron yields of >20MJ at a driven laser energy of ∼2MJ. The key factors behind the robust ignition and moderate energy gain of such GDPS implosions are as follows: (1) The high initial density of the high-Z pusher shell can be placed at a very high adiabat while the DT fuel is maintained at a relatively low-entropy state; therefore, such implosions can still provide enough compression ρR>1g/cm^{2} for sufficient confinement; (2) the high-Z layer significantly reduces heat-conduction loss from the hot spot since thermal conductivity scales as ∼1/Z; and (3) possible radiation trapping may offer an additional advantage for reducing energy loss from such high-Z targets.

Warm dense matter is a material state in the region of parameter space connecting condensed matter to classical plasma physics. In this intermediate regime, we investigate the significance of non-adiabatic electron-ion interactions upon ion dynamics. To disentangle non-adiabatic from adiabatic electron-ion interactions, we compare the ion self-diffusion coefficient from the non-adiabatic electron force field computational model with an adiabatic, classical molecular dynamics simulation. A classical pair potential developed through a force-matching algorithm ensures the only difference between the models is due to the electronic inertia. We implement this new method to characterize non-adiabatic effects on the self-diffusion of warm dense hydrogen over a wide range of temperatures and densities. Ultimately we show that the impact of non-adiabatic effects is negligible for equilibrium ion dynamics in warm dense hydrogen. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

In the dynamic-shell (DS) concept [V. N. Goncharov et al., Novel Hot-Spot Ignition Designs for Inertial Confinement Fusion with Liquid-Deuterium-Tritium Spheres, Phys. Rev. Lett. 125, 065001 (2020).PRLTAO0031-900710.1103/PhysRevLett.125.065001] for laser-driven inertial confinement fusion the deuterium-tritium fuel is initially in the form of a homogeneous liquid inside a wetted-foam spherical shell. This fuel is ignited using a conventional implosion, which is preceded by a initial compression of the fuel followed by its expansion and dynamic formation of a high-density fuel shell with a low-density interior. This Letter reports on a scaled-down, proof-of-principle experiment on the OMEGA laser demonstrating, for the first time, the feasibility of DS formation. A shell is formed by convergent shocks launched by laser pulses at the edge of a plasma sphere, with the plasma itself formed as a result of laser-driven compression and relaxation of a surrogate plastic-foam ball target. Three x-ray diagnostics, namely, 1D spatially resolved self-emission streaked imaging, 2D self-emission framed imaging, and backlighting radiography, have shown good agreement with the predicted evolution of the DS and its stability to low Legendre mode perturbations introduced by laser irradiation and target asymmetries.

Inverse bremsstrahlung absorption was measured based on transmission through a finite-length plasma that was thoroughly characterized using spatially resolved Thomson scattering. Expected absorption was then calculated using the diagnosed plasma conditions while varying the absorption model components. To match data, it is necessary to account for (i) the Langdon effect; (ii) laser-frequency (rather than plasma-frequency) dependence in the Coulomb logarithm, as is typical of bremsstrahlung theories but not transport theories; and (iii) a correction due to ion screening. Radiation-hydrodynamic simulations of inertial confinement fusion implosions have to date used a Coulomb logarithm from the transport literature and no screening correction. We anticipate that updating the model for collisional absorption will substantially revise our understanding of laser-target coupling for such implosions.

In plasmas, electronic states can be well-localized bound states or itinerant free states, or something in between. In self-consistent treatments of plasma electronic structure such as the average-atom model, all states must be accurately resolved in order to achieve a converged numerical solution. This is a challenging numerical and algorithmic problem in large part due to the continuum of free states which is relatively expensive and difficult to resolve accurately. Siegert states are an appealing alternative. They form a complete eigenbasis with a purely discrete spectrum while still being equivalent to a representation in terms of the usual bound states and free states. However, many of their properties are unintuitive, and it is not obvious that they are suitable for self-consistent plasma electronic structure calculations. Here it is demonstrated that Siegert states can be used to accurately solve an average-atom model and offer advantages over the traditional finite-difference approach, including a concrete physical picture of pressure ionization and continuum resonances.

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.

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., (© 2022. The Author(s).)

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.

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.

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.

The pressure and internal energy of an ultracold plasma in a state of quasiequilibrium are evaluated using classical molecular dynamics simulations. Coulomb collapse is avoided by modeling electron-ion interactions using an attractive Coulomb potential with a repulsive core. We present a method to separate the contribution of classical bound states, which form due to recombination, from the contribution of free charges when evaluating these thermodynamic state variables. It is found that the contribution from free charges is independent of the choice of repulsive core length scale when it is sufficiently short-ranged. The partial pressure associated with the free charges is found to closely follow that of the one-component plasma model, reaching negative values at strong coupling, while the total system pressure remains positive. This pseudopotential model is also applied to Debye-Hückel theory to describe the weakly coupled regime.

Self-diffusion and interdiffusion coefficients of binary ionic mixtures are evaluated using the effective potential theory (EPT), and the predictions are compared with the results of molecular dynamics simulations. We find that EPT agrees with molecular dynamics from weak coupling well into the strong-coupling regime, which is a similar range of coupling strengths as previously observed in comparisons with the one-component plasma. Within this range, typical relative errors of approximately 20% and worst-case relative errors of approximately 40% are observed. We also examine the Darken model, which approximates the interdiffusion coefficients based on the self-diffusion coefficients.