A fast GPU-accelerated Monte Carlo engine for calculation of MLC-collimated electron fields.

Although intensity-modulated radiation therapy and volumetric arc therapy have revolutionized photon external beam therapies, the technological advances associated with electron beam therapy have fallen behind. Modern linear accelerators contain technologies that would allow for more advanced forms of electron treatments, such as beam collimation, using the conventional photon multi-leaf collimator (MLC); however, no commercial solutions exist that calculate dose from such beam delivery modes. Additionally, for clinical adoption to occur, dose calculation times would need to be on par with that of modern dose calculation algorithms. This work developed a graphics processing unit (GPU)-accelerated Monte Carlo (MC) engine incorporating the Varian TrueBeam linac head geometry for a rapid calculation of electron beams collimated using the conventional photon MLC. A compute unified device architecture framework was created for the following: (1) transport of electrons and photons through the linac head geometry, considering multiple scattering, Bremsstrahlung, Møller, Compton, and pair production interactions; (2) electron and photon propagation through the CT geometry, considering all interactions plus the photoelectric effect; and (3) secondary particle cascades through the linac head and within the CT geometry. The linac head collimating geometry was modeled according to the specifications provided by the vendor, who also provided phase-space files. The MC was benchmarked against EGSnrc/DOSXYZnrc/GEANT by simulating individual interactions with simple geometries, pencil, and square beam dose calculations in various phantoms. MC-calculated dose distributions for MLC and jaw-collimated electron fields were compared to measurements in a water phantom and with radiochromic film. Pencil and square beam dose distributions are in good agreement with DOSXYZnrc. Angular and spatial distributions for multiple scattering and secondary particle production in thin slab geometries are in good agreement with EGSnrc and GEANT. Dose profiles for MLC and jaw-collimated 6-20-MeV electron beams showed an average absolute difference of 1.1 and 1.9 mm for the FWHM and 80%-20% penumbra from measured profiles. Percent depth doses showed differences of <5% for as compared to measurement. The computation time on an NVIDIA Tesla V100 card was 2.5 min to achieve a dose uncertainty of <1%, which is ∼300 times faster than published results in a similar geometry using a single-CPU core. The GPU-based MC can quickly calculate dose for electron fields collimated using the conventional photon MLC. The fast calculation times will allow for a rapid calculation of electron fields for mixed photon and electron particle therapy.

Brost EE ，Wan Chan Tseung H ，Antolak JA 《-》

Validation of a Monte Carlo model for multi leaf collimator based electron delivery.

To develop and validate a Monte Carlo model of the Varian TrueBeam to study electron collimation using the existing photon multi-leaf collimators (pMLC), instead of conventional electron applicators and apertures. A complete Monte Carlo model of the Varian TrueBeam was developed using Tool for particle simulation (TOPAS) (version 3.1.p3). Vendor-supplied information was used to model the treatment head components and the source parameters. A phase space plane was setup above the collimating jaws and captured particles were reused until a statistical uncertainty of 1% was achieved in the central axis. Electron energies 6, 9, 12, 16, and 20 MeV with a jaw-defined field of 20 × 20 cm2 at iso-center, pMLC-defined fields of 6.8 × 6.8 cm2 and 11.4 × 11.4 cm2 at 80 cm source-to-surface distance (SSD) and an applicator-defined field of 10 × 10 cm2 at iso-center were evaluated. All the measurements except the applicator-defined fields were measured using an ionization chamber in a water tank using 80 cm SSD. The dose difference, distance-to-agreement and gamma index were used to evaluate the agreement between the Monte Carlo calculations and measurements. Contributions of electron scattering off pMLC leaves and inter-leaf leakage on dose profiles were evaluated and compared with Monte Carlo calculations. Electron transport through a heterogeneous phantom was simulated and the resulting dose distributions were compared with film measurements. The validated Monte Carlo model was used to simulate several clinically motivated cases to demonstrate the benefit of pMLC-based electron delivery compared to applicator-based electron delivery. Calculated and measured percentage depth-dose (PDD) curves agree within 2% after normalization. The agreement between normalized percentage depth dose curves were evaluated using one-dimensional gamma analysis with a local tolerance of 2%/1 mm and the %points passing gamma criteria was 100% for all energies. For jaw-defined fields, calculated profiles agree with measurements with pass rates of >97% for 2%/2 mm gamma criteria. Calculated FWHM and penumbra width agree with measurements within 0.4 cm. For fields with tertiary collimation using an pMLC or applicator, the average gamma pass rate of compared profiles was 98% with 2%/2 mm gamma criteria. The profiles measured to evaluate the pMLC leaf scattering agreed with Monte Carlo calculations with an average gamma pass rate of 96.5% with 3%/2 mm gamma criteria. Measured dose profiles below the heterogenous phantom agreed well with calculated profiles and matched within 2.5% for most points. The calculated clinically applicable cases using TOPAS MC and Eclipse TPS for single enface electron beam, electron-photon mixed beam and a matched electron-electron beam exhibited a reasonable agreement in PDDs, profiles and dose volume histograms. We present a validation of a Monte Carlo model of Varian TrueBeam for pMLC-based electron delivery. Monte Carlo calculations agreed with measurements satisfying gamma criterion of 1%/1 mm for depth dose curves and 2%/1 mm for dose profiles. The simulation of clinically applicable cases demonstrated the clinical utility of pMLC-based electrons and the use of MC simulations for development of advanced radiation therapy techniques.

