Spectroscopic study of prompt-gamma emission for range verification in proton therapy

First experimental demonstration of real-time neutron capture discrimination in helium and carbon ion therapy

This work provides the first experimental proof of an increased neutron capture photon signal following the introduction of boron to a PMMA phantom during helium and carbon ion therapies in Neutron Capture Enhanced Particle Therapy (NCEPT). NCEPT leverages [Formula: see text]B neutron capture, leading to the emission of detectable 478 keV photons. Experiments were performed at the Heavy Ion Medical Accelerator in Chiba, Japan, with two Poly(methyl methacrylate) (PMMA) targets, one bearing a boron insert. The BeNEdiCTE gamma-ray detector measured an increase in the 478 keV signal of 45 ± 7% and 26 ± 2% for carbon and helium ion irradiation, respectively. Our Geant4 Monte Carlo simulation model, developed to investigate photon origins, found less than 30% of detected photons originated from the insert, while boron in the detector's circuit boards contributed over 65%. Further, the model investigated detector sensitivity, establishing its capability to record a 10% increase in 478 keV photon detection at a target [Formula: see text]B concentration of 500 ppm using spectral windowing alone, and 25% when combined with temporal windowing. The linear response extended to concentrations up to 20,000 ppm. The increase in the signal in all evaluated cases confirm the potential of the proposed detector design for neutron capture quantification in NCEPT.

Near-field coded-mask technique and its potential for proton therapy monitoring

Objective.Prompt-gamma imaging encompasses several approaches to the online monitoring of the beam range or deposited dose distribution in proton therapy. We test one of the imaging techniques - a coded mask approach - both experimentally and via simulations.Approach.Two imaging setups have been investigated experimentally. Each of them comprised a structured tungsten collimator in the form of a modified uniformly redundant array mask and a LYSO:Ce scintillation detector of fine granularity. The setups differed in detector dimensions and operation mode (1D or 2D imaging). A series of measurements with radioactive sources have been conducted, testing the performance of the setups for near-field gamma imaging. Additionally, Monte Carlo simulations of a larger setup of the same type were conducted, investigating its performance with a realistic gamma source distribution occurring during proton therapy.Main results.The images of point-like sources reconstructed from two small-scale prototypes' data using the maximum-likelihood expectation maximisation algorithm constitute the experimental proof of principle for the near-field coded-mask imaging modality, both in the 1D and the 2D mode. Their precision allowed us to calibrate out certain systematic offsets appearing due to the limited alignment accuracy of setup elements. The simulation of the full-scale setup yielded a mean distal falloff retrieval precision of 0.72 mm in the studies for beam energy range 89.5-107.9 MeV and with 1 × 108protons (a typical number for distal spots). The implemented algorithm of image reconstruction is relatively fast-a typical procedure needs several seconds.Significance.Coded-mask imaging appears a valid option for proton therapy monitoring. The results of simulations let us conclude that the proposed full-scale setup is competitive with the knife-edge-shaped and the multi-parallel slit cameras investigated by other groups.

Influence of the background in Compton camera images for proton therapy treatment monitoring

Objective. Background events are one of the most relevant contributions to image degradation in Compton camera imaging for hadron therapy treatment monitoring. A study of the background and its contribution to image degradation is important to define future strategies to reduce the background in the system.Approach. In this simulation study, the percentage of different kinds of events and their contribution to the reconstructed image in a two-layer Compton camera have been evaluated. To this end, GATE v8.2 simulations of a proton beam impinging on a PMMA phantom have been carried out, for different proton beam energies and at different beam intensities.Main results. For a simulated Compton camera made of Lanthanum (III) Bromide monolithic crystals, coincidences caused by neutrons arriving from the phantom are the most common type of background produced by secondary radiations in the Compton camera, causing between 13% and 33% of the detected coincidences, depending on the beam energy. Results also show that random coincidences are a significant cause of image degradation at high beam intensities, and their influence in the reconstructed images is studied for values of the time coincidence windows from 500 ps to 100 ns.Significance. Results indicate the timing capabilities required to retrieve the fall-off position with good precision. Still, the noise observed in the image when no randoms are considered make us consider further background rejection methods.

Prompt-gamma fall-off estimation with C-ion irradiation at clinical energies, using a knife-edge slit camera: A Monte Carlo study

PURPOSE

