Beam Position Projection Algorithms in Proton Pencil Beam Scanning.

Beam position uncertainties along the beam trajectory arise from the accelerator, beamline, and scanning magnets (SMs). They can be monitored in real time, e.g., through strip ionization chambers (ICs), and treatments can be paused if needed. Delivery is more reliable and accurate if the beam position is projected from monitored nozzle parameters to the isocenter, allowing for accurate online corrections to be performed. Beam position projection algorithms are also used in post-delivery log file analyses. In this paper, we investigate the four potential algorithms that can be applied to all pencil beam scanning (PBS) nozzles. For some combinations of nozzle configurations and algorithms, however, the projection uses beam properties determined offline (e.g., through beam tuning or technical commissioning). The best algorithm minimizes either the total uncertainty (i.e., offline and online) or the total offline uncertainty in the projection. Four beam position algorithms are analyzed (A1-A4). Two nozzle lengths are used as examples: a large nozzle (1.5 m length) and a small nozzle (0.4 m length). Three nozzle configurations are considered: IC after SM, IC before SM, and ICs on both sides. Default uncertainties are selected for ion chamber measurements, nozzle entrance beam position and angle, and scanning magnet angle. The results for other uncertainties can be determined by scaling these results or repeating the error propagation. We show the propagation of errors from two locations and the SM angle to the isocenter for all the algorithms. The best choice of algorithm depends on the nozzle length and is A1 and A3 for the large and small nozzles, respectively. If the total offline uncertainty is to be minimized (a better choice if the offline uncertainty is not stable), the best choice of algorithm changes to A1 for the small nozzle for some hardware configurations. Reducing the nozzle length can help to reduce the gantry size and make proton therapy more accessible. This work is important for designing smaller nozzles and, consequently, smaller gantries. This work is also important for log file analyses.

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.

Effects of spot parameters in pencil beam scanning treatment planning.

BACKGROUND: Spot size σ (in air at isocenter), interspot spacing d, and spot charge q influence dose delivery efficiency and plan quality in Intensity Modulated Proton Therapy (IMPT) treatment planning. The choice and range of parameters varies among different manufacturers. The goal of this work is to demonstrate the influence of the spot parameters on dose quality and delivery in IMPT treatment plans, to show their interdependence, and to make practitioners aware of the spot parameter values for a certain facility. Our study could help as a guideline to make the trade-off between treatment quality and time in existing PBS centers and in future systems. METHODS: We created plans for seven patients and a phantom, with different tumor sites and volumes, and compared the effect of small-, medium-, and large-spot widths (σ = 2.5, 5, and 10 mm) and interspot distances (1σ, 1.5σ, and 1.75σ) on dose, spot charge, and treatment time. Moreover, we quantified how postplanning charge threshold cuts affect plan quality and the total number of spots to deliver, for different spot widths and interspot distances. We show the effect of a minimum charge (or MU) cutoff value for a given proton delivery system. RESULTS: Spot size had a strong influence on dose: larger spots resulted in more protons delivered outside the target region. We observed dose differences of 2-13 Gy (RBE) between 2.5 mm and 10 mm spots, where the amount of extra dose was due to dose penumbra around the target region. Interspot distance had little influence on dose quality for our patient group. Both parameters strongly influence spot charge in the plans and thus the possible impact of postplanning charge threshold cuts. If such charge thresholds are not included in the treatment planning system (TPS), it is important that the practitioner validates that a given combination of lower charge threshold, interspot spacing, and spot size does not result in a plan degradation. Low average spot charge occurs for small spots, small interspot distances, many beam directions, and low fractional dose values. CONCLUSIONS: The choice of spot parameters values is a trade-off between accelerator and beam line design, plan quality, and treatment efficiency. We recommend the use of small spot sizes for better organ-at-risk sparing and lateral interspot distances of 1.5σ to avoid long treatment times. We note that plan quality is influenced by the charge cutoff. Our results show that the charge cutoff can be sufficiently large (i.e., 106 protons) to accommodate limitations on beam delivery systems. It is, therefore, not necessary per se to include the charge cutoff in the treatment planning optimization such that Pareto navigation (e.g., as practiced at our institution) is not excluded and optimal plans can be obtained without, perhaps, a bias from the charge cutoff. We recommend that the impact of a minimum charge cut impact is carefully verified for the spot sizes and spot distances applied or that it is accommodated in the TPS.

