Fast CPU-based Monte Carlo simulation for radiotherapy dose calculation

Performance assessment of variant UNet-based deep-learning dose engines for MR-Linac-based prostate IMRT plans

Objective. UNet-based deep-learning (DL) architectures are promising dose engines for traditional linear accelerator (Linac) models. Current UNet-based engines, however, were designed differently with various strategies, making it challenging to fairly compare the results from different studies. The objective of this study is to thoroughly evaluate the performance of UNet-based models on magnetic-resonance (MR)-Linac-based intensity-modulated radiation therapy (IMRT) dose calculations.Approach. The UNet-based models, including the standard-UNet, cascaded-UNet, dense-dilated-UNet, residual-UNet, HD-UNet, and attention-aware-UNet, were implemented. The model input is patient CT and IMRT field dose in water, and the output is patient dose calculated by DL model. The reference dose was calculated by the Monaco Monte Carlo module. Twenty training and ten test cases of prostate patients were included. The accuracy of the DL-calculated doses was measured using gamma analysis, and the calculation efficiency was evaluated by inference time.Results. All the studied models effectively corrected low-accuracy doses in water to high-accuracy patient doses in a magnetic field. The gamma passing rates between reference and DL-calculated doses were over 86% (1%/1 mm), 98% (2%/2 mm), and 99% (3%/3 mm) for all the models. The inference times ranged from 0.03 (graphics processing unit) to 7.5 (central processing unit) seconds. Each model demonstrated different strengths in calculation accuracy and efficiency; Res-UNet achieved the highest accuracy, HD-UNet offered high accuracy with the fewest parameters but the longest inference, dense-dilated-UNet was consistently accurate regardless of model levels, standard-UNet had the shortest inference but relatively lower accuracy, and the others showed average performance. Therefore, the best-performing model would depend on the specific clinical needs and available computational resources.Significance. The feasibility of using common UNet-based models for MR-Linac-based dose calculations has been explored in this study. By using the same model input type, patient training data, and computing environment, a fair assessment of the models' performance was present.

Assessing the Effect of Directional Bremsstrahlung Splitting on the Output Spectra and Parameters Using BEAMnrc Monte Carlo Simulation Package

Introduction:

EGSnrc software package is one of the computational packages for Monte Carlo simulation in radiation therapy and has several subset codes. Directional bremsstrahlung splitting (DBS) is a technique that applies braking radiations in interactions in this software. This study aimed to evaluate the effect of this technique on the simulation time, uncertainty, particle number of phase-space data, and photon beam spectrum resulting from a medical linear accelerator (LINAC).

Materials and methods:

The gantry of the accelerator, including the materials and geometries of different parts, was simulated using the BEAMnrc code (a subset code in the EGSnrc package). The phase-space data were recorded in different parts of the LINAC. The DBS values (1, 10, 100, and 1000) were changed, and their effects were evaluated on the simulation parameters and output spectra.

Results:

Increasing the DBS value from 1 to 1000 resulted in an increase in the simulation time from 1.778 to 11.310 hours, and increasing the number of particles in the phase-space plane (5 590 732-180 328 382). When the DBS had been picked up from 1 to 100, the simulation uncertainty decreased by about 1.29%. In addition, the DBS increment value from 100 to 1000 leads to an increase in uncertainty and simulation time of about 0.71% and 315%, respectively.

Conclusion:

Although using the DBS technique reduces the simulation time or uncertainty, increasing the DBS from a specific value, equal to 100 in our study, increases simulation uncertainties and times. Therefore, we propose considering a specific DBS value as we obtained for the Monte Carlo simulation of photon beams produced by linear accelerators.

Determination of dose distributions by monte-carlo simulation of 6 MV photon beam of varian vitalbeam accelerator using geant4 multithreaded code

Background:

Accuracy of dose delivery in radiation therapy is a primary requirement for effective cancer treatment. In practice, dose delivery accuracy of ±5% is desired. To achieve this accuracy level, an accurate method for calculating the dose distributions in the tumor volume is required. Monte-Carlo method is one of the methods considered to be the most accurate for calculating dose distributions.

Materials and Methods:

