Simulation of depth-dose curves and water equivalent ratios of energetic proton beams in cortical bone

TL;DR

This study uses advanced simulation techniques combining molecular dynamics and Monte Carlo methods to accurately determine proton beam depth-dose curves and water equivalent ratios in cortical bone, crucial for precise radiotherapy.

Contribution

It introduces a detailed simulation approach that integrates molecular dynamics with Monte Carlo methods, providing more accurate data for proton therapy in cortical bone.

Findings

Simulation results align well with scarce experimental data

Deviations from ICRU standard stopping powers highlight importance of accurate inputs

Emphasizes need for precise cortical bone stopping quantities in radiotherapy

Abstract

We have determined the depth-dose curve, the penetration range, and the water equivalent ratio (WER), for proton beams of clinical energies in cortical bone, by means of a detailed and accurate simulation that combines molecular dynamics and Monte Carlo techniques. The fundamental input quantities (stopping power and energy loss straggling) for the simulation were obtained from a reliable electronic excitation spectrum of the condensed-phase target, which takes into account the organic and mineral phases that form it. Our simulations with these inputs, that are in excellent agreement with the scarce data available for a cortical bone target, deviate from simulations performed using other stopping quantities, such as those provided by the International Commission on Radiation Units and Measurements (ICRU) in its widely used Report 49. The results of this work emphasize the importance of an accurate determination of the stopping quantities of cortical bone in order to advance towards the millimetric precision for the proton penetration ranges and deposited dose needed in radiotherapy.