Background and purpose Previous studies have demonstrated that mini-ridge filters (MRF) can reduce beam delivery time in pencil beam scanning (PBS) proton therapy, but the effect on plan quality has not been analyzed in detail. Furthermore, the range modulation produced by the MRF may lead to an added benefit of reduced dose-averaged linear energy transfer (LETd) on the distal edge of the beam. The aim of this study is to demonstrate a statistically significant reduction in beam delivery time and distal LETd with a MRF without any degradation in plan quality. Materials and methods Twenty patient plans from a range of treatment sites, including five deep inspiration breath hold (DIBH) patients, have been re-optimized with a Monte Carlo based MRF beam model and compared to the corresponding clinical plans in beam delivery time, clinical acceptability, and a range of dose and LETd metrics. Results The average delivery time per beam decreased from (96.7 ± 4.9) s to (75.3 ± 4.0) s (p < 0.001). DIBH patients saw a reduction from (82.8 ± 15.5) s to (70.4 ± 12.7) s. The average reduction in maximum LETd to any organ-at-risk (OAR) within the 500 cGyRBE isodose line was (0.43 ± 0.06) keV/μm (p < 0.001). There was no significant difference in target dose homogeneity, conformality, robust target coverage, OAR dose metrics or LETd within the target. Conclusions MRFs can be used in PBS proton therapy to achieve a clinically relevant reduction in beam delivery time, and distal LETd without sacrificing plan quality.

BACKGROUND:Neurosurgical cranial titanium mesh and screws are commonly encountered in postoperative radiation therapy. However, only a limited number of reports are available in the context of proton therapy, resulting in a lack of consensus among the proton centers regarding the protocol for handling the hardware. PURPOSE:This study is to examine the impact of the hardware in proton plans. The results serve as evidence for proton centers to generate standard operating procedures to manage the hardware in proton treatment. METHODS:Plans with different gantry angles and material overrides are generated on the CT images of a phantom made of the hardware. The dose distributions of the plans with and without material override, at different depths are compared. Films and ionization chambers are used to measure the plans and the measurements are compared to the treatment planning system (TPS) calculations by gamma analysis. RESULTS:There are some overdose and underdose regions downstream of the hardware. The overdose and underdose values are within a few percent of the prescribed dose when multiple fields with large hinge angles are used. The gamma analysis results show that the measurements agree with the TPS calculations within limits that are clinically relevant. CONCLUSION:The study has demonstrated the influence of the hardware on proton plans. Based on the result of this study, a standard operating procedure of managing the hardware has been implemented in our clinic.

Background Until today, the majority of ocular proton treatments worldwide were planned with the EYEPLAN treatment planning system (TPS). Recently, the commercial, computed tomography (CT)-based TPS for ocular proton therapy RayOcular was released, which follows the general concepts of model-based treatment planning approach in conjunction with a pencil-beam-type dose algorithm (PBA). Purpose To validate RayOcular with respect to two main features: accurate geometrical representation of the eye model and accuracy of its dose calculation algorithm in combination with an Ion Beam Applications (IBA) eye treatment delivery system. Methods Different 3D-printed eye-ball-phantoms were fabricated to test the geometrical representation of the corresponding CT-based model, both in orthogonal 2D images for X-ray image overlay and in fundus view overlaid with a funduscopy. For the latter, the phantom was equipped with a lens matching refraction of the human eye. Funduscopy was acquired in a Zeiss Claus 500 camera. Tantalum clips and fiducials attached to the phantoms were localized in the TPS model, and residual deviations to the actual position in X-ray images for various orientations of the phantom were determined, after the nominal eye orientation was corrected in RayOcular to obtain a best overall fit. In the fundus view, deviations between known and displayed distances were measured. Dose calculation accuracy of the PBA on a 0.2 mm grid was investigated by comparing between measured lateral and depth-dose profiles in water for various combinations of range, modulation, and field-size. Ultimately, the modeling of dose distributions behind wedges was tested. A 1D gamma-test was applied, and the lateral and distal penumbra were further compared. Results Average residuals between model clips and visible clips/fiducials in orthogonal X-ray images were within 0.3 mm, including different orientations of the phantom. The differences between measured distances on the registered funduscopy image in the RayOcular fundus view and the known ground-truth were within 1 mm up to 10.5 mm distance from the posterior pole. No clear benefit projection of either polar mode or camera mode could be identified, the latter mimicking camera properties. Measured dose distributions were reproduced with gamma-test pass-rates of >95% with 2%/0.3 mm for depth and lateral profiles in the middle of spread-out Bragg-peaks. Distal falloff and lateral penumbra were within 0.2 mm for fields without a wedge. For shallow depths, the agreement was worse, reaching pass-rates down to 80% with 5%/0.3 mm when comparing lateral profiles in air. This is caused by low-energy protons from a scatter source in the IBA system not modeled by RayOcular. Dose distributions modified by wedges were reproduced, matching the wedge-induced broadening of the lateral penumbra to within 0.4 mm for the investigated cases and showing the excess dose within the field due to wedge scatter. Conclusion RayOcular was validated for its use with an IBA single scattering delivery nozzle. Geometric modeling of the eye and representation of 2D projections fulfill clinical requirements. The PBA dose calculation reproduces measured distributions and allows explicit handling of wedges, overcoming approximations of simpler dose calculation algorithms used in other systems.

