Search

Purpose: This paper provides a comparison between a fast, commercial, in-patient Monte Carlo dose calculation algorithm (GPUMCD) and geant4. It also evaluates the dosimetric impact of the application of an external 1.5 T magnetic field., Methods: A stand-alone version of the Elekta™ GPUMCD algorithm, to be used within the Monaco treatment planning system to model dose for the Elekta™ magnetic resonance imaging (MRI) Linac, was compared against GEANT4 (v10.1). This was done in the presence or absence of a 1.5 T static magnetic field directed orthogonally to the radiation beam axis. Phantoms with material compositions of water, ICRU lung, ICRU compact-bone, and titanium were used for this purpose. Beams with 2 MeV monoenergetic photons as well as a 7 MV histogrammed spectrum representing the MRI Linac spectrum were emitted from a point source using a nominal source-to-surface distance of 142.5 cm. Field sizes ranged from 1.5 × 1.5 to 10 × 10 cm(2). Dose scoring was performed using a 3D grid comprising 1 mm(3) voxels. The production thresholds were equivalent for both codes. Results were analyzed based upon a voxel by voxel dose difference between the two codes and also using a volumetric gamma analysis., Results: Comparisons were drawn from central axis depth doses, cross beam profiles, and isodose contours. Both in the presence and absence of a 1.5 T static magnetic field the relative differences in doses scored along the beam central axis were less than 1% for the homogeneous water phantom and all results matched within a maximum of ±2% for heterogeneous phantoms. Volumetric gamma analysis indicated that more than 99% of the examined volume passed gamma criteria of 2%-2 mm (dose difference and distance to agreement, respectively). These criteria were chosen because the minimum primary statistical uncertainty in dose scoring voxels was 0.5%. The presence of the magnetic field affects the dose at the interface depending upon the density of the material on either sides of the interface. This effect varies with the field size. For example, at the water-lung interface a 33.94% increase in dose was observed (relative to the Dmax), by both GPUMCD and GEANT4 for the field size of 2 × 2 cm(2) (compared to no B-field case), which increased to 47.83% for the field size of 5 × 5 cm(2) in the presence of the magnetic field. Similarly, at the lung-water interface, the dose decreased by 19.21% (relative to Dmax) for a field size of 2 × 2 cm(2) and by 30.01% for 5 × 5 cm(2) field size. For more complex combinations of materials the dose deposition also becomes more complex., Conclusions: The GPUMCD algorithm showed good agreement against GEANT4 both in the presence and absence of a 1.5 T external magnetic field. The application of 1.5 T magnetic field significantly alters the dose at the interfaces by either increasing or decreasing the dose depending upon the density of the material on either side of the interfaces.

The aim of this study was to evaluate the potential of bGPUMCD, a Monte Carlo algorithm executed on Graphics Processing Units (GPUs), for fast dose calculations in permanent prostate implant dosimetry. It also aimed to validate a low dose rate brachytherapy source in terms of TG-43 metrics and to use this source to compute dose distributions for permanent prostate implant in very short times. The physics of bGPUMCD was reviewed and extended to include Rayleigh scattering and fluorescence from photoelectric interactions for all materials involved. The radial and anisotropy functions were obtained for the Nucletron SelectSeed in TG-43 conditions. These functions were compared to those found in the MD Anderson Imaging and Radiation Oncology Core brachytherapy source registry which are considered the TG-43 reference values. After appropriate calibration of the source, permanent prostate implant dose distributions were calculated for four patients and compared to an already validated Geant4 algorithm. The radial function calculated from bGPUMCD showed excellent agreement (differences within 1.3%) with TG-43 accepted values. The anisotropy functions at r = 1 cm and r = 4 cm were within 2% of TG-43 values for angles over 17.5°. For permanent prostate implants, Monte Carlo-based dose distributions with a statistical uncertainty of 1% or less for the target volume were obtained in 30 s or less for 1 × 1 × 1 mm(3) calculation grids. Dosimetric indices were very similar (within 2.7%) to those obtained with a validated, independent Monte Carlo code (Geant4) performing the calculations for the same cases in a much longer time (tens of minutes to more than a hour). bGPUMCD is a promising code that lets envision the use of Monte Carlo techniques in a clinical environment, with sub-minute execution times on a standard workstation. Future work will explore the use of this code with an inverse planning method to provide a complete Monte Carlo-based planning solution.

