pith. sign in

arxiv: 1907.04323 · v1 · pith:4MBBBY64new · submitted 2019-07-09 · ⚛️ physics.ins-det

Monte Carlo calculation and verification of the geometrical factors for the NPDGamma experiment

Pith reviewed 2026-05-25 00:02 UTC · model grok-4.3

classification ⚛️ physics.ins-det
keywords NPDGammaMonte Carlogeometrical factorsparity violationneutron captureCsI detectorsMCNPX
0
0 comments X

The pith

Monte Carlo modeling supplies per-detector corrections to the geometrical factors for the NPDGamma parity-violation measurement.

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

The NPDGamma experiment extracts a small parity-violating asymmetry from gamma rays emitted after polarized neutrons capture on liquid parahydrogen. Earlier work approximated detector sensitivities with a step-wise sinusoidal function that assumed cylindrical symmetry around the target. Real mechanical deviations break that symmetry, so the paper replaces the approximation with a Monte Carlo calculation that gives an individual sensitivity value for each cesium iodide detector. The calculation uses source-modified MCNPX and is checked by measuring the known parity-violating asymmetry from neutron capture on chlorine.

Core claim

The Monte Carlo model computes the sensitivity of each detector in the array to the physics asymmetry, providing corrections that replace the earlier cylindrical approximation, and these corrections are validated by reproducing the measured chlorine asymmetry.

What carries the argument

Source-modified MCNPX Monte Carlo simulation that tracks the response of each individual CsI detector to the parity-violating asymmetry signal.

If this is right

  • The extracted physics asymmetry for the parahydrogen target uses a unique sensitivity weight for every detector rather than a single functional form.
  • The same Monte Carlo geometry can be re-run for any change in detector placement or target position without re-deriving analytic approximations.
  • The chlorine measurement serves as an in-situ calibration that transfers directly to the parahydrogen data set.

Where Pith is reading between the lines

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

  • Future parity-violation experiments with large detector arrays may need similar Monte Carlo corrections whenever mechanical tolerances break perfect symmetry.
  • The validation approach of using a secondary target with a known asymmetry could be repeated for other neutron-capture targets to cross-check the model.

Load-bearing premise

The MCNPX model with the added source modifications accurately reproduces the physical geometry, materials, and neutron interactions of the actual NPDGamma apparatus.

What would settle it

A statistically significant mismatch between the Monte Carlo prediction and the measured parity-violating asymmetry for the chlorine target would show that the model does not correctly represent the detector responses.

Figures

Figures reproduced from arXiv: 1907.04323 by David Blyth, Elise Tang, Geoffrey L. Greene, James D. Bowman, Kyle B. Grammer, Matthew Musgrave, Nadia Fomin, Zhaowen Tang.

