A method for the optimization of radiation therapy treatment plans is disclosed. The disclosed method is equally-applicable to robotic radiosurgery as well as other types of radiosurgical delivery, intensity-modulated radiotherapy (IMRT), volumetric modulated arc therapy (VMAT), and three-dimensional conformal radiotherapy (3DCRT). A population-based heuristic approximation is used to perform a global search, and subsequently, a deterministic local trajectory search is employed to further refine the initial solution.