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 have conducted 0D, 1D, 2D, and 3D simulations of an idealized coaxial impedance-matched Marx generator (IMG) [Phys. Rev. Accel. Beams 20, 040402 (2017)]. The 0D calculations were conducted with a four-element circuit model; the 1D, 2D, and 3D calculations were conducted with highly resolved, fully electromagnetic representations. The IMG consists of 30 stages distributed axially and connected electrically in series. Each stage is powered by two bricks separated by 180 degrees and connected electrically in parallel. Each brick comprises two opposite-polarity capacitors in series with a single switch. The bricks drive an internal impedance-matched coaxial transmission line terminated by a resistive load. The simulations neglect effects due to the switch-triggering circuit, the capacitor-charging circuit, external conducting boundaries, and reactive components of the load. We find dimensionality does not significantly affect the electrical power delivered by the IMG to its load: peak load powers estimated by the 0D, 1D, 2D, and 3D simulations agree to within 1%. The 3D calculations demonstrate that electromagnetic power radiated by the bricks, and axial gaps between stages, reduces the peak load power by less than similar to 1%. Each simulation assumes the load impedance is 34% above that at which the load power is maximized. Operating an IMG with such an overmatched load offers several advantages while decreasing the peak load power by only 2%. The 0D, 1D, 2D, and 3D models outlined herein could be adapted to assess computationally competing IMG designs, and conduct a variety of numerical IMG experiments, before an IMG is constructed.

Magnetohydrodynamic (MHD) simulations of pulsed power experiments frequently result in unphysical runaway heating when extended to low values of mass density. Traditionally, this has been addressed by the use of conductivity floors below which the plasma is given an arbitrarily low conductivity value. Here, a low-density treatment is presented that allows for low-density material to carry current while maintaining a stable temperature. This treatment is implemented in the Ares multiphysics code. Numerical modifications to the conventional vacuum treatment include an energy-conserving density floor, a modified averaging procedure to determine the thermal conductivity at the edge of adjacent zones, and a modified averaging procedure for determining the thermal conductivity in mixed zones that contain material from differing regions. Additionally, anomalous resistivity and Bohm diffusion models are implemented as simplified models for microphysics-induced enhancement of collisional transport. The advantage of these various improvements are illustrated in a simple 1-D pulsed power target where the combined changes result in stable temperatures within the lower density regions.

Pulsed power generators create high-energy-density conditions by rapidly delivering an immense pulse of electrical current to a compact imploding load. Accurately measuring the shape and amplitude of this load current pulse is essential to understanding the behavior of all pulsed power experiments. At the Z Pulsed Power Facility, the closest-in load current measurements are provided by velocimetry techniques such as VISAR (velocity interferometer system for any reflector) and PDV (photonic Doppler velocimetry). Here, fiber-coupled interferometers measure the velocity history of an exploding metallic flyer plate that is embedded in the vertical walls of the current return can. The flyer plate is driven outward by the magnetic pressure from the load current such that magnetohydrodynamic modeling can be used to determine the load current waveform from the measured velocity history. In this paper, we present the first load current velocimetry measurements to be made from the horizontal top flyer plate that carries current radially inward from the return can to the load. These spatially resolved measurements, which span R = 5- 9 mm, are enabled by a transformative new velocimetry diagnostic-a line-imaging velocity interferometer called Z Line VISAR (ZLV)-whose optical performance overcomes the measurement challenges presented by the steep velocity gradients encountered on the top flyer plate. To validate ZLV's capabilities, a 14-MA, 100-ns experiment was conducted to losslessly couple current up the return can and radially inward across the top flyer plate. Comparisons between the ZLV data obtained from this experiment and two-dimensional magnetohydrodynamic simulations driven with the current measured on the return can indicate that the current delivery across the top flyer plate is indeed lossless to within the few-percent uncertainty of the ZLV data. Given that the current coupling is lossless, the experimental results are used to demonstrate that one-dimensional current unfold techniques can be applied to generate a radially resolved load current map from the ZLV velocity data. This analysis provides a template for how to use the ZLV diagnostic to determine the efficacy of current delivery in future experiments where losses may occur in close proximity to the load.

The hydrodynamic instability growth of a reshocked single-mode interface between high energy density fluids is studied. A laser-driven shock wave is used to drive an initially solid, sinusoidal interface between a dense plastic (1.43 g/cc) and a light foam (≈ 0.110 g/cc). After the interface has grown to a nonlinear state where the amplitude is of order of the wavelength, it is reshocked. The reshock compresses the nonlinear perturbation, which then grows at about twice the rate. While the pre-reshock growth rate is sensitive to the initial amplitude and wavelength of the perturbation, the post-reshock growth rate is comparatively insensitive to the initial condition. Qualitatively, we observe that the perturbations are less coherent after reshock, consistent with the idea that having a reshock accelerates the transition to turbulence. We find that some memory of the initial condition remains, even after reshock at late time: it appears if the initial perturbations have large enough wavelengths, and the flow structure of size comparable to the initial wavelength persists through reshock. Our results agree with design simulations and are consistent with the phenomenology of reshock studies in conventional gaseous shock tubes.

