Use of short-lived positron emitters for in-beam and real-time ß+ range monitoring in proton therapy.

AIM: The purpose of this work is to evaluate the precision with which the GEANT4 toolkit simulates the production of ß+ emitters relevant for in-beam and real-time PET in proton therapy. BACKGROUND: An important evolution in proton therapy is the implementation of in-beam and real-time verification of the range of protons by measuring the correlation between the activity of ß+ and dose deposition. For that purpose, it is important that the simulation of the various ß+ emitters be sufficiently realistic, in particular for the 12N short-lived emitter that is required for efficient in-beam and real-time monitoring. METHODS: The GEANT4 toolkit was used to simulate positron emitter production for a proton beam of 55 MeV in a cubic PMMA target and results are compared to experimental data. RESULTS: The three ß+ emitters with the highest production rates in the experimental data (11C, 15O and 12N) are also those with the highest production rate in the simulation. Production rates differ by 8% to 174%. For the 12N isotope, the ß+ spatial distribution in the simulation shows major deviations from the data. The effect of the long range (of the order of 20 mm) of the ß+ originating from 12N is also shown and discussed. CONCLUSIONS: At first order, the GEANT4 simulation of the ß+ activity presents significant deviations from the data. The need for precise cross-section measurements versus energy below 30 MeV is of first priority in order to evaluate the feasibility of in-beam and real-time PET.

A simulation study of gold nanoparticles localisation effects on radiation enhancement at the mitochondrion scale.

This paper presents a Monte-Carlo study focusing on the effects of gold nanoparticles on the energy deposition patterns produced by incident photons in the close vicinity of the mitochondrial network modeled as a tube. Spherical shaped gold nanoparticles of 30 nm diameter were placed in a micrometric (10 × 10 × 10 µm3) water phantom containing a tube of 300 nm diameter and 5 µm length. The tube represented a mitochondrial fragment and nanoparticles were distributed in the water phantom outside the tube. Photons of 120 keV were simulated using the Geant4 Livermore processes and the Geant4-DNA electron processes to account for secondary electrons collisions. The Livermore processes took into account the Auger cascade inside the gold material. A data mining algorithm was then used to analyze the energy deposition clusters inside the water phantom and the tube. A comparison was made between the results obtained for a uniform distribution of nanoparticles and a vesicle distribution model. The results including energy deposition clusters are also compared to dose enhancement ratios.

Microdosimetry in 3D realistic mitochondria phantoms: Geant4 Monte Carlo tracking of 250keV photons in phantoms reconstructed from microscopic images.

Mitochondria are considered to be sensitive radiation targets since they control processes vital to the cell's functioning. These organelles are starting to get attention and some studies are investigating the radiation dose inside them. In previous studies, mitochondria are represented as simple ellipsoids inside the cell not taking into consideration the complexity of their shape. In this study, realistic phantoms are built based on deconvolved widefield fluorescent microscopic images of the mitochondrial networks of fibroblast cells. The phantoms are imported into Geant4 as tessellated volumes taking into account the geometrical complexity of these organelles. Irradiation with 250keV photons is performed and the lineal energy is calculated. The lineal energy distributions inside the produced phantoms are compared with those calculated inside simple volumes, a sphere and an ellipsoid, where the effect of the shape and volume is clearly seen on lineal energies.

Measurement of carbon ion microdosimetric distributions with ultrathin 3D silicon diodes.

The commissioning of an ion beam for hadrontherapy requires the evaluation of the biologically weighted effective dose that results from the microdosimetric properties of the therapy beam. The spectra of the energy imparted at cellular and sub-cellular scales are fundamental to the determination of the biological effect of the beam. These magnitudes are related to the microdosimetric distributions of the ion beam at different points along the beam path. This work is dedicated to the measurement of microdosimetric spectra at several depths in the central axis of a (12)C beam with an energy of 94.98 AMeV using a novel 3D ultrathin silicon diode detector. Data is compared with Monte Carlo calculations providing an excellent agreement (deviations are less than 2% for the most probable lineal energy value) up to the Bragg peak. The results show the feasibility to determine with high precision the lineal energy transfer spectrum of a hadrontherapy beam with these silicon devices.

Distributions of secondary particles in proton and carbon-ion therapy: a comparison between GATE/Geant4 and FLUKA Monte Carlo codes.

Monte Carlo simulations play a crucial role for in-vivo treatment monitoring based on PET and prompt gamma imaging in proton and carbon-ion therapies. The accuracy of the nuclear fragmentation models implemented in these codes might affect the quality of the treatment verification. In this paper, we investigate the nuclear models implemented in GATE/Geant4 and FLUKA by comparing the angular and energy distributions of secondary particles exiting a homogeneous target of PMMA. Comparison results were restricted to fragmentation of (16)O and (12)C. Despite the very simple target and set-up, substantial discrepancies were observed between the two codes. For instance, the number of high energy (>1 MeV) prompt gammas exiting the target was about twice as large with GATE/Geant4 than with FLUKA both for proton and carbon ion beams. Such differences were not observed for the predicted annihilation photon production yields, for which ratios of 1.09 and 1.20 were obtained between GATE and FLUKA for the proton beam and the carbon ion beam, respectively. For neutrons and protons, discrepancies from 14% (exiting protons-carbon ion beam) to 57% (exiting neutrons-proton beam) have been identified in production yields as well as in the energy spectra for neutrons.

