Reduced models for approximating the impact of kinetic electron behavior on the transport of thermal energy and magnetic field are investigated. The thermal flux limiter has improved agreement with Vlasov-Fokker-Planck data when a harmonic form is used that adjusts the electron mean free path to account for electron-electron collisions; these results apply to both unmagnetized and magnetized plasmas. Once a magnetic field is incorporated, the mean free path should also be modified using the electron gyroradius. A flux limiter on Nernst advection of magnetic fields is also required; a form that limits Nernst by the same fraction as the thermal heat-flow best reproduces kinetic simulations. A flux limiter form for the cross terms (Righi-Leduc and cross-gradient-Nernst) is also suggested. Hohlraum simulations relevant to fusion experiments on the National Ignition Facility are found to be sensitive to all of these details.

We present simulations of heat flow relevant to gaspipe experiments on the National Ignition Facility (NIF) to investigate kinetic effects on transport phenomena. D2 and neopentane (C5H12) filled targets are used to study the laser preheat stage of a MagLIF scheme where anaxial magnetic field is sometimes applied to the target. Simulations were done with the radiation-MHD code HYDRA with a collision-dominated fluid model and the Schurtz nonlocal electron thermal conduction model. Using the Schurtz model to evolve the electron temperature increased the heat front propagation of neopentane gas targets compared to a local model by limiting radial heat flow. This increases electron temperature near the axis, which decreases laser absorption. We find the effect of heat flow models on temperature profiles and laser propagation is modest. Beyond the Schurtz model, we utilize HYDRA to initialize plasma conditions for the Vlasov Fokker-Planck K2 code. We run K2 until a quasi-steady state is reached and examine the impact of kinetic effects on heat transport. Although axial heat flow is well predicted by fluid models, the fluid model consistently over predicts radial heat flow up to 150 temperature gradient of D2 filled gaspipes. On the other hand, the Schurtz nonlocal electron conduction model is found to be adequate for capturing kinetic heat flow in gaspipes.

The accelerated electron spectrum from high-intensity laser–solid interaction is often conveniently described using a Boltzmann distribution, whose temperature is known within the field as the hot-electron temperature. The importance of the electron temperature is highlighted by the sheer number of experimental and simulation studies on the subject over the past three decades. Recently, multi-kJ, multi-ps pulses have yielded electron spectra with temperatures far beyond the expected ponderomotive result. Expressions that predict the electron temperature considering laser parameters beyond intensity and wavelength have been developed, albeit using small datasets. In this review, we present what is, to the best of our knowledge, the largest dataset of electron temperatures gathered from experimental measurements and particle-in-cell simulations. This dataset allows us to compare existing analytical and empirical hot-electron temperature scaling models over a wide parameter range. We also develop new scaling models that incorporate the laser pulse duration of the laser and the plasma scale length. Three models that include pulse-duration and scale length dependence are especially successful at predicting both simulated and experimental data. The dataset will soon be made publicly available to encourage further investigation.

Indirect Drive Inertial Confinement Fusion Experiments on the National Ignition Facility (NIF) have achieved a burning plasma state with neutron yields exceeding 170 kJ, roughly 3 times the prior record and a necessary stage for igniting plasmas. The results are achieved despite multiple sources of degradations that lead to high variability in performance. Results shown here, for the first time, include an empirical correction factor for mode-2 asymmetry in the burning plasma regime in addition to previously determined corrections for radiative mix and mode-1. Analysis shows that including these three corrections alone accounts for the measured fusion performance variability in the two highest performing experimental campaigns on the NIF to within error. Here we quantify the performance sensitivity to mode-2 symmetry in the burning plasma regime and apply the results, in the form of an empirical correction to a 1D performance model. Furthermore, we find the sensitivity to mode-2 determined through a series of integrated 2D radiation hydrodynamic simulations to be consistent with the experimentally determined sensitivity only when including alpha-heating. Recent improvements in the indirect-drive inertial confinement fusion experiments include the achievement of burning plasma state. Here the authors report the scaling of neutron yield in a burning plasma of Deuterium-Tritium fusion reaction by including the mode-2 asymmetry.