In-vivo range verification has been a hot topic in particle therapy since two decades. Many efforts have been done for proton therapy, while fewer studies were conducted considering a beam of carbon ions. In the present work, a simulation study was performed to show whether it is possible to measure the prompt-gamma fall-off inside the high neutron background typical of carbon-ion irradiation, using a knife-edge slit camera. In addition to this, we wanted to estimate the uncertainty in retrieving the particle range in the case of a pencil beam of C-ions at clinically relevant energy of 150 MeVu.

METHODS

For these purposes, the Monte Carlo code FLUKA was adopted for simulations and three different analytical methods were implemented to get the accuracy in the range retrieval of the simulated set-up.

RESULTS

The analysis of simulation data has brought to the promising and desired precision of about 4 mm in the determination of the dose profile fall-off in case of a spill irradiation, for which all the three cited methods were coherent in their predictions.

CONCLUSIONS

The Prompt Gamma Imaging technique should be further studied as a tool to reduce range uncertainties affecting carbon ion radiation therapy.

Machine learning-based event recognition in SiFi Compton camera imaging for proton therapy monitoring

Objective. Online monitoring of dose distribution in proton therapy is currently being investigated with the detection of prompt gamma (PG) radiation emitted from a patient during irradiation. The SiPM and scintillation Fiber based Compton Camera (SiFi-CC) setup is being developed for this aim. Approach. A machine learning approach to recognize Compton events is proposed, reconstructing the PG emission profile during proton therapy. The proposed method was verified on pseudo-data generated by a Geant4 simulation for a single proton beam impinging on a polymethyl methacrylate (PMMA) phantom. Three different models including the boosted decision tree (BDT), multilayer perception (MLP) neural network, and k-nearest neighbour (k-NN) were trained using 10-fold cross-validation and then their performances were assessed using the receiver operating characteristic (ROI) curves. Subsequently, after event selection by the most robust model, a software based on the List-Mode Maximum Likelihood Estimation Maximization (LM-MLEM) algorithm was applied for the reconstruction of the PG emission distribution profile. Main results. It was demonstrated that the BDT model excels in signal/background separation compared to the other two. Furthermore, the reconstructed PG vertex distribution after event selection showed a significant improvement in distal falloff position determination. Significance. A highly satisfactory agreement between the reconstructed distal edge position and that of the simulated Compton events was achieved. It was also shown that a position resolution of 3.5 mm full width at half maximum (FWHM) in distal edge position determination is feasible with the proposed setup.

Detection and discrimination of neutron capture events for NCEPT dose quantification

Neutron Capture Enhanced Particle Therapy (NCEPT) boosts the effectiveness of particle therapy by capturing thermal neutrons produced by beam-target nuclear interactions in and around the treatment site, using tumour-specific [Formula: see text]B or [Formula: see text]Gd-based neutron capture agents. Neutron captures release high-LET secondary particles together with gamma photons with energies of 478 keV or one of several energies up to 7.94 MeV, for [Formula: see text]B and [Formula: see text]Gd, respectively. A key requirement for NCEPT's translation is the development of in vivo dosimetry techniques which can measure both the direct ion dose and the dose due to neutron capture. In this work, we report signatures which can be used to discriminate between photons resulting from neutron capture and those originating from other processes. A Geant4 Monte Carlo simulation study into timing and energy thresholds for discrimination of prompt gamma photons resulting from thermal neutron capture during NCEPT was conducted. Three simulated [Formula: see text] mm[Formula: see text] cubic PMMA targets were irradiated by [Formula: see text]He or [Formula: see text]C ion beams with a spread out Bragg peak (SOBP) depth range of 60 mm; one target is homogeneous while the others include [Formula: see text] mm[Formula: see text] neutron capture inserts (NCIs) of pure [Formula: see text]B or [Formula: see text]Gd located at the distal edge of the SOBP. The arrival times of photons and neutrons entering a simulated [Formula: see text] mm[Formula: see text] ideal detector were recorded. A temporal mask of 50-60 ns was found to be optimal for maximising the discrimination of the photons resulting from the neutron capture by boron and gadolinium. A range of candidate detector and thermal neutron shielding materials were simulated, and detections meeting the proposed acceptance criteria (i.e. falling within the target energy window and arriving 60 ns post beam-off) were classified as true or false positives, depending on their origin. The ratio of true/false positives ([Formula: see text]) was calculated; for targets with [Formula: see text]B and [Formula: see text]Gd NCIs, the detector materials which resulted in the highest [Formula: see text] were cadmium-shielded CdTe and boron-shielded LSO, respectively. The optimal irradiation period for both carbon and helium ions was 1 µs for the [Formula: see text]B NCI and 1 ms for the [Formula: see text]Gd NCI.