Impact of spot charge inaccuracies in IMPT treatments.

BACKGROUND: Spot charge is one parameter of pencil-beam scanning dose delivery system whose accuracy is typically high but whose required value has not been investigated. In this work we quantify the dose impact of spot charge inaccuracies on the dose distribution in patients. Knowing the effect of charge errors is relevant for conventional proton machines, as well as for new generation proton machines, where ensuring accurate charge may be challenging. METHODS: Through perturbation of spot charge in treatment plans for seven patients and a phantom, we evaluated the dose impact of absolute (up to 5× 106 protons) and relative (up to 30%) charge errors. We investigated the dependence on beam width by studying scenarios with small, medium and large beam sizes. Treatment plan statistics included the Γ passing rate, dose-volume-histograms and dose differences. RESULTS: The allowable absolute charge error for small spot plans was about 2× 106 protons. Larger limits would be allowed if larger spots were used. For relative errors, the maximum allowable error size for small, medium and large spots was about 13%, 8% and 6% for small, medium and large spots, respectively. CONCLUSIONS: Dose distributions turned out to be surprisingly robust against random spot charge perturbation. Our study suggests that ensuring spot charge errors as small as 1-2% as is commonly aimed at in conventional proton therapy machines, is clinically not strictly needed.

Impact of Spot Size and Beam-Shaping Devices on the Treatment Plan Quality for Pencil Beam Scanning Proton Therapy.

PURPOSE: This study aimed to assess the clinical impact of spot size and the addition of apertures and range compensators on the treatment quality of pencil beam scanning (PBS) proton therapy and to define when PBS could improve on passive scattering proton therapy (PSPT). METHODS AND MATERIALS: The patient cohort included 14 pediatric patients treated with PSPT. Six PBS plans were created and optimized for each patient using 3 spot sizes (∼12-, 5.4-, and 2.5-mm median sigma at isocenter for 90- to 230-MeV range) and adding apertures and compensators to plans with the 2 larger spots. Conformity and homogeneity indices, dose-volume histogram parameters, equivalent uniform dose (EUD), normal tissue complication probability (NTCP), and integral dose were quantified and compared with the respective PSPT plans. RESULTS: The results clearly indicated that PBS with the largest spots does not necessarily offer a dosimetric or clinical advantage over PSPT. With comparable target coverage, the mean dose (Dmean) to healthy organs was on average 6.3% larger than PSPT when using this spot size. However, adding apertures to plans with large spots improved the treatment quality by decreasing the average Dmean and EUD by up to 8.6% and 3.2% of the prescribed dose, respectively. Decreasing the spot size further improved all plans, lowering the average Dmean and EUD by up to 11.6% and 10.9% compared with PSPT, respectively, and eliminated the need for beam-shaping devices. The NTCP decreased with spot size and addition of apertures, with maximum reduction of 5.4% relative to PSPT. CONCLUSIONS: The added benefit of using PBS strongly depends on the delivery configurations. Facilities limited to large spot sizes (>∼8 mm median sigma at isocenter) are recommended to use apertures to reduce treatment-related toxicities, at least for complex and/or small tumors.

PBS machine interlocks using EWMA.

