Achievement of Target Gain Larger than Unity in an Inertial Fusion Experiment

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.

Neutron time of flight (nToF) detectors for inertial fusion experiments

Neutrons generated in Inertial Confinement Fusion (ICF) experiments provide valuable information to interpret the conditions reached in the plasma. The neutron time-of-flight (nToF) technique is well suited for measuring the neutron energy spectrum due to the short time (100 ps) over which neutrons are typically emitted in ICF experiments. By locating detectors 10s of meters from the source, the neutron energy spectrum can be measured to high precision. We present a contextual review of the current state of the art in nToF detectors at ICF facilities in the United States, outlining the physics that can be measured, the detector technologies currently deployed and analysis techniques used.

Alpha heating of indirect-drive layered implosions on the National Ignition Facility

In order to understand how close current layered implosions in indirect-drive inertial confinement fusion are to ignition, it is necessary to measure the level of alpha heating present. To this end, pairs of experiments were performed that consisted of a low-yield tritium-hydrogen-deuterium (THD) layered implosion and a high-yield deuterium-tritium (DT) layered implosion to validate experimentally current simulation-based methods of determining yield amplification. The THD capsules were designed to reduce simultaneously DT neutron yield (alpha heating) and maintain hydrodynamic similarity with the higher yield DT capsules. The ratio of the yields measured in these experiments then allowed the alpha heating level of the DT layered implosions to be determined. The level of alpha heating inferred is consistent with fits to simulations expressed in terms of experimentally measurable quantities and enables us to infer the level of alpha heating in recent high-performing implosions.

Design of a multi-detector, single line-of-sight, time-of-flight system to measure time-resolved neutron energy spectra

In the dynamic environment of burning, thermonuclear deuterium-tritium plasmas, diagnosing the time-resolved neutron energy spectrum is of critical importance. Strategies exist for this diagnosis in magnetic confinement fusion plasmas, which presently have a lifetime of ∼10¹² longer than inertial confinement fusion (ICF) plasmas. Here, we present a novel concept for a simple, precise, and scale-able diagnostic to measure time-resolved neutron spectra in ICF plasmas. The concept leverages general tomographic reconstruction techniques adapted to time-of-flight parameter space, and then employs an updated Monte Carlo algorithm and National Ignition Facility-relevant constraints to reconstruct the time-evolving neutron energy spectrum. Reconstructed spectra of the primary 14.028 MeV n_DT peak are in good agreement with the exact synthetic spectra. The technique is also used to reconstruct the time-evolving downscattered spectrum, although the present implementation shows significantly more error.

Constraining time-dependent ion temperature measurements in inertial confinement fusion (ICF) implosions with an intermediate distance neutron time-of-flight (nToF) detector

A concept for using an intermediate distance (0.3-3.0 m) neutron time-of-flight (nToF) to provide a constraint on the measurement of the time-dependence of ion temperature in inertial confinement fusion implosions is presented. Simulated nToF signals at different distances are generated and, with a priori knowledge of the burn-averaged quantities and burn history, analyzed to determine requirements for a future detector. Results indicate a signal-to-noise ratio >50 and time resolution <20 ps to constrain the ion temperature gradient to ∼±25% (0.5 keV/100 ps).

Construction and study of instrument response functions for analysis of the National Ignition Facility (NIF) neutron time-of-flight detectors

The analysis of the National Ignition Facility (NIF) neutron time-of-flight (nToF) detectors uses a forward-fit routine that depends critically on the instrument response functions (IRFs) of the diagnostics. The details of the IRFs used can have large impacts on measurements such as ion temperature and down-scattered ratio (DSR). Here, we report on the recent steps taken to construct and validate nToF IRFs at the NIF to an increased degree of accuracy, as well as remove the need for fixed DSR baseline offsets. The IRF is treated in two parts: a "core," measured experimentally with an x-ray impulse source, and a "tail" that occurs later in time and has limited experimental data. The tail region is calibrated with the data from indirect drive exploding pusher shots, which have little neutron scattering and are traditionally assumed to have zero DSR. Using analytic modeling estimates, the non-zero DSR for these shots is estimated. The impact of varying IRF tail components on DSR is investigated with a systematic parameter study, and good agreement is found with the non-zero DSR estimates. These approaches will be used to improve the precision and uncertainty of NIF nToF DSR measurements.

Lawson Criterion for Ignition Exceeded in an Inertial Fusion Experiment

For more than half a century, researchers around the world have been engaged in attempts to achieve fusion ignition as a proof of principle of various fusion concepts. Following the Lawson criterion, an ignited plasma is one where the fusion heating power is high enough to overcome all the physical processes that cool the fusion plasma, creating a positive thermodynamic feedback loop with rapidly increasing temperature. In inertially confined fusion, ignition is a state where the fusion plasma can begin "burn propagation" into surrounding cold fuel, enabling the possibility of high energy gain. While "scientific breakeven" (i.e., unity target gain) has not yet been achieved (here target gain is 0.72, 1.37 MJ of fusion for 1.92 MJ of laser energy), this Letter reports the first controlled fusion experiment, using laser indirect drive, on the National Ignition Facility to produce capsule gain (here 5.8) and reach ignition by nine different formulations of the Lawson criterion.

