Experimental assessment of proton dose calculation accuracy in inhomogeneous media.

PURPOSE: Proton therapy with Pencil Beam Scanning (PBS) has the potential to improve radiotherapy treatments. Unfortunately, its promises are jeopardized by the sensitivity of the dose distributions to uncertainties, including dose calculation accuracy in inhomogeneous media. Monte Carlo dose engines (MC) are expected to handle heterogeneities better than analytical algorithms like the pencil-beam convolution algorithm (PBA). In this study, an experimental phantom has been devised to maximize the effect of heterogeneities and to quantify the capability of several dose engines (MC and PBA) to handle these. METHODS: An inhomogeneous phantom made of water surrounding a long insert of bone tissue substitute (1×10×10 cm3) was irradiated with a mono-energetic PBS field (10×10 cm2). A 2D ion chamber array (MatriXX, IBA Dosimetry GmbH) lied right behind the bone. The beam energy was such that the expected range of the protons exceeded the detector position in water and did not attain it in bone. The measurement was compared to the following engines: Geant4.9.5, PENH, MCsquare, as well as the MC and PBA algorithms of RayStation (RaySearch Laboratories AB). RESULTS: For a γ-index criteria of 2%/2mm, the passing rates are 93.8% for Geant4.9.5, 97.4% for PENH, 93.4% for MCsquare, 95.9% for RayStation MC, and 44.7% for PBA. The differences in γ-index passing rates between MC and RayStation PBA calculations can exceed 50%. CONCLUSION: The performance of dose calculation algorithms in highly inhomogeneous media was evaluated in a dedicated experiment. MC dose engines performed overall satisfactorily while large deviations were observed with PBA as expected.

Dose distribution for gynecological brachytherapy with dose accumulation between insertions: Feasibility study.

PURPOSE: For gynecological treatments, it is standard to acquire CT images and preferably also MR images before each treatment to calculate the dose of the day. The dose of the complete treatment is calculated by adding the dose metrics of each fraction. It makes the conservative assumption that the same part of the organs at risk always receives the highest dose. The dose calculated this way often limits the prescription dose or the target coverage. We investigated the use of deformable image registration (DIR) as an alternative method to assess the cumulative dose for a treatment course. METHODS AND MATERIALS: Rigid registration is preformed on CT images, followed by DIR. DIR can be based either solely on the three-dimensional images or combined with organ contours. To improve DIR in the pelvic region with low CT contrast, we propose (1) using contours drawn on CT or (2) modifying artificially the contrast in certain volumes. The dose matrix from fraction_n (n > 1) is deformed using a calculated deformation field. RESULTS: The use of the contrast-enhanced images or of contour information helps to guide the DIR. However, because of the very high dose gradients involved in brachytherapy, the uncertainty on the accumulated dose remains of the order of 5-10%. Even for good contour matching, a small local error in the deformation can have significant consequences for the dose distribution. CONCLUSIONS: Using DIR, based on image features and contours, allows to accumulate the dose from different brachytherapy fractions. A robust validation procedure should be developed.