Kaluarachchi MM ，Saleh ZH ，Schwer ML ，Klein EE ... -  《-》

New capabilities of the Monte Carlo dose engine ARCHER-RT: Clinical validation of the Varian TrueBeam machine for VMAT external beam radiotherapy.

The Monte Carlo radiation transport method is considered the most accurate approach for absorbed dose calculations in external beam radiation therapy. In this study, an efficient and accurate source model of the Varian TrueBeam 6X STx Linac is developed and integrated with a fast Monte Carlo photon-electron transport absorbed dose engine, ARCHER-RT, which is capable of being executed on CPUs, NVIDIA GPUs, and AMD GPUs. This capability of fast yet accurate radiation dose calculation is essential for clinical utility of this new technology. This paper describes the software and algorithmic developments made to the ARCHER-RT absorbed dose engine. AMD's Heterogeneous-Compute Interface for Portability (HIP) was implemented in ARCHER-RT to allow for device independent execution on NVIDIA and AMD GPUs. Architecture-specific atomic-add algorithms have been identified and both more accurate single-precision and double-precision computational absorbed dose calculation methods have been added to ARCHER-RT and validated through a test case to evaluate the accuracy and performance of the algorithms. The validity of the source model and the radiation transport physics were benchmarked against Monte Carlo simulations performed with EGSnrc. Secondary dose-check physics plans, and a clinical prostate treatment plan were calculated to demonstrate the applicability of the platform for clinical use. Absorbed dose difference maps and gamma analyses were conducted to establish the accuracy and consistency between the two Monte Carlo models. Timing studies were conducted on a CPU, an NVIDIA GPU, and an AMD GPU to evaluate the computational speed of ARCHER-RT. Percent depth doses were computed for different field sizes ranging from 1.5 cm2  × 1.5 cm2 to 22 cm2  × 40cm2 and the two codes agreed for all points outside high gradient regions within 3%. Axial profiles computed for a 10 cm2  × 10 cm2 field for multiple depths agreed for all points outside high gradient regions within 2%. The test case investigating the impact of native single-precision compared to double-precision showed differences in voxels as large as 71.47% and the implementation of KAS single-precision reduced the difference to less than 0.01%. The 3%/3mm gamma pass rates for an MPPG5a multileaf collimator (MLC) test case and a clinical VMAT prostate plan were 94.2% and 98.4% respectively. Timing studies demonstrated the calculation of a VMAT plan was completed in 50.3, 187.9, and 216.8 s on an NVIDIA GPU, AMD GPU, and Intel CPU, respectively. ARCHER-RT is capable of patient-specific VMAT external beam photon absorbed dose calculations and its potential has been demonstrated by benchmarking against a well validated EGSnrc model of a Varian TrueBeam. Additionally, the implementation of AMD's HIP has shown the flexibility of the ARCHER-RT platform for device independent calculations. This work demonstrates the significant addition of functionality added to ARCHER-RT framework which has marked utility for both research and clinical applications and demonstrates further that Monte Carlo-based absorbed dose engines like ARCHER-RT have the potential for widespread clinical implementation.

Adam DP ，Liu T ，Caracappa PF ，Bednarz BP ，Xu XG ... -  《-》

A GPU OpenCL based cross-platform Monte Carlo dose calculation engine (goMC).

Tian Z ，Shi F ，Folkerts M ，Qin N ，Jiang SB ，Jia X ... -  《-》

ARCHERRT - a GPU-based and photon-electron coupled Monte Carlo dose computing engine for radiation therapy: software development and application to helical tomotherapy.

Su L ，Yang Y ，Bednarz B ，Sterpin E ，Du X ，Liu T ，Ji W ，Xu XG ... -  《-》