pith. sign in

arxiv: 2506.06927 · v3 · submitted 2025-06-07 · ⚛️ physics.pop-ph · physics.hist-ph· physics.soc-ph

An elementary method to determine the critical mass of a sphere of fissile material based on a separation of neutron transport and nuclear reaction processes

Steven K. Lamoreaux This is my paper

Pith reviewed 2026-05-19 11:36 UTC · model grok-4.3

classification ⚛️ physics.pop-ph physics.hist-phphysics.soc-ph
keywords critical massfissile sphereneutron path lengthelementary calculationneutron multiplicationcritical radiusfission chain reaction
0
0 comments X

The pith

A statistical threshold on average neutron path length inside a sphere gives the critical radius for fissile material without solving the diffusion equation.

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

The paper shows that the critical mass of a pure fissile sphere follows directly from requiring that the average neutron path length produce enough new neutrons to offset absorption and surface leakage, with the path length linking the nuclear reaction rates to the sphere geometry. This separation lets the calculation use only elementary statistics and calculus instead of transport equations. The resulting radius matches tabulated critical masses to within a few percent for common fissile isotopes. The same approach extends immediately to materials with impurities or isotopic mixes and supplies quick order-of-magnitude checks for more elaborate simulations.

Core claim

By imposing the steady-state condition that the total number of neutrons inside the sphere remains constant, the average distance a neutron travels before escaping or being absorbed can be set equal to the path length required to generate replacement neutrons from fission; this length, together with the mean free path between scatterings, fixes the sphere radius directly.

What carries the argument

The threshold condition that neutron number does not change with time, which isolates the nuclear reaction rates from the geometric transport problem through the single shared quantity of total path length.

If this is right

  • An algebraic expression for critical radius is obtained directly from the path-length balance.
  • The method applies without change to impure or isotopically mixed fissile material.
  • It supplies a rapid design guide or consistency check for Monte Carlo neutronics runs.
  • Direct numerical comparison is possible with the Oppenheimer-Bethe formula and with diffusion-equation solutions under adjusted boundary conditions.

Where Pith is reading between the lines

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

  • The same path-length balance could be adapted to estimate critical sizes for non-spherical shapes by replacing the sphere escape probability with the appropriate geometric factor.
  • Because the derivation uses only average path length, it offers a transparent starting point for teaching how leakage and multiplication compete before students encounter the full transport equation.

Load-bearing premise

The nuclear reaction rates and the geometric transport can be connected solely through the total distance a neutron travels inside the material.

What would settle it

Apply the derived radius formula to a known fissile sphere such as 93-percent uranium-235 and check whether the predicted critical mass falls outside the few-percent agreement window with established experimental values.

Figures

Figures reproduced from arXiv: 2506.06927 by Steven K. Lamoreaux.

