Commissioning and initial experience with the first clinical gantry-mounted proton therapy system

Commissioning of a novel gantry-less proton therapy system

Purpose

The focus of this article is to describe the configuration, testing, and commissioning of a novel gantry-less synchrotron-based proton therapy (PT) facility.

Materials and methods

The described PT system delivers protons with a water equivalent range between 4 and 38 cm in 1800 energy layers. The fixed beam delivery permits a maximum field size of 28 × 30 cm2. The patient positioning and imaging system includes a six-degree-of-freedom robotic arm, a convertible patient chair, a vertical 4DCT, and an orthogonal 2D X-ray imaging system.

Results

The spot positioning reproducibility was consistent within ±1 mm. The width (σ) of the beam profile at the isocenter was energy dependent and ranged from 2.8 mm to 7.7 mm. Absolute dose reproducibility was measured and deviations were found to be <0.62% for all possible beam scenarios. The built-in dose monitoring system was successfully tested for its ability to generate interlocks under specific conditions (beam spot deviation ≥2 mm, individual spot dose ≥10% or ≥0.25 Gy, spot energy deviation ≥0.5 MeV). The robot positioning exhibited a consistent reproducibility within ±1 mm. All tested scenarios achieved laser-free initial 3D/3D image-guided positioning within ±5 mm. Subsequent 2D/3D positioning showed an accuracy of ±1 mm. A single 2D/3D image registration event corrected positions in all cases. Results of gamma analysis (3%, 3 mm) demonstrated pass rates greater than 95% for head and neck, thorax, abdomen treatment plans.

Conclusions

We report on the performance of a novel single-room gantry-less PT system comprised of a compact synchrotron and an adjustable (from nearly horizontal to almost vertical) patient positioning system. The commissioning results show high accuracy and reproducibility of the main proton beam parameters and the patient positioning system. The new PT facility started patient treatments in March 2023, which were the first in Israel and the Middle Eastern region.

Sectored single-energy volumetric-modulated proton arc therapy (VPAT): A preliminary multi-disease-site concept study

PURPOSE

To explore the feasibility of a novel intensity-modulated proton arc technique that uses a single-energy beam from the cyclotron. The beam energy is externally modulated at each gantry angle by a tertiary energy modulator (EM). We hypothesize that irradiating in an arc without requiring an energy change from the cyclotron will achieve a faster delivery (main advantage of our technique) while keeping clinically desirable dosimetric results.

METHODS

In a retrospective cohort of four patients with female pelvis, prostate, lung, and brain cancers, we investigated our volumetric-modulated proton arc therapy (VPAT) technique. Arcs were simulated by sectors of 1°-spaced static beams. Keeping the energy requested from the cyclotron the same for each entire arc was supported by a predesigned EM placed in front of the nozzle. As a feasibility measure, EM thicknesses were calculated. Delivery times and doses to targets and organs at risk (OARs) were compared to those of the clinical plans.

RESULTS

VPAT plans were comparable to their clinical counterparts in achieving target dose conformity, being robust to uncertainties, and meeting clinical dose-volume constraints. Cyclotron energies for the four cases were within 159-220 MeV, and energy modulation range was 69-100 MeV, equivalent to 13-19 cm of water-equivalent thickness (WET). Plan delivery times were reduced from > 5 min in our clinical practice to < 3.5 min in VPAT.

CONCLUSION

For the evaluated plans, the novel VPAT approach achieved shorter delivery times without sacrificing robustness, OAR sparing or target coverage. VPAT's EMs had WETs implementable in a clinical setup.

ESTRO-EPTN radiation dosimetry guidelines for the acquisition of proton pencil beam modelling data

Proton therapy (PT) is an advancing radiotherapy modality increasingly integrated into clinical settings, transitioning from research facilities to hospital environments. A critical aspect of the commissioning of a proton pencil beam scanning delivery system is the acquisition of experimental beam data for accurate beam modelling within the treatment planning system (TPS). These guidelines describe in detail the acquisition of proton pencil beam modelling data. First, it outlines the intrinsic characteristics of a proton pencil beam-energy distribution, angular-spatial distribution and particle number. Then, it lists the input data typically requested by TPSs. Finally, it describes in detail the set of experimental measurements recommended for the acquisition of proton pencil beam modelling data-integrated depth-dose curves, spot maps in air, and reference dosimetry. The rigorous characterization of these beam parameters is essential for ensuring the safe and precise delivery of proton therapy treatments.

Proton Therapy Equipment Installation, Upgrades, and Building Design

PURPOSE

This work aims at reviewing challenges and pitfalls in proton facility design related to equipment upgrade or replacement. Proton therapy was initially developed at research institutions in the 1950s which ushered in the use of hospital-based machines in 1990s. We are approaching an era where older commercial machines are reaching the end of their life and require replacement. The future widespread application of proton therapy depends on cost reduction; customized building design and installation are significant expenses.

METHODS AND MATERIALS

We take this opportunity to discuss how commercial proton machines have been installed and how buildings housing the equipment have been designed.

RESULTS

Data on dimensions and weights of the larger components of proton systems (cyclotron main magnet and gantries) are presented and innovative, non-gantry-based, patient positioning systems are discussed.

CONCLUSIONS

We argue that careful consideration of the building design to include larger elevators, hoistways from above, wide corridors and access slopes to below grade installations, generic vault and treatment room layouts to accommodate multiple vendor's equipment, and modular system design can provide specific benefits during planning, installation, maintenance, and replacement phases of the project. Room temperature magnet coils can be constructed in a more modular manner: a potential configuration is presented. There is scope for constructing gantries and magnet yokes from smaller modular sub-units. These considerations would allow a hospital to replace a commercial machine at its end of life in a manner similar to a linac.

A bi-directional beam-line energy ramping for efficient patient treatment with scanned proton therapy

