Calibration method and performance of a time-of-flight detector to measure absolute beam energy in proton therapy

BACKGROUND

The beam energy is one of the most significant parameters in particle therapy since it is directly correlated to the particles' penetration depth inside the patient. Nowadays, the range accuracy is guaranteed by offline routine quality control checks mainly performed with water phantoms, 2D detectors with PMMA wedges, or multi-layer ionization chambers. The latter feature low sensitivity, slow collection time, and response dependent on external parameters, which represent limiting factors for the quality controls of beams delivered with fast energy switching modalities, as foreseen in future treatments. In this context, a device based on solid-state detectors technology, able to perform a direct and absolute beam energy measurement, is proposed as a viable alternative for quality assurance measurements and beam commissioning, paving the way for online range monitoring and treatment verification.

PURPOSE

This work follows the proof of concept of an energy monitoring system for clinical proton beams, based on Ultra Fast Silicon Detectors (featuring tenths of ps time resolution in 50 μm active thickness, and single particle detection capability) and time-of-flight techniques. An upgrade of such a system is presented here, together with the description of a dedicated self-calibration method, proving that this second prototype is able to assess the mean particles energy of a monoenergetic beam without any constraint on the beam temporal structure, neither any a priori knowledge of the beam energy for the calibration of the system.

METHODS

A new detector geometry, consisting of sensors segmented in strips, has been designed and implemented in order to enhance the statistics of coincident protons, thus improving the accuracy of the measured time differences. The prototype was tested on the cyclotron proton beam of the Trento Protontherapy Center (TPC). In addition, a dedicated self-calibration method, exploiting the measurement of monoenergetic beams crossing the two telescope sensors for different flight distances, was introduced to remove the systematic uncertainties independently from any external reference.

RESULTS

The novel calibration strategy was applied to the experimental data collected at TPC (Trento) and CNAO (Pavia). Deviations between measured and reference beam energies in the order of a few hundreds of keV with a maximum uncertainty of 0.5 MeV were found, in compliance with the clinically required water range accuracy of 1 mm.

CONCLUSIONS

The presented version of the telescope system, minimally perturbative of the beam, relies on a few seconds of acquisition time to achieve the required clinical accuracy and therefore represents a feasible solution for beam commission, quality assurance checks, and online beam energy monitoring.

Uniformity calibration for large area XY strip parallel plate ionization chamber

Demand for large area parallel plate ionization chamber (PPIC) or large area ionization chamber (LAIC) has risen in recent years due to several advantages of the large effective area in monitoring therapeutic radiation beams. PPICs are designed for the measurements of beam profile and dosimetry in radiation therapy quality assurance (QA) procedures.Objective. Heterogeneous responses over the large sensitive area pose an undeniable concern for the straightforward applications of PPICs in clinical dosimetry. Uniformity calibration for the detector response is thus essential for the accurate performance of each PPIC unit.Approach.A large area XY strip PPIC, characterized by a large effective area of 345.44 × 345.44 mm²and 256 readout channels, was investigated in this study. A new systematic uniformity calibration is developed to improve the lateral response of the PPIC over the measurements for both narrow beams and large square field beams. A 2D response map of the PPIC was obtained by a spot-scanning method using a compact x-ray tube (mini x-ray). The mini x-ray, providing stable radiation (uncertainty <0.1%), was moved with a step size of 20 mm in 2 dimensions across the entire PPIC surface to collect a complete spot scan. Different uniformity calibration methods were introduced for the measurement of the PPIC by adopting the information from the detector 2D response map.Main results.Deviation of the detector response, before calibration, was observed to reach about 7% for the testing PPIC unit which is much higher than the recommended uniformity response of 1% (IAEA TRS-398). The uniformity response of the PPIC improved significantly to less than 1% across the detector surface after calibration.Significance.The proposed methods enable the practical application of PPIC in routine clinical dosimetry and can be reliably adopted by any radiation facility to perform daily and monthly QA.

Clinical dose assessment for scanned carbon-ion radiotherapy using linear energy transfer measurements and Monte Carlo simulations