It was recently shown that the use of Coulomb logarithms appropriate for bremsstrahlung radiation (rather than transport processes) along with corrections for the Langdon effect and ion screening reproduced measurements of collisional absorption in well-characterized underdense plasmas [D. Turnbull et al., Phys. Rev. Lett. 130, 145103 (2023)]. However, it was recognized at the time that the use of the standard absorption-reduction factor from Langdon's seminal paper was inconsistent with the use of Coulomb logarithms that are thermally averaged over a Maxwellian electron-velocity distribution function. A more accurate approach would be to average over the expected super-Gaussian distribution function while accounting for the Gaunt factor's velocity dependence, which somewhat mitigates the Langdon effect; however, at that time, this theory matched the data less well. This conflict is now eliminated with the additional insight that the ionization state of our mid-Z ion species (when present) was lower than had been assumed, as evidenced by the Thomson-scattering data and time-dependent Cretin simulations. We are now able to show that an improved treatment of the Langdon effect provides the best match to data. Otherwise, the prior conclusions remain unchanged. We also show an example of the substantial expected impact to the absorption rate in calculations of indirect-drive hohlraums.

On December 5, 2022, an indirect drive fusion implosion on the National Ignition Facility (NIF) achieved a target gain G_{target} of 1.5. This is the first laboratory demonstration of exceeding "scientific breakeven" (or G_{target}>1) where 2.05 MJ of 351 nm laser light produced 3.1 MJ of total fusion yield, a result which significantly exceeds the Lawson criterion for fusion ignition as reported in a previous NIF implosion [H. Abu-Shawareb et al. (Indirect Drive ICF Collaboration), Phys. Rev. Lett. 129, 075001 (2022)PRLTAO0031-900710.1103/PhysRevLett.129.075001]. This achievement is the culmination of more than five decades of research and gives proof that laboratory fusion, based on fundamental physics principles, is possible. This Letter reports on the target, laser, design, and experimental advancements that led to this result.

We provide analytic expressions for the effective Coulomb logarithm for inverse bremsstrahlung absorption which predict significant corrections to the Langdon effect and overall absorption rate compared to previous estimates. The calculation of the collisional absorption rate of laser energy in a plasma by the inverse bremsstrahlung mechanism usually makes the approximation of a constant Coulomb logarithm. We dispense with this approximation and instead take into account the velocity dependence of the Coulomb logarithm, leading to a more accurate expression for the absorption rate valid in both classical and quantum conditions. In contrast to previous work, the laser intensity enters into the Coulomb logarithm. In most laser-plasma interactions the electron distribution function is super-Gaussian [Langdon, Phys. Rev. Lett. 44, 575 (1980)], and we find the absorption rate under these conditions is increased by as much as similar to 30% compared to previous estimates at low density. In many cases of interest the correction to Langdon's predicted reduction in absorption is large; for example at Z = 6 and T-e = 400 eV the Langdon prediction for the absorption is in error by a factor of similar to 2. However, we also account for the additional effect of plasma screening, which predicts a reduction in absorption by a similar amount (up to approximate to 30%). These two effects compete to determine the overall absorption, which may be increased or decreased, depending on the conditions. The corrections can be incorporated into radiation-hydrodynamics simulation codes by replacing the familiar Coulomb logarithm with an analytic expression which depends on the super-Gaussian order "M" and the screening length.

We describe the fastVFP code for solution of the Vlasov-Fokker-Planck equation for non-local electron transport and the generation of magnetic field, especially for application to laser-produced plasmas. We describe the essential features of the code that make it fast and robust and suitable for inclusion as a transport package in a fluid simulation. We present a few sample results that demonstrate the abilities of the code.

