Benchmarking GEANT4 nuclear models for hadron therapy with 95 MeV/nucleon carbon ions

In carbon-therapy, the interaction of the incoming beam with human tissues may lead to the production of a large amount of nuclear fragments and secondary light particles. An accurate estimation of the biological dose deposited into the tumor and the surrounding healthy tissues thus requires sophisticated simulation tools based on nuclear reaction models. The validity of such models requires intensive comparisons with as many sets of experimental data as possible. Up to now, a rather limited set of double di erential carbon fragmentation cross sections have been measured in the energy range used in hadrontherapy (up to 400 MeV/A). However, new data have been recently obtained at intermediate energy (95 MeV/A). The aim of this work is to compare the reaction models embedded in the GEANT4 Monte Carlo toolkit with these new data. The strengths and weaknesses of each tested model, i.e. G4BinaryLightIonReaction, G4QMDReaction and INCL++, coupled to two di fferent de-excitation models, i.e. the generalized evaporation model and the Fermi break-up are discussed

Simulation study on light ions identification methods for carbon beams from 95 to 400 MeV/A

Monte Carlo simulations have been performed in order to evaluate the efficiencies of several light ions identification techniques. The detection system was composed with layers of scintillating material to measure either the deposited energy or the time-of-flight of ions produced by nuclear reactions between 12C projectiles and a PMMA target. Well known techniques such as (DELTA) E--Range, (DELTA) E--E--ToF and (DELTA)E--E are presented and their particle identification efficiencies are compared one to another regarding the generated charge and mass of the particle to be identified. The simulations allowed to change the beam energy matching the ones proposed in an hadron therapy facility, namely from 95 to 400 MeV/A

Limitation of energy deposition in classical N body dynamics

Energy transfers in collisions between classical clusters are studied with Classical N Body Dynamics calculations for different entrance channels. It is shown that the energy per particle transferred to thermalised classical clusters does not exceed the energy of the least bound particle in the cluster in its ``ground state''. This limitation is observed during the whole time of the collision, except for the heaviest system.Comment: 13 pages, 15 figures, 1 tabl

Double di ffential fragmentation cross sections measurements of 95 MeV/u 12C on thin targets for hadrontherapy

During therapeutic treatment with heavy ions like carbon, the beam undergoes nuclear fragmentation and secondary light charged particles, in particular protons and alpha particles, are produced. To estimate the dose deposited into the tumors and the surrounding healthy tissues, an accurate prediction on the fluences of these secondary fragments is necessary. Nowadays, a very limited set of double di ffential carbon fragmentation cross sections are being measured in the energy range used in hadrontherapy (40 to 400 MeV/u). Therefore, new measurements are performed to determine the double di ffential cross section of carbon on di erent thin targets. This work describes the experimental results of an experiment performed on May 2011 at GANIL. The double di ffential cross sections and the angular distributions of secondary fragments produced in the 12C fragmentation at 95 MeV/u on thin targets (C, CH2, Al, Al2O3, Ti and PMMA) have been measured. The experimental setup will be precisely described, the systematic error study will be explained and all the experimental data will be presented.Comment: Submitted to PR

New approach of fragment charge correlations in 129Xe+(nat)Sn central collisions

A previous analysis of the charge (Z) correlations in the $\Delta Z-$ plane for Xe+Sn central collisions at 32 MeV/u has shown an enhancement in the production of equally sized fragments (low $\Delta Z$) which was interpreted as an evidence for spinodal decomposition. However the signal is weak and rises the question of the estimation of the uncorrelated yield. After a critical analysis of its robustness, we propose in this paper a new technique to build the uncorrelated yield in the charge correlation function. The application of this method to Xe+Sn central collision data at 32, 39, 45 and 50 MeV/u does not show any particular enhancement of the correlation function in any $\Delta Z$ bin.Comment: 23 pages, 9 figures, revised version with an added figure and minor changes. To appear in Nuclear Physics

Mid-rapidity charge distribution in peripheral heavy ion collisions

The charge density distribution with respect to the velocity of matter produced in peripheral heavy ion reactions around Fermi energy is investigated. The experimental finding of enhancement of mid-rapidity matter shows the necessity to include correlations beyond BUU which was performed in the framework of nonlocal kinetic theory. Different theoretical improvements are discussed. While the in-medium cross section changes the number of collisions, it leads the transferred energy almost unchanged. In contrast the nonlocal scenario changes the energy transferred during collisions and leads to an enhancement of mid-rapidity matter. The renormalisation of quasiparticle energies is shown to be possible to include in nonlocal scenarios and and leads to a further enhancement of mid-rapidity matter distribution. This renormalisation is accompanied by a dynamical softening of the equation of state seen in longer oscillation periods of the excited compressional collective mode. We propose to include quasiparticle renormalization by using the Pauli-rejected collisions which circumvent the problem of backflows in Landau theory. Using the maximum relative velocity of projectile and target like fragments we associate experimental events with impact parameters of the simulations. For peripheral collisions we find a reasonable agreement between experiment and theory. For more central collisions the velocity damping is higher in one - body simulations than observed experimentally which is due to missing cluster formations in the used kinetic theory

Disappearance of Transverse Flow in Central Collisions for Heavier Nuclei

For the first time, mass dependence of balance energy only for heavier systems has been studied. Our results are in excellent agreement with the data which allow us to predict the balance energy of U+U, for the first time, around 37-39 MeV/nucleon. Also our results indicate a hard equation of state along with nucleon-nucleon cross-section around 40 mb.Comment: 4 pages, 3 figures, submitted to Phys. Rev. Let

Analysis of anisotropic flow with Lee-Yang zeroes

We present a new method to extract anisotropic flow in heavy ion collisions from the genuine correlation among a large number of particles. Anisotropic flow is obtained from the zeroes in the complex plane of a generating function of azimuthal correlations, in close analogy with the theory of phase transitions by Lee and Yang. Flow is first estimated globally, i.e., averaged over the phase space covered by the detector, and then differentially, as a function of transverse momentum and rapidity for identified particles. The corresponding estimates are less biased by nonflow correlations than with any other method. The practical implementation of the method is rather straightforward. Furthermore, it automatically takes into account most corrections due to azimuthal anisotropies in the detector acceptance. The main limitation of the method is statistical errors, which can be significantly larger than with the ``standard'' method of flow analysis if the flow and/or the event multiplicities are too small. In practice, we expect this to be the most accurate method to analyze directed and elliptic flow in fixed-target heavy-ion collisions between 100 MeV and 10 GeV per nucleon (at the Darmstadt SIS synchrotron and the Brookhaven Alternating Gradient Synchrotron), and elliptic flow at ultrarelativistic energies (at the Brookhaven Relativistic Heavy Ion Collider, and the forthcoming Large Hadron Collider at CERN).Comment: 32 pages, 7 eps figures, RevTe