Towards real-time PGS range monitoring in proton therapy of prostate cancer

Proton therapy of prostate cancer (PCPT) was linked with increased levels of gastrointestinal toxicity in its early use compared to intensity-modulated radiation therapy (IMRT). The higher radiation dose to the rectum by proton beams is mainly due to anatomical variations. Here, we demonstrate an approach to monitor rectal radiation exposure in PCPT based on prompt gamma spectroscopy (PGS). Endorectal balloons (ERBs) are used to stabilize prostate movement during radiotherapy. These ERBs are usually filled with water. However, other water solutions containing elements with higher atomic numbers, such as silicon, may enable the use of PGS to monitor the radiation exposure of the rectum. Protons hitting silicon atoms emit prompt gamma rays with a specific energy of 1.78 MeV, which can be used to monitor whether the ERB is being hit. In a binary approach, we search the silicon energy peaks for every irradiated prostate region. We demonstrate this technique for both single-spot irradiation and real treatment plans. Real-time feedback based on the ERB being hit column-wise is feasible and would allow clinicians to decide whether to adapt or continue treatment. This technique may be extended to other cancer types and organs at risk, such as the oesophagus.

Proton range verification with MACACO II Compton camera enhanced by a neural network for event selection

The applicability extent of hadron therapy for tumor treatment is currently limited by the lack of reliable online monitoring techniques. An active topic of investigation is the research of monitoring systems based on the detection of secondary radiation produced during treatment. MACACO, a multi-layer Compton camera based on LaBr₃ scintillator crystals and SiPMs, is being developed at IFIC-Valencia for this purpose. This work reports the results obtained from measurements of a 150 MeV proton beam impinging on a PMMA target. A neural network trained on Monte Carlo simulations is used for event selection, increasing the signal to background ratio before image reconstruction. Images of the measured prompt gamma distributions are reconstructed by means of a spectral reconstruction code, through which the 4.439 MeV spectral line is resolved. Images of the emission distribution at this energy are reconstructed, allowing calculation of the distal fall-off and identification of target displacements of 3 mm.

Prompt gamma spectroscopy for absolute range verification of 12C ions at synchrotron-based facilities

The physical range uncertainty limits the exploitation of the full potential of charged particle therapy. In this work, we face this issue aiming to measure the absolute Bragg peak position in the target. We investigate p, ⁴He, ¹²C and ¹⁶O beams accelerated at the Heidelberg Ion-Beam Therapy Center. The residual range of the primary ¹²C ions is correlated to the energy spectrum of the prompt gamma radiation. The prompt gamma spectroscopy method was demonstrated for proton beams accelerated by cyclotrons and is developed here for the first time for heavier ions accelerated by a synchrotron. We develop a detector system that includes (i) a spectroscopic unit based on cerium(III) bromide and bismuth germanium oxide scintillating crystals, (ii) a beam trigger based on an array of scintillating fibers and (iii) a data acquisition system based on a FlashADC. We test the system in two different scenarios. In the first series of experiments, we detect and identify 19 independent spectral lines over a wide gamma energy spectrum in the presence of the four ion species for different targets, including a water target with a titanium insert. In the second series of experiments, we introduce a collimator aiming to relate the spectral information to the range of the primary particles. We perform extensive measurements for a ¹²C beam and demonstrate submillimetric precision for the measurement of its Bragg peak position in the experimental setup. The features of the energy and time spectra for gamma radiation induced by p, ⁴He and ¹⁶O are investigated upstream and downstream from the Bragg peak position. We conclude the analysis by extrapolating the required future developments, which would be needed to achieve range verification with a 2 mm accuracy during a single fraction delivery of [Formula: see text] physical dose.

PIBS: Proton and ion beam spectroscopy for in vivo measurements of oxygen, carbon, and calcium concentrations in the human body