Observation of Hydrodynamic Flows in Imploding Fusion Plasmas on the National Ignition Facility

Inertial confinement fusion implosions designed to have minimal fluid motion at peak compression often show significant linear flows in the laboratory, attributable per simulations to percent-level imbalances in the laser drive illumination symmetry. We present experimental results which intentionally varied the mode 1 drive imbalance by up to 4% to test hydrodynamic predictions of flows and the resultant imploded core asymmetries and performance, as measured by a combination of DT neutron spectroscopy and high-resolution x-ray core imaging. Neutron yields decrease by up to 50%, and anisotropic neutron Doppler broadening increases by 20%, in agreement with simulations. Furthermore, a tracer jet from the capsule fill-tube perturbation that is entrained by the hot-spot flow confirms the average flow speeds deduced from neutron spectroscopy.

Three-dimensional diagnostics and measurements of inertial confinement fusion plasmas

Recent inertial confinement fusion measurements have highlighted the importance of 3D asymmetry effects on implosion performance. One prominent example is the bulk drift velocity of the deuterium-tritium plasma undergoing fusion ("hotspot"), v_HS. Upgrades to the National Ignition Facility neutron time-of-flight diagnostics now provide v_HS to better than 1 part in 10⁴ and enable cross correlations with other measurements. This work presents the impact of v_HS on the neutron yield, downscatter ratio, apparent ion temperature, electron temperature, and 2D x-ray emission. The necessary improvements to diagnostic suites to take these measurements are also detailed. The benefits of using cross-diagnostic analysis to test hotspot models and theory are discussed, and cross-shot trends are shown.

Understanding the effects of neutron scattering for neutron-yield-isotropy measurements at the NIF

Neutron-yield diagnostics at the NIF have been upgraded to include 48 detectors placed around the NIF target chamber to assess the DT-neutron-yield isotropy for inertial confinement fusion experiments. Real-time neutron-activation detectors are used to understand yield asymmetries due to Doppler shifts in the neutron energy attributed to hotspot motion, variations in the fuel and ablator areal densities, and other physics effects. In order to isolate target physics effects, we must understand the contribution due to neutron scattering associated with the different hardware configurations used for each experiment. We present results from several calibration experiments that demonstrate the ability to achieve our goal of 1% or better precision in determining the yield isotropy.

Interpolating individual line-of-sight neutron spectrometer measurements onto the "sky" at the National Ignition Facility (NIF)

Nuclear diagnostics provide measurements of inertial confinement fusion implosions used as metrics of performance for the shot. The interpretation of these measurements for shots with low mode asymmetries requires a way of combining the data to produce a "sky map" where the individual line-of-sight values are used to interpolate to other positions in the sky. These interpolations can provide information regarding the orientation of the low mode asymmetries. We describe the interpolation method, associated uncertainties, and correlations between different metrics, e.g., T_ion, down scatter ratio, and hot-spot velocity direction. This work is also related to recently reported studies [H. G. Rinderknecht et al., Phys. Rev. Lett. 124, 145002 (2020) and K. M. Woo et al., Phys. Plasmas 27, 062702 (2020)] of low mode asymmetries. We report an analysis that makes use of a newly commissioned line of sight, a scheme for incorporating multiple neutron spectrum measurement types, and recent work on the sources of implosion asymmetry to provide a more complete picture of implosion performance.

Real-time nuclear activation detectors for measuring neutron angular distributions at the National Ignition Facility (invited)

The Real Time Nuclear Activation Detector (RTNAD) array at NIF measures the distribution of 14 MeV neutrons emitted by deuterium-tritium (DT) fueled inertial confinement fusion implosions. The uniformity of the neutron distribution is an important indication of implosion symmetry and DT shell integrity. The array consists of 48 LaBr₃(Ce) crystal gamma-ray spectrometers mounted outside the NIF target chamber, which continuously monitor the slow decay of the 909 keV gamma-ray line from activated ⁸⁹Zr located in Zr cups surrounding each crystal. The measured decay rate dramatically increases during a DT implosion in proportion to the number of 14 MeV neutrons striking each Zr cup. The neutrons produce activated ⁸⁹Zr through an (n, 2n) reaction on ⁹⁰Zr, which is insensitive to low energy neutrons. The neutron flux along the detector line-of-sight at shot time is determined by extrapolating the fitted 909 keV decay curve back to shot time. Automatic analysis algorithms were developed to handle the non-stop data stream. The large number of detectors and the high statistical accuracy of the array enable the spherical harmonic modes of the neutron angular distribution to be measured up to L ≤ 4 to provide a better understanding of implosion dynamics. In addition, these data combined with measurements of the down-scattered neutrons can be used to derive fuel areal density distributions. This paper will describe the RTNAD hardware and analysis procedures.

Proof-of-concept of a neutron time-of-flight ellipsoidal detector

