Keywords:
Monte Carlo simulation, metropolis hasting sampling, path-by-path scheme, scattering calculation
Low computational efficiency is the largest obstacle hindering practical applications of Monte Carlo particle transport simulation. We develop a new MC scheme for photon transport that employs Metropolis-Hasting sampling algorithm to conduct the simulations via a path-by-path scheme in this GPU-based Metropolis MC (gMMC) package. By using the Metropolis algorithm, gMMC is able to sample an entire path of a photon each time, extending from the x-ray source to the detector. The sampled paths over a long run will follow a distribution governed by the particle transport physics. gMMC was benchmarked against an in-house developed GPU-based Monte Carol simulation tool gMCDRR that performs simulations in the conventional particle-by-particle scheme. An example problem of x-ray scatter calculation was studied. For pencil beam case, speed-up factors of 200 to 300 times for first order scatter and 12 to 18 times for second order scatter can be achieved, while the average differences comparing to the conventional approach are within 1%.
- Yusi Chen
- Southern Medical University
- Jianhui Ma
- Southern Medical University
- Bin Li
- Southern Medical University
- Zhen Tian
- UT Southern Medical Center
- Linghong Zhou
- Southern Medical University
- Xun Jia
- UT Southern Medical Center
- Yuan Xu
- Southern Medical University
- Y. Xu, T. Bai, H. Yan, L. Ouyang, A. Pompos, J. Wang, L. Zhou, S. B. Jiang, and X. Jia, "A practical cone-beam CT scatter correction method with optimized Monte Carlo simulations for image-guided radiation therapy," Physics in Medicine & Biology, vol. 60, pp. 3567-87, 2015.
- G. Pratx and L. Xing, "GPU computing in medical physics: a review," Medical Physics, vol. 38, pp. 2685-97, 2011.
- X. Jia, P. Ziegenhein, and S. B. Jiang, "GPU-based high-performance computing for radiation therapy,“ Physics in Medicine and Biology, vol. 59, p. R151, 2014.
- X. Jia, H. Yan, L. Cervino, M. Folkerts, and S. B. Jiang, "A GPU tool for efficient, accurate, and realistic simulation of cone beam CT projections," Med Phys, vol. 39, pp. 7368-78, Dec 2012.
- X. Jia, H. Yan, X. Gu, and S. B. Jiang, "Fast Monte Carlo simulation for patient-specific CT/CBCT imaging dose calculation," Phys Med Biol, vol. 57, pp. 577-90, Feb 07 2012.
- N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, "Equation of State Calculations by Fast Computing Machines," Journal of Biochemical & Biophysical Methods, vol. 21, pp. 1087-1092, 1952.
- W. K. Hastings, "Monte Carlo sampling methods using Markov chains and their applications," Biometrika, vol. 57, pp. 97-109, 1970.
- C. P. Robert and G. Casella, The Metropolis—Hastings Algorithm: Springer New York, 2004.