Scattered neutron spectroscopy is a diagnostic technique commonly used to measure areal density in inertial confinement fusion experiments. Deleterious areal density asymmetries modify the shape of the scattered neutron spectrum. In this work, a novel analysis is developed, which can be used to fit the shape change. This will allow experimental scattered neutron spectroscopy to directly infer the amplitude and mode of the areal density asymmetries, with little sensitivity to confounding factors that affect other diagnostics for areal density. The model is tested on spectra produced by a neutron transport calculation with both isotropic and anisotropic primary fusion neutron sources. Multiple lines of sight are required to infer the areal density distribution over the whole sphere—we investigate the error propagation and optimal detector arrangement associated with the inference of mode 1 asymmetries.

Low mode drive asymmetries are proposed as a major degradation mechanism in inertial confinement fusion (ICF) implosions. When mode 1 drive asymmetries are present, several experimental signatures are measured; these include hotspot bulk flow velocity,1–4 anisotropy and differences in DT and DD inferred temperatures,5 and asymmetric hotspot shape.6 Hydrodynamic simulations suggest that the drive asymmetries, which can be used to explain these observations, will also cause corresponding fuel areal density asymmetries that reduce confinement.1,7–10 Recent work at the National Ignition Facility (NIF) has added to this evidence by exploring the relationship between hotspot velocity and areal density asymmetries measured using activation diagnostics.11 Activation diagnostics, while sensitive to areal density asymmetries, are also affected by hotspot velocity. Additional spectroscopic measurements are required to approximately remove this degeneracy.12 Scattered neutron images require sophisticated tomographic reconstruction techniques to extract the fuel density distribution.13 Scattered neutron spectroscopy is another method commonly used to measure fuel areal densities.14,15 It can be used at lower areal densities at which activation measurements are difficult, as is the case at OMEGA. While different areal densities can be measured on different detector lines of sight, the effect of areal density asymmetries on the scattered neutron spectrum has not been thoroughly investigated.

In this paper, we explore the effects of fuel areal density asymmetries on the scattered neutron spectrum. The kinematic scattering relationships suggest that properties of the areal density asymmetries should be contained within a single measured spectrum. A physics-based fitting model for the scattered neutron spectrum in the presence of areal density asymmetries is constructed from the nuclear cross sections. This will allow experimental spectroscopic measurements to leverage their high spectral resolution to directly infer more about the areal density distribution than currently possible. By combining multiple lines of sight, the mode, amplitude, and direction of the asymmetries can be measured.

In Sec. II, we will outline the geometry involved in neutron scattering. This will allow the areal density measured in neutron spectroscopy to be precisely defined. A model for a single scatter neutron spectrum from a DT ICF implosion will be described in Sec. III. This will be used to fit results from a neutron transport calculation in Sec. IV. Finally in Sec. V, how to combine results from multiple lines of sight will be discussed. This will include an error analysis as well as a methodology to evaluate the optimal detector arrangement for general projection measurements, such as the hotspot velocity.2 

The link between areal density asymmetries and scattered neutron spectral shape arises due to kinematic energy-angle relationships. These relationships are one-to-one for elastic processes, such as nD and nT scattering. From classical kinematics, in an elastic collision of a neutron with a stationary ion of mass Amn, the incoming neutron energy, E, and outcoming neutron energy, E, can be directly related to the scattering cosine, μs,

(1)

Conversely, inelastic processes, such as n(D,2n)p and n(T,2n)D, produce a distribution of energies for a given scattering angle and vice versa. A feature that is shared amongst all nuclear interactions of interest is azimuthal symmetry. This allows us to define a “scattering cone” at a given scattering angle (or cosine) around which the differential cross section is only energy dependent.

To understand the geometric interaction of scattering cones and areal density asymmetries, we will consider the simplified case for which the birth spectra of primary fusion neutrons are isotropic. This allows separation of the spatial and spectral dimensions. The spatial effects on the neutron spectrum come through the areal density around the scattering cone. In ICF implosions, a central hotspot acts as an extended source of neutrons. Therefore, we must consider the neutron-averaged line-integrated density, or ρL, rather than the more common ρR. Working in the scattering cone geometry, we find ρLs.c. along chords at a fixed scattering angle (with respect to our detector). For simplicity, we will consider the detector to lie along the z-axis, i.e., θ = 0. We then consider the areal density seen by a beam of neutrons born at r, with initial direction, Ω̂, before they scatter into the detector line of sight,

(2)

To neutron-average, we must sum over the whole reacting volume. The ρL and ρLs.c. are then defined as follows:

(3)
(4)
(5)

where Yn is the fusion neutron yield, ρ is the mass density, and Rn is the neutron production rate. The geometry of these chord integrals is shown in Fig. 1.

FIG. 1.

Diagram showing the geometry of the chord integrals given in Eqs. (3)–(5). Line integrals of the density, ρ, along direction Ω̂ are taken from each point inside the central source region V. The chord integrals are weighted by the neutron production rate, Rn, at the source point.

FIG. 1.

Diagram showing the geometry of the chord integrals given in Eqs. (3)–(5). Line integrals of the density, ρ, along direction Ω̂ are taken from each point inside the central source region V. The chord integrals are weighted by the neutron production rate, Rn, at the source point.

Close modal

For a point source, ρL will coincide with the ρR. However, this is not a good approximation for a typical ICF implosion. With perfect spherical symmetry, ρL is proportional to ρR with a coefficient, which depends on hydrodynamic profiles.14,16 In asymmetric conditions, ρL is well-defined and measurable, while ρR is not.

As ρL is defined over the surface of a sphere, it is natural to expand it in terms of spherical harmonics, Ylm,

(6)

Since the scattering cone integrates over the azimuthal angle, ϕs, it is simple to show that

(7)