Objective. Dosimetric commissioning of treatment planning systems (TPS) focuses on validating the agreement of the physical dose with experimental data. For carbon-ion radiotherapy, the commissioning of the relative biological effectiveness (RBE) is necessary to predict the clinical outcome based on the radiation quality of the mixed radiation field. In this study, we proposed a approach for RBE commissioning using Monte Carlo (MC) simulations, which was further strengthen by RBE validation based on linear energy transfer (LET) measurements.Approach. First, we tuned the MC simulation based on the results of dosimetric experiments including the beam ranges, beam sizes, and MU calibrations. Furthermore, we compared simulated results to measured depth- and radial-LET distributions of the 430 MeV u^-1carbon-ion spot beam with a 1.5 mm², 36μm thick silicon detector. The measured dose-averaged LET (LET_d) and RBE were compared with the simulated results. The RBE was calculated based on the mixed beam model with linear-quadratic parameters depending on the LET. Finally, TPS-calculated clinical dose profiles were validated through the tuned MC-based calculations.Main results. A 10 keVμm^-1and 0.15 agreement for LET_dand RBE, respectively, were found between simulation and measurement results obtained for a 2σlateral size of 430 MeV u^-1carbon-ion spot beam in water. These results suggested that the tuned MC simulation can be used with acceptable precision for the RBE and LET calculations of carbon-ion spot beam within the clinical energy range. For physical and clinical doses, the TPS- and MC-based calculations showed good agreements within 1.0% at the centre of the spread-out Bragg peaks.Significance. The tuned MC simulation can accurately reproduce the actual carbon-ion beams, and it can be used to validate the physical and clinical dose distributions calculated by TPS. Moreover, the MC simulation can be used for dosimetric commissioning, including clinical doses, without LET measurements.

A high spatiotemporal resolution 2D strip ionization chamber array for proton pencil beam scanning FLASH radiotherapy

PURPOSE

Experimental measurements of 2D dose rate distributions in proton pencil beam scanning (PBS) FLASH radiation therapy (RT) are currently lacking. In this study, we characterize a newly designed 2D strip-segmented ionization chamber array (SICA) with high spatial and temporal resolution and demonstrate its applications in a modern proton PBS delivery system at both conventional and ultra-high dose rates.

METHODS

A dedicated research beamline of the Varian ProBeam system was employed to deliver a 250 MeV proton PBS beam with nozzle currents up to 215 nA. In the research and clinical beamlines, the spatial, temporal, and dosimetric performance of the SICA was characterized and compared with measurements using parallel-plate ion chambers (IBA PPC05 and PTW Advanced Markus chamber), a 2D scintillator camera (IBA Lynx), Gafchromic films (EBT-XD), and a Faraday Cup. A novel reconstruction approach was proposed to enable the measurement of 2D dose and dose rate distributions using such a strip-type detector.

RESULTS

The SICA demonstrated a position accuracy of 0.12 ± 0.02 mm at a 20 kHz sampling rate (50 μs per event) and a linearity of R² > 0.99 for both dose and dose rate with nozzle beam currents ranging from 1 nA to 215 nA. The 2D dose comparison to the film measurement resulted in a gamma passing rate of 99.8% (2 mm/2%). A measurement-based proton PBS 2D FLASH dose rate distribution was compared to simulation results and showed a gamma passing rate of 97.3% (2 mm/2%).

CONCLUSIONS

The newly designed SICA demonstrated excellent spatial, temporal, and dosimetric performance and is well suited for commissioning, quality assurance (QA), and a wide range of clinical applications in proton PBS clinical and FLASH radiotherapy. This article is protected by copyright. All rights reserved.

Proton therapy monitoring: spatiotemporal emission reconstruction with prompt gamma timing and implementation with PET detectors

