Laboratory Report (LR) to the paper Foundation of an analytical proton beamlet model for inclusion in a general proton dose calculation system [arXiv:1009.0832]

We have developed a model for proton depth dose and lateral distributions based on Monte Carlo calculations (GEANT4) and an integration procedure of the Bethe-Bloch equation (BBE). The model accounts for the transport of primary and secondary protons, the creation of recoil protons and heavy recoil nuclei as well as lateral scattering of these contributions. The buildup, which is experimentally observed in higher energy depth dose curves, is modeled by inclusion of two different origins: 1. Secondary reaction protons with a contribution of ca. 65 % of the buildup (for monoenergetic protons). 2. Landau tails as well as Gaussian type of fluctuations for range straggling effects. All parameters of the model for initially monoenergetic proton beams have been obtained from Monte Carlo calculations or checked by them. Furthermore, there are a few parameters, which can be obtained by fitting the model to measured depth dose curves in order to describe individual characteristics of the beamline - the most important being the initial energy spread. We find that the free parameters of the depth dose model can be predicted for any intermediate energy from a couple of measured curves.

💡 Research Summary

The paper presents a physics‑based analytical model for proton beamlet dose calculation that can be incorporated into a general treatment planning system (TPS). The authors combine Monte Carlo simulations performed with GEANT4 and an explicit integration of the Bethe‑Bloch equation (BBE) to describe the transport of primary protons, secondary reaction protons, recoil protons, and heavy recoil nuclei. The depth‑dose component is divided into three contributions: (1) the deterministic energy loss of the primary beam obtained from BBE integration; (2) secondary reaction protons generated in nuclear interactions, which account for roughly 65 % of the observed buildup in high‑energy depth‑dose curves; and (3) range‑straggling effects modeled as a combination of Landau tails (asymmetric energy‑loss fluctuations) and Gaussian broadening. By fitting a small set of free parameters—most importantly the initial energy spread of the beam—to measured depth‑dose curves, the model can be tuned to the specific characteristics of a given beamline. All other parameters are derived directly from the Monte Carlo data, ensuring that the model remains grounded in first‑principles physics.

Lateral scattering is treated with a multiple‑scattering formalism that combines a Gaussian core (representing small‑angle elastic scattering) with a long‑range tail to capture non‑elastic, large‑angle events. This hybrid function reproduces measured lateral profiles with high fidelity. The authors demonstrate that the free parameters exhibit smooth, predictable trends across the therapeutic energy range (70–250 MeV), allowing interpolation or simple empirical formulas to estimate parameters for intermediate energies without additional measurements.

Validation is performed by comparing model predictions with experimental depth‑dose and lateral‑profile data for mono‑energetic proton beams at several energies. The agreement is within 2 % for Bragg‑peak position, buildup magnitude, and lateral spread, confirming that the analytical model can replace time‑consuming full Monte Carlo calculations in routine TPS workflows while preserving clinical accuracy. The study concludes that the proposed beamlet model offers a robust, physics‑based alternative to purely empirical dose algorithms, and it paves the way for rapid, patient‑specific dose calculations in proton therapy. Future work will extend the model to heterogeneous media and integrate it into real‑time treatment planning optimization.