where Pl are the Legendre polynomials and μs is the scattering cosine [cos(θs)]. Thus, a Legendre expansion in ρLs.c. is general as the spherical harmonic expansion of ρL is complete. The areal density Legendre coefficients, ρLl, are directly proportional to the m = 0 spherical harmonic coefficients. Modes with m0 vanish during azimuthal integration so do not affect ρLs.c.. Therefore, scattered spectra are blind to m0 modes. Note that these spherical harmonic expansions are defined in reference to a coordinate system centered on the detector line of sight. The m modes within the experimental geometry, e.g., those in the equatorial plane at NIF, are measurable with an appropriately orientated detector. The coefficients of order l along one line of sight are a linear combination of the coefficients, also of order l, on a different line of sight. Additionally, due to the extended nature of the hotspot neutron source, ρL is likely to be dominated by low modes.17 

For mode l = 1 areal density asymmetries, this is particularly simple,

(8)

where Ω̂det is the detector line of sight and ρL1 is a vector that points along the axis of the mode 1 and has a magnitude given by the half the peak to trough areal density asymmetry. The addition of an isotropic and vector component resembles the formula for the DT primary centroid shift.18 Therefore, a similar 4 detector setup2–4 is required in order to constrain the mode 1 areal density asymmetry using scattered spectra. For every higher mode, additional 2l+1 detectors are required in order to back out the physical ρLlm coefficients.

The energy spectrum of the scattered neutrons depends on the spectrum of primary fusion neutrons and the areal density distribution. The nuclear interaction differential cross sections provide the link between these spectral and spatial dependencies. Given an birth neutron energy spectrum, Qb(E), the scattered neutron spectrum is found as follows:

(9)
(10)

where m¯ is the average ion mass and the double differential cross section is the first term inside the integrand in Eq. (9). The dN/dE can be measured via time-of-flight or magnetic recoil spectrometers. The primary neutron spectrum is also measured and this can be used to separately constrain the birth spectrum; this leaves the areal density distribution as the single unknown. These equations will stand as the basis for a model to describe the singly scattered neutron spectrum in asymmetric ICF implosions.

Using the scattering cone integration, one can construct a full singly scattered neutron spectrum model for asymmetric ICF implosions. The single scatter approximation limits the scope of this model to areal densities 200 mg/cm2. In this work, we will focus on this areal density regime. The single scatter approximation also implies that the level of attenuation is low since the total cross section is dominated by scattering rather than absorption. Multiple scattering becomes increasingly important at higher areal densities. One can extend the scattering cone method to multiple scattering, but many of the approximations made are invalid. More accurate models for higher areal densities will be the subject of future work. Working with higher energy scattered neutrons can also be used to reduce the impact of multiple scattering.17 At these higher energies, both the scattering cross section and the number of multiply scattered neutrons are lower. At 200 mg/cm2, the ratio of doubly to singly scattered neutrons is just ∼2% within 10–12 MeV compared to ∼10% below the nT edge, and these rise linearly with increased areal density. Hence, the model limit of 200 mg/cm2 ensures suitable accuracy across the whole range of neutron energies considered.

The set of reactions we will consider are D(T,n)α, D(D,n)3He, T(T,2n)α, nT elastic, nD elastic, n(D,2n)p, and n(T,2n)D. Cross-sectional data for these reactions were taken from the ENDF19 and CENDL20 nuclear data libraries. The Bosch–Hale21 DT and DD fusion reactivities and the Appelbe22 calculation of the TT temperature dependent spectral shape were used. The DT and DD primary moments were calculated using relativistically correct expressions derived by Ballabio.23 For stationary target ions, the differential cross sections can simply be taken from nuclear data libraries, as will be done in this model. The effects of scattering ion velocities can be included for elastic processes using the method outlined in Crilly et al.24 These ion velocity effects are particularly important for the spectral shape of the backscatter edges. It is worth noting that most of the scattering interactions have been well measured,25,26 but this is not the case for n(T,2n)D. Models disagree on the cross section at 14 MeV (Refs. 19 and 27) and experimental data are limited28—this prompts further theoretical and experimental effort to reduce these uncertainties.

In current indirect drive implosions, remaining ablator areal densities29 can be several hundreds of mg/cm2. The higher mass nuclei mean that these ablator areal densities are generally less effective scatterers than equivalent fuel areal densities. However, scattering from the remaining ablator can prove a non-negligible contribution to the scattered neutron spectrum.17 A minimal extension to the spectral model, described above, could include the cross sections for ablator scattering interactions in a similar manner to those of the fuel. Carbon γ-ray measurements29–31 could be used to infer the carbon areal density as an input or constraint to this extended model. In this work, we will focus on the dominant fuel contribution to the scattered neutron spectrum.

First, we will look at the form of the scattering cone neutron spectrum as given in Eq. (9). The double differential cross section term can be calculated independently of the areal density distribution. The birth neutron spectrum of the primary fusion reactions can be calculated using the models listed above. For a stationary uniform plasma, the spectral shapes and yields of the primaries are only dependent on a single parameter, the ion temperature. More complex source conditions will alter the birth spectrum, but the model will aim to describe the birth spectra using only a single temperature, taken as the burn averaged DT ion temperature, TiDT. Figure 2 shows the differential cross section terms for a 4 keV ion temperature case.

FIG. 2.

(a) The DT double differential cross section integrated over all primary fusion sources. The two high intensity lines correspond to nD and nT elastic scattering of DT fusion neutrons. The white areas of the plot are kinematically forbidden. The spectral shapes and yield ratios of the primary neutron sources were calculated assuming TiDT = 4 keV and 50/50 DT isotopic composition. (b) A plot of the differential cross section with respect to outgoing energy for each scattering component. (c) A plot of the differential cross section with respect to scattering cosine for each scattering component. These single differential cross sections are marginalizations of the double differential cross section.

FIG. 2.

