Neutron Wall Load in Fusion & UCN Systems

• NWL is a design metric that quantifies neutron power/flux on first-wall surfaces in fusion devices and analogous loss rates in UCN traps.
• It informs key engineering decisions such as component survivability, heat load management, and shielding placement using analytic geometric models.
• NWL analysis integrates neutron polarization, diffuse reflection laws, and transport modeling to balance reactor design and operational performance.

Searching arXiv for the specified papers on neutron wall loading and related neutron–wall interactions. Neutron wall load (NWL) denotes a wall-interaction quantity whose precise meaning depends on the neutron-physics subfield. In fusion-device design, NWL is the neutron power or flux incident on the first wall and is treated as a geometric, design-oriented metric for identifying where neutrons strike the machine boundary, with direct relevance to shielding, blanket design, magnet lifetime, component survivability, and peaking-factor control (Schwartz, 15 Jul 2025). In ultracold-neutron (UCN) material-trap physics, the closely analogous quantity is the wall-loss rate caused primarily by absorption at trap walls, with additional sensitivity to geometry, gravity, and diffuse reflection laws (Grigoriev et al., 2024). In confined UCN bottles, neutron–wall Casimir–Polder attraction contributes to the effective wall potential, even though that literature does not define NWL as a separate formula (Higa et al., 2018). In neutron-detector instrumentation, “neutron wall” refers instead to a large detector array, and the corresponding electronics literature addresses high-load operation operationally rather than through an explicit NWL definition (Qian-Shun et al., 2014).

1. Terminological scope

The term is not used uniformly across the cited literature. Its dominant formal usage is in fusion engineering, where it refers to neutron power or flux incident on a wall surface. In UCN lifetime work, the analogous wall quantity is a neutron-loss rate from material boundaries. In neutron-wall detector electronics, the expression “neutron wall” names the detector system itself, and the cited work does not define an NWL metric (Schwartz, 15 Jul 2025).

Context	Quantity of interest	Status of the term “NWL”
Spin-polarized fusion in axisymmetric tori	Neutron power/flux incident on first wall	Explicitly defined
UCN material traps	Wall-loss rate from absorption, with inelastic upscattering as an additional contribution	Analogue of NWL
UCN near conducting walls	Casimir–Polder wall contribution to effective trapping potential	Not a separate NWL formula
Neutron wall detector electronics	Operational detector loading, fast timing/charge readout under high-rate conditions	Not explicitly defined

This multiplicity of usage matters because superficially similar phrases refer to distinct observables: direct neutron power loading in fusion, neutron survival losses in UCN storage, wall-induced electromagnetic potentials in bottle confinement, and electronics burden in detector systems. A plausible implication is that “NWL” should always be interpreted together with the wall model, transport model, and experimental context.

2. Fusion NWL as a direct geometric wall-flux quantity

For a wall patch at position $\boldsymbol{r}_w$, the fusion paper defines neutron wall loading as neutron energy times a wall-flux integral over the source volume,

$$Q^{\mathrm{NWL}}(\boldsymbol{r}_w)=\frac{E_n}{4\pi}\int_{V_s} dV \left(\Phi^{\mathrm{source}}\;\hat{n}\cdot\frac{\boldsymbol{r}_s-\boldsymbol{r}_w}{|\boldsymbol{r}_s-\boldsymbol{r}_w|^3}\right),$$

where $E_n$ is the neutron energy per reaction, $V_s$ is the source volume, $\Phi^{\mathrm{source}}$ is the local neutron source density, and $\hat{n}(\boldsymbol{r}_w)$ is the outward wall normal (Schwartz, 15 Jul 2025). The kernel is a free-space geometric transport factor normalized by $1/(4\pi)$ for isotropic emission.

The same work is explicit that this NWL is a simplified direct source-to-wall quantity. It drops the energy spread by treating $E_n$ as constant for every neutron and does not include neutron scattering, blanket transport, secondary scattering in the walls, or Monte Carlo neutron propagation (Schwartz, 15 Jul 2025). Accordingly, the formalism is not a full neutronics calculation.