Figure 1
Figure 1. Figure 1: Neutrons enter the RFSR from the left after passing [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Detector numbering scheme with the neutron prop [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: GUD, GLR, and GZ ideal geometrical factors for the 35Cl target inside the downstream end of the RFSR. The three directionally weighted energy deposition tallies are calculated identically to a normal , F6 en￾ergy deposition tally but the energy deposition tally at each scattering event is weighted by the initial γ￾ray direction tags using the tallyx subroutine, T j x = 1 Nn X Nn k=1 X l E j kl( ˆkγ · xˆ)k,… view at source ↗
Figure 4
Figure 4. Figure 4: Grid scan data points using the 4 mCi 137Cs source are shown in blue. The MCNPX model grid points are shown in red. scan was simulated using MCNPX using a cylindri￾cal γ-ray source with energy 662 keV in order to de￾termine the ideal detector response. The grid scan patterns are shown in figure 4. The MCNPX grid represents the ideal response of the detectors to the cesium source and would differ from the s… view at source ↗
Figure 5
Figure 5. Figure 5: Extracted δφ correction angles with uncertainties determined from a fit of the measured detector response for a 137Cs source to the MCNPX calculated detector response using equation 22. order polynomial in x and y using GNU Scientific Library [14], P(x, y) = X 8 p=0 X 8−i q=0 kp,qx p y q . (22) The laboratory coordinate system differs from the model coordinate system with a rotation by δφ,  x ′ y ′  =  … view at source ↗
Figure 6
Figure 6. Figure 6: 35Cl asymmetry fitted to ring-specific geometrical factors (ie. amplitude for each ring varying with cos(θ)) and scaled by the resulting scaling parameter, APV, in equation 27. The reduced χ 2 is 0.88. 0 10 20 30 40 −1.0 −0.5 0.0 0.5 1.0 Detector Geometric factor A j raw/APV,const GPV [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: 35Cl asymmetry fitted to constant amplitude ge￾ometrical factors for each ring (analogous to the LANSCE procedure [8]) and scaled by the resulting scaling parameter, APV,const, in equation 27. The reduced χ 2 is 2.57. raw asymmetry and an uncertainty for each detector, which was then fit to the geometrical factors using A j raw = APCG j PC ′ + APVG j PV ′ . (27) The χ 2 per degree of freedom from this fit … view at source ↗
Figure 10
Figure 10. Figure 10: Ratio of the PV (top) and PC (bottom) geometrical [PITH_FULL_IMAGE:figures/full_fig_p009_10.png] view at source ↗
Figure 9
Figure 9. Figure 9: GLR and GPC ideal geometrical factors for the liquid hydrogen target. The liquid hydrogen geometrical factors are shown in figures 8 and 9. The neutron scattering correction to GLR is responsible for the significant difference be￾tween GLR and GPC. The orthohydrogen scattering cross section is a fac￾tor of 102 higher in magnitude than the absorption and parahydrogen scattering cross sections [10][19]. 0.99… view at source ↗
read the original abstract

The NPDGamma experiment measures the parity-violating asymmetry in $\gamma$-ray emission in the capture of polarized neutrons on liquid parahydrogen. The sensitivity to the asymmetry for each detector in the array is used as a parameter in the extraction of the physics asymmetry from the measured data. The detector array is approximately cylindrically symmetric around the target and a step-wise sinusoidal function has been used for the sensitivity in the previous iteration of the NPDGamma experiment, but deviations from cylindrical symmetry necessitate the use of a Monte Carlo model to determine corrections to the geometrical factors. For the calculations, source code modifications to MCNPX were done in order to calculate the sensitivity of each cesium iodide detector to the physics asymmetry. We describe the MCNPX model and results from calculations and how the results are validated through measurement of the parity violating asymmetry of $\gamma$-rays from neutron capture on chlorine.

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

0 major / 2 minor

Summary. The paper describes Monte Carlo calculations with MCNPX (including source code modifications) to compute per-detector sensitivities to the parity-violating asymmetry for the CsI array in the NPDGamma experiment. These sensitivities serve as parameters in extracting the physics asymmetry from polarized neutron capture on liquid parahydrogen. The authors argue that deviations from cylindrical symmetry around the target require a full Monte Carlo treatment rather than the step-wise sinusoidal approximation used previously, and they validate the model by comparing simulated and measured asymmetries from neutron capture on chlorine.

Significance. If the model and validation hold, the work supplies corrected geometrical factors that improve the precision of asymmetry extraction in a flagship parity-violation measurement. The chlorine asymmetry comparison provides an independent empirical anchor outside the parahydrogen data set itself. The approach of targeted MCNPX modifications to tally detector-specific responses is a standard but necessary engineering step for non-ideal geometries; documenting it with validation strengthens the reliability of the NPDGamma analysis pipeline.

minor comments (2)
  1. [Abstract / §3] The abstract states that source code modifications were made to MCNPX but does not specify the nature or location of those changes (e.g., which tally or source routine was altered). A brief description or pseudocode in §3 would allow readers to assess reproducibility.
  2. [Abstract] No quantitative comparison (e.g., χ² or fractional difference) between the Monte Carlo chlorine asymmetry and the measured value is provided in the abstract or summary; the validation is described only qualitatively.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading and positive assessment of the manuscript, including recognition of the significance of the MCNPX-based corrections to the geometrical factors and the independent validation via the chlorine asymmetry measurement. We appreciate the recommendation for minor revision.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper's central workflow is a standard Monte Carlo computation (MCNPX with targeted source modifications) of per-detector geometrical factors for the NPDGamma array, followed by empirical validation against a measured chlorine parity-violating asymmetry. This validation step is independent of the parahydrogen target result and supplies external grounding. No self-definitional equations, fitted-input predictions, load-bearing self-citations, or ansatz smuggling are present in the described derivation chain; the claim that cylindrical symmetry deviations require MC treatment follows directly from the geometry without circular reduction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based on abstract only; the central claim rests on the validity of the Monte Carlo simulation and the chlorine validation measurement. No explicit free parameters or invented entities are mentioned.

axioms (1)
  • domain assumption The MCNPX code with modifications accurately models neutron capture, gamma emission and detector response.
    Invoked implicitly as the basis for calculating sensitivities.

pith-pipeline@v0.9.0 · 5705 in / 1054 out tokens · 24860 ms · 2026-05-25T00:02:09.747517+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

19 extracted references · 19 canonical work pages

  1. [1]

    Fomin, G

    N. Fomin, G. Greene, R. Allen, V. Cianciolo, C. Crawford, T. Tito, P. Huffman, E. Iver- son, R. Mahurin, W. Snow, Fundamental neu- tron physics beamline at the spallation neutron source at ORNL, Nucl. Instrum. Methods Phys. Res., Sect. A 773 (2015) 45–51

  2. [2]

    Balascuta, R

    S. Balascuta, R. Alarcon, S. Baeßler, G. Greene, A. Mietke, C. Crawford, R. Milburn, S. Penttil¨ a, J. Prince, J. Sch¨ adler, S. Bae β ler, G. Greene, A. Mietke, C. Crawford, R. Milburn, S. Pent- tila, J. Prince, J. Sch¨ adler, Others, The imple- mentation of a super mirror polarizer at the SNS fundamental neutron physics beamline, Nuclear Instruments and...

  3. [3]

    P.-N. Seo, L. Barr´ on-Palos, J. Bowman, T. Chupp, C. Crawford, M. Dabaghyan, M. Dawkins, S. Freedman, T. Gentile, M. Ger- icke, R. C. Gillis, G. Greene, F. Hersman, G. Jones, M. Kandes, S. Lamoreaux, B. Lauss, M. Leuschner, R. Mahurin, M. Mason, J. Mei, G. Mitchell, H. Nann, S. Page, S. Pent- til¨ a, W. Ramsay, A. Salas-Bacci, S. Santra, M. Sharma, T. Sm...

  4. [4]

    Gericke, C

    M. Gericke, C. Blessinger, J. Bowman, R. C. Gillis, J. Hartfield, T. Ino, M. Leuschner, Y. Ma- suda, G. Mitchell, S. Muto, H. Nann, S. Page, S. Penttil¨ a, W. Ramsay, P.-N. Seo, W. Snow, J. Tasson, W. Wilburn, A current mode detector array for γ ray asymmetry measurements, Nu- clear Instruments and Methods in Physics Re- search, Section A: Accelerators, Sp...

  5. [5]

    Santra, L

    S. Santra, L. Barr´ on-Palos, C. Blessinger, J. Bowman, T. Chupp, S. Covrig, C. Crawford, M. Dabaghyan, J. Dadras, M. Dawkins, M. Ger- icke, W. Fox, R. Gillis, M. Leuschner, B. Lo- zowski, R. Mahurin, M. Mason, J. Mei, H. Nann, S. Penttil¨ a, A. Salas-Bacci, M. Sharma, W. Snow, W. Wilburn, L. Barr´ on Palos, C. Blessinger, J. Bowman, T. Chupp, S. Cov- rig...

  6. [6]

    J. Fry, R. Alarcon, R. Allen, E. Askanazi, S. Bal- ascuta, L. Barron-Palos, S. Baeßler, A. Barzilov, C. Blessinger, D. Blyth, J. D. Bowman, J. R. Calarco, T. E. Chupp, C. E. Coppola, C. Craw- ford, K. Craycraft, M. Dabaghyan, D. Evans, J. Favela, C. Fieseler, N. Fomin, W. Fox, S. Freedman, E. Frleˇ z, C. Fu, C. Garcia, I. Gar- ishvili, M. T. Gericke, R. C...

  7. [7]

    Musgrave, S

    M. Musgrave, S. Baeßler, S. Balascuta, L. Barr´ on-Palos, D. Blyth, J. Bowman, T. Chupp, V. Cianciolo, C. Crawford, K. Cray- craft, N. Fomin, J. Fry, M. Gericke, R. Gillis, K. Grammer, G. Greene, J. Hamblen, C. Hayes, P. Huffman, C. Jiang, S. Kucuker, M. McCrea, P. Mueller, S. Penttil¨ a, W. Snow, E. Tang, Z. Tang, X. Tong, W. Wilburn, Measurement of the a...

  8. [8]

    Gericke, R

    M. Gericke, R. Alarcon, S. Balascuta, L. Barr´ on- Palos, C. Blessinger, J. Bowman, R. Carlini, W. Chen, T. Chupp, C. Crawford, S. Covrig, M. Dabaghyan, N. Fomin, S. Freedman, T. Gen- tile, R. Gillis, G. Greene, F. Hersman, T. Ino, G. Jones, B. Lauss, M. Leuschner, W. Lozowski, R. Mahurin, Y. Masuda, J. Mei, G. Mitchell, S. Muto, H. Nann, S. Page, S. Pent...

  9. [9]

    D. B. Pelowitz, MCNPX User’s Manual, Los Alamos National Laboratory, 2011

  10. [10]

    Chadwick, M

    M. Chadwick, M. Herman, P. Obloˇ zinsk´ y, M. Dunn, Y. Danon, A. Kahler, D. Smith, B. Pritychenko, G. Arbanas, R. Arcilla, R. Brewer, D. A. Brown, R. Capote, A. D. Carl- son, Y. S. Cho, H. Derrien, K. Guber, G. M. Hale, S. Hoblit, S. Holloway, T. D. Johnson, T. Kawano, B. C. Kiedrowski, H. Kim, S. Ku- nieda, N. M. Larson, L. Leal, J. P. Lestone, R. C. Lit...

  11. [11]

    Berger, J

    M. Berger, J. Coursey, M. Zucker, J. Chang, STAR, PSTAR, and ASTAR: Computer Pro- grams for Calculating Stopping-Power and Range Tables for Electrons, Protons, and He- lium Ions (version 1.2.3) (2005)

  12. [12]

    Gericke, J

    M. Gericke, J. Bowman, M. Johnson, Mott- Schwinger scattering of polarized low energy neutrons up to thermal energies, Phys. Rev. C 78 (4) (2008) 044003

  13. [13]

    Csoto, B

    A. Csoto, B. F. Gibson, G. L. Payne, Parity con- serving gamma asymmetry in n-p radiative cap- ture, Phys. Rev. C 56 (2) (1997) 631–634

  14. [14]

    Galassi, J

    M. Galassi, J. Davies, J. Theiler, B. Gough, G. Jungman, M. Booth, F. Rossi, GNU Scien- tific Library Reference Manual, Vol. 954161734, 2009

  15. [15]

    Fomin, First results from the NPDGamma experiment at the Spallation Neutron Source, AIP Conf

    N. Fomin, First results from the NPDGamma experiment at the Spallation Neutron Source, AIP Conf. Proc. 1560 (2013) 145–148

  16. [16]

    Mitchell, C

    G. Mitchell, C. Blessinger, J. Bowman, T. Chupp, K. Coulter, M. Gericke, G. Jones, M. Leuschner, H. Nann, S. Page, S. Penttil¨ a, T. Smith, W. Snow, W. Wilburn, A measure- ment of parity-violating gamma-ray asymme- tries in polarized cold neutron capture on 35Cl, 113Cd, and 139La, Nuclear Instruments and Methods in Physics Research Section A: Acceler- ato...

  17. [17]

    Barr´ on-Palos, R

    L. Barr´ on-Palos, R. Alarcon, S. Balascuta, C. Blessinger, J. Bowman, T. Chupp, S. Covrig, C. Crawford, M. Dabaghyan, J. Dadras, M. Dawkins, W. Fox, M. Gericke, R. Gillis, B. Lauss, M. Leuschner, B. Lozowski, R. Mahurin, M. Mason, J. Mei, H. Nann, S. Penttil¨ a, W. Ramsay, A. Salas-Bacci, S. Santra, P.-N. Seo, M. Sharma, T. Smith, W. Snow, W. Wilburn, V....

  18. [18]

    K. B. Grammer, R. Alarcon, L. Barr´ on-Palos, D. Blyth, J. D. Bowman, J. Calarco, C. Craw- ford, K. Craycraft, D. Evans, N. Fomin, J. Fry, M. Gericke, R. C. Gillis, G. L. Greene, J. Ham- blen, C. Hayes, S. Kucuker, R. Mahurin, M. Maldonado-Vel´ azquez, E. Martin, M. Mc- Crea, P. E. Mueller, M. Musgrave, H. Nann, S. I. Penttil¨ a, W. M. Snow, Z. Tang, W. S...

  19. [19]

    Mughabghab, Atlas of Neutron Resonances, 5th Edition, Elsevier, New York, 2006

    S. Mughabghab, Atlas of Neutron Resonances, 5th Edition, Elsevier, New York, 2006. 12