G4 linac-MT code was used to simulate a 6 MV photon beam. The initial electron beam parameters were tuned to validate the beam modeling from depth doses and beam profile. The dose distributions measured in water phantom were compared to the calculated dose distributions based on gamma index criterion.

Results:

The beam tuning showed the initial electron energy, sigma and full width at half maximum of 6.2 MeV, 0.8 MeV, and 1.18 mm, respectively, best match the measured dose distributions. The gamma index tests showed the calculated depth doses and beam profile were generally comparable with measurements, passing the standard acceptance criterion of 2%/2 mm. The simulated photon beam was justified by the index of beam quality, which showed excellent agreement with measured doses with a discrepancy of 0.1%.

Conclusion:

The observed agreement confirm the accuracy of the simulated 6 MV photon beam. It can therefore be used as radiation source for calculating dose distributions and further investigations aimed at improving dose delivery and planning in cancer patients.

XIORT-MC: A real-time MC-based dose computation tool for low- energy X-rays intraoperative radiation therapy

PURPOSE

The INTRABEAM system is a miniature accelerator for low-energy X-ray Intra-Operative Radiation Therapy (IORT), and it could benefit from a fast and accurate dose computation tool. With regards to accuracy, dose computed with Monte Carlo (MC) simulations are the gold standard, however, they require a large computational effort and consequently they are not suitable for real-time dose planning. This work presents a comparison of the implementation on Graphics Processing Unit (GPU) of two different dose calculation algorithms based on MC phase-space (PHSP) information to compute dose distributions for the INTRABEAM device within seconds and with the accuracy of realistic MC simulations.

METHODS

The MC-based algorithms we present incorporate photoelectric, Compton and Rayleigh effects for the interaction of low-energy X-rays. XIORT-MC (X-ray Intra-Operative Radiation Therapy Monte Carlo) includes two dose calculation algorithms; a Woodcock-based MC algorithm (WC-MC) and a Hybrid MC algorithm (HMC), and it is implemented in CPU and in GPU. Detailed MC simulations have been generated to validate our tool in homogeneous and heterogeneous conditions with all INTRABEAM applicators, including three clinically realistic CT-based simulations. A performance study has been done to determine the acceleration reached with the code, in both CPU and GPU implementations.

RESULTS

Dose distributions were obtained with the HMC and the WC-MC and compared to standard reference MC simulations with more than 95% voxels fulfilling a 7%-0.5 mm gamma evaluation in all the cases considered. The CPU-HMC is 100 times more efficient than the reference MC, and the CPU-WC-MC is about 50 times more efficient. With the GPU implementation, the particle tracking of the WC-MC is faster than the HMC, with the extraction of the particle's information from the PHSP file taking a major part of the time. However, thanks to the variance reduction techniques implemented in the HMC, up to 400 times less particles are needed in the HMC to reach the same level of noise than the WC-MC. Therefore, in our implementation for INTRABEAM energies, the HMC is about 1.3 times more efficient than the WC-MC in an NVIDIA GeForce GTX 1080 Ti card and about 5.5 times more efficient in an NVIDIA GeForce RTX 3090. Dose with noise below 5% has been obtained in realistic situations in less than 5 s with the WC-MC and in less than 0.5 s with the HMC.

CONCLUSIONS

The XIORT-MC is a dose computation tool designed to take full advantage of modern GPUs, making possible to obtain MC-grade accurate dose distributions within seconds. Its high speed allows a real-time dose calculation that includes the realistic effects of the beam in voxelized geometries of patients. It can be used as a dose-planning tool in the operating room during a XIORT treatment with any INTRABEAM device.

Characterization of Extrafocal Dose Influence on the Out-of-Field Dose Distribution by Monte Carlo Simulations and Dose Measurements