An initially planar shock wave propagating into a medium of nonuniform density will be perturbed, leading to the generation of postshock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly nonuniform density (order-unity normalized variance, σ_{ρ}/ρ[over ¯]∼1) and strong shocks (shock Mach numbers M[over ¯]_{s}≳10). This leads to a highly disrupted shock and a turbulent postshock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the postshock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the preshock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The postshock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.

A transverse magnetic field can suppress the hydrodynamic instability growth of an unstable plasma interface. This effect, of interest to inertial confinement fusion and astrophysics, has mostly been studied via simulation. Here, we present the design of an experiment at the National Ignition Facility (NIF) to demonstrate this effect in a laboratory. Simulations indicate that the timescale for the diffusion of the magnetic field across the mixing region should be at least comparable to the timescale of the instability growth in order to have a measurable suppression effect. This motivates the use of lower density target materials than usual high energy density (HED) hydrodynamics experiments to permit faster hydrodynamics and higher plasma conductivities (through higher temperature), for a given laser drive and magnetic field. We discuss a target design for creating a Rayleigh–Taylor unstable HED plasma interface that uses 320 mg/cc iodine-doped carbon foam as the heavy material, 20 mg/cc carbon foam as the light material, and a 6 μm amplitude, a 120 μm wavelength ripple machined at the interface, which shows a measurable suppression effect with a nominal NIF drive and a 30 T magnetic field (the present facility limit). Models indicate lower density foams that may display even larger suppression effects, as the Hohlraum drive also radiatively preheats the foam to permit even higher temperatures and, hence, higher conductivities.

This paper describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer-Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy-dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Peclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Peclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

We describe a new radiography–based technique for measuring the width of the mixing layer that develops when high–energy–density (HED) fluids of different density mix at an interface. In a traditional radiography target, the denser material includes an embedded, density–matched, high opacity, tracer layer to focus the radiograph on a 2D slice of the 3D flow field. The high opacity layer enhances the contrast between the dense and light materials that permits a clear identification of the spike front, i.e. penetration of the dense material into the light material, but obscures the bubble front, i.e. penetration of the light material into the heavy when the fluids are highly interpenetrated, as in a turbulent mixing layer. On the other hand, removing the tracer layer makes it difficult to identify the interface at all, due to poor contrast and an increased susceptibility to edge effects. Our solution is a split tracer approach where half the target uses the traditional geometry while the other half utilizes a special tracer layer in the light material to produce a radiograph with inverse x-ray opacity characteristics when backlit at a specified x–ray energy. Hydrodynamic equivalency is maintained between the two halves of the target. By making a single target with the traditional and inverse material pairs set side-by-side, we are able to simultaneously measure the penetration of bubbles and spikes in a deeply nonlinear flow, yielding the full mix width. We discuss this new target concept and present our first experiments at the National Ignition Facility (NIF) demonstrating the viability of this technique.

A line VISAR (Velocity Interferometer System for Any Reflector) has been designed and commissioned at the Sandia National Laboratory's Z-machine. The instrument consists of an F/2 collection system, beam transport, and an interferometer table that contains two Mach-Zehnder type interferometers and an eight channel Gated Optical Imaging (GOI) system. The VISAR probe laser operates at the 532 nm wavelength, and the GOI bandpass is 540-600 nm. The output of each interferometer is passed to an optical streak camera with four selectable sweep speeds. The system is designed with three interchangeable optics modules to select a full field of view of 1 mm, 2 mm, or 4 mm. The optical beam transport system connects the target image plane to the interferometers and the gated optical imagers. The target is integrated into a sacrificial final optics assembly that is integral to the transport beamline.

Pulsed power accelerators compress electrical energy in space and time to provide versatile experimental platforms for high energy density and inertial confinement fusion science. The 80-TW "Z" pulsed power facility at Sandia National Laboratories is the largest pulsed power device in the world today. Z discharges up to 22MJ of energy stored in its capacitor banks into a current pulse that rises in 100ns and peaks at a current as high as 30 MA in low-inductance cylindrical targets. Considerable progress has been made over the past 15years in the use of pulsed power as a precision scientific tool. This paper reviews developments at Sandia in inertial confinement fusion, dynamic materials science, x-ray radiation science, and pulsed power engineering, with an emphasis on progress since a previous review of research on Z in Physics of Plasmas in 2005.

