Compton Camera and Prompt Gamma Ray Timing: Two Methods for In Vivo Range Assessment in Proton Therapy
- 1Institute of Radiooncology, Helmholtz-Zentrum Dresden – Rossendorf, Dresden, Germany
- 2OncoRay – National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Dresden, Germany
- 3Institute of Radiation Physics, Helmholtz-Zentrum Dresden – Rossendorf, Dresden, Germany
- 4German Cancer Consortium (DKTK), German Cancer Research Center (DKFZ), Heidelberg, Germany
Proton beams are promising means for treating tumors. Such charged particles stop at a defined depth, where the ionization density is maximum. As the dose deposit beyond this distal edge is very low, proton therapy minimizes the damage to normal tissue compared to photon therapy. Nevertheless, inherent range uncertainties cast doubts on the irradiation of tumors close to organs at risk and lead to the application of conservative safety margins. This constrains significantly the potential benefits of protons over photons. In this context, several research groups are developing experimental tools for range verification based on the detection of prompt gammas, a nuclear by-product of the proton irradiation. At OncoRay and Helmholtz-Zentrum Dresden-Rossendorf, detector components have been characterized in realistic radiation environments as a step toward a clinical Compton camera. On the one hand, corresponding experimental methods and results obtained during the ENTERVISION training network are reviewed. On the other hand, a novel method based on timing spectroscopy has been proposed as an alternative to collimated imaging systems. The first tests of the timing method at a clinical proton accelerator are summarized, its applicability in a clinical environment for challenging the current safety margins is assessed, and the factors limiting its precision are discussed.
In the first decades of the 20th century, during the rise of particle accelerators, physicists studied the interaction of fast charged particles with matter. The energy loss of heavy ions (as opposed to light electrons) within a target medium was described by Bethe’s stopping power formula (1). The ionization, namely, the Coulomb collisions where the accelerated ions strip out electrons of the atoms of the target, is the predominant loss mechanism for non-relativistic ion beams (2).
The engineering race toward high-energy accelerators endowed heavy charged particles a penetration depth in tissue comparable to the body dimensions. This opened up the possibility of using protons for medical applications, as neutrons, electrons, gamma, or X-rays had been applied before in the field of radiotherapy, which emerged after Röntgen’s X-ray discovery in 1895 (3).
In 1946, Wilson predicted the physical, in particular dosimetric, properties of a proton beam (4) for a therapeutic scenario and founded the field of proton therapy. The straight beam trajectory, the finite particle range, as well as the increase of the ionization density close to the stopping point, aroused the interest of the medical community. In the context of cancer treatment, this ionizing radiation was expected to damage the cells of the target tumor and eventually cause their death, while sparing most efficiently surrounding normal tissue.
The first experimental treatments were performed during the 1950s at Berkeley Radiation Laboratory, USA, and in Uppsala, Sweden (2, 5). However, it was not until 1990 that the first hospital-based proton facility in Loma Linda University (USA) was created. Since then, the number of therapy centers has increased steadily, and carbon or other ions have been also introduced. Nowadays, more than 15,000 patients are treated per year in around 50 facilities worldwide (5).
Several distinguishing features of accelerated protons are listed below:
• The particle range: the protons stop at a defined penetration depth, which depends on their initial velocity.
• The Bragg peak: the ionization cross section increases for low proton velocities, so that the dose deposition density is maximum close to the particle range (2). Conversely, the dose at the entrance point is minimum.
• The distal penumbra: beyond the particle range, the dose deposition falls steeply to zero.
• The lateral penumbra: the beam trajectory is straight and the spread in the transversal dimension due to multiple Coulomb scattering is small but increases steadily with depth [(6), p. 15].
• Tissue composition: the dose deposition curve and proton range are strongly dependent on the stopping power of the traversed tissue (on its density and composition).
• Secondary products: neutrons, annihilation photons, and prompt gammas are produced throughout the proton track. These secondary by-products release a small dose compared to the incident protons (2), and they can exit the patient. Neutron emission is focused in forward direction (7).
In theory, proton therapy has several advantages over conventional photon therapy:
• The distal edge: the steep dose gradient at the distal edge is promising for sparing critical organs close to the irradiated tumor. In contrast, the slowly falling depth dose curve of photons impedes that strategy.
• The integral dose: thanks to the Bragg peak, the dose can be focused on the tumor volume, and the damage on normal tissue before and beyond is minimized. For photons, the dose in healthy tissue can only be distributed on a larger volume by irradiating from several directions, but its integral is higher than for protons.
The main drawbacks compared to photons are:
• The lack of large clinical trials and evidence about the superiority of proton beam therapy for the majority of tumor entities (10). Whether or not the lower integral dose translates to better clinical outcome (11), e.g., less normal-tissue toxicity, is yet to be proven.
• Range uncertainties (12) of the proton beam due to intrinsic factors, such as patient or organ motion, as well as restricted knowledge about the tissue composition, are prone to yield severe differences with respect to the planned dose, specially at the distal edge, whereas photon plans are much less sensitive in this regard (13). To circumvent the risk, safety margins are applied and robust treatment plans are designed, at the price of an increase of the dose in the normal tissue compared to the dose-optimum (but more risky) plan.
These disadvantages question the cost-effectiveness of ion beam therapy and fuel the controversy about their clinical superiority [(14), chapters 2.11–2.13] over photon therapy. There is an urgent need for techniques that tackle one of the major weaknesses of proton therapy: the intrinsic range uncertainties, which limit the ultimate precision with which ion beams can be safely delivered. The most common sources of range errors are:
• Stopping power ambiguity due to degeneracy of Hounsfield values depending on tissue composition (15).
• Patient alignment errors.
• Anatomy changes between or during treatment fractions, as cavity filling, change of weight, tissue swelling, or tumor shrinkage.
• Organ motion in the thoracic and abdominal region.
• Biological factors (16).
The proton range is strongly dependent on the composition of the traversed tissue. Photons are less dependent on these factors, and the absence of a sharp edge constrains the maximum dose deviations due to target shifts or path composition variations. The absence of tools in clinical routine for measuring in vivo and in real time, the actual distal fall-off edge, together with the high sensitivity of the proton range to tissue composition, force medical physicists to add safety margins and apply field patching techniques in order to obtain a robust treatment plan (17). Notwithstanding the theoretically superior dose profile of ions, broad safety margins (Figure 1) waste substantially the outstanding traits of ion over photon beam therapy.
Figure 1. Safety margins applied at different clinical proton therapy facilities: (3.5% + 3-mm) at Loma Linda University Medical Center (LLUMC) (3.5% + 2-mm) at Universitäts Protonen Therapie Dresden (UPTD) (3.5% + 1-mm) at Massachusetts General Hospital (MGH) and (2.5% + 1.5-mm) at University of Florida Health Proton Therapy Institute (UFH). Range bonus refers to the margin added to the prescribed range to ensure full tumor coverage even in the case of an undershoot. These centers may apply bigger margins in specific treatment scenarios (62).
It should be emphasized that most cancers are treated successfully with surgery, electron or photon therapy, chemotherapy, brachytherapy, whereas proton therapy covers just a residual percentage (18). Still, the improvement in the accelerator technology, delivery systems, and the trend toward personalized medicine make proton beams an attractive alternative for certain patient ages and types of tumors. It is estimated that ~10% of cancer patients, especially children, would benefit from proton therapy (reduction of late side effects) compared to conventional techniques (18). Hence, proton therapy is still in the headlines, the number of facilities is increasing from year to year, and questions concerning the improvement of the technique and quality assurance are of great interest.
In this context, several groups across the world aim at an experimental device that measures the particle range and even the dose profile, preferably in real time (13). Numerous techniques have been proposed in the last two decades and are reviewed in Ref. (13, 19, 20). This paper fits into this context and summarizes two different methods for monitoring the dose delivery of proton beams in real time based on Prompt Gamma Imaging (PGI).
Prompt gamma rays, a by-product emitted in nuclear reactions along the proton track, cf. Figure 2, cover a broad energy spectrum with several prominent characteristic lines, cf. Ref. (21). The high gamma ray energy ensures that they can be detected outside the patient without severe attenuation. The spatial emission distribution correlates to the dose deposition map of the incident protons (21, 22) and provides an indirect measurement of the particle range. Such correlation is dependent on prompt gamma ray energy and tissue composition (23–26) and stems from the maximum of the nuclear cross section at low (~10-MeV) proton energies (27).
Figure 2. Schematic of a typical prompt gamma ray production process 12C(p,p′)12C*. Left: nuclear collision of a proton (p) with a 12C nucleus of the target. Center: the proton is scattered (p′) and the nucleus is left in an excited state 12C*. Right: the relaxation to the ground state 12C is accompanied by a 4.4-MeV prompt gamma emission γp.
These gamma rays are prompt, i.e., they are emitted almost instantaneously after the collision, which is interesting for real-time range verification. The gamma ray production over 1-MeV is considerable, around 0.16 per proton (on average) at 160-MeV beam energy (28). The gamma ray emission rate depends on the beam current, duty cycle, and micro-time structure of the considered accelerator. Taking as an example the Cyclone® 230 (C230) isochronous cyclotron of IBA (Louvain-la-Neuve, Belgium) and a realistic treatment plan with pencil beam scanning, the peak beam current is ≈2-nA, there are about 109 gamma ray emissions per second, and 106-s−1 events are registered in a ø2″ × 2″ LaBr3 scintillator at 30-cm distance (29). This large gamma ray rate, as well as the inherent neutron background, poses a serious challenge on the detector and electronics design. Note that the so-called neutron background is mostly indirect, due to the detection of gamma rays following neutron interactions or neutron captures in surrounding materials, rather than from the interactions of neutrons in the detector itself.
In the field of nuclear medicine, commercial gamma cameras are used in clinical routine to obtain images of gamma ray distributions. Hence, one may think that the imaging of prompt gamma rays is not an issue, as the technology is already established. However, together with the detection rate and background, the high gamma ray energies and polychromatic energy spectrum prevent the direct use of the gamma camera as PGI device. In comparison, the gamma ray energies in nuclear medicine range between 80 and 511-keV. This significant difference is outlined in Figure 3: larger collimators and detectors are needed to absorb high-energy prompt gamma rays, normally after multiple interaction processes. For example, a 2-mm layer of lead has a 99% attenuation power for 140-keV photons, but only 9% for 4.4-MeV gammas; a 1-cm thick CsI crystal has a detection efficiency of 98% for the first and just 15% for the latter. Thick collimators reduce the system efficiency and deteriorate the image quality, whereas large detectors increase critically the system price and enlarge the footprint. Hence, alternative concepts are needed.
Figure 3. Mean free path λ of photons of energy Eγ in different representative materials: water (body), CsI (detector), and Pb (collimator). λ is calculated as the inverse of the total attenuation factor μ that is extracted from (63). Explicit points at typical energies are drawn: 140-keV for 99Tc (nuclear medicine), 511-keV for PET, and 4-MeV for PGI.
Dedicated PGI detector systems have been designed and tested in the last decade based on active or passive collimation. A pin-hole camera (30) is the pioneer approach to scan the prompt gamma emission distribution in a right angle to the beam track. Many research groups have performed experiments based on slit cameras at proton or carbon beams (31–37). The knife-edge-shaped camera has demonstrated the feasibility of millimeter range verification at clinical current intensities (38) in real time on a spot basis with realistic treatment plans and heterogeneous phantoms (39).
Among actively collimated systems, most efforts are concentrated on the Compton camera (40). It comprises multiple position sensitive gamma ray detectors, which are arranged in one scatterer and one absorber, or in several scatter planes. The prompt gamma rays reach the detectors, and the energy deposit as well as the point of interaction in each plane are measured, cf. Figure 4 for the two-plane camera. The Compton equation (41) relates the scattering angle θ to the initial (Eγ) and final (Eγ′) photon energies:
where Ls and La are the energies released in scatterer and absorber, respectively, and mec2 = 511-keV is the electron rest energy. In contrast to a slit camera, no collimation is needed in order to reconstruct the angle of incidence of the gamma ray, and two-dimensional (2D) or even three-dimensional (3D) images instead of one-dimensional (1D) profiles may be obtained. More single gamma rays and directions can be detected, but the condition of simultaneous interaction in the different camera stages limits the overall efficiency. Furthermore, the instrumentation requirements in terms of spatial, time, and energy resolution for the detectors of a Compton camera are especially high, and the reconstruction algorithm is complex and computationally intensive, as the incident direction cannot be recovered univocally for each event. Nowadays, a PGI Compton camera prototype demonstrating range verification in a clinical scenario is still a challenge several institutes aim at (42–47), and the only published experimental results at a proton beam are constrained to <2-MeV gammas (48, 49) or to beam currents far below the clinical case (50). Technical complexity, electronics expense, low coincident efficiency, high detector load, radiation background, and the elevated percentage of random coincidences are intrinsic hurdles that cast doubts on the applicability of this concept (19).
Figure 4. Incoherent scattering event in a two-plane Compton camera. The cone surface contains the possible incidence directions (any generatrix) of the initial photon (γ). It interacts with the scatterer plane and deposits an energy Ls. The scattered photon (γ′) releases the rest of the energy La in the absorber. The line connecting both interaction points (in orange) is the propagation direction of γ′. This defines the axis (directrix) of the aforementioned cone, with half-opening (scattering) angle θ and vertex at the scatter point.
In the recent years, one can identify a trend toward less complicated PGI systems, at least concerning hardware. These may have a faster translation into clinical practice due to their lower price (35, 37, 51). The Prompt Gamma Ray Timing (PGT) method (28) is one of these novel approaches, which relies on a single monolithic detector with excellent timing resolution and no collimation. As a consequence of the measurable transit time of ions through matter, the detection times of prompt gamma rays encode essential information about their spatial emission point. Figure 5 illustrates this physical effect: the deeper the proton interaction (prompt gamma emission) point, the larger the proton transit time and time of flight of the gammas to the detector. Applying the Continuous Slowing Down Approximation (CSDA), the transit time can be derived mathematically (28) if the initial beam energy E0 and target composition are known. First, the proton energy variation per unit length yields
where E is the kinetic energy, ρ (z) the mass density of the target at a depth z, and S(E) the stopping power, that depends on target material and thus indirectly on the depth z. The kinetic energy E of the proton at any depth z > z0 is
Figure 5. Illustrative sketch of the PGT method. The therapeutic proton p (the projectile) slows down as it penetrates the target and interacts with a nucleus, which emits a prompt gamma ray. The time between the entrance of the proton to the target (the start flag) and the arrival of the gamma ray to the detector (the stop flag at the yellow rectangle) encodes the proton transit time and gamma time of flight, which can be correlated to the depth of interaction of the proton (gamma emission point).
The relativistic velocity v of the particle is a function of the kinetic energy E:
where mp is the proton rest mass. Finally, the equation of the proton transit time yields
where dz′ has been exchanged with dE using equation (2).
The low cost and small footprint of PGT makes this concept very tempting. A major limitation is that the time spectra are not only blurred by the resolution of the detector but also by the time width of the accelerated bunches. This implies that the PGT method is not applicable at all clinical accelerators: only to those with a specific micro-time structure. For the widespread accelerator C230 (5), the micro bunch time spread can reach up to 2-ns for clinical beam energies (29). Here, range shifts can be identified based on distribution momenta. It is under exploration whether PGT is only applicable for pencil beams or also for passively scattered ones. In order to know if other types of clinical accelerators are compatible with the PGT approach, the specific micro pulse structure has to be measured.
Rather than an in-depth review of the literature, this manuscript provides a summary of two separate topics developed within the collaborative framework of ENTERVISION (52): the characterization of detector components for the absorber plane of a clinical Compton camera and the first test of the PGT method with heterogeneous phantoms at a clinical proton center, corresponding to the publications (53) and (29), respectively.
2. Comparison of BGO and LSO Scintillation Detectors
Bi4Ge3O12 (BGO) and Lu2SiO5:Ce (LSO) scintillators are straightforward candidates to absorber detectors of a Compton camera aiming at PGI. These are used traditionally in Positron Emission Tomography (PET) scanners. Despite its higher price and 30% lower photoabsorption efficiency, LSO has gained importance due to its higher light yield and fast decay time. It is questionable if this conclusion can be transcribed to the PGI scenario. In order to assess the choice between the two options, benchmark experiments are conducted at different accelerators for comparing BGO and LSO detectors in terms of energy, spatial, and time resolution. Other factors, such as intrinsic radiation, absorption efficiency, and cost-effectiveness ratio, are also discussed.
The basic detection unit in commercial PET scanners is the block detector. It consists of a square matrix of segmented or pixellated scintillating crystals coupled to four light-sharing Photomultiplier Tubes (PMT), as depicted in Figure 6. The pixel where the gamma ray interacts can be calculated based on the ratio of light collected by each PMT. The block detectors used in this comparative study and their properties are listed in Table 1. They are named as LSO2 and BGO1 when referring to the concrete detector results, in contrast to LSO and BGO when speaking about general features of the scintillation materials.
Figure 6. Sketch of the LSO2 (left) and BGO1 (right) block detectors with the PMT numbering convention and crystal coordinate system, namely, the XBlock and YBlock axis (relative position between 0 and 1). Crystals are depicted in orange, PMTs in blue, and the light guide in yellow. Reproduced with permission from Ref. (53).
Table 1. Comparison of the properties of the different block detectors from Siemens Medical Solutions USA, Inc. Molecular Imaging Division, whose sketch is depicted in Figure 6.
In Figure 7, the relative energy resolution RE as Full Width at Half Maximum (FWHM) of the LSO2 and BGO1 detectors is compared. The energy resolution of LSO2 is better across the whole energy range. At 511-keV, RE,LSO2 ≈ 11% and RE,BGO1 ≈ 18%. Nonetheless, the differences are less pronounced for high-energy photopeaks (53, 54). At 4.4-MeV, RE,LSO2 ≈ 7% versus RE,BGO1 ≈ 8%. In other words, LSO2 excels at the PET scenario (below 1-MeV), whereas for the PGI scenario (above 2-MeV), the difference in performance is less significant, and the higher price of LSO does not imply a much better detector quality.
Figure 7. Relative energy resolution RE = ΣE/E of the LSO2 (black line) and BGO1 (red line) block detectors as a function of energy deposit E. ΣE refers to the FWHM of the photopeak. The empirical fit to the experimental points is for LSO2 and for BGO1. Reproduced with permission from Ref. (53).
The reason for the comparable energy resolution at high energies is the following. The relative energy resolution RE depends on two independent contributions: the statistical and the intrinsic one. The first one depends on the light yield and is proportional to the inverse square root of number of (collected) scintillation photons. The latter is due to non-linearity effects (55) and is dependent on the crystal structure. At low photon energies, i.e., the range of usual radioactive sources or in case of PET, the statistical contribution dominates over the intrinsic one. As LSO has a four times higher light yield than BGO, its energy resolution is significantly better. At high photon energies, i.e., the PGI energy range, the number of scintillation photons is larger, so that the statistical contribution is smaller and the intrinsic contribution starts to dominate. This intrinsic factor is comparable for BGO and LSO (54, 56, 57), which explains their similar performance concerning energy resolution at the PGI energy range.
With respect to the spatial resolution, one can conclude that the discrimination power between pixels of the flood map increases with the energy range for both block detectors, cf. Figure 8. This effect is also due to the lower statistical relative uncertainty of events with high-energy deposit. The spots in the flood map of BGO1 are broader than those of LSO2 in any case, but become very sharp in the PGI range. This points out that one could segment the BGO1 block detector in, e.g., 13 × 13 instead of 8 × 8 and achieve the spatial resolution of LSO2 without jeopardizing the pixel discrimination in the flood map of the PGI range.
Figure 8. Block detector flood map (top figure: LSO2 and bottom figure: BGO1) for different energy ranges at the ELBE accelerator (distortion correction applied for LSO2). Non-uniformities are due to the fact that the focused bremsstrahlung beam spot is smaller than the detector size as well as the different extension (number of bins) of each pixel spot in the map. Reproduced with permission from Ref. (53).
Regarding the time resolution, cf. Figure 9, the LSO2 detector beats BGO1 over the whole energy range thanks to the larger light yield and shorter decay time. A good time resolution is mandatory for a PGI Compton camera in order to suppress delayed radiation background (58). In order to analyze if the timing resolution of BGO1 is sufficient for this goal, we calculate the figure of merit FoMBSR (59):
where Σt,det is the detector time resolution as FWHM, Σt,bunch the bunch time spread (FWHM), and the bunch period (the inverse of the radio frequency). This ratio measures the amount of background that can be suppressed thanks to timing measurements in a pulsed accelerator. In the case of the C230 machine, where Tbunch = 9.4-ns and Σt,bunch ≈ 2-ns (for 100-MeV protons) (29), the background suppression ratios for 4-MeV prompt gammas are FoMBSR,LSO2 = 84% and FoMBSR,BGO1 = 64%. LSO2 has undoubtedly better performance, but BGO1 could be also acceptable.
Figure 9. Time resolution Σt,det (FWHM) of the LSO2 (black line) and BGO1 (red line) block detectors measured at the ELBE accelerator as a function of the released energy. The experimental points are approximately reproduced by the curves for LSO2 and for BGO1. Reproduced with permission from Ref. (53).
Other material features, such as the decay time and the intrinsic radioactivity, are worth to discuss. The decay time of BGO (7.5 times longer than LSO) implies a limit of around 300-kcps detector load. Taking into account the high rates expected in the PGI scenario, about 1-Mcps, one might be forced to reduce the area of the BGO block detectors or increase the distance to the beam axis. These rates are also quite challenging for the electronics and data processing. On the other hand, it is well known that LSO has a high intrinsic radioactivity below 1-MeV due to 176Lu, namely, through β− decay and subsequent gamma ray cascade. The simultaneous detection of the electron (stopped in the LSO crystal) and the gamma ray (in the scatterer plane) could produce a significant fraction of false coincidences in a Compton camera.
3. Prompt Gamma Ray Timing with Heterogeneous Targets
The PGT method is a promising and novel method for range verification proposed by Golnik et al. (28) based on experiments with homogeneous phantoms at a research accelerator. To further explore its potential and the limitations that may appear when translating the concept to a realistic irradiation scenario, specific experiments at a clinical proton center with heterogeneous targets are conducted. The concrete goals are to test the robustness of the PGT method, its precision, and limitations, as well as the capability of detecting range shifts due to heterogeneities. Furthermore, the next steps toward a clinical PGT prototype are identified.
The experiment is carried out at the Westdeutsches Protonentherapiezentrum Essen (WPE), Germany. This clinical proton center comprises a C230 cyclotron, with a radio frequency close to 106-MHz. A horizontal pencil beam (no scanning) in the gantry treatment room irradiates a cylindrical target, see Figure 10 (bottom). The inner shell contains slices of custom thickness and composition, so that different heterogeneous targets can be configured. Available materials are polymethyl methacrylate (PMMA), air (hollow slice), and cortical bone. The detectors listed in Table 2 are set up at an angle α and distance d, as described in Figure 10 (top). Whenever a target is thick enough to completely stop the impinging protons, we label it as full (for a given proton energy).
Figure 10. Top: sketch of the experimental setup at the WPE proton therapy facility and of the longitudinal cross section of the target. The accelerated protons bunched with 106 MHz collide with the PMMA cylindrical target. It contains a cavity or bone insert of thickness h located at a distance f from the beam entrance point. Bottom: photograph of the experimental setup with three detectors at different ring angles α. The linear stage on the center of the ring holds the two hollow joined half cylinders, in which PMMA, cavity, or bone slices can be inserted. The beam incidence is horizontal from the left, where the snout of the nozzle is seen. Reproduced with permission from Ref. (29).
The bunch time spread is an important limiting factor for the PGT method. In Figure 11 (left), the bunch width is characterized as a function of the proton energy, ranging from 350-ps at 230-MeV to 2-ns at 100-MeV. The time spread can be reduced up to a factor of two by adjusting the momentum spread limiting slits, cf. Figure 11 (right), the main component of the energy selection system of the C230 cyclotron (60).
Figure 11. Left: PGT spectra of the B1 detector with a thin PMMA target for three different proton energies and the usual slit closing (25-mm) at the WPE gantry. Right: bunch time spread Σt,bunch at 100-MeV proton energy and a full PMMA target for the B3 detector. Reproduced with permission from Ref. (29).
For the acquisition of PGT spectra, the detection time of the gamma rays with respect to the arrival of the protons to the target has to be measured. As usual in research accelerators, the radio frequency can be used as reference time for the bunch arrival (except for an offset). However, Figure 12 shows that this offset is not constant on a large time scale (29). These phase drifts of the proton bunch with respect to the RF signal may be caused due to temperature changes or main coil current instabilities, among other factors.
Figure 12. Left: PGT spectra of the B3 detector, a homogeneous PMMA target and a proton energy of 230-MeV. Independent redundant measurements of about 5-min with a separation of around 1-h are overlaid. Right: PGT spectra of the B3 detector, a full homogeneous PMMA target, and a proton energy of 160-MeV for different values of the main coil current of the cyclotron (detailed in the legend). Reproduced with permission from Ref. (29).
With regard to homogeneous targets, the stacked target experiment is accomplished as follows: homogeneous PMMA targets of various thicknesses are irradiated with 230-MeV protons. The increase of the target thickness correlates to an increase of the area of the PGT distribution and a shift to the right in its mean value, due to the enlarged region of prompt gamma emission, as observed in Figure 13 (left). As the bunch time spread is significantly lower than the proton transit time, one can resolve the prompt gamma emission density as a function of the (timewise) depth with much less blurring than for 160 or 100-MeV proton energies. In Figure 13 (right), we calculate the PGT distributions according to the analytical simple Box (simBox) model (28). The shape is qualitatively similar to the experimental spectra, but the model is too simple to reproduce, e.g., the fall-off edge corresponding to the Bragg peak or the background radiation.
Figure 13. Left: experimental PGT spectra of the B1 detector for 230-MeV protons and different PMMA target thicknesses. Right: spectra calculated with the analytical simBox model (28). Both: the two vertical dashed lines refer to the expected front face and proton range positions. Reproduced with permission from Ref. (29).
Concerning heterogeneous targets, air cavities of different thicknesses are placed at different depths inside a full PMMA target. The experimental results are shown in Figure 14. The deficit in the gamma ray production inside the air cavity can be identified as a dip in the PGT spectra at a time position and with a magnitude correlated to its location and thickness, respectively. The falling edge of the spectrum shifts steadily to the right according to the cavity thickness (beam overshoot).
Figure 14. Experimental PGT spectra of the L0 detector for 230-MeV protons and different air cavities inside the full (400-mm) PMMA target at f and h (front face position and thickness as described in Figure 10). Reproduced with permission from Ref. (29).
An analogous experiment with a bone insert (20-mm thick) at different positions is carried out (29). The resulting PGT spectra are depicted in Figure 15 (left). An increase of the gamma ray production due to the higher density of bone is visible in the PGT spectrum at a time correlated to the insert position. A shift to the left in the falling edge of the distribution can be identified (undershoot) with respect to the homogeneous case. In Figure 15 (right), the PGT spectra are converted to depth profiles by using the transit time equation, cf. equation (5), and applying detector solid angle and gamma time-of-flight corrections.
Figure 15. Left: experimental PGT spectra of the L0 detector for 230-MeV protons and a bone insert inside the full (400-mm) PMMA target at f and h (front face position and thickness as described in Figure 10). Right: time to distance axis conversion from the left PGT spectra after application of stopping power, sensitivity, and gamma time-of-flight corrections as well as background subtraction. z refers to the depth with respect to the target’s front face (beam entrance point). Vertical dashed lines mark the centroid of the bump according to the simBox model. Reproduced with permission from Ref. (29).
In Figure 16, for 230-MeV protons, the effect of a beam overshoot (air cavity) and undershoot (bone insert) on PGT spectra is compared with respect to a reference measurement (homogeneous PMMA target). Moreover, the detectability of the range shift based on the location of the trailing edge is analyzed as a function of the number of irradiated protons. For clinical treatment plans, the strongest spot in pencil beam scanning, which is usually at the distal edge, yields close to 108 protons. This hints that a single detector is able to recover 5-mm range errors of the distal spot based on the PGT method with realistic beam currents at the C230 accelerator.
Figure 16. Top row and bottom left plot: PGT spectra of the L0 detector for 230-MeV protons. The homogeneous case (red curve) corresponds to a full PMMA target (400-mm). A heterogeneous slice is placed inside the full PMMA target at f = 169-mm and h = 5-mm (air cavity – blue curve) or f = 169-mm and h = 20-mm (bone insert – black curve), where f and h are the front face position and thickness, as described in Figure 10. The legend header contains the number of protons associated with each spectrum. Bottom right: shift of the falling edge with respect to the homogeneous case depending on the number of protons. The dashed lines depict the expected shift of the trailing edge according to the simBox model. Reproduced with permission from Ref. (29).
Prompt gamma rays, produced in nuclear reactions of accelerated protons with tissue, are valid signatures for retrieving the range of therapeutic protons. Several imaging systems are under development in the scientific community. Among others, the Compton camera and the Prompt Gamma Ray Timing (PGT) method have been studied intensively during the last years in collaboration with the ENTERVISION project (52).
A Compton camera requires position-sensitive detectors with high resolution and efficiency. The characterization of different candidate detectors in the PGI energy range is mandatory for assessing the material choice based on measurements that complement previous simulations and textbook knowledge. BGO and LSO block detectors from commercial PET scanners are compared in terms of energy, spatial, and time resolution, as well as price, absorption efficiency, and intrinsic background. As expected, the overall performance of LSO is better, but BGO closes the gap in the PGI range. The reason is that the high gamma ray energies (compared to the PET scenario) and thus number of scintillation photons balances the lower light yield. In addition, BGO has a higher photoabsorption efficiency, no intrinsic radioactivity, and low cost. Hence, BGO is a competitive alternative for the absorber of a Compton camera, thanks to the superior cost-effectiveness ratio in the PGI field (53).
PGT is an innovative method for range assessment based on a low footprint detector setup at minimum expense. First tests at a clinical accelerator and with heterogeneous phantoms reveal the capability of measuring 5-mm range shifts (due to heterogeneities) for beam spots with 108 protons (29). The bunch time spread is a crucial factor that affects the resolution of the PGT method. It depends on the delivered proton energy and the settings of the energy selection system. Furthermore, bunch phase drifts are found throughout the experiment, which pose a challenge on the robustness of the PGT method on a large time scale. Hence, it is advisable to introduce a proton bunch monitor (60, 61) that measures the bunch time structure as well as the potential phase drifts. A larger detector load and acquisition throughput are mandatory to improve the number of gamma rays detected per proton, so that statistically significant conclusions on range errors can be drawn for more spots of the treatment plan. Quantification of the range shifts based on more sophisticated models or simulations are also necessary. Experiments with upgraded detectors and electronics, realistic treatment plans, and anthropomorphic phantoms are ongoing. These are the next steps toward clinical translation and development of a first robust prototype.
All authors were involved in the conceptual design, management, support, and accomplishment of the experiments reviewed in this paper as well as in the critical revision of the data analysis and obtained results.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors like to thank S. Akhmadaliev3, D. Bemmerer3, M. Berthel2, A. K. Biegun, J. v. Borany3, P. Dendooven, A. Dreyer2, A. Hartmann3, K. Heidel3, S. Helmbrecht3, G. Janssens, A. Junghans3, U. Just2, M. Kempe3, A. Laso-García3, D. Prieels, M. Priegnitz3, H. Rohling3, K. Schmidt3, S. Schöne3, R. Schwengner3, J. Smeets, M. Sobiella3, F. Vander Stappen, A. Wagner3, L. Wagner3, and D. Weinberger3 for the excellent assistance as well as the crews of the different accelerators for stable operations. (Superscripts refer to affiliation numbers).
The summarized work was supported by the German Federal Ministry of Education and Research (BMBF-03Z1NN12) and the European Commission (FP7 Grant Agreement No. 241851 and No. 264552, and IA-ENSAR FP7 Contract No. RII3-CT-2010-262010).
5. Particle Therapy Co-Operative Group. Patient Statistics. Villigen: Particle Therapy Co-Operative Group (2014). Available from: http://www.ptcog.ch/archive/patient_statistics/Patientstatistics-updateMar2014.pdf
6. Fiedler F. Anwendung des in-Beam PET Therapiemonitorings auf Präisionsbestrahlungen mit Helium-Ionen. Ph.D. thesis, Technische Universität Dresden, Dresden (2008). Available from: http://nbn-resolving.de/urn:nbn:de:bsz:14-ds-1206441845447-47490
7. Robert C, Dedes G, Battistoni G, Böhlen TT, Buvat I, Cerutti F, et al. Distributions of secondary particles in proton and carbon-ion therapy: a comparison between GATE/Geant4 and FLUKA Monte Carlo codes. Phys Med Biol (2013) 58(9):2879. doi:10.1088/0031-9155/58/9/2879
9. Peeters A, Grutters JP, Pijls-Johannesma M, Reimoser S, Ruysscher DD, Severens JL, et al. How costly is particle therapy? Cost analysis of external beam radiotherapy with carbon-ions, protons and photons. Radiother Oncol (2010) 95(1):45–53. doi:10.1016/j.radonc.2009.12.002
10. Freeman T. Proton Therapy: Randomized Trials or Model-Based Selection? Bristol: Medical Physics Web (2015). Available from: http://medicalphysicsweb.org/cws/article/opinion/61592
11. Freeman T. Will Protons Gradually Replace Photons? Bristol: Medical Physics Web (2012). Available from: http://medicalphysicsweb.org/cws/article/research/50584
14. Orton CG, Hendee WR. Controversies in Medical Physics: A Compendium of Point/Counterpoint Debates. College Park, MD: American Association of Physicists in Medicine (2008). p. 110–23. Available from: https://www.aapm.org/pubs/books/PointCounterpointCompendium.pdf
16. Grün R, Friedrich T, Krämer M, Zink K, Durante M, Engenhart-Cabillic R, et al. Physical and biological factors determining the effective proton range. Med Phys (2013) 40(11):111716. doi:10.1118/1.4824321
19. Smeets J. Prompt Gamma Imaging with a Slit Camera for Real Time Range Control in Particle Therapy. Ph.D. thesis, Université Claude Bernard Lyon 1, Lyon (2012). Available from: http://theses.ulb.ac.be/ETD-db/collection/available/ULBetd-10182012-123741/unrestricted/Thesis.pdf
23. Moteabbed M, España S, Paganetti H. Monte Carlo patient study on the comparison of prompt gamma and PET imaging for range verification in proton therapy. Phys Med Biol (2011) 56(4):1063. doi:10.1088/0031-9155/56/4/012
24. Polf JC, Panthi R, Mackin DS, McCleskey M, Saastamoinen A, Roeder BT, et al. Measurement of characteristic prompt gamma rays emitted from oxygen and carbon in tissue-equivalent samples during proton beam irradiation. Phys Med Biol (2013) 58(17):5821. doi:10.1088/0031-9155/58/17/5821
25. Gueth P, Dauvergne D, Freud N, Létang JM, Ray C, Testa E, et al. Machine learning-based patient specific prompt-gamma dose monitoring in proton therapy. Phys Med Biol (2013) 58(13):4563. doi:10.1088/0031-9155/58/13/4563
26. Janssen F, Landry G, Lopes PC, Dedes G, Smeets J, Schaart DR, et al. Factors influencing the accuracy of beam range estimation in proton therapy using prompt gamma emission. Phys Med Biol (2014) 59(15):4427. doi:10.1088/0031-9155/59/15/4427
28. Golnik C, Hueso-González F, Müller A, Dendooven P, Enghardt W, Fiedler F, et al. Range assessment in particle therapy based on prompt γ-ray timing measurements. Phys Med Biol (2014) 59(18):5399. doi:10.1088/0031-9155/59/18/5399
29. Hueso-González F, Enghardt W, Fiedler F, Golnik C, Janssens G, Petzoldt J, et al. First test of the prompt gamma ray timing method with heterogeneous targets at a clinical proton therapy facility. Phys Med Biol (2015) 60(16):6247. doi:10.1088/0031-9155/60/16/6247
30. Seo KS, Kim CH, Kim JW. Comparison of titanium hydride (TiH2) and paraffin as neutron moderator material in a prompt gamma scanning system. J Korean Phys Soc (2006) 48(4):855. doi:10.3938/jkps.48.855. Available from: http://www.kps.or.kr/jkps/abstract_view.asp?articleuid=CB6536D4-BA6B-4C17-A954-7B0A6A2566B7&globalmenu=3&localmenu=10
31. Testa E, Bajard M, Chevallier M, Dauvergne D, Le Foulher F, Freud N, et al. Monitoring the Bragg peak location of 73MeV/u carbon ions by means of prompt γ-ray measurements. Appl Phys Lett (2008) 93(9):093506. doi:10.1063/1.2975841
33. Pinto M, Dauvergne D, Freud N, Krimmer J, Letang JM, Ray C, et al. Design optimisation of a TOF-based collimated camera prototype for online hadrontherapy monitoring. Phys Med Biol (2014) 59(24):7653. doi:10.1088/0031-9155/59/24/7653
34. Cambraia Lopes P, Pinto M, Simoes H, Biegun A, Dendooven P, Oxley D, et al. Optimization of collimator designs for real-time proton range verification by measuring prompt gamma rays. Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC). Anaheim, CA: IEEE (2012). p. 3864–70. doi:10.1109/NSSMIC.2012.6551886
37. Testa M, Min CH, Verburg JM, Schümann J, Lu HM, Paganetti H. Range verification of passively scattered proton beams based on prompt gamma time patterns. Phys Med Biol (2014) 59(15):4181. doi:10.1088/0031-9155/59/15/4181
38. Perali I, Celani A, Bombelli L, Fiorini C, Camera F, Clementel E, et al. Prompt gamma imaging of proton pencil beams at clinical dose rate. Phys Med Biol (2014) 59(19):5849. doi:10.1088/0031-9155/59/19/5849
39. Priegnitz M, Helmbrecht S, Janssens G, Perali I, Smeets J, Stappen FV, et al. Measurement of prompt gamma profiles in inhomogeneous targets with a knife-edge slit camera during proton irradiation. Phys Med Biol (2015) 60(12):4849. doi:10.1088/0031-9155/60/12/4849
42. Seo H, Park JH, Ushakov A, Kim CH, Kim JK, Lee JH, et al. Experimental performance of double-scattering Compton camera with anthropomorphic phantom. J Instrum (2011) 6(01):C01024. doi:10.1088/1748-0221/6/01/C01024
43. Krimmer J, Ley J-L, Abellan C, Cachemiche J-P, Caponetto L, Chen X, et al. Development of a Compton camera for medical applications based on silicon strip and scintillation detectors. Nucl Instrum Methods Phys Res A (2014) 787:98–101. doi:10.1016/j.nima.2014.11.042
44. Llosá G, Cabello J, Callier S, Gillam J, Lacasta C, Rafecas M, et al. First Compton telescope prototype based on continuous LaBr3-SiPM detectors. Nucl Instrum Methods Phys Res A (2013) 718(0):130–3. doi:10.1016/j.nima.2012.08.074
45. Thirolf P, Lang C, Aldawood S, Van der Kolff H, Maier L, Schaart D, et al. Development of a Compton camera for online range monitoring of laser-accelerated proton beams via prompt-gamma detection. EPJ Web of Conferences. (Vol. 66), Florence: EDP Sciences (2014). 11036 p. doi:10.1051/epjconf/20146611036
46. Kormoll T, Fiedler F, Schöne S, Wüstemann J, Zuber K, Enghardt W. A Compton imager for in-vivo dosimetry of proton beams – a design study. Nucl Instrum Methods Phys Res A (2011) 626-627:114–9. doi:10.1016/j.nima.2010.10.031
47. McCleskey M, Kaye W, Mackin D, Beddar S, He Z, Polf J. Evaluation of a multistage CdZnTe Compton camera for prompt γ imaging for proton therapy. Nucl Instrum Methods Phys Res A (2015) 785(0):163–9. doi:10.1016/j.nima.2015.02.030
48. Kabuki S, Ueno K, Kurosawa S, Iwaki S, Kubo H, Miuchi K, et al. Study on the use of electron-tracking Compton gamma-ray camera to monitor the therapeutic proton dose distribution in real time. Nuclear Science Symposium Conference Record (NSS/MIC). (Vol. J05-13), Orlando, FL: IEEE (2009). p. 2437–40. doi:10.1109/NSSMIC.2009.5402130
49. Kurosawa S, Kubo H, Ueno K, Kabuki S, Iwaki S, Takahashi M, et al. Prompt gamma detection for range verification in proton therapy. Curr Appl Phys (2012) 12(2):364–8. doi:10.1016/j.cap.2011.07.027
50. Polf JC, Avery S, Mackin DS, Beddar S. Imaging of prompt gamma rays emitted during delivery of clinical proton beams with a Compton camera: feasibility studies for range verification. Phys Med Biol (2015) 60(18):7085. doi:10.1088/0031-9155/60/18/7085
51. Assmann W, Kellnberger S, Reinhardt S, Lehrack S, Edlich A, Thirolf PG, et al. Ionoacoustic characterization of the proton Bragg peak with submillimeter accuracy. Med Phys (2015) 42(2):567–74. doi:10.1118/1.4905047
52. ENTERVISION. Research Training in 3D Digital Imaging for Cancer Radiation Therapy. Geneva: CERN – European Laboratory for Particle Physics (2015). Available from: http://entervision.web.cern.ch/entervision/
53. Hueso-González F, Biegun AK, Dendooven P, Enghardt W, Fiedler F, Golnik C, et al. Comparison of LSO and BGO block detectors for prompt gamma imaging in ion beam therapy. J Instrum (2015) 10(9):09015. doi:10.1088/1748-0221/10/09/P09015
54. Roemer K, Pausch G, Bemmerer D, Berthel M, Dreyer A, Golnik C, et al. Characterization of scintillator crystals for usage as prompt gamma monitors in particle therapy. J Instrum (2015) 10(10):10033. doi:10.1088/1748-0221/10/10/P10033
55. Dorenbos P, de Haas J, van Eijk C. Non-proportionality in the scintillation response and the energy resolution obtainable with scintillation crystals. IEEE Trans Nucl Sci (1995) 42(6):2190–202. doi:10.1109/23.489415
56. Balcerzyk M, Moszynski M, Kapusta M, Wolski D, Pawelke J, Melcher C. YSO, LSO, GSO and LGSO. a study of energy resolution and nonproportionality. IEEE Trans Nucl Sci (2000) 47(4):1319–23. doi:10.1109/23.872971
57. Moszynski M, Balcerzyk M, Czarnacki W, Kapusta M, Klamra W, Syntfeld A, et al. Intrinsic energy resolution and light yield nonproportionality of BGO. IEEE Trans Nucl Sci (2004) 51(3):1074–9. doi:10.1109/TNS.2004.829491
58. Biegun AK, Seravalli E, Lopes PC, Rinaldi I, Pinto M, Oxley DC, et al. Time-of-flight neutron rejection to improve prompt gamma imaging for proton range verification: a simulation study. Phys Med Biol (2012) 57(20):6429. doi:10.1088/0031-9155/57/20/6429
59. Hueso-González F, Golnik C, Berthel M, Dreyer A, Enghardt W, Fiedler F, et al. Test of Compton camera components for prompt gamma imaging at the ELBE bremsstrahlung beam. J Instrum (2014) 9(05):05002. doi:10.1088/1748-0221/9/05/P05002
60. Petzoldt J, Roemer KE, Enghardt W, Fiedler F, Golnik C, Hueso-González F, et al. Characterization of the microbunch time structure of proton pencil beams at a clinical treatment facility. Phys Med Biol (2016) 61(6):2432. doi:10.1088/0031-9155/61/6/2432
61. Werner T, Dreyer A, Golnik C, Hueso-Gonzalez F, Enghardt W, Straessner A, et al. Studies of a proton phase beam monitor for range verification in proton therapy. ANIMMA Conference Record 2015. Lisbon (2015).
63. Hubbell JH, Seltzer SM. Tables of X-Ray Mass Attenuation Coefficients and Mass Energy-Absorption Coefficients from 1 keV to 20 MeV for Elements Z = 1 to 92 and 48 Additional Substances of Dosimetric Interest. Gaithersburg, MD: NIST (2004). Available from: http://www.nist.gov/pml/data/xraycoef/
Keywords: proton therapy, range verification, in vivo dosimetry, Compton imaging, block detector, scintillation, prompt gamma ray timing
Citation: Hueso-González F, Fiedler F, Golnik C, Kormoll T, Pausch G, Petzoldt J, Römer KE and Enghardt W (2016) Compton Camera and Prompt Gamma Ray Timing: Two Methods for In Vivo Range Assessment in Proton Therapy. Front. Oncol. 6:80. doi: 10.3389/fonc.2016.00080
Received: 29 September 2015; Accepted: 21 March 2016;
Published: 12 April 2016
Edited by:Marco Durante, GSI Helmholtz Centre for Heavy Ion Research, Germany
Reviewed by:Nitin Ohri, Albert Einstein College of Medicine, USA
Paulo A. V. Crespo, University of Coimbra, Portugal
Copyright: © 2016 Hueso-González, Fiedler, Golnik, Kormoll, Pausch, Petzoldt, Römer and Enghardt. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fernando Hueso-González, firstname.lastname@example.org