Optimal radiotherapy treatment planning using minimum entropy models

We study the problem of finding an optimal radiotherapy treatment plan. A time-dependent Boltzmann particle transport model is used to model the interaction between radiative particles with tissue. This model allows for the modeling of inhomogeneities in the body’s tissues. It also allows for anisotropic distributed sources of radiation—as in brachytherapy—and external beam sources—as in teletherapy. We study two optimization problems: minimizing the deviation from a spatially-dependent prescribed dose through a quadratic tracking functional; and minimizing the survival of tumor cells through the use of the linear-quadratic model of radiobiological cell response. For each problem, we derive the optimality systems. In order to solve the state and adjoint equations, we use the minimum entropy approximation; the advantages of this method are discussed. Numerical results are then presented.