Monte Carlo investigation of collimator scatter of proton-therapy beams produced using the passive scattering method

Modelling small block aperture in an in-house developed GPU-accelerated Monte Carlo-based dose engine for pencil beam scanning proton therapy

Purpose. To enhance an in-house graphic-processing-unit accelerated virtual particle (VP)-based Monte Carlo (MC) proton dose engine (VPMC) to model aperture blocks in both dose calculation and optimization for pencil beam scanning proton therapy (PBSPT)-based stereotactic radiosurgery (SRS).Methods and materials. A module to simulate VPs passing through patient-specific aperture blocks was developed and integrated in VPMC based on simulation results of realistic particles (primary protons and their secondaries). To validate the aperture block module, VPMC was first validated by an opensource MC code, MCsquare, in eight water phantom simulations with 3 cm thick brass apertures: four were with aperture openings of 1, 2, 3, and 4 cm without a range shifter, while the other four were with same aperture opening configurations with a range shifter of 45 mm water equivalent thickness. Then, VPMC was benchmarked with MCsquare and RayStation MC for 10 patients with small targets (average volume 8.4 c.c. with range of 0.4-43.3 c.c.). Finally, 3 typical patients were selected for robust optimization with aperture blocks using VPMC.Results. In the water phantoms, 3D gamma passing rate (2%/2 mm/10%) between VPMC and MCsquare was 99.71 ± 0.23%. In the patient geometries, 3D gamma passing rates (3%/2 mm/10%) between VPMC/MCsquare and RayStation MC were 97.79 ± 2.21%/97.78 ± 1.97%, respectively. Meanwhile, the calculation time was drastically decreased from 112.45 ± 114.08 s (MCsquare) to 8.20 ± 6.42 s (VPMC) with the same statistical uncertainties of ~0.5%. The robustly optimized plans met all the dose-volume-constraints (DVCs) for the targets and OARs per our institutional protocols. The mean calculation time for 13 influence matrices in robust optimization by VPMC was 41.6 s and the subsequent on-the-fly 'trial-and-error' optimization procedure took only 71.4 s on average for the selected three patients.Conclusion. VPMC has been successfully enhanced to model aperture blocks in dose calculation and optimization for the PBSPT-based SRS.

Optimization of proton pencil beam positioning in collimated fields

BACKGROUND

The addition of static or dynamic collimator systems to the pencil beam scanning delivery technique increases the number of options for lateral field shaping. The collimator shape needs to be optimized together with the intensity modulation of spots.

PURPOSE

To minimize the proton field's lateral penumbra by investigating the fundamental relations between spot and collimating aperture edge position.

METHODS

Analytical approaches describing the effect of spot position on the resulting spot profile are presented. The theoretical description is then compared with Monte Carlo simulations in TOPAS and in the RayStation treatment planning system, as well as with radiochromic film measurements at a clinical proton therapy facility. In the model, one single spot profile is analyzed for various spot positions in air. Further, irradiation setups in water with different energies, the combination with a range shifter, and two-dimensional proton fields were investigated in silico.

RESULTS

The further the single spot is placed beyond the collimating aperture edge ('overscanning'), the sharper the relative lateral dose fall-off and thus the lateral penumbra. Overscanning up to 5 mm $5\,\text{mm}$ reduced the lateral penumbra by about 20% on average after a propagation of 13 cm $13\,\text{cm}$ in air. This benefit from overscanning is first predicted by the analytical proofs and later verified by simulations and measurements. Corresponding analyses in water confirm the benefit in lateral penumbra with spot position optimization as observed theoretically and in air. The combination of spot overscanning with fluence modulation facilitated an additional improvement.

CONCLUSIONS

The lateral penumbra of single spots in collimated scanned proton fields can be improved by the method of spot overscanning. This suggests a better sparing of proximal organs at risk in smaller water depths at higher energies, especially in the plateau of the depth dose distribution. All in all, spot overscanning in collimated scanned proton fields offers particular potential in combination with techniques such as fluence modulation or dynamic collimation for optimizing the lateral penumbra to spare normal tissue.

Evaluation of the dosimetric effect of scattered protons in clinical practice in passive scattering proton therapy

The present study verified and evaluated the dosimetric effects of protons scattered from a snout and an aperture in clinical practice, when a range compensator was included. The dose distribution calculated by a treatment planning system (TPS) was compared with the measured dose distribution and the dose distribution calculated by Monte Carlo simulation at several depths. The difference between the measured and calculated results was analyzed using Monte Carlo simulation with filtration of scattering in the snout and aperture. The dependence of the effects of scattered protons on snout size, beam range, and minimum thickness of the range compensator was also investigated using the Monte Carlo simulation. The simulated and measured results showed that the additional dose compared with the results calculated by the TPS at shallow depths was mainly due to protons scattered by the snout and aperture. This additional dose was filtered by the structure of the range compensator so that it was observed under the thin region of the range compensator. The maximum difference was measured at a depth of 16 mm (8.25%), with the difference decreasing with depth. Analysis of protons contributing to the additional dose showed that the contribution of protons scattered from the snout was greater than that of protons scattered from the aperture when a narrow snout was used. In the Monte Carlo simulation, this effect of scattered protons was reduced when wider snouts and longer‐range proton beams were used. This effect was also reduced when thicker range compensator bases were used, even with a narrow snout. This study verified the effect of scattered protons even when a range compensator was included and emphasized the importance of snout‐scattered protons when a narrow snout is used for small fields. It indicated that this additional dose can be reduced by wider snouts, longer range proton beams, and thicker range compensator bases. These results provide a better understanding of the additional dose from scattered protons in clinical practice.

Reduction of superficial radiation dose with bolus in passive scattering proton beam therapy

Purpose

In passive scattering proton beam therapy, scattered protons from the snout and aperture increase the superficial dose, however, treatment planning systems (TPSs) based on analytic algorithms (such as proton convolution superposition) are often inaccurate in this aspect. This additional dose can cause permanent alopecia or severe radiation dermatitis. This study aimed to evaluate the effect of bolus on the superficial radiation dose in passive scattering proton beam therapy.

Methods

We drew a clinical target volume (CTV) and a scalp‐p (phantom), and created plans using a TPS for a solid water phantom with and without bolus. We calculated the dose distribution in the established plans independently with Monte Carlo (MC) simulation and measured the actual dose distribution with an array of ion chambers and radiochromic films. To assess the clinical impact of bolus on scalp dose, we conducted independent dose verification using MC simulation in a clinical case.

