Polynomial modelling of proton trajectories in homogeneous media for fast most likely path estimation and trajectory simulation.

Protons undergo many small angle deflections when traversing a medium, such as a patient. This effect, known as multiple Coulomb scattering (MCS), leads to degraded image resolution in proton radiography and computed tomography (CT) and to lateral spreading of the dose distribution in proton therapy. To optimally account for MCS in proton imaging, the most likely path (MLP) of a proton is estimated based on its position and propagation angle measured in front of and behind the object. In this work, we propose a functional which quantifies the likelihood of a proton trajectory and study how it can be used to model proton trajectories in a homogeneous medium. We focus on two aspects: first, we present an analytical method to quickly generate proton trajectories in a homogeneous medium based on the likelihood functional and validate it through Monte Carlo simulations. It could be used for fast generation of proton CT images without a full Monte Carlo simulation, or potentially to complement the components in a treatment planning Monte Carlo which simulate MCS. Second, by maximising the likelihood functional, we derive an expression for the MLP which is equivalent to the conventional ones reported in the literature yet computationally more convenient. Moreover, we show that the MLP is strictly a polynomial function if the protons' energy loss in the medium is approximated as a polynomial and that the orders of both are linked. We validate our MLP through Monte Carlo simulations and compare proton CT images reconstructed with our expression and with the conventional one. We find that an MLP polynomial of orders larger than five do not lead to increased spatial resolution compared to lower order expressions.