Objective.The treatment of mobile tumours using Pencil Beam Scanning (PBS) has become more prevalent in the last decade. However, to achieve the same beam delivery quality as for static tumours, treatments have to be combined with motion mitigation techniques, not limited but including, breath hold, gating and re-scanning, which typically prolong treatment time. In this article we present a novel method of bi-directional energy modulation and demonstrate our initial experience in improvement of treatment efficiency. Approach.At Paul Scherrer Institute Gantry 2 mobile tumours are treated by combining PBS with gating and volumetric re-scanning (VR), where the target volume is irradiated multiple times. Initial implementation of VR used only descending beam energies, creating a substantial dead time due to the beam-line initialization (ramping) before each re-scan. In 2019 we commissioned an energy meandering strategy that allows us to avoid beam line ramping in-between energy series while maintaining beam delivery quality.Main results.The measured beam parameters difference for both energy sequence are in the order of the typical daily variations: 0.2 mm in beam position and 0.2 mm in range. Using machine log files, we performed point-to-point dose difference calculations between original and new applications where we observed dose differences of less than 2%. After three years of operation employing bi-directional energy modulation, we have analysed the individual beam delivery time for 181 patients and have compared this to simulations of the timing behaviour assuming uni-directional energy sequence application. Depending on treatment complexity, we obtained plan delivery time reductions of up to 55%, with a median time gain of 17% for all types of treatments.Significance. Bi-directional energy modulation can help improving patient treatment efficiency by reducing delivery times especially for complex and specialised irradiations. It could be implemented in many existing facilities without significant additional hardware upgrades.

Synchronized high-speed scintillation imaging of proton beams, generated by a gantry-mounted synchrocyclotron, on a pulse-by-pulse basis

BACKGROUND

With the emergence of more complex and novel proton delivery techniques, there is a need for quality assurance (QA) tools with high spatiotemporal resolution to conveniently measure the spatial and temporal properties of the beam. In this context, scintillation-based dosimeters, if synchronized with the radiation beam and corrected for ionization quenching, are appealing.

PURPOSE

To develop a synchronized high-speed scintillation imaging system for characterization and verification of the proton therapy beams on a pulse-by-pulse basis.

MATERIALS AND METHODS

A 30 cm × 30 cm × 5 cm block of BC-408 plastic scintillator placed in a light-tight housing was irradiated by proton beams generated by a Mevion S250^TM proton therapy synchrocyclotron. A high-speed camera system, placed perpendicular to the beam direction and facing the scintillator, was synchronized to the accelerator's pulses to capture images. Opening and closing of the camera's shutter was controlled by setting a proper time delay and exposure time, respectively. The scintillation signal was recorded as a set of two-dimensional (2D) images. Empirical correction factors were applied to the images to correct for the non-uniformity of the pixel sensitivity and quenching of the scintillator. Proton range and modulation were obtained from the corrected images.

RESULTS

The camera system was able to capture all data on a pulse-by-pulse basis at a rate of ∼504 frames per second. The applied empirical correction method for ionization quenching was effective and the corrected composite image provided a 2D map of dose distribution. The measured range (depth of distal 90%) through scintillation imaging agreed within 1.2 mm with that obtained from ionization chamber measurement.

CONCLUSION

A high-speed camera system capable of capturing scintillation signals from individual proton pulses was developed. The scintillation imaging system is promising for rapid proton beam characterization and verification. This article is protected by copyright. All rights reserved.

Brainstem Toxicity in Pediatric Patients Treated with Protons Using a Single-vault Synchrocyclotron System

Purpose

Cranial radiation therapy remains an integral component of curative treatment for pediatric patients with brain tumors. Proton beam radiation therapy (PBT) can limit collateral radiation dose to surrounding normal tissue, thus reducing off-target exposure while maintaining appropriate tumor coverage. While PBT offers significant advantages over photon therapy for pediatric patients with intracranial malignancies, cases of brainstem necrosis after PBT have raised concerns that PBT may pose an increased risk of necrosis over photon therapy. We investigated the incidence of brainstem necrosis at our institution in children treated with PBT for intracranial malignancies.

Patients and Methods

Patients with pediatric brain tumor treated with passively scattered PBT, using a gantry-mounted, synchrocyclotron single-vault system between 2013 and 2018, were retrospectively reviewed. Inclusion criteria included patients 21 years of age or younger who received a minimum 0.1 cm³ maximum brainstem dose of 50 Gray relative biological effectiveness (GyRBE). Patients were assessed for “central nervous system necrosis” in the brainstem per the Common Terminology Criteria for Adverse Events (CTCAE), version 5.0 (US National Cancer Institute, Bethesda, Maryland) criteria.

Results

Fifty-eight patients were included for analysis. The median age was 10.3 years. Twenty-one (36.2%) patients received craniospinal irradiation. Thirty-four (58.6%) patients received chemotherapy. The median prescription radiation dose was 54 GyRBE. Regarding published dosimetric constraints used at 3 separate proton centers, the goal brainstem D50% <52 GyRBE was exceeded in 23 (40%) patients, but the brainstem Dmax <58 GyRBE was not exceeded in any patients. No patient experienced grade ≥2 brainstem injury. One patient demonstrated radiographic changes consistent with grade 1 toxicity. This patient had myeloablative chemotherapy with tandem stem cell rescue before PBT.

Conclusion

Our data demonstrates a low risk of any brainstem injury in children treated with passively scattered PBT using a single-vault synchrocyclotron.

Towards subpercentage uncertainty proton stopping-power mapping via dual-energy CT: Direct experimental validation and uncertainty analysis of a statistical iterative image reconstruction method

PURPOSE

To assess the potential of a joint dual-energy computerized tomography (CT) reconstruction process (statistical image reconstruction method built on a basis vector model (JSIR-BVM)) implemented on a 16-slice commercial CT scanner to measure high spatial resolution stopping-power ratio (SPR) maps with uncertainties of less than 1%.

METHODS

