Archer(Rt) - A Gpu-Based And Photon-Electron Coupled Monte Carlo Dose Computing Engine For Radiation Therapy: Software Development And Application To Helical Tomotherapy
Lin Su,Youming Yang,Bryan Bednarz,Edmond Sterpin,Xining Du,Tianyu Liu,Wei Ji,X George Xu
DOI: https://doi.org/10.1118/1.4884229
IF: 4.506
2014-01-01
Medical Physics
Abstract:Purpose: Using the graphical processing units (GPU) hardware technology, an extremely fast Monte Carlo (MC) code ARCHER(RT) is developed for radiation dose calculations in radiation therapy. This paper describes the detailed software development and testing for three clinical TomoTherapy (R) cases: the prostate, lung, and head & neck.Methods: To obtain clinically relevant dose distributions, phase space files (PSFs) created from optimized radiation therapy treatment plan fluence maps were used as the input to ARCHERRT. Patient-specific phantoms were constructed from patient CT images. Batch simulations were employed to facilitate the time-consuming task of loading large PSFs, and to improve the estimation of statistical uncertainty. Furthermore, two different Woodcock tracking algorithms were implemented and their relative performance was compared. The dose curves of an Elekta accelerator PSF incident on a homogeneous water phantom were benchmarked against DOSXYZnrc. For each of the treatment cases, dose volume histograms and isodose maps were produced from ARCHER(RT) and the general-purpose code, GEANT4. The gamma index analysis was performed to evaluate the similarity of voxel doses obtained from these two codes. The hardware accelerators used in this study are one NVIDIA K20 GPU, one NVIDIA K40 GPU, and six NVIDIA M2090 GPUs. In addition, to make a fairer comparison of the CPU and GPU performance, a multithreaded CPU code was developed using OpenMP and tested on an Intel E5-2620 CPU.Results: For the water phantom, the depth dose curve and dose profiles from ARCHER(RT) agree well with DOSXYZnrc. For clinical cases, results from ARCHER(RT) are compared with those from GEANT4 and good agreement is observed. Gamma index test is performed for voxels whose dose is greater than 10% of maximum dose. For 2%/2mm criteria, the passing rates for the prostate, lung case, and head & neck cases are 99.7%, 98.5%, and 97.2%, respectively. Due to specific architecture of GPU, modified Woodcock tracking algorithm performed inferior to the original one. ARCHER(RT) achieves a fast speed for PSF-based dose calculations. With a single M2090 card, the simulations cost about 60, 50, 80 s for three cases, respectively, with the 1% statistical error in the PTV. Using the latest K40 card, the simulations are 1.7-1.8 times faster. More impressively, six M2090 cards could finish the simulations in 8.9-13.4 s. For comparison, the same simulations on Intel E5-2620 (12 hyperthreading) cost about 500-800 s.Conclusions: ARCHER(RT) was developed successfully to perform fast and accurate MC dose calculation for radiotherapy using PSFs and patient CT phantoms. c 2014 American Association of Physicists in Medicine.