Bragg peak

Computational science

Engineering

Statistics

Computational physics

Proton

Nuclear physics

Operating system

Intensity modulation

Medical physics

Dosimetry

Parallel computing

Computer science

Sobp

Statistical physics

Phase noise

Physics

Xeon

Code of practice

Mathematics

Optics

Phase modulation

Monte Carlo method

Medicine

Proton therapy

Nuclear medicine

Intensity (physics)

Engineering ethics

Medical Physics

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 &lt;2.5% (&lt;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.

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