Delivery of pencil beam scanning (PBS) requires the on-line measurement of several beam parameters. If the measurement is outside of specified tolerances and a binary threshold algorithm is used, the beam will be paused. Given instrumentation and statistical noise such a system can lead to many pauses which could increase the treatment time. Statistical quality control methods are typically used on manufacturing lines to monitor a process and give early detection of a gradual problem and stop the process if a deviation is statistically significant. These methods can be used to develop a more intuitive algorithm for (PBS) delivery systems that is robust and safe and leads to decreased treatment times. The Exponentially Weighted Moving Average (EWMA) control scheme monitors deviations in beam properties which are averaged over a specified number of measurements with greater weight applied to the more recent ones. Simulation of an EWMA-style algorithm safely detected shifts in random and systematic delivery errors without false alarms. Binary and EWMA methods can be combined for improved reliability without sacrificing patient safety. In the EWMA method, the mean of a beam property can be related to systematic uncertainties and the standard deviation can be related to random uncertainties. This method allows one to have separate interlock levels for each type of uncertainty and to detect systematic trends.

Shortening delivery times of intensity modulated proton therapy by reducing proton energy layers during treatment plan optimization.

PURPOSE: To shorten delivery times of intensity modulated proton therapy by reducing the number of energy layers in the treatment plan. METHODS AND MATERIALS: We have developed an energy layer reduction method, which was implemented into our in-house-developed multicriteria treatment planning system "Erasmus-iCycle." The method consisted of 2 components: (1) minimizing the logarithm of the total spot weight per energy layer; and (2) iteratively excluding low-weighted energy layers. The method was benchmarked by comparing a robust "time-efficient plan" (with energy layer reduction) with a robust "standard clinical plan" (without energy layer reduction) for 5 oropharyngeal cases and 5 prostate cases. Both plans of each patient had equal robust plan quality, because the worst-case dose parameters of the standard clinical plan were used as dose constraints for the time-efficient plan. Worst-case robust optimization was performed, accounting for setup errors of 3 mm and range errors of 3% + 1 mm. We evaluated the number of energy layers and the expected delivery time per fraction, assuming 30 seconds per beam direction, 10 ms per spot, and 400 Giga-protons per minute. The energy switching time was varied from 0.1 to 5 seconds. RESULTS: The number of energy layers was on average reduced by 45% (range, 30%-56%) for the oropharyngeal cases and by 28% (range, 25%-32%) for the prostate cases. When assuming 1, 2, or 5 seconds energy switching time, the average delivery time was shortened from 3.9 to 3.0 minutes (25%), 6.0 to 4.2 minutes (32%), or 12.3 to 7.7 minutes (38%) for the oropharyngeal cases, and from 3.4 to 2.9 minutes (16%), 5.2 to 4.2 minutes (20%), or 10.6 to 8.0 minutes (24%) for the prostate cases. CONCLUSIONS: Delivery times of intensity modulated proton therapy can be reduced substantially without compromising robust plan quality. Shorter delivery times are likely to reduce treatment uncertainties and costs.

Dose uncertainties in IMPT for oropharyngeal cancer in the presence of anatomical, range, and setup errors.

PURPOSE: Setup, range, and anatomical uncertainties influence the dose delivered with intensity modulated proton therapy (IMPT), but clinical quantification of these errors for oropharyngeal cancer is lacking. We quantified these factors and investigated treatment fidelity, that is, robustness, as influenced by adaptive planning and by applying more beam directions. METHODS AND MATERIALS: We used an in-house treatment planning system with multicriteria optimization of pencil beam energies, directions, and weights to create treatment plans for 3-, 5-, and 7-beam directions for 10 oropharyngeal cancer patients. The dose prescription was a simultaneously integrated boost scheme, prescribing 66 Gy to primary tumor and positive neck levels (clinical target volume-66 Gy; CTV-66 Gy) and 54 Gy to elective neck levels (CTV-54 Gy). Doses were recalculated in 3700 simulations of setup, range, and anatomical uncertainties. Repeat computed tomography (CT) scans were used to evaluate an adaptive planning strategy using nonrigid registration for dose accumulation. RESULTS: For the recalculated 3-beam plans including all treatment uncertainty sources, only 69% (CTV-66 Gy) and 88% (CTV-54 Gy) of the simulations had a dose received by 98% of the target volume (D98%) >95% of the prescription dose. Doses to organs at risk (OARs) showed considerable spread around planned values. Causes for major deviations were mixed. Adaptive planning based on repeat imaging positively affected dose delivery accuracy: in the presence of the other errors, percentages of treatments with D98% >95% increased to 96% (CTV-66 Gy) and 100% (CTV-54 Gy). Plans with more beam directions were not more robust. CONCLUSIONS: For oropharyngeal cancer patients, treatment uncertainties can result in significant differences between planned and delivered IMPT doses. Given the mixed causes for major deviations, we advise repeat diagnostic CT scans during treatment, recalculation of the dose, and if required, adaptive planning to improve adequate IMPT dose delivery.