JSIR-BVM was used to reconstruct images of effective electron density and mean excitation energy from dual-energy CT (DECT) sinograms for 10 high-purity samples of known density and atomic composition inserted into head and body phantoms. The measured DECT data consisted of 90 and 140 kVp axial sinograms serially acquired on a Philips Brilliance Big Bore CT scanner without beam-hardening corrections. The corresponding SPRs were subsequently measured directly via ion chamber measurements on a MEVION S250 superconducting synchrocyclotron and evaluated theoretically from the known sample compositions and densities. Deviations of JSIR-BVM SPR values from their theoretically calculated and directly measured ground-truth values were evaluated for our JSIR-BVM method and our implementation of the Hünemohr-Saito (H-S) DECT image-domain decomposition technique for SPR imaging. A thorough uncertainty analysis was then performed for five different scenarios (comparison of JSIR-BVM stopping-power ratio/stopping power (SPR/SP) to International Commission on Radiation Measurements and Units benchmarks; comparison of JSIR-BVM SPR to measured benchmarks; and uncertainties in JSIR-BVM SPR/SP maps for patients of unknown composition) per the Joint Committee for Guides in Metrology and the Guide to Expression of Uncertainty in Measurement, including the impact of uncertainties in measured photon spectra, sample composition and density, photon cross section and I-value models, and random measurement uncertainty. Estimated SPR uncertainty for three main tissue groups in patients of unknown composition and the weighted proportion of each tissue type for three proton treatment sites were then used to derive a composite range uncertainty for our method.

RESULTS

Mean JSIR-BVM SPR estimates deviated by less than 1% from their theoretical and directly measured ground-truth values for most inserts and phantom geometries except for high-density Delrin and Teflon samples with SPR error relative to proton measurements of 1.1% and -1.0% (head phantom) and 1.1% and -1.1% (body phantom). The overall root-mean-square (RMS) deviations over all samples were 0.39% and 0.52% (head phantom) and 0.43% and 0.57% (body phantom) relative to theoretical and directly measured ground-truth SPRs, respectively. The corresponding RMS (maximum) errors for the image-domain decomposition method were 2.68% and 2.73% (4.68% and 4.99%) for the head phantom and 0.71% and 0.87% (1.37% and 1.66%) for the body phantom. Compared to H-S SPR maps, JSIR-BVM yielded 30% sharper and twofold sharper images for soft tissues and bone-like surrogates, respectively, while reducing noise by factors of 6 and 3, respectively. The uncertainty (coverage factor k = 1) of the DECT-to-benchmark values comparison ranged from 0.5% to 1.5% and is dominated by scanning-beam photon-spectra uncertainties. An analysis of the SPR uncertainty for patients of unknown composition showed a JSIR-BVM uncertainty of 0.65%, 1.21%, and 0.77% for soft-, lung-, and bony-tissue groups which led to a composite range uncertainty of 0.6-0.9%.

CONCLUSIONS

Observed JSIR-BVM SPR estimation errors were all less than 50% of the estimated k = 1 total uncertainty of our benchmarking experiment, demonstrating that JSIR-BVM high spatial resolution, low-noise SPR mapping is feasible and is robust to variations in the geometry of the scanned object. In contrast, the much larger H-S SPR estimation errors are dominated by imaging noise and residual beam-hardening artifacts. While the uncertainties characteristic of our current JSIR-BVM implementation can be as large as 1.5%, achieving < 1% total uncertainty is feasible by improving the accuracy of scanner-specific scatter-profile and photon-spectrum estimates. With its robustness to beam-hardening artifact, image noise, and variations in phantom size and geometry, JSIR-BVM has the potential to achieve high spatial-resolution SPR mapping with subpercentage accuracy and estimated uncertainty in the clinical setting.

Future Developments in Charged Particle Therapy: Improving Beam Delivery for Efficiency and Efficacy

The physical and clinical benefits of charged particle therapy (CPT) are well recognized. However, the availability of CPT and complete exploitation of dosimetric advantages are still limited by high facility costs and technological challenges. There are extensive ongoing efforts to improve upon these, which will lead to greater accessibility, superior delivery, and therefore better treatment outcomes. Yet, the issue of cost remains a primary hurdle as utility of CPT is largely driven by the affordability, complexity and performance of current technology. Modern delivery techniques are necessary but limited by extended treatment times. Several of these aspects can be addressed by developments in the beam delivery system (BDS) which determines the overall shaping and timing capabilities enabling high quality treatments. The energy layer switching time (ELST) is a limiting constraint of the BDS and a determinant of the beam delivery time (BDT), along with the accelerator and other factors. This review evaluates the delivery process in detail, presenting the limitations and developments for the BDS and related accelerator technology, toward decreasing the BDT. As extended BDT impacts motion and has dosimetric implications for treatment, we discuss avenues to minimize the ELST and overview the clinical benefits and feasibility of a large energy acceptance BDS. These developments support the possibility of advanced modalities and faster delivery for a greater range of treatment indications which could also further reduce costs. Further work to realize methodologies such as volumetric rescanning, FLASH, arc, multi-ion and online image guided therapies are discussed. In this review we examine how increased treatment efficiency and efficacy could be achieved with improvements in beam delivery and how this could lead to faster and higher quality treatments for the future of CPT.

Development of a storage phosphor imaging system for proton pencil beam spot profile determination

PURPOSE

Accurate two-dimensional (2D) profile measurements at submillimeter precision are necessary for proton beam commissioning and periodic quality assurance (QA) purposes and are currently performed at our institution with a commercial scintillation detector (Lynx PT) with limited means for independent checks. The purpose of this work was to create an independent dosimetry system consisting of an in-house optical scanner and a BaFBrI:Eu²⁺ storage phosphor dosimeter by: (a) determining the optimal settings for the optical scanner, (b) measuring 2D proton spot profiles with the storage phosphors, and (c) comparing them to similar measurements using a commercial scintillation detector.

METHODS

