Monte Carlo vs. Pencil Beam based optimization of stereotactic lung IMRT

<p>Abstract</p> <p>Background</p> <p>The purpose of the present study is to compare finite size pencil beam (fsPB) and Monte Carlo (MC) based optimization of lung intensity-modulated stereotactic radiotherapy (lung IMSRT).</p> <p>Materials and methods</p> <p>A fsPB and a MC algorithm as implemented in a biological IMRT planning system were validated by film measurements in a static lung phantom. Then, they were applied for static lung IMSRT planning based on three different geometrical patient models (one phase static CT, density overwrite one phase static CT, average CT) of the same patient. Both 6 and 15 MV beam energies were used. The resulting treatment plans were compared by how well they fulfilled the prescribed optimization constraints both for the dose distributions calculated on the static patient models and for the accumulated dose, recalculated with MC on each of 8 CTs of a 4DCT set.</p> <p>Results</p> <p>In the phantom measurements, the MC dose engine showed discrepancies < 2%, while the fsPB dose engine showed discrepancies of up to 8% in the presence of lateral electron disequilibrium in the target. In the patient plan optimization, this translates into violations of organ at risk constraints and unpredictable target doses for the fsPB optimized plans. For the 4D MC recalculated dose distribution, MC optimized plans always underestimate the target doses, but the organ at risk doses were comparable. The results depend on the static patient model, and the smallest discrepancy was found for the MC optimized plan on the density overwrite one phase static CT model.</p> <p>Conclusions</p> <p>It is feasible to employ the MC dose engine for optimization of lung IMSRT and the plans are superior to fsPB. Use of static patient models introduces a bias in the MC dose distribution compared to the 4D MC recalculated dose, but this bias is predictable and therefore MC based optimization on static patient models is considered safe.</p

Optimization of extracranial stereotactic radiation therapy of small lung lesions using accurate dose calculation algorithms

BACKGROUND: The aim of this study was to compare and to validate different dose calculation algorithms for the use in radiation therapy of small lung lesions and to optimize the treatment planning using accurate dose calculation algorithms. METHODS: A 9-field conformal treatment plan was generated on an inhomogeneous phantom with lung mimics and a soft tissue equivalent insert, mimicking a lung tumor. The dose distribution was calculated with the Pencil Beam and Collapsed Cone algorithms implemented in Masterplan (Nucletron) and the Monte Carlo system XVMC and validated using Gafchromic EBT films. Differences in dose distribution were evaluated. The plans were then optimized by adding segments to the outer shell of the target in order to increase the dose near the interface to the lung. RESULTS: The Pencil Beam algorithm overestimated the dose by up to 15% compared to the measurements. Collapsed Cone and Monte Carlo predicted the dose more accurately with a maximum difference of -8% and -3% respectively compared to the film. Plan optimization by adding small segments to the peripheral parts of the target, creating a 2-step fluence modulation, allowed to increase target coverage and homogeneity as compared to the uncorrected 9 field plan. CONCLUSION: The use of forward 2-step fluence modulation in radiotherapy of small lung lesions allows the improvement of tumor coverage and dose homogeneity as compared to non-modulated treatment plans and may thus help to increase the local tumor control probability. While the Collapsed Cone algorithm is closer to measurements than the Pencil Beam algorithm, both algorithms are limited at tissue/lung interfaces, leaving Monte-Carlo the most accurate algorithm for dose prediction

Photon beam relative dose validation of the DPM Monte Carlo code in lung‐equivalent media

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/134925/1/mp5671.pd

Evaluation of skin dose associated with different frequencies of bolus applications in post-mastectomy three-dimensional conformal radiotherapy

<p>Abstract</p> <p>Background</p> <p>The study aimed to calculate chest-wall skin dose associated with different frequencies of bolus applications in post-mastectomy three-dimensional conformal radiotherapy (3D-CRT) and to provide detailed information in the selection of an appropriate bolus regimen in this clinical setting.</p> <p>Methods</p> <p>CT-Simulation scans of 22 post-mastectomy patients were used. Chest wall for clinical target volume (CTV) and a volume including 2-mm surface thickness of the chest wall for skin structures were delineated. Precise PLAN 2.11 treatment planning system (TPS) was used for 3D-CRT planning. 50 Gy in 25 fractions were prescribed using tangential fields and 6-MV photons. Six different frequencies of bolus applications (0, 5, 10, 15, 20, and 25) were administered. Cumulative dose-volume histograms were generated for each bolus regimen. The minimum, maximum and mean skin doses associated with the bolus regimens were compared. To test the accuracy of TPS dose calculations, experimental measurements were performed using EBT gafchromic films.</p> <p>Results</p> <p>The mean, minimum and maximum skin doses were significantly increased with increasing days of bolus applications (p < 0.001). The minimum skin doses for 0, 5, 10, 15, 20, and 25 days of bolus applications were 73.0% ± 2.0%, 78.2% ± 2.0%, 83.3% ± 1.7%, 88.3% ± 1.6%, 92.2% ± 1.7%, and 93.8% ± 1.8%, respectively. The minimum skin dose increments between 20 and 25 (1.6% ± 1.0%), and 15 and 20 (4.0% ± 1.0%) days of bolus applications were significantly lower than the dose increments between 0 and 5 (5.2% ± 0.6%), 5 and 10 (5.1% ± 0.8%), and 10 and 15 (4.9% ± 0.8%) days of bolus applications (p < 0.001). The maximum skin doses for 0, 5, 10, 15, 20, and 25 days of bolus applications were 110.1% ± 1.1%, 110.3% ± 1.1%, 110.5% ± 1.2%, 110.8% ± 1.3%, 111.2% ± 1.5%, and 112.2% ± 1.7%, respectively. The maximum skin dose increments between 20 and 25 (1.0% ± 0.6%), and 15 and 20 (0.4% ± 0.3%) days of bolus applications were significantly higher than the dose increments between 0 and 5 (0.2% ± 0.2%), 5 and 10 (0.2% ± 0.2%), and 10 and 15 (0.2% ± 0.2%) days of bolus applications (p ≤ 0.003). The TPS overestimated the near-surface dose 10.8% at 2-mm below the skin surface.</p> <p>Conclusion</p> <p>In post-mastectomy 3D-CRT, using a 1-cm thick bolus in up to 15 of the total 25 fractions increased minimum skin doses with a tolerable increase in maximum doses.</p