(a) The DT double differential cross section integrated over all primary fusion sources. The two high intensity lines correspond to nD and nT elastic scattering of DT fusion neutrons. The white areas of the plot are kinematically forbidden. The spectral shapes and yield ratios of the primary neutron sources were calculated assuming TiDT = 4 keV and 50/50 DT isotopic composition. (b) A plot of the differential cross section with respect to outgoing energy for each scattering component. (c) A plot of the differential cross section with respect to scattering cosine for each scattering component. These single differential cross sections are marginalizations of the double differential cross section.

Close modal

While limiting the energy range of a spectral measurement might be dominated by a particular angular range, there are still contributions from other scattering angles. This becomes increasingly important at lower energies where many scattering sources have similar amplitude. Areal density measurements found using a down-scatter ratio (DSR) technique will, therefore, find the following average measurement of the areal density distribution:

(11)
(12)

where the notation on the double differential cross sections indicates that they have been integrated over the DT peak only. Note we convert the DSR to a ρL, not a ρR. The conversion to ρR is sensitive to hydrodynamic profiles and, therefore, varies between implosions. Estimates for this geometry effect exist in the literature.14,16Figure 3 shows the angular weighting function for various DSR ranges. Inferring areal density asymmetries using DSR techniques from multiple lines of sight will, therefore, have to account for this averaging effect as well as the line of sight projection effect of the ρL distribution. Therefore, it may be advantageous to perform a fit to scattered spectrum over a large energy range to infer the ρLl coefficients directly. Fitting the spectrum also has the additional advantage that a single line of sight measurement can be used to quantify areal density asymmetries, whereas DSRs require multiple lines of sight to do this. This will be investigated further in the following sections of this work. Finally, we will consider the effect of cross section uncertainty on areal density inference. This can be estimated by calculating the fractional uncertainty in the average scattering cross section, σ. The uncertainty on the inelastic processes is likely the dominant source of error, leading to

(13)

where Δ indicates a fractional error and subscript i indicates the inelastic component of the cross section. The ratio σi/σ quantifies the fraction of scattered neutrons within the DSR range, which have undergone inelastic processes. This fraction is 2%, 11%, and 38% for 10–12, 6–12, and 3.5–4 MeV, respectively. We can, therefore, see that 3.5–4 MeV is the most sensitive to inelastic processes. For this range, fractional errors in inelastic cross sections of ∼10% (Ref. 26) produce ∼4% error in areal density inference.

FIG. 3.

The angular weighting function for various down scatter ratio ranges in MeV. The solid lines show the total, the dashed lines show the D contribution (from all scattering processes), and the dotted lines show the T contribution. In the 3.5–4 MeV range, the (n,2n), nD, and nT operate in distinct angular ranges, moving from small (μs1) to large (μs1) scattering angles, respectively.

FIG. 3.

The angular weighting function for various down scatter ratio ranges in MeV. The solid lines show the total, the dashed lines show the D contribution (from all scattering processes), and the dotted lines show the T contribution. In the 3.5–4 MeV range, the (n,2n), nD, and nT operate in distinct angular ranges, moving from small (μs1) to large (μs1) scattering angles, respectively.

Close modal

Given the scattering cone neutron spectrum, one can compute the resultant total neutron spectrum by integration over scattering cosine [as shown in Eq. (10)]. Figure 4 shows the spectrum for various l = 1 and l = 2 areal density asymmetry amplitudes. Deviations from the symmetric case are notable and different behavior is observed for l = 1 and l = 2. A simple interpretation based on the elastic energy-angle relationship in Eq. (1) fails to predict some of the behavior. For example, for l = 1, one might expect a positive mode 1 to be lower amplitude compared to the symmetric case in the backscatter regions. This is not the case due to the (n,2n) reactions, which have a very forwarded peaked differential cross section (c.f., Fig. 2). Their contribution to the elastic backscatter regions samples areal density from the forward direction—opposite to that sampled by the elastic interactions (c.f., the 3.5–4 MeV result in Fig. 3). A diagrammatic representation of this effect is given in Fig. 5. This demonstrates the need for the model outlined in this section in order to understand the effects of areal density asymmetries on neutron spectra.

FIG. 4.

The full primary and scattered neutron spectrum for various l = 1 (top) and l = 2 (bottom) areal density asymmetry amplitudes. The primary spectrum shapes correspond to a 4 keV 50/50 DT plasma. The ρL0 was taken to be 150 mg/cm2 for both cases.

FIG. 4.

The full primary and scattered neutron spectrum for various l = 1 (top) and l = 2 (bottom) areal density asymmetry amplitudes. The primary spectrum shapes correspond to a 4 keV 50/50 DT plasma. The ρL0 was taken to be 150 mg/cm2 for both cases.

Close modal
FIG. 5.

Diagram showing some example scattering trajectories for different nuclear interactions, which produce outgoing neutrons at 3.5 MeV. The hotspot and shell configuration, shown by the red and blue regions, respectively, results in a mode 1 areal density asymmetry in the direction indicated in the diagram.

FIG. 5.

Diagram showing some example scattering trajectories for different nuclear interactions, which produce outgoing neutrons at 3.5 MeV. The hotspot and shell configuration, shown by the red and blue regions, respectively, results in a mode 1 areal density asymmetry in the direction indicated in the diagram.

Close modal

The model laid out in Sec. III will be compared to results from a more complete neutron transport description of an ICF implosion. The inverse ray trace method described in Crilly et al.17 allows the calculation of neutron spectra given 3D grids of hydrodynamic conditions. As a suitable test problem, we consider an ice block model for a mode 1 asymmetry where a stationary spherical hotspot is displaced within a dense cold fuel layer, as in Fig. 5. The hotspot was chosen to have a parabolic temperature profile to test the spectral model's reliance on a single average temperature. The ice block configuration was then chosen such that TiDT = 3.9 keV, ρL0 = 150 mg/cm2, and ρL1/ρL0 = 30%. The detector line of sight was then rotated with respect to the mode 1 axis to explore the relationship derived in Eq. (8). Neutron spectra were calculated for these different detectors using the inverse ray trace method. These spectra were fitted using the model outlined in Sec. III to obtain ρL0 and ρL1 for each detector.