An in-house 2D laboratory optical scanner was constructed and spatially calibrated for accurate 2D photostimulated luminescence (PSL) dosimetry. Square 5 × 5 cm² BaFBrI:Eu²⁺ dosimeter samples were uniformly irradiated with line scans performed to determine the physical and electronic scanner settings resulting in the highest signal-to-noise ratios (SNR) at a sub-millimeter spatial resolution. The resultant spatial resolution of the scanner was then quantitatively assessed by measuring (a) line pairs on a standard X-ray lead bar phantom and (b) modulation transfer functions. Following this, 2D proton spot profiles from a Mevion S250i Hyperscan proton unit were obtained at 1, 10, 20, 30, 40, and 50 monitor unit (MU) settings at maximum energy (E₀  = 227.1 MeV) and compared to baseline profiles from a commercial scintillation detector, where 1 MU is calibrated to deliver 1 Gy absolute proton dose-to-water under reference conditions, that is, 41 × 41 proton spots uniformly spaced by 0.25 cm within a 10 × 10 cm² square field size at maximum energy (227.1 MeV) in water at depth of 5 cm at isocenter. The dosimetric system's sensitivities to (a) ±1 mm positional shifts and (b) ±0.3 mm beam lateral spread changes were quantitatively evaluated through a Gaussian fitting of the crossline and inline plots of the respective artificially shifted beam profiles.

RESULTS

The physical scanner settings of (a) Δτ = 27 ms time interval between data samples, (b) v_x  = 1.235 cm/s scanning speed, (c) 1% laser transmission (0.02 mW power) and (d) (Δx, Δy) = (0.33, 0.50 mm) pixel sizes with electronic settings of (a) 300 microseconds time constant, (b) normal dynamic reserve, (c) 24 dB/oct low pass filter slope, and (d) 160 Hz chopping frequency resulted in the highest SNR while maintaining sub-millimeter spatial resolution. The BaFBr_0.85 I_0.15 :Eu²⁺ storage phosphor dosimeters were linear from 1 to 50 MU and their profiles did not saturate up to 150 MU. The scanner was able to detect lateral displacements of ±1 mm in both the crossline and inline directions and ±0.3 mm beam spread changes that were artificially introduced by varying the incident proton energy. Specific to our proton unit, proton energy changes of ±1 MeV can also be detected indirectly via beam spread measurements.

CONCLUSION

Our combined dosimetric system including an in-house laboratory optical scanner and reusable BaFBr_0.85 I_0.15 :Eu²⁺ storage phosphors demonstrated a sufficient spatial resolution and dosimetric accuracy to support its use as an independent proton spot measurement dosimeter system. Its wide dynamic range allows for other versatile applications such as proton halo measurements.

Future technological developments in proton therapy - A predicted technological breakthrough

In the current spectrum of cancer treatments, despite high costs, a lack of robust evidence based on clinical outcomes or technical and radiobiological uncertainties, particle therapy and in particular proton therapy (PT) is rapidly growing. Despite proton therapy being more than fifty years old (first proposed by Wilson in 1946) and more than 220,000 patients having been treated with in 2020, many technological challenges remain and numerous new technical developments that must be integrated into existing systems. This article presents an overview of on-going technical developments and innovations that we felt were most important today, as well as those that have the potential to significantly shape the future of proton therapy. Indeed, efforts have been done continuously to improve the efficiency of a PT system, in terms of cost, technology and delivery technics, and a number of different developments pursued in the accelerator field will first be presented. Significant developments are also underway in terms of transport and spatial resolution achievable with pencil beam scanning, or conformation of the dose to the target: we will therefore discuss beam focusing and collimation issues which are important parameters for the development of these techniques, as well as proton arc therapy. State of the art and alternative approaches to adaptive PT and the future of adaptive PT will finally be reviewed. Through these overviews, we will finally see how advances in these different areas will allow the potential for robust dose shaping in proton therapy to be maximised, probably foreshadowing a future era of maturity for the PT technique.

Simulation study of proton arc therapy with the compact single-room MEVION-S250 proton therapy system

Aim:

The purpose of this study is to investigate the feasibility of proton arc therapy (PAT) using the double-scattering MEVION-S250 proton system. The treatment planning and dose delivery parameters from PAT were compared with conventional treatment planning techniques.

Materials and methods:

PAT was simulated with multiple conformal and fixed-aperture beams (5–15) using the MEVION-S250-double-scattering proton system. Conformal apertures were simulated with the Eclipse-treatment-planning system: (a) using a static single aperture that provides the best average conformal circular or rectangular apertures to cover the tumour from different angular views (SPAT), and (b) dynamic conformal apertures of the tumour shape at each irradiation angle that can be obtained from a multi-leaf-collimator system (DPAT).

Results:

The DPAT and SPAT plans provided superior dose coverage and sparing of normal tissues in comparison with conventional plans (CPT). The entrance normal tissue and skin doses (<40%) were lowered significantly by delivering dose from different directions over a wider angular view compared to conventional plans that have large entrance dose from only two fields. While the mean and minimum doses from PAT and CPT were comparable, the maximum doses from arc plans were lower than the maximum doses in conventional plans. The SPAT and DPAT plans had comparable dose parameters for regularly shaped targets. The heterogeneity index (HI) was superior in PAT plans which improved with increasing number of beams in arc plans for the different treatment sites. The conformality index (CI) depends on the treatment site and complexity of the shape of the planning target volume where for brain, pancreatic and lung tumours, PAT plans conformality was comparable and sometimes superior to CPT; and HI and CI were generally better in DPAT compared to SPAT.

Conclusions:

PAT plans have superior dose coverage and sparing of normal tissues compared to CPT plans using the MEVION double-scattering system as shown in this simulation study. Ideally, conformal proton arcs require beam shaping and dose delivery with the gantry moving; however, the MEVION double-scattering system lacks a multi-leaf collimator system and cannot deliver dose during gantry rotation. The single aperture conformal proton therapy technique is more time and cost effective compared with conventional techniques that are used currently with the MEVION proton therapy system because of the elimination of the need for patient-specific compensators. In present study, PAT was simulated with the MEVION double-scattering proton therapy system; however, it can be performed also with other proton therapy systems.

A Monte Carlo based analytic model of the in-room neutron ambient dose equivalent for a Mevion gantry-mounted passively scattered proton system