The time-resolved measurement of neutrons emitted from nuclear implosions at inertial confinement fusion facilities is used to characterize the fusing plasma. Several significant quantities are routinely measured by neutron time-of-flight (nToF) detectors in these experiments. Current nToF detectors use scintillators as well as solid-state Cherenkov radiators. The latter has an inherently faster time response and can provide a co-registered γ-ray measurement as well as improved precision in the bulk hot-spot velocity. This work discusses a nToF ellipsoidal detector that also utilizes a solid-state Cherenkov radiator. The detector has the potential to achieve a fast instrument response function allowing for characterization of the γ-ray burn history as well as the ability to field the detector closer to the fusion source. Proof-of-concept testing of the nToF ellipsoidal detector has been conducted at the National Ignition Facility using commercial optics. A time-resolved neutron signal has been measured from the diagnostic. Preliminary simulations corroborate the results.

Optimal choice of multiple line-of-sight measurements determining plasma hotspot velocity at the National Ignition Facility

The measurement of plasma hotspot velocity provides an important diagnostic of implosion performance for inertial confinement fusion experiments at the National Ignition Facility. The shift of the fusion product neutron mean kinetic energy as measured along multiple line-of-sight time-of-flight spectrometers provides velocity vector components from which the hotspot velocity is inferred. Multiple measurements improve the hotspot velocity inference; however, practical considerations of available space, operational overhead, and instrumentation costs limit the number of possible line-of-sight measurements. We propose a solution to this classical "experiment design" problem that optimizes the precision of the velocity inference for a limited number of measurements.

The five line-of-sight neutron time-of-flight (nToF) suite on the National Ignition Facility (NIF)

Measurement of the neutron spectrum from inertial confinement fusion implosions is one of the primary diagnostics of implosion performance. Analysis of the spectrum gives access to quantities such as neutron yield, hot-spot velocity, apparent ion temperature, and compressed fuel ρr through measurement of the down-scatter ratio. On the National Ignition Facility, the neutron time-of-flight suite has been upgraded to include five independent, collimated lines of sight, each comprising a high dynamic range bibenzyl/diphenylacetylene-stilbene scintillator [R. Hatarik et al., Plasma Fusion Res. 9, 4404104 (2014)] and high-speed fused silica Cherenkov detectors [A. S. Moore et al., Rev. Sci. Instrum. 89, 10I120 (2018)].

Time-Resolved Fuel Density Profiles of the Stagnation Phase of Indirect-Drive Inertial Confinement Implosions

The implosion efficiency in inertial confinement fusion depends on the degree of stagnated fuel compression, density uniformity, sphericity, and minimum residual kinetic energy achieved. Compton scattering-mediated 50-200 keV x-ray radiographs of indirect-drive cryogenic implosions at the National Ignition Facility capture the dynamic evolution of the fuel as it goes through peak compression, revealing low-mode 3D nonuniformities and thicker fuel with lower peak density than simulated. By differencing two radiographs taken at different times during the same implosion, we also measure the residual kinetic energy not transferred to the hot spot and quantify its impact on the implosion performance.

Azimuthal Drive Asymmetry in Inertial Confinement Fusion Implosions on the National Ignition Facility

Data from nuclear diagnostics present correlated signatures of azimuthal implosion asymmetry in recent indirect-drive inertial confinement fusion (ICF) implosion campaigns performed at the National Ignition Facility (NIF). The mean hot-spot velocity, inferred from the Doppler shift of 14 MeV neutrons produced by deuterium-tritium (DT) fusion, is systematically directed toward one azimuthal half of the NIF target chamber, centered on ϕ≈70°. Areal density (ρR) asymmetry of the converged DT fuel, inferred from nuclear activation diagnostics, presents a minimum ρR in the same direction as the hot-spot velocity and with ΔρR amplitude correlated with velocity magnitude. These two correlated observations, which are seen in all recent campaigns with cryogenic layers of DT fuel, are a known signature of asymmetry in the fuel convergence, implying a systematic azimuthal drive asymmetry across a wide range of shot and target configurations. The direction of the implied radiation asymmetry is observed to cluster toward the hohlraum diagnostic windows. This low-mode asymmetry degrades hot-spot conditions at peak convergence and limits implosion performance and yield.

Neutron Time-of-Flight Measurements of Charged-Particle Energy Loss in Inertial Confinement Fusion Plasmas

Neutron spectra from secondary ^{3}H(d,n)α reactions produced by an implosion of a deuterium-gas capsule at the National Ignition Facility have been measured with order-of-magnitude improvements in statistics and resolution over past experiments. These new data and their sensitivity to the energy loss of fast tritons emitted from thermal ^{2}H(d,p)^{3}H reactions enable the first statistically significant investigation of charged-particle stopping via the emitted neutron spectrum. Radiation-hydrodynamic simulations, constrained to match a number of observables from the implosion, were used to predict the neutron spectra while employing two different energy loss models. This analysis represents the first test of stopping models under inertial confinement fusion conditions, covering plasma temperatures of k_{B}T≈1-4  keV and particle densities of n≈(12-2)×10^{24}  cm^{-3}. Under these conditions, we find significant deviations of our data from a theory employing classical collisions whereas the theory including quantum diffraction agrees with our data.

Generating keV ion distributions for nuclear reactions at near solid-density using intense short-pulse lasers

