分类: 核科学技术 >> 辐射物理与技术 提交时间: 2023-06-01
摘要: In recent years, graphics processing units (GPUs) have been applied to accelerate Monte Carlo (MC) simulations for proton dose calculation in radiotherapy. Nonetheless, current GPU platforms, such as CUDA and OpenCL, suffer from cross-platform limitation or relatively high programming barrier. However, the Taichi toolkit, which was developed to overcome these difficulties, has been successfully applied to high-performance numerical computations. Based on the class II condensed history simulation scheme with various proton-nucleus interactions, we developed a proton MC transport GPU-accelerated engine using the Taichi toolkit. Dose distributions in homogeneous and heterogeneous materials were calculated for 110, 160, and 200 MeV protons and were compared with those obtained by full MC simulations using Topas. The gamma passing rates were greater than 0.99 and 0.95 with criteria of 2 mm, 2% and 1 mm, 1%, respectively, in all the tested conditions. Moreover, the calculation speed was at least 5800 times faster than that of Topas, and the number of lines of code was approximately 10 times lesser than those of CUDA or OpenCL. Our study provides a highly accurate, efficient, and easy-to-use proton dose calculation engine for algorithm developers, students, and medical physicists.
分类: 物理学 >> 核物理学 提交时间: 2025-01-06
摘要: Dose calculation is the foundation of Boron Neutron Capture Therapy (BNCT). MagicDose, a dose calcula- tion program for the BNCT treatment planning system, is developed based on the Monte Carlo method. Firstly, the voxel phantom of the modified Snyder head with 16, 8 mm are constructed, and the deviation of each result on y=x and the calculation time are statistical.The modified Snyder head phantom with tumor at three different spatial resolutions of 16, 8, and 1 mm are constructed, and the depth-dose-rate curves and spatial distribution maps are analyzed. Finally, the patient's head CT data is used as an application. The results show that the results calculated by MagicDose and MCNP are in good consistency, demonstrating that the computational efficiency of MagicDose is also better than that of MCNP. As the spatial resolution increases, the variability of the dose rate results is smaller. The voxel size and the number of threads are both inversely proportional to the time. For the CT model, the voxel phantom is successfully constructed and the calculation results are reasonable. The above results verify the correctness of MagicDose, which also provides a reference for optimizing the design of the voxel phantom in clinical treatment.