Out-of-field scattered and transmitted extrafocal radiation may induce secondary cancer in long-term survivors of external radiotherapy. Pediatric patients have higher life expectancy and tend to receive higher secondary radiation damage due to geometric and biological factors. The goal of this study is to characterize the location and the magnitude of extrafocal dose regions in the case of three-dimensional conformal radiotherapy and volumetric arc therapy, to apply this information to clinical treatment cases, and to provide mitigation strategies. Extrafocal dose has been investigated in a Varian TrueBeam linac equipped with a high-definition 120 multileaf collimator using different physical and virtual phantoms, dose calculation (including Monte Carlo techniques), and dose measurement methods. All Monte Carlo calculations showed excellent agreement with measurements. Treatment planning system calculations failed to provide reliable results out of the treatment field. Both Monte Carlo calculations and dose measurements showed regions with higher dose (extrafocal dose areas) when compared to the background. These areas start to be noticeable beyond 11 cm from the isocenter in the direction perpendicular to the multileaf collimator leaves' travel direction. Out-of-field extrafocal doses up to 160% of the mean dose transmitted through the closed multileaf collimator were registered. Two overlapping components were observed in the extrafocal distribution: the first is an almost elliptical blurred dose distribution, and the second is a well-defined rectangular dose distribution. Extra precautions should be taken into consideration when treating pediatric patients with a high-definition 120 multileaf collimator to avoid directing the extrafocal radiation into a radiosensitive organ during external beam therapy.

SIMU-RAD programme: a learning tool for radiation (photons and charged particles) interaction

Radiation education is necessary for a wide variety of people, such as radiation workers particularly for students of secondary school and higher education institution who learn radiation sciences. The fact that we could not see or feel radiation makes it difficult to understand it. The use of radiation trajectories shown on a personal computer should be useful to overcome this difficulty. In order to understand radiation behaviour inside the material, we have developed a Simu-Rad (Copyright: LY2018002738) by using Monte Carlo simulation programme. One who has no programming knowledge is able to simulate photons in a material through the developed programme. The program could become a computer aided learning tool for radiation related courses. We aim to facilitate lecturer from ‘The Traditional Classroom’ to ‘The Flipped Classroom’ for radiation education concerning in the era of IR 4.0. To validate our radiation simulator, we calculate photon linear attenuation coefficient (µ) of an aluminium material which commonly used as a filter in diagnostic radiology. µ is one of the main characteristics to understand how the radiation attenuated inside the materials. We calculate at energy photon of 662 keV (Cs-137 radiation source) to compare our results of µ with the XCOM database. Consequently, the results from the developed simulator comparable with the database verified our programme to be used for radiation study.

Technical Note: An hybrid parallel implementation for EGSnrc Monte Carlo user codes

PURPOSE

The purpose of this study was to present a parallel solution for the EGSnrc Monte Carlo code system combining MPI and OpenMP programming models as an alternative to the provided implementation, based on the use of a batch-queueing system (BQS).

METHODS

Relying on a previous implementation based on OpenMP by E. Doerner and P. Caprile [Med. Phys. 44, 6672 (2017)], this work incorporates MPI features to efficiently distribute the simulation on current high-performance computing (HPC) systems. These features are introduced through properly defined macros, which are enabled depending on the compilation flags given by the user. The presented solution was benchmarked using the DOSXYZnrc code for a 6 MV clinical photon beam impinging on an homogeneous water phantom.

RESULTS

The platform validation against the serial run results confirmed that the introduction of new features does not modify the final dose distribution. The performance tests indicated that the new implementation was able to handle efficiently the workload distribution among the computing units available. Using all the resources available, the hybrid simulation was 10% faster than the MPI only solution and 30% faster than the BQS implementation.

CONCLUSIONS

The hybrid method presented is a viable solution to parallelize MC simulations using the EGSnrc codes in distributed computing systems in an simple and efficient way, taking advantage of the available resources and giving the user the possibility of choosing between different parallelization schemes (only OpenMP/MPI or a combination of both).

Performance of a hybrid Monte Carlo-Pencil Beam dose algorithm for proton therapy inverse planning

PURPOSE

Analytical algorithms have a limited accuracy when modeling very heterogeneous tumor sites. This work addresses the performance of a hybrid dose optimizer that combines both Monte Carlo (MC) and pencil beam (PB) dose engines to get the best trade-off between speed and accuracy for proton therapy plans.

METHODS

The hybrid algorithm calculates the optimal spot weights (w) by means of an iterative optimization process where the dose at each iteration is computed by using a precomputed dose influence matrix based on the conventional PB plus a correction term c obtained from a MC simulation. Updates of c can be triggered as often as necessary by calling the MC dose engine with the last corrected values of w as input. In order to analyze the performance of the hybrid algorithm against dose calculation errors, it was applied to a simplistic water phantom for which several test cases with different errors were simulated, including proton range uncertainties. Afterwards, the algorithm was used in three clinical cases (prostate, lung, and brain) and benchmarked against full MC-based optimization. The influence of different stopping criteria in the final results was also investigated.