Within those assumptions, NWL is presented as a first-order design quantity. It identifies where neutrons hit the first wall, supports assessment of heat and damage loading, informs first-wall and blanket placement, and motivates use of the peaking factor defined as max load / average load (Schwartz, 15 Jul 2025). The cited design contexts include tokamaks, spherical tokamaks, and related axisymmetric devices, with particular concern for wall regions such as the center stack and divertor-like structures. This framing makes NWL primarily a geometry-sensitive engineering observable rather than a reaction-rate quantity alone.

3. Spin polarization and anisotropic neutron emission

The fusion analysis uses the standard polarized DT formalism. Deuterium has three spin projections and tritium has two, combined into three effective modes,

$a=d_+t_+ + d_-t_-,$

$b=d_0t_+ + d_0t_-,$

$$Q^{\mathrm{NWL}}(\boldsymbol{r}_w)=\frac{E_n}{4\pi}\int_{V_s} dV \left(\Phi^{\mathrm{source}}\;\hat{n}\cdot\frac{\boldsymbol{r}_s-\boldsymbol{r}_w}{|\boldsymbol{r}_s-\boldsymbol{r}_w|^3}\right),$$0

with $$Q^{\mathrm{NWL}}(\boldsymbol{r}_w)=\frac{E_n}{4\pi}\int_{V_s} dV \left(\Phi^{\mathrm{source}}\;\hat{n}\cdot\frac{\boldsymbol{r}_s-\boldsymbol{r}_w}{|\boldsymbol{r}_s-\boldsymbol{r}_w|^3}\right),$$1 and $$Q^{\mathrm{NWL}}(\boldsymbol{r}_w)=\frac{E_n}{4\pi}\int_{V_s} dV \left(\Phi^{\mathrm{source}}\;\hat{n}\cdot\frac{\boldsymbol{r}_s-\boldsymbol{r}_w}{|\boldsymbol{r}_s-\boldsymbol{r}_w|^3}\right),$$2 for non-polarized plasma (Schwartz, 15 Jul 2025).

The corresponding differential cross section is

$$Q^{\mathrm{NWL}}(\boldsymbol{r}_w)=\frac{E_n}{4\pi}\int_{V_s} dV \left(\Phi^{\mathrm{source}}\;\hat{n}\cdot\frac{\boldsymbol{r}_s-\boldsymbol{r}_w}{|\boldsymbol{r}_s-\boldsymbol{r}_w|^3}\right),$$3

where $$Q^{\mathrm{NWL}}(\boldsymbol{r}_w)=\frac{E_n}{4\pi}\int_{V_s} dV \left(\Phi^{\mathrm{source}}\;\hat{n}\cdot\frac{\boldsymbol{r}_s-\boldsymbol{r}_w}{|\boldsymbol{r}_s-\boldsymbol{r}_w|^3}\right),$$4 is the angle relative to the local magnetic field $$Q^{\mathrm{NWL}}(\boldsymbol{r}_w)=\frac{E_n}{4\pi}\int_{V_s} dV \left(\Phi^{\mathrm{source}}\;\hat{n}\cdot\frac{\boldsymbol{r}_s-\boldsymbol{r}_w}{|\boldsymbol{r}_s-\boldsymbol{r}_w|^3}\right),$$5 (Schwartz, 15 Jul 2025). NWL therefore depends not only on source placement and wall geometry but also on emission anisotropy induced by fuel polarization.