Our understanding of a large range of astrophysical phenomena depends on a precise knowledge of charged particle nuclear reactions that occur at very low rates, which are difficult to measure under relevant plasma conditions. Here, we describe a method for generating dense plasmas at effective ion temperatures >20 keV, sufficient to induce measurable charged particle nuclear reactions. Our approach uses ultra-intense lasers to drive micron-sized, encapsulated nanofoam targets. Energetic electrons generated in the intense laser interaction pass through the foam, inducing a rapid expansion of the foam ions; this results in a hot, near-solid density plasma. We present the laser and target conditions necessary to achieve these conditions and illustrate the system performance using three-dimensional particle-in-cell simulations, outline potential applications and calculate expected nuclear reaction rates in the D(d,n) and 12C(p,γ) systems assuming CD, or CH aerogel foams.

Beam Energy and Centrality Dependence of Direct-Photon Emission from Ultrarelativistic Heavy-Ion Collisions

The PHENIX collaboration presents first measurements of low-momentum (0.4<p_{T}<3  GeV/c) direct-photon yields from Au+Au collisions at sqrt[s_{NN}]=39 and 62.4 GeV. For both beam energies the direct-photon yields are substantially enhanced with respect to expectations from prompt processes, similar to the yields observed in Au+Au collisions at sqrt[s_{NN}]=200. Analyzing the photon yield as a function of the experimental observable dN_{ch}/dη reveals that the low-momentum (>1  GeV/c) direct-photon yield dN_{γ}^{dir}/dη is a smooth function of dN_{ch}/dη and can be well described as proportional to (dN_{ch}/dη)^{α} with α≈1.25. This scaling behavior holds for a wide range of beam energies at the Relativistic Heavy Ion Collider and the Large Hadron Collider, for centrality selected samples, as well as for different A+A collision systems. At a given beam energy, the scaling also holds for high p_{T} (>5  GeV/c), but when results from different collision energies are compared, an additional sqrt[s_{NN}]-dependent multiplicative factor is needed to describe the integrated-direct-photon yield.

Ab initio response functions for Cherenkov-based neutron detectors

Neutron time-of-flight diagnostics at the NIF were recently outfitted with Cherenkov detectors. A fused silica radiator delivers sub-nanosecond response time and is optically coupled to a microchannel plate photomultiplier tube with gain from ∼1 to 10⁴. Capitalizing on fast time response and gamma-ray sensitivity, these systems can provide better than 30 ps precision for measuring first moments of neutron distributions. Generation of ab initio instrument response functions (IRFs) is critical to meet the <1% uncertainty needed. A combination of Monte Carlo modeling, benchtop characterization, and in situ comparison is employed. Close agreement is shown between the modeled IRFs and in situ measurements using the NIF's short-pulse advanced radiographic capability beams. First and second moments of neutron spectra calculated using ab initio IRFs agree well with established scintillator measurements. Next-step designs offer increased sensitivity and time-response.

Uncertainty analysis of response functions and γ -backgrounds on Tion and t0 measurements from Cherenkov neutron detectors at the National Ignition Facility (NIF)

Cherenkov radiators deployed to measure the neutron time-of-flight spectrum have response times associated with the neutron transit across the detector and are free from long time response tails characteristic of scintillation detectors. The Cherenkov radiation results from simple physical processes which makes them amenable to high fidelity Monte Carlo simulation. The instrument response function of neutron time-of-flight systems is a major contributor to both the systematic and statistical uncertainties of the parameters used to describe these spectra; in particular, the first and second moments of these distributions are associated with arrival time, t₀, and ion temperature, T_ion. We present the results of uncertainty analysis showing the significant reduction of the uncertainty in determining these quantities in the Cherenkov detector system recently deployed at NIF. The increased sensitivity to gamma radiation requires additional consideration of the effect of this background to the uncertainties in both t₀ and T_ion.

Using multiple neutron time of flight detectors to determine the hot spot velocity

An important diagnostic value of a shot at the National Ignition Facility is the resultant center-of-mass motion of the imploding capsule. This residual velocity reduces the efficiency of converting laser energy into plasma temperature. A new analysis method extracts the effective hot spot motion by using information from multiple neutron time-of-flight (nToF) lines-of-sight (LoSs). This technique fits a near Gaussian spectrum to the nToF scope traces and overcomes reliance on models to relate the plasma temperature to the mean energy of the emitted neutrons. This method requires having at least four nToF LoSs. The results of this analysis will be compared to an approach where each LoS is analyzed separately and a model is used to infer the mean energy of the emitted neutrons.

A fused silica Cherenkov radiator for high precision time-of-flight measurement of DT γ and neutron spectra (invited)

A fused silica Cherenkov radiator has been implemented at the National Ignition Facility to provide a new high precision measurement of the time-of-flight spectrum of 14.1 MeV DT fusion neutrons. This detector enables a high precision (<30 ps) co-registered measurement of both a thresholded γ-ray and a neutron spectrum on a single record. Other methods typically require γ and neutron signals to be co-registered via other diagnostics and/or dedicated timing experiments. Analysis of the co-registered γ and neutron signals allows precise extraction of the mean neutron energy and bulk hot-spot velocity, both of which were not possible with prior scintillator technologies. Initial measurements demonstrate the feasibility of this measurement and indicate that combined detection of neutrons and γ-rays on multiple lines-of-sight should enable the bulk vector velocity of the implosion hot-spot to be determined to ≈5 km/s and reduced uncertainty in the spectral width ≈0.1 keV.