The goal of this study was to develop a Monte Carlo (MC)-based analytical model that can predict the in-room ambient dose equivalent from a Mevion gantry-mounted passively scattered proton system. The Mevion S250 and treatment vault were simulated using the MCNPX MC code. The results of the in-room neutron dose measurements, using an FHT 762 WENDI-II detector, were employed to benchmark the MC-derived values. After tuning the MCNPX MC code, for the same beam delivery parameters, the code was used to calculate the neutron spectra and ambient dose equivalent in the vault and at varying angles from the isocenter. Then, based on the calculations, an analytical model was reconstructed and data were fitted to derive the model parameters at 95% confidence intervals (CI). The MCNPX codes were tuned to within about 19% of the measured values for most of the measurements in the vault. For the maze, up to 0.08 mSv Gy-1 discrepancies were found between the experimental measurements and MCNPX calculated results. The analytical model showed up to 18% discrepancy for distances between 100 and 600 cm from the isocenter compared to the MC calculations. The model may underestimate the neutron ambient dose equivalent up to 21% for distances less than 100 cm from the isocenter. The proposed analytical model can be used to estimate the contribution of the secondary neutron dose from the Mevion S250 for the design of local shielding inside the proton therapy treatment vault.

Particle therapy in the future of precision therapy

The first hospital-based treatment facilities for particle therapy started operation about thirty years ago. Since then, the clinical experience with protons and carbon ions has grown continuously and more than 200,000 patients have been treated to date. The promising clinical results led to a rapidly increasing number of treatment facilities and many new facilities are planned or under construction all over the world. An inverted depth-dose profile combined with potential radiobiological advantages make charged particles a precious tool for the treatment of tumours that are particularly radioresistant or located nearby sensitive structures. A rising number of trials have already confirmed the benefits of particle therapy in selected clinical situations and further improvements in beam delivery, image guidance and treatment planning are expected. This review summarises some physical and biological characteristics of accelerated charged particles and gives some examples of their clinical application. Furthermore, challenges and future perspectives of particle therapy will be discussed.

Commissioning and beam characterization of the first gantry-mounted accelerator pencil beam scanning proton system

PURPOSE

To present and discuss beam characteristics and commissioning process of the first gantry-mounted accelerator single room pencil beam scanning (PBS) proton system.

METHODS

The Mevion HYPERSCAN employs a design configuration with a synchrocyclotron mounted on the gantry to eliminate the traditional beamline and a nozzle that contains the dosimetry monitoring chambers, the energy modulator (Energy Selector (ES)), and an Adaptive Aperture (AA). To characterize the beam, we measured the integrated depth dose (IDDs) for 12 energies, from highest energy of 227 MeV down to 28 MeV with a range difference ~ 2 cm between the adjacent energies, using a large radius Bragg peak chamber; single-spot profiles in air at five locations along the beam central axis using radiochromic EBT3 film and cross compared with a scintillation detector; and determined the output using a densely packed spot map. To access the performance of AA, we measured interleaf leakage and the penumbra reduction effect. Monte Carlo simulation using TOPAS was performed to study spot size variation along the beam path, beam divergence, and energy spectrum.

RESULTS

This proton system is calibrated to deliver 1 Gy dose at 5 cm depth in water using the highest beam energy by delivering 1 MU/spot to a 10 × 10 cm² map with a 2.5 mm spot spacing. The spot size in air varies from 4 mm to 26 mm from 227 MeV to 28 MeV at the isocenter plane with the nozzle retracted 23.6 cm from isocenter. The beam divergence of 28 MeV beam is ~ 52.7 mrad, which is nearly 22 times that of 227 MeV proton beam. The binary design of the ES has resulted in shifts of the effective SSD toward the isocenter as the energy is modulated lower. The peaks of IDD curves have a constant 80%-80% width of 8.4 mm at all energies. The interleaf leakage of the AA is less than 1.5% at the highest energy; and the AA can reduce the penumbra by 2 mm to 13 mm for the 227 and 28 MeV energies at isocenter plane in air.

CONCLUSIONS

The unique design of the HYPERSCAN proton system has yielded beam characteristics significantly different from that of other proton systems in terms of the Bragg peak shapes, spot sizes, and the penumbra sharpening effect of the AA. The combination of the ES and AA has made PBS implementation possible without using beam transport line and range shifter devices. Different considerations may be required in treatment planning optimization to account for different design and beam characteristics.

A Monte Carlo-based analytic model of neutron dose equivalent for a mevion gantry-mounted passively scattered proton system for craniospinal irradiation

PURPOSE

To calculate in- and out-of-field neutron spectra and dose equivalent, using Monte Carlo (MC) simulation, for a Mevion gantry-mounted passively scattered proton system in craniospinal irradiation. An analytical model based on the MC calculations that estimates in- and out-of-field neutron dose equivalent from proton Craniospinal irradiation (CSI) was also developed.

METHODS

The MCNPX MC code was used to simulate a Mevion S250 proton therapy system. The simulated proton depth doses and profiles for pristine and spread-out Bragg peaks were benchmarked against the measured data. Previous measurements using extended-range Bonner spheres were used to verify the calculated neutron spectra and dose equivalent. Using the benchmarked results as a reference condition, a correction-based analytical model was reconstructed by fitting the data to derive model parameters at 95% confidence interval. Sensitivity analysis of brass aperture opening, thickness of the Lucite (PMMA) range compensator, and modulation width was performed to obtain correction parameters for nonreference conditions.

RESULTS

For the neutron dose equivalent per therapeutic proton dose, the MCNPX calculated dose equivalent matched the measured values to within 8%. The benchmarked neutron dose equivalent at the isocenter was 41.2 and 20.8 mSv/Gy, for cranial and spinal fields, respectively. For in- and out-of-field neutron dose calculations, the correction-based analytical model showed up to 17% discrepancy compared to the MC calculations. The correction factors may provide a conservative estimation of neutron dose, especially for depth ≤ 5 cm and regions underneath the brass aperture.

CONCLUSION

The proposed analytical model can be used to estimate the contribution of the neutron dose to the overall CSI treatment dose. Moreover, the model can be employed to estimate the neutron dose to the implantable cardiac electronic devices.

Sensitivity analysis of Monte Carlo model of a gantry-mounted passively scattered proton system