RESULTS

The hybrid algorithm achieved excellent results provided that the estimated range in a homogeneous material is the same for the two dose engines involved, i.e., PB and MC. For the three patient cases, the hybrid plans were clinically equivalent to those obtained with full MC-based optimization. Only a single update of c was needed in the hybrid algorithm to fulfill the clinical dose constraints, which represents an extra computation time to obtain c that ranged from 1 (brain) to 4 min (lung) with respect to the conventional PB-based optimization, and an estimated average gain factor of 14 with respect to full MC-based optimization.

CONCLUSION

The hybrid algorithm provides an improved trade-off between accuracy and speed. This algorithm can be immediately considered as an option for improving dose calculation accuracy of commercial analytical treatment planning systems, without a significant increase in the computation time (≪5 min) with respect to current PB-based optimization.

MRI-guided lung SBRT: Present and future developments

Stereotactic body radiotherapy (SBRT) is rapidly becoming an alternative to surgery for the treatment of early-stage non-small cell lung cancer patients. Lung SBRT is administered in a hypo-fractionated, conformal manner, delivering high doses to the target. To avoid normal-tissue toxicity, it is crucial to limit the exposure of nearby healthy organs-at-risk (OAR). Current image-guided radiotherapy strategies for lung SBRT are mostly based on X-ray imaging modalities. Although still in its infancy, magnetic resonance imaging (MRI) guidance for lung SBRT is not exposure-limited and MRI promises to improve crucial soft-tissue contrast. Looking beyond anatomical imaging, functional MRI is expected to inform treatment decisions and adaptations in the future. This review summarises and discusses how MRI could be advantageous to the different links of the radiotherapy treatment chain for lung SBRT: diagnosis and staging, tumour and OAR delineation, treatment planning, and inter- or intrafractional motion management. Special emphasis is placed on a new generation of hybrid MRI treatment devices and their potential for real-time adaptive radiotherapy.

Towards real-time photon Monte Carlo dose calculation in the cloud

Near real-time application of Monte Carlo (MC) dose calculation in clinic and research is hindered by the long computational runtimes of established software. Currently, fast MC software solutions are available utilising accelerators such as graphical processing units (GPUs) or clusters based on central processing units (CPUs). Both platforms are expensive in terms of purchase costs and maintenance and, in case of the GPU, provide only limited scalability. In this work we propose a cloud-based MC solution, which offers high scalability of accurate photon dose calculations. The MC simulations run on a private virtual supercomputer that is formed in the cloud. Computational resources can be provisioned dynamically at low cost without upfront investment in expensive hardware. A client-server software solution has been developed which controls the simulations and transports data to and from the cloud efficiently and securely. The client application integrates seamlessly into a treatment planning system. It runs the MC simulation workflow automatically and securely exchanges simulation data with the server side application that controls the virtual supercomputer. Advanced encryption standards were used to add an additional security layer, which encrypts and decrypts patient data on-the-fly at the processor register level. We could show that our cloud-based MC framework enables near real-time dose computation. It delivers excellent linear scaling for high-resolution datasets with absolute runtimes of 1.1 seconds to 10.9 seconds for simulating a clinical prostate and liver case up to 1% statistical uncertainty. The computation runtimes include the transportation of data to and from the cloud as well as process scheduling and synchronisation overhead. Cloud-based MC simulations offer a fast, affordable and easily accessible alternative for near real-time accurate dose calculations to currently used GPU or cluster solutions.

Expert system classifier for adaptive radiation therapy in prostate cancer

