Analytic and Monte Carlo calculations of dose-mean lineal energy for 1 MeV-1 GeV protons with application to radiation protection quality factor.

Abstract

Radiation quality for determining biological effects is commonly linked to the microdosimetric quantity lineal energy ( $y$ ) and to the dose-mean lineal energy ( $y_{\text{D}}$ ). Calculations of $y_{\text{D}}$ are typically performed by specialised Monte Carlo track-structure (MCTS) codes, which can be time-intensive. Thus, microdosimetry-based analytic models are potentially useful for practical calculations. Analytic model calculations of proton $y_{\text{D}}$ and radiation protection quality factor ( $Q$ ) values in sub-micron liquid water spheres (diameter 10-1000 nm) over a broad energy range (1 MeV-1 GeV) are compared against MCTS simulations by PHITS, RITRACKS, and Geant4-DNA. Additionally, an improved analytic microdosimetry model is proposed. The original analytic model of Xapsos is refined and model parameters are updated based on Geant4-DNA physics model. Direct proton energy deposition is described by an alternative energy-loss straggling distribution and the contribution of secondary electrons is calculated using the dielectric formulation of the relativistic Born approximation. MCTS simulations of proton $y_{\text{D}}$ values using the latest versions of the PHITS, RITRACKS, and Geant4-DNA are reported along with the Monte Carlo Damage Simulation (MCDS) algorithm. The $y_{\text{D}}$ datasets are then used within the Theory of Dual Radiation Action (TDRA) to illustrate variations in $Q$ with proton energy. By a careful selection of parameters, overall differences at the ~ 10% level between the proposed analytic model and the MCTS codes can be attained, significantly improving upon existing models. MCDS estimates of $y_{\text{D}}$ are generally much lower than estimates from MCTS simulations. The differences of $Q$ among the examined methods are somewhat smaller than those of $y_{\text{D}}$ . Still, estimates of proton $Q$ values by the present model are in better agreement with MCTS-based estimates than the existing analytic models. An improved microdosimetry-based analytic model is presented for calculating proton $y_{\text{D}}$ values over a broad range of proton energies (1 MeV-1 GeV) and target sizes (10-1000 nm) in very good agreement with state-of-the-art MCTS simulations. It is envisioned that the proposed model might be used as an alternative to CPU-intensive MCTS simulations and advance practical microdosimetry and quality factor calculations in medical, accelerator, and space radiation applications.