Radiotherapy dose calculations in the presence of hip prostheses

The high density and atomic number of hip prostheses for patients undergoing pelvic radiotherapy challenge our ability to accurately calculate dose. A new clinical dose calculation algorithm, Monte Carlo, will allow accurate calculation of the radiation transport both within and beyond hip prostheses. The aim of this research was to investigate, for both phantom and patient geometries, the capability of various dose calculation algorithms to yield accurate treatment plans. Dose distributions in phantom and patient geometries with high atomic number prostheses were calculated using Monte Carlo, superposition, pencil beam, and no-heterogeneity correction algorithms. The phantom dose distributions were analyzed by depth dose and dose profile curves. The patient dose distributions were analyzed by isodose curves, dose-volume histograms (DVHs) and tumor control probability/normal tissue complication probability (TCP/NTCP) calculations. Monte Carlo calculations predicted the dose enhancement and reduction at the proximal and distal prosthesis interfaces respectively, whereas superposition and pencil beam calculations did not. However, further from the prosthesis, the differences between the dose calculation algorithms diminished. Treatment plans calculated with superposition showed similar isodose curves, DVHs, and TCP/NTCP as the Monte Carlo plans, except in the bladder, where Monte Carlo predicted a slightly lower dose. Treatment plans calculated with either the pencil beam method or with no heterogeneity correction differed significantly from the Monte Carlo plans.