N. Freud

Monitoring the Bragg peak location of 73MeV∕u carbon ions by means of prompt γ-ray measurements

By means of a time-of-flight technique, we measured the longitudinal profile of prompt γ-rays emitted by 73MeV∕u C13 ions irradiating a polymethyl methacrylate target. This technique allowed us to minimize the shielding against neutrons and scattered γ-rays, and to correlate prompt gamma emission to the ion path. This correlation, together with a high counting rate, paves the way toward real-time monitoring of the longitudinal dose profile during ion therapy treatments. Moreover, the time correlation between the prompt gamma detection and the transverse position of the incident ions measured by a beam monitor can provide real-time three dimensional control of the irradiation.

A computer code to simulate X-ray imaging techniques

Abstract A computer code was developed to simulate the operation of radiographic, radioscopic or tomographic devices. The simulation is based on ray-tracing techniques and on the X-ray attenuation law. The use of computer-aided drawing (CAD) models enables simulations to be carried out with complex three-dimensional (3D) objects and the geometry of every component of the imaging chain, from the source to the detector, can be defined. Geometric unsharpness, for example, can be easily taken into account, even in complex configurations. Automatic translations or rotations of the object can be performed to simulate radioscopic or tomographic image acquisition. Simulations can be carried out with monochromatic or polychromatic beam spectra. This feature enables, for example, the beam hardening phenomenon to be dealt with or dual energy imaging techniques to be studied. The simulation principle is completely deterministic and consequently the computed images present no photon noise. Nevertheless, the variance of the signal associated with each pixel of the detector can be determined, which enables contrast-to-noise ratio (CNR) maps to be computed, in order to predict quantitatively the detectability of defects in the inspected object. The CNR is a relevant indicator for optimizing the experimental parameters. This paper provides several examples of simulated images that illustrate some of the rich possibilities offered by our software. Depending on the simulation type, the computation time order of magnitude can vary from 0.1 s (simple radiographic projection) up to several hours (3D tomography) on a PC, with a 400 MHz microprocessor. Our simulation tool proves to be useful in developing new specific applications, in choosing the most suitable components when designing a new testing chain, and in saving time by reducing the number of experimental tests.

Design Guidelines for a Double Scattering Compton Camera for Prompt-

In hadrontherapy in order to fully take advantage of the assets of the ion irradiation, the position of the Bragg peak has to be monitored accurately. Here, we investigate a monitoring method relying on the detection in real time of the prompt γ emitted quasi instantaneously during the nuclear fragmentation processes. Our detection system combines a beam hodoscope and a double scattering Compton camera. The prompt-γ emission points are reconstructed by intersecting the ion trajectories given by the hodoscope and the Compton cones reconstructed with the camera. We propose here to study in terms of point spread function and efficiency the theoretical feasibility of the emission points reconstruction with our set-up in the case of a photon point source in air. First we analyze the nature of all the interactions which are likely to produce an energy deposit in the three detectors of the camera. It is underlined that upper energy thresholds in both scatter detectors are required in order to select mainly Compton events (one Compton interaction in each scatter detector and one interaction in the absorber detector). Then, we study the influence of various parameters such as the photon energy and the inter-detector distances on the Compton camera response. These studies are carried out by means of Geant4 simulations. We use a source with a spectrum corresponding to the prompt-γ spectrum emitted during the carbon ion irradiation of a water phantom. In the current configuration, the spatial resolution of the Compton camera is about 6 mm (Full Width at Half Maximum) and the detection efficiency 10-5. Finally, provided the detection efficiency is increased, the clinical applicability of our system is considered.

\gamma

PURPOSE Proton CT (pCT) has the potential to accurately measure the electron density map of tissues at low doses but the spatial resolution is prohibitive if the curved paths of protons in matter is not accounted for. The authors propose to account for an estimate of the most likely path of protons in a filtered backprojection (FBP) reconstruction algorithm. METHODS The energy loss of protons is first binned in several proton radiographs at different distances to the proton source to exploit the depth-dependency of the estimate of the most likely path. This process is named the distance-driven binning. A voxel-specific backprojection is then used to select the adequate radiograph in the distance-driven binning in order to propagate in the pCT image the best achievable spatial resolution in proton radiographs. The improvement in spatial resolution is demonstrated using Monte Carlo simulations of resolution phantoms. RESULTS The spatial resolution in the distance-driven binning depended on the distance of the objects from the source and was optimal in the binned radiograph corresponding to that distance. The spatial resolution in the reconstructed pCT images decreased with the depth in the scanned object but it was always better than previous FBP algorithms assuming straight line paths. In a water cylinder with 20 cm diameter, the observed range of spatial resolutions was 0.7 - 1.6 mm compared to 1.0 - 2.4 mm at best with a straight line path assumption. The improvement was strongly enhanced in shorter 200° scans. CONCLUSIONS Improved spatial resolution was obtained in pCT images with filtered backprojection reconstruction using most likely path estimates of protons. The improvement in spatial resolution combined with the practicality of FBP algorithms compared to iterative reconstruction algorithms makes this new algorithm a candidate of choice for clinical pCT.

Imaging During Ion Beam Therapy: A Monte Carlo Simulation Study

Prompt-gamma profile was measured at WPE-Essen using 160 MeV protons impinging a movable PMMA target. A single collimated detector was used with time-of-flight (TOF) to reduce the background due to neutrons. The target entrance rise and the Bragg peak falloff retrieval precision was determined as a function of incident proton number by a fitting procedure using independent data sets. Assuming improved sensitivity of this camera design by using a greater number of detectors, retrieval precisions of 1 to 2 mm (rms) are expected for a clinical pencil beam. TOF improves the contrast-to-noise ratio and the performance of the method significantly.