In-beam quality assurance using induced ß(+) activity in hadrontherapy: a preliminary physical requirements study using Geant4.

Light and heavy ions particle therapy, mainly by means of protons and carbon ions, represents an advantageous treatment modality for deep-seated and/or radioresistant tumours. An in-beam quality assurance principle is based on the detection of secondary particles induced by nuclear fragmentations between projectile and target nuclei. Three different strategies are currently under investigation: prompt γ rays imaging, proton interaction vertex imaging and in-beam positron emission tomography. Geant4 simulations have been performed first in order to assess the accuracy of some hadronic models to reproduce experimental data. Two different kinds of data have been considered: ß(+)-emitting isotopes and prompt γ-ray production rates. On the one hand simulations reproduce experimental ß(+) emitting isotopes production rates to an accuracy of 24%. Moreover simulated ß(+) emitting nuclei production rate as a function of depth reproduce well the peak-to-plateau ratio of experimental data. On the other hand by tuning the tolerance factor of the photon evaporation model available in Geant4, we reduce significantly prompt γ-ray production rates until a very good agreement is reached with experimental data. Then we have estimated the total amount of induced annihilation photons and prompt γ rays for a simple treatment plan of ∼1 physical Gy in a homogenous equivalent soft tissue tumour (6 cm depth, 4 cm radius and 2 cm wide). The average annihilation photons emitted during a 45 s irradiation in a 4 π solid angle are ∼2 × 10(6) annihilation photon pairs and 10(8) single prompt γ whose energy ranges from a few keV to 10 MeV.

Molecular scale track structure simulations in liquid water using the Geant4-DNA Monte-Carlo processes.

This paper presents a study of energy deposits induced by ionising particles in liquid water at the molecular scale. Particles track structures were generated using the Geant4-DNA processes of the Geant4 Monte-Carlo toolkit. These processes cover electrons (0.025 eV-1 MeV), protons (1 keV-100 MeV), hydrogen atoms (1 keV-100 MeV) and alpha particles (10 keV-40 MeV) including their different charge states. Electron ranges and lineal energies for protons were calculated in nanometric and micrometric volumes.

Proton production in relativistic heavy ion collisions; comparison with a thermodynamical model.

Experimental results concerning proton production in nuclear collisions, obtained at Saturne with the Diogene 4 pi facility, are compared with the predictions of a thermodynamical model, using collective velocity distributions combined with a statistical thermodynamics in local rest frames. Experimental differential cross sections for alpha + nucleus and Neon + nucleus central collisions at incident energies between 200 and 800 MeV per nucleon are well reproduced by the model, for an angular range 30-110 degrees in the laboratory system. Extracted values of the temperatures are compared with those given by other authors.

Nuclear collective flow and charged-pion emission in Ne-nucleus collisions at E/A = 800 MeV.

Triple-differential cross sections of charged pions were measured for collisions of Ne projectiles at E/A = 800 MeV with NaF, Nb, and Pb targets. The reaction plane was estimated event by event from the light-baryon momentum distribution. For heavy targets, preferential emission of charged pions away from the interaction zone towards the projectile side was observed in the transverse direction. Such a preferential emission, which is not predicted by cascade calculations, may be attributed to a stronger pion absorption by the heavier spectator remnant.

Proton-proton correlations at small relative momentum in neon-nucleus collisions at E/A=400 and 800 MeV.

Proton-proton small angle correlations have been measured in neon-nucleus collisions, using the 4 pi detector Diogene, at 400 and 800 MeV per nucleon incident energies. Values of the size of the emitting region are obtained by comparison with the Koonin formula, taking into account the biases of the apparatus. The dependence of the density on target mass and incident energy is also analysed.

Exclusive measurements of mean pion multiplicities in 4He-nucleus reactions from 200 to 800 MeV/nucleon.

Mean multiplicities of pi+ and pi- in 4He collisions with C, Cu, and Pb at 200, 600, and 800 MeV/u, and with C and Pb at 400 MeV/u have been measured using the large solid angle detector Diogene. The independence of pion multiplicity on projectile incident energy, target mass and proton multiplicity is studied in comparison with intra-nuclear cascade predictions. The discrepancy between experimental results and theory is pointed out and discussed.

The Diogene 4 pi detector at Saturne.

Diogene, an electronic 4 pi detector, has been built and installed at the Saturne synchrotron in Saclay. The forward angular range (0 degree-6 degrees) is covered by 48 time-of-flight scintillator telescopes that provide charge identification. The trajectories of fragments emitted at larger angles are recorded in a cylindrical 0.4-m3 Pictorial Drift Chamber (PDC) surrounding the target. The PDC is inside a 1-T magnetic field; the axis of the PDC cylinder and the magnetic field are parallel to the beam. Good identification has been obtained for both positive and negative pi mesons and for hydrogen and helium isotopes. Multiplicities in relativistic nucleus-nucleus reactions up to 40 have been detected, limited mainly by the present electronics.