Testing a Cherenkov neutron time-of-flight detector on OMEGA

A Cherenkov neutron time-of-flight (nTOF) detector developed and constructed at Lawrence Livermore National Laboratory was tested at 13 m from the target in a collimated line of sight (LOS) and at 5.3 m from the target in the open space inside the OMEGA Target Bay. Neutrons interacting with the quartz rod generate gammas, which through Compton scattering produce relativistic electrons that give rise to Cherenkov light. A photomultiplier tube (PMT) transferred the Cherenkov light into an amplified electrical signal. The Cherenkov nTOF detector consists of an 8-mm-diam, 25-cm quartz hexagonal prism coupled with a Hamamatsu gated PMT R5916U-52. The tests were performed with DT direct-drive implosions with cryogenic and room-temperature targets, producing a wide range of neutron yields and ion temperatures. The results of the tests and comparison with other nTOF detectors on OMEGA are presented. In the collimated LOS at 13 m from the target, the Cherenkov nTOF detector demonstrated good precision measurement in both the yield and ion temperature for DT yields above 3 × 10¹³.

Characterization of photodetector temporal response for neutron time-of-flight (nToF) diagnostics at the National Ignition Facility

The temporal response of a microchannel plate photomultiplier tube used in the suite of neutron time of flight (nToF) diagnostics at the National Ignition Facility has been characterized to reduce uncertainty in, and understanding of, shot parameters obtained from nTOF data. A short pulse laser, neutral density glass filters, and electrical attenuators were used to gather statistically significant samples of photodetector impulse response functions (IRF) in rapid succession. Individual components have been absolutely calibrated to minimize systematic uncertainties. The zeroth (collected charge), first (transit time), and second central moments (transit time spread) of the IRF were calculated as either the bias voltage or the amount of light incident on the detector was varied. Timing reference was provided by a monitor photodiode viewing a pickoff of the incident laser pulse. The primary sources of uncertainty are jitter in the monitor photodiode and the statistical variation across our measurement period. The spreads in the first moment, with respect to the timing photodiode, and the square root of the second central moment were found to be less than 50 ps and 150 ps, respectively.

Thermal Temperature Measurements of Inertial Fusion Implosions

Accurate measurement of the thermal temperature in inertially confined fusion plasmas is essential for characterizing ignition performance and validating the basic physics understanding of the stagnation conditions. We present experimental results from cryogenic deuterium-tritium implosions on the National Ignition Facility using a differential filter spectrometer designed to measure the thermal electron temperature from x-ray continuum emission from the stagnated plasma. Furthermore, electron temperature measurements, used in conjunction with the Doppler-broadened DT neutron spectra, allow one to infer the partition of energy in the hot spot between internal energy and unconverted kinetic energy.

Calibration of scintillation-light filters for neutron time-of-flight spectrometers at the National Ignition Facility

Sixty-four neutral density filters constructed of metal plates with 88 apertures of varying diameter have been radiographed with a soft x-ray source and CCD camera at National Security Technologies, Livermore. An analysis of the radiographs fits the radial dependence of the apertures' image intensities to sigmoid functions, which can describe the rapidly decreasing intensity towards the apertures' edges. The fitted image intensities determine the relative attenuation value of each filter. Absolute attenuation values of several imaged filters, measured in situ during calibration experiments, normalize the relative quantities which are now used in analyses of neutron spectrometer data at the National Ignition Facility.

Uncertainty analysis of signal deconvolution using a measured instrument response function

A common analysis procedure minimizes the ln-likelihood that a set of experimental observables matches a parameterized model of the observation. The model includes a description of the underlying physical process as well as the instrument response function (IRF). In the case investigated here, the National Ignition Facility (NIF) neutron time-of-flight (nTOF) spectrometers, the IRF is constructed from measurements and models. IRF measurements have a finite precision that can make significant contributions to determine the uncertainty estimate of the physical model's parameters. We apply a Bayesian analysis to properly account for IRF uncertainties in calculating the ln-likelihood function used to find the optimum physical parameters.

High-resolution measurements of the DT neutron spectrum using new CD foils in the Magnetic Recoil neutron Spectrometer (MRS) on the National Ignition Facility

The Magnetic Recoil neutron Spectrometer (MRS) on the National Ignition Facility measures the DT neutron spectrum from cryogenically layered inertial confinement fusion implosions. Yield, areal density, apparent ion temperature, and directional fluid flow are inferred from the MRS data. This paper describes recent advances in MRS measurements of the primary peak using new, thinner, reduced-area deuterated plastic (CD) conversion foils. The new foils allow operation of MRS at yields 2 orders of magnitude higher than previously possible, at a resolution down to ∼200 keV FWHM.

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.

Measurements of Elliptic and Triangular Flow in High-Multiplicity 3He+Au Collisions at √(s(NN))=200 GeV