Figure 1
Figure 1. Figure 1: FIG. 1: Critical radius for [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Number of first-generation neutrons produced [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Geometrical picture for calculating the mean [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
read the original abstract

A simplified method to calculate the critical mass of a fissile material sphere is presented. This is a purely pedagogical study, in part to elucidate the historical evolution of criticality calculations. This method employs only elementary calculus and straightforward statistical arguments by formulating the problem in terms of the threshold condition that the number of neutrons in the sphere does not change with time; the average neutron path length in the material must be long enough to produce enough fission neutrons to balance losses by absorption due to nuclear reactions and leakage through the surface. This separates the nuclear reaction part of the problem from the geometry and mechanics of neutron transport, the only connection being the total path length which together with the distance between scatterings determines the sphere radius. This leads to an expression for the critical radius without the need to solve the diffusion equation. Comparison with known critical masses shows agreement at the few-percent level. The analysis can also be applied to impure materials, isotopically or otherwise, and can be extended to general neutronics estimations as a design guide or for order-of-magnitude checking of Monte Carlo N-Particle (MCNP) simulations. A comparison is made with the Oppenheimer-Bethe criticality formula, with the results of other calculations, and with the diffusion equation approach via a new treatment of the boundary conditions.

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 / 2 minor

Summary. The manuscript presents a simplified pedagogical method to calculate the critical mass of a fissile material sphere. It formulates the problem via the stationary-neutron-number condition, using elementary calculus and statistical path-length arguments to balance fission neutron production against absorption and leakage losses. Nuclear reactions are separated from transport geometry, with the sole connection being the total average neutron path length that determines the critical radius without solving the diffusion or transport equation. The resulting explicit formula is compared to known critical masses (few-percent agreement), extended to impure materials, and contrasted with the Oppenheimer-Bethe formula and a diffusion-equation treatment employing new boundary conditions.

Significance. If the central approximations hold, the work supplies a transparent, elementary route to criticality estimates that separates reaction physics from geometry. This is valuable for pedagogy, historical context, order-of-magnitude design checks, and cross-validation of Monte Carlo codes such as MCNP. The reported few-percent agreement with tabulated critical masses and the extension to impure fissile materials demonstrate practical reach for a simplified approach. The explicit formula and direct comparison to historical methods are strengths.

major comments (2)
  1. [Formulation of the threshold condition and derivation of critical radius] The separation of nuclear reaction rates from neutron transport geometry via total path length alone (invoked in the threshold condition for stationary neutron number) is load-bearing for the claim that no diffusion equation is needed. The average path length is an integral over the actual flux and angular distribution; the paper's auxiliary statistical model (random-walk or chord-length averaging) therefore requires explicit justification or independent validation against exact transport solutions to confirm that the few-percent agreement is not an artifact of the model choice.
  2. [Comparison with known critical masses and other calculations] In the comparison with known critical masses, the effective neutron multiplication per path length appears as a free parameter drawn from literature values. Clarify whether this parameter is held fixed from independent data or adjusted to reproduce the tabulated masses; if the latter, the validation partially re-uses the same information that informed the inputs and the independence of the test should be qualified.
minor comments (2)
  1. [Abstract] The abstract states that the analysis 'can be extended to general neutronics estimations'; a brief concrete example of such an extension would strengthen the claim.
  2. [Throughout the equations and text] Notation for mean free paths, cross sections, and the effective multiplication factor should be introduced once and used consistently; a short nomenclature table would aid readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments and the recommendation for minor revision. Our responses to the major comments are provided below, along with indications of how we plan to revise the manuscript.

read point-by-point responses

  1. Referee: [Formulation of the threshold condition and derivation of critical radius] The separation of nuclear reaction rates from neutron transport geometry via total path length alone (invoked in the threshold condition for stationary neutron number) is load-bearing for the claim that no diffusion equation is needed. The average path length is an integral over the actual flux and angular distribution; the paper's auxiliary statistical model (random-walk or chord-length averaging) therefore requires explicit justification or independent validation against exact transport solutions to confirm that the few-percent agreement is not an artifact of the model choice.

    Authors: We appreciate this insightful comment. Our pedagogical approach deliberately employs a simplified statistical model for the average neutron path length, based on chord-length averaging in a sphere, which yields an average path length of 4R/3. This is a standard result in integral geometry and is justified in the context of the elementary method presented. While we recognize that this does not capture the full details of the neutron flux and angular distribution as in exact transport theory, the method's strength lies in its transparency and the observed agreement with established critical masses at the few-percent level. To strengthen the manuscript, we will include additional justification referencing the mean chord length and discuss the approximation's scope and limitations relative to diffusion or Monte Carlo methods. revision: yes

  2. Referee: [Comparison with known critical masses and other calculations] In the comparison with known critical masses, the effective neutron multiplication per path length appears as a free parameter drawn from literature values. Clarify whether this parameter is held fixed from independent data or adjusted to reproduce the tabulated masses; if the latter, the validation partially re-uses the same information that informed the inputs and the independence of the test should be qualified.

    Authors: The effective neutron multiplication per path length is calculated using fixed nuclear parameters (such as fission cross-sections, absorption cross-sections, and average neutrons per fission) taken directly from independent literature sources, without any adjustment to match the critical mass data. The comparison with tabulated critical masses therefore represents an a posteriori validation of the model's predictive capability. We will revise the manuscript to explicitly state the origin of these parameters and emphasize their independence from the critical mass values being compared. revision: yes

Circularity Check

0 steps flagged

Derivation self-contained; comparison to external benchmarks is validation

full rationale

The paper derives an approximate critical-radius expression from the stationary-neutron-number condition using only elementary statistical averaging over total path length, with nuclear cross-sections and mean free paths taken from independent prior literature. The resulting formula is then compared to tabulated experimental critical masses for agreement at the few-percent level. This comparison constitutes an external check rather than a fit by construction, and no equation reduces the output to the inputs. The separation of reaction rates from geometry is presented as an ansatz whose accuracy is tested against known results, not assumed to be exact. No self-citation chain or uniqueness theorem is invoked to force the result.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The derivation assumes standard nuclear cross-sections and mean free paths are known inputs; the balance condition is enforced by setting average path length equal to a multiple of the fission mean free path, with the multiplier treated as an adjustable statistical factor.

free parameters (1)
  • effective neutron multiplication per path length
    Statistical factor chosen so that fission production balances absorption plus leakage; calibrated against known critical masses.
axioms (2)
  • domain assumption Neutron number is constant at criticality when average path length satisfies the production-loss balance
    Invoked in the threshold condition that opens the derivation.
  • domain assumption Nuclear reaction rates and transport geometry can be separated except for total path length
    Stated as the central methodological step.

pith-pipeline@v0.9.0 · 5771 in / 1374 out tokens · 62987 ms · 2026-05-19T11:36:43.197102+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.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    the threshold condition that the number of neutrons in the sphere does not change with time; the average neutron path length in the material must be long enough to produce enough fission neutrons to balance losses by absorption due to nuclear reactions and leakage through the surface. This separates the nuclear reaction part of the problem from the geometry and mechanics of neutron transport, the only connection being the total path length

  • IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    R_c = ε ℓ_s ℓ_c with ℓ_c = −1/(n σ0) ln(1 − σ0/(ν σ_f − σ0)) and ε ≈ √(10/(4π))

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.

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages

  1. [1]

    Fissions will occur, producing a number of neutrons Nf “ νN ˆ σf σ0 ˙ ` 1 ´ e´nσ 0ℓ˘

    work page

  2. [2]

    Neutrons will be absorbed, losing a number N0 “ N P0 “ N ` 1 ´ e´nσ 0ℓ˘

    work page

  3. [3]

    Iτw as 9N {N “ 1{τw. The presence of other processes does not affect this “leakage

    N Neutrons will leave the sphere, according to the definition of τw and ℓ. The threshold of criticality is defined by N being con- stant, we must have an ensemble average such that with each period τw, the number of neutrons produced by fis- sion equals the loss by emission (neutrons exiting through the sphere surface) and by absorption, νNf ´ N0 ´ N “ 0. (3...

    work page 2006

  4. [4]

    Nsℓs “ ℓspR2{ℓ2 sq “ R2{ℓs “ nσtrR2, (13) where ℓs is again the mean free path for scattering. The number of neutrons N produced per neutron is then N “ ˆ σf σ0 ν ´ 1 ˙

    72 and the fission threshold is quite high at 10 MeV. IV. NUMBER OF FIRST-GENERATION NEUTRONS PRODUCED AS A FUNCTION OF SPHERE RADIUS In this section we consider the number of first gener- ation neutrons that are produced and on average reach the mean distance ℓ0 to the edge of the sphere after a single neutron is placed randomly within the sphere. The path...

    work page

  5. [5]

    Ni η . (16) Studies of this effect were pioneered by Feynman et al. during the Manhattan projects[14] The fact that multi- plication measurements can provide a “nondestructive

    42, and NN “ 1. 25, which sum to N “ 10. As produc- tion in the sphere and leakage from the sphere are sepa- rate random processes that are coupled only through N , 7 there are independent chances of upward vs. downward fluctuations for both production and loss; they might sometimes coincide while sometimes going in the oppo- site direction. This implies t...

    work page

  6. [6]

    π ? 3 d 1 rpv ´ 1qn2σtrσf s ˆ f. The factor f is a correction derived from a “more accu- rate integral theory

    02 ˆ 1023, or Ng „ 72. The mean distance trav- eled between fissions is 1 {nσf “ 11 cm. (Scattering occurs over a smaller distance, but this does not mat- ter so far as computing a timescale is concerned.) The total distance traveled is consequently 72 ˆ 11 „ 800 cm for the uncompressed metal. The velocity of a 1 MeV neutron is 1 . 4 ˆ 109 cm/s, so the tim...

    work page 1943

  7. [7]

    Chyba and Caroline R

    Christopher F. Chyba and Caroline R. Milne, Simple calculation of the critical mass for highly enriched ura- nium and plutonium-239, American Journal of Physics 82, 977 (2014); doi: 10.1119/1.4885379 View on- line: https://doi.org/10.1119/1.4885379, and reference s therein

    work page doi:10.1119/1.4885379 2014

  8. [8]

    Derringh, Estimate of the critical mass of a fissionable isotope, Am

    E. Derringh, Estimate of the critical mass of a fissionable isotope, Am. J. Phys. 58, 363–364 (1990).https://doi.org/10.1119/1.16172

    work page doi:10.1119/1.16172 1990

  9. [9]

    B. C. Reed, The Physics of the Manhattan Project, 4th ed. (Springer, Heidelberg, 2021)

    work page 2021

  10. [10]

    Chadwick, Nuclear Science for the Manhattan Project and Comparison to Today’s ENDF Data , Nuclear Technology 207, Supplement 1, pp

    Mark B. Chadwick, Nuclear Science for the Manhattan Project and Comparison to Today’s ENDF Data , Nuclear Technology 207, Supplement 1, pp. S24-S61 (2021). DOI: https://doi.org/10.1080/00295450.2021.1901002

    work page doi:10.1080/00295450.2021.1901002 2021

  11. [11]

    K. M. Case and P. F. Zweifel, Linear Transport Theory (Addison-Wesley, New York, 1967). See also, J. J. Duder- stadt and W. R. Martin, Transport Theory (Wiley, New York, 1979)

    work page 1967

  12. [12]

    H. D. Seldby et al., Fission Product Data Measured at Los Alamos for Fission Spectrum and Thermal Neutrons on 239Pu, 235U, 238U , Nucl. Data Sheets 111, 2891 (2010); https://doi.org/10.1016/j.nds.2010.11.002

    work page doi:10.1016/j.nds.2010.11.002 2010

  13. [13]

    Chadwick, Private Commumication (2025)

    M.B. Chadwick, Private Commumication (2025)

    work page 2025

  14. [14]

    https://wwwndc.jaea.go.jp/cgi- bin/Tab80WWW.cgi?lib=J40&iso=U238

    work page

  15. [15]

    Neudecker, A.E

    D. Neudecker, A.E. Lovell, T. Kawano, P. Talou, Evaluated 238U(n,f ) Average Prompt Fission Neutron Multiplicities Including the CGMF Model , Los Alamos National Laboratory, Los Alamos, NM, 87545, USA, https://www.osti.gov/servlets/purl/1889944, Sept. 28, 2022

    work page arXiv 2022

  16. [16]

    Paxton and N.L

    H.C. Paxton and N.L. Pruvost, Critical Di- mensions of Systems Containing 235U, 239Pu, and 233U, Los Alamos Report LA-10860- MS (UC 46), July 1987. Also available at https://www.nuclearweaponarchive.org/Nwfaq/Nfaq4- 1.html#Nfaq4.1.7.1

    work page 1987

  17. [17]

    01 for A « 240, see A.M

    The elastic scattering causes a fractional energy loss due to recoil, of approximately 1 ´ E E1 « 1 ´ A ´ 1 A ` 1 “ 0. 01 for A « 240, see A.M. Weinberg and E.P. Wigner, The Physical Theory of Neutron Chain Reactors (University of Chicago Press, Chicago, 1958), Eq. (10.1a) and (10.6). For inelastic scattering, the incident neutron excites in- ternal state...

    work page 1958

  18. [18]

    Glaser, On the Proliferation Potential of Uranium Fuel for Research Reactors at Various Enrichment Lev- els, Science and Global Security 14(1) 1-24 (2006)

    A. Glaser, On the Proliferation Potential of Uranium Fuel for Research Reactors at Various Enrichment Lev- els, Science and Global Security 14(1) 1-24 (2006)

    work page 2006

  19. [19]

    doi:10.1007/978-0- 387-87857-7

    Hans Fischer, A History of the Central Limit Theorem: From Classical to Modern Probability Theory, Sources and Studies in the History of Mathematics and Physical Sciences (New York: Springer, 2011). doi:10.1007/978-0- 387-87857-7

    work page doi:10.1007/978-0- 2011

  20. [20]

    Peter Galison, Feynman ’s War: Modelling Weapons, Modelling Nature, Stud. Hist. Phil. Mod. Phys., Vol. 29, No. 3, 391—434, 1998

    work page 1998

  21. [21]

    Templin et al., Reactor Physics Constants: Ar- gonne National Laboratory Report ANL-5800 2nd ed

    L.J. Templin et al., Reactor Physics Constants: Ar- gonne National Laboratory Report ANL-5800 2nd ed. (1963) https://www.osti.gov/servlets/purl/4620873 DOI:https://doi.org/10.2172/4620873 Fast neutron reac- tions are described in Sec. 7. It should be noted that many parameters are given to 4 or 5 digits of precision; in re- ality, achieving better than 5%...

    work page doi:10.2172/4620873 1963

  22. [22]

    Johnson, S

    N. Johnson, S. Kotz, and N. Balakrishnan, Continuous Univariate Distributions, Vol. 1, 2nd ed. (Boston, MA: Houghton Mifflin, 1994) p. 421

    work page 1994

  23. [23]

    J. M. Pearson and B. C. Reed, Remarks on the yield of fission bombs , Am. J. Phys. 92(9) 680-685 (2024). 13

    work page 2024

  24. [24]

    Henry (https://math.stackexchange.com/users/6460 /henry), Expected Value of Random Walk, URL (version: 2012- 01-28): https://math.stackexchange.com/q/103170

    work page 2012

  25. [25]

    Johnson, T. (2024). The first indication of neutrino-induced nuclear fission (Order No. 31632293). Available from ProQuest Disserta- tions & Theses Global. (3156296008). Retrieved from https://www.proquest.com/dissertations- theses/first-indication-neutrino-induced-nuclear- fission/docview/3156296008/se-2; see also E. Kolbe, K. Langanke, and G. M. Fuller, Neu...

    work page arXiv 2024

  26. [26]

    (https://www.sciencedirect.com/science /article/pii/S0370269323001661)

    Yonghao Chen, et al., Measurement of neutron- induced fission cross sections of 232Th from 1 to 300 MeV at CSNS Back-n, Physics Let- ters B, 839, 137832 (2023) ISSN 0370-2693, https://doi.org/10.1016/j.physletb.2023.137832. (https://www.sciencedirect.com/science /article/pii/S0370269323001661)

    work page doi:10.1016/j.physletb.2023.137832 2023

  27. [27]

    https://wwwndc.jaea.go.jp/cgi- bin/Tab80WWW.cgi?/data/JENDL/JENDL-4- prc/intern/Th232.intern The number of prompt neutrons per fission is around 2.4, however, there is 10% uncertainty in this, compare A. B. Smith, R.G. Nobles and S.A. Cox, Phys. Rev. 115, 1242 (1959), and D.S. Mather, P. Fieldhouse, A. Moat, Nucl. Phys. 66, 149 (1965)

    work page 1959