A pencil beam algorithm for intensity modulated proton therapy derived from Monte Carlo simulations

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.

[1]  A. Lomax,et al.  Intensity modulation methods for proton radiotherapy. , 1999, Physics in medicine and biology.

[2]  C Nauraye,et al.  Experimental determination and verification of the parameters used in a proton pencil beam algorithm. , 2001, Medical physics.

[3]  P. Petti,et al.  Differential-pencil-beam dose calculations for charged particles. , 1992, Medical physics.

[4]  M. Fippel,et al.  A Monte Carlo dose calculation algorithm for proton therapy. , 2004, Medical physics.

[5]  F Nüsslin,et al.  An objective function for radiation treatment optimization based on local biological measures. , 1999, Physics in medicine and biology.

[6]  T Bortfeld,et al.  An analytical approximation of the Bragg curve for therapeutic proton beams. , 1997, Medical physics.

[7]  Michael Lee,et al.  An empirical method to build up a model of proton dose distribution for a radiotherapy treatment-planning package , 1993 .

[8]  U Isacsson,et al.  Implementation of pencil kernel and depth penetration algorithms for treatment planning of proton beams. , 2000, Physics in medicine and biology.

[9]  E Pedroni,et al.  Monte Carlo dose calculations for spot scanned proton therapy , 2005, Physics in medicine and biology.

[10]  E Pedroni,et al.  Experimental characterization and physical modelling of the dose distribution of scanned proton pencil beams , 2005, Physics in medicine and biology.

[11]  R. Jeraj,et al.  The effect of dose calculation accuracy on inverse treatment planning. , 2002, Physics in medicine and biology.

[12]  B. Rossi,et al.  Cosmic-Ray Theory , 1941 .

[13]  P Andreo,et al.  Monte Carlo and analytical calculation of proton pencil beams for computerized treatment plan optimization , 1997, Physics in medicine and biology.

[14]  T Kanai,et al.  Spot scanning system for proton radiotherapy. , 1980, Medical physics.

[15]  Fridtjof Nüsslin,et al.  Hyperion — An integrated IMRT planning tool , 2000 .

[16]  M Goitein,et al.  A pencil beam algorithm for proton dose calculations. , 1996, Physics in medicine and biology.

[17]  I. Kawrakow,et al.  3D electron dose calculation using a Voxel based Monte Carlo algorithm (VMC). , 1996, Medical physics.

[18]  Robert J. Schneider,et al.  Multiple Coulomb scattering of 160 MeV protons , 1993 .

[19]  G Kraft,et al.  Design and construction of a ripple filter for a smoothed depth dose distribution in conformal particle therapy. , 1999, Physics in medicine and biology.

[20]  E. Pedroni,et al.  Dose calculation models for proton treatment planning using a dynamic beam delivery system: an attempt to include density heterogeneity effects in the analytical dose calculation. , 1999, Physics in medicine and biology.

[21]  I. Kawrakow,et al.  Investigation of variance reduction techniques for Monte Carlo photon dose calculation using XVMC , 2000, Physics in medicine and biology.

[22]  Alessandra Bolsi,et al.  Treatment planning and verification of proton therapy using spot scanning: initial experiences. , 2004, Medical physics.

[23]  U. Oelfke,et al.  Two-dimensional pencil beam scaling: an improved proton dose algorithm for heterogeneous media. , 2002, Physics in medicine and biology.

[24]  J O Deasy,et al.  A proton dose calculation algorithm for conformal therapy simulations based on Molière's theory of lateral deflections. , 1998, Medical physics.

[25]  C. Ma,et al.  A particle track-repeating algorithm for proton beam dose calculation , 2005, Physics in medicine and biology.

[26]  M. Fippel.  Fast Monte Carlo dose calculation for photon beams based on the VMC electron algorithm. , 1999, Medical physics.