pith. sign in

arxiv: 2605.20432 · v2 · pith:WUDTRJXEnew · submitted 2026-05-19 · ⚛️ physics.plasm-ph · nucl-ex

NeSST: A Python Tool for Neutron Spectra and Synthetic Diagnostics in Inertial Confinement Fusion

Aidan Crilly This is my paper

Pith reviewed 2026-05-21 07:09 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph nucl-ex
keywords neutron spectrainertial confinement fusionsynthetic diagnosticsENDF cross sectionsneutron time-of-flightICF implosionsscattering kernelsLegendre expansions
0
0 comments X

The pith

NeSST is a Python tool that constructs primary and singly scattered neutron spectra from ICF implosions using ENDF cross sections, relativistic kinematics, Legendre expansions for asymmetries, and ion-velocity kernels.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents NeSST, an open-source Python package for rapidly building primary and scattered neutron spectra in inertial confinement fusion implosions. It evaluates spectra for DT, DD, and TT reactions plus scattering from ablator materials by directly reading ENDF libraries for elastic and inelastic cross sections. Features include relativistic corrections to scattering kinematics, Legendre mode expansions to capture areal density asymmetries from implosion non-uniformities, and pre-computed kernels that account for ion velocities shaping the backscatter edge. A synthetic neutron time-of-flight module then converts the energy spectra into detector signals that include configurable instrument responses, scintillator sensitivity, and beamline attenuation. This matters because neutron spectra serve as primary diagnostics for fusion performance, and a fast, configurable tool supports analysis and fitting without full transport simulations.

Core claim

NeSST reads differential and total nuclear cross sections from ENDF libraries to compute elastic scattering such as nD and nT along with inelastic channels like D(n,2n)p and T(n,2n)D for DT fuel and additional ablator materials such as 12C. The code applies relativistic corrections to elastic scattering kinematics, incorporates areal density asymmetries through a Legendre mode expansion of the neutron-averaged projected areal density, handles the effect of scattering ion velocities on the neutron backscatter edge via pre-computed velocity-dependent kernels, and converts the resulting energy spectra into synthetic nToF signals using a full forward model with instrument response functions,能量 -

What carries the argument

NeSST itself, which combines ENDF-sourced cross sections with a single-scattering framework, Legendre expansions for areal-density asymmetries, and pre-computed ion-velocity kernels to generate spectra and synthetic diagnostics.

If this is right

  • Spectral signatures of implosion non-uniformities can be computed and fitted directly from measured neutron spectra.
  • Synthetic nToF detector signals can be generated that include full instrument response, energy-dependent scintillator sensitivity, and beamline attenuation.
  • Scattering contributions from ablator materials such as carbon can be treated in the same framework as fuel scattering.
  • Relativistic kinematics improve the modeling of high-energy neutron scattering edges.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The open-source release allows rapid iteration on diagnostic interpretation for ongoing ICF experiments at facilities that record nToF data.
  • Integration with hydrodynamic simulation outputs could enable parameter studies of how implosion asymmetries translate into observable neutron spectra.
  • The modular structure supports extension to additional reaction channels or multi-scatter corrections if experimental data require them.

Load-bearing premise

Single scattering dominates the relevant spectral features and Legendre expansions of neutron-averaged projected areal density plus pre-computed velocity kernels sufficiently capture non-uniformity and backscatter edge effects without requiring full Monte Carlo transport.

What would settle it

A direct comparison of NeSST-generated backscatter edges and asymmetry signatures against full Monte Carlo neutron transport results for the same asymmetric ICF implosion conditions would show clear discrepancies if the single-scattering and Legendre approximations are insufficient.

Figures