A classifier-based expert system was developed to compare delivered and planned radiation therapy in prostate cancer patients. Its aim is to automatically identify patients that can benefit from an adaptive treatment strategy. The study predominantly addresses dosimetric uncertainties and critical issues caused by motion of hollow organs. 1200 MVCT images of 38 prostate adenocarcinoma cases were analyzed. An automatic daily re-contouring of structures (i.e. rectum, bladder and femoral heads), rigid/deformable registration and dose warping was carried out to simulate dose and volume variations during therapy. Support vector machine, K-means clustering algorithms and similarity index analysis were used to create an unsupervised predictive tool to detect incorrect setup and/or morphological changes as a consequence of inadequate patient preparation due to stochastic physiological changes, supporting clinical decision-making. After training on a dataset that was considered sufficiently dosimetrically stable, the system identified two equally sized macro clusters with distinctly different volumetric and dosimetric baseline properties and defined thresholds for these two clusters. Application to the test cohort resulted in 25% of the patients located outside the two macro clusters thresholds and which were therefore suspected to be dosimetrically unstable. In these patients, over the treatment course, mean volumetric changes of 30 and 40% for rectum and bladder were detected which possibly represents values justifying adjustment of patient preparation, frequent re-planning or a plan-of-the-day strategy. Based on our research, by combining daily IGRT images with rigid/deformable registration and dose warping, it is possible to apply a machine learning approach to the clinical setting obtaining useful information for a decision regarding an individualized adaptive strategy. Especially for treatments influenced by the movement of hollow organs, this could reduce inadequate treatments and possibly reduce toxicity, thereby increasing overall RT efficacy.

Fast multipurpose Monte Carlo simulation for proton therapy using multi- and many-core CPU architectures

PURPOSE

Accuracy in proton therapy treatment planning can be improved using Monte Carlo (MC) simulations. However the long computation time of such methods hinders their use in clinical routine. This work aims to develop a fast multipurpose Monte Carlo simulation tool for proton therapy using massively parallel central processing unit (CPU) architectures.

METHODS

A new Monte Carlo, called MCsquare (many-core Monte Carlo), has been designed and optimized for the last generation of Intel Xeon processors and Intel Xeon Phi coprocessors. These massively parallel architectures offer the flexibility and the computational power suitable to MC methods. The class-II condensed history algorithm of MCsquare provides a fast and yet accurate method of simulating heavy charged particles such as protons, deuterons, and alphas inside voxelized geometries. Hard ionizations, with energy losses above a user-specified threshold, are simulated individually while soft events are regrouped in a multiple scattering theory. Elastic and inelastic nuclear interactions are sampled from ICRU 63 differential cross sections, thereby allowing for the computation of prompt gamma emission profiles. MCsquare has been benchmarked with the gate/geant4 Monte Carlo application for homogeneous and heterogeneous geometries.

RESULTS

Comparisons with gate/geant4 for various geometries show deviations within 2%-1 mm. In spite of the limited memory bandwidth of the coprocessor simulation time is below 25 s for 10(7) primary 200 MeV protons in average soft tissues using all Xeon Phi and CPU resources embedded in a single desktop unit.

CONCLUSIONS

MCsquare exploits the flexibility of CPU architectures to provide a multipurpose MC simulation tool. Optimized code enables the use of accurate MC calculation within a reasonable computation time, adequate for clinical practice. MCsquare also simulates prompt gamma emission and can thus be used also for in vivo range verification.

Real-time 4D dose reconstruction for tracked dynamic MLC deliveries for lung SBRT

PURPOSE

This study provides a proof of concept for real-time 4D dose reconstruction for lung stereotactic body radiation therapy (SBRT) with multileaf collimator (MLC) tracking and assesses the impact of tumor tracking on the size of target margins.

METHODS

The authors have implemented real-time 4D dose reconstruction by connecting their tracking and delivery software to an Agility MLC at an Elekta Synergy linac and to their in-house treatment planning software (TPS). Actual MLC apertures and (simulated) target positions are reported to the TPS every 40 ms. The dose is calculated in real-time from 4DCT data directly after each reported aperture by utilization of precalculated dose-influence data based on a Monte Carlo algorithm. The dose is accumulated onto the peak-exhale (reference) phase using energy-mass transfer mapping. To investigate the impact of a potentially reducible safety margin, the authors have created and delivered treatment plans designed for a conventional internal target volume (ITV) + 5 mm, a midventilation approach, and three tracking scenarios for four lung SBRT patients. For the tracking plans, a moving target volume (MTV) was established by delineating the gross target volume (GTV) on every 4DCT phase. These were rigidly aligned to the reference phase, resulting in a unified maximum GTV to which a 1, 3, or 5 mm isotropic margin was added. All scenarios were planned for 9-beam step-and-shoot IMRT to meet the criteria of RTOG 1021 (3 × 18 Gy). The GTV 3D center-of-volume shift varied from 6 to 14 mm.

