Treatment plans optimized for intensity modulated proton therapy (IMPT) may be very sensitive to setup errors and range uncertainties. If these errors are not accounted for during treatment planning, the dose distribution realized in the patient may by strongly degraded compared to the planned dose distribution. The authors implemented the probabilistic approach to incorporate uncertainties directly into the optimization of an intensity modulated treatment plan. Following this approach, the dose distribution depends on a set of random variables which parameterize the uncertainty, as does the objective function used to optimize the treatment plan. The authors optimize the expected value of the objective function. They investigate IMPT treatment planning regarding range uncertainties and setup errors. They demonstrate that incorporating these uncertainties into the optimization yields qualitatively different treatment plans compared to conventional plans which do not account for uncertainty. The sensitivity of an IMPT plan depends on the dose contributions of individual beam directions. Roughly speaking, steep dose gradients in beam direction make treatment plans sensitive to range errors. Steep lateral dose gradients make plans sensitive to setup errors. More robust treatment plans are obtained by redistributing dose among different beam directions. This can be achieved by the probabilistic approach. In contrast, the safety margin approach as widely applied in photon therapy fails in IMPT and is neither suitable for handling range variations nor setup errors.

A pencil beam algorithm as a component of an optimization algorithm for intensity modulated proton therapy (IMPT) is presented. The pencil beam algorithm is tuned to the special accuracy requirements of IMPT, where in heterogeneous geometries both the position and distortion of the Bragg peak and the lateral scatter pose problems which are amplified by the spot weight optimization. Heterogeneity corrections are implemented by a multiple raytracing approach using fluence-weighted sub-spots. In order to derive nuclear interaction corrections, Monte Carlo simulations were performed. The contribution of long ranged products of nuclear interactions is taken into account by a fit to the Monte Carlo results. Energy-dependent stopping power ratios are also implemented. Scatter in optional beam line accessories such as range shifters or ripple filters is taken into account. The collimator can also be included, but without additional scattering. Finally, dose distributions are benchmarked against Monte Carlo simulations, showing 3%/1 mm agreement for simple heterogeneous phantoms. In the case of more complicated phantoms, principal shortcomings of pencil beam algorithms are evident. The influence of these effects on IMPT dose distributions is shown in clinical examples.

A Monte Carlo (MC) code (VMCpro) for treatment planning in proton beam therapy of cancer is introduced. It is based on ideas of the Voxel Monte Carlo algorithm for photons and electrons and is applicable to human tissue for clinical proton energies. In the present paper the implementation of electromagnetic and nuclear interactions is described. They are modeled by a Class II condensed history algorithm with continuous energy loss, ionization, multiple scattering, range straggling, delta-electron transport, nuclear elastic proton nucleus scattering and inelastic proton nucleus reactions. VMCpro is faster than the general purpose MC codes FLUKA by a factor of 13 and GEANT4 by a factor of 35 for simulations in a phantom with inhomogeneities. For dose calculations in patients the speed improvement is larger, because VMCpro has only a weak dependency on the heterogeneity of the calculation grid. Dose distributions produced with VMCpro are in agreement with GEANT4 results. Integrated or broad beam depth dose curves show maximum deviations not larger than 1% or 0.5 mm in regions with large dose gradients for the examples presented here.

ABSTRACT

The conversion of computed tomography (CT) numbers into material composition and mass density data influences the accuracy of patient dose calculations in Monte Carlo treatment planning (MCTP). The aim of our work was to develop a CT conversion scheme by performing a stoichiometric CT calibration. Fourteen dosimetrically equivalent tissue subsets (bins), of which ten bone bins, were created. After validating the proposed CT conversion scheme on phantoms, it was compared to a conventional five bin scheme with only one bone bin. This resulted in dose distributions D(14) and D(5) for nine clinical patient cases in a European multi-centre study. The observed local relative differences in dose to medium were mostly smaller than 5%. The dose-volume histograms of both targets and organs at risk were comparable, although within bony structures D(14) was found to be slightly but systematically higher than D(5). Converting dose to medium to dose to water (D(14) to D(14wat) and D(5) to D(5wat)) resulted in larger local differences as D(5wat) became up to 10% higher than D(14wat). In conclusion, multiple bone bins need to be introduced when Monte Carlo (MC) calculations of patient dose distributions are converted to dose to water.

No abstract

This article reports on a 4D-treatment planning workshop (4DTPW), held on 7-8 December 2009 at the Paul Scherrer Institut (PSI) in Villigen, Switzerland. The participants were all members of institutions actively involved in particle therapy delivery and research. The purpose of the 4DTPW was to discuss current approaches, challenges, and future research directions in 4D-treatment planning in the context of actively scanned particle radiotherapy. Key aspects were addressed in plenary sessions, in which leaders of the field summarized the state-of-the-art. Each plenary session was followed by an extensive discussion. As a result, this article presents a summary of recommendations for the treatment of mobile targets (intrafractional changes) with actively scanned particles and a list of requirements to elaborate and apply these guidelines clinically.

No abstract