Purpose

This study aimed to present guidance on the correlation between treatment nozzle and proton source parameters, and dose distribution of a passive double scattering compact proton therapy unit, known as Mevion S250.

Methods

All 24 beam options were modeled using the MCNPX MC code. The calculated physical dose for pristine peak, profiles, and spread out Bragg peak (SOBP) were benchmarked with the measured data. Track‐averaged LET (LET_t) and dose‐averaged LET (LET_d) distributions were also calculated. For the sensitivity investigations, proton beam line parameters including Average Energy (AE), Energy Spread (ES), Spot Size (SS), Beam Angle (BA), Beam Offset (OA), and Second scatter Offset (SO) from central Axis, and also First Scatter (FS) thickness were simulated in different stages to obtain the uncertainty of the derived results on the physical dose and LET distribution in a water phantom.

Results

For the physical dose distribution, the MCNPX MC model matched measurements data for all the options to within 2 mm and 2% criterion. The Mevion S250 was found to have a LET_t between 0.46 and 8.76 keV.μm^–1 and a corresponding LET_d between 0.84 and 15.91 keV.μm^–1. For all the options, the AE and ES had the greatest effect on the resulting depth of pristine peak and peak‐to‐plateau ratio respectively. BA, OA, and SO significantly decreased the flatness and symmetry of the profiles. The LETs were found to be sensitive to the AE, ES, and SS, especially in the peak region.

Conclusions

This study revealed the importance of considering detailed beam parameters, and identifying those that resulted in large effects on the physical dose distribution and LETs for a compact proton therapy machine.

Beam commissioning of the first compact proton therapy system with spot scanning and dynamic field collimation

OBJECTIVES

To describe the measurements and to present the results of the beam commissioning and the beam model validation of a compact, gantry-mounted, spot scanning proton accelerator system with dynamic layer-by-layer field collimation.

METHODS

We performed measurements of depth dose distributions in water, spot and scanned field size in air at different positions from the isocenter plane, spot position over the 20 × 20 cm² scanned area, beam monitor calibration in terms of absorbed dose to water and specific field collimation measurements at different gantry angles to commission the system. To validate the beam model in the treatment planning system (TPS), we measured spot profiles in water at different depths, absolute dose in water of single energy layers of different field sizes and inversely optimised spread-out Bragg peaks (SOBP) under normal and oblique beam incidence, field size and penumbra in water of SOBPs, and patient treatment specific quality assurance in homogeneous and heterogeneous phantoms.

RESULTS

Energy range, spot size, spot position and dose output were consistent at all gantry angles with 0.3 mm, 0.4 mm, 0.6 mm and 0.5% maximum deviations, respectively. Uncollimated spot size (one sigma) in air with an air-gap of 10 cm ranged from 4.1 to 16.4 mm covering a range from 32.2 to 1.9 cm in water, respectively. Absolute dose measurements were within 3% when comparing TPS and experimental data. Gamma pass rates >98% and >96% at 3%/3 mm were obtained when performing 2D dose measurements in homogeneous and in heterogeneous media, respectively. Leaf position was within ±1 mm at all gantry angles and nozzle positions.

CONCLUSIONS

Beam characterisation and machine commissioning results, and the exhaustive end-to-end tests performed to assess the proper functionality of the system, confirm that it is safe and accurate to treat patients.

ADVANCES IN KNOWLEDGE

This is the first paper addressing the beam commissioning and the beam validation of a compact, gantry-mounted, pencil beam scanning proton accelerator system with dynamic layer-by-layer multileaf collimation.

A novel approach to Verify air gap and SSD for proton radiotherapy using surface imaging

PURPOSE

To develop a novel approach to accurately verify patient set up in proton radiotherapy, especially for the verification of the nozzle - body surface air gap and source-to-skin distance (SSD), the consistency and accuracy of which is extremely important in proton treatment.

METHODS

Patient body surfaces can be captured and monitored with the optical surface imaging system during radiation treatment for improved intrafraction accuracy. An in-house software package was developed to reconstruct the patient body surface in the treatment position from the optical surface imaging reference capture and to calculate the corresponding nozzle - body surface air gap and SSD. To validate this method, a mannequin was scanned on a CT simulator and proton plans were generated for a Mevion S250 Proton machine with 20 gantry/couch angle combinations, as well as two different snout sizes, in the Varian Eclipse Treatment Planning Systems (TPS). The surface generated in the TPS from the CT scan was imported into the optical imaging system as an RT Structure for the purpose of validating and establishing a benchmark for ground truth comparison. The optical imaging surface reference capture was acquired at the treatment setup position after orthogonal kV imaging to confirm the positioning. The air gaps and SSDs calculated with the developed method from the surface captured at the treatment setup position (VRT surface) and the CT based surface imported from the TPS were compared to those calculated in TPS. The same approach was also applied to 14 clinical treatment fields for 10 patients to further validate the methodology.

RESULTS

The air gaps and SSDs calculated from our program agreed well with the corresponding values derived from the TPS. For the phantom results, using the CT surface, the absolute differences in the air gap were 0.45 mm ± 0.33 mm for the small snout, and 0.51 mm ± 0.49 mm for the large snout, and the absolute differences in SSD were 0.68 mm ± 0.42 mm regardless of snout size. Using the VRT surface, the absolute differences in air gap were 1.17 mm ± 1.17 mm and 2.1 mm ± 3.09 mm for the small and large snouts, respectively, and the absolute differences in SSD were 0.81 mm ± 0.45 mm. Similarly, for patient data, using the CT surface, the absolute differences in air gap were 0.42 mm ± 0.49 mm, and the absolute differences in SSD were 1.92 mm ± 1.4 mm. Using the VRT surface, the absolute differences in the air gap were 2.35 mm ± 2.3 mm, and the absolute differences in SSD were 2.7 mm ± 2.17 mm.

CONCLUSION

These results showed the feasibility and robustness of using an optical surface imaging approach to conveniently determine the air gap and SSD in proton treatment, providing an accurate and efficient way to confirm the target depth at treatment.