Results

In the solid water phantom without bolus, the calculated scalp‐p volume receiving 190 cGy was 20% with TPS but 80% with MC simulation when the CTV received 200 cGy. With 2 cm bolus, this decreased from 80% to 10% in MC simulation. With the measurements, average superficial dose to the scalp‐p was reduced by 5.2% when 2 cm bolus was applied. In the clinical case, the scalp‐c (clinical) volume receiving 3000 cGy decreased from 74% to 63% when 2 cm bolus was applied.

Conclusion

This study revealed that bolus can reduce radiation dose at the superficial body area and alleviate toxicity in passive scattering proton beam therapy.

DICOM-RT Ion interface to utilize MC simulations in routine clinical workflow for proton pencil beam radiotherapy

To adopt Monte Carlo (MC) simulations as an independent dose calculation method for proton pencil beam radiotherapy, an interface that converts the plan information in DICOM format into MC components such as geometries and beam source is a crucial element. For this purpose, a DICOM-RT Ion interface (https://github.com/topasmc/dicom-interface) has been developed and integrated into the TOPAS MC code to perform such conversions on-the-fly. DICOM-RT objects utilized in this interface include Ion Plan (RTIP), Ion Beams Treatment Record (RTIBTR), CT image, and Dose. Beamline geometries, gantry and patient coordinate systems, and fluence maps are determined from RTIP and/or RTIBTR. In this interface, DICOM information is processed and delivered to a MC engine in two steps. A MC model, which consists of beamline geometries and beam source, to represent a treatment machine is created by a DICOM parser of the interface. The complexities from different DICOM types, various beamline configurations and source models are handled in this step. Next, geometry information and beam source are transferred to TOPAS on-the-fly via the developed TOPAS extensions. This interface with two treatment machines was successfully deployed into our automated MC workflow which provides simulated dose and LET distributions in a patient or a water phantom automatically when a new plan is identified. The developed interface provides novel features such as handling multiple treatment systems based on different DICOM types, DICOM conversions on-the-fly, and flexible sampling methods that significantly reduce the burden of handling DICOM based plan or treatment record information for MC simulations.

THE TOP-IMPLART PROTON LINEAR ACCELERATOR: INTERIM CHARACTERISTICS OF THE 35 MEV BEAM

In the framework of the Italian TOP-IMPLART project (Regione Lazio), ENEA-Frascati, ISS and IFO are developing and constructing the first proton linear accelerator based on an actively scanned beam for tumor radiotherapy with final energy of 150 MeV. An important feature of this accelerator is modularity: an exploitable beam can be delivered at any stage of its construction, which allows for immediate characterization and virtually continuous improvement of its performance. Currently, a sequence of 3 GHz accelerating modules combined with a commercial injector operating at 425 MHz delivers protons up to 35 MeV. Several dosimetry systems were used to obtain preliminary characteristics of the 35-MeV beam in terms of stability and homogeneity. Short-term stability and homogeneity better than 3% and 2.6%, respectively, were demonstrated; for stability an improvement with respect to the respective value obtained for the previous 27 MeV beam.

Physical and biological impacts of collimator-scattered protons in spot-scanning proton therapy

To improve the penumbra of low‐energy beams used in spot‐scanning proton therapy, various collimation systems have been proposed and used in clinics. In this paper, focused on patient‐specific brass collimators, the collimator‐scattered protons' physical and biological effects were investigated. The Geant4 Monte Carlo code was used to model the collimators mounted on the scanning nozzle of the Hokkaido University Hospital. A systematic survey was performed in water phantom with various‐sized rectangular targets; range (5–20 cm), spread‐out Bragg peak (SOBP) (5–10 cm), and field size (2 × 2–16 × 16 cm²). It revealed that both the range and SOBP dependences of the physical dose increase had similar trends to passive scattering methods, that is, it increased largely with the range and slightly with the SOBP. The physical impact was maximized at the surface (3%–22% for the tested geometries) and decreased with depth. In contrast, the field size (FS) dependence differed from that observed in passive scattering: the increase was high for both small and large FSs. This may be attributed to the different phase‐space shapes at the target boundary between the two dose delivery methods. Next, the biological impact was estimated based on the increase in dose‐averaged linear energy transfer (LET_d) and relative biological effectiveness (RBE). The LET_d of the collimator‐scattered protons were several keV/μm higher than that of unscattered ones; however, since this large increase was observed only at the positions receiving a small scattered dose, the overall LET_d increase was negligible. As a consequence, the RBE increase did not exceed 0.05. Finally, the effects on patient geometries were estimated by testing two patient plans, and a negligible RBE increase (0.9% at most in the critical organs at surface) was observed in both cases. Therefore, the impact of collimator‐scattered protons is almost entirely attributed to the physical dose increase, while the RBE increase is negligible.

Effect of Irradiation Time on Biological Effectiveness and Tumor Control Probability in Proton Therapy

PURPOSE

The biological effectiveness of proton beams may decrease with irradiation time because of sublethal damage repair (SLDR). The purpose of this study is to systematically evaluate this effect in hypofractionated proton therapy for various target sizes, depths, and prescribed doses per fraction.

METHODS AND MATERIALS

Plans with a single spread-out Bragg peak beam were created using a constant relative biological effectiveness (RBE) of 1.1 to cover targets of 6 different sizes located at 3 different depths in water. Biological doses of 2, 3, 5, 10, and 20 Gy (RBE) were prescribed to the targets. First, to investigate the depth variation of the biological effectiveness, the biological dose in instantaneous irradiation was recalculated based on the microdosimetric kinetic model. SLDR was then taken into account in the microdosimetric kinetic model during treatments to obtain the irradiation time-dependent biological effectiveness for irradiation time T of 5 to 60 minutes and beam interruption time τ of 0 to 60 minutes. The tumor control probabilities were calculated for single-fraction proton therapy fields of different Ts and τs, and the curative doses were evaluated at a tumor control probability of 90%.

RESULTS

The biological effectiveness decreased with longer T and τ and higher prescribed dose. The maximum decrease in the biological effectiveness was 21% with a 20 Gy (RBE) prescribed dose. In single-fraction proton therapy, the curative dose increased linearly by approximately 33% to 35% with the increase of T from 0 to 60 minutes.

CONCLUSIONS

The biological effectiveness varies largely with T and τ because of SLDR during treatments. This effect was pronounced for high prescribed doses per fraction. Thus, the effect of SLDR needs to be considered in hypofractionated and single-fraction proton therapies in relation to size and depth of the target.

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.

A REVIEW OF ANALYTICAL MODELS OF STRAY RADIATION EXPOSURES FROM PHOTON- AND PROTON-BEAM RADIOTHERAPIES

External-beam radiation therapy is safe, effective and widely used to treat cancer. With 5-year cancer survival for adults above 70%, increasingly research is focusing on quantifying and reducing treatment-related morbidity. Reducing exposures to healthy tissues is one strategy, which can be accomplished with advanced-technology radiotherapies, such as intensity-modulated photon therapy and proton therapy. Both of these modalities provide good conformation of the therapeutic dose to the tumor volume, but they also deliver stray radiation to the whole body that increases the risk of radiogenic second cancers. To minimize these risks, one needs to create and compare candidate treatment plans that explicitly take into account these risks. Currently, clinical practice does not include routine calculation of stray radiation exposure and, consequently, the assessment of corresponding risks is difficult. In this article, we review recent progress toward stray dose algorithms that are suitable for large-scale clinical use. In particular, we emphasize the current state of physics-based dose algorithms for intensity-modulated photon radiotherapy and proton therapy.

A correction scheme for a simplified analytical random walk model algorithm of proton dose calculation in distal Bragg peak regions

The lateral homogeneity assumption is used in most analytical algorithms for proton dose, such as the pencil-beam algorithms and our simplified analytical random walk model. To improve the dose calculation in the distal fall-off region in heterogeneous media, we analyzed primary proton fluence near heterogeneous media and propose to calculate the lateral fluence with voxel-specific Gaussian distributions. The lateral fluence from a beamlet is no longer expressed by a single Gaussian for all the lateral voxels, but by a specific Gaussian for each lateral voxel. The voxel-specific Gaussian for the beamlet of interest is calculated by re-initializing the fluence deviation on an effective surface where the proton energies of the beamlet of interest and the beamlet passing the voxel are the same. The dose improvement from the correction scheme was demonstrated by the dose distributions in two sets of heterogeneous phantoms consisting of cortical bone, lung, and water and by evaluating distributions in example patients with a head-and-neck tumor and metal spinal implants. The dose distributions from Monte Carlo simulations were used as the reference. The correction scheme effectively improved the dose calculation accuracy in the distal fall-off region and increased the gamma test pass rate. The extra computation for the correction was about 20% of that for the original algorithm but is dependent upon patient geometry.

Simplified estimation method for dose distributions around field junctions in proton craniospinal irradiation

In radiotherapy involving craniospinal irradiation (CSI), field junctions of therapeutic beams are necessary, because a CSI target is generally several times larger than the maximum field size of the beams. The purpose of this study was to develop a simplified method for estimating dose uniformity around the field junctions in proton CSI. We estimated the dose profiles around the field junctions of proton beams using a simplified field-junction model, in which partial lateral dose distributions around the field edge were assumed to be approximated using the error function. We measured the lateral dose distributions of the proton beams planned for the CSI treatment using a two-dimensional (2D) ionization chamber array. Although dose hot spots and cold spots tend to be underestimated by a chamber array because of the partial volume effect of the sensitive volume and discrete chamber positions, the model estimation results were fairly consistent with the measurements obtained using a 2D chamber array subjected to CSI-simulated serial irradiation. The simplified junction model enabled us to estimate the dose distributions and dependence of the setup position gap on the dose uniformity around the field junctions on the basis of the field-by-field dose profiles measured using the 2D chamber array.

A method for modeling laterally asymmetric proton beamlets resulting from collimation

PURPOSE

To introduce a method to model the 3D dose distribution of laterally asymmetric proton beamlets resulting from collimation. The model enables rapid beamlet calculation for spot scanning (SS) delivery using a novel penumbra-reducing dynamic collimation system (DCS) with two pairs of trimmers oriented perpendicular to each other.

METHODS

Trimmed beamlet dose distributions in water were simulated with MCNPX and the collimating effects noted in the simulations were validated by experimental measurement. The simulated beamlets were modeled analytically using integral depth dose curves along with an asymmetric Gaussian function to represent fluence in the beam's eye view (BEV). The BEV parameters consisted of Gaussian standard deviations (sigmas) along each primary axis (σ(x1),σ(x2),σ(y1),σ(y2)) together with the spatial location of the maximum dose (μ(x),μ(y)). Percent depth dose variation with trimmer position was accounted for with a depth-dependent correction function. Beamlet growth with depth was accounted for by combining the in-air divergence with Hong's fit of the Highland approximation along each axis in the BEV.

RESULTS

The beamlet model showed excellent agreement with the Monte Carlo simulation data used as a benchmark. The overall passing rate for a 3D gamma test with 3%/3 mm passing criteria was 96.1% between the analytical model and Monte Carlo data in an example treatment plan.

CONCLUSIONS

The analytical model is capable of accurately representing individual asymmetric beamlets resulting from use of the DCS. This method enables integration of the DCS into a treatment planning system to perform dose computation in patient datasets. The method could be generalized for use with any SS collimation system in which blades, leaves, or trimmers are used to laterally sharpen beamlets.

A Review of Radiotherapy-Induced Late Effects Research after Advanced Technology Treatments

The number of incident cancers and long-term cancer survivors is expected to increase substantially for at least a decade. Advanced technology radiotherapies, e.g., using beams of protons and photons, offer dosimetric advantages that theoretically yield better outcomes. In general, evidence from controlled clinical trials and epidemiology studies are lacking. To conduct these studies, new research methods and infrastructure will be needed. In the paper, we review several key research methods of relevance to late effects after advanced technology proton-beam and photon-beam radiotherapies. In particular, we focus on the determination of exposures to therapeutic and stray radiation and related uncertainties, with discussion of recent advances in exposure calculation methods, uncertainties, in silico studies, computing infrastructure, electronic medical records, and risk visualization. We identify six key areas of methodology and infrastructure that will be needed to conduct future outcome studies of radiation late effects.

An empirical model for calculation of the collimator contamination dose in therapeutic proton beams

Collimators are used as lateral beam shaping devices in proton therapy with passive scattering beam lines. The dose contamination due to collimator scattering can be as high as 10% of the maximum dose and influences calculation of the output factor or monitor units (MU). To date, commercial treatment planning systems generally use a zero-thickness collimator approximation ignoring edge scattering in the aperture collimator and few analytical models have been proposed to take scattering effects into account, mainly limited to the inner collimator face component. The aim of this study was to characterize and model aperture contamination by means of a fast and accurate analytical model. The entrance face collimator scatter distribution was modeled as a 3D secondary dose source. Predicted dose contaminations were compared to measurements and Monte Carlo simulations. Measurements were performed on two different proton beam lines (a fixed horizontal beam line and a gantry beam line) with divergent apertures and for several field sizes and energies. Discrepancies between analytical algorithm dose prediction and measurements were decreased from 10% to 2% using the proposed model. Gamma-index (2%/1 mm) was respected for more than 90% of pixels. The proposed analytical algorithm increases the accuracy of analytical dose calculations with reasonable computation times.

Use of diverging apertures to minimize the edge scatter in passive scattering proton therapy

The purpose of this study was to evaluate the use of diverging-cut aperture to minimize collimator contamination in proton therapy. Two sets of apertures with nondivergent and divergent edge were fabricated to produce a 10 cm × 10 cm field at the radiation isocenter of a single-room proton therapy unit. Transverse profiles were acquired in a scanning water tank with both aperture sets. Up to 9.5% extra dose was observed from aperture scattering near the field edges with the nondivergent aperture set at 2 cm above the water surface and remained 3.0% at depth of 10 cm. For the divergent set, the contamination was reduced to less than 3.5% and 1.3%, respectively. Our study demonstrated that scattering from apertures contaminated the dose distribution near the field edge at shallow depth. A diverging-cut aperture was capable of reducing the contamination and is recommended for use in passive scattering proton therapy, especially when critical organs are lateral and proximal to the target at shallow depth.

Evaluation of monitor unit calculation based on measurement and calculation with a simplified Monte Carlo method for passive beam delivery system in proton beam therapy

Calibrating the dose per monitor unit (DMU) for individual patients is important to deliver the prescribed dose in radiation therapy. We have developed a DMU calculation method combining measurement data and calculation with a simplified Monte Carlo method for the double scattering system in proton beam therapy at the National Cancer Center Hospital East in Japan. The DMU calculation method determines the clinical DMU by the multiplication of three factors: a beam spreading device factor FBSD, a patient‐specific device factor FPSD, and a field‐size correction factor FFS(A). We compared the calculated and the measured DMU for 75 dose fields in clinical cases. The calculated DMUs were in agreement with measurements in ±1.1% for all of 25 fields in prostate cancer cases, and in ±3% for 94% of 50 fields in head and neck (H&N) and lung cancer cases, including irregular shape fields and small fields. Although the FBSD in the DMU calculations is dominant as expected, we found that the patient‐specific device factor and field‐size correction also contribute significantly to the calculated DMU. This DMU calculation method will be able to substitute the conventional DMU measurement for the majority of clinical cases with a reasonable calculation time required for clinical use.

PACS number: 87.55.kh

Semi-analytical model for output factor calculations in proton beam therapy with consideration for the collimator aperture edge

In the development of an external radiotherapy treatment planning system, the output factor (OPF) is an important value for the monitor unit calculations. We developed a proton OPF calculation model with consideration for the collimator aperture edge to account for the dependence of the OPF on the collimator aperture and distance in proton beam therapy. Five parameters in the model were obtained by fitting with OPFs measured by a pinpoint chamber with the circular radiation fields of various field radii and collimator distances. The OPF model calculation using the fitted model parameters could explain the measurement results to within 1.6% error in typical proton treatment beams with 6- and 12 cm SOBP widths through a range shifter and a circular aperture more than 10.6 mm in radius. The calibration depth dependences of the model parameters were approximated by linear or quadratic functions. The semi-analytical OPF model calculation was tested with various MLC aperture shapes that included circles of various sizes as well as a rectangle, parallelogram, and L-shape for an intermediate proton treatment beam condition. The pre-calculated OPFs agreed well with the measured values, to within 2.7% error up to 620 mm in the collimator distance, though the maximum difference was 5.1% in the case of the largest collimator distance of 740 mm. The OPF calculation model would allow more accurate monitor unit calculations for therapeutic proton beams within the expected range of collimator conditions in clinical use.

Reducing the cost of proton radiation therapy: the feasibility of a streamlined treatment technique for prostate cancer

Proton radiation therapy is an effective modality for cancer treatments, but the cost of proton therapy is much higher compared to conventional radiotherapy and this presents a formidable barrier to most clinical practices that wish to offer proton therapy. Little attention in literature has been paid to the costs associated with collimators, range compensators and hypofractionation. The objective of this study was to evaluate the feasibility of cost-saving modifications to the present standard of care for proton treatments for prostate cancer. In particular, we quantified the dosimetric impact of a treatment technique in which custom fabricated collimators were replaced with a multileaf collimator (MLC) and the custom range compensators (RC) were eliminated. The dosimetric impacts of these modifications were assessed for 10 patients with a commercial treatment planning system (TPS) and confirmed with corresponding Monte Carlo simulations. We assessed the impact on lifetime risks of radiogenic second cancers using detailed dose reconstructions and predictive dose-risk models based on epidemiologic data. We also performed illustrative calculations, using an isoeffect model, to examine the potential for hypofractionation. Specifically, we bracketed plausible intervals of proton fraction size and total treatment dose that were equivalent to a conventional photon treatment of 79.2 Gy in 44 fractions. Our results revealed that eliminating the RC and using an MLC had negligible effect on predicted dose distributions and second cancer risks. Even modest hypofractionation strategies can yield substantial cost savings. Together, our results suggest that it is feasible to modify the standard of care to increase treatment efficiency, reduce treatment costs to patients and insurers, while preserving high treatment quality.

The physics of proton therapy

The physics of proton therapy has advanced considerably since it was proposed in 1946. Today analytical equations and numerical simulation methods are available to predict and characterize many aspects of proton therapy. This article reviews the basic aspects of the physics of proton therapy, including proton interaction mechanisms, proton transport calculations, the determination of dose from therapeutic and stray radiations, and shielding design. The article discusses underlying processes as well as selected practical experimental and theoretical methods. We conclude by briefly speculating on possible future areas of research of relevance to the physics of proton therapy.

A method to select aperture margin in collimated spot scanning proton therapy

The use of collimator or aperture may sharpen the lateral dose gradient for spot scanning proton therapy. However, to date, there has not been a standard method to determine the aperture margin for a single field in collimated spot scanning proton therapy. This study describes a theoretical framework to select the optimal aperture margin for a single field, and also presents the spot spacing limit required such that the optimal aperture margin exists. Since, for a proton pencil beam partially intercepted by collimator, the maximum point dose (spot center) shifts away from the original pencil beam central axis, we propose that the optimal margin should be equal to the maximum pencil beam center shift under the condition that spot spacing is small with respect to the maximum pencil beam center shift, which can be numerically determined based on beam modeling data. A test case is presented which demonstrates agreement with the prediction made based on the proposed methods. When apertures are applied in a commercial treatment planning system this method may be implemented.

Calculations and measurements of the scintillator-to-water stopping power ratio of liquid scintillators for use in proton radiotherapy

Liquid scintillators are a promising detector for high-resolution three-dimensional proton therapy dosimetry. Because the scintillator comprises both the active volume of the detector and the phantom material, an ideal scintillator will exhibit water equivalence in its radiological properties. One of the most fundamental of these is the scintillator's stopping power. The objective of this study was to compare calculations and measurements of scintillator-to-water stopping power ratios to evaluate the suitability of the liquid scintillators BC-531 and OptiPhase HiSafe 3 for proton dosimetry. We also measured the relative scintillation output of the two scintillators. Both calculations and measurements show that the linear stopping power of OptiPhase is significantly closer to water than that of BC-531. BC-531 has a somewhat higher scintillation output. OptiPhase can be mixed with water at high concentrations, which further improves its scintillator-to-water stopping power ratio. However, this causes the solution to become cloudy, which has a negative impact on the scintillation output and spatial resolution of the detector. OptiPhase is preferred over BC-531 for proton dosimetry because its density and scintillator-to-water stopping power ratio are more water equivalent.

Radiation burden from secondary doses to patients undergoing radiation therapy with photons and light ions and radiation doses from imaging modalities

Ionising radiation is increasingly used for the treatment of cancer, being the source of a considerable fraction of the medical irradiation to patients. With the increasing success rate of cancer treatments and longer life expectancy of the treated patients, the issue of secondary cancer incidence is of growing concern, especially for paediatric patients who may live long after the treatment and be more susceptible to carcinogenesis. Also, additional imaging procedures like computed tomography, kilovoltage and megavoltage imaging and positron emission tomography, alone or in conjunction with radiation therapy, may add to the radiation burden associated with the risk of occurrence of secondary cancers. This work has been based on literature studies and is focussed on the assessment of secondary doses to healthy tissues that are delivered by the use of modern radiation therapy and diagnostic imaging modalities in the clinical environment.

Implementation of an improved dose-per-MU model for double-scattered proton beams to address interbeamline modulation width variability

Because treatment planning systems (TPSs) generally do not provide monitor units (MUs) for double‐scattered proton plans, models to predict MUs as a function of the range and the nominal modulation width requested of the beam delivery system, such as the one developed by the MGH group, have been proposed. For a given nominal modulation width, however, the measured modulation width depends on the accuracy of the vendor's calibration process and may differ from this nominal value, and also from one beamline to the next. Although such a difference can be replicated in our TPS, the output dependence on range and modulation width for each beam option or suboption has to be modeled separately for each beamline in order to achieve maximal 3% inaccuracy. As a consequence, the MGH output model may not be directly transferable. This work, therefore, serves to extend the model to more general clinic situations. In this paper, a parameterized linear‐quadratic transformation is introduced to convert the nominal modulation width to the measured modulation width for each beam option or suboption on a per‐beamline basis. Fit parameters are derived for each beamline from measurements of 60 reference beams spanning the minimum and maximum ranges, and modulation widths from 2 cm to full range per option or suboption. Using the modeled modulation width, we extract the MGH parameters for the output dependence on range and modulation width. Our method has been tested with 1784 patient‐specific fields delivered across three different beamlines at our facility. For these fields, all measured outputs fall within 3%, and 64.4% fall within 1%, of our model. Using a parameterized linear‐quadratic modulation width, MU calculation models can be established on a per‐beamline basis for each double scattering beam option or suboption.

PACS number: 87.53.Qc

An algorithm to assess the need for clinical Monte Carlo dose calculation for small proton therapy fields based on quantification of tissue heterogeneity

PURPOSE

In proton therapy, complex density heterogeneities within the beam path constitute a challenge to dose calculation algorithms. This might question the reliability of dose distributions predicted by treatment planning systems based on analytical dose calculation. For cases in which substantial dose errors are expected, resorting to Monte Carlo dose calculations might be essential to ensure a successful treatment outcome and therefore the benefit is worth a presumably long computation time. The aim of this study was to define an indicator for the accuracy of dose delivery based on analytical dose calculations in treatment planning systems for small proton therapy fields to identify those patients for which Monte Carlo dose calculation is warranted.

METHODS

Fourteen patients treated at our facility with small passively scattered proton beams (apertures diameters below 7 cm) were selected. Plans were generated in the XiO treatment planning system in combination with a pencil beam algorithm developed at the Massachusetts General Hospital and compared to Monte Carlo dose calculations. Differences in the dose to the 50% of the gross tumor volume (D50, GTV) were assessed in a field-by-field basis. A simple and fast methodology was developed to quantify the inhomogeneity of the tissue traversed by a single small proton beam using a heterogeneity index (HI)-a concept presented by Plugfelder et al. [Med. Phys. 34, 1506-1513 (2007)] for scanned proton beams. Finally, the potential correlation between the error made by the pencil beam based treatment planning algorithm for each field and the level of tissue heterogeneity traversed by the proton beam given by the HI was evaluated.

RESULTS

Discrepancies up to 5.4% were found in D50 for single fields, although dose differences were within clinical tolerance levels (<3%) when combining all of the fields involved in the treatment. The discrepancies found for each field exhibited a strong correlation to their associated HI-values (Spearman's ρ=0.8, p<0.0001); the higher the level of tissue inhomogeneities for a particular field, the larger the error by the analytical algorithm. With the established correlation a threshold for HI can be set by choosing a tolerance level of 2-3%-commonly accepted in radiotherapy.

CONCLUSIONS

The HI is a good indicator for the accuracy of proton field delivery in terms of GTV prescription dose coverage when small fields are delivered. Each HI-value was obtained from the CT image in less than 3 min on a computer with 2 GHz CPU allowing implementation of this methodology in clinical routine. For HI-values exceeding the threshold, either a change in beam direction (if feasible) or a recalculation of the dose with Monte Carlo would be highly recommended.

Quantification of dose perturbations induced by external and internal accessories in ocular proton therapy and evaluation of their dosimetric impact

PURPOSE

Proton scattering on beam shaping devices and protons slowing down on media with different densities within the treatment volume may produce dose perturbations and range variations that are not predicted by treatment planning systems. The aim of this work was to assess the dosimetric impact of elements present in ocular proton therapy treatments that may disturb the prescribed treatment plan. Both distal beam shaping devices and intraocular elements were considered.

METHODS

A wedge filter, tantalum fiducial marker, hemispherical compensator, two intraocular endotamponades (densities 0.97 and 1.92 g cm(-3)) and an intraocular eye lens (IOL) were considered in the study. For these elements, longitudinal dose distributions were measured and∕or calculated in water in beam alignment for a clinical spread-out Bragg peak. Under the same conditions, the unperturbed dose distributions were similarly measured and∕or calculated in the absence of the element. The dosimetric impact was assessed by comparison of unperturbed and perturbed dose distributions. Measurements and calculations were carried out with two methods. Measurements are based on EBT3 films with dedicated software, which makes use of a calibration curve and correction for the quenching effect. Calculations are based on the Monte Carlo (MC) code MCNPX and reproduce the experimental conditions. Both dose maps are obtained with a resolution of 300 dpi.

RESULTS

The degree of disturbance of distal beam shaping devices is low for the wedge filter (2% overdose ripple all along the central axis) and moderate for the hemispherical compensator (two bands of variable overdose of up to 10% downstream the compensator lateral edges and -5% underdose on the plateau at off-axis distance of 5 cm). Tantalum clips produce important dose shadows (-20% behind the clip parallel to the beam and range reduction of 1.1 mm) and bands of overdose (15%). The presence of endotamponades modifies the dose distribution very significantly (-5% underdose on the plateau and 3 mm range prolongation for the tamponade with density 0.97 g cm(-3) and -15% underdose on plateau and 8 mm range reduction for that with density 1.92 g cm(-3)). No dose perturbations were found for the IOL. The high performance of EBT3 film and MC tools used was confirmed and good agreement was found between them (percentage of pixels passing the gamma test >87%).

CONCLUSIONS

The degree of disturbance by external beam shaping devices remains low in ocular proton therapy and can be reduced by bringing accessories closer to the eye. Tantalum fiducial markers must be located in such a way that dose perturbation is not projected on the tumor. The treatment of patients with intraocular endotamponades must be carefully managed. It is fundamental that radiation oncologists and medical physicists are informed about the presence of such substances prior to the treatment.

Experimental verification of dose calculation using the simplified Monte Carlo method with an improved initial beam model for a beam-wobbling system

A beam delivery system using a single-radius-beam-wobbling method has been used to form a conformal irradiation field for proton radiotherapy in Japan. A proton beam broadened by the beam-wobbling system provides a non-Gaussian distribution of projection angle different in two mutually orthogonal planes with a common beam central axis, at a certain position. However, the conventional initial beam model for dose calculations has been using an approximation of symmetric Gaussian angular distribution with the same variance in both planes (called here a Gaussian model with symmetric variance (GMSV)), instead of the accurate one. We have developed a more accurate initial beam model defined as a non-Gaussian model with asymmetric variance (NonGMAV), and applied it to dose calculations using the simplified Monte Carlo (SMC) method. The initial beam model takes into account the different distances of two beam-wobbling magnets from the iso-center and also the different amplitudes of kick angle given by each magnet. We have confirmed that the calculation using the SMC with NonGMAV reproduced the measured dose distribution formed in air by a mono-energetic proton beam passing through a square aperture collimator better than with the GMSV and with a Gaussian model with asymmetric variance (GMAV) in which different variances of angular distributions are used in the two mutually orthogonal planes. Measured dose distributions in a homogeneous phantom formed by a modulated proton beam passing through a range shifter and an L-shaped range compensator, were consistent with calculations using the SMC with GMAV and NonGMAV, but in disagreement with calculations using the SMC with GMSV. Measured lateral penumbrae in a lateral direction were reproduced better by calculations using the SMC with NonGMAV than by those with GMAV, when an aperture collimator with a smaller opening was used. We found that such a difference can be attributed to the non-Gaussian angular distribution of the initial beam at a lateral position for the beam-wobbling system. Calculations using the SMC with NonGMAV are effective to reproduce dose distributions formed by a beam-wobbling system more accurately than that with GMSV or that with GMAV.

Determination of the quenching correction factors for plastic scintillation detectors in therapeutic high-energy proton beams

Plastic scintillation detectors (PSDs) have many advantages over other detectors in small field dosimetry due to their high spatial resolution, excellent water equivalence and instantaneous readout. However, in proton beams, the PSDs undergo a quenching effect which makes the signal level reduced significantly when the detector is close to the Bragg peak where the linear energy transfer (LET) for protons is very high. This study measures the quenching correction factor (QCF) for a PSD in clinical passive-scattering proton beams and investigates the feasibility of using PSDs in depth-dose measurements in proton beams. A polystyrene-based PSD (BCF-12, ϕ0.5 mm × 4 mm) was used to measure the depth-dose curves in a water phantom for monoenergetic unmodulated proton beams of nominal energies 100, 180 and 250 MeV. A Markus plane-parallel ion chamber was also used to get the dose distributions for the same proton beams. From these results, the QCF as a function of depth was derived for these proton beams. Next, the LET depth distributions for these proton beams were calculated by using the MCNPX Monte Carlo code, based on the experimentally validated nozzle models for these passive-scattering proton beams. Then the relationship between the QCF and the proton LET could be derived as an empirical formula. Finally, the obtained empirical formula was applied to the PSD measurements to get the corrected depth-dose curves and they were compared to the ion chamber measurements. A linear relationship between the QCF and LET, i.e. Birks' formula, was obtained for the proton beams studied. The result is in agreement with the literature. The PSD measurements after the quenching corrections agree with ion chamber measurements within 5%. PSDs are good dosimeters for proton beam measurement if the quenching effect is corrected appropriately.

Feasibility of using PRESAGE® for relative 3D dosimetry of small proton fields

Small field dosimetry is challenging due to the finite size of the conventional detectors that underestimate the dose distribution. With the fast development of the dynamic proton beam delivery system, it is essential to find a dosimeter which can be used for 3D dosimetry of small proton fields. We investigated the feasibility of using a proton formula PRESAGE® for 3D dosimetry of small fields in a uniform scanning proton beam delivery system with dose layer stacking technology. The relationship between optical density and the absorbed dose was found to be linear through small volume cuvette studies for both photon and proton irradiation. Two circular fields and three patient-specific fields were used for proton treatment planning calculation and beam delivery. The measured results were compared with the calculated results in the form of lateral dose profiles, depth dose, isodose plots and gamma index analysis. For the circular field study, lateral dose profile comparison showed that the relative PRESAGE® profile falls within ± 5% from the calculated profile for most of the spatial range. For unmodulated depth dose comparison, the agreement between the measured and calculated results was within 3% in the beam entrance region before the Bragg peak. However, at the Bragg peak, there was about 20% underestimation of the absorbed dose from PRESAGE®. For patient-specific field 3D dosimetry, most of the data points within the target volume passed gamma analysis for 3% relative dose difference and 3 mm distance to agreement criteria. Our results suggest that this proton formula PRESAGE® dosimeter has the potential for 3D dosimetry of small fields in proton therapy, but further investigation is needed to improve the dose under-response of the PRESAGE® in the Bragg peak region.

Range uncertainties in proton therapy and the role of Monte Carlo simulations

The main advantages of proton therapy are the reduced total energy deposited in the patient as compared to photon techniques and the finite range of the proton beam. The latter adds an additional degree of freedom to treatment planning. The range in tissue is associated with considerable uncertainties caused by imaging, patient setup, beam delivery and dose calculation. Reducing the uncertainties would allow a reduction of the treatment volume and thus allow a better utilization of the advantages of protons. This paper summarizes the role of Monte Carlo simulations when aiming at a reduction of range uncertainties in proton therapy. Differences in dose calculation when comparing Monte Carlo with analytical algorithms are analyzed as well as range uncertainties due to material constants and CT conversion. Range uncertainties due to biological effects and the role of Monte Carlo for in vivo range verification are discussed. Furthermore, the current range uncertainty recipes used at several proton therapy facilities are revisited. We conclude that a significant impact of Monte Carlo dose calculation can be expected in complex geometries where local range uncertainties due to multiple Coulomb scattering will reduce the accuracy of analytical algorithms. In these cases Monte Carlo techniques might reduce the range uncertainty by several mm.

Range uncertainties in proton therapy and the role of Monte Carlo simulations

The main advantages of proton therapy are the reduced total energy deposited in the patient as compared to photon techniques and the finite range of the proton beam. The latter adds an additional degree of freedom to treatment planning. The range in tissue is associated with considerable uncertainties caused by imaging, patient setup, beam delivery and dose calculation. Reducing the uncertainties would allow a reduction of the treatment volume and thus allow a better utilization of the advantages of protons. This paper summarizes the role of Monte Carlo simulations when aiming at a reduction of range uncertainties in proton therapy. Differences in dose calculation when comparing Monte Carlo with analytical algorithms are analyzed as well as range uncertainties due to material constants and CT conversion. Range uncertainties due to biological effects and the role of Monte Carlo for in vivo range verification are discussed. Furthermore, the current range uncertainty recipes used at several proton therapy facilities are revisited. We conclude that a significant impact of Monte Carlo dose calculation can be expected in complex geometries where local range uncertainties due to multiple Coulomb scattering will reduce the accuracy of analytical algorithms. In these cases Monte Carlo techniques might reduce the range uncertainty by several mm.

Uncertainties and correction methods when modeling passive scattering proton therapy treatment heads with Monte Carlo

Nowadays, Monte Carlo models of proton therapy treatment heads are being used to improve beam delivery systems and to calculate the radiation field for patient dose calculations. The achievable accuracy of the model depends on the exact knowledge of the treatment head geometry and time structure, the material characteristics, and the underlying physics. This work aimed at studying the uncertainties in treatment head simulations for passive scattering proton therapy. The sensitivities of spread-out Bragg peak (SOBP) dose distributions on material densities, mean ionization potentials, initial proton beam energy spread and spot size were investigated. An improved understanding of the nature of these parameters may help to improve agreement between calculated and measured SOBP dose distributions and to ensure that the range, modulation width, and uniformity are within clinical tolerance levels. Furthermore, we present a method to make small corrections to the uniformity of spread-out Bragg peaks by utilizing the time structure of the beam delivery. In addition, we re-commissioned the models of the two proton treatment heads located at our facility using the aforementioned correction methods presented in this paper.

Dosimetric accuracy of planning and delivering small proton therapy fields

A detailed Monte Carlo model of a proton therapy treatment head was commissioned to simulate the delivery of small field proton treatments. Small fields are challenging in the planning and quality assurance process because of aperture scattering and dosimetric disequilibrium. Four patients with small fields used in all or parts of their treatment course were studied, including two stereotactic patients and two fractionated patients. For the two stereotactic patients the overall difference of the dose covering at least 95% of the gross tumor volume between the Monte Carlo calculations and the delivery was -0.2% and -1.6%, respectively. For the two fractionated patients the overall difference of the dose covering at least 95% of the clinical target volume was -3.0% and 1.0%, respectively. We have thus confirmed that our current planning and delivery process for small proton fields is accurate enough to treat small lesions in the patient. Furthermore, we studied the impact of field size corrections and identified limitations of the pencil beam algorithm for predicting hot and cold spots and range degradation in the target due scattering in heterogeneities. For the four cases studied in this paper, these limitations appear to impact individual field calculations, but do not have a significant impact on the prescribed dose over multiple fields.

An MCNPX Monte Carlo model of a discrete spot scanning proton beam therapy nozzle

PURPOSE

The purposes of this study were to validate a discrete spot scanning proton beam nozzle using the Monte Carlo (MC) code MCNPX and use the MC validated model to investigate the effects of a low-dose envelope, which surrounds the beam's central axis, on measurements of integral depth dose (IDD) profiles.

METHODS

An accurate model of the discrete spot scanning beam nozzle from The University of Texas M. D. Anderson Cancer Center (Houston, Texas) was developed on the basis of blueprints provided by the manufacturer of the nozzle. The authors performed simulations of single proton pencil beams of various energies using the standard multiple Coulomb scattering (MCS) algorithm within the MCNPX source code and a new MCS algorithm, which was implemented in the MCNPX source code. The MC models were validated by comparing calculated in-air and in-water lateral profiles and percentage depth dose profiles for single pencil beams with their corresponding measured values. The models were then further tested by comparing the calculated and measured three-dimensional (3-D) dose distributions. Finally, an IDD profile was calculated with different scoring radii to determine the limitations on the use of commercially available plane-parallel ionization chambers to measure IDD.

RESULTS

The distance to agreement, defined as the distance between the nearest positions of two equivalent distributions with the same value of dose, between measured and simulated ranges was within 0.13 cm for both MCS algorithms. For low and intermediate pencil beam energies, the MC simulations using the standard MCS algorithm were in better agreement with measurements. Conversely, the new MCS algorithm produced better results for high-energy single pencil beams. The IDD profile calculated with cylindrical tallies with an area equivalent to the area of the largest commercially available ionization chamber showed up to 7.8% underestimation of the integral dose in certain depths of the IDD profile.

CONCLUSIONS

The authors conclude that a combination of MCS algorithms is required to accurately reproduce experimental data of single pencil beams and 3-D dose distributions for the scanning beam nozzle. In addition, the MC simulations showed that because of the low-dose envelope, ionization chambers with radii as large as 4.08 cm are insufficient to accurately measure IDD profiles for a 221.8 MeV pencil beam in the scanning beam nozzle.

Development and verification of an analytical algorithm to predict absorbed dose distributions in ocular proton therapy using Monte Carlo simulations

Proton beam radiotherapy is an effective and non-invasive treatment for uveal melanoma. Recent research efforts have focused on improving the dosimetric accuracy of treatment planning and overcoming the present limitation of relative analytical dose calculations. Monte Carlo algorithms have been shown to accurately predict dose per monitor unit (D/MU) values, but this has yet to be shown for analytical algorithms dedicated to ocular proton therapy, which are typically less computationally expensive than Monte Carlo algorithms. The objective of this study was to determine if an analytical method could predict absolute dose distributions and D/MU values for a variety of treatment fields like those used in ocular proton therapy. To accomplish this objective, we used a previously validated Monte Carlo model of an ocular nozzle to develop an analytical algorithm to predict three-dimensional distributions of D/MU values from pristine Bragg peaks and therapeutically useful spread-out Bragg peaks (SOBPs). Results demonstrated generally good agreement between the analytical and Monte Carlo absolute dose calculations. While agreement in the proximal region decreased for beams with less penetrating Bragg peaks compared with the open-beam condition, the difference was shown to be largely attributable to edge-scattered protons. A method for including this effect in any future analytical algorithm was proposed. Comparisons of D/MU values showed typical agreement to within 0.5%. We conclude that analytical algorithms can be employed to accurately predict absolute proton dose distributions delivered by an ocular nozzle.

A sector-integration method for dose/MU calculation in a uniform scanning proton beam

An accurate, simple and time-saving sector integration method for calculating the proton output (dose/monitor unit, MU) is presented based on the following treatment field parameters: aperture shape, aperture size, measuring position, beam range and beam modulation. The model is validated with dose/MU values for 431 fields previously measured at our center. The measurements were obtained in a uniform scanning proton beam with a parallel plate ionization chamber in a water phantom. For beam penetration depths of clinical interest (6-27 cm water), dose/MU values were measured as a function of spread-out Bragg peak (SOBP) extent and aperture diameter. First, 90 randomly selected fields were used to derive the model parameters, which were used to compute the dose/MU values for the remaining 341 fields. The min, max, average and the standard deviation of the difference between the calculated and the measured dose/MU values of the 341 fields were used to evaluate the accuracy and stability, for different energy ranges, aperture sizes, measurement positions and SOBP values. The experimental results of the five different functional sets showed that the calculation model is accurate with calculation errors ranging from -2.4% to 3.3%, and 99% of the errors are less than +/-2%. The accuracy increases with higher energy, larger SOBP and bigger aperture size. The average error in the dose/MU calculation for small fields (field size <25 cm(2)) is 0.31 +/- 0.96 (%).

REDUCING STRAY RADIATION DOSE FOR A PEDIATRIC PATIENT RECEIVING PROTON CRANIOSPINAL IRRADIATION

The aim of this study was to quantify stray radiation dose from neutrons emanating from a proton treatment unit and to evaluate methods of reducing this dose for a pediatric patient undergoing craniospinal irradiation. The organ equivalent doses and effective dose from stray radiation were estimated for a 30.6-Gy treatment using Monte Carlo simulations of a passive scattering treatment unit and a patient-specific voxelized anatomy. The treatment plan was based on computed tomography images of a 10-yr-old male patient. The contribution to stray radiation was evaluated for the standard nozzle and for the same nozzle but with modest modifications to suppress stray radiation. The modifications included enhancing the local shielding between the patient and the primary external neutron source and increasing the distance between them. The effective dose from stray radiation emanating from the standard nozzle was 322 mSv; enhancements to the nozzle reduced the effective dose by as much as 43%. These results add to the body of evidence that modest enhancements to the treatment unit can reduce substantially the effective dose from stray radiation.

Field size dependence of the output factor in passively scattered proton therapy: influence of range, modulation, air gap, and machine settings

At the Francis H. Burr Proton Therapy Center field specific output factors (i.e., dose per monitor unit) for patient treatments were modeled for all beamlines (two gantries, fixed stereotactic, and fixed eye beamline). The authors evaluated the accuracy of dose calculation and output model for small fields. Measurements in a water phantom were performed in three of our beamlines quantifying the dependency of the output factor on the field size for a variety of proton ranges. The influence of snout size, air gap, modulation, and second scatterer was investigated. The impact of field size on output depends strongly on the depth of interest. The air gap has a notable influence on small field outputs. A field size specific correction factor to the output is necessary if the latter was modeled or measured without the custom hardware in place. The output was shown to be field size dependent even for large fields, indicating an effect beyond charged particle disequilibrium caused by lateral scatter.

Commissioning a passive-scattering proton therapy nozzle for accurate SOBP delivery

Proton radiotherapy centers that currently use passively scattered proton beams do field specific calibrations for a non-negligible fraction of treatment fields, which is time and resource consuming. Our improved understanding of the passive scattering mode of the IBA universal nozzle, especially of the current modulation function, allowed us to re-commission our treatment control system for accurate delivery of SOBPs of any range and modulation, and to predict the output for each of these fields. We moved away from individual field calibrations to a state where continued quality assurance of SOBP field delivery is ensured by limited system-wide measurements that only require one hour per week. This manuscript reports on a protocol for generation of desired SOBPs and prediction of dose output.