The results of this analysis are shown in Fig. 6. Every fit retrieved a value of ρL0 within 4 mg/cm2 of the true value of 150 mg/cm2. This error originates in the approximations of the model, namely, single temperature primary spectra and no attenuation. For ρL1, it was found that a fit to the spectrum between 3 and 6 MeV produced the best results. This range of energies is sensitive to all scattering angles (c.f., Fig. 2) but also narrow enough that the effect of differential attenuation of the scattered neutrons is small. Narrow high energy ranges are more robust to attenuation effects but have lower angular sensitivity.

FIG. 6.

(Top) Plot showing the neutron spectra from a mode 1 ice block model at various angles to the asymmetry axis. (Bottom) Results from fitting the spectra in two different energy ranges. These fits were used to infer the magnitude of mode 1 areal density asymmetry. The theoretical result from Eq. (8) is shown by the solid black line as a function of the cosine between detector and mode 1 direction, μ.

FIG. 6.

(Top) Plot showing the neutron spectra from a mode 1 ice block model at various angles to the asymmetry axis. (Bottom) Results from fitting the spectra in two different energy ranges. These fits were used to infer the magnitude of mode 1 areal density asymmetry. The theoretical result from Eq. (8) is shown by the solid black line as a function of the cosine between detector and mode 1 direction, μ.

Close modal

It is common for the primary spectra to exhibit anisotropy due to bulk fluid flow of the fusing plasma. For a mode 1 drive asymmetry, the hotspot velocity is aligned with the areal density asymmetry.1,7,8,11 According to the experimental data analysis by Rinderknecht et al.,11 the following linear relationship between mode 1 areal density asymmetry and hotspot velocity matches experimental data:

(14)

Thus, the 30% asymmetry presented here corresponds to a hotspot velocity ∼77 km/s. To see the effect of a spectral anisotropy of this magnitude, a uniform hotspot velocity (of 77 km/s) is introduced to the ice block model. For the fitting method, the spectral anisotropy was included through the birth spectrum term, Qb(E,μs,vHS), in Eq. (9) with the assumption that the mode 1 areal density asymmetry and hotspot flow velocity were aligned. It was found that the same results as shown for the stationary hotspot case were obtained when the hotspot flow velocity was properly accounted for. If the anisotropy of the birth spectrum was not included, then an additional absolute 1% error was introduced to the inferred ρL1/ρL0.

Finally, an areal density inference based on DSRs was performed. For a pure mode 1, the ρLDSR as defined in Eq. (12) can be simplified as

(15)
(16)

where μsDSR is the average scattering cosine being sampled in the DSR range. The 10–12 MeV (Ref. 14) and 3.5–4 MeV (Ref. 32) down scatter ranges have μsDSR values of 0.70 and −0.22, respectively (based on the results of Fig. 3). Between 3.5 and 4 MeV, if the nD and (n,2n) backgrounds under the nT backscatter edge can be successfully subtracted, this range becomes much more backward focused with a μsDSR = −0.89. However, the subtraction method must be suitably general as the scattering backgrounds are themselves sensitive to areal density asymmetries.

Unlike the fitting method, values of ρL0 and ρL1 cannot be inferred separately on a single line of sight when using DSRs. Figure 7 shows the inferred ρLDSR for the ice block spectra presented in Fig. 6. Without taking into account the average scattering cosine, μsDSR, the DSRs as a function of viewing angle are a poor representation of the underlying areal density asymmetry. This is due to the fact that the areal density probed by the DSR measurement is not along the detector line of sight. Because of this, interpreting the ρL through measured DSRs requires additional care. If the underlying areal density distribution is assumed to be mode 1, a linear fit can be performed once the measurements have been projected using the appropriate angle, μsDSR. For higher modes, Eq. (15) must be expanded to higher order introducing higher moments of μs, complicating further the relationship between DSR and underlying physical areal density asymmetries. The results of a mode 1 fit to the data in Fig. 7 are values of ρL0 = 148 mg/cm2 and ρL1/ρL0 = 30%, within a few percent of the true values. Here, we have used the fact we know the direction of the mode 1 a priori to analyze the data, and in the next section, we will consider the general case of a completely unknown ρL1.

FIG. 7.

Down scatter ratio inferred areal density [using Eq. (15)] as a function of the angle between the detector line of sight and the mode 1 axis. The results are projected using the average scattering cosine as given in Eq. (16), shown as hollow circles. The black dashed shows a linear fit to the projected areal densities.

FIG. 7.

Down scatter ratio inferred areal density [using Eq. (15)] as a function of the angle between the detector line of sight and the mode 1 axis. The results are projected using the average scattering cosine as given in Eq. (16), shown as hollow circles. The black dashed shows a linear fit to the projected areal densities.

Close modal

As shown in Sec. II, a single measurement of the scattered neutron spectrum contains incomplete information about the physical ρL(θ,ϕ) distribution. Multiple measurements are required in order to infer the underlying spherical harmonic moments. The inference of hotspot flows2,4 and apparent ion temperature anisotropy33 from primary neutron spectra requires similar analysis.

In this section, we will consider a more general approach to inferring the mode 1 areal density asymmetry vector, ρL1, and the error analysis. We will approach this as a general projection measurement first. We will consider the case where we have nd detectors and ndnp, where np is the number of physical parameters we are inferring. The measurements can then be related to the underlying physical parameters with the following system of linear equations:

(17)

where m is the length nd vector of measurements, A is the projection nd×np matrix, and p is the length np vector of physical parameters. The normal equation can be used to find the least squares distance between measurement and inferred physical parameters. This leads to the definition of the pseudo-inverse, A+,

(18)
(19)
(20)

We are also interested in how error propagates from measurement to inferred values. By considering the covariance matrix of p, which is denoted Σp, one can find a detector arrangement that minimizes this error,

(21)