Proton and ion beam therapy has proven to benefit tumour control with lower side-effects, mostly in paediatrics. Here we demonstrate a feasible technique for proton and ion beam spectroscopy (PIBS) capable of determining the elemental compositions of the irradiated tissues during particle therapy. This follows the developments in prompt gamma imaging for online range verification and the inheritance from prompt gamma neutron activation analysis. Samples of water solutions were prepared to emulate varying oxygen and carbon concentrations. The irradiation of those samples and other tissue surrogate inserts by protons and ion beams under clinical conditions clearly showed a logarithmic relationship between the target elemental composition and the prompt gamma production. This finding is in line with the known logarithmic dependence of the pH with the proton molar concentration. Elemental concentration changes of 1% for calcium and 2% for oxygen in adipose, brain, breast, liver, muscle and bone-related tissue surrogates were clearly identified. Real-time in vivo measurements of oxygen, carbon and calcium concentrations will be evaluated in a pre-clinical and clinical environment. This technique should have an important impact in the assessment of tumour hypoxia over the course of several treatment fractions and the tracking of calcifications in brain metastases.

Application of the phase-space distribution approach of Monte Carlo for radiation contamination dose estimation from the (n,γ), (γ,n) nuclear reactions and linac leakage photons in the megavoltage radiotherapy facility

Aim

The aim of this study was to characterize the radiation contamination inside and outside the megavoltage radiotherapy room.

Background

Radiation contamination components in the 18 MV linac room are the secondary neutron, prompt gamma ray, electron and linac leakage radiation.

Materials and Methods

An 18 MV linac modeled in a typical bunker employing the MCNPX code of Monte Carlo. For fast calculation, phase-space distribution (PSD) file modeling was applied and the calculations were conducted for the radiation contamination components dose and spectra at 6 locations inside and outside the bunker.

Results

The results showed that the difference of measured and calculated percent depth-dose (PDD) and photo beam-profile (PBP) datasets were lower than acceptable values. At isocenter, the obtained photon dose and neutron fluence were 2.4 × 10^-14 Gy/initial e° and 2.22 × 10^-8 n°/cm², respectively. Then, neutron apparent source strength (Q_N) value was found as 1.34 × 10¹² n°/Gy X at isocenter and the model verified to photon and neutron calculations. A surface at 2 cm below the flattening filter was modeled as phase-space (PS) file for PDD and PBP calculations. Then by use of a spherical cell in the center of the linac target as a PS surface, contaminant radiations dose, fluence and spectra were estimated at 6 locations in a considerably short time, using the registered history of all particles and photons in the 13GB PSD file as primary source in the second step.

Conclusion

Designing the PSD file in MC modeling helps user to solve the problems with complex geometry and physics precisely in a shorter run-time.

MC safe bunker designing for an 18 MV linac with nanoparticles included primary barriers and effect of the nanoparticles on the shielding aspects

Aim

The aim of this study was to design a safe bunker for an 18 MV linac in to configuration; primary barriers made from nanoparticle-containing concrete and pure concrete.

Background

Application of some nanoparticles in the shielding materials has been studied and it was shown that the presence of some nanoparticles improved radiation shielding properties.

Materials and methods

Some percentage of different nanoparticles were modeled by the MCNP5 code of MC in the megavoltage radiotherapy treatment room's primary barriers. Other parts of the designed room, such as secondary barriers and maze door, were modeled as ordinary pure concrete. A safe bunker was designed according to the MC derived spectra at primary and secondary barriers location using a modeled and benchmarked 18 MV linac in free air. Then, the thickness of the required shielding materials for the door and also concrete for the walls and primary barriers were calculated separately.

Results

According to the results, required concrete thickness in primary and secondary barriers was reduced by around 0.8% compared to pure concrete application. Additionally, required lead and BPE decreased by 25% and 15%, respectively, due to primary barriers nanoparticles.

Conclusions

It was concluded that application of some nanoparticles in the shielding materials structures in megavoltage radiotherapy can make the shielding effective.

Range optimization for target and organs at risk in dynamic adaptive passive scattering proton beam therapy - A proof of concept

PURPOSE

The purpose of this study was to design and develop a new range optimization for target and organs at risk (OARs) in dynamic adaptive proton beam therapy (PBT).

METHODS

The new range optimization for target and OARs (RO-TO) was optimized to maintain target dose coverage but not to increase the dose exposure of OARs, while the other procedure, range optimization for target (RO-T), only focused on target dose coverage. A retrospective analysis of a patient who received PBT for abdominal lymph node metastases was performed to show the effectiveness of our new approach. The original plan (OP), which had a total dose of 60 Gy (relative biological effectiveness; RBE), was generated using six treatment fields. Bone-based registration (BR) and tumor-based registration (TR) were performed on each pretreatment daily CT image dataset acquired once every four fractions, to align the isocenter.

RESULTS

