The microdosimetric kinetic model (MKM) is widely used for estimating relative biological effectiveness (RBE)-weighted doses for various radiotherapies because it can determine the surviving fraction of irradiated cells based on only the lineal energy distribution, and it is independent of the radiation type and ion species. However, the applicability of the method to proton therapy has not yet been investigated thoroughly. In this study, we validated the RBE-weighted dose calculated by the MKM in tandem with the Monte Carlo code PHITS for proton therapy by considering the complete simulation geometry of the clinical proton beam line. The physical dose, lineal energy distribution, and RBE-weighted dose for a 155 MeV mono-energetic and spread-out Bragg peak (SOBP) beam of 60 mm width were evaluated. In estimating the physical dose, the calculated depth dose distribution by irradiating the mono-energetic beam using PHITS was consistent with the data measured by a diode detector. A maximum difference of 3.1% in the depth distribution was observed for the SOBP beam. In the RBE-weighted dose validation, the calculated lineal energy distributions generally agreed well with the published measurement data. The calculated and measured RBE-weighted doses were in excellent agreement, except at the Bragg peak region of the mono-energetic beam, where the calculation overestimated the measured data by ~15%. This research has provided a computational microdosimetric approach based on a combination of PHITS and MKM for typical clinical proton beams. The developed RBE-estimator function has potential application in the treatment planning system for various radiotherapies.

In proton therapy, the Bragg peak of a proton beam reportedly deteriorates when passing though heterogeneous structures such as human lungs. Previous studies have used heterogeneous random voxel phantoms, in which soft tissues and air are randomly allotted to render the phantoms the same density as human lungs, for conducting Monte Carlo (MC) simulations. However, measurements of these phantoms are complicated owing to their difficult‐to‐manufacture shape. In the present study, we used Voronoi tessellation to design a phantom that can be manufactured, and prepared a Voronoi lung phantom for which both measurement and MC calculations are possible. Our aim was to evaluate the effectiveness of this phantom as a new lung phantom for investigating proton beam Bragg peak deterioration. For this purpose, we measured and calculated the percentage depth dose and the distal falloff widths (DFW) passing through the phantom. For the 155 MeV beam, the measured and calculated DFW values with the Voronoi lung phantom were 0.40 and 0.39 cm, respectively. For the 200 MeV beam, the measured and calculated DFW values with the Voronoi lung phantom were both 0.48 cm. Our results indicate that both the measurements and MC calculations exhibited high reproducibility with plastinated lung sample from human body in previous studies. We found that better results were obtained using the Voronoi lung phantom than using other previous phantoms. The designed phantom may contribute significantly to the improvement of measurement precision. This study suggests that the Voronoi lung phantom is useful for simulating the effects of the heterogeneous structure of lungs on proton beam deterioration.