During the last decade, studies have shown that 3D list-mode ordered-subset expectation-maximization (LM-OSEM) algorithms for positron emission tomography (PET) reconstruction could be effectively computed and considerably accelerated by graphics processing unit (GPU) devices. However, most of these studies rely on pre-calculated sensitivity matrices. In many cases, the time required to compute this matrix can be longer than the reconstruction time itself. In fact, the relatively long time required for the calculation of the patient-specific sensitivity matrix is considered as one of the main obstacle in introducing a list-mode PET reconstruction algorithm for routine clinical use. The objective of this work is to accelerate a fully 3D LM-OSEM algorithm, including the calculation of the sensitivity matrix that accounts for the patient-specific attenuation and normalization corrections. For this purpose, sensitivity matrix calculations and list-mode OSEM reconstructions were implemented on GPUs, using the geometry of a commercial PET system. The system matrices were built on-the-fly by using an approach with multiple rays per detector pair. The reconstructions were performed for a volume of 188 × 188 × 57 voxels of 2 × 2 × 3.15 mm(3) and for another volume of 144 × 144 × 57 voxels of 4 × 4 × 3.15 mm(3). The time to compute the sensitivity matrix for the 188 × 188 × 57 array was 9 s while the LM-OSEM algorithm performed at a rate of 1.1 millions of events per second. For the 144 × 144 × 57 array, the respective numbers are 8 s for the sensitivity matrix and 0.8 million of events per second for the LM-OSEM step. This work lets envision fast reconstructions for advanced PET applications such as real time dynamic studies and parametric image reconstructions.

Purpose: To establish the accuracy and speed of bGPUMCD, a GPU-oriented Monte Carlo code used for high dose rate brachytherapy dose calculations. The first objective is to evaluate the time required for dose calculation when full Monte Carlo generated dose distribution kernels are used for plan optimization. The second objective is to assess the accuracy and speed when recalculating pre-optimized plans, consisting of many dwell positions., Methods: bGPUMCD is tested with three clinical treatment plans : one prostate case, one breast case, and one rectum case with a shielded applicator. Reference distributions, generated with GEANT4, are used as a basis of comparison. Calculations of full dose distributions of pre-optimized treatment plans as well as single dwell dosimetry are performed. Single source dosimetry, based on TG-43 parameters reproduction, is also presented for the microSelectron V2 (Nucletron, Veenendaal, The Netherlands)., Results: In timing experiments, the computation of single dwell position dose kernels takes between 0.25 and 0.5 s. bGPUMCD can compute full dose distributions of previously optimized plans in ∼2 s. bGPUMCD is capable of computing pre-optimized brachytherapy plans within 1% for the prostate case and 2% for the breast and shielded applicator cases, when comparing the dosimetric parameters D90 and V100 of the reference (GEANT4) and bGPUMCD distributions. For all voxels within the target, an absolute average difference of approximately 1% is found for the prostate case, less than 2% for the breast case and less than 2% for the rectum case with shielded applicator. Larger point differences (>5%) are found within bony regions in the prostate case, where bGPUMCD underdoses compared to GEANT4. Single source dosimetry results are mostly within 2% for the radial function and within 1%-4% for the anisotropic function., Conclusions: bGPUMCD has the potential to allow for fast MC dose calculation in a clinical setting for all phases of HDR treatment planning, from dose kernel calculations for plan optimization to plan recalculation.

Purpose: To validate GPUMCD, a new package for fast Monte Carlo dose calculations based on the GPU (graphics processing unit), as a tool for low-energy single seed brachytherapy dosimetry for specific seed models. As the currently accepted method of dose calculation in low-energy brachytherapy computations relies on severe approximations, a Monte Carlo based approach would result in more accurate dose calculations, taking in to consideration the patient anatomy as well as interseed attenuation. The first step is to evaluate the capability of GPUMCD to reproduce low-energy, single source, brachytherapy calculations which could ultimately result in fast and accurate, Monte Carlo based, brachytherapy dose calculations for routine planning., Methods: A mixed geometry engine was integrated to GPUMCD capable of handling parametric as well as voxelized geometries. In order to evaluate GPUMCD for brachytherapy calculations, several dosimetry parameters were computed and compared to values found in the literature. These parameters, defined by the AAPM Task-Group No. 43, are the radial dose function, the 2D anisotropy function, and the dose rate constant. These three parameters were computed for two different brachytherapy sources: the Amersham OncoSeed 6711 and the Imagyn IsoStar IS-12501., Results: GPUMCD was shown to yield dosimetric parameters similar to those found in the literature. It reproduces radial dose functions to within 1.25% for both sources in the 0.5< r <10 cm range. The 2D anisotropy function was found to be within 3% at r =5 cm and within 4% at r = 1 cm. The dose rate constants obtained were within the range of other values reported in the literature., Conclusion: GPUMCD was shown to be able to reproduce various TG-43 parameters for two different low-energy brachytherapy sources found in the literature. The next step is to test GPUMCD as a fast clinical Monte Carlo brachytherapy dose calculations with multiple seeds and patient geometry, potentially providing more accurate results than the TG-43 formalism while being much faster than calculations using general purpose Monte Carlo codes.