Intensity modulated proton therapy for postmastectomy radiation of bilateral implant reconstructed breasts: a treatment planning study.

BACKGROUND AND PURPOSE: Delivery of post-mastectomy radiation (PMRT) in women with bilateral implants represents a technical challenge, particularly when attempting to cover regional lymph nodes. Intensity modulated proton therapy (IMPT) holds the potential to improve dose delivery and spare non-target tissues. The purpose of this study was to compare IMPT to three-dimensional (3D) conformal radiation following bilateral mastectomy and reconstruction. MATERIALS AND METHODS: Ten IMPT, 3D conformal photon/electron (P/E), and 3D photon (wide tangent) plans were created for 5 patients with breast cancer, all of whom had bilateral breast implants. Using RTOG guidelines, a physician delineated contours for both target volumes and organs-at-risk. Plans were designed to achieve 95% coverage of all targets (chest wall, IMN, SCV, axilla) to a dose of 50.4 Gy or Gy (RBE) while maximally sparing organs-at-risk. RESULTS: IMPT plans conferred similar target volume coverage with enhanced homogeneity. Both mean heart and lung doses using IMPT were significantly decreased compared to both P/E and wide tangent planning. CONCLUSIONS: IMPT provides improved homogeneity to the chest wall and regional lymphatics in the post-mastectomy setting with improved sparing of surrounding normal structures for woman with reconstructed breasts. IMPT may enable women with mastectomy to undergo radiation therapy without the need for delay in breast reconstruction.

Interpolation of tabulated proton Bragg peaks.

Treatment planning databases for pencil beam scanning can be large, difficult to manage and problematic for quality assurance when they contain tabulated Bragg peaks at small range resolution. Smaller range resolution, in the absence of an accurate interpolation method, improves the accuracy in dose calculations. In this work, we derive an approximate scaling function to interpolate between tabulated Bragg peaks, and determine the accuracy of this interpolation technique and the minimum number of tabulated peaks in a treatment planning database. With the new interpolation technique, three tabulated mono-energetic Bragg peaks (N = 3) are a suitable lower limit for N to achieve interpolation accuracy better than ±1% of the maximum dose in pristine and spread out Bragg peaks for ranges between 6.8 and 32.1 cm of water.

Numerical solutions of the γ-index in two and three dimensions.

The γ-index is used routinely to establish correspondence between two dose distributions. The definition of the γ-index can be written with a single equation but solving this equation at millions of points is computationally expensive, especially in three dimensions. Our goal is to extend the vector-equation method in Bakai et al (2003 Phys. Med. Biol.48 3543-53) to higher order for better accuracy and, as important, to determine the magnitude of accuracy in a higher order solution. We construct a numerical framework for calculating the γ-index in two and three dimensions and present an efficient method for calculating the γ-index with zeroth-, first- and second-order methods using tricubic spline interpolation. For an intensity-modulated radiation therapy example with 1.78 × 106 voxels, the zeroth-order, first-order, first-order iterations and semi-second-order methods calculate the three-dimensional γ-index in 1.5, 4.7, 34.7 and 35.6 s with 36.7%, 1.1%, 0.2% and 0.8% accuracy, respectively. The accuracy of linear interpolation with this example is 1.0%. We present efficient numerical methods for calculating the three-dimensional γ-index with tricubic spline interpolation. The first-order method with iterations is the most accurate and fastest choice of the numerical methods if the dose distributions may have large second-order gradients. Furthermore, the difference between iterations can be used to determine the accuracy of the method.

