Treatment planning optimization using constrained simulated annealing.

A variation of simulated annealing optimization called 'constrained simulated annealing' is used with a simple annealing schedule to optimize beam weights and angles in radiation therapy treatment planning. Constrained simulated annealing is demonstrated using two contrasting objective functions which incorporate both biological response and dose-volume considerations. The first objective function maximizes the probability of a complication-free treatment (PCFT) by minimizing the normal tissue complications subject to the constraint that the entire target volume receives a prescribed minimum turmourcidal dose with a specified dose homogeneity. Probabilities of normal tissue complication are based on published normal tissue complication probability functions and computed from dose-volume histograms. The second objective function maximizes the isocentre dose subject to a set of customized normal tissue dose-volume and target volume dose homogeneity constraints (MVDL). Although the PCFT objective function gives consistently lower estimates of normal tissue complication probabilities, the ability to specify individualized dose-volume limits, and therefore the individualized probability of complication, for an individual organ makes the MDVL objective function more useful for treatment planning.