Abstract:Purpose Pre-calculation of accurate dose deposition kernels for treatment planning of spot-based radiotherapies, such as Gamma Knife (GK) and Gamma Pod (GP), can be very time-consuming and may require large data storage with an enormous number of possible spots. We proposed a novel kernel decomposition (KD) model to address accurate and fast (real-time) dose calculation with reduced data storage requirements for spot-based treatment planning. The application of the KD model was demonstrated for clinical GK and GP radiotherapy platforms. Methods The dose deposition kernel at each spot (shot position) is modeled as the product of a shift-invariant kernel based on a reference kernel and spatially variant scale factor. The reference kernel, one for each collimator, is defined at the center of the commissioning phantom for GK and at the center of the treatment target for GP and calculated using the Monte Carlo (MC) method. The spatially variant scale factor is defined as the ratio of the mean tissue maximum ratio (TMR) at the candidate shot position to that at the reference kernel position, and the mean TMR map is calculated within the entire volume through parallel beam ray tracing on the density image followed by averaging over all source directions. The proposed KD dose calculations were compared with the MC method and with the GK and GP treatment planning system (TPS) computations for various shot positions and collimator sizes utilizing a phantom and 14 and 12 clinical plans for GK and GP, respectively. Results For the phantom study, the KD Gamma index (3%/1 mm) passing rates were greater than 99% (median 100%) relative to the MC doses, except for the shots close to the boundary. The passing rates dropped below 90% for 8 mm (16 mm) shots positioned within similar to 1 cm (similar to 2 cm) of the boundary. For the clinical GK plans, the KD Gamma passing rates were greater than 99% (median 100%) compared to the MC and greater than 92% (median 99%) compared to the TPS. For the clinical GP plans, the KD Gamma passing rates were greater than 95% (median 98%) compared to the MC and greater than 91% (median 97%) compared to the TPS. The scale factors were calculated in sub-seconds with GPU implementation and only need to be calculated once before treatment plan optimization. The calculation of the dose kernel was also within sub-seconds without requiring beam-by-beam calculation commonly done in the TPS. Conclusion The proposed model can provide an accurate dose and enables real-time dose and derivative calculations by kernel shifting and scaling without pre-calculating or requiring large data storage for GK and GP dose deposition kernels during treatment planning. This model could be useful for spot-based radiotherapy treatment planning by allowing an efficient global fine search for optimal spots.

Stochastic Versus Deterministic Kernel-Based Superposition Approaches For Dose Calculation Of Intensity-Modulated Arcs

Application of dose kernel calculation using a simplified Monte Carlo method to treatment plan for scanned proton beams

SU‐C‐BRC‐03: Development of a Novel Strategy for On‐Demand Monte Carlo and Deterministic Dose Calculation Treatment Planning and Optimization for External Beam Photon and Particle Therapy

A New Monte Carlo-based Treatment Plan Optimization Approach for Intensity Modulated Radiation Therapy

Comparing radiation treatments using intensity-modulated beams, multiple arcs, and single arcs.

Request of Dose Calculation Method in Arc-modulated Radiation Therapy

Leaf-Sequencing For Intensity-Modulated Arc Therapy Using Graph Algorithms

A New Pencil Beam Model for Photon Dose Calculations in Heterogeneous Media

Implementation of convolution/superposition model of photon dose calculation

A Pencil Beam Dose Calculation Model for Cyberknife System

Research on the Convolution/Superposition Algorithm Based on Monte Carlo Calculation in IMRT

SU-E-T-252: Developing a Pencil Beam Dose Calculation Algorithm for CyberKnife System.

Real-time dose computation: GPU-accelerated source modeling and superposition/convolution

Moving Gpu-Opencl-Based Monte Carlo Dose Calculation Toward Clinical Use: Automatic Beam Commissioning and Source Sampling for Treatment Plan Dose Calculation

SU‐E‐T‐175: Clinical Evaluations of Monte Carlo‐Based Inverse Treatment Plan Optimization for Intensity Modulated Radiotherapy

Automatic Commissioning of a GPU-based Monte Carlo Radiation Dose Calculation Code for Photon Radiotherapy

Dose Kernel Decomposition for Spot-Based Radiotherapy Treatment Planning.

A Monte-Carlo-based and GPU-accelerated 4D-dose calculator for a pencil-beam scanning proton therapy system

An efficient dose calculation strategy for intensity modulated proton therapy.

GPU-accelerated Monte Carlo convolution/superposition implementation for dose calculation

A Two-Stage Sequential Linear Programming Approach to Imrt Dose Optimization