Both range adaptive approaches achieved better coverage (D_95%) and homogeneity (D_5%-D_95%) than BR and TR only. However, RO-T showed the greatest increases in D_2cc and D_mean values of the small intestine and stomach and exceeded the limitations of dose exposure for those OARs. RO-TO showed comparable or superior dose sparing compared with the OP for all OARs.

CONCLUSIONS

Our results suggest that BR and TR alone may reduce target dose coverage, and that RO-T may increase the dose exposure to the OARs. RO-TO may achieve the planned dose delivery to the target and OARs more efficiently than the OP. The technique requires testing on a large clinical dataset.

Precision imaging of 4.4 MeV gamma rays using a 3-D position sensitive Compton camera

Imaging of nuclear gamma-ray lines in the 1–10 MeV range is far from being established in both medical and physical applications. In proton therapy, 4.4 MeV gamma rays are emitted from the excited nucleus of either ¹²C* or ¹¹B* and are considered good indicators of dose delivery and/or range verification. Further, in gamma-ray astronomy, 4.4 MeV gamma rays are produced by cosmic ray interactions in the interstellar medium, and can thus be used to probe nucleothynthesis in the universe. In this paper, we present a high-precision image of 4.4 MeV gamma rays taken by newly developed 3-D position sensitive Compton camera (3D-PSCC). To mimic the situation in proton therapy, we first irradiated water, PMMA and Ca(OH)2 with a 70 MeV proton beam, then we identified various nuclear lines with the HPGe detector. The 4.4 MeV gamma rays constitute a broad peak, including single and double escape peaks. Thus, by setting an energy window of 3D-PSCC from 3 to 5 MeV, we show that a gamma ray image sharply concentrates near the Bragg peak, as expected from the minimum energy threshold and sharp peak profile in the cross section of ¹²C(p,p)¹²C*.

Exradin W1 plastic scintillation detector for in vivo skin dosimetry in passive scattering proton therapy

In vivo skin dosimetry is desirable in passive scattering proton therapy because of the possibility of high entrance dose with a small number of fields. However, suitable detectors are needed to determine skin dose in proton therapy. Plastic scintillation detectors (PSDs) are particularly well suited for applications in proton therapy because of their water equivalence, small size, and ease of use. We investigated the utility of the Exradin W1, a commercially available PSD, for in vivo skin dosimetry during passive scattering proton therapy. We evaluated the accuracy of the Exradin W1 in six patients undergoing proton therapy for prostate cancer, as part of an Institutional Review Board-approved protocol. Over 22 weeks, we compared in vivo PSD measurements with in-phantom ionization chamber measurements and doses from the treatment planning system, resulting in 96 in vivo measurements. Temperature and ionization quenching correction factors were applied on the basis of the dose response of the PSD in a phantom. The calibrated PSD exhibited an average 7.8% under-response (±1% standard deviation) owing to ionization quenching. We observed 4% under-response at 37 °C relative to the calibration-temperature response. After temperature and quenching corrections were applied, the overall PSD dose response was within ±1% of the expected dose for all patients. The dose differences between the PSD and ionization chamber measurements for all treatment fields were within ±2% (standard deviation 0.67%). The PSD was highly accurate for in vivo skin dosimetry in passively scattered proton beams and could be useful in verifying proton therapy delivery.

Imaging of monochromatic beams by measuring secondary electron bremsstrahlung for carbon-ion therapy using a pinhole x-ray camera

A feasibility study on the imaging of monochromatic carbon-ion beams for carbon-ion therapy was performed. The evaluation was based on Monte Carlo simulations and beam-irradiation experiments, using a pinhole x-ray camera, which measured secondary electron bremsstrahlung (SEB). The simulation results indicated that the trajectories of the carbon-ion beams with injection energies of 278, 249 and 218 MeV/u in a water phantom, were clearly imaged by measuring the SEB with energies from 30 to 60 keV, using a pinhole camera. The Bragg-peak positions for these three injection energies were located at the positions where the ratios of the counts of SEB acquisitions to the maximum counts were approximately 0.23, 0.26 and 0.29, respectively. Moreover, we experimentally demonstrated that it was possible to identify the Bragg-peak positons, at the positions where the ratios coincided with the simulation results. However, the estimated Bragg-peak positions for the injection energies of 278 and 249 MeV/u were slightly deeper than the expected positions. In conclusion, for both the simulations and experiments, we found that the 25 mm shifts in the Bragg-peak positions can be observed by this method.