{"title":"一种gpu加速的蒙特卡罗剂量引擎用于外射束放射治疗。","authors":"Zihao Liu, Yuxiang Wang, Yiqun Han, Panpan Hu, Cheng Zheng, Bing Yan, Yidong Yang","doi":"10.1002/mp.17899","DOIUrl":null,"url":null,"abstract":"<div>\n \n \n <section>\n \n <h3> Background</h3>\n \n <p>Accurate dose computation is crucial in intensity-modulated radiation therapy. Owing to its high accuracy, Monte Carlo method is considered the gold standard for radiation dose computation. Its efficiency, however, demands continuous improvement.</p>\n </section>\n \n <section>\n \n <h3> Purpose</h3>\n \n <p>This study aims to develop a GPU-accelerated Monte Carlo radiation dose engine (GARDEN) for fast and accurate dose computation in external beam radiotherapy.</p>\n </section>\n \n <section>\n \n <h3> Methods</h3>\n \n <p>In GARDEN simulation, photon and electron transport were modeled using Woodcock tracking and Class II condensed history technique, respectively. To enhance GPU computational efficiency, warp convergence optimization and coalesced access methods were employed. A novel linear accelerator (Linac) head model was established by incorporating a virtual source and a digital collimator system. The physics was verified against GEANT4 in both homogeneous and heterogeneous phantoms. The Linac head model was commissioned using data measured in a water tank and validated by comparing simulation with film doses for two alternating open and closed MLC patterns. Finally, computational efficiency and accuracy were further evaluated in clinical IMRT and VMAT treatment plans.</p>\n </section>\n \n <section>\n \n <h3> Results</h3>\n \n <p>GARDEN was more than 2500 times faster than GEANT4, with dose differences less than 1% in both homogeneous water and heterogeneous water-lung-bone phantoms. Compared to commission data, the average differences in percentage depth dose curves were less than 1%, and the penumbra differences in lateral dose profiles were less than 1 mm for various radiation field sizes. For two MLC patterns, the gamma pass rates between GARDEN simulations and films were 98.78% and 97.30% at 2%/2 mm criteria, respectively. Both IMRT and VMAT treatment plans achieved gamma pass rates exceeding 99.23% at 3%/3 mm criteria compared to GEANT4 results, with GARDEN completing the dose calculations within 3 s at ∼1% uncertainty on an i9-13900K CPU and NVIDIA 4080 GPU.</p>\n </section>\n \n <section>\n \n <h3> Conclusion</h3>\n \n <p>The accuracy and efficiency of GARDEN has been benchmarked against GEANT4 and validated in both phantoms and clinical treatment plans. With its capability for fast and accurate dose computation, GARDEN shows strong potential for applications in treatment planning and quality assurance.</p>\n </section>\n </div>","PeriodicalId":18384,"journal":{"name":"Medical physics","volume":"52 7","pages":""},"PeriodicalIF":3.2000,"publicationDate":"2025-05-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"A GPU-accelerated Monte Carlo dose engine for external beam radiotherapy\",\"authors\":\"Zihao Liu, Yuxiang Wang, Yiqun Han, Panpan Hu, Cheng Zheng, Bing Yan, Yidong Yang\",\"doi\":\"10.1002/mp.17899\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div>\\n \\n \\n <section>\\n \\n <h3> Background</h3>\\n \\n <p>Accurate dose computation is crucial in intensity-modulated radiation therapy. Owing to its high accuracy, Monte Carlo method is considered the gold standard for radiation dose computation. Its efficiency, however, demands continuous improvement.</p>\\n </section>\\n \\n <section>\\n \\n <h3> Purpose</h3>\\n \\n <p>This study aims to develop a GPU-accelerated Monte Carlo radiation dose engine (GARDEN) for fast and accurate dose computation in external beam radiotherapy.</p>\\n </section>\\n \\n <section>\\n \\n <h3> Methods</h3>\\n \\n <p>In GARDEN simulation, photon and electron transport were modeled using Woodcock tracking and Class II condensed history technique, respectively. To enhance GPU computational efficiency, warp convergence optimization and coalesced access methods were employed. A novel linear accelerator (Linac) head model was established by incorporating a virtual source and a digital collimator system. The physics was verified against GEANT4 in both homogeneous and heterogeneous phantoms. The Linac head model was commissioned using data measured in a water tank and validated by comparing simulation with film doses for two alternating open and closed MLC patterns. Finally, computational efficiency and accuracy were further evaluated in clinical IMRT and VMAT treatment plans.</p>\\n </section>\\n \\n <section>\\n \\n <h3> Results</h3>\\n \\n <p>GARDEN was more than 2500 times faster than GEANT4, with dose differences less than 1% in both homogeneous water and heterogeneous water-lung-bone phantoms. Compared to commission data, the average differences in percentage depth dose curves were less than 1%, and the penumbra differences in lateral dose profiles were less than 1 mm for various radiation field sizes. For two MLC patterns, the gamma pass rates between GARDEN simulations and films were 98.78% and 97.30% at 2%/2 mm criteria, respectively. Both IMRT and VMAT treatment plans achieved gamma pass rates exceeding 99.23% at 3%/3 mm criteria compared to GEANT4 results, with GARDEN completing the dose calculations within 3 s at ∼1% uncertainty on an i9-13900K CPU and NVIDIA 4080 GPU.</p>\\n </section>\\n \\n <section>\\n \\n <h3> Conclusion</h3>\\n \\n <p>The accuracy and efficiency of GARDEN has been benchmarked against GEANT4 and validated in both phantoms and clinical treatment plans. With its capability for fast and accurate dose computation, GARDEN shows strong potential for applications in treatment planning and quality assurance.</p>\\n </section>\\n </div>\",\"PeriodicalId\":18384,\"journal\":{\"name\":\"Medical physics\",\"volume\":\"52 7\",\"pages\":\"\"},\"PeriodicalIF\":3.2000,\"publicationDate\":\"2025-05-29\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Medical physics\",\"FirstCategoryId\":\"3\",\"ListUrlMain\":\"https://onlinelibrary.wiley.com/doi/10.1002/mp.17899\",\"RegionNum\":2,\"RegionCategory\":\"医学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"RADIOLOGY, NUCLEAR MEDICINE & MEDICAL IMAGING\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Medical physics","FirstCategoryId":"3","ListUrlMain":"https://onlinelibrary.wiley.com/doi/10.1002/mp.17899","RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"RADIOLOGY, NUCLEAR MEDICINE & MEDICAL IMAGING","Score":null,"Total":0}
A GPU-accelerated Monte Carlo dose engine for external beam radiotherapy
Background
Accurate dose computation is crucial in intensity-modulated radiation therapy. Owing to its high accuracy, Monte Carlo method is considered the gold standard for radiation dose computation. Its efficiency, however, demands continuous improvement.
Purpose
This study aims to develop a GPU-accelerated Monte Carlo radiation dose engine (GARDEN) for fast and accurate dose computation in external beam radiotherapy.
Methods
In GARDEN simulation, photon and electron transport were modeled using Woodcock tracking and Class II condensed history technique, respectively. To enhance GPU computational efficiency, warp convergence optimization and coalesced access methods were employed. A novel linear accelerator (Linac) head model was established by incorporating a virtual source and a digital collimator system. The physics was verified against GEANT4 in both homogeneous and heterogeneous phantoms. The Linac head model was commissioned using data measured in a water tank and validated by comparing simulation with film doses for two alternating open and closed MLC patterns. Finally, computational efficiency and accuracy were further evaluated in clinical IMRT and VMAT treatment plans.
Results
GARDEN was more than 2500 times faster than GEANT4, with dose differences less than 1% in both homogeneous water and heterogeneous water-lung-bone phantoms. Compared to commission data, the average differences in percentage depth dose curves were less than 1%, and the penumbra differences in lateral dose profiles were less than 1 mm for various radiation field sizes. For two MLC patterns, the gamma pass rates between GARDEN simulations and films were 98.78% and 97.30% at 2%/2 mm criteria, respectively. Both IMRT and VMAT treatment plans achieved gamma pass rates exceeding 99.23% at 3%/3 mm criteria compared to GEANT4 results, with GARDEN completing the dose calculations within 3 s at ∼1% uncertainty on an i9-13900K CPU and NVIDIA 4080 GPU.
Conclusion
The accuracy and efficiency of GARDEN has been benchmarked against GEANT4 and validated in both phantoms and clinical treatment plans. With its capability for fast and accurate dose computation, GARDEN shows strong potential for applications in treatment planning and quality assurance.
期刊介绍:
Medical Physics publishes original, high impact physics, imaging science, and engineering research that advances patient diagnosis and therapy through contributions in 1) Basic science developments with high potential for clinical translation 2) Clinical applications of cutting edge engineering and physics innovations 3) Broadly applicable and innovative clinical physics developments
Medical Physics is a journal of global scope and reach. By publishing in Medical Physics your research will reach an international, multidisciplinary audience including practicing medical physicists as well as physics- and engineering based translational scientists. We work closely with authors of promising articles to improve their quality.