A semi‐empirical model for the therapeutic range shift estimation caused by inhomogeneities in proton beam therapy

The purpose of this study was to devise a simple semi‐empirical model to estimate the range shift in clinical practices with high‐Z inhomogeneity in proton beam. A semi‐empirical model utilizing the logarithmic dependence on Z in stopping power from Bohr's classical approach has been developed to calculate the range shift due to the presence of inhomogeneity. Range shift from metallic plates of atomic number Z of various thicknesses were measured in water using a parallel plate ionization chamber and calculated with the FLUKA Monte Carlo code. The proton range shifts for bone and polymethyl methacrylate (PMMA) were estimated using the semi‐empirical model and compared with Monte Carlo calculation. The semi‐empirical equation to determine range shift and water equivalent thickness is presented. The model predicts a shift of the therapeutic range to within 2.5% accuracy for initial proton energies of 50 to 250 MeV and atomic numbers from 3.3 (effective Z for water) to 82. This equation is independent of beam energy, and thus provides range shift from high‐Z materials without the knowledge of proton energy. The proposed method of calculating the therapeutic range shift accurately requires only knowledge of the effective or actual atomic number of the inhomogeneity and the thickness of the inhomogeneity along the beam direction. The model generalizes the range shift calculation for any material based on its effective atomic number, and permits reliable prediction of the range shift for material combinations where no data is currently available. The proposed model can be readily implemented in routine clinical practice for proton range shift estimation and quality assurance on the treatment planning. PACS numbers: 87.53.‐j, 87.55.‐x, 87.53.Bn, 87.55.D‐, 87.55.Qr, 87.55.K‐