Purpose: The primary fluence of a proton pencil beam exiting the accelerator is enveloped by a region of secondaries, commonly called "spray". Although small in magnitude, this spray may affect dose distributions in pencil beam scanning mode e.g., in the calculation of the small field output, if not modelled properly in a treatment planning system (TPS). The purpose of this study was to dosimetrically benchmark the Monte Carlo (MC) dose engine of the RayStation TPS (v.10A) in small proton fields and systematically compare single Gaussian (SG) and double Gaussian (DG) modeling of initial proton fluence providing a more accurate representation of the nozzle spray.Methods: The initial proton fluence distribution for SG/DG beam modeling was deduced from two-dimensional measurements in air with a scintillation screen with electronic readout. The DG model was either based on direct fits of the two Gaussians to the measured profiles, or by an iterative optimization procedure, which uses the measured profiles to mimic in-air scan-field factor (SF) measurements. To validate the DG beam models SFs, i.e. relative doses to a 10 x 10 cm(2) field, were measured in water for three different initial proton energies (100MeV, 160MeV, 226.7MeV) and square field sizes from 1x1cm(2) to 10x10cm(2) using a small field ionization chamber (IBA CC01) and an IBA ProteusPlus system (universal nozzle). Furthermore, the dose to the center of spherical target volumes (diameters: 1cm to 10cm) was determined using the same small volume ionization chamber (IC). A comprehensive uncertainty analysis was performed, including estimates of influence factors typical for small field dosimetry deduced from a simple two-dimensional analytical model of the relative fluence distribution. Measurements were compared to the predictions of the RayStation TPS.Results: SFs deviated by more than 2% from TPS predictions in all fields <4x4cm(2) with a maximum deviation of 5.8% for SG modeling. In contrast, deviations were smaller than 2% for all field-sizes and proton energies when using the directly fitted DG model. The optimized DG model performed similarly except for slightly larger deviations in the 1x1cm(2) scan-fields. The uncertainty estimates showed a significant impact of pencil beam size variations (+/- 5%) resulting in up to 5.0% standard uncertainty. The point doses within spherical irradiation volumes deviated from calculations by up to 3.3% for the SG model and 2.0% for the DG model.Conclusion: Properly representing nozzle spray in RayStation's MC-based dose engine using a DG beam model was found to reduce the deviation to measurements in small spherical targets to below 2%. A thorough uncertainty analysis shows a similar magnitude for the combined standard uncertainty of such measurements.