We present the first measurement of elliptic (v(2)) and triangular (v(3)) flow in high-multiplicity (3)He+Au collisions at √(s(NN))=200  GeV. Two-particle correlations, where the particles have a large separation in pseudorapidity, are compared in (3)He+Au and in p+p collisions and indicate that collective effects dominate the second and third Fourier components for the correlations observed in the (3)He+Au system. The collective behavior is quantified in terms of elliptic v(2) and triangular v(3) anisotropy coefficients measured with respect to their corresponding event planes. The v(2) values are comparable to those previously measured in d+Au collisions at the same nucleon-nucleon center-of-mass energy. Comparisons with various theoretical predictions are made, including to models where the hot spots created by the impact of the three (3)He nucleons on the Au nucleus expand hydrodynamically to generate the triangular flow. The agreement of these models with data may indicate the formation of low-viscosity quark-gluon plasma even in these small collision systems.

First high-convergence cryogenic implosion in a near-vacuum hohlraum

Recent experiments on the National Ignition Facility [M. J. Edwards et al., Phys. Plasmas 20, 070501 (2013)] demonstrate that utilizing a near-vacuum hohlraum (low pressure gas-filled) is a viable option for high convergence cryogenic deuterium-tritium (DT) layered capsule implosions. This is made possible by using a dense ablator (high-density carbon), which shortens the drive duration needed to achieve high convergence: a measured 40% higher hohlraum efficiency than typical gas-filled hohlraums, which requires less laser energy going into the hohlraum, and an observed better symmetry control than anticipated by standard hydrodynamics simulations. The first series of near-vacuum hohlraum experiments culminated in a 6.8 ns, 1.2 MJ laser pulse driving a 2-shock, high adiabat (α∼3.5) cryogenic DT layered high density carbon capsule. This resulted in one of the best performances so far on the NIF relative to laser energy, with a measured primary neutron yield of 1.8×10(15) neutrons, with 20% calculated alpha heating at convergence ∼27×.

Azimuthal-angle dependence of charged-pion-interferometry measurements with respect to second- and third-order event planes in Au+Au collisions at √[S(NN)]=200 GeV

Charged-pion-interferometry measurements were made with respect to the second- and third-order event plane for Au+Au collisions at sqrt[s_{NN}]=200  GeV. A strong azimuthal-angle dependence of the extracted Gaussian-source radii was observed with respect to both the second- and third-order event planes. The results for the second-order dependence indicate that the initial eccentricity is reduced during the medium evolution, which is consistent with previous results. In contrast, the results for the third-order dependence indicate that the initial triangular shape is significantly reduced and potentially reversed by the end of the medium evolution, and that the third-order oscillations are largely dominated by the dynamical effects from triangular flow.

Performance of high-convergence, layered DT implosions with extended-duration pulses at the National Ignition Facility

Radiation-driven, low-adiabat, cryogenic DT layered plastic capsule implosions were carried out on the National Ignition Facility (NIF) to study the sensitivity of performance to peak power and drive duration. An implosion with extended drive and at reduced peak power of 350 TW achieved the highest compression with fuel areal density of ~1.3±0.1 g/cm2, representing a significant step from previously measured ~1.0 g/cm2 toward a goal of 1.5 g/cm2. Future experiments will focus on understanding and mitigating hydrodynamic instabilities and mix, and improving symmetry required to reach the threshold for thermonuclear ignition on NIF.

Medium modification of jet fragmentation in Au+Au collisions at √[s(NN)]=200 GeV measured in direct photon-hadron correlations

The jet fragmentation function is measured with direct photon-hadron correlations in p+p and Au+Au collisions at √[s(NN)]=200 GeV. The p(T) of the photon is an excellent approximation to the initial p(T) of the jet and the ratio z(T)=p(T)(h)/p(T)(γ) is used as a proxy for the jet fragmentation function. A statistical subtraction is used to extract the direct photon-hadron yields in Au+Au collisions while a photon isolation cut is applied in p+p. I(AA), the ratio of hadron yield opposite the photon in Au+Au to that in p+p, indicates modification of the jet fragmentation function. Suppression, most likely due to energy loss in the medium, is seen at high z(T). The associated hadron yield at low z(T) is enhanced at large angles. Such a trend is expected from redistribution of the lost energy into increased production of low-momentum particles.

Measurement of direct photons in Au+Au collisions at √(s(NN))=200 GeV

We report the measurement of direct photons at midrapidity in Au+Au collisions at √(s(NN))=200 GeV. The direct photon signal was extracted for the transverse momentum range of 4 GeV/c<p(T)<22 GeV/c, using a statistical method to subtract decay photons from the inclusive photon sample. The direct photon nuclear modification factor R(AA) was calculated as a function of p(T) for different Au+Au collision centralities using the measured p+p direct photon spectrum and compared to theoretical predictions. R(AA) was found to be consistent with unity for all centralities over the entire measured p(T) range. Theoretical models that account for modifications of initial direct photon production due to modified parton distribution functions in Au and the different isospin composition of the nuclei predict a modest change of R(AA) from unity. They are consistent with the data. Models with compensating effects of the quark-gluon plasma on high-energy photons, such as suppression of jet-fragmentation photons and induced-photon bremsstrahlung from partons traversing the medium, are also consistent with this measurement.

