Integrated Simulation of Fragmentation, Evaporation, and Gamma-decay Processes in the Interaction of Cosmic-ray Heavy Ions with the Atmosphere using PHITS


General-purpose Monte-Carlo radiation tranport calculation code PHITS is applied to calculate prompt gamma-ray emission from cosmic-ray heavy ions fragmented in the atmosphere. Event-by-event simulation of spallation reactions by cosmic-ray heavy ions was performed by combination of three reaction models, responsible for different reaction phases.

[1] S. Agostinelli et al., “Geant4-a simulation toolkit”, Nuclear Instruments and Methods in Physics Research A, 506, 250-303 (2003)

[2] F. Ballarini et al., “The physics of the FLUKA code: Recent development”, Advances in Space Research, 40, 1339-1349 (2002).

[3] Laurie S. Waters et al., “MCNPX User’s Manual”, LA-UR, 02-2607 (2002).

[4] T.Sato et al., “Particle and Heavy Ion Transport code System, PHITS, version 2.52”, Journal of Nuclear Science and Technology, 50, 913-923 (2013)

[5] T.Wilcox et al,“Correlated gammas using CGM and MCNPX,”Progress in Nuclear Energy, 63, 1-6(2013)

[6] A. Ferrari et al., “The production of residual nuclei in peripheral high energy nucleusnucleus interactions”, Zeitschrift fur Physik C, 71, 75-86 (1996).

[7] T. Ogawa et al., “Development of gamma de-excitation model for prediction of prompt gamma-rays and isomer production based on energy-dependent level structure treatment”, Nuclear Instruments and Methods in Physics Research A, 325,35- 42 (2014)

[8] K. Niita et al., ”Analysis of the (N,xN’) reactions by quantum molecular dynamics plus statistical decay model”, Physical Review C, 52, 2620-2635 (1995).

[9] T.Ogawa et al., ”Energy-dependent fragmentation cross sections of relativistic 12C”, Physical Review C, 92, 024614 (2015)

[10] S. Furihata, ”Statistical analysis of light fragment production from medium energyproton-induced reactions”, Nuclear Instruments and Methods in Physics Research Section B, 171, 3, 251-258 (2000).

[11] Jag Tuli et al., “Evaluated Nuclear Structure Data File”,

[12] J.Aichelin et al., ”Quantum Molecular Dynamics”, PHYSICS REPORTS, 202, 233-360 (1991)

[13] V.F. Weisskopf, “Statistics and Nuclear Reactions”, Physical Review, 52, 295 (1937).

[14] V. F. Weisskopf et al., Physical Review, 57, 472 (1940)

[15] A. Gilbert et al., “A composit nuclear-level density formula with shell corrections”, Canadian Journal of Physics, 43, 1446-1496 (1965).

[16] I. Dostrovsky et al., “Monte Carlo Calculations of Nuclear Evaporation Processes. III. Applications to Low-Energy Reactions”, Physical Review, 116, 683 (1959)