BACKGROUND:Pencil beam scanning (PBS) proton therapy allows for far superior dose conformality compared with passive scattering techniques. However, one drawback of PBS is that the beam delivery time can be long, particularly when treating superficial disease. Minimizing beam delivery time is important for patient comfort and precision of treatment delivery. Mini-ridge filters (MRF) have been shown to reduce beam delivery time for synchrotron-based PBS. Given that cyclotron systems are widely used in proton therapy it is necessary to investigate the potential clinical benefit of mini-ridge filters in such systems. PURPOSE:To demonstrate the clinical benefit of using a MRF to reduce beam delivery time for patients with large target volumes and superficial disease in cyclotron-based PBS proton therapy. METHODS:A MRF beam model was generated by simulating the effect of a MRF on our clinical beam data assuming a fixed snout position relative to the isocenter. The beam model was validated with a series of measurements. The model was used to optimize treatment plans in a water phantom and on six patient DICOM datasets to further study the effect of the MRF and for comparison with physician-approved clinical treatment plans. Beam delivery time was measured for six plans with and without the MRF to demonstrate the reduction achievable. Plans with and without MRF were reviewed to confirm clinical acceptability by a radiation oncologist. Patient-specific QA measurements were carried out with a two-dimensional ionization chamber array detector for one representative patient's plan optimized with the MRF beam model. RESULTS:Results show good agreement between the simulated beam model and measurements with mean and maximum deviations of 0.06 mm (0.45%) and 0.61 mm (4.9%). The increase in Bragg peak width (FWHM) ranged from 2.7 mm at 226 MeV to 6.1 mm at 70 MeV. The mean and maximum reduction in beam delivery time observed per field was 29.1 s (32.2%) and 79.7 s (55.3%). CONCLUSION:MRFs can be used to reduce treatment time in cyclotron-based PBS proton therapy without sacrificing plan quality. This is particularly beneficial for patients with large targets and superficial disease such as in breast cancer where treatment times are generally long, as well as patients treated with deep inspiration breath hold (DIBH).

The use of scanned proton beams in external beam radiation therapy has seen a rapid development over the past decade. This technique places new demands on treatment planning, as compared to conventional photon-based radiation therapy. In this article, several proton specific functions as implemented in the treatment planning system RayStation are presented. We will cover algorithms for energy layer and spot selection, basic optimization including the handling of spot weight limits, optimization of the linear energy transfer (LET) distribution, robust optimization including the special case of 4D optimization, proton arc planning, and automatic planning using deep learning. We will further present the Monte Carlo (MC) proton dose engine in RayStation to some detail, from the material interpretation of the CT data, through the beam model parameterization, to the actual MC transport mechanism. Useful tools for plan evaluation, including robustness evaluation, and the versatile scripting interface are also described. The overall aim of the paper is to give an overview of some of the key proton planning functions in RayStation, with example usages, and at the same time provide the details about the underlying algorithms that previously have not been fully publicly available.

Purpose To evaluate the impact of beam quality in terms of distal fall-off (DFO, 90%-10%) and lateral penumbra (LP, 80%-20%) of single beam ocular proton therapy (OPT) and to derive resulting ideal requirements for future systems. Methods Nine different beam models with DFO varying between 1 and 4 mm and LP between 1 and 4 mm were created. Beam models were incorporated into the RayStation with RayOcular treatment planning system version 10 B (RaySearch Laboratories, Stockholm, Sweden). Each beam model was applied for eight typical clinical cases, covering different sizes and locations of uveal melanoma. Plans with and without an additional wedge were created, resulting in 117 plans with a total prescribed median dose of 60 Gy(RBE) to the clinical target volume. Treatment plans were analyzed in terms of V20-V80 penumbra volume, D1 (dose to 1% of the volume) for optic disc and macula, optic nerve V30 (volume receiving 30 Gy(RBE), i.e., 50% of prescription), as well as average dose to lens and ciliary body. An LP-dependent aperture margin was based on estimated uncertainties, ranging from 1.7 to 4.0 mm. Results V20-V80 showed a strong influence by LP, while DFO was less relevant. The optic disc D1 reached an extra dose of up to 3000 cGy(RBE), comparing the defined technical limit of DFO = LP = 1 mm with DFO = 3 mm/LP = 4 mm. The latter may result from a pencil-beam scanning (PBS) system with static apertures. Plans employing a wedge showed an improvement for organs at risk sparing. Conclusion Plan quality is strongly influenced by initial beam parameters. The impact of LP is more pronounced when compared to DFO. The latter becomes important in the treatment of posterior tumors near the macula, optic disc or optic nerve. The plan quality achieved by dedicated OPT nozzles in single- or double-scattering design might not be achievable with modified PBS systems.