Monte Carlo study of the potential reduction in out-of-field dose using a patient-specific aperture in pencil beam scanning proton therapy.

This study is aimed at identifying the potential benefits of using a patient-specific aperture in proton beam scanning. For this purpose, an accurate Monte Carlo model of the pencil beam scanning (PBS) proton therapy (PT) treatment head at Massachusetts General Hospital (MGH) was developed based on an existing model of the passive double-scattering (DS) system. The Monte Carlo code specifies the treatment head at MGH with sub-millimeter accuracy. The code was configured based on the results of experimental measurements performed at MGH. This model was then used to compare out-of-field doses in simulated DS treatments and PBS treatments. For the conditions explored, the penumbra in PBS is wider than in DS, leading to higher absorbed doses and equivalent doses adjacent to the primary field edge. For lateral distances greater than 10 cm from the field edge, the doses in PBS appear to be lower than those observed for DS. We found that placing a patient-specific aperture at nozzle exit during PBS treatments can potentially reduce doses lateral to the primary radiation field by over an order of magnitude. In conclusion, using a patient-specific aperture has the potential to further improve the normal tissue sparing capabilities of PBS.

Golden beam data for proton pencil-beam scanning.

Proton, as well as other ion, beams applied by electro-magnetic deflection in pencil-beam scanning (PBS) are minimally perturbed and thus can be quantified a priori by their fundamental interactions in a medium. This a priori quantification permits an optimal reduction of characterizing measurements on a particular PBS delivery system. The combination of a priori quantification and measurements will then suffice to fully describe the physical interactions necessary for treatment planning purposes. We consider, for proton beams, these interactions and derive a 'Golden' beam data set. The Golden beam data set quantifies the pristine Bragg peak depth-dose distribution in terms of primary, multiple Coulomb scatter, and secondary, nuclear scatter, components. The set reduces the required measurements on a PBS delivery system to the measurement of energy spread and initial phase space as a function of energy. The depth doses are described in absolute units of Gy(RBE) mm² Gp⁻¹, where Gp equals 109 (giga) protons, thus providing a direct mapping from treatment planning parameters to integrated beam current. We used these Golden beam data on our PBS delivery systems and demonstrated that they yield absolute dosimetry well within clinical tolerance.

A fast optimization algorithm for multicriteria intensity modulated proton therapy planning.