At National Ignition Facility (NIF), yield amplification due to alpha particle heating approached ~3 in the highest performing inertial confinement fusion (ICF) implosions, while yield amplification of ~15-30 is needed for ignition. Hydrodynamic instabilities are a major factor in degradation of implosions while understanding and mitigation of the instabilities are critical to achieving ignition. This article describes recent and planned hydrodynamic instability experiments with several focused platforms that have been developed to directly measure these instabilities in all phases of ICF implosions. Measurements of ripple-shock generation at OMEGA laser have indicated initial seeds for the instabilities in three ablators - plastic (CH), beryllium, and high-density carbon (HDC). Hydrodynamic Growth Radiography (HGR) platform was used to measure instability growth at the ablation front in the acceleration phase of implosions. This platform used pre-imposed 2-D perturbations for growth factor measurements at different perturbation wavelengths and was also used to measure growth of “native roughness” modulations, fill tubes, and capsule support membranes or “tents”. Also, in the acceleration phase several new experimental platforms have been or are being developed to measure instability growth at the ablator-ice interface. In the deceleration phase of implosions, “self-emission” and “self-backlighting” platforms were developed to measure perturbations near peak compression. This article reviews recent progress and results.

Hydrodynamic instabilities are a major factor in degradation of inertial confinement fusion (ICF) implosions. In the highest performing implosions on National Ignition Facility, yield amplification (YA) due to alpha particle heating approached ∼3, while YA of ∼15–30 is needed for ignition. Understanding and mitigation of the instabilities are critical to achieving ignition. This article reviews several experimental platforms that have been developed to directly measure these instabilities in all phases of ICF implosions. Measurements of ripple-shock propagation at OMEGA laser has provided results on initial seeds for the instabilities in three ablators—plastic (CH), beryllium, and high-density carbon. At the ablation front, instability growth of pre-imposed modulations was measured in the linear regime using the hydrodynamic growth radiography platform. This platform was extended for modulation growth of ‘native roughness’ modulations and engineering features (fill tubes and capsule support membranes or ‘tents’). Several new experimental platforms have or are being developed to measure instability growth at the ablator–ice interface. In the deceleration phase of implosions, complementary ‘self-emission’ and ‘self-backlighting’ platforms were developed to measure low-mode asymmetries and high-mode perturbations near peak compression.

The Rayleigh-Taylor (RT) instability occurs at an interface between two fluids of differing density during an acceleration. These instabilities can occur in very diverse settings, from inertial confinement fusion (ICF) implosions over spatial scales of [Formula: see text] cm (10-1,000 μm) to supernova explosions at spatial scales of [Formula: see text] cm and larger. We describe experiments and techniques for reducing ("stabilizing") RT growth in high-energy density (HED) settings on the National Ignition Facility (NIF) at Lawrence Livermore National Laboratory. Three unique regimes of stabilization are described: (i) at an ablation front, (ii) behind a radiative shock, and (iii) due to material strength. For comparison, we also show results from nonstabilized "classical" RT instability evolution in HED regimes on the NIF. Examples from experiments on the NIF in each regime are given. These phenomena also occur in several astrophysical scenarios and planetary science [Drake R (2005) Plasma Phys Controlled Fusion 47:B419-B440; Dahl TW, Stevenson DJ (2010) Earth Planet Sci Lett 295:177-186].

The Rayleigh-Taylor (RT) instability is a common occurrence in nature, notably in astrophysical systems like supernovae, where it serves to mix the dense layers of the interior of an exploding star with the low-density stellar wind surrounding it, and in inertial confinement fusion experiments, where it mixes cooler materials with the central hot spot in an imploding capsule and stifles the desired nuclear reactions. In both of these examples, the radiative flux generated by strong shocks in the system may play a role in partially stabilizing RT instabilities. Here, we present experiments performed on the National Ignition Facility, designed to isolate and study the role of radiation and heat conduction from a shock front in the stabilization of hydrodynamic instabilities. By varying the laser power delivered to a shock-tube target with an embedded, unstable interface, the radiative fluxes generated at the shock front could be controlled. We observe decreased RT growth when the shock significantly heats the medium around it, in contrast to a system where the shock did not produce significant heating. Both systems are modeled with a modified set of buoyancy-drag equations accounting for ablative stabilization, and the experimental results are consistent with ablative stabilization when the shock is radiative. This result has important implications for our understanding of astrophysical radiative shocks and supernova radiative hydrodynamics [Kuranz et al., Nature Communications 9(1), 1564 (2018)].