Self-Similar Structure and Experimental Signatures of Suprathermal Ion Distribution in Inertial Confinement Fusion Implosions

Impact of mid-Z gas fill on dynamics and performance of shock-driven implosions at the OMEGA laser

Shock-driven implosions with 100% deuterium (D_{2}) gas fill compared to implosions with 50:50 nitrogen-deuterium (N_{2}D_{2}) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-Z fill gas on implosion performance. Ion temperature (T_{ion}) as inferred from the width of measured DD-neutron spectra is seen to be 34%±6% higher for the N_{2}D_{2} implosions than for the D_{2}-only case, while the DD-neutron yield from the D_{2}-only implosion is 7.2±0.5 times higher than from the N_{2}D_{2} gas fill. The T_{ion} enhancement for N_{2}D_{2} is observed in spite of the higher Z, which might be expected to lead to higher radiative loss, and higher shock strength for the D_{2}-only versus N_{2}D_{2} implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-Z material. This picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the N_{2}D_{2} implosions relative to the pure D_{2}-fill implosion is determined to result from the reduced amount of D_{2} in the fill (fourfold effect on yield) combined with a lower fraction of the D_{2} fuel being hot enough to burn in the N_{2}D_{2} case. The experimental yield and T_{ion} ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot D_{2} in the N_{2}D_{2} relative to the D_{2}-only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D xrage simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased T_{ion} and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the N_{2}D_{2} versus D_{2}-only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

Experimental Measurements of Ion Diffusion Coefficients and Heating in a Multi-Ion-Species Plasma Shock

Collisional plasma shocks generated from supersonic flows are an important feature in many astrophysical and laboratory high-energy-density plasmas. Compared to single-ion-species plasma shocks, plasma shock fronts with multiple ion species contain additional structure, including interspecies ion separation driven by gradients in species concentration, temperature, pressure, and electric potential. We present time-resolved density and temperature measurements of two ion species in collisional plasma shocks produced by head-on merging of supersonic plasma jets, allowing determination of the ion diffusion coefficients. Our results provide the first experimental validation of the fundamental inter-ion-species transport theory. The temperature separation, a higher-order effect reported here, is valuable for advancements in modeling HED and ICF experiments.

Thermal decoupling of deuterium and tritium during the inertial confinement fusion shock-convergence phase

A series of thin glass-shell shock-driven DT gas-filled capsule implosions was conducted at the OMEGA laser facility. These experiments generate conditions relevant to the central plasma during the shock-convergence phase of ablatively driven inertial confinement fusion (ICF) implosions. The spectral temperatures inferred from the DTn and DDn spectra are most consistent with a two-ion-temperature plasma, where the initial apparent temperature ratio, T_{T}/T_{D}, is 1.5. This is an experimental confirmation of the long-standing conjecture that plasma shocks couple energy directly proportional to the species mass in multi-ion plasmas. The apparent temperature ratio trend with equilibration time matches expected thermal equilibration described by hydrodynamic theory. This indicates that deuterium and tritium ions have different energy distributions for the time period surrounding shock convergence in ignition-relevant ICF implosions.

Species Separation and Hydrogen Streaming upon Shock Release from Polystyrene under Inertial Confinement Fusion Conditions

Shock release from inertial confinement fusion (ICF) shells poses a great challenge to single-fluid hydrodynamic equations, especially for describing materials composed of different ion species. This has been evidenced by a recent experiment [Haberberger et al., Phys. Rev. Lett. 123, 235001 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.235001], in which low-density plasmas (10^{19} to 10^{20}  cm^{-3}) are measured to move far ahead of what radiation-hydrodynamic simulations predict. To understand such experimental observations, we have performed large-scale nonequilibrium molecular-dynamics simulations of shock release in polystyrene (CH) at experimental conditions. These simulations revealed that upon shock releasing from the back surface of a CH foil, hydrogen can stream out of the bulk of the foil due to its mass being lighter than carbon. This released hydrogen, exhibiting a much broader velocity distribution than carbon, forms low-density plasmas moving in nearly constant velocities ahead of the in-flight shell, which is in quantitative agreement with the experimental measurements. Such kinetic effect of species separation is currently missing in single-fluid radiation-hydrodynamics codes for ICF simulations.

Kinetic simulations of fusion ignition with hot-spot ablator mix

Inertial confinement fusion fuel suffers increased x-ray radiation losses when carbon from the capsule ablator mixes into the hot-spot. Here, we present one- and two-dimensional ion Vlasov-Fokker-Planck simulations that resolve hot-spot self-heating in the presence of a localized spike of carbon mix, totalling 1.9% of the hot-spot mass. The mix region cools and contracts over tens of picoseconds, increasing its α particle stopping power and radiative losses. This makes a localized mix region more severe than an equal amount of uniformly distributed mix. There is also a purely kinetic effect that reduces fusion reactivity by several percent, since faster ions in the tail of the distribution are absorbed by the mix region. Radiative cooling and contraction of the spike induces fluid motion, causing neutron spectrum broadening. This artificially increases the inferred experimental ion temperatures and gives line of sight variations.

Two-dimensional implosion simulations with a kinetic particle code

We perform two-dimensional implosion simulations using a Monte Carlo kinetic particle code. The application of a kinetic transport code is motivated, in part, by the occurrence of nonequilibrium effects in inertial confinement fusion capsule implosions, which cannot be fully captured by hydrodynamic simulations. Kinetic methods, on the other hand, are able to describe both continuum and rarefied flows. We perform simple two-dimensional disk implosion simulations using one-particle species and compare the results to simulations with the hydrodynamics code rage. The impact of the particle mean free path on the implosion is also explored. In a second study, we focus on the formation of fluid instabilities from induced perturbations. We find good agreement with hydrodynamic studies regarding the location of the shock and the implosion dynamics. Differences are found in the evolution of fluid instabilities, originating from the higher resolution of rage and statistical noise in the kinetic studies.

Indications of flow near maximum compression in layered deuterium-tritium implosions at the National Ignition Facility

An accurate understanding of burn dynamics in implosions of cryogenically layered deuterium (D) and tritium (T) filled capsules, obtained partly through precision diagnosis of these experiments, is essential for assessing the impediments to achieving ignition at the National Ignition Facility. We present measurements of neutrons from such implosions. The apparent ion temperatures T_{ion} are inferred from the variance of the primary neutron spectrum. Consistently higher DT than DD T_{ion} are observed and the difference is seen to increase with increasing apparent DT T_{ion}. The line-of-sight rms variations of both DD and DT T_{ion} are small, ∼150eV, indicating an isotropic source. The DD neutron yields are consistently high relative to the DT neutron yields given the observed T_{ion}. Spatial and temporal variations of the DT temperature and density, DD-DT differential attenuation in the surrounding DT fuel, and fluid motion variations contribute to a DT T_{ion} greater than the DD T_{ion}, but are in a one-dimensional model insufficient to explain the data. We hypothesize that in a three-dimensional interpretation, these effects combined could explain the results.