RESULTS

Real-time dose reconstruction at 25 Hz could be realized on a single workstation due to the highly efficient implementation of dose calculation and dose accumulation. Decreased PTV margins resulted in inadequate target coverage during untracked deliveries for patients with substantial tumor motion. MLC tracking could ensure the GTV target dose for these patients. Organ-at-risk (OAR) doses were consistently reduced by decreased PTV margins. The tracked MTV + 1 mm deliveries resulted in the following OAR dose reductions: lung V20 up to 3.5%, spinal cord D2 up to 0.9 Gy/Fx, and proximal airways D2 up to 1.4 Gy/Fx.

CONCLUSIONS

The authors could show that for patient data at clinical resolution and realistic motion conditions, the delivered dose could be reconstructed in 4D for the whole lung volume in real-time. The dose distributions show that reduced margins yield lower doses to healthy tissue, whilst target dose can be maintained using dynamic MLC tracking.

A machine learning tool for re-planning and adaptive RT: A multicenter cohort investigation

PURPOSE

To predict patients who would benefit from adaptive radiotherapy (ART) and re-planning intervention based on machine learning from anatomical and dosimetric variations in a retrospective dataset.

MATERIALS AND METHODS

90 patients (pts) treated for head-neck cancer (H&N) formed a multicenter data-set. 41 H&N pts (45.6%) were considered for learning; 49 pts (54.4%) were used to test the tool. A homemade machine-learning classifier was developed to analyze volume and dose variations of parotid glands (PG). Using deformable image registration (DIR) and GPU, patients' conditions were analyzed automatically. Support Vector Machines (SVM) was used for time-series evaluation. "Inadequate" class identified patients that might benefit from replanning. Double-blind evaluation by two radiation oncologists (ROs) was carried out to validate day/week selected for re-planning by the classifier.

RESULTS

The cohort was affected by PG mean reduction of 23.7±8.8%. During the first 3weeks, 86.7% cases show PG deformation aligned with predefined tolerance, thus not requiring re-planning. From 4th week, an increased number of pts would potentially benefit from re-planning: a mean of 58% of cases, with an inter-center variability of 8.3%, showed "inadequate" conditions. 11% of cases showed "bias" due to DIR and script failure; 6% showed "warning" output due to potential positioning issues. Comparing re-planning suggested by tool with recommended by ROs, the 4th week seems the most favorable time in 70% cases.

CONCLUSIONS

SVM and decision-making tool was applied to overcome ART challenges. Pts would benefit from ART and ideal time for re-planning intervention was identified in this retrospective analysis.

Assessment of MLC tracking performance during hypofractionated prostate radiotherapy using real-time dose reconstruction

By adapting to the actual patient anatomy during treatment, tracked multi-leaf collimator (MLC) treatment deliveries offer an opportunity for margin reduction and healthy tissue sparing. This is assumed to be especially relevant for hypofractionated protocols in which intrafractional motion does not easily average out. In order to confidently deliver tracked treatments with potentially reduced margins, it is necessary to monitor not only the patient anatomy but also the actually delivered dose during irradiation. In this study, we present a novel real-time online dose reconstruction tool which calculates actually delivered dose based on pre-calculated dose influence data in less than 10 ms at a rate of 25 Hz. Using this tool we investigate the impact of clinical target volume (CTV) to planning target volume (PTV) margins on CTV coverage and organ-at-risk dose. On our research linear accelerator, a set of four different CTV-to-PTV margins were tested for three patient cases subject to four different motion conditions. Based on this data, we can conclude that tracking eliminates dose cold spots which can occur in the CTV during conventional deliveries even for the smallest CTV-to-PTV margin of 1 mm. Changes of organ-at-risk dose do occur frequently during MLC tracking and are not negligible in some cases. Intrafractional dose reconstruction is expected to become an important element in any attempt of re-planning the treatment plan during the delivery based on the observed anatomy of the day.