Reproducing the NIRS-QST Clinical Dose Calculations for Carbon Ion Radiotherapy Using Microdosimetric Probability Density Distributions

Ion radiotherapy requires accurate relative biological effectiveness (RBE) calculations to account for the markedly different biological effects of ions compared to photons. Microdosimetric RBE models rely on descriptions of the energy deposition at the microscopic scale, either through radial dose distributions (RDDs) or microdosimetric probability density distributions. While RDD approaches focus on the theoretical description of the energy deposition around the ion track, microdosimetric distributions offer the advantage of being experimentally measurable, which is crucial for quality assurance programs. As the results of microdosimetric RBE models depend on whether RDD or microdosimetric distributions are used, the model parameters are not interchangeable between these approaches. This study presents and validates a method to reproduce the published reference biological and clinical dose calculations at NIRS-QST for only carbon ion radiotherapy by using the modified microdosimetric kinetic model (MKM) alongside microdosimetric distributions instead of the reference RDD approach. To achieve this, Monte Carlo simulations were performed to estimate the variation of the radiation quality within and outside the field of pristine and spread-out Bragg peaks. By appropriately optimizing the modified MKM parameters for microdosimetric distributions assessed within water spheres, we successfully reproduced the results of calculations using the reference NIRS-QST RDD, generally within 2%.