Three limiting polarization modes are emphasized. In A mode, $$Q^{\mathrm{NWL}}(\boldsymbol{r}_w)=\frac{E_n}{4\pi}\int_{V_s} dV \left(\Phi^{\mathrm{source}}\;\hat{n}\cdot\frac{\boldsymbol{r}_s-\boldsymbol{r}_w}{|\boldsymbol{r}_s-\boldsymbol{r}_w|^3}\right),$$6, the total cross section is $$Q^{\mathrm{NWL}}(\boldsymbol{r}_w)=\frac{E_n}{4\pi}\int_{V_s} dV \left(\Phi^{\mathrm{source}}\;\hat{n}\cdot\frac{\boldsymbol{r}_s-\boldsymbol{r}_w}{|\boldsymbol{r}_s-\boldsymbol{r}_w|^3}\right),$$7, there is a 50% enhancement relative to unpolarized DT, and neutron and alpha emission are proportional to $$Q^{\mathrm{NWL}}(\boldsymbol{r}_w)=\frac{E_n}{4\pi}\int_{V_s} dV \left(\Phi^{\mathrm{source}}\;\hat{n}\cdot\frac{\boldsymbol{r}_s-\boldsymbol{r}_w}{|\boldsymbol{r}_s-\boldsymbol{r}_w|^3}\right),$$8, hence preferentially perpendicular to $$Q^{\mathrm{NWL}}(\boldsymbol{r}_w)=\frac{E_n}{4\pi}\int_{V_s} dV \left(\Phi^{\mathrm{source}}\;\hat{n}\cdot\frac{\boldsymbol{r}_s-\boldsymbol{r}_w}{|\boldsymbol{r}_s-\boldsymbol{r}_w|^3}\right),$$9. In B mode, $E_n$0, the total cross section is $E_n$1, equal to unpolarized DT, while the emission pattern is proportional to $E_n$2, preferentially parallel to $E_n$3. In C mode, $E_n$4, the angular distribution is the same as in B mode but the total cross section is reduced to $E_n$5 (Schwartz, 15 Jul 2025).

The paper also identifies two isotropic-emission mixtures: $E_n$6, with a 20% enhanced total cross section and isotropic emission, and $E_n$7, with 75% of the unpolarized cross section and isotropic emission (Schwartz, 15 Jul 2025). These combinations show that changes in total reaction rate and changes in neutron directionality need not coincide. A plausible implication is that polarization choice is intrinsically system-dependent because fusion gain, neutron directionality, alpha deposition and heating, and wall-load constraints must be balanced simultaneously.

4. Analytic construction in axisymmetric geometries

The core analytic strategy is to represent the source as infinitesimal filamentary ring sources in an axisymmetric torus. This is motivated by the fact that axisymmetry reduces the source to a ring at fixed radius and height, while more complicated volumetric sources can be approximated as weighted sums of rings whose individual contributions can be computed analytically (Schwartz, 15 Jul 2025).

For a source ring at radius $E_n$8, height $E_n$9, and toroidal angle $V_s$0,

$V_s$1

and for a wall target $V_s$2 the source-to-target vector is $V_s$3 (Schwartz, 15 Jul 2025). The wall contribution involves the geometric dilution factor $V_s$4, the incidence factor $V_s$5, and, for polarized sources, the angular factor $V_s$6.

The paper first derives intuition-building results for a rectangular-cross-section torus, including the inboard wall, floor and ceiling, and outboard wall. For an inboard-wall target, source visibility is restricted to the arc $V_s$7 with $V_s$8, and the resulting isotropic and polarized reduced intensity kernels are given in closed form using incomplete elliptic integrals $V_s$9 and $\Phi^{\mathrm{source}}$0 (Schwartz, 15 Jul 2025). The B/C kernel is constructed from an isotropic part and a $\Phi^{\mathrm{source}}$1-weighted auxiliary kernel.

The analysis is then generalized to arbitrary axisymmetric magnetic-field direction using a polar angle $\Phi^{\mathrm{source}}$2 away from toroidal direction and an azimuthal angle $\Phi^{\mathrm{source}}$3 describing the tilt plane. Although the integrands become more complicated, they remain analytically reducible to sums of standard component integrals (Schwartz, 15 Jul 2025). The main geometric generalization is to a torus with convex poloidal cross section. Convexity ensures that each source ring is visible from every wall point only over a limited arc, bounded by line-of-sight tangency to an outward-facing wall segment. The limiting angle is obtained from a tangency condition,

$\Phi^{\mathrm{source}}$4

with associated contact fraction and contact point definitions involving $\Phi^{\mathrm{source}}$5, $\Phi^{\mathrm{source}}$6, and $\Phi^{\mathrm{source}}$7 (Schwartz, 15 Jul 2025). The paper also specifies failure cases in which a segment does not limit the view, namely negative discriminant, $\Phi^{\mathrm{source}}$8, or contact height outside the segment.