Neutron spectrometry--an essential tool for diagnosing implosions at the National Ignition Facility (invited)

DT neutron yield (Y(n)), ion temperature (T(i)), and down-scatter ratio (dsr) determined from measured neutron spectra are essential metrics for diagnosing the performance of inertial confinement fusion (ICF) implosions at the National Ignition Facility (NIF). A suite of neutron-time-of-flight (nTOF) spectrometers and a magnetic recoil spectrometer (MRS) have been implemented in different locations around the NIF target chamber, providing good implosion coverage and the complementarity required for reliable measurements of Y(n), T(i), and dsr. From the measured dsr value, an areal density (ρR) is determined through the relationship ρR(tot) (g∕cm(2)) = (20.4 ± 0.6) × dsr(10-12 MeV). The proportionality constant is determined considering implosion geometry, neutron attenuation, and energy range used for the dsr measurement. To ensure high accuracy in the measurements, a series of commissioning experiments using exploding pushers have been used for in situ calibration of the as-built spectrometers, which are now performing to the required accuracy. Recent data obtained with the MRS and nTOFs indicate that the implosion performance of cryogenically layered DT implosions, characterized by the experimental ignition threshold factor (ITFx), which is a function of dsr (or fuel ρR) and Y(n), has improved almost two orders of magnitude since the first shot in September, 2010.

Observation of direct-photon collective flow in Au + Au collisions at √s(NN)] = 200 GeV

The second Fourier component v(2) of the azimuthal anisotropy with respect to the reaction plane is measured for direct photons at midrapidity and transverse momentum (p(T)) of 1-12 GeV/c in Au + Au collisions at √s(NN)] = 200 GeV. Previous measurements of this quantity for hadrons with p(T) < 6 GeV/c indicate that the medium behaves like a nearly perfect fluid, while for p(T) > 6 GeV/c a reduced anisotropy is interpreted in terms of a path-length dependence for parton energy loss. In this measurement with the PHENIX detector at the Relativistic Heavy Ion Collider we find that for p(T) > 4 GeV/c the anisotropy for direct photons is consistent with zero, which is as expected if the dominant source of direct photons is initial hard scattering. However, in the p(T) < 4 GeV/c region dominated by thermal photons, we find a substantial direct-photon v(2) comparable to that of hadrons, whereas model calculations for thermal photons in this kinematic region underpredict the observed v(2).

Assembly of high-areal-density deuterium-tritium fuel from indirectly driven cryogenic implosions

The National Ignition Facility has been used to compress deuterium-tritium to an average areal density of ~1.0±0.1 g cm(-2), which is 67% of the ignition requirement. These conditions were obtained using 192 laser beams with total energy of 1-1.6 MJ and peak power up to 420 TW to create a hohlraum drive with a shaped power profile, peaking at a soft x-ray radiation temperature of 275-300 eV. This pulse delivered a series of shocks that compressed a capsule containing cryogenic deuterium-tritium to a radius of 25-35 μm. Neutron images of the implosion were used to estimate a fuel density of 500-800 g cm(-3).

Measurements of higher order flow harmonics in Au+Au collisions at √s(NN)=200 GeV

Flow coefficients ν(n) for n=2, 3, 4, characterizing the anisotropic collective flow in Au+Au collisions at √s(NN)=200 GeV, are measured relative to event planes Ψ(n), determined at large rapidity. We report ν(n) as a function of transverse momentum and collision centrality, and study the correlations among the event planes of different order n. The ν(n) are well described by hydrodynamic models which employ a Glauber Monte Carlo initial state geometry with fluctuations, providing additional constraining power on the interplay between initial conditions and the effects of viscosity as the system evolves. This new constraint can serve to improve the precision of the extracted shear viscosity to entropy density ratio η/s.

Cold nuclear matter effects on J/ψ yields as a function of rapidity and nuclear geometry in d+A collisions at sqrt[s(NN)]=200 GeV

We present measurements of J/ψ yields in d+Au collisions at sqrt[s(NN)]=200  GeV recorded by the PHENIX experiment and compare them with yields in p+p collisions at the same energy per nucleon-nucleon collision. The measurements cover a large kinematic range in J/ψ rapidity (-2.2<y<2.4) with high statistical precision and are compared with two theoretical models: one with nuclear shadowing combined with final state breakup and one with coherent gluon saturation effects. In order to remove model dependent systematic uncertainties we also compare the data to a simple geometric model. The forward rapidity data are inconsistent with nuclear modifications that are linear or exponential in the density weighted longitudinal thickness, such as those from the final state breakup of the bound state.

Azimuthal anisotropy of π⁰ production in Au+Au collisions at sqrt((s)NN)=200 GeV: path-length dependence of jet quenching and the role of initial geometry