Purpose: Monte Carlo methods are considered as the gold standard for dosimetric computations in radiotherapy. Their execution time is, however, still an obstacle to the routine use of Monte Carlo packages in a clinical setting. To address this problem, a completely new, and designed from the ground up for the GPU, Monte Carlo dose calculation package for voxelized geometries is proposed: GPUMCD., Method: GPUMCD implements a coupled photon-electron Monte Carlo simulation for energies in the range of 0.01-20 MeV. An analog simulation of photon interactions is used and a class II condensed history method has been implemented for the simulation of electrons. A new GPU random number generator, some divergence reduction methods, as well as other optimization strategies are also described. GPUMCD was run on a NVIDIA GTX480, while single threaded implementations of EGSnrc and DPM were run on an Intel Core i7 860., Results: Dosimetric results obtained with GPUMCD were compared to EGSnrc. In all but one test case, 98% or more of all significant voxels passed the gamma criteria of 2%-2 mm. In terms of execution speed and efficiency, GPUMCD is more than 900 times faster than EGSnrc and more than 200 times faster than DPM, a Monte Carlo package aiming fast executions. Absolute execution times of less than 0.3 s are found for the simulation of 1M electrons and 4M photons in water for monoenergetic beams of 15 MeV, including GPU-CPU memory transfers., Conclusion: GPUMCD, a new GPU-oriented Monte Carlo dose calculation platform, has been compared to EGSnrc and DPM in terms of dosimetric results and execution speed. Its accuracy and speed make it an interesting solution for full Monte Carlo dose calculation in radiation oncology.

Purpose: Graphic processing units (GPUs) are increasingly used for scientific applications, where their parallel architecture and unprecedented computing power density can be exploited to accelerate calculations. In this paper, a new GPU implementation of a convolution/superposition (CS) algorithm is presented., Methods: This new GPU implementation has been designed from the ground-up to use the graphics card's strengths and to avoid its weaknesses. The CS GPU algorithm takes into account beam hardening, off-axis softening, kernel tilting, and relies heavily on raytracing through patient imaging data. Implementation details are reported as well as a multi-GPU solution., Results: An overall single-GPU acceleration factor of 908x was achieved when compared to a nonoptimized version of the CS algorithm implemented in PlanUNC in single threaded central processing unit (CPU) mode, resulting in approximatively 2.8 s per beam for a 3D dose computation on a 0.4 cm grid. A comparison to an established commercial system leads to an acceleration factor of approximately 29x or 0.58 versus 16.6 s per beam in single threaded mode. An acceleration factor of 46x has been obtained for the total energy released per mass (TERMA) calculation and a 943x acceleration factor for the CS calculation compared to PlanUNC. Dose distributions also have been obtained for a simple water-lung phantom to verify that the implementation gives accurate results., Conclusions: These results suggest that GPUs are an attractive solution for radiation therapy applications and that careful design, taking the GPU architecture into account, is critical in obtaining significant acceleration factors. These results potentially can have a significant impact on complex dose delivery techniques requiring intensive dose calculations such as intensity-modulated radiation therapy (IMRT) and arc therapy. They also are relevant for adaptive radiation therapy where dose results must be obtained rapidly.

The numerical calculation of dose is central to treatment planning in radiation therapy and is at the core of optimization strategies for modern delivery techniques. In a clinical environment, dose calculation algorithms are required to be accurate and fast. The accuracy is typically achieved through the integration of patient-specific data and extensive beam modeling, which generally results in slower algorithms. In order to alleviate execution speed problems, the authors have implemented a modern dose calculation algorithm on a massively parallel hardware architecture. More specifically, they have implemented a convolution-superposition photon beam dose calculation algorithm on a commodity graphics processing unit (GPU). They have investigated a simple porting scenario as well as slightly more complex GPU optimization strategies. They have achieved speed improvement factors ranging from 10 to 20 times with GPU implementations compared to central processing unit (CPU) implementations, with higher values corresponding to larger kernel and calculation grid sizes. In all cases, they preserved the numerical accuracy of the GPU calculations with respect to the CPU calculations. These results show that streaming architectures such as GPUs can significantly accelerate dose calculation algorithms and let envision benefits for numerically intensive processes such as optimizing strategies, in particular, for complex delivery techniques such as IMRT and are therapy.