Objective. In this study we introduce spatiotemporal emission reconstruction prompt gamma timing (SER-PGT), a new method to directly reconstruct the prompt photon emission in the space and time domains inside the patient in proton therapy.Approach. SER-PGT is based on the numerical optimisation of a multidimensional likelihood function, followed by a post-processing of the results. The current approach relies on a specific implementation of the maximum-likelihood expectation maximisation algorithm. The robustness of the method is guaranteed by the complete absence of any information about the target composition in the algorithm.Main results. Accurate Monte Carlo simulations indicate a range resolution of about 0.5 cm (standard deviation) when considering 10⁷primary protons impinging on an homogeneous phantom. Preliminary results on an anthropomorphic phantom are also reported.Significance. By showing the feasibility for the reconstruction of the primary particle range using PET detectors, this study provides significant basis for the development of an hybrid in-beam PET and prompt photon device.

Gaussian fitting algorithm with multi-geometric parameters for rotated elliptical beam profiling using pixel ion chamber

PURPOSE

A high-precision rotated elliptical beam profiling method based on pixel ion chamber is proposed in this paper. This method aims to improve the accuracy by modeling the transverse profile of rotated beam as an ellipse with additional correlation coefficient and eliminating the fitting error due to the volume averaging effect of pixel ion chamber.

METHODS

In pencil beam scanning (PBS) proton therapy systems, the transverse beam profile model is generally represented as a standard Gaussian distribution. Considering the elliptical spots, two-dimensional (2D) joint Gaussian distribution characterized with the correlation coefficient ρ is adopted in this study. Gaussian-type particle distribution with white noise was generated and processed in MATLAB to simulate the secondary particle collection in the pixel ion chamber. The simulated pixel ion chamber is a commercially available ion chamber which consists of 12 × 12 small square pixels (3.75 × 3.75 mm² ) with a 0.05 mm interval. The simulated signals were preprocessed by filtering with the noise threshold and extracting the maximum simply connected domain (MSCD) of the signal. Then, five geometric parameters that identify the transverse beam profiles were fitted under different signal-to-noise ratio (SNR) conditions: the center of the beam (x₀ , y₀ ), the spot size (σ_major , σ_minor ), and the rotation angle θ formed between the major axes of elliptical spot and the x axes of the ion chamber. First, the simulated signals were preprocessed by filtering with the noise threshold and extracting the MSCD of the signal. Second, a rectification curve of systematic error in fitted spot size versus the prescribed spot size was used to predict the systematic error due to the volume averaging effect. Finally, the effects of fitting errors on therapeutic dose were evaluated in terms of gamma index and relative dose difference.

RESULTS

When the SNR is not less than 20 dB, the relative fitting error of spot size and the absolute fitting error of angle θ are less than 1% and 6.1°, respectively. The fitting error of beam center increases with spot size and will not exceed 0.22 mm when spot size reaches up to 12 mm. At a SNR equal to 20 dB, neither cold nor hot spots were presented in dose distribution calculated with the fitted spot parameters.

CONCLUSION

The improved Gaussian fitting algorithm performs well when SNR is not less than 20 dB. This method can effectively distinguish the nominal beam and rotated elliptical beam. An ideal systematic error curve can be predicted and used to correct the fitted spot size, thus eliminating the systematic error due to the volume averaging effect of the pixel ion chamber. The fitting error of spot size cannot be fully corrected, but it is negligible and shows little effect on the overall therapeutic dose.

Characterization of a high-resolution 2D transmission ion chamber for independent validation of proton pencil beam scanning of conventional and FLASH dose delivery

INTRODUCTION

Ultra-high dose rate (FLASH) radiotherapy has become a popular research topic with the potential to reduce normal tissue toxicities without losing the benefit of tumor control. The development of FLASH proton pencil beam scanning (PBS) delivery requires accurate dosimetry despite high beam currents with correspondingly high ionization densities in the monitoring chamber. In this study, we characterized a newly designed high-resolution position sensing transmission ionization chamber with a purpose-built multichannel electrometer for both conventional and FLASH dose rate proton radiotherapy.

METHODS