In order to proceed with the analysis, we will assume that the detectors are independent and have the same error, σm. Then, the covariance matrix of the measurements Σm=σm2I. Using the properties of symmetric matrices for ATA, this simplifies Eq. (21),

(22)

The determinant of Σp is a measure of the total error in p. Therefore, the determinant of ATA must be maximized to minimize this error. This defines a method to evaluate error propagation and optimal detector setup for measurements which can be written as linear projection problems.2,33

We will investigate the optimal detector arrangement for multiple DSR measurements of a mode 1 areal density asymmetry as a worked example of this methodology. The projection matrix, A, for the DSR method can be found by combining Eqs. (8) and (15),

(23)
(24)

where the bracketed superscripts denote the detector index. If one sets μsDSR(i)=1, then this projection matrix also applies for the measurement of vector hotspot velocity and isotropic Gamow shifts for primary DT and DD neutron spectra.2,4 Evaluating ATA gives the following 4 × 4 matrix:

(25)
(26)

In order for the inferred quantities to have no covariances, then this matrix must be diagonal. Setting the off diagonal components to zero will set requirements on the detector lines of sight. As discussed earlier, the minimal error detector arrangement is the one that maximizes the determinant,

(27)
(28)

where Σu is the 3 × 3 covariance matrix of detector projection vectors, u. We will investigate the case where every detector uses the same DSR range such that μsDSR(i)=μsDSR. The angular factors can then be factored out,

(29)

leaving only dependence on the (unit length) line of sight vectors, Ω̂(i). Therefore, we can now use statistical results for the analysis of spherical data. The normalized orientation matrix, T=Ω̂Ω̂T, has positive eigenvalues, τi, that sum to unity.34 Its determinant is then maximized when all these eigenvalues are equal. This describes a uniform isotropic distribution for which

(30)

leading to

(31)

These requirements also lead to a diagonal ATA with elements (nd,nd/3,nd/3,nd/3). For nd = 4, the optimum is found for a tetrahedral arrangement. The solution to this problem is analogous to arranging masses on the unit sphere such that the center of gravity is at the origin and all principal moments of inertia are equal. This second constraint is satisfied if there are two axes with n-fold (n 3) rotational symmetry.35 

Using the analysis outlined in this section, we can calculate new optimal lines of sight for mode 1 projection problems given current experimental detector arrangements. This will apply for hotspot velocity inference and mode 1 areal density measurements, which use the same DSR range on all detectors. Given the current (nd = 4) NIF nToF arrangement,2 the next optimal line of sight would be along a (θϕ) of (77–68). At OMEGA, for the current nToF arrangement4 (nd = 5 unique lines of sight), the next optimal line of sight is along (135–301). The analysis is more involved if multiple different DSR ranges are used. This is because the scattering angle terms cannot be factored out [c.f., Eq. (29)]. The optimization problem must then be solved separately for any given combination of DSR ranges. Previous work in optimizing new lines of sight has been carried out for scattered neutron imaging using the NIF neutron imaging system (NIS).36 The optimal new detector location was taken as the one which would maximize coverage of the sky. Using the analysis in this section produces a different optimal line of sight to that which maximizes coverage. This highlights the different considerations required for different diagnostic techniques.

Scattered neutron spectroscopy contains valuable information on areal density asymmetries in ICF implosions. Areal density asymmetries lead to a reduction in confinement, and so accurate measurement of such asymmetries is key in understanding current perturbation sources and improving experimental performance.

In this work, we developed a model for the singly scattered neutron spectrum, which can fully account for the effects of an asymmetric areal density distribution. This required a more general description of the scattering geometry, which defined the measurable neutron averaged areal density, ρL. The model outlined in this work will allow fits to spectroscopic data to infer the amplitude and mode number of the areal density asymmetries. The model was tested on neutron transport results for a simple mode 1 test case. Fitting the spectra from multiple lines of sight reproduced the theoretical results, even in the presence of non-uniform temperature and anisotropic birth spectra. This analysis provides the most direct measurement of areal density with little sensitivity to other factors. However, going to higher areal densities presents challenges due to the increasing level of multiple scattering.

Due to the azimuthal symmetry of scattering, a single spectroscopic line of sight is insensitive to variation in ρL along the azimuthal direction. Therefore, combining measurements from a sufficient number of detectors is required to obtain a complete picture of the ρL distribution across the whole sphere. The error propagation for combining the measurements from multiple lines of sight was explored. This led to a methodology for calculating the optimal detector arrangement, which minimized the error in the inferred physical quantities. This error analysis is also relevant to the measurement of primary DT and DD spectra, which involve similar line of sight dependencies.2,4

This work was supported by the Lawrence Livermore National Laboratory through the Academic Partnership Program.

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