eIMRT: a web platform for the verification and optimization of radiation treatment plans

The eIMRT platform is a remote distributed computing tool that provides users with Internet access to three different services: Monte Carlo optimization of treatment plans, CRT & IMRT treatment optimization, and a database of relevant radiation treatments/clinical cases. These services are accessible through a user-friendly and platform independent web page. Its flexible and scalable design focuses on providing the final users with services rather than a collection of software pieces. All input and output data (CT, contours, treatment plans and dose distributions) are handled using the DICOM format. The design, implementation, and support of the verification and optimization algorithms are hidden to the user. This allows a unified, robust handling of the software and hardware that enables these computation-intensive services. The eIMRT platform is currently hosted by the Galician Supercomputing Center (CESGA) and may be accessible upon request (there is a demo version at http://eimrt.cesga.es:8080/ eIMRT2/demo; request access in http://eimrt.cesga.es/signup.html). This paper describes all aspects of the eIMRT algorithms in depth, its user interface, and its services. Due to the flexible design of the platform, it has numerous applications including the intercenter comparison of treatment planning, the quality assurance of radiation treatments, the design and implementation of new approaches to certain types of treatments, and the sharing of information on radiation treatment techniques. In addition, the web platform and software tools developed for treatment verification and optimization have a modular design that allows the user to extend them with new algorithms. This software is not a commercial product. It is the result of the collaborative effort of different public research institutions and is planned to be distributed as an open source project. In this way, it will be available to any user; new releases will be generated with the new implemented codes or upgradesThis work was financed by Xunta de Galicia of Spain through grant PGIDT05SIN00101CT and by the European Community through the BeInGrid projectS

First steps towards a fast-neutron therapy planning program

<p>Abstract</p> <p>Background</p> <p>The Monte Carlo code GEANT4 was used to implement first steps towards a treatment planning program for fast-neutron therapy at the FRM II research reactor in Garching, Germany. Depth dose curves were calculated inside a water phantom using measured primary neutron and simulated primary photon spectra and compared with depth dose curves measured earlier. The calculations were performed with GEANT4 in two different ways, simulating a simple box geometry and splitting this box into millions of small voxels (this was done to validate the voxelisation procedure that was also used to voxelise the human body).</p> <p>Results</p> <p>In both cases, the dose distributions were very similar to those measured in the water phantom, up to a depth of 30 cm. In order to model the situation of patients treated at the FRM II MEDAPP therapy beamline for salivary gland tumors, a human voxel phantom was implemented in GEANT4 and irradiated with the implemented MEDAPP neutron and photon spectra. The 3D dose distribution calculated inside the head of the phantom was similar to the depth dose curves in the water phantom, with some differences that are explained by differences in elementary composition. The lateral dose distribution was studied at various depths. The calculated cumulative dose volume histograms for the voxel phantom show the exposure of organs at risk surrounding the tumor.</p> <p>Conclusions</p> <p>In order to minimize the dose to healthy tissue, a conformal treatment is necessary. This can only be accomplished with the help of an advanced treatment planning system like the one developed here. Although all calculations were done for absorbed dose only, any biological dose weighting can be implemented easily, to take into account the increased radiobiological effectiveness of neutrons compared to photons.</p

Theoretical methods for the calculation of Bragg curves and 3D distributions of proton beams

The well-known Bragg-Kleeman rule RCSDA = A dot E0p has become a pioneer work in radiation physics of charged particles and is still a useful tool to estimate the range RCSDA of approximately monoenergetic protons with initial energy E0 in a homogeneous medium. The rule is based on the continuous-slowing-down-approximation (CSDA). It results from a generalized (nonrelativistic) Langevin equation and a modification of the phenomenological friction term. The complete integration of this equation provides information about the residual energy E(z) and dE(z)/dz at each position z (0 <= z <= RCSDA). A relativistic extension of the generalized Langevin equation yields the formula RCSDA = A dot (E0 +E02/2M dot c2)p. The initial energy of therapeutic protons satisfies E0 << 2M dot c2 (M dot c2 = 938.276 MeV), which enables us to consider the relativistic contributions as correction terms. Besides this phenomenological starting-point, a complete integration of the Bethe-Bloch equation (BBE) is developed, which also provides the determination of RCSDA, E(z) and dE(z)/dz and uses only those parameters given by the BBE itself (i.e., without further empirical parameters like modification of friction). The results obtained in the context of the aforementioned methods are compared with Monte-Carlo calculations (GEANT4); this Monte-Carlo code is also used with regard to further topics such as lateral scatter, nuclear interactions, and buildup effects. In the framework of the CSDA, the energy transfer from protons to environmental atomic electrons does not account for local fluctuations.Comment: 97 pages review pape