Matsuya, Yusuke

Abstract Objective . Estimation of the probability density of the microdosimetric quantities in macroscopic matter is indispensable for applying the concept of microdosimetry to medical physics and radiological protection. The Particle and Heavy Ion Transport code System (PHITS) enables estimating the microdosimetric probability densities due to its unique hybrid modality between the Monte Carlo and analytical approaches called the microdosimetric function. It can convert the deposition energies calculated by the macroscopic Monte Carlo radiation transport simulation to microdosimetric probability densities in water using an analytical function based on the track-structure simulations. Approach . In this study, we improved this function using the latest track-structure simulation codes implemented in PHITS. The improved function is capable of calculating the probability densities of not only the conventional microdosimetric quantities such as lineal energy but also the number of ionization events occurring in a target site, the so-called ionization cluster size distribution, for arbitrary site diameters from 3 nm to 1 μ m. Main results . The accuracy of the improved function was well verified by comparing the microdosimetric probability densities measured by tissue-equivalent proportional counters with the corresponding data calculated in this study. Test calculations for clonogenic cell survival using the improved function coupled with the modified microdosimetric kinetic model suggested a slight increase of its relative biological effectiveness compared with our previous estimations. As a new application of the improved function, we calculated the relative biological effectiveness of the single-strand break and double-strand break yields for proton irradiations using the updated PHITS coupled with the simplified DNA damage estimation model, and confirmed its equivalence in accuracy and its superiority in computational time compared to our previously proposed method based on the track-structure simulation. Significance . From these features, we concluded that the improved function could expand the application fields of PHITS by bridging the gap between microdosimetry and macrodosimetry.

Abstract Proton beam therapy allows irradiating tumor volumes with reduced side effects on normal tissues with respect to conventional x-ray radiotherapy. Biological effects such as cell killing after proton beam irradiations depend on the proton kinetic energy, which is intrinsically related to early DNA damage induction. As such, DNA damage estimation based on Monte Carlo simulations is a research topic of worldwide interest. Such simulation is a mean of investigating the mechanisms of DNA strand break formations. However, past modellings considering chemical processes and DNA structures require long calculation times. Particle and heavy ion transport system (PHITS) is one of the general-purpose Monte Carlo codes that can simulate track structure of protons, meanwhile cannot handle radical dynamics simulation in liquid water. It also includes a simple model enabling the efficient estimation of DNA damage yields only from the spatial distribution of ionizations and excitations without DNA geometry, which was originally developed for electron track-structure simulations. In this study, we investigated the potential application of the model to protons without any modification. The yields of single-strand breaks, double-strand breaks (DSBs) and the complex DSBs were assessed as functions of the proton kinetic energy. The PHITS-based estimation showed that the DSB yields increased as the linear energy transfer (LET) increased, and reproduced the experimental and simulated yields of various DNA damage types induced by protons with LET up to about 30 keV μ m −1 . These results suggest that the current DNA damage model implemented in PHITS is sufficient for estimating DNA lesion yields induced after protons irradiation except at very low energies (below 1 MeV). This model contributes to evaluating early biological impacts in radiation therapy.