KnE Energy | The 3rd International Conference on Particle Physics and Astrophysics (ICPPA) | pages: 391–398

, , , , and

1. Introduction

The GeV-class gamma-rays from astronomical objects such as blazers and pulsars are observed by satellites (e.g. Fermi-LAT) or telescope arrays to study cosmic-ray acceleration mechanism, evolution of galaxies, etc. On the other hand, the cosmic ray heavy ions interacting with the nuclei in the atmosphere produce various secondary particles such as nucleons, nuclear clusters, fragments, mesons, and prompt gamma-rays. The gamma-rays from energetic projectile fragments are boosted by Doppler Effect and observed as GeV-range gamma-rays in the laboratory frame. The energetic gamma-rays from cosmic-ray heavy ions should be considered in the experiments.

Three different calculation steps, dynamic phase, evaporation and prompt gamma-decay, are necessary to simulate production of such high energy gamma-rays (Fig.1). It is also vital to simulate reactions on event-by-event basis to deduce the excitation energy of each fragment. A number of simulation codes have been developed for prediction of cosmic ray reactions; however, very few can simulate event-by-event prompt gamma-ray emission[1-4] because it is necessary to determine the isotopic species produced in each reaction and to calculate the internal transition intensity between discrete levels based on the nuclear structure data [5-7].

The general purpose radiation transport code PHITS (Particle and Heavy Ion Transport code System) incorporates a heavy ion reaction model JQMD (JAERI Quantum Molecular Dynamics Model) [8,9], statistical decay model GEM (Generalized Evaporation Model) [10], and prompt gamma-ray production model EBITEM (ENSDF[11]-Based Isomeric Transition and isomEr production Model) [7] and therefore PHITS is capable of simulating prompt gamma-rays from heavy ion fragments.

Figure 1

Calculation scheme of prompt gamma-rays from cosmic heavy ion fragments.


Prompt gamma-ray emission was simulated in typical conditions such as 10 AGeV Oxygen interactions with nitrogen and 1 AGeV boron interactions with nitrogen. The calculated gamma-ray energy spectra showed sharp peaks corresponding to particular internal transition paths.

2. Method

In this study, cosmic ray transport and interaction simulation is performed by using PHITS. PHITS is a general-purpose radiation transport simulation code which can simulate transport and reactions of nearly all kinds of particles from thermal energies to 1 TeV. In PHITS, the 3 reaction phases (dynamic phase, evaporation and prompt gamma-decay) were simulated by the JQMD, GEM, and EBITEM, respectively. Because all these models are event generators, secondary particle emission was calculated in event-by-event basis. The three models are briefly explained below.

Dynamic phase model JQMD

JQMD is one of the reaction event generators based on the quantum molecular dynamics (QMD) approach [12], which is often used for simulation of nucleus-nucleus reaction dynamic phase. In comparison with other QMD models, JQMD can simulate various perspectives, such as secondary particle production and fragment yields, by using fewer parameters.

In JQMD, nucleons follow a Lorentz-covariant equation of motion. The long range interaction between nucleons is described by the potential,


where A is one of the Skyrme force constant (-219.4MeV), ρ s is the saturation density (0.168 fm -3 ), τ is another Skyrme force constant (4/3), B is the other Skyrme force constant (165.3 MeV), ρ i is the local nucleon density at the position of i-th particle, i is the spatial coordinate of the centroid of the i-th particle, c i is charge of i-th particle, e is elementary charge, L is the width of wave packet (2 fm 2 ), C s is the symmetry force constant (25 MeV) and ρ ij is the wave function overlap of i-th particle and j-th particle. Strong repulsion attributed to the hard core is not included in this potential. Instead, the nucleon pairs close to each other undergo elastic or inelastic scattering based on their cross sections. Production of baryon resonance particles such as Δand N* as well as absorption of pions by nucleons are also considered.

After 150fm/c of time evolution, nucleons close to each other are bound to form clusters. Mass, charge, kinetic energy, and excitation energy are then calculated by summing the contribution from the constituent nucleons. More details on JQMD are available elsewhere [7,8].

Evaporation model GEM

