A computational study on different penalty approaches for constrained optimization in radiation therapy treatment planning with a simulated annealing algorithm

In intensity modulated radiation therapy (IMRT) a treatment plan is a high dimensionality optimization problem with the goal to give the prescribed radiation dose to the Planning Target Volume (PTV) while sparing critical organs. A clinically acceptable plan is usually generated by a numerical optimization process in pursuit of attaining the above mentioned goal. Incorporation of dose volume constraints (DVCs) for the OARs introduce an additional degree of impediment to the optimization task. Heuristic algorithms have been ascertained in the past as a powerful tool for various problems in radiation therapy, such as beam angle optimization (BAO) and IMRT treatment planning. Simulated Annealing (SA) algorithm has the capability to find global minima for bound-constrained optimization problems. However, its performance depends on the penalty method which is utilized for the constraints. In the current study we investigate the performance of three penalty methods for IMRT treatment planning for five prostate cases with dose volume constraints. In addition, sensitivity analysis was performed in order to study the effect of the penalty variables on the treatment plan. Finally the merits and demerits of each method have been discussed.

[1]  Sandro Ridella,et al.  Minimizing multimodal functions of continuous variables with the “simulated annealing” algorithmCorrigenda for this article is available here , 1987, TOMS.

[2]  Yixin Chen,et al.  Simulated annealing with asymptotic convergence for nonlinear constrained optimization , 2007, J. Glob. Optim..

[3]  Xin-She Yang,et al.  Nature-Inspired Metaheuristic Algorithms: Second Edition , 2010 .

[4]  D. Mitra,et al.  Convergence and finite-time behavior of simulated annealing , 1985, 1985 24th IEEE Conference on Decision and Control.

[5]  Sandro Ridella,et al.  Corrigenda: "Minimizing Multimodal Functions of Continuous Variables with the 'Simulated Annealing' Algorithm" , 1989, ACM Trans. Math. Softw..

[6]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[7]  E. Bernhard,et al.  How does radiation kill cells? , 1999, Current opinion in chemical biology.

[8]  Xin-She Yang,et al.  Review of Metaheuristics and Generalized Evolutionary Walk Algorithm , 2011, 1105.3668.

[9]  Christian Blum,et al.  Metaheuristics in combinatorial optimization: Overview and conceptual comparison , 2003, CSUR.

[10]  Joseph O Deasy,et al.  CERR: a computational environment for radiotherapy research. , 2003, Medical physics.

[11]  T Bortfeld,et al.  Optimized planning using physical objectives and constraints. , 1999, Seminars in radiation oncology.

[12]  Yu Lei,et al.  A self-tuned bat algorithm for optimization in radiation therapy treatment planning , 2014, 15th IEEE/ACIS International Conference on Software Engineering, Artificial Intelligence, Networking and Parallel/Distributed Computing (SNPD).

[13]  Michel Gendreau,et al.  Metaheuristics in Combinatorial Optimization , 2022 .

[14]  T. Bortfeld.  IMRT: a review and preview , 2006, Physics in medicine and biology.