Vlasov–Fokker–Planck simulation codes occupy an important niche in modeling laser-produced plasmas, since they are well suited to studying the effect of collisions on electron kinetic phenomena, especially energy transport. One of the most important elements of energy transport is the absorption of laser light by the plasma; however, simulating this in detail requires resolving oscillations of the laser light, whose characteristic timescale is orders of magnitude shorter than the simulation time needed to study transport physics. For this reason, most Vlasov–Fokker–Planck codes used to study electron transport in laser plasmas rely on simplified models of the laser–plasma coupling. Their underlying assumptions nominally preclude their use for modeling laser light having short-scale structure in space or time, such as broadband lasers. In this work, we derive a more general computational framework suitable for arbitrarily structured laser fields. Our approach is based on an extended set of Vlasov–Fokker–Planck equations that separately solve for the low- and high-frequency plasma response. We implement these extended Vlasov–Fokker–Planck equations in the spherical harmonic code K2 and demonstrate the performance of the method on several laser absorption test problems, with particular attention to the judicious selection of time steps, time integrators, and spherical harmonic truncation, according to the intensity and spectrum of the laser light under consideration. Comparison with the widely used Langdon absorption operator shows the Langdon operator performs remarkably well for predicting laser heating in the simple cases considered here, even in situations that would seem to violate its underlying assumptions.

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.

Strong shocks are essential components in many high-energy-density environments such as inertial confinement fusion implosions. However, the experimental measurements of the spatial structures of such shocks are sparse. In this paper, the soft x-ray emission of a shock front in a helium gas mixture (90% helium, 10% neon) and a pure neon gas was spatially resolved using an imaging spectrometer. We observe that the shock width in the helium mixture gas is about twice as large as in the pure neon gas. Moreover, they exhibit different precursor layers, where electron temperature greatly exceeds ion temperature, extending for more than ∼350µm with the helium gas mixture but less than 30µm in the pure neon. At the shock front, calculations show that the electrons are strongly collisional with mean-free path two orders of magnitude shorter than the characteristic length of the shock. However, the helium ions can reach a kinetic regime as a consequence of their mean-free path being comparable to the scale of the shock. A radiation-hydrodynamic simulation demonstrates the impact of thermal conduction on the formation of the precursors with charge state, Z, playing a major role in heat flow and the precursor formation in both the helium mixture and the pure neon gases. Particle-in-cell simulations are also performed to study the ion kinetic effects on the formation of the observed precursors. A group of fast-streaming ions is observed leading the shock only in the helium gas mixture. Both effects explain the longer precursor layer in the helium shock.

We show that the ponderomotive force associated with laser speckles can scatter electrons in a laser-produced plasma in a manner similar to Coulomb scattering. Analytic expressions for the effective collision rates are given. The electron-speckle collisions become important at high laser intensity or during filamentation, affecting both long- and short-pulse laser intensity regimes. As an example, we find that the effective collision rate in the laser-overlap region of hohlraums on the National Ignition Facility is expected to exceed the Coulomb collision rate by 1 order of magnitude, leading to a fundamental change to the electron transport properties. At the high intensities characteristic of short-pulse laser-plasma interactions (I≳10^{17}  W cm^{-2}), the scattering is strong enough to cause the direct absorption of laser energy, generating hot electrons with energy scaling as E≈1.44(I/10^{18}  W cm^{-2})^{1/2} MeV, close to experimentally observed results.

Current and future applications of intense proton sources abound, including radiography, cancer therapy, warm dense matter generation, and inertial confinement fusion. With increasingly efficient acceleration and focusing mechanisms, proton current densities may soon approach and exceed 1010 A/cm2, e.g., via intense laser drivers. Simulations have previously shown that in this current density regime, beam-induced field generation plays a significant role in beam transport through dense plasmas. Here, we present a theoretical model for the generation of resistive magnetic fields by intense proton beam transport through solid density plasmas. The theoretical evolution of the magnetic field profile is calculated using an analytic model for aluminum resistivity, heat capacity, and stopping power, applicable from cold matter to hot plasma. The effects of various beam and material parameters on the field are investigated and explained for both monoenergetic and Maxwellian proton beams. For a proton beam with Maxwellian temperature 5 MeV and total energy 10 J, the model calculates resistive magnetic fields up to 150 T in aluminum. The calculated field profiles from several beam cases are compared with 2D hybrid particle-in-cell simulations, with good agreement found in magnitude and time scale.

