Date of Award


Degree Type


Degree Name

Doctor of Philosophy


Nuclear Engineering

Major Professor

Lawrence H. Heilbronn

Committee Members

Lawrence Townsend, Jason Hayward, Yanwen Zhang


The study of secondary particles produced from heavy-ion interactions is important in heavy ion radiotherapy, space radiation protection, and shielding at accelerator facilities. This dissertation focuses on the study of secondary neutron production as they are of special concern among all secondary particles.

The first part of this dissertation is the measurement of secondary neutrons created from 4He [helium] stopped in various target materials together with the model calculations accomplished by PHITS, FLUKA, and MCNP transport codes. The comparison results show that the physics models need improvements particularly in the predictions of 1) neutrons created from the 4He interactions at the high-energy end of the spectra at each angle for FLUKA’s and PHITS’s models, 2) the high-energy peaks at 0degree for all systems and all models, and 3) the low-energy neutrons at small angles for 230-MeV/nucleon [megaelectron volt per nucleon] 4He stopping in the light targets. However, the model calculations agree with the experiment data well at intermediate and large angles in intermediate and low energy regions.

The second part is the benchmark of the neutron production cross section data with model calculations fulfilled by PHITS, FLUKA, and MCNP. The studied cases cover wide ranges of projectile species, beam energy and target nuclei mass. Some significant differences do appear not only among model calculations but also between measured data and calculations. In particular, LAQGSM03.03 implemented in MCNP6 significantly overestimates the high-energy peak in the forward direction in the light and very light system at 400 MeV/nucleon. RQMD implemented in FLUKA 2011.2c overestimates the neutron cross sections at intermediate energies in nearly all systems expect the lightest targets in our studies cases. The greatest inter-model difference appear on low-energy neutrons at forward angles in the system of 400-MeV/nucleon 132Xe (xenon) and copper target, and it is inferred that GEM implemented in PHITS 2.73 over-predicts neutrons produced from evaporation.

The results of both experimental study and model calculations provide critical information for validation and verification of the current radiation transport codes used for simulating heavy-ion interactions and help lead to improvements in the physics models.


Figure 44 in Appendix 1 is a duplicate of Figure 43 due to an editing mistake. The data for Figure 44 can be referred to in the second column of Figure 39. The MCNP6 calculation for the 400-MeV/nucleon carbon ions on the lithium target, which is shown in the first column of Figure 39 and Figure 43 in Appendix 1, is in error due to a mistake in the input to the code. The corrected MCNP calculations now show better agreement to the data.

Files over 3MB may be slow to open. For best results, right-click and select "save as..."