PURPOSE:To perform the validation of the GPU-based (Graphical Processing Unit based) proton Monte Carlo (MC) dose engine implemented in a commercial TPS (RayStation 10B) and to report final dose calculation times for clinical cases. MATERIALS AND METHODS:440 patients treated at the Proton Therapy Center of Trento, Italy, between 2018 and 2019 were selected for this study. 636 approved plans with 3361 beams computed with the clinically implemented CPU-MC dose engine (version 4.2 and 4.5), were used for the validation of the new algorithm. For each beam, the dose was recalculated using the new GPU-MC dose engine with the initial CPU computation settings and compared to the original CPU-MC dose. Beam dose difference distributions were studied to ensure that the two dose distributions were equal within the expected fluctuations of the MC statistical uncertainty (s) of each computation. Plan dose distributions were compared with respect to the dosimetric indices D98, D50 and D1 of all ROIs defined as targets. A complete assessment of the computation time as a function of s and dose grid voxel size was done. RESULTS:The median over all mean beam dose differences between CPU- and GPU-MC was -0.01% and the median of the corresponding standard deviations was close to (√2s) both for simulations with an s of 0.5% and 1.0% per beam. This shows that the two dose distributions can be considered equal. All the DVH indices showed an average difference below 0.04%. About half of the plans were computed with 1.0% statistical uncertainty on a 2 mm dose calculation grid, for which the median computation time was 5.2 s. The median computational speed for all plans in the study was 8.4 million protons/second. CONCLUSION:A validation of a clinical MC algorithm running on GPU was performed on a large pool of patients treated with pencil beam scanning proton therapy. We demonstrated that the differences with the previous CPU-based MC were only due to the intrinsic statistical fluctuations of the MC method, which translated to insignificant differences on plan dose level. The significant increase in dose calculation speed is expected to facilitate new clinical workflows.

Background and purpose: Monte Carlo simulations as well as analytical computations of proton transport in material media require accurate values of multiple Coulomb scattering (MCS) angles. High-quality experimental data on MCS angles in the energy range for proton therapy are, however, sparse. In this work, MCS modeling in proton transport was evaluated employing an experimental method to measure these angles on a medical proton beamline in clinically relevant materials. Results are compared to Monte Carlo simulations and analytical models. Materials and methods: Aluminum, brass, and lucite (PMMA) scatterers of clinically relevant thicknesses were irradiated with protons at 100, 160, and 220 MeV. Resulting spatial distributions of individual pencil beams were measured with a scintillating screen. The MCS angles were determined by deconvolution and a virtual point source approach. Results were compared to those obtained with the Monte Carlo codes PENH, TOPAS, and RayStation Monte Carlo, as well as the analytical models RayStation Pencil Beam Algorithm and the Moliere/Fano/Hanson variant of the Moliere theory. Results: Experimental data obtained with the presented methodology agree with previously published results within 6%, with an average deviation of 3%. The combined average uncertainty of the experimental data yielded 1.8%, while the combined maximum uncertainty was below 4%. The obtained Monte Carlo results for PENH, TOPAS, and RayStation deviate on average for all considered energies, materials and thicknesses, by 2.5%, 3.4%, and 2.8% from the experimental data, respectively. For the analytical models, the average deviations amount to 4.5% and 2.9% for the RayStation Pencil Beam Algorithm and the MoliAre/Fano/Hanson model, respectively. Conclusion: The experimental method developed for the present work allowed to measure MCS angles in clinical proton facilities with good accuracy. The presented method permits to extend the database on experimental MCS angles which is rather limited. This work further provides benchmark data for lucite in thicknesses relevant for clinical applications. The data may serve to validate dose engines of treatment planning systems and secondary dose check software. The Monte Carlo and analytical algorithms studied are capable of reproducing MCS data within the required accuracy for clinical applications. (C) 2021 The Authors. Medical Physics published by Wiley Periodicals LLC on behalf of American Association of Physicists in Medicine

