A benchmarking method to evaluate the accuracy of a commercial proton monte carlo pencil beam scanning treatment planning system.

AcurosPT is a Monte Carlo algorithm in the Eclipse 13.7 treatment planning system, which is designed to provide rapid and accurate dose calculations for proton therapy. Computational run-time in minimized by simplifying or eliminating less significant physics processes. In this article, the accuracy of AcurosPT was benchmarked against both measurement and an independent MC calculation, TOPAS. Such a method can be applied to any new MC calculation for the detection of potential inaccuracies. To validate multiple Coulomb scattering (MCS) which affects primary beam broadening, single spot profiles in a Solidwater® phantom were compared for beams of five selected proton energies between AcurosPT, measurement and TOPAS. The spot Gaussian sigma in AcurosPT was found to increase faster with depth than both measurement and TOPAS, suggesting that the MCS algorithm in AcurosPT overestimates the scattering effect. To validate AcurosPT modeling of the halo component beyond primary beam broadening, field size factors (FSF) were compared for multi-spot profiles measured in a water phantom. The FSF for small field sizes were found to disagree with measurement, with the disagreement increasing with depth. Conversely, TOPAS simulations of the same FSF consistently agreed with measurement to within 1.5%. The disagreement in absolute dose between AcurosPT and measurement was smaller than 2% at the mid-range depth of multi-energy beams. While AcurosPT calculates acceptable dose distributions for typical clinical beams, users are cautioned of potentially larger errors at distal depths due to overestimated MCS and halo implementation.

Selective robust optimization: A new intensity-modulated proton therapy optimization strategy.

PURPOSE: To develop a new robust optimization strategy for intensity-modulated proton therapy as an important step in translating robust proton treatment planning from research to clinical applications. METHODS: In selective robust optimization, a worst-case-based robust optimization algorithm is extended, and terms of the objective function are selectively computed from either the worst-case dose or the nominal dose. Two lung cancer cases and one head and neck cancer case were used to demonstrate the practical significance of the proposed robust planning strategy. The lung cancer cases had minimal tumor motion less than 5 mm, and, for the demonstration of the methodology, are assumed to be static. RESULTS: Selective robust optimization achieved robust clinical target volume (CTV) coverage and at the same time increased nominal planning target volume coverage to 95.8%, compared to the 84.6% coverage achieved with CTV-based robust optimization in one of the lung cases. In the other lung case, the maximum dose in selective robust optimization was lowered from a dose of 131.3% in the CTV-based robust optimization to 113.6%. Selective robust optimization provided robust CTV coverage in the head and neck case, and at the same time improved controls over isodose distribution so that clinical requirements may be readily met. CONCLUSIONS: Selective robust optimization may provide the flexibility and capability necessary for meeting various clinical requirements in addition to achieving the required plan robustness in practical proton treatment planning settings.

Testing of the analytical anisotropic algorithm for photon dose calculation.

The analytical anisotropic algorithm (AAA) was implemented in the Eclipse (Varian Medical Systems) treatment planning system to replace the single pencil beam (SPB) algorithm for the calculation of dose distributions for photon beams. AAA was developed to improve the dose calculation accuracy, especially in heterogeneous media. The total dose deposition is calculated as the superposition of the dose deposited by two photon sources (primary and secondary) and by an electron contamination source. The photon dose is calculated as a three-dimensional convolution of Monte-Carlo precalculated scatter kernels, scaled according to the electron density matrix. For the configuration of AAA, an optimization algorithm determines the parameters characterizing the multiple source model by optimizing the agreement between the calculated and measured depth dose curves and profiles for the basic beam data. We have combined the acceptance tests obtained in three different departments for 6, 15, and 18 MV photon beams. The accuracy of AAA was tested for different field sizes (symmetric and asymmetric) for open fields, wedged fields, and static and dynamic multileaf collimation fields. Depth dose behavior at different source-to-phantom distances was investigated. Measurements were performed on homogeneous, water equivalent phantoms, on simple phantoms containing cork inhomogeneities, and on the thorax of an anthropomorphic phantom. Comparisons were made among measurements, AAA, and SPB calculations. The optimization procedure for the configuration of the algorithm was successful in reproducing the basic beam data with an overall accuracy of 3%, 1 mm in the build-up region, and 1%, 1 mm elsewhere. Testing of the algorithm in more clinical setups showed comparable results for depth dose curves, profiles, and monitor units of symmetric open and wedged beams below dmax. The electron contamination model was found to be suboptimal to model the dose around dmax, especially for physical wedges at smaller source to phantom distances. For the asymmetric field verification, absolute dose difference of up to 4% were observed for the most extreme asymmetries. Compared to the SPB, the penumbra modeling is considerably improved (1%, 1 mm). At the interface between solid water and cork, profiles show a better agreement with AAA. Depth dose curves in the cork are substantially better with AAA than with SPB. Improvements are more pronounced for 18 MV than for 6 MV. Point dose measurements in the thoracic phantom are mostly within 5%. In general, we can conclude that, compared to SPB, AAA improves the accuracy of dose calculations. Particular progress was made with respect to the penumbra and low dose regions. In heterogeneous materials, improvements are substantial and more pronounced for high (18 MV) than for low (6 MV) energies.