FLUKA simulation of target fragmentation in proton therapy.

In proton therapy, secondary fragments are created in nuclear interactions of the beam with the target nuclei. The secondary fragments have low kinetic energies and high atomic numbers as compared to primary protons. Fragments have a high LET and deposit all their energy close to the generation point. For their characteristics, secondary fragments can alter the dose distribution and lead to an increase of RBE for the same delivered physical dose. Moreover, the radiobiological impact of target fragmentation is significant mostly in the region before the Bragg peak, where generally healthy tissues are present, and immediately after Bragg peak. Considering the high biological impact of those particles, especially in the case of healthy tissues or organs at risk, the inclusion of target fragmentation processes in the dose calculation of a treatment planning system can be relevant to improve the treatment accuracy and for this reason it is one of the major tasks of the MoVe IT project. In this study, Monte Carlo simulations were employed to fully characterize the mixed radiation field generated by target fragmentation in proton therapy. The dose averaged LET has been evaluated in case of a Spread Out Bragg Peak (SOBP). Starting from LET distribution, RBE has been evaluated with two different phenomenological models. In order to characterize the mixed radiation field, the production cross section has been evaluated by means of the FLUKA code. The future development of present work is to generate a MC database of fragments fluence to be included in TPS.

On Molière and Fermi-Eyges scattering theories in hadrontherapy.

We consider the evaluation of lateral spread distributions of charged particle beams at therapeutic energies, due to an absorber in the form of a homogeneous slab or of a stack. We show that the Molière theory has the same degree of flexibility as the Fermi-Eyges, but is much more accurate and does not present particular computing difficulties with the energy loss formula we have employed. It is also shown that the non-Gaussian shape of the projected one dimensional (1D) distributions is not a complication for passing from the projected to the spatial two-dimensional (2D) distribution, if one assumes circular symmetry. All the calculations are compared with the results of the FLUKA code. The nuclear interaction is not considered here, because it is outside of the scope of this work.

A model for the accurate computation of the lateral scattering of protons in water.

A pencil beam model for the calculation of the lateral scattering in water of protons for any therapeutic energy and depth is presented. It is based on the full Molière theory, taking into account the energy loss and the effects of mixtures and compounds. Concerning the electromagnetic part, the model has no free parameters and is in very good agreement with the FLUKA Monte Carlo (MC) code. The effects of the nuclear interactions are parametrized with a two-parameter tail function, adjusted on MC data calculated with FLUKA. The model, after the convolution with the beam and the detector response, is in agreement with recent proton data in water from HIT. The model gives results with the same accuracy of the MC codes based on Molière theory, with a much shorter computing time.