A decades-long quest to achieve fusion energy target gain and ignition in a controlled laboratory experiment, dating back to 1962, has been realized at the National Ignition Facility (NIF) on December 5, 2022 [Abu-Shawareb et al., Phys. Rev. Lett. 132, 065102 (2024)] where an imploded pellet of deuterium and tritium (DT) fuel generated more fusion energy (3.15 MJ) than laser energy incident on the target (2.05 MJ). In these experiments, laser beams incident on the inside of a cylindrical can (Hohlraum) generate an intense 3 × 106 million degree x-ray radiation bath that is used to spherically implode 2 mm diameter pellets containing frozen deuterium and tritium. The maximum fusion energy produced in this configuration to date is 3.88 MJ using 2.05 MJ of incident laser energy and 5.2 MJ using 2.2 MJ of incident laser energy, producing a new record target gain of 2.4×. This paper describes the physics (target and laser) design of this platform and follow-on experiments that show increased performance. We show robust megajoule fusion energy output using this design as well as explore design modification using radiation hydrodynamic simulations benchmarked against experimental data, which can further improve the performance of this platform.

Fusion offers the promise of clean and limitless energy, which could be the long-term answer to rapidly growing global energy demands. However, net fusion energy gain has not yet been realized in the laboratory due to the difficulties in creating and confining the extreme plasma conditions needed for fusion to occur. Overcoming the electrostatic repulsion of the fusion reactants requires tens of millions of degrees, and creating these conditions in the laboratory has previously required more heating energy, or driver energy, than is produced from fusion. Generating more fusion energy delivered to the target is an important milestone in the quest for fusion energy and has been a long-standing goal of all fusion energy approaches. On December 5, 2022, we have finally achieved that goal in an experiment on the National Ignition Facility (NIF)1 in Livermore, CA, which was also the last remaining milestone for the National Ignition Facility (NIF).2 Since then, we have repeated this result and achieved target gain exceeding unity five times to date at the NIF.

In the laser indirect-drive approach,3–8 lasers are shined onto the inside of a Au-lined Depleted Uranium (DU) cylindrical cavity (Hohlraum) to generate a 3 million degree x-ray radiation bath. A hollow diamond sphere containing a cryogenically frozen layer of DT, with a DT gas core, is suspended inside the Hohlraum using a thin membrane (capsule support “tent”). The x-ray radiation cavity heats the outside of the diamond shell causing it to expand radially outward which launches (implodes) the fuel and a fraction of the diamond mass radially inward, doing mechanical PdV work on the fuel and squeezing it to smaller volumes, which results in an increase in compression and temperature, see Fig. 1.

FIG. 1.

(left) Schematic of the target configuration showing the Au-lined DU hohlraum and spherical diamond capsule that contains DT fuel, with an expanded pie-diagram of the capsule and fuel for N221204 compared to N210808. Radiation hydrodynamics simulations of the integrated hohlraum and capsule (top) show the material positions and Au plasma filling into the hohlraum at time = 6 ns for N221204 (right) compared to N210808 (left). The radiation drive symmetry was controlled by transferring energy between laser beams through wavelength detuning (Δλ) and by precisely defining the balance of laser powers between the inner (23° and 30°) and outer (44° and 50°) beams throughout the duration of the laser pulse (top right). The amount of transfer was increased for N221204 compared to N210808 via increasing Δλ from 1.8 Angstroms to 2.75 Angstroms to maintain symmetry with the longer laser pulse and thicker ablator. The “Peak” of the laser power was extended by about 300 ps to drive the thicker diamond capsule and the “trough” was lengthened to preserve the shock timing. (bottom right) Simulated hohlraum internal radiation temperature (Trad) as a function of time for N221204, together with the simulated scaled capsule radii vs time for the N221204 radiation drive with the thinner (N210808) and thicker (N221204) diamond capsules.

FIG. 1.

(left) Schematic of the target configuration showing the Au-lined DU hohlraum and spherical diamond capsule that contains DT fuel, with an expanded pie-diagram of the capsule and fuel for N221204 compared to N210808. Radiation hydrodynamics simulations of the integrated hohlraum and capsule (top) show the material positions and Au plasma filling into the hohlraum at time = 6 ns for N221204 (right) compared to N210808 (left). The radiation drive symmetry was controlled by transferring energy between laser beams through wavelength detuning (Δλ) and by precisely defining the balance of laser powers between the inner (23° and 30°) and outer (44° and 50°) beams throughout the duration of the laser pulse (top right). The amount of transfer was increased for N221204 compared to N210808 via increasing Δλ from 1.8 Angstroms to 2.75 Angstroms to maintain symmetry with the longer laser pulse and thicker ablator. The “Peak” of the laser power was extended by about 300 ps to drive the thicker diamond capsule and the “trough” was lengthened to preserve the shock timing. (bottom right) Simulated hohlraum internal radiation temperature (Trad) as a function of time for N221204, together with the simulated scaled capsule radii vs time for the N221204 radiation drive with the thinner (N210808) and thicker (N221204) diamond capsules.

Close modal

A main challenge is doing sufficient work on the central plasma to ignite the fuel. For this to occur, extreme temperatures are required to balance energy losses in the system, e.g., Bremsstrahlung x-ray losses, conduction losses, and explosion phase losses. Self-heating from the alpha particles that are born in the fusion reactions is required to increase the plasma temperature to ignition conditions [D + T n (14.1 MeV)+4He (3.5 MeV)]9 which requires achieving sufficient hotspot densities to stop the alpha particles in the central hot plasma. Once the central hot plasma ignites, a burn wave propagates through a cooler denser shell of DT fuel that both initially does work on the implosion and confines the plasma long enough to enable “burning-up” of the remaining cooler denser DT fuel. The plasma is confined for 100 ps via its own inertia and then explodes under the extreme pressure of the fusing plasma. On December 5, 2022, 4% of the DT fuel present in the experiment was “burned,” up, or “fused,” creating energetic neutron particles. This paper describes the design that enabled this result together with simulated projections for future improvements.

Early designs at the NIF8 using plastic (CH) capsules to hold the DT fuel suffered from hydrodynamic instabilities seeded by the “capsule support tent”11 which caused capsule material to enter the hot-fusing plasma and radiate energy away. Increasing the laser power in the early part of the laser pulse (“CH HiFoot”) resulted in increased stability12 but was ultimately also limited by the tent perturbation. Introducing diamond (High Density Carbon (HDC)) ablators13–17 reduced sensitivity to the tent perturbation due to its higher density (3.5 g/cm3) and also made symmetry control easier due to shorter laser pulses. Due to the higher density of HDC, thinner ablators could be used, which reduces the transit time of multiple shocks moving through the ablator and fuel, resulting in shorter laser pulses.

The early hohlraum designs were also limited by laser plasma interactions18 which resulted in inefficiency of driver coupling to the hohlraum. This was mainly a result of a high level of helium gas fill that was used to hold back the expanding wall plasma, but resulted in increased levels of Stimulated Raman Scattering (SRS). Reducing the helium gas fill density inside the hohlraum19–21 resulted in lower levels of laser backscatter and an increase in hohlraum efficiency, which enabled fielding larger case-to-capsule-ratio (CCR) hohlraums and making symmetry easier to control without the use of intentional cross-beam-energy-transfer (CBET).18,22

While imploding the capsules, it is important to control the time-varying radiation flux symmetry within the hohlraum to 1%, to efficiently do work on the implosion.23–26 The laser–wall interaction causes the hohlraum wall to ingress into the path of laser beams which makes controlling symmetry a main challenge. A main method for controlling symmetry is using CBET to move energy between laser beams by tuning the frequency of the density perturbation caused by the ponder-motive pressure from the interference pattern of the laser beams. Changing the relative wavelength of the overlapping laser beams changes the frequency of this perturbation and creates a refractive index modulation which acts as a Bragg scatter and redirects the incoming laser light. Previously, using this technique in hohlraums with high levels of He gas fill resulted in increased levels of laser-backscatter and loss of symmetry control.27 Development of CBET in low-gas-filled hohlraums28,29 has been instrumental in enabling symmetry control for reduced CCR target designs.

The HYBRID-E design30–35 was one of a series of designs36 aimed to increase the energy delivered to the hotspot plasma by increasing the scale (S) of the implosion 10%–15% while maintaining the hotspot pressure required for inertial confinement fusion and other critical implosion properties, within the current limitations of the NIF. Driving larger scale implosions to similar velocities as smaller scale implosions while using the same amount of laser energy requires developing more efficient smaller hohlraums with less radiation loss. This results in a smaller case to capsule ratios (CCR), or smaller hohlraum diameter for a given capsule size, which is more challenging for controlling implosion symmetry. Advancements in symmetry control using CBET in low-gas-filled hohlraums28,29 and development of a semi-empirical model for late-time hohlraum symmetry37 calibrated against years of experimental data were instrumental in the design of HYBRID-E. Even with a more efficient hohlraum, design tradeoffs were necessary between key implosion metrics and what could be achieved using a laser-driven hohlraum. It was expected that the simulations (HYDRA38) coupled with past experimental data could get close in determining the scale of implosion that would best balance these tradeoffs. Since simulations cannot predict to the degree of precision required to delicately balance the hohlraum and implosion design tradeoffs and achieve ignition, testing two implosion scales (initial inner radius of 1050 and 1100  μm) was part of the initial HYBRID-E experimental plan.