GEM is one of the event generators to simulate statistical decay of excited nuclei. Nuclei undergo sequential particle emission based on the Weisskopf-Ewing approach [13,14] in GEM. Particle emission probability Γis calculated by,

σinvε=whiteπRb2c11+bε for neutron whiteπRb2cj1Vε for others

where subscript i denotes initial nucleus, subscript j denotes ejectile particle, and subscript d denotes daughter nucleus. S j is ejectile spin, m j is ejectile mass,ρ d is level density of daughter nucleus calculated by Fermi gas model or by the model proposed by Gilbert and Cameron [15], Q is Q value, V is Coulomb potential between the projectile and the ejectile, R b , c 1 , c j , and b are model parameters given in [16], respectively.

The initial state of statistical decay is specified by the clustering phase of JQMD. The mass, charge, kinetic energy, and excitation energy are passed to GEM.

Prompt gamma-decay model EBITEM

A theoretical model to simulate gamma de-excitation of excited nuclei, EBITEM (ENSDF-Based Isomeric Transition and isomEr production Model), is based on the Evaluated Nuclear Structure Data File (ENSDF), and theories [15]. When excitation energy became smaller than the nucleon separation energy, evaporation by GEM was stopped and the residue, characterized by mass, charge, angular momentum and excitation energy, starts gamma de-excitation. If the excitation energy was smaller than 3 MeV or the nuclear mass was smaller than 40 amu, the level structure taken from ENSDF was employed. Otherwise, the level density was calculated by the Gilbert-Cameron formula expressed in Eq. (3),

ρU,Jπa14U54 Exp (2aU)2J+12σ2πσ2 Exp (J(J+1)2σ2),

where U is excitation energy, a is level density parameter, J is angular momentum, and σ is spin cut-off parameter.

Figure 2

Secondary particle trajectories in one reaction.


The probability of transition from one state to another state was taken from ENSDF if the nucleus was at an excitation state given in ENSDF; otherwise it was calculated by the single-particle model. The transition probability was calculated as a function of total angular momentum transfer and gamma-ray energy. By sampling the de-excitation path at random, the next level characterized by excitation energy and total angular momentum was determined. By repeating this process corresponding to de-excitation, the nuclei reach either the ground state or the metastable states.

The detailed description of the theory and validation of EBITEM is provided elsewhere [7].

3. Results and Discussion

Some case studies on the application of PHITS are discussed.

  • Event-by-event analysis

Fig. 2 shows particle trajectories after an interaction of 12 C with 14 N at 100 AGeV in the atmosphere. The air density averaged from the ground level to 10 km in the sky was adopted as the air density in this calculation. Various particles such as muons, pions, photons, neutrons, and neutrinos are seen as particle tracks. High energy gamma-rays are also observed in this event. In the laboratory frame, two photons with energies 410 MeV and 646 MeV were detected in coincidence and their emission angle was shifted by 0.17 degrees from each other. This gamma emission event is attributed to internal transition of 12 C emitting 3.12 MeV (7.56MeV - > 4.44MeV) and 4.44 MeV (4.44MeV- > Ground) gamma-rays. Gamma-ray energy in the laboratory frame is shifted by Doppler effect taking into account for the projectile energy (100 AGeV) and the angle between gamma emission and projectile momentum.

Figure 3

Energy spectrum of gamma-rays from 11B interacting with nitrogen.


In this way, PHITS can be used to reproduce one event starting from one heavy ion interacting with the target nucleus.

  • Average over events

Fig. 3 shows energy spectrum of gamma-rays from the fragments produced by 10 AGeV 11 B interacting with nitrogen. The gamma-ray energies are calculated in the rest frame of the projectile.

Gamma-ray peaks attributed to various fragment species are seen in the spectrum. It is possible to simulate correlation between the fragment and its gamma-ray. In reality, this peak structure is blurred because the kinetic energy of fragments is distributed and gamma-ray energies are shifted by Doppler Effect.

4. Conclusion

In this study, it is shown that PHITS is capable of simulating interactions of cosmic-ray heavy ions. In particular, for predicting prompt gamma-rays from the fragments of cosmic-ray heavy ions, combination of a non-equilibrium reaction model, a statistical decay model and a prompt-gamma decay model is essential.