1.
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
(
6
),
062704
(
2020
).
2.
R.
Hatarik
,
R. C.
Nora
,
B. K.
Spears
,
M. J.
Eckart
,
G. P.
Grim
,
E. P.
Hartouni
,
A. S.
Moore
, and
D. J.
Schlossberg
, “
Using multiple neutron time of flight detectors to determine the hot spot velocity
,”
Rev. Sci. Instrum.
89
(
10
),
10I138
(
2018
).
3.
O. M.
Mannion
,
V. Y.
Glebov
,
C. J.
Forrest
,
J. P.
Knauer
,
V. N.
Goncharov
,
S. P.
Regan
,
T. C.
Sangster
,
C.
Stoeckl
, and
M.
Gatu Johnson
, “
Calibration of a neutron time-of-flight detector with a rapid instrument response function for measurements of bulk fluid motion on omega
,”
Rev. Sci. Instrum.
89
(
10
),
10I131
(
2018
).
4.
O. M.
Mannion
,
J. P.
Knauer
,
V. Y.
Glebov
,
C. J.
Forrest
,
A.
Liu
,
Z. L.
Mohamed
,
M. H.
Romanofsky
,
T. C.
Sangster
,
C.
Stoeckl
, and
S. P.
Regan
, “
A suite of neutron time-of-flight detectors to measure hot-spot motion in direct-drive inertial confinement fusion experiments on omega
,”
Nucl. Instrum. Methods Phys. Res., Sect. A
964
,
163774
(
2020
).
5.
M.
Gatu Johnson
,
J. P.
Knauer
,
C. J.
Cerjan
,
M. J.
Eckart
,
G. P.
Grim
,
E. P.
Hartouni
,
R.
Hatarik
,
J. D.
Kilkenny
,
D. H.
Munro
,
D. B.
Sayre
,
B. K.
Spears
,
R. M.
Bionta
,
E. J.
Bond
,
J. A.
Caggiano
,
D.
Callahan
,
D. T.
Casey
,
T.
Döppner
,
J. A.
Frenje
,
V. Y.
Glebov
,
O.
Hurricane
,
A.
Kritcher
,
S.
LePape
,
T.
Ma
,
A.
Mackinnon
,
N.
Meezan
,
P.
Patel
,
R. D.
Petrasso
,
J. E.
Ralph
,
P. T.
Springer
, and
C. B.
Yeamans
, “
Indications of flow near maximum compression in layered deuterium-tritium implosions at the National Ignition Facility
,”
Phys. Rev. E
94
(
2
),
021202
(
2016
).
6.
J. J.
Ruby
,
A.
Pak
,
J. E.
Field
,
T.
Ma
,
B. K.
Spears
,
L. R.
Benedetti
,
D. K.
Bradley
,
L. F.
Berzak Hopkins
,
D. T.
Casey
,
T.
Döppner
,
D.
Eder
,
D.
Fittinghoff
,
G.
Grim
,
R.
Hatarik
,
D. E.
Hinkel
,
N.
Izumi
,
J. D.
Kilkenny
,
S. F.
Khan
,
J. P.
Knauer
,
A. L.
Kritcher
,
F. E.
Merrill
,
J. D.
Moody
,
S. R.
Nagel
,
H.-S.
Park
,
J. D.
Salmonson
,
D. B.
Sayre
,
D. A.
Callahan
,
W. W.
Hsing
,
O. A.
Hurricane
,
P. K.
Patel
, and
M. J.
Edwards
, “
Spatially resolved x-ray emission measurements of the residual velocity during the stagnation phase of inertial confinement fusion implosion experiments
,”
Phys. Plasmas
23
(
7
),
072701
(
2016
).
7.
B. K.
Spears
,
M. J.
Edwards
,
S.
Hatchett
,
J.
Kilkenny
,
J.
Knauer
,
A.
Kritcher
,
J.
Lindl
,
D.
Munro
,
P.
Patel
,
H. F.
Robey
 et al., “
