Dose optimization for intensity modulated radiotherapy (IMRT) using small field elements (beamlets) requires the computation of a large number of very small, often only virtual fields of typically a few mm to 1 cm in size. The primary requirements for a suitable dose computation algorithm are (1) speed and (2) proper consideration of the penumbra of the fields which are composed of these beamlets. Here, a finite size pencil beam (fsPB) algorithm is proposed which was specifically designed for the purpose of beamlet-based IMRT. The algorithm employs an analytical function for the cross-profiles of the beamlets which is based on the assumption of self-consistency, i.e. the requirement that an arbitrary superposition of abutting beamlets should add up to a homogeneous field. The depth dependence is stored in tables derived from Monte Carlo computed dose distributions. It is demonstrated that the algorithm produces accurately the output factors and cross-profiles of typical multi-leaf-shaped segments. Due to the accurate penumbra model, the dose distribution features physically feasible gradients at any stage of the iterative optimization, which eliminates the problem of large discrepancies in normal tissue dose due to misaligned gradients between optimized and recomputed treatment plans.

No abstract

For beamlet-based IMRT optimization, fast and less accurate dose computation algorithms are frequently used, while more accurate algorithms are needed to recompute the final dose for verification. In order to speed up the optimization process and ensure close proximity between dose in optimization and verification, proper consideration of dose gradients and tissue inhomogeneity effects should be ensured at every stage of the optimization. Due to their speed, pencil beam algorithms are often used for precalculation of beamlet dose distributions in IMRT treatment planning systems. However, accounting for tissue heterogeneities with these models requires the use of approximate rescaling methods. Recently, a finite size pencil beam (fsPB) algorithm, based on a simple and small set of data, was proposed which was specifically designed for the purpose of dose pre-computation in beamlet-based IMRT. The present work describes the incorporation of 3D density corrections, based on Monte Carlo simulations in heterogeneous phantoms, into this method improving the algorithm accuracy in inhomogeneous geometries while keeping its original speed and simplicity of commissioning. The algorithm affords the full accuracy of 3D density corrections at every stage of the optimization, hence providing the means for density related fluence modulation like penumbra shaping at field edges.

In this ROCOCO in silico trial, a reduction in mean dose to OARs was achieved using particle therapy compared to photons in the re-irradiation of HNSCC. There was a dosimetric benefit favouring carbon-ions above proton therapy. These dose reductions may potentially translate into lower severe complication rates related to the re-irradiation.