Four-dimensional radiotherapeutic dose calculation using biomechanical respiratory motion description

Organ motion due to patient breathing introduces a technical challenge for dosimetry and lung tumor treatment by hadron therapy. Accurate dose distribution estimation requires patient-specific information on tumor position, size, and shape as well as information regarding the material density and stopping power of the media along the beam path. A new 4D dosimetry method was developed, which can be coupled to any motion estimation method. As an illustration, the new method was implemented and tested with a biomechanical model and clinical data.

First, an anatomical model of the lung and tumor was synthesized with deformable tetrahedral grids using computed tomography (CT) images. The CT attenuation values were estimated at the grid vertices. Respiratory motion was simulated biomechanically based on nonlinear finite element analysis. Contrary to classical image-based methods where motion is described using deformable image registration algorithms, the dose distribution was accumulated over tetrahedral meshes that are deformed using biomechanical modeling based on finite element analysis.

The new method preserves the mass of the objects during simulation with an error between 1.6 and 3.6%. The new method was compared to an existing dose calculation method demonstrating significant differences between the two approaches and overall superior performance using the new method.

A unified model of 4D radiotherapy respiratory effects was developed where biomechanical simulations are coupled with dose calculations. Promising results demonstrate that this approach has significant potential for the treatment for moving tumors.