Filtered backprojection proton CT reconstruction along most likely paths

Monte Carlo simulations are nowadays essential tools for a wide range of research topics in the field of radiotherapy. They also play an important role in the effort to develop a real-time monitoring system for quality assurance in proton and carbon ion therapy, by means of prompt-gamma detection. The internal theoretical nuclear models of Monte Carlo simulation toolkits are of decisive importance for the accurate description of neutral or charged particle emission, produced by nuclear interactions between beam particles and target nuclei. We assess the performance of Geant4 nuclear models in the context of prompt-gamma emission, comparing them with experimental data from proton and carbon ion beams. As has been shown in the past and further indicated in our study, the prompt-gamma yields are consistently overestimated by Geant4 by a factor of about 100% to 200% over an energy range from 80 to 310 MeV/u for the case of (12)C, and to a lesser extent for 160 MeV protons. Furthermore, we focus on the quantum molecular dynamics (QMD) modeling of ion-ion collisions, in order to optimize its description of light nuclei, which are abundant in the human body and mainly anticipated in hadrontherapy applications. The optimization has been performed by benchmarking QMD free parameters with well established nuclear properties. In addition, we study the effect of this optimization on charged particle emission. With the usage of the proposed parameter values, discrepancies reduce to less than 70%, with the highest values being attributed to the nucleon-ion induced prompt-gammas. This conclusion, also confirmed by the disagreement we observe in the case of proton beams, indicates the need for further investigation on nuclear models which describe proton and neutron induced nuclear reactions.

Real-time proton beam range monitoring by means of prompt-gamma detection with a collimated camera

In hadron therapy in order to fully take advantage of the assets of the ion irradiation, the position of the Bragg peak has to be monitored accurately. Here, we propose a monitoring method relying on the detection in real time of the prompt γ emitted quasi instantaneously during the nuclear fragmentation processes. Our detection system combines a beam hodoscope and a double scattering Compton camera. The prompt γ emission points are reconstructed by intersecting the ion trajectories given by the hodoscope and the Compton cones reconstructed with the camera. We studied the influence of various parameters such as the photon energy and the inter-detector distances on the Compton camera response to a photon point source. This study was carried out by means of Geant4 simulations. In the current configuration, for a photon source with a typical prompt γ spectrum, the spatial resolution of the Compton camera is about 5.6 mm and the detection efficiency 10-5.

Assessment and improvements of Geant4 hadronic models in the context of prompt-gamma hadrontherapy monitoring

Medical imaging and radiation therapy are widely used synchrotron-based techniques which have one thing in common: a significant dose delivery to typically biological samples. Among the ways to provide the experimenters with image guidance techniques indicating optimization strategies, Monte Carlo simulation has become the gold standard for accurately predicting radiation dose levels under specific irradiation conditions. A highly important hampering factor of this method is, however, its slow statistical convergence. A track length estimator (TLE) module has been coded and implemented for the first time in the open-source Monte Carlo code GATE/Geant4. Results obtained with the module and the procedures used to validate them are presented. A database of energy-absorption coefficients was also generated, which is used by the TLE calculations and is now also included in GATE/Geant4. The validation was carried out by comparing the TLE-simulated doses with experimental data in a synchrotron radiation computed tomography experiment. The TLE technique shows good agreement versus both experimental measurements and the results of a classical Monte Carlo simulation. Compared with the latter, it is possible to reach a pre-defined statistical uncertainty in about two to three orders of magnitude less time for complex geometries without loss of accuracy.

Design study of a Compton camera for prompt γ imaging during ion beam therapy

The goal of this study is to tune the design of the absorber detector of a Compton camera for prompt γ-ray imaging during ion beam therapy. The response of the Compton camera to a photon point source with a realistic energy spectrum (corresponding to the prompt γ-ray spectrum emitted during the carbon irradiation of a water phantom) is studied by means of Geant4 simulations. Our Compton camera consists of a stack of 2 mm thick silicon strip detectors as a scatter detector and of a scintillator plate as an absorber detector. Four scintillators are considered: LYSO, NaI, LaBr3 and BGO. LYSO and BGO appear as the most suitable materials, due to their high photo-electric cross-sections, which leads to a high percentage of fully absorbed photons. Depth-of-interaction measurements are shown to have limited influence on the spatial resolution of the camera. In our case, the thickness which gives the best compromise between a high percentage of photons that are fully absorbed and a low parallax error is about 4 cm for the LYSO detector and 4.5 cm for the BGO detector. The influence of the width of the absorber detector on the spatial resolution is not very pronounced as long as it is lower than 30 cm.

An efficient numerical tool for dose deposition prediction applied to synchrotron medical imaging and radiation therapy.

Longitudinal prompt-gamma ray profiles have been measured with a multi-slit multi-detector configuration at a 75 MeV/u 13C beam and with a PMMA target. Selections in time-of-flight and energy have been applied in order to discriminate prompt-gamma rays produced in the target from background events. The ion ranges which have been extracted from each individual detector module agree amongst each other and are consistent with theoretical expectations. In a separate dedicated experiment with 200 MeV/u 12C ions the fraction of inter-detector scattering has been determined to be on the 10%-level via a combination of experimental results and simulations. At the same experiment different collimator configurations have been tested and the shielding properties of tungsten and lead for prompt-gamma rays have been measured.