In addition, needful functions and data such as transformation between the frames, event-by-event correlation and nuclear level structure are implemented to help users to obtain desired quantities in useful conditions.

These facts indicate that PHITS is one of the most robust tools to plan and analyze cosmic ray experiments. In particular, PHITS offers a unique capability to simulate prompt gamma-rays from cosmic ray heavy ion fragments.


We express our gratitude to the operation team of the Center for Computational Science and E-systems (CCSE) at the Japan Atomic Energy Agency. The simulations reported in this paper were executed on a PC-cluster system located at CCSE. This work was supported by JSPS Grants-in-Aid for Scientific Research (Number 26790072).



Geant4 Collaboration Geant4##hssm###8212;a simulation toolkit Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 20035063250303


F. Ballarini The physics of the FLUKA code: Recent development , Advances in Space Research 200240Recent developmentThe physics of the FLUKA code


S. Laurie Waters. MCNPX Users Manual MCNPX Users Manual 2002022607


T. Sato K. Niita N. Matsuda S. Hashimoto Y. Iwamoto S. Noda T. Ogawa H. Iwase H. Nakashima T. Fukahori K. Okumura T. Kai S. Chiba T. Furuta L. Sihver Particle and heavy ion transport code system, PHITS, version 2.52 Journal of Nuclear Science and Technology 20135099139232-s2.0-8488244550010.1080/00223131.2013.814553


T. Wilcox T. Kawano G. W. McKinney J. S. Hendricks Correlated gammas using CGM and MCNPX Progress in Nuclear Energy 201363162-s2.0-8487016433110.1016/j.pnucene.2012.10.002


A. Ferrari J. Ranft S. Roesler P. R. Sala The production of residual nuclei in peripheral high energy nucleus-nucleus interactions Zeitschrift fur Physik C-Particles and Fields 199671175862-s2.0-3374912292110.1007/s002880050149


T. Ogawa S. Hashimoto T. Sato K. Niita 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 Section B: Beam Interactions with Materials and Atoms 201432535422-s2.0-8489767766210.1016/j.nimb.2014.02.007


K. Niita S. Chiba T. Maruyama T. Maruyama H. Takada T. Fukahori Y. Nakahara A. Iwamoto Analysis of the (N,xN) reactions by quantum molecular dynamics plus statistical decay model Physical Review C: Nuclear Physics 1995525262026352-s2.0-000082635410.1103/PhysRevC.52.2620


T. Ogawa T. Sato S. Hashimoto D. Satoh S. Tsuda K. Niita Erratum: Energy-dependent fragmentation cross sections of relativistic C 12 (Physical Review C - Nuclear Physics 92, 024614 (2015)) Physical Review C nuclear physics 20159222-s2.0-8494085467410.1103/PhysRevC.92.029904029904


S. Furihata Statistical analysis of light fragment production from medium energy proton-induced reactions Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 200017132512582-s2.0-003432088510.1016/S0168-583X(00)00332-3


Jag. Tuli Evaluated Nuclear Structure Data File,


J. Aichelin ‘Quantum’ molecular dynamics−a dynamical microscopic n-body approach to investigate fragment formation and the nuclear equation of state in heavy ion collisions Physics Reports 199120223336010.1016/0370-1573(91)90094-3


V. Weisskopf Statistics and nuclear reactions Physical Review A: Atomic, Molecular and Optical Physics 19375242953032-s2.0-3614901880210.1103/PhysRev.52.295


V. F. Weisskopf D. H. Ewing On the yield of nuclear reactions with heavy elements Physical Review A: Atomic, Molecular and Optical Physics 19405764724852-s2.0-1884439947410.1103/PhysRev.57.472


A. Gilbert A composit nuclear-level density formula with shell corrections Canadian Journal of Physics 1965431446149610.1139/p65-139


I. Dostrovsky Z. Fraenkel G. Friedlander Monte Carlo Calculations of Nuclear Evaporation Processes. III. Applications to Low-Energy Reactions Physical Review A: Atomic, Molecular and Optical Physics 195911636837022-s2.0-3614901683510.1103/PhysRev.116.683



  • Downloads 15
  • Views 80



ISSN: 2413-5453