Thermal transport of uniformly laser-irradiated spheres of various materials is investigated computationally. One-dimensional simulations of low- to mid-Z materials (Be, Al, and Cu) are performed to evaluate the impact of nonlocal electron transport on experimental observables under laser intensities of relevance to direct-drive inertial confinement fusion. We compare thermal transport models of different levels of fidelity: flux-limited Spitzer–Harm diffusion, the Schurtz–Nicolai–Busquet (SNB) reduced-order nonlocal model, and a Fokker–Planck description. Spitzer–Harm diffusion with different flux-limiter factors are compared with different implementations of the SNB model in the HYDRA radiation hydrodynamics code. Under the conditions of interest, the peak heat flux in the thermal front with the SNB model shows good agreement with Fokker–Planck calculations, with the largest errors below 10% at 1015 W/cm2 laser intensity. From HYDRA-SNB simulations, two experimentally relevant effects are observed from nonlocal heat transport when compared to flux-limited Spitzer–Harm modeling: coronal temperatures are cooler due to reduced heat fluxes in the expanding plasma and (for mid-Z materials) x-ray emissions are enhanced due to preheating in the dense plasma.

This paper demonstrates the capability of optical Thomson scattering (OTS) to measure thermal transport, and provides support to radiation hydrodynamic and kinetic simulations of electron thermal transport and plasma evolution. OTS theory and plasma simulations are applied to the interpretation of experimental measurements of laser-produced plasma from spherical gold targets on the OMEGA laser facility. The dynamical form factor, S(k→,ω), of electron density fluctuations that is used in the fitting of Thomson scattering spectra includes ion–ion collisions and effects of non-Maxwellian distribution functions. OTS measurements and their interpretation are consistent with the nonlocal transport model in radiation hydrodynamic simulations as well as with kinetic simulations in the second half of the probe pulse duration. In particular, the reversal of heat transport during cooling is observed to be consistent with simulations, while some discrepancies are noted during the initial heating of the Au targets.

Laser-foil interactions provide a simple platform for studying physics related to hohlraums and direct-drive ICF. Proton radiography of these systems allows for probing to densities greater than critical, where the temperature gradients are steep and classical transport breaks down [1] . Quantification of magnetic flux generation around a laser-target interaction confirms that the plasma is non-local; a Biermann suppression term must be included in Gorgon simulations [3] , [4] to match the experimental magnetic flux [2] . Higher magnetic flux generation is found as the target Z is increased, due to shorter electron mean-free paths. These results allow for improved estimates for the impact of self-generated fields on both direct-drive ablation fronts and hohlraums. A direct measurement of Nernst in a laser absorption region is also sought, as this would further constrain the impact of selfgenerated magnetic fields.

Laser-plasma interactions (such as those used to drive ICF capsules) are still poorly described by radiation hydrodynamics modelling. Traditional techniques used to diagnose these plasmas are often limited to probing the under-dense regions. Proton radiography, however, can probe electric and magnetic fields in the dense plasma that transports heat into the ablation surface.

Measurements were made of the return current instability growth rate, demonstrating its concurrence with nonlocal transport. Thomson scattering was used to measure a maximum growth rate of 5.1×10^{9}  Hz, which was 3 times less than classical Spitzer-Härm theory predicts. The measured plasma conditions indicate the heat flux was nonlocal, and Vlasov-Fokker-Planck simulations that account for nonlocality reproduce the measured growth rates. Furthermore, the threshold for the return current instability was measured (δ_{T}=0.017±0.002) to be in good agreement with previous theoretical models.

We present a computational study of laser-driven ion acceleration that optimizes a combination of target transparency and extended field acceleration at moderate relativistic intensity. Our scheme applies two sequential laser pulses irradiating a thin target foil along the same direction: The first pulse drives a rapid expansion of the target, while the second one drives a quasistatic electric field in the expanding target with increasing electron temperature. In our particle-in-cell simulations we observe proton peak energies and numbers enhanced by factors of up to 3 compared with regular target-normal sheath acceleration.