A dosimetric analysis of proton beam therapy using different snouts

Purpose

This exploration is intended to analyse the dosimetric characteristics of proton beams of multiple energies using different snout sizes.

Materials and methods

A synchrotron was used for the extraction of eight proton beam energies (100–250 MeV). Dosimetric measurements were taken in a water phantom that was irradiated with a proton beam emanating from the gantry system at angles 0, 90, 180 and 270 degree using a large and a medium snout. The range of beam energies in the phantom, their corresponding centre modulation depth (CMD) and the width of spread out Bragg peak (SOBP) were measured by Markus chamber. Double scattering technique was employed for the creation of SOBPs.

Results

The range of proton beams varied from 4·3 cm for 100 MeV beam to 28·5 cm for 250 MeV beam with the medium snout and from 4·3 cm for 100 MeV to 25 cm for 250 MeV beam with large snout in the water phantom. SOBP width showed a variation from 4 to 10 cm with medium and large snout. While determining the output with medium snout, the discrepancy of 1·1% was observed between the maximum and minimum mean values of output for all the given set of energies and angles. There occurred a difference of 0·9% between the maximum and minimum mean values of output with the large snout. Beam output at SOBP centre was 12% higher with large snout as compared to that with medium snout for all the given beam energies. Flatness and symmetry were found within ±2·5% tolerance limits with medium and large snouts.

Conclusion

Flatness and symmetry were found within explicit limits with both medium and large snouts. Large snout produced higher beam output than that of medium snout at the centre of SOBP. This exploration can be extended to the determination of beam output, flatness and symmetry with a small snout.

Commissioning of a fluoroscopic-based real-time markerless tumor tracking system in a superconducting rotating gantry for carbon-ion pencil beam scanning treatment

PURPOSE

To perform the final quality assurance of our fluoroscopic-based markerless tumor tracking for gated carbon-ion pencil beam scanning (C-PBS) radiotherapy using a rotating gantry system, we evaluated the geometrical accuracy and tumor tracking accuracy using a moving chest phantom with simulated respiration.

METHODS

The positions of the dynamic flat panel detector (DFPD) and x-ray tube are subject to changes due to gantry sag. To compensate for this, we generated a geometrical calibration table (gantry flex map) in 15° gantry angle steps by the bundle adjustment method. We evaluated five metrics: (a) Geometrical calibration was evaluated by calculating chest phantom positional error using 2D/3D registration software for each 5° step of the gantry angle. (b) Moving phantom displacement accuracy was measured (±10 mm in 1-mm steps) with a laser sensor. (c) Tracking accuracy was evaluated with machine learning (ML) and multi-template matching (MTM) algorithms, which used fluoroscopic images and digitally reconstructed radiographic (DRR) images as training data. The chest phantom was continuously moved ±10 mm in a sinusoidal path with a moving cycle of 4 s and respiration was simulated with ±5 mm expansion/contraction with a cycle of 2 s. This was performed with the gantry angle set at 0°, 45°, 120°, and 240°. (d) Four types of interlock function were evaluated: tumor velocity, DFPD image brightness variation, tracking anomaly detection, and tracking positional inconsistency in between the two corresponding rays. (e) Gate on/off latency, gating control system latency, and beam irradiation latency were measured using a laser sensor and an oscilloscope.

RESULTS

By applying the gantry flex map, phantom positional accuracy was improved from 1.03 mm/0.33° to <0.45 mm/0.27° for all gantry angles. The moving phantom displacement error was 0.1 mm. Due to long computation time, the tracking accuracy achieved with ML was <0.49 mm (=95% confidence interval [CI]) for imaging rates of 15 and 7.5 fps; those at 30 fps were decreased to 1.84 mm (95% CI: 1.79 mm-1.92 mm). The tracking positional accuracy with MTM was <0.52 mm (=95% CI) for all gantry angles and imaging frame rates. The tumor velocity interlock signal delay time was 44.7 ms (=1.3 frame). DFPD image brightness interlock latency was 34 ms (=1.0 frame). The tracking positional error was improved from 2.27 ± 2.67 mm to 0.25 ± 0.24 mm by the tracking anomaly detection interlock function. Tracking positional inconsistency interlock signal was output within 5.0 ms. The gate on/off latency was <82.7 ± 7.6 ms. The gating control system latency was <3.1 ± 1.0 ms. The beam irradiation latency was <8.7 ± 1.2 ms.

CONCLUSIONS

Our markerless tracking system is now ready for clinical use. We hope to shorten the computation time needed by the ML algorithm at 30 fps in the future.

Experimental Assessment of Proton Dose Calculation Accuracy in Small-Field Delivery Using a Mevion S250 Proton Therapy System

Purpose:

Dose calculation accuracy of the Varian Eclipse treatment planning system (TPS) is empirically assessed for small-aperture fields using a Mevion S250 double scattering proton therapy system.

Materials and Methods:

Five spherical pseudotumors were modeled in a RANDO head phantom. Plans were generated for the targets with apertures of 1, 2, 3, 4, or 5 cm diameter using one, two, and three beams. Depth-dose curves and lateral profiles of the beams were taken with the planned blocks and compared to Eclipse calculations. Dose distributions measured with EBT3 films in the phantom were also compared to Eclipse calculations. Film quenching effect was simulated and considered.

Results:

Depth-dose scans in water showed a range pullback (up to 2.0 mm), modulation widening (up to 9.5 mm), and dose escalation in proximal end and sub-peak region (up to 15.5%) when compared to the Eclipse calculations for small fields. Measured full width at half maximum and penumbrae for lateral profiles differed <1.0 mm from calculations for most comparisons. In the phantom study, Eclipse TPS underestimated sub-peak dose. Gamma passing rates improved with each beam added to the plans. Greater range pullback and modulation degradation versus water scans were observed due to film quenching, which became more noticeable as target size increased.

Conclusions:

Eclipse TPS generates acceptable target coverage for small targets with carefully arranged multiple beams despite relatively large dose discrepancy for each beam. Surface doses higher than Eclipse calculations can be mitigated with multiple beams. When using EBT3 film, the quenching effect should be considered.