The dosimetry and positioning accuracies of the ion chamber were fully characterized with a clinical scanning beam. On the FLASH proton beamline, the cyclotron output current reached up to 350 nA with a maximum energy of 226.2 MeV, with 210 ± 3 nA nozzle pencil beam current. The ion recombination effect was characterized under various bias voltages up to 1000 V and different beam intensities. The charge collected by the transmission ion chamber was compared with the measurements from a Faraday cup.

RESULTS

Cross-calibrated with an Advanced Markus chamber (PTW, Freiburg, Germany) in a uniform PBS proton beam field at clinical beam setting, the ion chamber calibration was 38.0 and 36.7 GyE·mm² /nC at 100 and 226.2 MeV, respectively. The ion recombination effect increased with larger cyclotron current at lower bias voltage while remaining ≤0.5 ± 0.5% with ≥200 V of bias voltage. Above 200 V, the normalized ion chamber readings demonstrated good linearity with the mass stopping power in air for both clinical and FLASH beam intensities. The spot positioning accuracy was measured to be 0.10 ± 0.08 mm in two orthogonal directions.

CONCLUSION

We characterized a transmission ion chamber system under both conventional and FLASH beam current densities and demonstrated its suitability for use as a proton pencil beam dose and spot position delivery monitor under FLASH dose rate conditions.

Monte Carlo computation of 3D distributions of stopping power ratios in light ion beam therapy using GATE-RTion

PURPOSE

This paper presents a novel method for the calculation of three-dimensional (3D) Bragg-Gray water-to-detector stopping power ratio (s_w,det ) distributions for proton and carbon ion beams.

METHODS

Contrary to previously published fluence-based calculations of the stopping power ratio, the s_w,det calculation method used in this work is based on the specific way GATE/Geant4 scores the energy deposition. It only requires the use of the so-called DoseActor, as available in GATE, for the calculation of the s_w,det at any point of a 3D dose distribution. The simulations are performed using GATE-RTion v1.0, a dedicated GATE release that was validated for the clinical use in light ion beam therapy.

RESULTS

The Bragg-Gray water-to-air stopping power ratio (s_w,air ) was calculated for monoenergetic proton and carbon ion beams with the default stopping power data in GATE-RTion v1.0 and the new ICRU90 recommendation. The s_w,air differences between the use of the default and the ICRU90 configuration were 0.6% and 5.4% at the physical range (R₈₀ - 80% dose level in the distal dose fall-off) for a 70 MeV proton beam and a 120 MeV/u carbon ion beam, respectively. For protons, the s_w,det results for lithium fluoride, silicon, gadolinium oxysulfide, and the active layer material of EBT2 (radiochromic film) were compared with the literature and a reasonable agreement was found. For a real patient treatment plan, the 3D distributions of s_w,det in proton beams were calculated.

CONCLUSIONS

Our method was validated by comparison with available literature data. Its equivalence with Bragg-Gray cavity theory was demonstrated mathematically. The capability of GATE-RTion v1.0 for the s_w,det calculation at any point of a 3D dose distribution for simple and complex proton and carbon ion plans was presented.

A new detector for the beam energy measurement in proton therapy: a feasibility study

The proof of concept of a new device, capable of determining in a few seconds the energy of clinical proton beams by measuring the time of flight (ToF) of protons, is presented. The prototype consists of two thin ultra fast silicon detector (UFSD) pads, aligned along the beam direction in a telescope configuration and readout by a digitizer. The method developed for extracting the energy at the isocenter from the measured ToF, validated by Monte Carlo simulations, and the procedure used to calibrate the system are also presented and discussed in detail. The prototype was tested at the Centro Nazionale di Adroterapia Oncologica (CNAO, Pavia, Italy), at several beam energies, covering the entire clinical range, and using different distances between the sensors. The measured beam energies were benchmarked against the nominal CNAO energy values, obtained during the commissioning of the centre from the measured ranges in water. Deviations of few hundreds of keV have been achieved for all considered proton beam energies for distances between the two sensors larger than 60 cm, indicating a sensitivity to the corresponding beam range in water smaller than the clinical tolerance of 1 mm. Moreover, few seconds of irradiation were necessary to collect the required statistics. These preliminary results indicate that a telescope of UFSDs could achieve in a short time the accuracy required for the clinical application and therefore encourage further investigations towards the improvement and the optimization of the present prototype.

