Validation of Geant4 nuclear reaction models for hadrontherapy and preliminary results with SMF and BLOB

Reliable nuclear fragmentation models are of utmost importance in hadrontherapy, where Monte Carlo (MC) simulations are used to compute the input parameters of the treatment planning software, to validate the deposited dose calculation, to evaluate the biological effectiveness of the radiation, to correlate the bþ emitters production in the patient body with the delivered dose, and to allow a non- invasive treatment verification. Despite of its large use, the models implemented in Geant4 have shown severe limitations in reproducing the measured secondaries yields in ions interaction below 100 MeV/A, in term of production rates, angular and energy distributions [1–3]. We will present a benchmark of the Geant4 models with double-differential cross sec- tion and angular distributions of the secondary fragments produced in the 12C fragmentation at 62 MeV/A on thin carbon target, such a benchmark includes the recently implemented model INCL++ [4,5]. Moreover, we will present the preliminary results, obtained in simulating the same interaction, with SMF [6] and BLOB [7]. Both, SMF and BLOB are semiclassical one-body approaches to solve the Boltzmann-Langevin equation. They include an identical treatment of the mean-field propagation, on the basis of the same effective interaction, but they differ in the way fluctuations are included. In particular, while SMF employs a Uehling-Uhlenbeck collision term and introduces fluctuations as projected on the density space, BLOB introduces fluctuations in full phase space through a modified collision term where nucleon-nucleon correlations are explicitly involved. Both of them, SMF and BLOB, have been developed to sim- ulate the heavy ion interactions in the Fermi-energy regime. We will show their capabilities in describing 12C fragmentation foreseen their implementation in Geant4.