Validation of the relative biological effectiveness of active-energy scanning carbon-ion radiotherapy on a commercial treatment planning system with a microdosimetic kinetic model

Abstract

Abstract Background The study objective was to validate the relative biological effectiveness (RBE) calculated by the modified microdosimetric kinetic model in RayStation (Ray-MKM) for active-energy scanning carbon-ion radiotherapy. Methods The Ray-MKM was benchmarked using a spread-out Bragg-peak (SOBP) plan, which was suggested in literature from the National Institute of Radiobiological Science (NIRS) in Japan. The residual RBE differences from the MKM at NIRS (NIRS-MKM) were derived using several SOBP plans with different ranges, SOBP widths, and prescriptions. To investigate the origins of the differences, we compared the saturation-corrected dose-mean specific energy $$Z_{1D}^{*}$$ Z 1 D ∗ of the aforementioned SOBPs. Furthermore, we converted the RBE-weighted doses with the Ray-MKM to those with local effect model I (LEM doses). The purpose was to investigate whether the Ray-MKM could reproduce the RBE-weighted conversion study. Results The benchmark determined the value of the clinical dose scaling factor, $$F_{clin}$$ F clin , as 2.40. The target mean RBE deviations between the Ray-MKM and NIRS-MKM were median: 0.6 (minimum: 0.0 to maximum: 1.69) %. The $$Z_{1D}^{*}$$ Z 1 D ∗ difference in-depth led to the RBE difference in-depth and was remarkable at the distal end. The converted LEM doses from the Ray-MKM doses were comparable (the deviation being − 1.8–0.7%) to existing literature. Conclusion This study validated the Ray-MKM based on our active-energy scanning carbon-ion beam via phantom studies. The Ray-MKM could generate similar RBEs as the NIRS-MKM after benchmarking. Analysis based on $$Z_{1D}^{*}$$ Z 1 D ∗ indicated that the different beam qualities and fragment spectra caused the RBE differences. Since the absolute dose differences at the distal end were small, we neglected them. Furthermore, each centre may determine its centre-specific $$F_{clin}$$ F clin based on this approach.