Image guidance in proton therapy for lung cancer

Proton therapy is a promising but challenging treatment modality for the management of lung cancer. The technical challenges are due to respiratory motion, low dose tolerance of adjacent normal tissue and tissue density heterogeneity. Different imaging modalities are applied at various steps of lung proton therapy to provide information on target definition, target motion, proton range, patient setup and treatment outcome assessment. Imaging data is used to guide treatment design, treatment delivery, and treatment adaptation to ensure the treatment goal is achieved. This review article will summarize and compare various imaging techniques that can be used in every step of lung proton therapy to address these challenges.

The Mobius AIRO mobile CT for image-guided proton therapy: Characterization & commissioning

PURPOSE

The purpose of this study was to characterize the Mobius AIRO Mobile CT System for localization and image-guided proton therapy. This is the first known application of the AIRO for proton therapy.

METHODS

Five CT images of a Catphan® 504 phantom were acquired on the AIRO Mobile CT System, Varian EDGE radiosurgery system cone beam CT (CBCT), Philips Brilliance Big Bore 16 slice CT simulator, and Siemens SOMATOM Definition AS 20 slice CT simulator. DoseLAB software v.6.6 was utilized for image quality analysis. Modulation transfer function, scaling discrepancy, geometric distortion, spatial resolution, overall uniformity, minimum uniformity, contrast, high CNR, and maximum HU deviation were acquired. Low CNR was acquired manually using the CTP515 module. Localization accuracy and CT Dose Index were measured and compared to reported values on each imaging device. For treatment delivery systems (Edge and Mevion), the localization accuracy of the 3D imaging systems were compared to 2D imaging systems on each system.

RESULTS

The AIRO spatial resolution was 0.21 lp mm-1 compared with 0.40 lp mm-1 for the Philips CT Simulator, 0.37 lp mm-1 for the Edge CBCT, and 0.35 lp mm-1 for the Siemens CT Simulator. AIRO/Siemens and AIRO/Philips differences exceeded 100% for scaling discrepancy (191.2% and 145.8%). The AIRO exhibited higher dose (>27 mGy) than the Philips CT Simulator. Localization accuracy (based on the MIMI phantom) was 0.6° and 0.5 mm. Localization accuracy (based on Stereophan) demonstrated maximum AIRO-kV/kV shift differences of 0.1 mm in the x-direction, 0.1 mm in the y-direction, and 0.2 mm in the z-direction.

CONCLUSIONS

The localization accuracy of AIRO was determined to be within 0.6° and 0.5 mm despite its slightly lower image quality overall compared to other CT imaging systems at our institution. Based on our study, the Mobile AIRO CT system can be utilized accurately and reliably for image-guided proton therapy.

TOPAS Simulation of the Mevion S250 compact proton therapy unit

As proton therapy becomes increasingly popular, so does the need for Monte Carlo simulation studies involving accurate beam line modeling of proton treatment units. In this study, the 24 beam configurations of the Mevion S250 proton therapy system installed recently at our institution were modeled using the TOolkit for PArticle Simulation (TOPAS) code. Pristine Bragg peak, spread out Bragg peak (SOBP), and lateral beam profile dose distributions were simulated and matched to the measurements taken during commissioning of the unit. Differences in the range for all Percent Depth Dose (PDD) curves between measured and simulated data agreed to within 0.1 cm. For SOBP scans, the SOBP widths all agreed to within 0.3 cm. With regards to lateral beam profile comparisons between the measured and simulated data, the penumbras differed by less than 1 mm and the flatness differed by less than 1% in nearly all cases. This study shows that Monte Carlo simulation studies involving the Mevion S250 proton therapy unit can be a viable tool in commissioning and verification of the proton treatment planning system.

Comparability of three output prediction models for a compact passively double-scattered proton therapy system

The purpose of this study was to investigate comparability of three output prediction models for a compact double‐scattered proton therapy system. Two published output prediction models are commissioned for our Mevion S250 proton therapy system. Model A is a correction‐based model (Sahoo et al., Med Phys, 2008;35(11):5088–5097) and model B is an analytical model which employs a function of r = (R’‐M’)/M’ (Kooy et al., Phys Med Biol, 2005;50:5487–5456) where R’ is defined as depth of distal 100% dose with straggling and M’ is the width between distal 100% dose and proximal 100% dose with straggling instead of the theoretical definition due to more accurate output prediction. The r is converted to ((R‐0.31)‐0.81 × M)/(0.81 × M) with the vendor definition of R (distal 90% dose) and M (distal 90% dose‐to‐proximal 95% dose), where R’ = R‐0.31 (g cm⁻²) and M’ = 0.81 × M (g cm⁻²). In addition, a quartic polynomial fit model (model C) mathematically converted from model B is studied. The outputs of 272 sets of R and M covering the 24 double scattering options are measured. Each model's predicted output is compared to the measured output. For the total dataset, the percent difference between predicted (P) and measured (M) outputs ((P‐M)/M × 100%) were within ±3% using the three different models. The average differences (±standard deviation) were −0.13 ± 0.94%, −0.13 ± 1.20%, and −0.22 ± 1.11% for models A, B, and C, respectively. The p‐values of the t‐test were 0.912 (model A vs. B), 0.061 (model A vs. C), and 0.136 (model B vs. C). For all the options, all three models have clinically acceptable predictions. The differences between models A, B, and C are statistically insignificant; however, model A generally has the potential to more accurately predict the output if a larger dataset for commissioning is used. It is concluded that the models can be comparably used for the compact proton therapy system.

[Technical aspects of protontherapy: Setup, equipment and radioprotection]

The number of protontherapy facilities is still increasing rapidly with more than 30 ongoing projects and close to 60 currently under operation. Although the technology is now validated and robust, a proton facility cannot be considered as a standard radiation therapy equipment: its constraints in terms of building, services, project management are of paramount impact at the level of the hospital. Therefore, a protontherapy project must be carefully considered and prepared, which is mandatory for further fluid and efficient clinical operation.