Purpose: The purpose of this study is to improve upon a recently introduced bi-objective treatment planning method for prostate high-dose-rate (HDR) brachytherapy (BT), both in terms of resulting plan quality and runtime requirements, to the extent that its execution time is clinically acceptable. Methods: Bi-objective treatment planning is done using a state-of-the-art multiobjective evolutionary algorithm, which produces a large number of potential treatment plans with different trade-offs between coverage of the target volumes and sparing organs at risk. A graphics processing unit (GPU) is used for large-scale parallelization of dose calculations and the calculation of the dose-volume (DV) indices of potential treatment plans. Moreover, the objectives of the previously used bi-objective optimization model are modified to produce better results. Results: We applied the GPU-accelerated bi-objective treatment planning method to a set of 18 patients, resulting in a set containing a few hundred potential treatment plans with different trade-offs for each of these patients. Due to accelerations introduced in this article, results previously achieved after 1 hour are now achieved within 30 seconds of optimization. We found plans satisfying the clinical protocol for 15 of 18 patients, whereas this was the case for only 4 of 18 clinical plans. Higher quality treatment plans are obtained when the accuracy of DV index calculation is increased using more dose calculation points, requiring still no more than 3 minutes of optimization for 100 000 points. Conclusions: Large sets of high-quality treatment plans that trade-off coverage and sparing are now achievable within 30 seconds, due to the GPU-acceleration of a previously introduced bi-objective treatment planning method for prostate HDR brachytherapy. Higher quality plans can be achieved when optimizing for 3 minutes, which we still consider to be clinically acceptable. This allows for more insightful treatment plan selection in a clinical setting.

Purpose: To investigate the variation in computed dose-volume (DV) indices for high-dose-rate (HDR) prostate brachytherapy that can result from typical differences in computation settings in treatment planning systems (TPSs). Methods: Five factors were taken into account: number of dose-calculation points, radioactive source description, interpolation between delineated contours, intersections between delineated organ contours, and organ shape at the top and bottom contour using either full or partial slice thickness. Using in-house developed software, the DV indices of the treatment plans of 26 patients were calculated with different settings, and compared to a baseline setting that closely followed the default settings of the TPS used in our medical center. Studied organs were prostate and seminal vesicles, denoted as targets, and bladder, rectum, and urethra, denoted as organs at risk (OARs), which were delineated on MRI scans with a 3.3 mm slice thickness. Results: When sampling a fixed number of points in each organ, in order to achieve a width of the 95% confidence interval over all patients of the DV indices of 1% or less, only 32,000 points had to be sampled per target, but 256,000 points had to be sampled per OAR. For the remaining factors, DV indices changed up to 0.4% for rectum, 1.3% for urethra, and 2.6% for prostate. DV indices of the bladder changed especially if the high-dose-region was (partly) located at the most caudal contour, up to 8.5%, and DV indices of the vesicles changed especially if there were few delineated contours, up to 9.8%, both due to the use of full slice thickness for the top and bottom contour. Conclusions: The values of DV indices used in prostate HDR brachytherapy treatment planning are influenced by the computation settings in a TPS, especially at the most caudal part of the bladder, as well as in the seminal vesicles.