⁹Be and ¹²C Thick Targets Neutron Yields from Particle Accelerator-Induced Nuclear Reactions by PHITS Code

We calculated the double-differential cross-sections and neutron yields for the ⁹Be(⁴He, 𝑛)¹²C and ¹²C(𝑑, 𝑛)¹³N reactions at various energies using the PHITS code, which employs a Monte Carlo method. Our results were validated against existing experimental data, showing good agreement up to 30 MeV. Based on our findings, we identified optimal reactions at specific energies that maximize neutron yields, highlighting their potential for enhancing radiotherapy applications.