Monte Carlo Simulations of Particle Interactions with Tissue in Carbon Ion Therapy

MONTE-CARLO SIMULATION USING PHITS OF SECONDARY NEUTRONS PRODUCED IN-PATIENT DURING 16O ION THERAPY

In hadrontherapy, oxygen ions 16O can be currently considered as an alternative to carbon ions 12C designed specifically for the treatment of deep and radioresistant tumors. Secondary particles, particularly neutrons constitute a serious problem of undesirable additional irradiation to surrounding healthy tissue. The objective of this study is to evaluate, by Monte-Carlo simulation [code Particle and Heavy Ion Transport code System (PHITS)], the contribution in terms of dose of secondary neutrons produced during interaction 16O ion of 300 MeV u-1 in a soft tissue phantom. The dose of 16O ion, secondary particles and neutrons is evaluated, as well as the particle fluence and energy spectra of neutrons. The contribution to the total dose of secondary neutrons in a soft tissue phantom represents 0.1%. This dose, although apparently insignificant, is essential to conduct even more in-depth studies to understand the long-term effects of these secondary neutrons on the patient's body especially in pediatric case.

Dose calculation accuracy in particle therapy: Comparing carbon ions with protons

Purpose

This work presents the validation of an analytical pencil beam dose calculation algorithm in a commercial treatment planning system (TPS) for carbon ions by measurements of dose distributions in heterogeneous phantom geometries. Additionally, a comparison study of carbon ions versus protons is performed considering current best solutions in commercial TPS.

Methods

All treatment plans were optimized and calculated using the RayStation TPS (RaySearch, Sweden). The dose distributions calculated with the TPS were compared with measurements using a 24‐pinpoint ionization chamber array (T31015, PTW, Germany). Tissue‐like inhomogeneities (bone, lung, and soft tissue) were embedded in water, while a target volume of 4 x 4 x 4 cm³ was defined at two different depths behind the heterogeneities. In total, 10 different test cases, with and without range shifter as well as different air gaps, were investigated. Dose distributions inside as well as behind the target volume were evaluated.

Results

Inside the target volume, the mean dose difference between calculations and measurements, averaged over all test cases, was 1.6% for carbon ions. This compares well to the final agreement of 1.5% obtained in water at the commissioning stage of the TPS for carbon ions and is also within the clinically acceptable interval of 3%. The mean dose difference and maximal dose difference obtained outside the target area were 1.8% and 13.4%, respectively. The agreement of dose distributions for carbon ions in the target volumes was comparable or better to that between Monte Carlo (MC) dose calculations and measurements for protons. Percentage dose differences of more than 10% were present outside the target area behind bone–lung structures, where the carbon ion calculations systematically over predicted the dose. MC dose calculations for protons were superior to carbon ion beams outside the target volumes.

Conclusion

The pencil beam dose calculations for carbon ions in RayStation were found to be in good agreement with dosimetric measurements in heterogeneous geometries for points of interest located within the target. Large local discrepancies behind the target may contribute to incorrect dose predictions for organs at risk.

Automatic phase space generation for Monte Carlo calculations of intensity modulated particle therapy

Monte Carlo (MC) is generally considered as the most accurate dose calculation tool for particle therapy. However, a proper description of the beam particle kinematics is a necessary input for a realistic simulation. Such a description can be stored in phase space (PS) files for different beam energies. A PS file contains kinetic information such as energies, positions and travelling directions for particles traversing a plane perpendicular to the beam direction. The accuracy of PS files plays a critical role in the performance of the MC method for dose calculations. A PS file can be generated with a set of parameters describing analytically the beam kinematics. However, determining such parameters can be tedious and time consuming. Thus, we have developed an algorithm to obtain those parameters automatically and efficiently. In this paper, we presented such an algorithm and compared dose calculations using PS automatically generated for the Shanghai Proton and Heavy Ion Center (SPHIC) with measurements. The gamma-index for comparing calculated depth dose distributions (DDD) with measurements are above 96.0% with criterion 0.6%/0.6 mm. For each single energy, the mean difference percentage between calculated lateral spot sizes at 5 different locations along beam direction and measurements are below 3.5%.