Figures reproduced from arXiv: 2605.20432 by Aidan Crilly.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) The DT spectrum using various spectral shape [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Scattered spectrum components for an equimolar [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) Scattered spectra for symmetric ( [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Comparison of NeSST (dashed, N) with Minotaur [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. (a) Transmission of 15-m dry air at STP as a func [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: shows the effect on the time-resolved nToF signals of a burn duration of 100 ps with a linearly de￾creasing ion temperature with a central ion temperature of 5 keV. These finite burn duration effects can be in￾cluded in the forward model and used to fit time-resolved plasma conditions, including sensitivity and IRF effects, [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Synthetic code demonstration: forward-model self [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. (a) Craun, Smith & Bethe sensitivity models includ [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
read the original abstract

We present the Neutron Scattered Spectra Tool (NeSST), an open-source Python package for rapidly constructing primary and singly scattered neutron spectra from inertial confinement fusion (ICF) implosions. NeSST evaluates primary spectra for deuterium-tritium (DT), deuterium-deuterium (DD) and tritium-tritium (TT) reactions. Differential and total nuclear cross sections are read directly from Evaluated Nuclear Data File (ENDF) libraries. This enables elastic ($n$D, $n$T) and inelastic [D$(n,2n)$p, T$(n,2n)$D] scattering from DT fuel, as well as scattering from additional ablator materials such as $^{12}$C, to be treated within a common framework. Relativistic corrections to elastic scattering kinematics are included. Areal density asymmetries are incorporated through a Legendre mode expansion of the neutron-averaged projected areal density, allowing the spectral signatures of implosion non-uniformities to be computed and fitted. The effect of scattering ion velocities on the neutron backscatter edge shape is handled through pre-computed ion-velocity-dependent scattering kernels. A synthetic neutron time-of-flight (nToF) module converts energy spectra into detector signals with a full forward model that includes configurable instrument response functions (IRFs), energy-dependent scintillator sensitivity, and beamline attenuation. The code is publicly available at https://github.com/aidancrilly/NeSST

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript presents NeSST, an open-source Python package for rapidly constructing primary and singly scattered neutron spectra from ICF implosions. It supports DT, DD, and TT reactions plus ablator materials by reading differential and total cross sections directly from ENDF libraries, incorporates relativistic kinematics for elastic scattering, uses a Legendre mode expansion of the neutron-averaged projected areal density to model asymmetries, employs pre-computed ion-velocity-dependent scattering kernels for the backscatter edge, and includes a configurable synthetic nToF forward model with instrument response functions, scintillator sensitivity, and beamline attenuation.

Significance. If the implementation is accurate within its stated approximations, NeSST would provide a practical, efficient tool for generating synthetic neutron spectra and diagnostics in the ICF community. The open-source release on GitHub and reliance on standard ENDF libraries support reproducibility and community adoption. The handling of implosion non-uniformities via Legendre expansions and velocity effects on the backscatter edge addresses relevant experimental needs for nToF analysis at facilities such as NIF.

major comments (2)
  1. Scattering model (Methods): The central construction relies on the single-scattering approximation plus Legendre kernels and pre-computed velocity kernels, but no quantitative bound or test is given for the ρR range (e.g., ≳100 mg/cm²) where double scattering begins to fill the backscatter edge and alter down-scattered yield at the few-percent level required for diagnostic use.
  2. Validation and results sections: The manuscript contains no benchmark comparisons to full Monte Carlo transport codes, no error propagation analysis, and no direct comparison to experimental nToF data, leaving the claim that the generated spectra are 'usable' for diagnostic inference only partially supported.
minor comments (1)
  1. Code availability statement: While the GitHub link is given, the text would benefit from a short paragraph describing the module structure and providing a minimal usage example for constructing a spectrum with asymmetry modes.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive review and recommendation for minor revision. We address each major comment below and have revised the manuscript to strengthen the presentation of the scattering model and validation aspects while remaining within the scope of a tool-description paper.

read point-by-point responses

  1. Referee: Scattering model (Methods): The central construction relies on the single-scattering approximation plus Legendre kernels and pre-computed velocity kernels, but no quantitative bound or test is given for the ρR range (e.g., ≳100 mg/cm²) where double scattering begins to fill the backscatter edge and alter down-scattered yield at the few-percent level required for diagnostic use.

    Authors: We agree that an explicit discussion of the single-scattering approximation's range of validity is valuable. In the revised manuscript we have added a dedicated paragraph in the Methods section that cites prior Monte Carlo studies on multiple-scattering contributions in ICF neutron spectra. We state that for the areal-density range typical of current experiments (ρR ≲ 100 mg/cm²) the double-scattering correction to the backscatter edge remains below a few percent, consistent with the diagnostic precision requirements, and we note the approximation's breakdown at higher ρR. This addition provides the requested quantitative context without new calculations. revision: yes

  2. Referee: Validation and results sections: The manuscript contains no benchmark comparisons to full Monte Carlo transport codes, no error propagation analysis, and no direct comparison to experimental nToF data, leaving the claim that the generated spectra are 'usable' for diagnostic inference only partially supported.

    Authors: The manuscript's primary purpose is to document the NeSST implementation and its underlying physics models rather than to perform an exhaustive validation campaign. We have revised the Results section to include a brief discussion of accuracy expectations based on the ENDF libraries and relativistic kinematics, together with references to earlier literature that benchmarked comparable single-scattering models against Monte Carlo transport codes. Error propagation is facilitated by the code's modular structure, which allows users to propagate uncertainties from the input cross sections; we have clarified this in the text. Direct experimental nToF comparisons lie outside the present scope but are enabled by the synthetic diagnostic module; we have added a forward-looking statement indicating that such comparisons will be pursued in follow-on work. revision: partial

Circularity Check

0 steps flagged

NeSST tool implements standard external models and data without circular derivations

full rationale

The paper presents NeSST as a forward-modeling Python package that reads differential and total cross sections directly from ENDF libraries, applies relativistic elastic scattering kinematics, uses Legendre expansions of neutron-averaged projected areal density for asymmetries, and employs pre-computed ion-velocity-dependent scattering kernels. All core components are drawn from established nuclear data libraries and standard physics treatments rather than from any self-fitted parameters, internal definitions, or self-citation chains. The synthetic nToF module is likewise a configurable forward operator built on these external inputs. No load-bearing step reduces by construction to its own outputs, satisfying the criteria for a self-contained, non-circular implementation.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The tool depends on external nuclear data libraries and standard domain assumptions in ICF neutron modeling rather than introducing new free parameters or entities.

axioms (2)
  • domain assumption ENDF libraries contain accurate differential and total nuclear cross sections for elastic and inelastic scattering on D, T, and C
    Cross sections are read directly from ENDF for all scattering calculations.
  • domain assumption Singly scattered neutrons dominate the spectral signatures of interest in the relevant energy range
    Tool focuses exclusively on primary and singly scattered spectra.

pith-pipeline@v0.9.0 · 5786 in / 1487 out tokens · 35853 ms · 2026-05-21T07:09:14.853620+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.