The larger scale diamond capsule was initially tested in HYBRID-E but suffered from simultaneously trying to control the implosion velocity, coast-time, and symmetry.31 Following these experiments, a burning plasma32,33 was achieved by reducing the initial scale of the diamond capsule by only 50  μm, enabling more optimal balancing of these tradeoffs. Then, further reducing the size of the laser entrance holes, holes at the ends of the hohlraum where the laser beams enter, resulted in an even more efficient hohlraum which produced the same radiation temperature at lower laser power due to lower radiation losses. This enabled increasing the laser pulse length to improve the implosion properties39 and couple more driver energy to the implosion which led to the first experiment to achieve scientific ignition40 on August 8, 2021 [NIF shot No. N210808]34,41,42 and the first to produce >1 MJ of fusion energy, see Fig. 2.

FIG. 2.

Fusion energy produced as a function of date for Deuterium Tritium (DT) layered experiments on the NIF, color coded by design with HYBRID-E in magenta. Inset are measured neutron images10 of the hot-spot plasma burn-averaged over peak neutron production for shots N220919 (September 19, 2022), N221204 (December 5, 2022), and N230729 (July 29, 2023). The symmetry was adjusted between N220919 and N221204 which enabled the fusion energy produced to exceed the amount of laser energy into the target (blue dashed line). Repeat experiment N230729, despite being less symmetric, resulted in higher yield (3.88 MJ) and target gain (1.9) as a result of improved diamond capsule quality.

FIG. 2.

Fusion energy produced as a function of date for Deuterium Tritium (DT) layered experiments on the NIF, color coded by design with HYBRID-E in magenta. Inset are measured neutron images10 of the hot-spot plasma burn-averaged over peak neutron production for shots N220919 (September 19, 2022), N221204 (December 5, 2022), and N230729 (July 29, 2023). The symmetry was adjusted between N220919 and N221204 which enabled the fusion energy produced to exceed the amount of laser energy into the target (blue dashed line). Repeat experiment N230729, despite being less symmetric, resulted in higher yield (3.88 MJ) and target gain (1.9) as a result of improved diamond capsule quality.

Close modal

