Dosimetric evaluation and benchmarking of proton beam delivery methods in heterogeneous phantoms using the PHITS Monte Carlo code.

Proton therapy is a precise radiation technology; this study focuses on simulating dose distributions in pelvic and thoracic geometries representative of complex clinical scenarios. Two prominent delivery methods, Pencil Beam (PB) and Uniform Beam (UB), are employed to maximize damage to target cancer cells while minimizing harm to surrounding healthy tissues. This work proposes an efficient Monte Carlo simulation approach using PHITS to model energy deposition in anatomical regions involving high-density (bone) and low-density (lung) tissue interfaces. The simulation was designed to deliver clinically prescribed doses of 30 Gy (RBE) for pelvic model and 60 Gy (RBE) for Lung Cancer. To target volume at specific depths, beam energies were optimized to 100-120 MeV (PB) and 101-119 MeV (UB) for high-density pelvic region, and 78-100 MeV (PB) and 80-96 MeV (UB) for Lung Cancer, utilizing a global simulation factor of 108 proton histories. In the investigation of the pelvic (splenic) phantom model (18 fractions), the calculated doses for PB and UB were 29.81, 23.56, and 24.41 Gy (W) and 29.88, 26.07, and 24.58 Gy (W) in the GTV, CTV, and PTV regions, respectively. The absorbed doses for bone and skin were lower with the Pencil Beam (17.95 and 14.86 Gy) compared to the Uniform Beam (19.42 and 16.02 Gy). In the case of Lung Cancer (30 fractions), we estimated doses of 59.91, 42.34, and 26.33 Gy (W) for PB and 59.94, 45.73, and 29.28 Gy (W) for UB in the GTV, CTV, and PTV regions, respectively. Regarding Organs at Risk (OARs), the lung, bone, skin, and heart received 21.71, 22.44, 21.71, and 0.03 Gy (W) using the Pencil Beam, whereas the Uniform Beam resulted in higher OAR exposure of 23.32, 24.75, 25.11, and 0.04 Gy (W).