Out-of-field measurements and simulations of a proton pencil beam in a wide range of dose rates using a Timepix3 detector: Dose rate, flux and LET

Experimental validation of LET in intensity-modulated proton therapy with a miniaturized pixel detector

Objective.Given the increased interest in incorporating linear energy transfer (LET) as an optimization parameter in intensity-modulated proton therapy (IMPT), a solution for experimental validation of simulations and patient-specific quality assurance (PSQA) in terms of proton LET is needed. Here, we present the methodology and results of LET spectra measurements for spread-out Bragg peak (SOBP) and IMPT plans using a miniaturized pixel detector Timepix3.Approach.We used a MiniPIX Timepix3 detector that provides single-particle tracking, type-resolving power, and spectral information while allowing measurement in quasi-continuous mode. We performed measurements for SOBP and IMPT plans in homogeneous RW3 and heterogeneous CIRS head phantoms with reduced beam current. An artificial intelligence-based model was applied for proton identification and a GPU-accelerated FRED Monte Carlo (MC) code was applied for corresponding MC simulations.Main results.We compared the deposited energy and LET spectra obtained in mixed radiation fields from measurements and MC simulations. The peak positions of deposited energy and LET spectra for the SOBP and IMPT plans agree within the error bars. Discrepancies exceeding the error bars are only visible in the logarithmic scale in high-energy deposition and high-LET tails of the distributions. The mean relative difference of dose-averaged LET values between measurements and MC simulations for individual energy layers is about 5.1%.Significance.This study presents a methodology for assessing radiation quality in proton therapy through energy deposition and LET spectra measurements in uniform and clinical IMPT fields. Findings show an agreement between experimental data and MC simulations, validating our approach. The presented results demonstrate the feasibility of a commercially available Timepix3 detector to validate LET computations in IMPT fields and perform PSQA in terms of LET. This will support the implementation of LET in treatment planning, which will ultimately increase the effectiveness of the treatment.

Configuration of Timepix3 read-out parameters for spectral measurements in proton therapy applications

PURPOSE

With the increasing use of proton therapy, there is a growing emphasis on including radiation quality, often quantified by linear energy transfer, as a treatment plan optimization factor. The Timepix detectors offer energy-sensitive particle tracking useful for the characterization of proton linear energy transfer. To improve the detector's performance in mixed radiation fields produced in proton therapy, we customized the detector settings and performed the per-pixel energy calibration.

METHODS

The detection threshold and per-pixel signal shaping time (IKrum current) were customized, and energy calibration was performed for MiniPIX Timepix3. The detector calibration was verified using α source and clinical proton beams, as well as Monte Carlo simulations. The effects on the detector's performance, in terms of spectral saturation and pixel occupancy, were evaluated.

RESULTS

Measurements with proton beams showed a good agreement with simulations. With the customized settings, the measurable energy range in the detector data-driven mode was extended, and the signal duration time was reduced by 80%, while the yield of pixel time occupancy reduction depends on the number of occupied pixels. For performed measurements with proton beams, the number of occupied pixels was further reduced up to 40% due to the increased threshold.

CONCLUSIONS

Customized detector configuration of the Timepix3 detector allowed for reduced pixel occupancy and mitigation of signal saturation in a data-driven mode without significantly interfering with the energy deposition measurement. The presented approach enables the extension of the operational range, including higher intensities and mixed-radiation fields in particle radiotherapy environments.

Dosimetric and temporal beam characterization of individual pulses in FLASH radiotherapy using Timepix3 pixelated detector placed out-of-field

BACKGROUND

FLASH radiotherapy necessitates the development of advanced Quality Assurance methods and detectors for accurate monitoring of the radiation field. This study introduces enhanced time-resolution detection systems and methods used to measure the delivered number of pulses, investigate temporal structure of individual pulses and dose-per-pulse (DPP) based on secondary radiation particles produced in the experimental room.

METHODS

A 20 MeV electron beam generated from a linear accelerator (LINAC) was delivered to a water phantom. Ultra-high dose-per-pulse electron beams were used with a dose-per-pulse ranging from ̴ 1 Gy to over 7 Gy. The pulse lengths ranged from 1.18 µs to 2.88 µs at a pulse rate frequency of 5 Hz. A semiconductor pixel detector Timepix3 was used to track single secondary particles. Measurements were performed in the air, while the detector was positioned out-of-field at a lateral distance of 200 cm parallel with the LINAC exit window. The dose deposited was measured along with the pulse length and the nanostructure of the pulse.

RESULTS

The time of arrival (ToA) of single particles was measured with a resolution of 1.56 ns, while the deposited energy was measured with a resolution of several keV based on the Time over Threshold (ToT) value. The pulse count measured by the Timepix3 detector corresponded with the delivered values, which were measured using an in-flange integrating current transformer (ICT). A linear response (R2 = 0.999) was established between the delivered beam current and the measured dose at the detector position (orders of nGy). The difference between the average measured and delivered pulse length was ∼0.003(30) μs.