Purpose: This work aims to validate new 6D couch features and their implementation for seated radiotherapy in RayStation (RS) treatment planning system (TPS). Materials and methods: In RS TPS, new 6D couch features are (i) chair support device, (ii) patient treatment option of "Sitting: face towards the front of the chair", and (iii) patient support pitch and roll capabilities. The validation of pitch and roll was performed by comparing TPS generated DRRs with planar x-rays. Dosimetric tests through measurement by 2D ion chamber array were performed for beams created with varied scanning and treatment orientation and 6D couch rotations. For the implementation of 6D couch features for treatments in a seated position, the TPS and oncology information system (Mosaiq) settings are described for a commercial chair. An end-to-end test using an anthropomorphic phantom was performed to test the complete workflow from simulation to treatment delivery. Results: The 6D couch features were found to have a consistent implementation that met IEC 61712 standard. The DRRs were found to have an acceptable agreement with planar x-rays based on visual inspection. For dose map comparison between measured and calculated, the gamma index analysis for all the beams was >95% at a 3% dose-difference and 3 mm distance-to -agreement tolerances. For an end-to end-testing, the phantom was successfully set up at isocenter in the seated position and treatment was delivered. Conclusions: Chair-based treatments in a seated position can be implemented in RayStation through the use of newly released 6D couch features.

AbstractPurposeThe aim of this study is to validate the RayStation Monte Carlo (MC) dose algorithm using animal tissue neck phantoms and a water breast phantom.MethodsThree anthropomorphic phantoms were used in a clinical setting to test the RayStation MC dose algorithm. We used two real animal necks that were cut to a workable shape while frozen and then thawed before being CT scanned. Secondly, we made a patient breast phantom using a breast prosthesis filled with water and placed on a flat surface. Dose distributions in the animal and breast phantoms were measured using the MatriXX PT device.ResultsThe measured doses to the neck and breast phantoms compared exceptionally well with doses calculated by the analytical pencil beam (APB) and MC algorithms. The comparisons between APB and MC dose calculations and MatriXX PT measurements yielded an average depth difference for best gamma agreement of <1 mm for the neck phantoms. For the breast phantom better average gamma pass rates between measured and calculated dose distributions were observed for the MC than for the APB algorithms.ConclusionsThe MC dose calculations are more accurate than the APB calculations for the static phantoms conditions we evaluated, especially in areas where significant inhomogeneous interfaces are traversed by the beam.

Some clinical indications require small fields with sharp lateral dose gradients, which is technically challenging in proton beam therapy. This holds especially true for low-range fields applied with the spot scanning technique, where large beam profiles entering from the beam-line or the insertion of range shifting blocks lead to large lateral gradients. We regard the latter case and solve it by shifting the range shifting block far upstream in conjunction with a collimating aperture dose to the patient. The experiments of the current work are based on a commercial proton therapy treatment head designed for several delivery modes. In a research environment of the spot-scanning delivery mode a range shifter is inserted downstream of the scanning magnets in a slot which is usually employed only in a scattering delivery mode. This configuration is motivated by equations assuming a simple model of proton transport. In the experiments lateral dose planes are acquired with a scintillation screen and radiochromic films. Dose distributions are calculated with the Monte Carlo dose engine of the RayStation treatment planning system. We demonstrate that proton fields with 80%-20% lateral dose fall-off values between 1.4 mm and 4.0 mm can be achieved for water equivalent depths between 0 cm and 10 cm. The simulated lateral dose profiles agree with the experimental dose profiles. The sharpening of the field edges is set off by a broadening of the proton spots towards the center of the fields. This limits the clinical application mainly to small fields for which the distal and proximal conformality is of minor importance.

The objective of this study was to examine the use of proton pencil beam scanning for the treatment of moving lung tumors. A single-field uniform dose proton pencil beam scanning (PBS) plan was generated for the standard thorax phantom designed by the Imaging and Radiation Oncology Core (IROC) Houston QA Center. Robust optimization, including range and setup uncertainties as well as volumetric repainting, was used for the plan. Patient-specific quality assurance (QA) measurements were performed using both a water tank and a custom heterogeneous QA phantom. A custom moving phantom was used to find the optimal number of volumetric repainting. Both analytical and Monte Carlo (MC) algorithms were used for dose calculation and their accuracies were compared with actual measurements. A single ionization chamber, a 2-dimensional ionization chamber array, thermoluminescent dosimeters (TLDs), and films were used for dose measurements. The optimal number of volumetric repainting was found to be 4 times in our system. The mean dose overestimations on a moving target by analytical and MC algorithms based on a time-averaged computed tomography (CT) image of the phantom were found to be 4.8% and 2.4%, respectively. The mean gamma indexes for analytical and MC algorithms were 91% and 96%, respectively. The MC dose algorithm calculation was found to have a better agreement with measurements compared with the analytical algorithm. When treating moving lung tumors using proton PBS, the techniques of robust optimization, volumetric repainting, and MC dose calculation were found effective. Extra care needs to be taken when an analytical dose calculation algorithm is used.