PURPOSE: To describe a fast projection algorithm for optimizing intensity modulated proton therapy (IMPT) plans and to describe and demonstrate the use of this algorithm in multicriteria IMPT planning. METHODS: The authors develop a projection-based solver for a class of convex optimization problems and apply it to IMPT treatment planning. The speed of the solver permits its use in multicriteria optimization, where several optimizations are performed which span the space of possible treatment plans. The authors describe a plan database generation procedure which is customized to the requirements of the solver. The optimality precision of the solver can be specified by the user. RESULTS: The authors apply the algorithm to three clinical cases: A pancreas case, an esophagus case, and a tumor along the rib cage case. Detailed analysis of the pancreas case shows that the algorithm is orders of magnitude faster than industry-standard general purpose algorithms (MOSEK'S interior point optimizer, primal simplex optimizer, and dual simplex optimizer). Additionally, the projection solver has almost no memory overhead. CONCLUSIONS: The speed and guaranteed accuracy of the algorithm make it suitable for use in multicriteria treatment planning, which requires the computation of several diverse treatment plans. Additionally, given the low memory overhead of the algorithm, the method can be extended to include multiple geometric instances and proton range possibilities, for robust optimization.

A case study in proton pencil-beam scanning delivery.

PURPOSE: We completed an implementation of pencil-beam scanning (PBS), a technology whereby a focused beam of protons, of variable intensity and energy, is scanned over a plane perpendicular to the beam axis and in depth. The aim of radiotherapy is to improve the target to healthy tissue dose differential. We illustrate how PBS achieves this aim in a patient with a bulky tumor. METHODS AND MATERIALS: Our first deployment of PBS uses "broad" pencil-beams ranging from 20 to 35 mm (full-width-half-maximum) over the range interval from 32 to 7 g/cm(2). Such beam-brushes offer a unique opportunity for treating bulky tumors. We present a case study of a large (4,295 cc clinical target volume) retroperitoneal sarcoma treated to 50.4 Gy relative biological effectiveness (RBE) (presurgery) using a course of photons and protons to the clinical target volume and a course of protons to the gross target volume. RESULTS: We describe our system and present the dosimetry for all courses and provide an interdosimetric comparison. DISCUSSION: The use of PBS for bulky targets reduces the complexity of treatment planning and delivery compared with collimated proton fields. In addition, PBS obviates, especially for cases as presented here, the significant cost incurred in the construction of field-specific hardware. PBS offers improved dose distributions, reduced treatment time, and reduced cost of treatment.

Incorporation of the aperture thickness in proton pencil-beam dose calculations.

Field-specific apertures, of sufficient range-absorbing thickness, are used in the majority of proton-therapy treatments today. In current practice, these apertures are modelled as objects of infinitesimal thickness. Such an approximation, however, is not accurate if the aperture edge is close to, or extends over, the beam axis. Practical situations in which this occurs include off-axis patch fields, small apertures, and fields shaped with a multileaf collimator. We develop an extension of the pencil-beam dose model to incorporate the aperture thickness. We derive an exact solution as well as a computationally simpler approximate implementation. The model is validated using measurements of the lateral penumbra. For a set-up with a source size of 2.76 cm, a source-to-axis distance of 227 cm, and a aperture-to-axis distance of 35 cm, the maximum increase in penumbra for a 6 cm thick aperture compared to the thin-aperture model is about 2 mm. The maximum shift in the 95% isodose contour line is larger. The overall effect depends on the aperture thickness, the position of the aperture edge and the intrinsic source size and SAD, but is fairly insensitive to aperture-to-skin distance and depth in patient.

Four-dimensional proton treatment planning for lung tumors.

PURPOSE: In proton radiotherapy, respiration-induced variations in density lead to changes in radiologic path lengths and will possibly result in geometric misses. We compared different treatment planning strategies for lung tumors that compensate for respiratory motion. METHODS AND MATERIALS: Particle-specific treatment planning margins were applied to standard helical computed tomography (CT) scans as well as to "representative" CT scans. Margins were incorporated beam specific laterally by aperture widening and longitudinally by compensator smearing. Furthermore, treatment plans using full time-resolved 4D-computed tomography data were generated. RESULTS: Four-dimensional treatment planning guaranteed target coverage throughout a respiratory cycle. Use of a standard helical CT data set resulted in underdosing the target volume to 36% of the prescribed dose. For CT data representing average target positions, coverage can be expected but not guaranteed. In comparison to this strategy, 4D planning decreased the mean lung dose by up to 16% and the lung volume receiving 20 Gy (prescribed target dose 72 Gy) by up to 15%. CONCLUSION: When the three planning strategies are compared, only 4D proton treatment planning guarantees delivery of the prescribed dose throughout a respiratory cycle. Furthermore, the 4D planning approach results in equal or reduced dose to critical structures; even the ipsilateral lung is spared.

Intra- and interfractional patient motion for a variety of immobilization devices.

The magnitude of inter- and intrafractional patient motion has been assessed for a broad set of immobilization devices. Data was analyzed for the three ordinal directions--left-right (x), sup-inf (y), and ant-post (z)--and the combined spatial displacement. We have defined "rigid" and "non-rigid" immobilization devices depending on whether they could be rigidly and reproducibly connected to the treatment couch or not. The mean spatial displacement for intrafractional motion for rigid devices is 1.3 mm compared to 1.9 mm for nonrigid devices. The modified Gill-Thomas-Cosman frame performed best at controlling intrafractional patient motion, with a 95% probability of observing a three-dimensional (3D) vector length of motion (v95) of less than 1.8 mm, but could not be evaluated for interfractional motion. All other rigid and nonrigid immobilization devices had a v95 of more than 3 mm for intrafractional patient motion. Interfractional patient motion was only evaluated for the rigid devices. The mean total interfractional displacement was at least 3.0 mm for these devices while v95 was at least 6.0 mm.

The prediction of output factors for spread-out proton Bragg peak fields in clinical practice.

The reliable prediction of output factors for spread-out proton Bragg peak (SOBP) fields in clinical practice remained unrealized due to a lack of a consistent theoretical framework and the great number of variables introduced by the mechanical devices necessary for the production of such fields. These limitations necessitated an almost exclusive reliance on manual calibration for individual fields and empirical, ad hoc, models. We recently reported on a theoretical framework for the prediction of output factors for such fields. In this work, we describe the implementation of this framework in our clinical practice. In our practice, we use a treatment delivery nozzle that uses a limited, and constant, set of mechanical devices to produce SOBP fields over the full extent of clinical penetration depths, or ranges, and modulation widths. This use of a limited set of mechanical devices allows us to unfold the physical effects that affect the output factor. We describe these effects and their incorporation into the theoretical framework. We describe the calibration and protocol for SOBP fields, the effects of apertures and range-compensators and the use of output factors in the treatment planning process.

Stereotactic radiotherapy for vestibular schwannomas: favorable outcome with minimal toxicity.

OBJECTIVE: To determine the outcome and toxicity in patients with vestibular schwannomas treated with conventionally fractionated stereotactic radiotherapy (SRT) and to identify prognostic factors that are predictive of outcome. METHODS: Between 1992 and 2001, 70 patients with vestibular schwannomas were treated with linear accelerator-based SRT in our institutions. Eleven patients had neurofibromatosis Type II (NF2). The median age was 53 years (range, 17-82 yrs). The median tumor volume was 2.4 cm3 (range, 0.05-21.1 cm3). The indications for SRT were distributed as follows: 47% newly diagnosed, 31% progressive tumors after watchful waiting, 3% adjuvant postoperative radiation, and 19% recurrent tumors after surgical resection. The median dose was 54 Gy in 1.8 Gy per fraction, prescribed to 95% of the isodose line. Relocatable stereotactic frames were used for daily treatments. The median follow-up was 45.3 months. RESULTS: Tumor recurrence was defined as progressive enlargement of tumor on follow-up magnetic resonance imaging studies. One patient had a tumor recurrence at 38 months after SRT. The actuarial tumor control rates were 100 and 98% at 3 and 5 years, respectively. Three patients with a median tumor volume of 16.2 cm3 required surgical resection for persistent or increasing symptoms at a median of 37 months. The actuarial freedom from resection rates were 98 and 92% at 3 and 5 years, respectively. In multivariate analysis, tumor volume at time of treatment was predictive for neurosurgical intervention (surgical resection or shunt placement) after SRT (P = 0.001). The 3- and 5-year actuarial rates of freedom from any neurosurgical intervention were 100 and 97% for patients with tumor volume less than 8 cm3 and 74 and 47% respectively for patients with tumor of at least 8 cm3 (P < 0.0001). The 3-year actuarial rates of facial and trigeminal nerve preservation were 99 and 96%, respectively. Surgery before SRT was predictive of posttreatment trigeminal neuropathy. The 3-year actuarial rates of freedom from trigeminal neuropathy were 86 and 98% for patients with and without previous resection, respectively (P = 0.04). There was no difference in tumor control and cranial nerve function preservation rates seen in NF2 patients compared with non-NF2 patients. No second primary cancer or malignant transformation was observed. CONCLUSION: SRT in the conventionally fractionated approach results in a very favorable outcome with minimal toxicity, with results comparable to those of the best of the radiosurgery series. Patients with large tumors are more likely to undergo neurosurgical interventions after SRT. Patients who have undergone previous surgery are at increased risk of developing trigeminal neuropathy.