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

We have determined the depth-dose curve, penetration range, and water equivalent ratio 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, which 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 in its widely used Report No. 49 [M. J. Berger et al., Stopping powers and ranges for protons and alpha particles, International Commission on Radiation Units and Measurements, Bethesda, Maryland, 1993]. The results of this paper emphasize the importance of an accurate determination of the stopping quantities of cortical bone to advance towards the millimetric precision for the proton penetration ranges and deposited dose needed in radiotherapy.