CONCLUSION

This simple non-invasive method exhibits no limitations on the delivered DPP within the range used during this investigation.

Characterizing devices for validation of dose, dose rate, and LET in ultra high dose rate proton irradiations

BACKGROUND

Ultra high dose rate (UHDR) radiotherapy using ridge filter is a new treatment modality known as conformal FLASH that, when optimized for dose, dose rate (DR), and linear energy transfer (LET), has the potential to reduce damage to healthy tissue without sacrificing tumor killing efficacy via the FLASH effect.

PURPOSE

Clinical implementation of conformal FLASH proton therapy has been limited by quality assurance (QA) challenges, which include direct measurement of UHDR and LET. Voxel DR distributions and LET spectra at planning target margins are paramount to the DR/LET-related sparing of organs at risk. We hereby present a methodology to achieve experimental validation of these parameters.

METHODS

Dose, DR, and LET were measured for a conformal FLASH treatment plan involving a 250-MeV proton beam and a 3D-printed ridge filter designed to uniformly irradiate a spherical target. We measured dose and DR simultaneously using a 4D multi-layer strip ionization chamber (MLSIC) under UHDR conditions. Additionally, we developed an "under-sample and recover (USRe)" technique for a high-resolution pixelated semiconductor detector, Timepix3, to avoid event pile-up and to correct measured LET at high-proton-flux locations without undesirable beam modifications. Confirmation of these measurements was done using a MatriXX PT detector and by Monte Carlo (MC) simulations.

RESULTS

MC conformal FLASH computed doses had gamma passing rates of >95% (3 mm/3% criteria) when compared to MatriXX PT and MLSIC data. At the lateral margin, DR showed average agreement values within 0.3% of simulation at 100 Gy/s and fluctuations ∼10% at 15 Gy/s. LET spectra in the proximal, lateral, and distal margins had Bhattacharyya distances of <1.3%.

CONCLUSION

Our measurements with the MLSIC and Timepix3 detectors shown that the DR distributions for UHDR scenarios and LET spectra using USRe are in agreement with simulations. These results demonstrate that the methodology presented here can be used effectively for the experimental validation and QA of FLASH treatment plans.

Time structures of proton pencil beam scanning delivery on a microsecond scale measured with a pixelated semiconductor detector Timepix3

PURPOSE

The time structures of proton spot delivery in proton pencil beam scanning (PBS) radiation therapy are essential in many clinical applications. This study aims to characterize the time structures of proton PBS delivered by both synchrotron and synchrocyclotron accelerators using a non-invasive technique based on scattered particle tracking.

METHODS

A pixelated semiconductor detector, AdvaPIX-Timepix3, with a temporal resolution of 1.56 ns, was employed to measure time of arrival of secondary particles generated by a proton beam. The detector was placed laterally to the high-flux area of the beam in order to allow for single particle detection and not interfere with the treatment. The detector recorded counts of radiation events, their deposited energy and the timestamp associated with the single events. Individual recorded events and their temporal characteristics were used to analyze beam time structures, including energy layer switch time, magnet switch time, spot switch time, and the scanning speeds in the x and y directions. All the measurements were repeated 30 times on three dates, reducing statistical uncertainty.

RESULTS

The uncertainty of the measured energy layer switch times, magnet switch time, and the spot switch time were all within 1% of average values. The scanning speeds uncertainties were within 1.5% and are more precise than previously reported results. The measurements also revealed continuous sub-milliseconds proton spills at a low dose rate for the synchrotron accelerator and radiofrequency pulses at 7 µs and 1 ms repetition time for the synchrocyclotron accelerator.

CONCLUSION

The AdvaPIX-Timepix3 detector can be used to directly measure and monitor time structures on microseconds scale of the PBS proton beam delivery. This method yielded results with high precision and is completely independent of the machine log files.

Particle tracking, recognition and LET evaluation of out-of-field proton therapy delivered to a phantom with implants

Objective.This study aims to assess the composition of scattered particles generated in proton therapy for tumors situated proximal to some titanium (Ti) dental implants. The investigation involves decomposing the mixed field and recording Linear Energy Transfer (LET) spectra to quantify the influence of metallic dental inserts located behind the tumor.Approach.A therapeutic conformal proton beam was used to deliver the treatment plan to an anthropomorphic head phantom with two types of implants inserted in the target volume (made of Ti and plastic, respectively). The scattered radiation resulted during the irradiation was detected by a hybrid semiconductor pixel detector MiniPIX Timepix3 that was placed distal to the Spread-out Bragg peak. Visualization and field decomposition of stray radiation were generated using algorithms trained in particle recognition based on artificial intelligence neural networks (AI NN). Spectral sensitive aspects of the scattered radiation were collected using two angular positions of the detector relative to the beam direction: 0° and 60°.Results.Using AI NN, 3 classes of particles were identified: protons, electrons & photons, and ions & fast neutrons. Placing a Ti implant in the beam's path resulted in predominantly electrons and photons, contributing 52.2% of the total number of detected particles, whereas for plastic implants, the contribution was 65.4%. Scattered protons comprised 45.5% and 31.9% with and without metal inserts, respectively. The LET spectra were derived for each group of particles identified, with values ranging from 0.01 to 7.5 keVμm-1for Ti implants/plastic implants. The low-LET component was primarily composed of electrons and photons, while the high-LET component corresponded to protons and ions.Significance.This method, complemented by directional maps, holds the potential for evaluating and validating treatment plans involving stray radiation near organs at risk, offering precise discrimination of the mixed field, and enhancing in this way the LET calculation.

