Implementation of a dose calculation algorithm based on Monte Carlo simulations for treatment planning towards MRI guided ion beam therapy

Magnetic resonance guidance in particle therapy has the potential to improve the current performance of clinical workflows. However, the presence of magnetic fields challenges the current algorithms for treatment planning. To ensure proper dose calculations, compensation methods are required to guarantee that the maximum deposited energy of deflected beams lies in the target volume. In addition, proper modifications of the intrinsic dose calculation engines, accounting for magnetic fields, are needed. In this work, an algorithm for proton treatment planning in magnetic fields was implemented in a research treatment planning system (TPS), matRad. Setup-specific look up tables were generated using a validated MC model for a clinical proton beamline (62.4 - 215.7 MeV) interacting with a dipole magnet (B = 0-1 T). The algorithm was successfully benchmarked against MC simulations in water, showing gamma index (2%/2mm) global pass rates higher than 96% for different plan configurations. Additionally, absorbed depth doses were compared with experimental measurements in water. Differences within 2% and 3.5% in the Bragg peak and entrance regions, respectively, were found. Finally, treatment plans were generated and optimized for magnetic field strengths of 0 and 1 T to assess the performance of the proposed model. Equivalent treatment plans and dose volume histograms were achieved, independently of the magnetic field strength. Differences lower than 1.5% for plan quality indicators (D-2%, D-50%, D-90%, V-95% and V-105%) in water, a TG119 phantom and an exemplary prostate patient case were obtained. More complex treatment planning studies are foreseen to establish the limits of applicability of the proposed model.