We have measured the azimuthal anisotropy of π⁰ production for 1<p(T)<18  GeV/c for Au+Au collisions at sqrt((s)NN)=200  GeV. The observed anisotropy shows a gradual decrease for 3≲p(T)≲7-10  GeV/c, but remains positive beyond 10  GeV/c. The magnitude of this anisotropy is underpredicted, up to at least ∼10  GeV/c, by current perturbative QCD (PQCD) energy-loss model calculations. An estimate of the increase in anisotropy expected from initial-geometry modification due to gluon saturation effects and fluctuations is insufficient to account for this discrepancy. Calculations that implement a path-length dependence steeper than what is implied by current PQCD energy-loss models show reasonable agreement with the data.

Elliptic and hexadecapole flow of charged hadrons in Au+Au collisions at sq.rt(s(NN))=200 GeV

Differential measurements of the elliptic (v(2)) and hexadecapole (v(4)) Fourier flow coefficients are reported for charged hadrons as a function of transverse momentum (p(T)) and collision centrality or number of participant nucleons (N(part)) for Au+Au collisions at sq.rt(s(NN))=200  GeV. The v(2,4) measurements at pseudorapidity |η|≤0.35, obtained with four separate reaction-plane detectors positioned in the range 1.0<|η|<3.9, show good agreement, indicating the absence of significant Δη-dependent nonflow correlations. Sizable values for v(4)(p(T)) are observed with a ratio v(4)(p(T),N(part))/v(2)(2)(p(T),N(part))≈0.8 for 50≲N(part)≲200, which is compatible with the combined effects of a finite viscosity and initial eccentricity fluctuations. For N(part)≳200 this ratio increases up to 1.7 in the most central collisions.

Transition in yield and azimuthal shape modification in dihadron correlations in relativistic heavy ion collisions

Hard-scattered parton probes produced in collisions of large nuclei indicate large partonic energy loss, possibly with collective produced-medium response to the lost energy. We present measurements of π^{0} trigger particles at transverse momenta p{T}{t}=4-12  GeV/c and associated charged hadrons (p{T}{a}=0.5-7  GeV/c) vs relative azimuthal angle Δϕ in Au+Au and p+p collisions at sqrt[s{NN}]=200  GeV. The Au+Au distribution at low p{T}{a}, whose shape has been interpreted as a medium effect, is modified for p{T}{t}<7  GeV/c. At higher p{T}{t}, the data are consistent with unmodified or very weakly modified shapes, even for the lowest measured p{T}{a}, which quantitatively challenges some medium response models. The associated yield of hadrons opposing the trigger particle in Au+Au relative to p+p (I{AA}) is suppressed at high p{T} (I{AA}≈0.35-0.5), but less than for inclusive suppression (R{AA}≈0.2).

Enhanced production of direct photons in Au + Au collisions at square root(S(NN)) = 200 GeV and implications for the initial temperature

The production of e+ e- pairs for m(e+ e-)<0.3 GeV/c2 and 1<p(T)<5 GeV/c is measured in p+p and Au+Au collisions at square root(S(NN))=200 GeV. An enhanced yield above hadronic sources is observed. Treating the excess as photon internal conversions, the invariant yield of direct photons is deduced. In central Au+Au collisions, the excess of the direct photon yield over p+p is exponential in transverse momentum, with an inverse slope T=221+/-19(stat)+/-19(syst) MeV. Hydrodynamical models with initial temperatures ranging from T(init) approximately 300-600 MeV at times of approximately 0.6-0.15 fm/c after the collision are in qualitative agreement with the data. Lattice QCD predicts a phase transition to quark gluon plasma at approximately 170 MeV.

Charged Kaon interferometric probes of space-time evolution in Au+Au collisions at sqrt[S(NN)]=200 GeV

Bose-Einstein correlations of charged kaons are used to probe Au+Au collisions at sqrt[S(NN)]=200 GeV and are compared to charged pion probes, which have a larger hadronic scattering cross section. Three-dimensional Gaussian source radii are extracted, along with a one-dimensional kaon emission source function. The centrality dependences of the three Gaussian radii are well described by a single linear function of N(part)1/3 with a zero intercept. Imaging analysis shows a deviation from a Gaussian tail at r greater than or approximately equal to 10 fm, although the bulk emission at lower radius is well described by a Gaussian. The presence of a non-Gaussian tail in the kaon source reaffirms that the particle emission region in a heavy-ion collision is extended, and that similar measurements with pions are not solely due to the decay of long-lived resonances.

Measurement of bottom versus charm as a function of transverse momentum with electron-hadron correlations in p + p collisions at square root of s = 200 GeV

The momentum distribution of electrons from semileptonic decays of charm and bottom quarks for midrapidity |y|<0.35 in p+p collisions at square root of s=200 GeV is measured by the PHENIX experiment at the Relativistic Heavy Ion Collider over the transverse momentum range 2<pT<7 GeV/c. The ratio of the yield of electrons from bottom to that from charm is presented. The ratio is determined using partial D/D-->e(+/-)K(-/+)X (K unidentified) reconstruction. It is found that the yield of electrons from bottom becomes significant above 4 GeV/c in pT. A fixed-order-plus-next-to-leading-log perturbative quantum chromodynamics calculation agrees with the data within the theoretical and experimental uncertainties. The extracted total bottom production cross section at this energy is sigma(bb)=3.2(-1.1)(+1.2)(stat)(-1.3)(+1.4)(syst)mub.