Author(s)

Li Liao¹, Gino J. Lim¹, Xiaodong Zhang^2*

¹Department of Industrial Engineering, University of Houston, Houston, TX, USA.
²Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, TX, USA.

This paper uses a molecular dynamics (MD) method for intensity-modulated proton therapy treatment planning optimization to overcome the problem of gradient-based optimization methods such as quasi-Newton, which is sensitive to starting conditions and is easily trapped in local minima. We implemented a molecular dynamics (MD) method and a quasi-Newton method for plan optimization. Three types of cancer cases, prostate, head-and-neck and lung cancer, were tested with three starting different initial conditions. Overall, the MD method consistently resulted in solutions with lower objective function values (OFVs) compared to those from the quasi-Newton algorithm. Furthermore, the MD method converged on the same OFV regardless of its initial starting points used for the prostate cancer case.

Molecular Dynamics, Intensity-Modulated Proton Therapy, Quasi-Newton

Liao, L. , Lim, G. and Zhang, X. (2019) A Molecular Dynamics Approach for Optimizing Beam Intensities in IMPT Treatment Planning. Journal of Applied Mathematics and Physics, 7, 2130-2047. doi: 10.4236/jamp.2019.79146.