Mode 1 drive asymmetry in inertial confinement fusion implosions on the National Ignition Facility
,”
Phys. Plasmas
21
(
4
),
042702
(
2014
).
8.
J. P.
Chittenden
,
B. D.
Appelbe
,
F.
Manke
,
K.
McGlinchey
, and
N. P. L.
Niasse
, “
Signatures of asymmetry in neutron spectra and images predicted by three-dimensional radiation hydrodynamics simulations of indirect drive implosions
,”
Phys. Plasmas
23
(
5
),
052708
(
2016
).
9.
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
,
M. M.
Marinak
,
L. P.
Masse
,
D. H.
Munro
,
M. V.
Patel
,
P. K.
Patel
,
H. F.
Robey
,
C. R.
Schroeder
,
S. M.
Sepke
, and
M. J.
Edwards
, “
Three-dimensional modeling and hydrodynamic scaling of National Ignition Facility implosions
,”
Phys. Plasmas
26
(
5
),
050601
(
2019
).
10.
K.
McGlinchey
,
B. D.
Appelbe
,
A. J.
Crilly
,
J. K.
Tong
,
C. A.
Walsh
, and
J. P.
Chittenden
, “
Diagnostic signatures of performance degrading perturbations in inertial confinement fusion implosions
,”
Phys. Plasmas
25
(
12
),
122705
(
2018
).
11.
H. G.
Rinderknecht
,
D. T.
Casey
,
R.
Hatarik
,
R. M.
Bionta
,
B. J.
MacGowan
,
P.
Patel
,
O. L.
Landen
,
E. P.
Hartouni
, and
O. A.
Hurricane
, “
Azimuthal drive asymmetry in inertial confinement fusion implosions on the National Ignition Facility
,”
Phys. Rev. Lett.
124
,
145002
(
2020
).
12.
H. G.
Rinderknecht
,
R.
Bionta
,
G.
Grim
,
R.
Hatarik
,
H.
Khater
,
D.
Schlossberg
, and
C.
Yeamans
, “
Velocity correction for neutron activation diagnostics at the NIF
,”
Rev. Sci. Instrum.
89
(
10
),
10I125
(
2018
).
13.
P. L.
Volegov
,
S. H.
Batha
,
V.
Geppert-Kleinrath
,
C. R.
Danly
,
F. E.
Merrill
,
C. H.
Wilde
,
D. C.
Wilson
,
D. T.
Casey
,
D.
Fittinghoff
,
B.
Appelbe
,
J. P.
Chittenden
,
A. J.
Crilly
, and
K.
McGlinchey
, “
Density determination of the thermonuclear fuel region in inertial confinement fusion implosions
,”
J. Appl. Phys.
127
(
8
),
083301
(
2020
).
14.
M.
Gatu Johnson
,
J. A.
Frenje
,
D. T.
Casey
,
C. K.
Li
,
F. H.
Séguin
,
R.
Petrasso
,
R.
Ashabranner
,
R. M.
Bionta
,
D. L.
Bleuel
,
E. J.
Bond
,
J. A.
Caggiano
,
A.
Carpenter
,
C. J.
Cerjan
,
T. J.
Clancy
,
T.
Doeppner
,
M. J.
Eckart
,
M. J.
Edwards
,
S.
Friedrich
,
S. H.
Glenzer
,
S. W.
Haan
,
E. P.
Hartouni
,
R.
Hatarik
,
S. P.
Hatchett
,
O. S.
Jones
,
G.
Kyrala
,
S. L.
Pape
,
R. A.
Lerche
,
O. L.
Landen
,
T.
Ma
,
A. J.
MacKinnon
,
M. A.
McKernan
,
M. J.
Moran
,
E.
Moses
,
D. H.
Munro
,
J.
McNaney
,
H. S.
Park
,
J.
Ralph
,
B.
Remington
,
J. R.
Rygg
,
S. M.
Sepke
,
V.
Smalyuk
,
B.
Spears
,
P. T.
Springer
,
C. B.
Yeamans
,
M.
Farrell
,
D.
Jasion
,
J. D.
Kilkenny
,
A.
Nikroo
,
R.
Paguio
,
J. P.
Knauer
,
V. Y.
Glebov
,
T. C.
Sangster
,
R.
Betti
,
C.
Stoeckl
,
J.
Magoon
,
M. J.
Shoup
,
G. P.
Grim
,
J.
Kline
,
G. L.
Morgan
,
T. J.
Murphy
,
R. J.
Leeper
,
C. L.
Ruiz
,
G. W.
Cooper
, and
A. J.
Nelson
, “
Neutron spectrometry—An essential tool for diagnosing implosions at the National Ignition Facility (invited)
,”
Rev. Sci. Instrum.
83
(
10
),
10D308
(
2012
).
15.
V. Y.
Glebov
,
C. J.
Forrest
,
K. L.
Marshall
,
M.
Romanofsky
,
T. C.
Sangster
,
M. J.
Shoup
, and
C.
Stoeckl
, “
A new neutron time-of-flight detector for fuel-areal-density measurements on omega
,”
Rev. Sci. Instrum.
85
(
11
),
11E102
(
2014
).
16.
J. A.
Frenje
,
R.
Bionta
,
E. J.
Bond
,
J. A.
Caggiano
,
D. T.
Casey
,
C.
Cerjan
,
J.
Edwards
,
M.
Eckart
,
D. N.
Fittinghoff
,
S.
Friedrich
,
V. Y.
Glebov
,
S.
Glenzer
,
G.
Grim
,
S.
Haan
,
R.
Hatarik
,
S.
Hatchett
,
M.
Gatu Johnson
,
O. S.
Jones
,
J. D.
Kilkenny
,
J. P.
Knauer
,
O.
Landen
,
R.
Leeper
,
S. L.
Pape
,
R.
Lerche
,
C. K.
Li
,
A.
Mackinnon
,
J.
McNaney
,
F. E.
Merrill
,
M.
Moran
,
D. H.
Munro
,
T. J.
Murphy
,
R. D.
Petrasso
,
R.
Rygg
,
T. C.
Sangster
,
F. H.
Séguin
,
S.
Sepke
,
B.
Spears
,
P.
Springer
,
C.
Stoeckl
, and
D. C.
Wilson
, “
Diagnosing implosion performance at the National Ignition Facility (NIF) by means of neutron spectrometry
,”
Nucl. Fusion
53
(
4
),
043014
(
2013
).
17.
A. J.
Crilly
,
B. D.
Appelbe
,
K.
McGlinchey
,
C. A.
Walsh
,
J. K.
Tong
,
A. B.
Boxall
, and
J. P.
Chittenden
, “
Synthetic nuclear diagnostics for inferring plasma properties of inertial confinement fusion implosions
,”
Phys. Plasmas
25
(
12
),
122703
(
2018
).
18.
D. H.
Munro
, “
Interpreting inertial fusion neutron spectra
,”
Nucl. Fusion
56
(
3
),
036001
(
2016
).
19.
M. B.
Chadwick
,
P.
Obložinský
,
M.
Herman
,
N. M.
Greene
,
R. D.
McKnight
,
D. L.
Smith
,
P. G.
Young
,
R. E.
MacFarlane
,
G. M.
Hale
,
S. C.
Frankle
,
A. C.
Kahler
,
T.
Kawano
,
R. C.
Little
,
D. G.
Madland
,
P.
Moller
,
R. D.
Mosteller
,
P. R.
Page
,
P.
Talou
,
H.
Trellue
,
M. C.
White
,
W. B.
Wilson
,
R.
Arcilla
,
C. L.
Dunford
,
S. F.
Mughabghab
,
B.
Pritychenko
,
D.
Rochman
,
A. A.
Sonzogni
,
C. R.
Lubitz
,
T. H.
Trumbull
,
J. P.
Weinman
,
D. A.
Brown
,
D. E.
Cullen
,
D. P.
Heinrichs
,
D. P.
McNabb
,
H.
Derrien
,
M. E.
Dunn
,
N. M.
Larson
,
L. C.
Leal
,
A. D.
Carlson
,
R. C.
Block
,
J. B.
Briggs
,
E. T.
Cheng
,
H. C.
Huria
,
M. L.
Zerkle
,
K. S.
Kozier
,
A.
Courcelle
,
V.
Pronyaev
, and
S. C.
van der Marck
, “
ENDF/B-VII.0: Next generation evaluated nuclear data library for nuclear science and technology
,”
Nucl. Data Sheets
107
(
12
),
2931
(
2006
).
20.
Z. G.
Ge
,
Z. X.
Zhao
,
H. H.
Xia
,
Y. X.
Zhuang
,
T. J.
Liu
,
J. S.
Zhang
, and
H. C.
Wu
, “
The updated version of Chinese evaluated nuclear data library (CENDL-3.1)
,”
J. Korean Phys. Soc.
59
,
1052
(
2011
).
21.
H.-S.
Bosch
and
G. M.
Hale
, “
Improved formulas for fusion cross-sections and thermal reactivities
,”
Nucl. Fusion
32
(
4
),
611
(
1992
).
22.
B.
Appelbe
and
J.
Chittenden
, “
The effects of ion temperature on the energy spectra of t + t → 2n + α reaction products
,”
High Energy Density Phys.
19
,
29
37
(
2016
).
23.
L.
Ballabio
,
J. K.
Allne
, and
G.
Gorini
, “
Relativistic calculation of fusion product spectra
,”
Nucl. Fusion
38
(
11
),
1723
(
1998
).
24.
A. J.
Crilly
,
B. D.
Appelbe
,
O. M.
Mannion
,
C. J.
Forrest
,
V.
Gopalaswamy
,
C. A.
Walsh
, and
J. P.
Chittenden
, “
Neutron backscatter edge: A measure of the hydrodynamic properties of the dense DT fuel at stagnation in ICF experiments
,”
Phys. Plasmas
27
(
1
),
012701
(
2020
).
25.
J. A.
Frenje
,
C. K.
Li
,
F. H.
Seguin
,
D. T.
Casey
,
R. D.
Petrasso
,
D. P.
McNabb
,
P.
Navratil
,
S.
Quaglioni
,
T. C.
Sangster
,
V.
Yu Glebov
, and
D. D.
Meyerhofer
, “
Measurements of the differential cross sections for the elastic n3H and n2H scattering at 14.1 MeV by using an inertial confinement fusion facility
,”
Phys. Rev. Lett.
107
,
122502
(
2011
).
26.
C. J.
Forrest
,
A.
Deltuva
,
W. U.
Schröder
,
A. V.
Voinov
,
J. P.
Knauer
,
E. M.
Campbell
,
G. W.
Collins
,
V. Y.
Glebov
,
O. M.
Mannion
,
Z. L.
Mohamed
,
P. B.
Radha
,
S. P.
Regan
,
T. C.
Sangster
, and
C.
Stoeckl
, “
Deuteron breakup induced by 14-MeV neutrons from inertial confinement fusion
,”
Phys. Rev. C
100
,
034001
(
2019
).
27.
A.
Deltuva
and
A. C.
Fonseca
, “
Neutron-3H scattering above the four-nucleon breakup threshold
,”
Phys. Rev. C
86
,
011001
(
2012
).
28.
D. S.
Mather
and
L. F.
Pain
, “
Measurement of (n, 2n) and (n, 3n) cross sections at 14-MeV incident energy
,”
Technical Report No. AWRE-O-47/69 Atomic Weapons Research Establishment
, Aldermaston (England),
1969
.
29.
K. D.
Meaney
,
Y.
Kim
,
H.
Geppert-Kleinrath
,
H. W.
Herrmann
,
L.
Berzak Hopkins
,
N. M.
Hoffman
,
C.
Cerjan
,
O. L.
Landen
,
K.
Baker
,
J.
Carrera
 et al., “