A further simplification is the decomposition of the wall response into horizontal and vertical normal components. For a wall element with normal angle $\Phi^{\mathrm{source}}$9,

$\hat{n}(\boldsymbol{r}_w)$0

the total kernel is written as

$\hat{n}(\boldsymbol{r}_w)$1

with analogous generalized formulas for isotropic, A-mode, and B/C-mode emission (Schwartz, 15 Jul 2025). This separates geometry from wall orientation and is central to rapid evaluation.

The closed forms share a common dependence on algebraic prefactors, trigonometric functions of $\hat{n}(\boldsymbol{r}_w)$2, $\hat{n}(\boldsymbol{r}_w)$3, and $\hat{n}(\boldsymbol{r}_w)$4, and incomplete elliptic integrals

$\hat{n}(\boldsymbol{r}_w)$5

Because the formulas are analytic, fully differentiable, and implemented in the JAX-compatible Python package anarr, they are intended for fast design exploration, forward-mode and reverse-mode differentiation, and just-in-time compilation in optimization workflows (Schwartz, 15 Jul 2025). They are explicitly presented as suitable for scoping studies and first-wall shaping, including construction of shells with nearly uniform NWL, but not as a replacement for full Monte Carlo neutron transport or detailed blanket and shield engineering.

5. UCN wall-loss rate as an NWL analogue

In material UCN traps, the corresponding wall quantity is the neutron loss rate caused mainly by absorption by trap walls, with inelastic upscattering as an additional contribution. This loss rate is central to neutron lifetime experiments because the measured storage lifetime is shorter than the true beta-decay lifetime by an amount set by wall losses (Grigoriev et al., 2024).

The microscopic input is the absorption probability per wall collision for a neutron with normal incident velocity component $\hat{n}(\boldsymbol{r}_w)$6,

$\hat{n}(\boldsymbol{r}_w)$7

where $\hat{n}(\boldsymbol{r}_w)$8 is the wall-loss coefficient and $\hat{n}(\boldsymbol{r}_w)$9 is the limiting velocity corresponding to the optical potential $1/(4\pi)$0 (Grigoriev et al., 2024). Under the usual isotropic-incidence angular average,

$1/(4\pi)$1

or equivalently

$1/(4\pi)$2

with $1/(4\pi)$3 (Grigoriev et al., 2024).

In the gravity-neglecting limit and for isotropic velocity distributions, the wall-loss rate reduces to the familiar surface-to-volume estimate

$1/(4\pi)$4

which is the standard NWL-like scaling with wall area $1/(4\pi)$5, neutron speed $1/(4\pi)$6, average absorption probability, and volume $1/(4\pi)$7 (Grigoriev et al., 2024). The paper argues, however, that this standard picture is incomplete because gravity makes the UCN angular distribution height-dependent. A central distinction is drawn between the number distribution $1/(4\pi)$8 at the trap bottom and the density distribution $1/(4\pi)$9 at height $E_n$0; these are not the same in gravity, even if the bottom distribution is isotropic (Grigoriev et al., 2024).

For rectangular and cylindrical traps, the paper derives analytic “exact” methods based on the isotropic bottom distribution, including exact collision counts and survival probabilities. In a rectangular trap, for example,

$E_n$1

which already distinguishes bottom-wall losses from side-wall losses (Grigoriev et al., 2024). The standard gravity-including method, the oversimplified no-gravity method, and the exact bottom-distribution treatment can differ substantially. For side walls and total losses the reported differences can exceed 10% in experimentally relevant parameter ranges, and the paper states that wall-loss modeling can shift the extracted neutron lifetime by seconds, potentially contributing to the approximately 4 s tension between magnetic-trap and material-trap measurements (Grigoriev et al., 2024).

A major controversy concerns diffuse elastic reflections. The diffuse scattering probability is stated to be small, $E_n$2, but much larger than the absorption probability $E_n$3, so it can strongly reshape the phase-space distribution over many collisions (Grigoriev et al., 2024). Standard Monte Carlo practice often assumes Lambert’s cosine law,

$E_n$4