RaySearch Americas Inc. (NY) has introduced a commercial Monte Carlo dose algorithm (RS-MC) for routine clinical use in proton spot scanning. In this report, we provide a validation of this algorithm against phantom measurements and simulations in the GATE software package. We also compared the performance of the RayStation analytical algorithm (RS-PBA) against the RS-MC algorithm. A beam model (G-MC) for a spot scanning gantry at our proton center was implemented in the GATE software package. The model was validated against measurements in a water phantom and was used for benchmarking the RS-MC. Validation of the RS-MC was performed in a water phantom by measuring depth doses and profiles for three spread-out Bragg peak (SOBP) beams with normal incidence, an SOBP with oblique incidence, and an SOBP with a range shifter and large air gap. The RS-MC was also validated against measurements and simulations in heterogeneous phantoms created by placing lung or bone slabs in a water phantom. Lateral dose profiles near the distal end of the beam were measured with a microDiamond detector and compared to the G-MC simulations, RS-MC and RS-PBA. Finally, the RS-MC and RS-PBA were validated against measured dose distributions in an Alderson-Rando (AR) phantom. Measurements were made using Gafchromic film in the AR phantom and compared to doses using the RS-PBA and RS-MC algorithms. For SOBP depth doses in a water phantom, all three algorithms matched the measurements to within  ±3% at all points and a range within 1 mm. The RS-PBA algorithm showed up to a 10% difference in dose at the entrance for the beam with a range shifter and  >30 cm air gap, while the RS-MC and G-MC were always within 3% of the measurement. For an oblique beam incident at 45°, the RS-PBA algorithm showed up to 6% local dose differences and broadening of distal fall-off by 5 mm. Both the RS-MC and G-MC accurately predicted the depth dose to within  ±3% and distal fall-off to within 2 mm. In an anthropomorphic phantom, the gamma index (dose tolerance  =  3%, distance-to-agreement  =  3 mm) was greater than 90% for six out of seven planes using the RS-MC, and three out seven for the RS-PBA. The RS-MC algorithm demonstrated improved dosimetric accuracy over the RS-PBA in the presence of homogenous, heterogeneous and anthropomorphic phantoms. The computation performance of the RS-MC was similar to the RS-PBA algorithm. For complex disease sites like breast, head and neck, and lung cancer, the RS-MC algorithm will provide significantly more accurate treatment planning.

OBJECTIVES:To report complications after colonoscopy and surgery in patients with neoplasia detected through a population based colorectal cancer (CRC) screening programme in the capital region of Sweden.METHODS:Patients who after a positive FOBT screening result underwent colonoscopy from 1 January 2008 to 30 June 2012 were included. Mortality and complications within 30 days after colonoscopy or subsequent surgery were identified through national registers, and complications were assessed through review of medical charts. Complications were graded using the Clavien-Dindo classification.RESULTS:After 2984 colonoscopies, the complication rate was 1%. The risk of post-polypectomy bleeding was 14/1000. The risk of perforation was 1/1000 after a diagnostic colonoscopy and 2.5/1000 after a colonoscopy with polypectomy. One patient developed a post-polypectomy syndrome. There was one death which was not related to the colonoscopy. After surgery for 37 adenomas and 155 CRCs, the total complication rates were 27% and 50%, respectively. The rate of anastomotic leakage was 13% and 12% after surgery for adenomas and CRC, respectively. There were no deaths after surgery. The overall complication rate after colonoscopy and surgery for adenomas and cancer was 4%.CONCLUSIONS:Overall complication rates were acceptable and mortality low; however, the rate of anastomotic leakage after surgery for both adenomas and CRC was higher than expected.