Carbon ablator areal density at fusion burn: Observations and trends at the National Ignition Facility
,”
Phys. Plasmas
27
(
5
),
052702
(
2020
).
30.
N. M.
Hoffman
,
H. W.
Herrmann
,
Y. H.
Kim
,
H. H.
Hsu
,
C. J.
Horsfield
,
M. S.
Rubery
,
E. K.
Miller
,
E.
Grafil
,
W.
Stoeffl
,
J. A.
Church
,
C. S.
Young
,
J. M.
Mack
,
D. C.
Wilson
,
J. R.
Langenbrunner
,
S. C.
Evans
,
T. J.
Sedillo
,
V. Y.
Glebov
, and
T.
Duffy
, “
Measurement of areal density in the ablators of inertial-confinement-fusion capsules via detection of ablator (n, n′γ) gamma-ray emission
,”
Phys. Plasmas
20
(
4
),
042705
(
2013
).
31.
K. D.
Meaney
,
Y. H.
Kim
,
H.
Geppert-Kleinrath
,
H. W.
Herrmann
,
L. B.
Hopkins
, and
N. M.
Hoffman
, “
Diagnostic signature of the compressibility of the inertial-confinement-fusion pusher
,”
Phys. Rev. E
101
,
023208
(
2020
).
32.
C. J.
Forrest
,
P. B.
Radha
,
V. Y.
Glebov
,
V. N.
Goncharov
,
J. P.
Knauer
,
A.
Pruyne
,
M.
Romanofsky
,
T. C.
Sangster
,
M. J.
Shoup
 III
,
C.
Stoeckl
,
D. T.
Casey
,
M.
Gatu-Johnson
, and
S.
Gardner
, “
High-resolution spectroscopy used to measure inertial confinement fusion neutron spectra on omega (invited)
,”
Rev. Sci. Instrum.
83
(
10
),
10D919
(
2012
).
33.
K. M.
Woo
,
R.
Betti
,
D.
Shvarts
,
O. M.
Mannion
,
D.
Patel
,
V. N.
Goncharov
,
K. S.
Anderson
,
P. B.
Radha
,
J. P.
Knauer
,
A.
Bose
,
V.
Gopalaswamy
,
A. R.
Christopherson
,
E. M.
Campbell
,
J.
Sanz
, and
H.
Aluie
, “
Impact of three-dimensional hot-spot flow asymmetry on ion-temperature measurements in inertial confinement fusion experiments
,”
Phys. Plasmas
25
(
10
),
102710
(
2018
).
34.
N. I.
Fisher
,
T.
Lewis
, and
B. J. J.
Embleton
,
Statistical Analysis of Spherical Data
(
Cambridge University Press
,
1993
).
35.
P. K.
Aravind
, “
A comment on the moment of inertia of symmetrical solids
,”
Am. J. Phys.
60
(
8
),
754
755
(
1992
).
36.
S. H.
Batha
,
P. L.
Volegov
,
V. E.
Fatherley
,
V.
Geppert-Kleinrath
, and
C. H.
Wilde
, “
Optimizing neutron imaging line of sight locations for maximizing sampling of the cold fuel density in inertial confinement fusion implosions at the National Ignition Facility
,”
Rev. Sci. Instrum.
89
(
10
),
10I147
(
2018
).