On December 5, 2022 (shot #N221204),35,43–45 fusion energy target gain exceeding unity was achieved, where the fusion energy produced was greater for the first time than laser energy on target, after a symmetry correction from September 19, 2022 (fusion energy of 1.22 MJ and target gain of 60%), by making additional design changes to the HYBRID-E platform.35,43–45 A thicker diamond ablator was fielded to couple 8% more available laser energy46,47 to the implosion and increase the areal density of the DT fuel and capsule at maximum compression which leads to higher yield amplification and fuel burn-up fraction. The extra laser energy was enabled through upgrades to the optics damage mitigation and reprocessing. The new capability of NIF delivered 150 kJ additional laser energy (2.05 MJ) compared to prior experiments (1.9 MJ) at a maximum laser power of 440 TW, see Fig. 2 for laser power as a function of time for N210808 compared to N221204. The fusion energy produced was exceeded on July 29, 2023 (N230729) in this design using a higher quality diamond capsule, achieving 3.88 MJ of fusion energy and a target gain of 1.9×.

Figure 1 shows a schematic of the target and laser design. Laser beams are pointed onto the inside of a gold-lined Depleted Uranium (DU) cylindrical can (hohlraum), 11.24 mm length × 6.4 mm diameter, in four distinct laser sets defined by their angle of incidence (inner beams: 23° and 30° and outer beams: 44° and 50°). The laser heating of the hohlraum creates a radiation oven that is 3 × 106 K (>300 eV), see Fig. 1(bottom right). Compared to N210808, the laser energy was increased from 1.92 to 2.05 MJ at constant peak laser power (440 TW), see Fig. 1 (top right), and the early part of the laser pulse was extended to maintain the shock timing for the thicker capsule.

The radiation drive is absorbed in the outer region of a hollow diamond spherical capsule (1050  μm inner radius and 85  μm thick) sitting in the center of the hohlraum, which contains the DT fuel (65  μm thick ice layer with a DT gas core). The x-ray absorption heats and ablates the diamond which expands radially outward, causing the remaining diamond capsule and DT fuel to accelerate inward at extreme implosion velocities (vimp 380–400 km/s). Compared to the N210808 HYBRID-E design, the diamond ablator thickness was increased by 5.75  μm (N221204) and 6.9  μm (N230729) to provide higher energy coupling of the radiation drive to the hotspot, higher total areal density at peak compression, and improved stability. See Ref. 35 for additional design details. The thicker diamond ablator was used to make optimal use of the extra drive at the end of the pulse which would have otherwise not coupled significantly more driver energy to the implosion. Figure 1 (bottom right) shows the simulated implosion radius vs time for the 2.05 MJ laser pulse and the 6  μm thinner N210808 diamond ablator compared to the thicker (“as-shot”) N221204 diamond ablator. The thicker ablator causes the capsule to remain larger longer and absorb more of the extra drive at the end of the pulse, which would not have coupled significantly more to the thinner shell implosion. A thickness increase in 6  μm was initially chosen to balance the improvement in both ignition margin and total areal density. See Table I for a comparison of simulated key implosion metrics for the extended 2.05 MJ drive and same capsule thickness as N210808 vs using a 6  μm thicker capsule. Future experiments will test the limits of thicker capsule implosions in this platform to further increase yield-amplification (extra yield from self-heating).

TABLE I.

Key simulated implosion metrics for 2.05 MJ laser energy EL with and without thickening the capsule by  6  μm, normalized to the simulated implosion metrics for N210808.

Values rel. to N210808 2.05 MJ EL with the same capsule 2.05 MJ EL with  6  μm thicker capsule
Yield  1.1×  2.5× 
Diamond mass remaining  1.0×  +1.5% 
Implosion velocity (km/s)  1.0×  −15 km/s 
Coupled energy to the implosion  1.0×  +20 kJ 
Ignition margin (EP2 1.0×  >1.2× 
Total areal density  1.01×  >1.18× 
Values rel. to N210808 2.05 MJ EL with the same capsule 2.05 MJ EL with  6  μm thicker capsule
Yield  1.1×  2.5× 
Diamond mass remaining  1.0×  +1.5% 
Implosion velocity (km/s)  1.0×  −15 km/s 
Coupled energy to the implosion  1.0×  +20 kJ 
Ignition margin (EP2 1.0×  >1.2× 
Total areal density  1.01×  >1.18× 

The symmetry of the radiation drive becomes more difficult to control with longer laser pulses and more time for the plasma to fill the hohlraum which reduces the energy coupling to the hotspot.23 To adjust the symmetry, the amount of laser energy transfer from the “outer” to the “inner” beams (CBET) was increased by further detuning the relative wavelengths of these beams (Δλ = Inner beam—Outer beam before laser frequency tripling),18,22,28,29 where Δλ = 1.8  Å for N210808 and Δλ = 2.75  Å for N221204. Figure 1 (left) shows radiation hydrodynamics simulations of N221204 (right) and N210808 (left) at 6 ns after the start of the laser pulse, together with the calculated positions of the material boundaries [gold-lined depleted uranium hohlraum (orange), HDC ablator (light gray), and the position of the ice-gas interface (red)]. Simulated laser rays are also shown, where the color corresponds to normalized laser powers after CBET. Precisely adjusting the powers on the inner and outer laser cones during the entire laser pulse [Fig. 1(top right)] enabled controlling fluctuations in the radiation drive uniformity which can induce ρR variations in the compressed fuel. Detailed radiation hydrodynamic simulations using HYDRA38 with inline CBET48 were used to design these inner and outer powers. These simulations were benchmarked against a series of “tuning” experiments, see the  Appendix.

The first experimental test of this new design with 2.05 MJ of laser energy (N220919) was limited in the amount of CBET “allowed” due to risk of laser damage from Stimulated Brillouin Scattering (SBS)49 of the incident laser light back out of the hohlraum into the optics assembly. This resulted in a drive deficit at the waist of the hohlraum and an oblate hotspot, see Fig. 2. Increasing the amount of transfer through wavelength separation by 0.25 Angstroms from 2.5 Angstroms to 2.75 Angstroms and readjusting the time-dependent laser powers enabled achieving a round hotspot on N221204. On July 29, 2023 (N230729), a new diamond capsule “batch” was tested which had a lower amount of tungsten dopant,  0.4% vs  0.62% tested on N221204, see Fig. 3 for a pie-diagram of the diamond capsules noting the differences for these experiments. Due to the thicker ablator the early part of the laser pulse was extended by  30 ps and a slightly higher wavelength separation of 2.9 Angstroms was used to correct symmetry for the extra diamond “payload.” Thicker capsules take longer to implode and sample a longer duration of the pole-hot radiation flux asymmetry which drives the hotspot more oblate. Increasing the wavelength separation counteracts this by making the drive more waist-hot earlier in the pulse. Small extrapolations in symmetry are predicted for these changes using radiation-hydrodynamic simulations and experimentally measured sensitivities. However, there is a larger uncertainty in extrapolating to higher levels of CBET in this regime and the chosen wavelength separation was too large, producing a prolate implosion, see Figs. 2 and 3.

FIG. 3.

(left) Pie diagrams (not drawn to scale) noting the differences between diamond capsule “batches” shot on December 5, 2022 and July 29, 2023. Main differences include the level of tungsten dopant and target quality. (middle) Radiography of the two diamond capsule batches showing a high number of tungsten-carbide “chunk” features in the tungsten-doped region of the diamond ablator for the December 5, 2022 capsules compared to the highly uniform July 29, 2023 capsules. Both shells had similar quality in the outer un-doped diamond layer. (right) Low mode asymmetry of the July 29, 2023 experiment (top: measured mode one of the hot spot velocity (km/s) and bottom: measured mode two of the hot spot shape via 3D neutron reconstruction.10 Despite simulations predicting lower performance for the lower tungsten concentration July 29, 2023 capsules and worse low mode symmetry compared to the December 5, 2022 experiment, the July 29, 2023 experiment achieved record performance for 2.05 MJ of laser energy (3.88 MJ fusion energy and target gain 1.9×) currently hypothesized to be due to the improved capsule uniformity.

FIG. 3.

(left) Pie diagrams (not drawn to scale) noting the differences between diamond capsule “batches” shot on December 5, 2022 and July 29, 2023. Main differences include the level of tungsten dopant and target quality. (middle) Radiography of the two diamond capsule batches showing a high number of tungsten-carbide “chunk” features in the tungsten-doped region of the diamond ablator for the December 5, 2022 capsules compared to the highly uniform July 29, 2023 capsules. Both shells had similar quality in the outer un-doped diamond layer. (right) Low mode asymmetry of the July 29, 2023 experiment (top: measured mode one of the hot spot velocity (km/s) and bottom: measured mode two of the hot spot shape via 3D neutron reconstruction.10 Despite simulations predicting lower performance for the lower tungsten concentration July 29, 2023 capsules and worse low mode symmetry compared to the December 5, 2022 experiment, the July 29, 2023 experiment achieved record performance for 2.05 MJ of laser energy (3.88 MJ fusion energy and target gain 1.9×) currently hypothesized to be due to the improved capsule uniformity.

Close modal

Shot N230729 also had a worse mode one symmetry compared to N221204 (107 vs 36 km/s), largely due to an imbalance in laser delivery, aligned with other sources of mode one asymmetry (capsule non-uniformity, DT thickness non-uniformity, non-axially symmetric diagnostic windows). The impact of the combined sources of mode one asymmetry is measured as a net hotspot velocity, or speed at which the hotspot is moving in one direction. The measured hotspot velocity is amplified as the temperature rises for igniting implosions that are expanding during peak neutron production. We estimate this effect as well as the impact on performance for the measured level of mode one on N230729, see Fig. 5.

The lower tungsten dopant capsule fielded on N230729 was predicted to have lower performance (by up to 40%) compared to the higher tungsten dopant capsule of N221204 primarily due to the lower calculated total areal density which causes enhanced expansion through the thin spots following ignition which limits confinement. Thin spots in the dense shell are a result of low-mode asymmetries. Therefore, the level of reduction in the performance for the lower tungsten doped capsules increases with the level of low mode asymmetry. This can be seen in Fig. 4 (top) which shows the performance of the N210808 design (black), N221204 design (red), and N221204 design with lower tungsten concentration diamond capsules (blue stars) as a function residual kinetic energy, the amount of left over work that was not done on the hotspot as a result of the low mode asymmetries creating an inefficient piston. The higher optical depth capsules have more margin against RKE in the yield range where we are currently operating, with the symmetry-performance “cliff” at larger levels of RKE.

FIG. 4.

(top) Simulated neutron yield as a function of residual kinetic energy (RKE), or amount of energy not coupled to the hot spot as a result of low-mode asymmetries from simulations which do not include alpha-heating (burn-off). The red points and curve correspond to the N221204 higher-dopant capsule design, the black points and curve corresponding to the N210808 design, and the blue stars correspond to the N221204 design with lower optical depth (N230729). (bottom) Simulated sensitivity of neutron yield hot spot oblateness (-P2) from α-off simulations for the N221204 design (red points), where a time-varying flux asymmetry was adjusted in the simulations to change hot spot P2. The measured hot spot P2 was -12  μm for N220919 and 0  μm for N221204. The range of yields at a given level of hot spot P2 results from different assumptions about the magnitude of the time-varying P2 of the hot spot and DT shell, which will be benchmarked against data in upcoming experiments.

FIG. 4.

(top) Simulated neutron yield as a function of residual kinetic energy (RKE), or amount of energy not coupled to the hot spot as a result of low-mode asymmetries from simulations which do not include alpha-heating (burn-off). The red points and curve correspond to the N221204 higher-dopant capsule design, the black points and curve corresponding to the N210808 design, and the blue stars correspond to the N221204 design with lower optical depth (N230729). (bottom) Simulated sensitivity of neutron yield hot spot oblateness (-P2) from α-off simulations for the N221204 design (red points), where a time-varying flux asymmetry was adjusted in the simulations to change hot spot P2. The measured hot spot P2 was -12  μm for N220919 and 0  μm for N221204. The range of yields at a given level of hot spot P2 results from different assumptions about the magnitude of the time-varying P2 of the hot spot and DT shell, which will be benchmarked against data in upcoming experiments.

Close modal

The amount of RKE in Fig. 4 (top) includes the simulated time-varying modes one, two, and four of the radiation drive but may overestimate the level of fluctuation in the P2 flux asymmetry, or P2 “swing” [e.g., Fig. 4 (bottom)]. Historically, simulations are benchmarked against in-flight radiography measurements (at implosion radii of 200  μm) which calibrate the level of P2 swing before the experiment, but at the time of writing this paper, a calibration experiment for N221204 had not yet been completed. Figure 4 (bottom) shows the simulated neutron yield (red points) as a function of hotspot oblateness (P2) for the N221204 design. The range of simulated neutron yields at a given level of hotspot oblateness is from taking a plausible range for the P2 swing. The insets are simulated contours at peak compression from α-off simulations (simulations where α-deposition is artificially turned off (“burn off” or α-off)) of the hotspot temperature (left) and density (right) for example points. For the first two experiments (N220919 and N221204) the same diamond capsule batch was used, with the same quality and the measured yield increase as a function of change in hotspot oblateness (blue points) fell within the simulated range.

Post-shot simulations under-predicted the performance of N230729 given the expected degradation from the lower tungsten concentration when matching the measured hotspot P2 asymmetry and mode one hotspot velocity. To better match the observed performance, the level of P2 swing (which is currently unknown) was reduced in simulations. Upcoming tuning experiments will validate the time-varying P2 swing for this design. Using this model for the reduced swing that enables us to match performance for N230729 with the correct observed mode two and one asymmetries, we estimate the impact of the observed mode one.

Figure 5 (middle) shows the simulated fusion energy as a function of simulated hotspot velocity for α-off simulations and simulations that include α-deposition (“burn on” or α-on). These simulations include surface roughness, tent and fill tube perturbation features, and modes one, two, and four of the radiation drive. For higher yield and temperature implosions, the hotspot velocity is amplified as the implosion is burning on expansion, see Fig. 5 (right) for ratio of α-off/ α-on vs fusion energy. Matching the measured hotspot velocity in burn-on simulations gives a burn-off hotspot velocity (at maximum compression) of 70–75 km/s. This curve shows that reducing the mode one asymmetry could further to routinely achievable levels of < 60 km/s burn-on, increase the yield of N230729 by 30% (>5 MJ). The sensitivity of performance to mode one was benchmarked against the N210808 design, Fig. 5 (left), which experienced large mode one asymmetries from laser fluctuations in repeat experiments. The simulated impact of mode one on the performance for this design (red) was consistent with the measured impact (black points). Note, the mode one present for the N210808 design was also calculated to impact the performance by 25%.

FIG. 5.

(left) Neutron yield as a function of the mode one of the hot spot velocity (km/s) for the N210808 and repeat experiments, where the red points are normalized simulations and the black points are normalized measured values The experimental yields are normalized to N210808 with the normalization factor being the simulated reduction in yield for N210808 due to the observed hot spot velocity. (middle) Simulated fusion energy as a function of mode one velocity of the hot spot (km/s) for burn-on (α-on) and burn-off (α-off) simulations of N230729 which include surface roughness, tent and fill tube features, and low mode asymmetries (P1 (varying), P2, P4) to match the measured shape signatures and the performance of N230729. Here, an applied mode one of  0.7% reproduces the “Burn-on” measured hot spot velocity of >100 km/s. (right) Ratio of burn-on to burn-off hot spot velocity as a function of yield. For high fusion yields the burn-on hot spot velocity is increased by a larger fraction compared to the “burn-off” hot-spot velocity due to burning upon expansion.

FIG. 5.

(left) Neutron yield as a function of the mode one of the hot spot velocity (km/s) for the N210808 and repeat experiments, where the red points are normalized simulations and the black points are normalized measured values The experimental yields are normalized to N210808 with the normalization factor being the simulated reduction in yield for N210808 due to the observed hot spot velocity. (middle) Simulated fusion energy as a function of mode one velocity of the hot spot (km/s) for burn-on (α-on) and burn-off (α-off) simulations of N230729 which include surface roughness, tent and fill tube features, and low mode asymmetries (P1 (varying), P2, P4) to match the measured shape signatures and the performance of N230729. Here, an applied mode one of  0.7% reproduces the “Burn-on” measured hot spot velocity of >100 km/s. (right) Ratio of burn-on to burn-off hot spot velocity as a function of yield. For high fusion yields the burn-on hot spot velocity is increased by a larger fraction compared to the “burn-off” hot-spot velocity due to burning upon expansion.

Close modal

Experiment N230729 is also simulated to have generated more fusion energy with a corrected P2 that is less prolate. Figure 6 shows radiation hydrodynamic simulations of N230729 which match the level of performance (assuming the mode one was anomalous, and reducing it simulations) by picking a different level of P2 swing and matching the measured hotspot shape (left hand side). The simulations show contours at peak compression from α-off simulations of the hotspot temperature (left hand side) and density (right hand side). Reducing the level of prolate by decreasing the level of early-time waist-hot drive (a result of the increased wavelength separation compared to N221204) results in a 90% increase in performance with this correction alone, and >7MJ.

FIG. 6.

(left) Radiation Hydrodynamic simulations of N230729 which roughly captured the measured P2 and P4 symmetry and performance. (right) Reducing the level of mode two asymmetry by 5  μm results in a 90% increase in performance. The images include simulated contours at peak compression from α-off simulations of the hot spot temperature (left) and density (right) for example points.

FIG. 6.

(left) Radiation Hydrodynamic simulations of N230729 which roughly captured the measured P2 and P4 symmetry and performance. (right) Reducing the level of mode two asymmetry by 5  μm results in a 90% increase in performance. The images include simulated contours at peak compression from α-off simulations of the hot spot temperature (left) and density (right) for example points.

Close modal

Applying this model for the reduced swing, and accounting for as-shot differences in the P2 symmetry for N221204 and N230729, also gives a simulated performance for N221204 of >7MJ with higher dopant tungsten capsules and improved target quality, as well as lower residual kinetic energy (RKE). The difference between this prediction and the measured fusion energy of 3.15 MJ could be due to the reduced diamond capsule quality present for N221204. Experiments to study the impact of capsule quality vs lower tungsten dopant between N221204 and N230729 have bene proposed.

Following these experiments, further extending the laser pulse to 2.2 MJ of laser energy at the same peak power of 440TW (N231029 and N240210), resulted in a record fusion energy produced of 5.2 MJ, record target gain of 2.4×, and record fuel burn-up fraction of 7%. These experiments also increased the diamond ablator thickness by 5  μm compared to N230729. Ongoing work to improve the compressed shell symmetry in this platform, additionally increase the diamond ablator thickness, and increase the fuel burn-up fraction are ongoing.

The design presented in this paper provided more margin for achieving high fusion yields > 1 MJ in the presence of perturbations from non-ideal experimental conditions that reduced performance in the preceding N210808 design. These perturbations are predominantly unintentional low mode asymmetries from capsule non-uniformity or laser non-uniformity and reduced capsule quality that can seed hydrodynamic instabilities and mix ablator material into the hotspot. The hotspot energy density, a metric for ignition margin, was intentionally increased in this design by 20% through an increase in both the hotspot pressure and coupled energy.35 

Additional margin to these perturbations can be seen in Fig. 7 which shows the fusion energy produced for N210808 (black) and four experiments using the N221204 design (red). The solid bars denote experimental fielding conditions where the yield was not degraded as significantly by low mode asymmetries or high-Z mixing into the hotspot as the repeat experiments of N210808. The dashed red bars denote experiments in the N221204 design which experienced higher levels of low mode asymmetry (P1 or P2), degraded capsule quality (high-Z inclusions), or particles on the capsule. Note, the hotspot velocity being nearly the same for a significantly lower yield, N230618 vs N230729 (see Table II), represents a higher reduction in performance due to the fact that there is amplification of the hotspot velocity at higher yields. The hotspot velocity observed on N230618 was mainly due to non-uniformities in the initial capsule thickness and laser energy delivery, and the hotspot velocity observed on N230729 was mainly a result of an imbalance in laser energy delivery. Even with these perturbations, the N221204 design produced > 1 MJ of fusion energy, whereas the same perturbations applied to the N210808 design (black dashed bars) are simulated to have produced < MJ of fusion energy. The insets above the bars illustrate the main type of perturbation corresponding to that experiment. Table II shows variation in the experimental fielding conditions for the five experiments in the N221204 design shown in Fig. 7.

FIG. 7.

Fusion energy produced for the August 8, N210808 design (black) and the December 5, N221204 design (red). The crossed red bars denote fielded experiments in the N221204 design that suffered worse low mode asymmetries or capsule quality (particles on the capsule at shot-time). The crossed black bars denote the simulated yield for the N210808 design, hypothetically exposed to the same perturbations as the crossed red bars, which would have produced less energy than the N221204 design and <1 MJ of fusion yield.

FIG. 7.

Fusion energy produced for the August 8, N210808 design (black) and the December 5, N221204 design (red). The crossed red bars denote fielded experiments in the N221204 design that suffered worse low mode asymmetries or capsule quality (particles on the capsule at shot-time). The crossed black bars denote the simulated yield for the N210808 design, hypothetically exposed to the same perturbations as the crossed red bars, which would have produced less energy than the N221204 design and <1 MJ of fusion yield.

Close modal
TABLE II.

Variations in experimental fielding conditions for four experiments in the N221204 design. Larger values of hot spot velocity and P2 reduce the fusion yield by a greater amount.

N220919 N221204 N230618 N230729 N230904
Fusion energy (MJ)  1.22  3.15  1.6  3.88  2.09 
Hot spot velocity (km/s)  51  ± 9 (40 no- α 36  ± 6 (27 no- α 107  ± 7 (84 no- α 104  ± 9 (75 no- α 82  ± 17 (63 no- α
P2 PNIa (μm)  −10.5 (oblate)   0   0  +6 (prolate)  +6 (prolate) 
Capsule quality  >3000 high-Z inclusions  >3000 high-Z inclusions  <100 high-Z inclusions  <100 high-Z inclusions  <100 high-Z inclusions 
N220919 N221204 N230618 N230729 N230904
Fusion energy (MJ)  1.22  3.15  1.6  3.88  2.09 
Hot spot velocity (km/s)  51  ± 9 (40 no- α 36  ± 6 (27 no- α 107  ± 7 (84 no- α 104  ± 9 (75 no- α 82  ± 17 (63 no- α
P2 PNIa (μm)  −10.5 (oblate)   0   0  +6 (prolate)  +6 (prolate) 
Capsule quality  >3000 high-Z inclusions  >3000 high-Z inclusions  <100 high-Z inclusions  <100 high-Z inclusions  <100 high-Z inclusions 
a

PNI is the primary neutron image.

Ongoing work aims to further increase the compression or coupling of these designs to increase fuel-burn up fraction. This section presents two example methods to improve the performance of the N221204 design, in addition to the suggested higher capsule quality vs N221204 (high-W%  0.5% dopant) test and the mode two correction test of N230729.

Lowering the central gas fill density by lowering the “quench”50 temperature is predicted to increase convergence and performance of the N230729 design by more than 2×, see Fig. 8. These simulations include the “kitchen sink” (roughness, tent and fill-tube features, and low mode asymmetries) and match the performance of N230729 (left). Lowering the central gas fill density (right) from 0.000 43 to 0.0003 mg/cc increases convergence by  4.5%, the ignition margin by 20%, the total areal density by 10%, and the fusion yield by 2.35× (>9 MJ). These simulations also include differences in the DT gas composition with the change in initial temperature. Previous attempts to reduce the central fill density were performed in the N210808 (less stable) design with diamond capsules that had a significant number of defects, which resulted in lower than expected performance and complicated the understanding. Here we propose reducing the central gas fill density in the more stable N230729 design35 which has increased ignition margin and has shown less sensitivity to variations in diamond ablator capsule quality. This adjustment uses the same asymmetries that are present on N230729.

FIG. 8.

Radiation hydrodynamic simulations of the impact of reduced central DT vapor density at t = 0 by lowering the “quench” temperature. This reduces the central gas fill density from 0.000 43 to 0.0003 mg/cc which results in 4.5% higher convergence, 10% higher total areal density, and 1.2× increase in ignition margin (EP2). The performance is predicted to increase by 2.35x compared to N230729.

FIG. 8.

Radiation hydrodynamic simulations of the impact of reduced central DT vapor density at t = 0 by lowering the “quench” temperature. This reduces the central gas fill density from 0.000 43 to 0.0003 mg/cc which results in 4.5% higher convergence, 10% higher total areal density, and 1.2× increase in ignition margin (EP2). The performance is predicted to increase by 2.35x compared to N230729.

Close modal

Simulations suggest that the yield can also be substantially increased by further increasing the diamond thickness of the N221204 design using the same hohlraum configuration, laser energy, and adjusted shock timing. Once ignited, having a high total areal density of the fuel and ablator is advantageous for increasing yield amplification and fuel burn-up fraction. The first tests of the thicker ablator higher energy design took a substantial first step but started conservatively in thickness to balance the margin to ignite, increase in total areal density, and the symmetry risk of fielding thicker ablators. Exploring the limits of increase areal density could enhance robustness as well as performance (fusion energy produced and compression at a given yield) by several times, see Fig. 9 and Table III.

FIG. 9.

(left) Simulated neutron yield as a function of increased roughness of the diamond ablator interfaces and ice layer for future designs that further increase the diamond ablator thickness (green points) compared to N221204 and N230729 (red curve) by an additional 4, 6, 8, and 10  μm using the same 2.05 MJ laser energy and hohlraum geometry with adjusted shock timing. The black curves are for the N210808 design with (filled points) and without (open points) increased roughness on the DT ice layer. A roughness of 5× nominal matches the measured yield for N230729 without including low-mode asymmetries. For this level of roughness, the yield increases significantly by increasing ablator thickness. (right) Simulated increase in the Down-Scattered-Ratio (DSR) as a function of yield for the thicker diamond ablator simulations shown on the left hand side (green curve). The roll-back to higher DSR as yield increases is a result of burning closer to maximum compression vs the expansion phase.

FIG. 9.

(left) Simulated neutron yield as a function of increased roughness of the diamond ablator interfaces and ice layer for future designs that further increase the diamond ablator thickness (green points) compared to N221204 and N230729 (red curve) by an additional 4, 6, 8, and 10  μm using the same 2.05 MJ laser energy and hohlraum geometry with adjusted shock timing. The black curves are for the N210808 design with (filled points) and without (open points) increased roughness on the DT ice layer. A roughness of 5× nominal matches the measured yield for N230729 without including low-mode asymmetries. For this level of roughness, the yield increases significantly by increasing ablator thickness. (right) Simulated increase in the Down-Scattered-Ratio (DSR) as a function of yield for the thicker diamond ablator simulations shown on the left hand side (green curve). The roll-back to higher DSR as yield increases is a result of burning closer to maximum compression vs the expansion phase.

Close modal
TABLE III.

Key simulated implosion metrics for 2.05 MJ laser energy EL with an additional 4 and 10  μm thicker HDC ablator compared to baseline experiment N230729. Increasing the ablator thickness with the same laser driver has the potential to increase fusion yield and gain by more than a factor of three.

N230729 2.05 MJ + 4  μm HDC 2.05 MJ + 10  μm HDC
Implosion velocity (km/s)  386  383  374 
Ablator mass Remaining (%)  5.9  6.6 
Capsule absorbed energy (kJ)  245  257  274 
Norm. total ρRBT (g/cm2) α-off  1.0  1.08  1.16 
Norm. DSR (%) α-on  1.0  1.0  1.12 
Yield (MJ)  4.1  9.2  15.8 
N230729 2.05 MJ + 4  μm HDC 2.05 MJ + 10  μm HDC
Implosion velocity (km/s)  386  383  374 
Ablator mass Remaining (%)  5.9  6.6 
Capsule absorbed energy (kJ)  245  257  274 
Norm. total ρRBT (g/cm2) α-off  1.0  1.08  1.16 
Norm. DSR (%) α-on  1.0  1.0  1.12 
Yield (MJ)  4.1  9.2  15.8 

Figure 9 shows the simulated neutron yield as a function of increased roughness of the diamond ablator interfaces and ice layer with further increasing the diamond ablator thickness (green points) compared to N221204 and N230729 (red curve) by an additional 4, 6, 8, and 10  μm using the same 2.05 MJ laser energy and hohlraum geometry with adjusted shock timing. A roughness of 5× nominal was applied to match the measured yield for N230729 without including degradation from low-mode asymmetries. For this level of roughness, the yield increases significantly by increasing ablator thickness. While these simulations match yield by increasing interface roughness without including low mode asymmetries, additional studies including low mode asymmetries with nominal roughness also show the same level of impact for thicker ablators at these yield levels. The increase in compression at a given yield is also shown in Fig. 9 (right) via the Down-Scattered-Ratio (DSR) for the thicker diamond ablator simulations (green curve). The increase in DSR as yield increases for the green curve is a result of the thicker ablator designs burning closer to maximum compression in contrast to the expansion phase.

The design changes presented in this paper resulted in the first ever controlled fusion experiment (N221204) to exceed target gain greater > 1, where the fusion energy produced exceeded the laser energy used to drive the target by 1.54×. This design used a thicker diamond capsule (additional  7.6% in thickness) together with an extended higher energy laser drive (2.05 MJ of laser energy compared to 1.9 MJ)46,47 to increase the margin for ignition and enable higher DT fuel burn up fractions compared to N210808. A repeat experiment N230729 achieved the highest yield and gain to date for 2.05 MJ of laser driver energy (3.88 MJ and 1.9× target gain). Further extending the laser drive to 2.2 MJ and increasing the HDC ablator thickness by an additional  5  μm has achieved the maximum yield and target gain to date of 5.2 MJ and G  2.4×. We have also achieved target gain > 1 on five experiments, demonstrating robust and repeatable target gain. These proof-of-principle experiments demonstrate that there is nothing fundamentally limiting fusion energy gain in the laboratory.

The authors have no conflicts to disclose.

A. L. Kritcher: Conceptualization (equal); Investigation (equal); Software (equal); Supervision (equal); Writing – original draft (equal). D. J. Schlossberg: Data curation (equal); Supervision (equal). C. R. Weber: Methodology (supporting); Software (equal). C. V. Young: Formal analysis (supporting); Software (supporting). O. A. Hurricane: Conceptualization (equal); Supervision (supporting). E. Dewald: Data curation (supporting); Formal analysis (supporting). A. B. Zylstra: Data curation (lead); Supervision (equal). A. Allen: Data curation (supporting). B. Bachmann: Data curation (supporting). K. Baker: Data curation (supporting). S. Baxamusa: Data curation (supporting). T. Braun: Data curation (supporting). G. Brunton: Funding acquisition (supporting); Resources (supporting). D. A. Callahan: Conceptualization (supporting); Methodology (supporting); Software (supporting). D. T. Casey: Data curation (supporting); Formal analysis (supporting); Methodology (supporting). T. Chapman: Data curation (supporting). C. Choate: Data curation (supporting); Resources (equal). D. S. Clark: Formal analysis (supporting); Software (supporting). J.-M. G. Di Nicola: Resources (equal). L. Divol: Formal analysis (supporting). M. J. Edwards: Conceptualization (supporting); Project administration (supporting); Supervision (supporting). S. Haan: Formal analysis (supporting). T. Fehrenbach: Resources (equal). S. Hayes: Resources (equal). D. Hinkel: Methodology (equal); Supervision (supporting). M. Hohenberger: Data curation (supporting). K. Humbird: Software (supporting). N. Izumi: Data curation (supporting). O. Jones: Software (supporting). E. Kur: Software (supporting). B. Kustowski: Software (supporting). C. Kong: Resources (equal). O. L. Landen: Conceptualization (supporting); Data curation (supporting); Formal analysis (supporting); Investigation (supporting); Methodology (supporting). D. Larson: Resources (equal). X. Lepro-Chavez: Resources (equal). J. D. Lindl: Conceptualization (supporting); Formal analysis (supporting); Methodology (supporting). B. J. MacGowan: Formal analysis (supporting). S. MacLaren: Methodology (supporting). M. Marinak: Software (supporting). P. Michel: Formal analysis (supporting). M. Millot: Data curation (supporting); Formal analysis (supporting). A. Nikroo: Resources (supporting). R. Nora: Formal analysis (supporting). A. Pak: Data curation (supporting); Formal analysis (supporting). P. K. Patel: Formal analysis (supporting). J. E. Ralph: Data curation (supporting). M. Ratledge: Resources (supporting). M. S. Rubery: Data curation (supporting); Formal analysis (supporting). N. W. Ruof: Data curation (supporting); Formal analysis (supporting). S. M. Sepke: Software (supporting). M. Stadermann: Resources (equal). D. J. Strozzi: Formal analysis (supporting). T. I. Suratwala: Resources (equal). R. Tommasini: Data curation (supporting); Formal analysis (supporting). R. Town: Resources (equal); Supervision (supporting). Brandon Woodworth: Project administration (equal). B. Van Wonterghem: Resources (equal); Supervision (supporting). C. Wild: Resources (equal).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Radiation hydrodynamics simulations (HYDRA)38 were used to model the radiation drive created by the laser-hohlraum interaction as well as the physics of the capsule implosion. These simulations were performed in a two-step process: (1) integrated simulations of the hohlraum and capsule implosion which are benchmarked against experimental data from focused tuning experiments determine the spatially, temporally, and frequency resolved radiation drive surrounding the capsule implosion, and (2) the radiation drive is extracted and applied to higher resolution capsule-only kitchen sink simulations of the implosion to model engineering features with higher fidelity, such as material roughness of the capsule interfaces, a model for the capsule support tent, a model for the DT fill tube,51,52 non-uniformities in ablator thickness, and DT fuel layer thickness together with low mode non-uniformities of the radiation drive (including mode one, two, and four); see Fig. 10 for a comparison of the simulations to the experimental data.

FIG. 10.

(left) Schematic of the simulation methodology where the integrated simulations calibrated against experimental data determine the time-and spatially dependent radiation drive. The drive is extracted and applied to capsule-only simulations where the tent and fill-tube perturbation features, and roughness can be better modeled. (middle) Measured (black) and simulated (red) shock velocities for keyhole experiments of the N210808 and N221204 design, which were used (together with the time of maximum fusion production) to calibrate the radiation drive magnitude and symmetry. (right) simulations compared to experimental measurements of the inflight and hoy spot symmetry for N210808 which were used to calibrate the symmetry model.

FIG. 10.

(left) Schematic of the simulation methodology where the integrated simulations calibrated against experimental data determine the time-and spatially dependent radiation drive. The drive is extracted and applied to capsule-only simulations where the tent and fill-tube perturbation features, and roughness can be better modeled. (middle) Measured (black) and simulated (red) shock velocities for keyhole experiments of the N210808 and N221204 design, which were used (together with the time of maximum fusion production) to calibrate the radiation drive magnitude and symmetry. (right) simulations compared to experimental measurements of the inflight and hoy spot symmetry for N210808 which were used to calibrate the symmetry model.

Close modal

Simulations include detailed equations of state,53,54 radiation particle and neutron transport models,55,56 opacity models,57 and electron-ion coupling.58,59 The radiation drive is modeled with a flux-limited electron heat transport with a limiter of 0.15. An in-line cross beam energy transfer model was used with a saturation clamp on the electron density fluctuations of the ion acoustic wave (δn/n = 0.002 in until the start of the rise and 0.008 for the rest of the laser pulse). Post-shot simulations use the measured input target conditions and as-delivered laser powers vs time for all beams which can vary from the requested input conditions.

The simulations were adjusted in a common way between experiments to match tuning data and then applied to simulate the DT layered experiments. These adjustments include artificial multipliers on the input laser powers to match the in-flight symmetry, in-flight capsule velocity, and shock velocities along two lines of sight (pole and equator). Shock timing data were obtained immediately prior to N221204 which improved the symmetry model for temporal variations in the radiation drive. Back-scattered laser light due to the interaction of the lasers with the hohlraum wall plasma was energetically low on (N210808 and N221204) and mostly due to Stimulated Brillouin Scattering (SBS)49 of the inner beams early in peak laser power.

1.
E. I.
Moses
,
J.
Atherton
,
L.
Lagin
,
D.
Larson
,
C.
Keane
,
B.
MacGowan
,
R.
Patterson
,
M.
Spaeth
,
B.
Van Wonterghem
,
P.
Wegner
, and
R.
Kauffman
, “
The National Ignition Facility: Transition to a user facility
,”
J. Phys.: Conf. Ser.
688
,
012073
(
2016
).
2.
S. E.
Koonin
,
G. F.
Carrier
,
R.
Christy
,
R.
Conn
,
R.
Davidson
,
J.
Dawson
,
A.
DeMaria
,
P.
Doty
,
W.
Happer
, Jr.
,
G. L.
Kulcinski
et al, “
Review of the department of energy's inertial confinement fusion program
,”
J. Fusion Energy
10
,
157
(
1991
).
3.
R. E.
Kidder
, “
Laser fusion: The first ten years 1962–1972
,”
Report Nos. UCRL-BOOK-222681 and TRN: US0900468
(
Lawrence Livermore National Laboratory
,
Livermore
,
CA
,
2006
).
4.
J.
Nuckolls
,
L.
Wood
,
A.
Thiessen
, and
G.
Zimmerman
, “
Laser compression of matter to super-high densities: Thermonuclear (CTR) applications
,”
Nature
239
,
139
142
(
1972
).
5.
J. S.
Clarke
,
H. N.
Fisher
, and
R. J.
Mason
, “
Laser-driven implosion of spherical DT targets to thermonuclear burn conditions
,”
Phys. Rev. Lett.
30
,
249
249
(
1973
).
6.
J.
Lindl
, “
Development of the indirect-drive approach to inertial confinement fusion and the target physics basis for ignition and gain
,”
Phys. Plasmas
2
,
3933
4024
(
1995
).
7.
J. D.
Lindl
,
P.
Amendt
,
R. L.
Berger
,
S. G.
Glendinning
,
S. H.
Glenzer
,
S. W.
Haan
,
R. L.
Kauffman
,
O. L.
Landen
, and
L. J.
Suter
, “
The physics basis for ignition using indirect-drive targets on the National Ignition Facility
,”
Phys. Plasmas
11
,
339
491
(
2004
).
8.
S. W.
Haan
,
J. D.
Lindl
,
D. A.
Callahan
,
D. S.
Clark
,
J. D.
Salmonson
,
B. A.
Hammel
,
L. J.
Atherton
,
R. C.
Cook
,
M. J.
Edwards
,
S. H.
Glenzer
et al, “
Point design targets, specifications, and requirements for the 2010 ignition campaign on the National Ignition Facility
,”
Phys. Plasmas
18
,
051001
(
2011
).
9.
R. F.
Post
, “
Controlled fusion research—An application of the physics of high temperature plasmas
,”
Rev. Mod. Phys.
28
,
338
362
(
1956
).
10.
P.
Volegov
,
C. R.
Danly
,
D. N.
Fittinghoff
,
G. P.
Grim
,
N.
Guler
,
N.
Izumi
,
T.
Ma
,
F. E.
Merrill
,
A. L.
Warrick
,
C. H.
Wilde
, and
D. C.
Wilson
,
Rev. Sci. Instrum.
85
,
023508
(
2014
).
11.
S. R.
Nagel
,
S. W.
Haan
,
J. R.
Rygg
,
M.
Barrios
,
L. R.
Benedetti
,
D. K.
Bradley
,
J. E.
Field
,
B. A.
Hammel
,
N.
Izumi
,
O. S.
Jones
et al, “
Effect of the mounting membrane on shape in inertial confinement fusion implosions
,”
Phys. Plasmas
22
,
022704
(
2015
).
12.
T. R.
Dittrich
,
O. A.
Hurricane
,
D. A.
Callahan
,
E. L.
Dewald
,
T.
Döppner
,
D. E.
Hinkel
,
L. F. B.
Hopkins
,
S.
Le Pape
,
T.
Ma
,
J. L.
Milovich
et al, “
Design of a high-foot high-adiabat ICF capsule for the National Ignition Facility
,”
Phys. Rev. Lett.
112
(
5
),
055002
(
2014
).
13.
S.
Le Pape
,
L. F.
Berzak Hopkins
,
L.
Divol
,
A.
Pak
,
E. L.
Dewald
,
S.
Bhandarkar
,
L. R.
Bennedetti
,
T.
Bunn
,
J.
Biener
,
J.
Crippen
et al, “
Fusion energy output greater than the kinetic energy of an imploding shell at the National Ignition Facility
,”
Phys. Rev. Lett.
120
,
245003
(
2018
).
14.
L. B.
Hopkins
,
S.
LePape
,
L.
Divol
,
A.
Pak
,
E.
Dewald
,
D. D.
Ho
,
N.
Meezan
,
S.
Bhandarkar
,
L. R.
Benedett
,
T.
Bunn
et al, “
Toward a burning plasma state using diamond ablator Inertially Confined Fusion (ICF) implosions on the National Ignition Facility (NIF)
,”
Plasma Phys. Controlled Fusion
61
,
014023
(
2018
).
15.
K. L.
Baker
,
C. A.
Thomas
,
D. T.
Casey
,
M.
Hohenberger
,
S.
Khan
,
B. K.
Spears
,
O. L.
Landen
,
R.
Nora
,
D. T.
Woods
,
J. L.
Milovich
et al, “
Hotspot parameter scaling with velocity and yield for high-adiabat layered implosions at the National Ignition Facility
,”
Phys. Rev. E
102
,
023210
(
2020
).
16.
D. T.
Casey
,
C.
Thomas
,
K. L.
Baker
,
B. K.
Spears
,
M.
Hohenberger
,
S. F.
Khan
,
R. C.
Nora
,
C. R.
Weber
,
D. T.
Woods
,
O.
Hurricane
et al,
Phys. Plasmas
25
,
056308
(
2018
).
17.
M.
Hohenberger
,
D. T.
Casey
,
A. L.
Kritcher
,
A.
Pak
,
A. B.
Zylstra
,
C. A.
Thomas
,
K. L.
Baker
,
S. L.
Pape
,
B.
Bachmann
,
R. L.
Berger
et al, “
Integrated performance of large HDC-capsule implosions on the National Ignition Facility
,”
Phys. Plasmas
27
,
112704
(
2020
).
18.
P.
Michel
,
L.
Divol
,
E. A.
Williams
,
S.
Weber
,
C. A.
Thomas
,
D. A.
Callahan
,
S. W.
Haan
,
J. D.
Salmonson
,
S.
Dixit
,
D. E.
Hinkel
et al, “
Tuning the implosion symmetry of ICF targets via controlled crossed-beam energy transfer
,”
Phys. Rev. Lett.
102
,
025004
(
2009
).
19.
S. L.
Pape
,
L. F. B.
Hopkins
,
L.
Divol
,
N.
Meezan
,
D.
Turnbull
,
A. J.
Mackinnon
,
D.
Ho
,
J. S.
Ross
,
S.
Khan
,
A.
Pak
et al, “
The near vacuum hohlraum campaign at the NIF: A new approach
,”
Phys. Plasmas
23
,
056311
(
2016
).
20.
L. F.
Berzak Hopkins
,
N. B.
Meezan
,
S.
Le Pape
,
L.
Divol
,
A. J.
Mackinnon
,
D. D.
Ho
,
M.
Hohenberger
,
O. S.
Jones
,
G.
Kyrala
,
J. L.
Milovich
et al, “
First high-convergence cryogenic implosion in a near-vacuum hohlraum
,”
Phys. Rev. Lett.
114
,
175001
(
2015
).
21.
L.
Divol
,
A.
Pak
,
L. F. B.
Hopkins
,
S. L.
Pape
,
N. B.
Meezan
,
E. L.
Dewald
,
D. D.-M.
Ho
,
S. F.
Khan
,
A. J.
Mackinnon
,
J. S.
Ross
et al, “
Symmetry control of an indirectly driven high-density-carbon implosion at high convergence and high velocity
,”
Phys. Plasmas
24
,
056309
(
2017
).
22.
J. D.
Moody
,
P.
Michel
,
L.
Divol
,
R. L.
Berger
,
E.
Bond
,
D. K.
Bradley
,
D. A.
Callahan
,
E. L.
Dewald
,
S.
Dixit
,
M. J.
Edwards
et al, “
Multistep redirection by cross-beam power transfer of ultrahigh-power lasers in a plasma
,”
Nat. Phys.
8
,
344
349
(
2012
).
23.
A. L.
Kritcher
,
R.
Town
,
D.
Bradley
,
D.
Clark
,
B.
Spears
,
O.
Jones
,
S.
Haan
,
P. T.
Springer
,
J.
Lindl
,
R. H. H.
Scott
et al, “
Metrics for long wavelength asymmetries in inertial confinement fusion implosions on the National Ignition Facility
,”
Phys. Plasmas
21
,
042708
(
2014
).
24.
O. A.
Hurricane
,
D. T.
Casey
,
O.
Landen
,
A. L.
Kritcher
,
R.
Nora
,
P. K.
Patel
,
J. A.
Gaffney
,
K. D.
Humbird
,
J. E.
Field
,
M. K. G.
Kruse
,
J. L.
Peterson
, and
B. K.
Spears
, “
An analytic asymmetric-piston model for the impact of mode-1 shell asymmetry on ICF implosions
,”
Phys. Plasmas
27
,
062704
(
2020
).
25.
B.
MacGowan
,
O.
Landen
,
D.
Casey
,
C.
Young
,
D.
Callahan
,
E.
Hartouni
,
R.
Hatarik
,
M.
Hohenberger
,
T.
Ma
,
D.
Mariscal
et al, “
Trending low mode asymmetries in NIF capsule drive using a simple viewfactor metric
,”
High Energy Density Phys.
40
,
100944
(
2021
).
26.
D. T.
Casey
,
B. J.
MacGowan
,
J. D.
Sater
,
A. B.
Zylstra
,
O. L.
Landen
,
J.
Milovich
,
O. A.
Hurricane
,
A. L.
Kritcher
,
M.
Hohenberger
,
K.
Baker
et al, “
Evidence of three-dimensional asymmetries seeded by high-density carbon-ablator nonuniformity in experiments at the National Ignition Facility
,”
Phys. Rev. Lett.
126
,
025002
(
2021
).
27.
A. L.
Kritcher
,
D. E.
Hinkel
,
D. A.
Callahan
,
O. A.
Hurricane
,
D.
Clark
,
D. T.
Casey
,
E. L.
Dewald
,
T. R.
Dittrich
,
T.
Döppner
,
M. A. B.
Garcia
,
S.
Haan
,
L. F. B.
Hopkins
,
O.
Jones
,
O.
Landen
,
T.
Ma
,
N.
Meezan
,
J. L.
Milovich
,
A. E.
Pak
,
H. S.
Park
,
P. K.
Patel
,
J.
Ralph
,
H. F.
Robey
,
J. D.
Salmonson
,
S.
Sepke
,
B.
Spears
,
P. T.
Springer
,
C. A.
Thomas
,
R.
Town
,
P. M.
Celliers
, and
M. J.
Edwards
, “
Integrated modeling of cryogenic layered highfoot experiments at the NIF
,”
Phys. Plasmas
23
,
052709
(
2016
).
28.
A. L.
Kritcher
,
J.
Ralph
,
D. E.
Hinkel
,
T.
Döppner
,
M.
Millot
,
D.
Mariscal
,
R.
Benedetti
,
D. J.
Strozzi
,
T.
Chapman
,
C.
Goyon
et al, “
Energy transfer between lasers in low-gas-fill-density hohlraums
,”
Phys. Rev. E
98
,
053206
(
2018
).
29.
L. A.
Pickworth
,
T.
Doeppner
,
D. E.
Hinkel
,
J. E.
Ralph
,
B.
Bachmann
,
L. P.
Masse
,
L.
Divol
,
L. R.
Benedetti
,
P. M.
Celliers
,
H.
Chen
et al,
Phys. Plasmas
27
,
102702
(
2020
).
30.
A. B.
Zylstra
,
A. L.
Kritcher
,
O. A.
Hurricane
,
D. A.
Callahan
,
K.
Baker
,
T.
Braun
,
D. T.
Casey
,
D.
Clark
,
K.
Clark
,
T.
Döppner
et al, “
Record energetics for an inertial fusion implosion at NIF
,”
Phys. Rev. Lett.
126
,
025001
(
2021
).
31.
A. L.
Kritcher
,
A. B.
Zylstra
,
D. A.
Callahan
,
O. A.
Hurricane
,
C.
Weber
,
J.
Ralph
,
D. T.
Casey
,
A.
Pak
,
K.
Baker
,
B.
Bachmann
et al, “
Achieving record hot spot energies with large HDC implosions on NIF in HYBRID-E
,”
Phys. Plasmas
28
,
072706
(
2021
).
32.
A.
Kritcher
,
C.
Young
,
H.
Robey
,
C.
Weber
,
A.
Zylstra
,
O.
Hurricane
,
D.
Callahan
,
J.
Ralph
,
J.
Ross
,
K.
Baker
et al, “
Design of implosions reaching the burning plasma regime
,”
Nat. Phys.
18
,
251
(
2022
).
33.
A.
Zylstra
,
O.
Hurricane
,
D.
Callahan
,
A.
Kritcher
,
J.
Ralph
,
H. F.
Robey
,
J. S.
Ross
,
C. V.
Young
,
K.
Baker
,
D. T.
Casey
et al, “
Burning plasma achieved in inertial fusion
,”
Nature
601
,
542
548
(
2022
).
34.
A. L.
Kritcher
,
A. B.
Zylstra
,
D. A.
Callahan
,
O. A.
Hurricane
,
C. R.
Weber
,
D. S.
Clark
,
C. V.
Young
,
J. E.
Ralph
,
D. T.
Casey
,
A.
Pak
,
O. L.
Landen
et al,
Phys. Rev. E
106
,
025201
(
2022
).
35.
A. L.
Kritcher
,
A. B.
Zylstra
,
C. R.
Weber
,
O. A.
Hurricane
,
D. A.
Callahan
,
D. S.
Clark
,
L.
Divol
,
D. E.
Hinkel
,
K.
Humbird
,
O.
Jones
,
J. D.
Lindl
et al, “
Design of the first fusion experiment to achieve target energy gain G > 1
,”
Phys. Rev. E
109
,
025204
(
2024
).
36.
O. A.
Hurricane
,
P. T.
Springer
,
P. K.
Patel
,
D. A.
Callahan
,
K.
Baker
,
D. T.
Casey
,
L.
Divol
,
T.
Döppner
,
D. E.
Hinkel
,
M.
Hohenberger
et al, “
Approaching a burning plasma on the NIF
,”
Phys. Plasmas
26
,
052704
(
2019
).
37.
D. A.
Callahan
,
O. A.
Hurricane
,
J. E.
Ralph
,
C. A.
Thomas
,
K. L.
Baker
,
L. R.
Benedetti
,
L. F.
Berzak Hopkins
,
D. T.
Casey
,
T.
Chapman
,
C. E.
Czajka
et al, “
Exploring the limits of case-to-capsule ratio, pulse length, and picket energy for symmetric hohlraum drive on the National Ignition Facility Laser
,”
Phys. Plasmas
25
,
5
(
2018
).
38.
M. M.
Marinak
,
G. D.
Kerbel
,
N. A.
Gentile
,
O.
Jones
,
D.
Munro
,
S.
Pollaine
,
T. R.
Dittrich
, and
S. W.
Haan
, “
Three-dimensional hydra simulations of National Ignition Facility targets
,”
Phys. Plasmas
8
,
2275
2280
(
2001
).
39.
O. A.
Hurricane
,
A.
Kritcher
,
D. A.
Callahan
,
O.
Landen
,
P. K.
Patel
,
P. T.
Springer
,
D. T.
Casey
,
E. L.
Dewald
,
T. R.
Dittrich
,
T.
Döppner
et al, “
On the importance of minimizing ‘coast-time’ in x-ray driven inertially confined fusion implosions
,”
Phys. Plasmas
24
,
092706
(
2017
).
40.
J. D.
Lawson
, “
Some criteria for a power producing thermonuclear reactor
,”
Proc. Phys. Soc. Sect. B
70
,
6
10
(
1957
).
41.
H.
Abu-Shawareb
et al, “
Lawson criterion for ignition achieved in an inertial fusion experiment
,”
Phys. Rev. Lett.
129
,
075001
(
2022
).
42.
A.
Zylstra
,
A. L.
Kritcher
,
O. A.
Hurricane
,
D. A.
Callahan
,
J. E.
Ralph
,
D. T.
Casey
,
A.
Pak
,
O. L.
Landen
,
B.
Bachmann
,
K. L.
Baker
et al,
Phys. Rev. E
106
,
025202
(
2022
).
43.
A.
Pak
,
A.
Zylstra
,
K. L.
Baker
,
D. T.
Casey
,
E.
Dewald
,
L.
Divol
,
M.
Hohenberger
,
A. S.
Moore
,
J. E.
Ralph
,
D. J.
Schlossberg
et al, “
Observations and properties of the first laboratory fusion experiment to exceed target gain of 1
,”
Phys. Rev. E
109
,
025203
(
2024
).
44.
H.
Abu-Shawareb
et al, “
Achievement of target gain larger than unity in an inertial fusion experiment
,”
Phys. Rev. Lett.
132
,
065102
(
2024
).
45.
O. A.
Hurricane
,
D. A.
Callahan
,
D. T.
Casey
,
A. R.
Christopherson
,
A. L.
Kritcher
,
O. L.
Landen
,
S. A.
Maclaren
,
R.
Nora
,
P.
Patel
,
J.
Ralph
,
D.
Schlossberg
,
P. T.
Springer
,
C. V.
Young
, and
A. B.
Zylstra
, “
Energy principles of scientific breakeven in an inertial fusion experiment
,”
Phys. Rev. Lett.
132
,
065103
(
2024
).
46.
J.-M. G.
Di Nicola
, “
First demonstration of fusion ignition by inertial confinement fusion (ICF) at the National Ignition Facility (NIF) at LLNL
,”
Proc. SPIE
PC12401
,
PC1240101
(
2023
).
47.
J. D.
Nicola
,
T.
Bond
,
M.
Bowers
,
L.
Chang
,
M.
Hermann
,
R.
House
,
T.
Lewis
,
K.
Manes
,
G.
Mennerat
,
B.
MacGowan
,
R.
Negres
,
B.
Olejniczak
,
C.
Orth
,
T.
Parham
,
S.
Rana
,
B.
Raymond
,
M.
Rever
,
S.
Schrauth
,
M.
Shaw
,
M.
Spaeth
,
B. V.
Wonterghem
,
W.
Williams
,
C.
Widmayer
,
S.
Yang
,
P.
Whitman
, and
P.
Wegner
, “
The National Ignition Facility: Laser performance status and performance quad results at elevated energy
,”
Nucl. Fusion
59
,
032004
(
2018
).
48.
D. J.
Strozzi
,
D. S.
Bailey
,
P.
Michel
,
L.
Divol
,
S. M.
Sepke
,
G. D.
Kerbel
,
C. A.
Thomas
,
J. E.
Ralph
,
J. D.
Moody
, and
M. B.
Schneider
, “
Interplay of laser-plasma interactions and inertial fusion hydrodynamics
,”
Phys. Rev. Lett.
118
,
025002
(
2017
).
49.
R.
Kirkwood
,
D.
Strozzi
,
P.
Michel
,
D.
Callahan
,
B.
Raymond
,
G.
Gururangan
,
B.
MacGowan
,
NIF Team
, “
Laser backscatter damage risk assessments of NIF target experiments
,” in APS Meeting Abstracts,
2014
.
50.
B. J.
Kozioziemski
,
E. R.
Mapoles
,
J. D.
Sater
,
A. A.
Chernov
,
J. D.
Moody
,
J. B.
Lugten
, and
M. A.
Johnson
, “
Deuterium-tritium fuel layer formation for the National Ignition Facility
,”
Fusion Sci. Technol.
59
(
1
),
14
25
(
2011
).
51.
D. S.
Clark
,
C. R.
Weber
,
J. L.
Milovich
,
J. D.
Salmonson
,
A. L.
Kritcher
,
S. W.
Haan
,
B. A.
Hammel
,
D. E.
Hinkel
,
O. A.
Hurricane
,
O. S.
Jones
,
M. M.
Marinak
,
P. K.
Patel
,
H. F.
Robey
,
S. M.
Sepke
,
M. J.
Edwards
et al, “
Three-dimensional simulations of low foot and high foot implosion experiments on the National Ignition Facility
,”
Phys. Plasmas
23
,
056302
(
2016
).
52.
D. S.
Clark
,
C. R.
Weber
,
J. L.
Milovich
,
A. E.
Pak
,
D. T.
Casey
,
B. A.
Hammel
,
D. D.
Ho
,
O. S.
Jones
,
J. M.
Koning
,
A. L.
Kritcher
et al, “
Three-dimensional modeling and hydrodynamic scaling of National Ignition Facility implosions
,”
Phys. Plasmas
26
,
050601
(
2019
).
53.
L. X. L.
Benedict
,
K. P.
Driver
,
S.
Hamel
,
B.
Militzer
,
T.
Qi
,
A. A.
Correa
,
A.
Saul
, and
E.
Schwegler
, “
Multiphase equation of state for carbon addressing high pressures and temperatures
,”
Phys. Rev. B
89
,
224109
(
2014
).
54.
J.
Gaffney
,
S.
Hu
,
P.
Arnault
,
A.
Becker
,
L.
Benedict
,
T.
Boehly
,
P.
Celliers
,
D.
Ceperley
,
O.
Certik
,
J.
Clerouin
et al, “
A review of equation-of-state models for inertial confinement fusion materials
,”
High Energy Density Phys.
28
,
7
24
(
2018
).
55.
M. P.
Desjarlais
,
C. R.
Scullard
,
L. X.
Benedict
,
H. D.
Whitley
, and
R.
Redmer
, “
Density-functional calculations of transport properties in the nondegenerate limit and the role of electron-electron scattering
,”
Phys. Rev. E
95
,
033203
(
2017
).
56.
L. G.
Stanton
and
M. S.
Murillo
, “
Ionic transport in high-energy-density matter
,”
Phys. Rev. E
93
,
043203
(
2016
).
57.
C. A.
Iglesias
and
F. J.
Rogers
, “
Updated opal opacities
,”
Astrophys. J.
464
,
943
(
1996
).
58.
Y.
Lee
and
R.
More
, “
An electron conductivity model for dense plasmas
,”
Phys. Fluids
27
,
1273
1286
(
1984
).
59.
C. R.
Scullard
,
S.
Serna
,
L. X.
Benedict
,
C. L.
Ellison
, and
F. R.
Graziani
, “
Analytic expressions for electron-ion temperature equilibration rates from the Lenard-Balescu equation
,”
Phys. Rev. E
97
,
013205
(
2018
).