IMRT beam angle optimization using electromagnetism-like algorithm

The selection of appropriate beam irradiation directions in radiotherapy – beam angle optimization (BAO) problem – is very impor- tant for the quality of the treatment, both for improving tumor irradia- tion and for better organs sparing. However, the BAO problem is still not solved satisfactorily and, most of the time, beam directions continue to be manually selected in clinical practice which requires many trial and error iterations between selecting beam angles and computing fluence patterns until a suitable treatment is achieved. The objective of this pa- per is to introduce a new approach for the resolution of the BAO problem, using an hybrid electromagnetism-like algorithm with descent search to tackle this highly non-convex optimization problem. Electromagnetism- like algorithms are derivative-free optimization methods with the ability to avoid local entrapment. Moreover, the hybrid electromagnetism-like algorithm with descent search has a high ability of producing descent directions. A set of retrospective treated cases of head-and-neck tumors at the Portuguese Institute of Oncology of Coimbra is used to discuss the benefits of the proposed algorithm for the optimization of the BAO problem.Fundação para a Ciência e a Tecnologia (FCT

Theoretical methods for the calculation of Bragg curves and 3D distributions of proton beams

The well-known Bragg-Kleeman rule RCSDA = A dot E0p has become a pioneer work in radiation physics of charged particles and is still a useful tool to estimate the range RCSDA of approximately monoenergetic protons with initial energy E0 in a homogeneous medium. The rule is based on the continuous-slowing-down-approximation (CSDA). It results from a generalized (nonrelativistic) Langevin equation and a modification of the phenomenological friction term. The complete integration of this equation provides information about the residual energy E(z) and dE(z)/dz at each position z (0 <= z <= RCSDA). A relativistic extension of the generalized Langevin equation yields the formula RCSDA = A dot (E0 +E02/2M dot c2)p. The initial energy of therapeutic protons satisfies E0 << 2M dot c2 (M dot c2 = 938.276 MeV), which enables us to consider the relativistic contributions as correction terms. Besides this phenomenological starting-point, a complete integration of the Bethe-Bloch equation (BBE) is developed, which also provides the determination of RCSDA, E(z) and dE(z)/dz and uses only those parameters given by the BBE itself (i.e., without further empirical parameters like modification of friction). The results obtained in the context of the aforementioned methods are compared with Monte-Carlo calculations (GEANT4); this Monte-Carlo code is also used with regard to further topics such as lateral scatter, nuclear interactions, and buildup effects. In the framework of the CSDA, the energy transfer from protons to environmental atomic electrons does not account for local fluctuations.Comment: 97 pages review pape