Characterization of a Low-Cost Plastic Fiber Array Detector for Proton Beam Dosimetry

The Pencil Beam Scanning (PBS) technique in proton therapy uses fast magnets to scan the tumor volume rapidly. Changing the proton energy allows changing to layers in the third dimension, hence scanning the same volume several times. The PBS approach permits adapting the speed and/or current to modulate the delivered dose. We built a simple prototype that measures the dose distribution in a single step. The active detection material consists of a single layer of scintillating fibers (i.e., 1D) with an active length of 100 mm, a width of 18.25 mm, and an insignificant space (20 μm) between them. A commercial CMOS-based camera detects the scintillation light. Short exposure times allow running the camera at high frame rates, thus, monitoring the beam motion. A simple image processing method extracts the dose information from each fiber of the array. The prototype would allow scaling the concept to multiple layers read out by the same camera, such that the costs do not scale with the dimensions of the fiber array. Presented here are the characteristics of the prototype, studied under two modalities: spatial resolution, linearity, and energy dependence, characterized at the Center for Proton Therapy (Paul Scherrer Institute); the dose rate response, measured at an electron accelerator (Swiss Federal Institute of Metrology).

Accuracy assessment of the CNAO dose delivery system in the initial period of clinical activity and impact of later improvements on delivered dose distributions

PURPOSE

A retrospective analysis of the dose delivery system (DDS) performances of the initial clinical operation at CNAO (Centro Nazionale di Adroterapia Oncologica) is reported, and compared with the dose delivery accuracy following the implementation of a position feedback control.

METHODS

Log files and raw data of the DDS were analyzed for every field of patients treated with protons and carbon ions between January 2012 and April 2013 (~3800 fields). To investigate the DDS accuracy, the spot positions and the number of particles per spot measured by the DDS and prescribed by the treatment planning system were compared for each field. The impact of deviations on dose distributions was studied by comparing, through the gamma-index method, 2 three-dimensional (3D) physical dose maps (one for prescribed, one for measured data), generated by a validated dose computation software. The maximum gamma and the percentage of points with gamma ≤ 1 (passing volume) were studied as a function of the treatment day, and correlated with the deviations from the prescription in the measured number of particles and spot positions. Finally, delivered dose distributions of same treatment plans were compared before and after the implementation of a feedback algorithm for the correction of small position deviations, to study the effect on the delivery quality. A double comparison of prescribed and measured 3D maps, before and after feedback implementation, is reported and studied for a representative treatment delivered in 2012, redelivered on a polymethyl methacrylate (PMMA) block in 2018.

RESULTS

Systematic deviations of spot positions, mainly due to beam lateral offsets, were always found within 1.5 mm, with the exception of the initial clinical period. The number of particles was very stable, as possible deviations are exclusively related to the quantization error in the conversion from monitor counts to particles. For the chosen representative patient treatment, the gamma-index evaluation of prescribed and measured dose maps, before and after feedback implementation, showed a higher variability of maximum gamma for the 2012 irradiation, with respect to the reirradiation of 2018. However, the 2012 passing volume is >99.8% for the sum of all fields, which is comparable to the value of 2018, with the exception of one day with 98.2% passing volume, probably related to an instability of the accelerating system.

CONCLUSIONS

A detailed retrospective analysis of the DDS performances in the initial period of CNAO clinical activity is reported. The spot position deviations are referable to beam lateral offset fluctuations, while almost no deviation was found in the number of particles. The impact of deviations on dose distributions showed that the position feedback implementation and the increased beam control capability acquired after the first years of clinical experience led to an evident improvement in the DDS stability, evaluated in terms of gamma-index as a measure of the impact on dose distributions. However, the clinical effect of the maximum gamma variability found in the 2012 representative irradiation is mitigated by averaging along the number of fractions, and the high percentage of passing volumes confirmed the accuracy of the delivery even before the feedback implementation.