Determination of fast neutron RBE using a fully mechanistic computational model.

This work presents a model previously developed for estimating relative biological effectiveness (RBE) associated with high-LET particles. It is based on the combination of Monte Carlo simulations of particle interactions when traversing an atomic resolution DNA geometrical model. In addition, the model emulates the induction of lethal damage from the interaction of two sublethal lesions, taken as double-strand breaks. The Geant4-DNA package was used for simulations with liquid water as the transport medium. The RBE of neutron beams with energies ranging from 0.1 MeV up to 14 MeV was studied. The model succeeded in reproducing the general behavior of RBE as a function of neutron energy, including the RBE peak reported by experiments at approximately 0.4 MeV. Furthermore, the results of the model agree rather well with some experimental works. However, our results underestimate RBE for neutron energies above approximately 5 MeV due to the current limitations of Geant4-DNA for the tracking of heavy ions below 0.5 MeV/u.

Computational approach to determine the relative biological effectiveness of fast neutrons using the Geant4-DNA toolkit and a DNA atomic model from the Protein Data Bank.

This study proposes an innovative approach to estimate relative biological effectiveness (RBE) of fast neutrons using the Geant4 toolkit. The Geant4-DNA version cannot track heavy ions below 0.5 MeV/nucleon. In order to explore the impact of this issue, secondary particles are simulated instead of the primary low-energy neutrons. The Evaluated Nuclear Data File library is used to determine the cross sections for the elastic and inelastic interactions of neutrons with water and to find the contribution of each secondary particle spectrum. Two strategies are investigated in order to find the best possible approach and results. The first one takes into account only light particles, protons produced from elastic scattering, and α particles from inelastic scattering. Geantino particles are shot instead of heavy ions; hence all heavy ions are considered in the simulations, though their physical effects on DNA not. The second strategy takes into account all the heavy and light ions, although heavy ions cannot be tracked down to very low energies (E<0.5 MeV/nucleon). Our model is based on the combination of an atomic resolution DNA geometrical model and a Monte Carlo simulation toolkit for tracking particles. The atomic coordinates of the DNA double helix are extracted from the Protein Data Bank. Since secondary particle spectra are used instead of simulating the interaction of neutrons explicitly, this method reduces the computation times dramatically. Double-strand break induction is used as the end point for the estimation of the RBE of fast neutrons. ^{60}Co Î³ rays are used as the reference radiation quality. Both strategies succeed in reproducing the behavior of the RBE_{max} as a function of the incident neutron energy ranging from 0.1 to 14 MeV, including the position of its peak. A comparison of the behavior of the two strategies shows that for neutrons with energies less than 0.7 MeV, the effect of heavy ions would not be very significant, but above 0.7 MeV, heavy ions have an important role in neutron RBE.