Comprehensive Analysis of Particle Production, Transport, and Spatial Dose Distribution in Two Multilayer Shielding Designs for a 50 MeV Proton Accelerator Using FLUKA Monte Carlo Simulations
Abstract
This study investigates secondary particle production and transport processes in two different multilayer shielding configurations (Design 1 and Design 2) for a 50 MeV proton accelerator using the FLUKA Monte Carlo (MC) code. A monoenergetic proton beam interacting with a copper target generates a complex radiation field composed of photons, neutrons, protons, and light nuclear particles. Analysis of inelastic interactions shows that photon production is dominant in both designs at approximately 50%, followed by neutrons (~27%) and protons (~17%), indicating the significant role of electromagnetic processes in the low-energy regime. Stopping particle statistics reveal that nearly 99% of primary protons are absorbed within the shielding structures, while the contribution of light ions remains negligible. The investigation of low-energy neutron interactions demonstrates a strong neutron multiplication behavior, with more than 97% of the generated secondary particles being neutrons. Although both shielding designs exhibit similar overall particle production characteristics, Design 1 yields a slightly higher neutron multiplication factor (~4.16 particles per proton) compared to Design 2 (~3.74 particles per proton), suggesting that Design 2 provides relatively more effective moderation for low-energy neutrons. These findings indicate that while initial particle production is governed by electromagnetic interactions, the evolution of the radiation field is primarily controlled by neutron transport and multiplication mechanisms. Furthermore, one-dimensional ambient dose equivalent H*(10) profiles along the x and y directions clearly demonstrate the directional attenuation performance of the shielding configurations. In addition, relative error distributions obtained from MC simulations were analyzed to assess computational reliability, confirming that the majority of dose results fall within acceptable uncertainty limits (generally below 10%). Overall, the study provides a detailed physical interpretation of particle production mechanisms and dose distribution characteristics, offering valuable insights for the optimization of multilayer shielding designs in low-energy proton accelerator environments.
Keywords
Proton accelerator, Multilayer shielding, Secondary particle production, Neutron transport, Monte Carlo simulation
Supporting Institution
Ethical Statement
Thanks
References
- Agosteo, S., Magistris, M., Mereghetti, A., Silari, M., & Zajacova, Z. (2007). Shielding data for 100–250 MeV proton accelerators: Double differential neutron distributions and attenuation in concrete. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 265(2), 581–598. https://doi.org/10.1016/j.nimb.2007.09.046
- Agosteo, S., Magistris, M., Mereghetti, A., Silari, M., & Zajacova, Z. (2008). Shielding data for 100–250 MeV proton acceleratos: Attenuation of secondary radiation in thick iron and concrete/iron shields. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 266(15), 3406–3416. https://doi.org/10.1016/j.nimb.2008.05.002
- Alhazmi, N. K., Alanazi, S. D., Alrowily, H. S. A., Alkhamsan, K. F., Al-Rashidi, F. A. D., Mohammed, A. S. S., Khalaf, A. M., Alshammari, S. F. N., Alzaidi, A. A., Alshammari, F. T. N., & AlQusayri, A. A. K. (2024). Radiation safety compliance and protective practices among radiologic technologists. International Journal of Computational and Experimental Science and Engineering, 10(4), 3569–3577. https://doi.org/10.22399/ijcesen.4664
- Aydın, H., Süsoy Doğan, G., Tekin, H. O., Şen Baykal, D., & İltuş, Y. C. (2025). Examining the neutron and gamma attenuation characteristics of various amorphous structures with bioactive properties. International Journal of Computational and Experimental Science and Engineering, 11(3), 6973–6986. https://doi.org/10.22399/ijcesen.2176
- Battistoni, G., Boehlen, T., Cerutti, F., Chin, P. W., Esposito, L. S., Fassò, A., Ferrari, A., Lechner, A., Empl, A., Mairani, A., Mereghetti, A., Ortega, P. G., Ranft, J., Roesler, S., Sala, P. R., Vlachoudis, V., & Smirnov, G. (2015). Overview of the FLUKA code. Annals of Nuclear Energy, 82, 10–18. https://doi.org/10.1016/j.anucene.2014.11.007
- Böhlen, T. T., Cerutti, F., Chin, M. P. W., Fassò, A., Ferrari, A., Ortega, P. G., Mairani, A., Sala, P. R., Smirnov, G., & Vlachoudis, V. (2014). The FLUKA code: Developments and challenges for high energy and medical applications. Nuclear Data Sheets, 120, 211–214. https://doi.org/10.1016/j.nds.2014.07.049
- De Saint-Hubert, M., Farah, J., Klodowska, M., Romero-Expósito, M. T., Tyminska, K., Mares, V., Olko, P., Stolarczyk, L., & Trinkl, S. (2022). The influence of nuclear models and Monte Carlo radiation transport codes on stray neutron dose estimations in proton therapy. Radiation Measurements, 150, Article 106693. https://doi.org/10.1016/j.radmeas.2021.106693
- Emikönel, S., & Akkurt, İ. (2025). Radiation shielding properties of B₂O₃-Bi₂O₃ glass. International Journal of Computational and Experimental Science and Engineering, 11(2), 3461–3464. https://doi.org/10.22399/ijcesen.2157
- Ferrari, A., Sala, P. R., Fassò, A., & Ranft, J. (2005). FLUKA: A multi-particle transport code (program version 2005) (CERN-2005-010). CERN. https://doi.org/10.5170/CERN-2005-010
- Gökcan, A. O., & Sarıyer, D. (2026). High-energy neutron shielding analysis using curve fitting and machine learning: A Monte Carlo-based study in concrete and Ferro-Boron configurations. Annals of Nuclear Energy, 232, Article 112244. https://doi.org/10.1016/j.anucene.2026.112244