Linear energy transfer (LET) distribution outside small radiotherapy field edges produced by 6 MV X-rays

In modern radiotherapy with photons, the absorbed dose outside the radiation field is generally investigated. But it is well known that the biological damage depends not only on the absorbed dose but also on LET. This work investigated the dose-average LET (LΔ,D) outside several small radiotherapy fields to provide information that can help for better evaluating the biological effect in organs at risk close to the tumour volume. The electron fluences produced in liquid water by a 6 MV X-rays Varian iX linac were calculated using the EGSnrc Monte Carlo code. With the electron spectra, LΔ,D calculations were made for eight open small square fields and the reference field at water depths of 0.15 cm, 1.35 cm, 9.85 cm and 19.85 cm and several off-axis distances. The variation of LΔ,D from the centre of the beam to 2 cm outside the field's edge depends on the field size and water depth. Using radiobiological data reported in the literature for chromosomal aberrations as an endpoint for the induction of dicentrics determined in Human Lymphocytes, we estimated the maximum low-dose relative biological effectiveness, (RBEM) finding an increase of up to 100% from the centre of the beam to 2 cm from the field's edge.

Thermal neutron detection and track recognition method in reference and out-of-field radiotherapy FLASH electron fields using Timepix3 detectors

Objective.This work presents a method for enhanced detection, imaging, and measurement of the thermal neutron flux.Approach. Measurements were performed in a water tank, while the detector is positioned out-of-field of a 20 MeV ultra-high pulse dose rate electron beam. A semiconductor pixel detector Timepix3 with a silicon sensor partially covered by a6LiF neutron converter was used to measure the flux, spatial, and time characteristics of the neutron field. To provide absolute measurements of thermal neutron flux, the detection efficiency calibration of the detectors was performed in a reference thermal neutron field. Neutron signals are recognized and discriminated against other particles such as gamma rays and x-rays. This is achieved by the resolving power of the pixel detector using machine learning algorithms and high-resolution pattern recognition analysis of the high-energy tracks created by thermal neutron interactions in the converter.Main results. The resulting thermal neutrons equivalent dose was obtained using conversion factor (2.13(10) pSv·cm2) from thermal neutron fluence to thermal neutron equivalent dose obtained by Monte Carlo simulations. The calibrated detectors were used to characterize scattered radiation created by electron beams. The results at 12.0 cm depth in the beam axis inside of the water for a delivered dose per pulse of 1.85 Gy (pulse length of 2.4μs) at the reference depth, showed a contribution of flux of 4.07(8) × 103particles·cm-2·s-1and equivalent dose of 1.73(3) nSv per pulse, which is lower by ∼9 orders of magnitude than the delivered dose.Significance. The presented methodology for in-water measurements and identification of characteristic thermal neutrons tracks serves for the selective quantification of equivalent dose made by thermal neutrons in out-of-field particle therapy.

A Monte-Carlo-based study of a single-2D-detector proton-radiography system

PURPOSE

To assess the feasibility of a proton radiography (pRG) system based on a single thin pixelated detector for water-equivalent path length (WEPL) and relative stopping power (RSP) measurements.

METHODS

A model of a pRG system consisting of a single pixelated detector measuring energy deposition and proton fluence was investigated in a Geant4-based Monte Carlo study. At the position directly after an object traversed by a broad proton beam, spatial 2D distributions are calculated of the energy deposition in, and the number of protons entering the detector. Their ratio relates to the 2D distribution of the average stopping power of protons in the detector. The system response is calibrated against the residual range in water of the protons to provide the 2D distribution of the WEPL of the object. The WEPL distribution is converted into the distribution of the RSP of the object. Simulations have been done, where the system has been tested on 13 samples of homogeneous materials of which the RSPs have been calculated and compared with RSPs determined from simulations of residual-range-in-water, which we refer to as reference RSPs.

RESULTS

For both human-tissue- and non-human-tissue-equivalent materials, the RSPs derived with the detector agree with the reference values within 1%.

CONCLUSION

The study shows that a pRG system based on one thin pixelated detection screen has the potential to provide RSP predictions with an accuracy of 1%.