Technical Note: Integrating an open source Monte Carlo code “MCsquare” for clinical use in intensity-modulated proton therapy Journal Article


Authors: Deng, W.; Younkin, J. E.; Souris, K.; Huang, S.; Augustine, K.; Fatyga, M.; Ding, X.; Cohilis, M.; Bues, M.; Shan, J.; Stoker, J.; Lin, L.; Shen, J.; Liu, W.
Article Title: Technical Note: Integrating an open source Monte Carlo code “MCsquare” for clinical use in intensity-modulated proton therapy
Abstract: Purpose: To commission an open source Monte Carlo (MC) dose engine, “MCsquare” for a synchrotron-based proton machine, integrate it into our in-house C++-based I/O user interface and our web-based software platform, expand its functionalities, and improve calculation efficiency for intensity-modulated proton therapy (IMPT). Methods: We commissioned MCsquare using a double Gaussian beam model based on in-air lateral profiles, integrated depth dose of 97 beam energies, and measurements of various spread-out Bragg peaks (SOBPs). Then we integrated MCsquare into our C++-based dose calculation code and web-based second check platform “DOSeCHECK.” We validated the commissioned MCsquare based on 12 different patient geometries and compared the dose calculation with a well-benchmarked GPU-accelerated MC (gMC) dose engine. We further improved the MCsquare efficiency by employing the computed tomography (CT) resampling approach. We also expanded its functionality by adding a linear energy transfer (LET)-related model-dependent biological dose calculation. Results: Differences between MCsquare calculations and SOBP measurements were <2.5% (<1.5% for ~85% of measurements) in water. The dose distributions calculated using MCsquare agreed well with the results calculated using gMC in patient geometries. The average 3D gamma analysis (2%/2 mm) passing rates comparing MCsquare and gMC calculations in the 12 patient geometries were 98.0 ± 1.0%. The computation time to calculate one IMPT plan in patients’ geometries using an inexpensive CPU workstation (Intel Xeon E5-2680 2.50 GHz) was 2.3 ± 1.8 min after the variable resolution technique was adopted. All calculations except for one craniospinal patient were finished within 3.5 min. Conclusions: MCsquare was successfully commissioned for a synchrotron-based proton beam therapy delivery system and integrated into our web-based second check platform. After adopting CT resampling and implementing LET model-dependent biological dose calculation capabilities, MCsquare will be sufficiently efficient and powerful to achieve Monte Carlo-based and LET-guided robust optimization in IMPT, which will be done in the future studies. © 2020 American Association of Physicists in Medicine
Keywords: intensity modulated radiation therapy; computerized tomography; geometry; energy transfer; software; radiation dose distribution; dose calculation; monte carlo method; optimization; proton therapy; gaussian beams; linear energy transfer; user interfaces; monte carlo methods; websites; open source software; human; article; proton beams; spread out bragg peaks; robust optimization; proton beam therapy; open systems; engines; intensity-modulated proton therapy (impt); mcsquare; c++ (programming language); calculation efficiency; intensity modulated proton therapies; re-sampling approach; variable resolution; web-based softwares
Journal Title: Medical Physics
Volume: 47
Issue: 6
ISSN: 0094-2405
Publisher: American Association of Physicists in Medicine  
Date Published: 2020-06-01
Start Page: 2558
End Page: 2574
Language: English
DOI: 10.1002/mp.14125
PUBMED: 32153029
PROVIDER: scopus
DOI/URL:
Notes: Article -- Export Date: 3 August 2020 -- Source: Scopus
Altmetric
Citation Impact
BMJ Impact Analytics
MSK Authors
  1. Sheng Huang
    13 Huang