whereas the paper also studies an isotropic diffuse law. Its conclusion is not that one method is universally wrong, but that the standard gravity-including formula agrees with Monte Carlo only if the diffuse-wall physics closely resembles Lambert’s law. This suggests that accurate UCN wall-load estimation requires more than an $E_n$5 parameterization: it also requires the correct phase-space distribution, collision kinematics in gravity, and diffuse-reflection law.

6. Neutron–wall Casimir–Polder interaction in confined UCN systems

A separate but related wall effect arises from the neutron’s electromagnetic structure. Even though the neutron is neutral, it has electric and magnetic dipole polarizabilities, so vacuum fluctuations induce Casimir–Polder (CP) interactions with conductors. The relevant neutron responses are the dynamical dipole polarizabilities $E_n$6 and $E_n$7, fitted to chiral EFT results up to the pion-production threshold and at the onset of the $E_n$8 resonance, and used through the imaginary-frequency continuation $E_n$9 and $a=d_+t_+ + d_-t_-,$0 (Higa et al., 2018).

For a neutron at distance $a=d_+t_+ + d_-t_-,$1 from a conducting wall, the CP potential is

$a=d_+t_+ + d_-t_-,$2

with

$a=d_+t_+ + d_-t_-,$3

This interaction is attractive overall, falls off as a long-range interaction, and is controlled by the neutron’s electric polarizability together with retardation encoded in the exponential and polynomial factor (Higa et al., 2018).

For one neutron between two conducting walls separated by distance $a=d_+t_+ + d_-t_-,$4, with neutron position $a=d_+t_+ + d_-t_-,$5 measured from the midpoint, the potential becomes

$a=d_+t_+ + d_-t_-,$6

This is the bottle-like geometry of direct relevance to confined UCN, with an overall $a=d_+t_+ + d_-t_-,$7 scale and explicit dependence on neutron position $a=d_+t_+ + d_-t_-,$8 (Higa et al., 2018).

The paper emphasizes the standard crossover from shorter-distance van der Waals-like behavior to asymptotic Casimir–Polder tails and reports that the dynamic-polarizability curve has smaller magnitude than the static one at comparable distances (Higa et al., 2018). For bottle confinement, this means that the actual wall contribution to the neutron’s effective potential is weaker than an overly strong static estimate would suggest. The same work explicitly notes that these attractive CP forces can compete with the usual repulsive Fermi pseudopotential of the wall material, for example about 252 neV for Ni, making them relevant to the practical confinement conditions of ultracold neutrons inside bottles (Higa et al., 2018). A plausible implication is that any wall-load analysis for UCN bottles that focuses only on the repulsive optical potential omits a geometry-dependent long-range attraction.

7. Detector-wall usage and operational loading

In detector instrumentation, the phrase “neutron wall” refers to a large array detector rather than a first-wall or trap-wall boundary. The cited electronics paper presents a single-width NIM module with eight channels of TAC and QAC designed for the neutron wall detector in the HIRFL-CSR External Target Facility, where the detector measures neutrons in the approximate energy range of several tens of MeV up to 1 GeV in a time-of-flight scheme (Qian-Shun et al., 2014).

The module converts the START–STOP time interval to a voltage through a TAC and measures detector pulse charge through a QAC, enabling simultaneous time and charge readout from the same event. The TAC input range is 30 ns to 1 $a=d_+t_+ + d_-t_-,$9s, and the QAC input range is 40 pC to 600 pC. Reported performance includes TAC linearity error lower than 1.28% and time resolution less than 0.871%, and QAC linearity error lower than 0.81% with resolution better than 0.936% (Qian-Shun et al., 2014). The module also features discharge time below 200 ns, low output noise and offset, fast current splitting, gated integration, simple circuit structure, and low power dissipation.

This paper does not define or calculate NWL explicitly (Qian-Shun et al., 2014). Its relevance is indirect: it addresses fast neutron-detector signals, simultaneous timing and charge extraction, high count-rate or high-load conditions, and front-end readout for large array detectors. That is an operational notion of loading on the neutron-wall readout chain, not a formal neutron wall loading metric in the sense used in fusion engineering. The distinction is terminologically important because identical words refer here to detector occupancy and electronics burden rather than to neutron power incident on a confining wall.