4D Monte Carlo Dose Calculations for Particle Therapy Combined with the Spring Network Model of Lung Motion

To enable the evaluation of the impact of respiratory motion on charged particle therapy and to realize 4D treatment planning while keeping CT exposure as low as possible, we are developing a Monte Carlo dose calculation system combined with a computational biomechanical model of lung motion. The human lung is CT scanned for a single phase of a respiratory motion. The CT images are transformed into tetrahedral elements by automated segmentation. Then, the respiratory motion is simulated using a computational biomechanical model called the spring network model, and the calculated 3D shape of the lung for a given phase is transformed to a voxel data set. For each phase, assuming carbon-ion beam irradiation, biological dose distribution is calculated using the Monte Carlo particle and heavy ion transport code PHITS coupled with a microdosimetric kinetic model. The dose is mapped onto the reference data set to obtain the accumulated dose. The first version of the 4D dose calculation system we have developed so far can realistically reproduce the lung motion, successfully read the data set for each phase and calculate the accumulated dose. The number of the phases to be sampled and their weights can be set arbitrarily, without need of additional CT scanning. Our first simulations for a 70 MeV/u carbon ion beam with a diameter of 2 cm indicate that the dose distribution can significantly change with phase and that many data sets may be needed to accurately evaluate the dose to the surrounding normal tissue. The impact of the system we developed is two-fold: in the short term, it can be used to investigate different issues of 4D treatment planning; our goal is an entirely simulation-based 4D planning from a single CT scanning. The system is being developed further.