Fast Monte Carlo simulation for patient-specific CT/CBCT imaging dose calculation

Technical note: A GPU-based shared Monte Carlo method for fast photon transport in multi-energy x-ray exposures

BACKGROUND

The Monte Carlo (MC) method is an accurate technique for particle transport calculation due to the precise modeling of physical interactions. Nevertheless, the MC method still suffers from the problem of expensive computational cost, even with graphics processing unit (GPU) acceleration. Our previous works have investigated the acceleration strategies of photon transport simulation for single-energy CT. But for multi-energy CT, conventional individual simulation leads to unnecessary redundant calculation, consuming more time.

PURPOSE

This work proposes a novel GPU-based shared MC scheme (gSMC) to reduce unnecessary repeated simulations of similar photons between different spectra, thereby enhancing the efficiency of scatter estimation in multi-energy x-ray exposures.

METHODS

The shared MC method selects shared photons between different spectra using two strategies. Specifically, we introduce spectral region classification strategy to select photons with the same initial energy from different spectra, thus generating energy-shared photon groups. Subsequently, the multi-directional sampling strategy is utilized to select energy-and-direction-shared photons, which have the same initial direction, from energy-shared photon groups. Energy-and-direction-shared photons perform shared simulations, while others are simulated individually. Finally, all results are integrated to obtain scatter distribution estimations for different spectral cases.

RESULTS

The efficiency and accuracy of the proposed gSMC are evaluated on the digital phantom and clinical case. The experimental results demonstrate that gSMC can speed up the simulation in the digital case by ∼37.8% and the one in the clinical case by ∼20.6%, while keeping the differences in total scatter results within 0.09%, compared to the conventional MC package, which performs an individual simulation.

CONCLUSIONS

The proposed GPU-based shared MC simulation method can achieve fast photon transport calculation for multi-energy x-ray exposures.

Monte Carlo methods for medical imaging research

In radiation-based medical imaging research, computational modeling methods are used to design and validate imaging systems and post-processing algorithms. Monte Carlo methods are widely used for the computational modeling as they can model the systems accurately and intuitively by sampling interactions between particles and imaging subject with known probability distributions. This article reviews the physics behind Monte Carlo methods, their applications in medical imaging, and available MC codes for medical imaging research. Additionally, potential research areas related to Monte Carlo for medical imaging are discussed.

Development of a TOPAS Monte Carlo (MC) model extension to simulate the automatic exposure control function of a C-arm CBCT

PURPOSE

Accurate simulation of organ doses in C-arm CBCT is critical for estimating personalised patient dosimetry. However, system complexities such as automatic exposure control (AEC) and the incorporation of DICOM images into simulations are challenging. The aim of this study was to develop a model for mimicking the operation of an AEC system, which maintains a constant dose to the detector through mA modulation in order to facilitate more accurate MC dosimetry models for C-arm CBCT.

METHODS

A Siemens Artis Q Interventional Radiology (IR) C-arm system [Siemens, Erlangen, Germany] was modelled in TOol for PArticle Simulation (TOPAS) by incorporating system specifications such as rotational speed, number of projections and exam protocol parameters. A novel threshold scorer, AECScorer, was developed to model the AEC functionality. MC simulations were performed using a variety of imaged volumes including a CTDI phantom, an anthropomorphic phantom and a patient DICOM dataset.

RESULTS

The AECScorer extension provides a framework for a conditional scoring function within TOPAS which allows for the simulation of an AEC system. The AECScorer successfully equalises the dose to the detector for simple phantoms and DICOM imaging datasets by adjusting the number of histories simulated at each CBCT projection. This AECSCorer tool is applicable to other medical imaging systems requiring AEC simulation.

CONCLUSIONS

We demonstrate a novel threshold scorer in TOPAS for a C-arm CBCT setup. The presented AECScorer is the first step towards providing a system-, patient- and protocol-specific dose estimates from CBCT dosimetry applications.

Dose, dose, dose, but where is the patient dose?

The article reviews the historical developments in radiation dose metrices in medical imaging. It identifies the good, the bad, and the ugly aspects of current-day metrices. The actions on shifting focus from International Commission on Radiological Protection (ICRP) Reference-Man-based population-average phantoms to patient-specific computational phantoms have been proposed and discussed. Technological developments in recent years involving AI-based automatic organ segmentation and 'near real-time' Monte Carlo dose calculations suggest the feasibility and advantage of obtaining patient-specific organ doses. It appears that the time for ICRP and other international organizations to embrace 'patient-specific' dose quantity representing risk may have finally come. While the existing dose metrices meet specific demands, emphasis needs to be also placed on making radiation units understandable to the medical community.

An indirect estimation of x-ray spectrum via convolutional neural network and transmission measurement

Objective.In this work, we aim to propose an accurate and robust spectrum estimation method by synergistically combining x-ray imaging physics with a convolutional neural network (CNN).Approach.The approach relies on transmission measurements, and the estimated spectrum is formulated as a convolutional summation of a few model spectra generated using Monte Carlo simulation. The difference between the actual and estimated projections is utilized as the loss function to train the network. We contrasted this approach with the weighted sums of model spectra approach previously proposed. Comprehensive studies were performed to demonstrate the robustness and accuracy of the proposed approach in various scenarios.Main results.The results show the desirable accuracy of the CNN-based method for spectrum estimation. The ME and NRMSE were -0.021 keV and 3.04% for 80 kVp, and 0.006 keV and 4.44% for 100 kVp, superior to the previous approach. The robustness test and experimental study also demonstrated superior performances. The CNN-based approach yielded remarkably consistent results in phantoms with various material combinations, and the CNN-based approach was robust concerning spectrum generators and calibration phantoms.Significance. We proposed a method for estimating the real spectrum by integrating a deep learning model with real imaging physics. The results demonstrated that this method was accurate and robust in estimating the spectrum, and it is potentially helpful for broad x-ray imaging tasks.

Organ dose prediction for patients undergoing radiotherapy CBCT chest examinations using artificial intelligence

PURPOSE

To propose an artificial intelligence (AI)-based method for personalized and real-time dosimetry for chest CBCT acquisitions.

METHODS

CT images from 113 patients who underwent radiotherapy treatment were collected for simulating thorax examinations using cone-beam computed tomography (CBCT) with the Monte Carlo technique. These simulations yielded organ dose data, used to train and validate specific AI algorithms. The efficacy of these AI algorithms was evaluated by comparing dose predictions with the actual doses derived from Monte Carlo simulations, which are the ground truth, utilizing Bland-Altman plots for this comparative analysis.

RESULTS

The absolute mean discrepancies between the predicted doses and the ground truth are (0.9 ± 1.3)% for bones, (1.2 ± 1.2)% for the esophagus, (0.5 ± 1.3)% for the breast, (2.5 ± 1.4)% for the heart, (2.4 ± 2.1)% for lungs, (0.8 ± 0.6)% for the skin, and (1.7 ± 0.7)% for integral. Meanwhile, the maximum discrepancies between the predicted doses and the ground truth are (14.4 ± 1.3)% for bones, (12.9 ± 1.2)% for the esophagus, (9.4 ± 1.3)% for the breast, (14.6 ± 1.4)% for the heart, (21.2 ± 2.1)% for lungs, (10.0 ± 0.6)% for the skin, and (10.5 ± 0.7)% for integral.

CONCLUSIONS

AI models that can make real-time predictions of the organ doses for patients undergoing CBCT thorax examinations as part of radiotherapy pre-treatment positioning were developed. The results of this study clearly show that the doses predicted by analyzed AI models are in close agreement with those calculated using Monte Carlo simulations.

Virtual dosimetry study with three cone-beam breast computed tomography scanners using a fast GPU-based Monte Carlo code

Objective. To compare the dosimetric performance of three cone-beam breast computed tomography (BCT) scanners, using real-time Monte Carlo-based dose estimates obtained with the virtual clinical trials (VCT)-BREAST graphical processing unit (GPU)-accelerated platform dedicated to VCT in breast imaging. Approach. A GPU-based Monte Carlo (MC) code was developed for replicatingin silicothe geometric, x-ray spectra and detector setups adopted, respectively, in two research scanners and one commercial BCT scanner, adopting 80 kV, 60 kV and 49 kV tube voltage, respectively. Our cohort of virtual breasts included 16 anthropomorphic voxelized breast phantoms from a publicly available dataset. For each virtual patient, we simulated exams on the three scanners, up to a nominal simulated mean glandular dose of 5 mGy (primary photons launched, in the order of 1011-1012per scan). Simulated 3D dose maps (recorded for skin, adipose and glandular tissues) were compared for the same phantom, on the three scanners. MC simulations were implemented on a single NVIDIA GeForce RTX 3090 graphics card.Main results.Using the spread of the dose distribution as a figure of merit, we showed that, in the investigated phantoms, the glandular dose is more uniform within less dense breasts, and it is more uniformly distributed for scans at 80 kV and 60 kV, than at 49 kV. A realistic virtual study of each breast phantom was completed in about 3.0 h with less than 1% statistical uncertainty, with 109primary photons processed in 3.6 s computing time.Significance. We reported the first dosimetric study of the VCT-BREAST platform, a fast MC simulation tool for real-time virtual dosimetry and imaging trials in BCT, investigating the dose delivery performance of three clinical BCT scanners. This tool can be adopted to investigate also the effects on the 3D dose distribution produced by changes in the geometrical and spectrum characteristics of a cone-beam BCT scanner.

Roadmap for precision preclinical x-ray radiation studies

This Roadmap paper covers the field of precision preclinical x-ray radiation studies in animal models. It is mostly focused on models for cancer and normal tissue response to radiation, but also discusses other disease models. The recent technological evolution in imaging, irradiation, dosimetry and monitoring that have empowered these kinds of studies is discussed, and many developments in the near future are outlined. Finally, clinical translation and reverse translation are discussed.

Small animal photon counting cone-beam CT on a preclinical radiation research platform to improve radiation dose calculation accuracy

Objective.Cone beam CT (CBCT) in preclinical small animal irradiation platforms provides essential information for image guidance and radiation dose calculation for experiment planning. This project developed a photon-counting detector (PCD)-based multi(3)-energy (ME-)CBCT on a small animal irradiator to improve the accuracy of material differentiation and hence dose calculation, and compared to conventional flat panel detector (FPD)-based CBCT.Approach.We constructed a mechanical structure to mount a PCD to an existing preclinical irradiator platform and built a data acquisition pipeline to acquire x-ray projection data with a 100 kVp x-ray beam using three different energy thresholds in a single gantry rotation. We implemented an energy threshold optimization scheme to determine optimal thresholds to balance signal-to-noise ratios (SNRs) among energy channels. Pixel-based detector response calibration was performed to remove ring artifacts in reconstructed CBCT images. Feldkamp-Davis-Kress method was employed to reconstruct CBCT images and a total-variance regularization-based optimization model was used to decompose CBCT images into bone and water material images. We compared dose calculation results using PCD-based ME-CBCT with that of FPD-based CBCT.Main results.The optimal nominal energy thresholds were determined as 26, 56, and 90 keV, under which SNRs in a selected region-of-interest in the water region were 6.11, 5.91 and 5.93 in the three energy channels, respectively. Compared with dose calculation results using FPD-based CBCT, using PCD-based ME-CBCT reduced the mean relative error from 49.5% to 16.4% in bone regions and from 7.5% to 6.9% in soft tissue regions.Significance.PCD-based ME-CBCT is beneficial in improving radiation dose calculation accuracy in experiment planning of preclinical small animal irradiation researches.

Real-time estimation of patient-specific dose distributions for medical CT using the deep dose estimation

PURPOSE

With the rising number of computed tomography (CT) examinations and the trend toward personalized medicine, patient-specific dose estimates are becoming more and more important in CT imaging. However, current approaches are often too slow or too inaccurate to be applied routinely. Therefore, we propose the so-called deep dose estimation (DDE) to provide highly accurate patient dose distributions in real time METHODS: To combine accuracy and computational performance, the DDE algorithm uses a deep convolutional neural network to predict patient dose distributions. To do so, a U-net like architecture is trained to reproduce Monte Carlo simulations from a two-channel input consisting of a CT reconstruction and a first-order dose estimate. Here, the corresponding training data were generated using CT simulations based on 45 whole-body patient scans. For each patient, simulations were performed for different anatomies (pelvis, abdomen, thorax, head), different tube voltages (80 kV, 100 kV, 120 kV), different scan trajectories (circle, spiral), and with and without bowtie filtration and tube current modulation. Similar simulations were performed using a second set of eight whole-body CT scans from the Visual Concept Extraction Challenge in Radiology (Visceral) project to generate testing data. Finally, the DDE algorithm was evaluated with respect to the generalization to different scan parameters and the accuracy of organ dose and effective dose estimates based on an external organ segmentation.

RESULTS

DDE dose distributions were quantified in terms of the mean absolute percentage error (MAPE) and a gamma analysis with respect to the ground truth Monte Carlo simulation. Both measures indicate that DDE generalizes well to different scan parameters and different anatomical regions with a maximum MAPE of 6.3% and a minimum gamma passing rate of 91%. Evaluating the organ dose values for all organs listed in the International Commission on Radiological Protection (ICRP) recommendation, shows an average error of 3.1% and maximum error of 7.2% (bone surface).

CONCLUSIONS

The DDE algorithm provides an efficient approach to determine highly accurate dose distributions. Being able to process a whole-body CT scan in about 1.5 s, it provides a valuable alternative to Monte Carlo simulations on a graphics processing unit (GPU). Here, the main advantage of DDE is that it can be used on top of any existing Monte Carlo code such that real-time performance can be achieved without major adjustments. Thus, DDE opens up new options not only for dosimetry but also for scan and protocol optimization.

Virtual Clinical Trials in 2D and 3D X-ray Breast Imaging and Dosimetry: Comparison of CPU-Based and GPU-Based Monte Carlo Codes

Simple Summary

Virtual clinical trials in X-ray breast imaging may permit substantial reduction of the costs, times, and exposure risk to patient of clinical trials. Monte Carlo simulation techniques are increasingly adopted for VCT in breast imaging and dosimetry studies. This work aims to compare three different platforms for breast VCT studies, to develop real-time virtual DM, DBT and BCT examinations, where the in-silico image acquisition process takes a computational time comparable to that typical of a corresponding real clinical examination.

Abstract

Computational reproductions of medical imaging tests, a form of virtual clinical trials (VCTs), are increasingly being used, particularly in breast imaging research. The accuracy of the computational platform that is used for the imaging and dosimetry simulation processes is a fundamental requirement. Moreover, for practical usage, the imaging simulation computation time should be compatible with the clinical workflow. We compared three different platforms for in-silico X-ray 3D breast imaging: the Agata (University & INFN Napoli) that was based on the Geant4 toolkit and running on a CPU-based server architecture; the XRMC Monte Carlo (University of Cagliari) that was based on the use of variance reduction techniques, running on a CPU hardware; and the Monte Carlo code gCTD (University of Texas Southwestern Medical Center) running on a single GPU platform with CUDA environment. The tests simulated the irradiation of cylindrical objects as well as anthropomorphic breast phantoms and produced 2D and 3D images and 3D maps of absorbed dose. All the codes showed compatible results in terms of simulated dose maps and imaging values within a maximum discrepancy of 3%. The GPU-based code produced a reduction of the computation time up to factor 10⁴, and so permits real-time VCT studies for X-ray breast imaging.

Improving dose calculation accuracy in preclinical radiation experiments using multi-energy element resolved cone beam CT

Cone-beam CT (CBCT) in modern pre-clinical small-animal radiation research platforms provides volumetric images for image guidance and experiment planning purposes. In this work, we implemented multi-energy element-resolved (MEER) CBCT using three scans with different kVps on a SmART platform (Precision X-ray Inc.) We performed comprehensive calibration tasks achieve sufficient accuracy for this quantitative imaging purpose. For geometry calibration, we scanned a ball bearing phantom and used an analytical method together with an optimization approach to derive gantry-angle specific geometry parameters. Intensity calibration and correction included the corrections for detector lag, glare, and beam hardening. The corrected CBCT projection images acquired at 30, 40 and 60 kVp in multiple scans were used to reconstruct CBCT images using the Feldkamp-Davis-Kress reconstruction algorithm. After that, an optimization problem was solved to determine images of relative electron density (rED) and elemental composition (EC) that are needed for Monte Carlo-based radiation dose calculation. We demonstrated effectiveness of our CBCT calibration steps by showing improvements in image quality and successful material decomposition in cases with a small animal CT calibration phantom and a plastinated mouse phantom. It was found that artifacts induced by geometry inaccuracy, detector lag, glare and beam hardening were visually reduced. CT number mean errors were reduced from 19\% to 5\%. In the CT calibration phantom case, median errors in H, O, and Ca fractions for all the inserts were below 1\%, 2\%, and 4\% respectively, and median error in rED was less than 5\%. Compared to standard approach deriving material type and rED via CT number conversion, our approach improved Monte Carlo simulation-based dose calculation accuracy in bone regions. Mean dose error was reduced from 47.5\% to 10.9\%.

Quasi-Monte Carlo method for calculating X-ray scatter in CT

In this paper we transform the trajectories of X-ray as it interacts with a phantom into a high-dimensional integration problem and give the integral formula for the probability of photons emitted from the X-ray source through the phantom to reach the detector. We propose a superior algorithm called gQMCFRD, which combines GPU-based quasi-Monte Carlo (gQMC) method with forced random detection (FRD) technique to simulate this integral. QMC simulation is deterministic versions of Monte Carlo (MC) simulation, which uses deterministic low discrepancy points (such as Sobol' points) instead of the random points. By using the QMC and FRD technique, the gQMCFRD greatly increases the simulation convergence rate and efficiency. We benchmark gQMCFRD, GPU based MC tool (gMCDRR), which performs conventional simulations, a GPU-based Metropolis MC tool (gMMC), which uses the Metropolis-Hasting algorithm to sample the entire photon path from the X-ray source to the detector and gMCFRD, that uses random points for sampling against PENELOPE subroutines: MC-GPU. The results are in excellent agreement and the Efficiency Improvement Factor range 27 ∼ 37 (or 1.09 ∼ 1.16, or 0.12 ∼ 0.15, or 3.62 ∼ 3.70) by gQMCFRD (or gMCDRR, or gMMC, or gMCFRD) with comparison to MC-GPU in all cases. It shows that gQMCFRD is more effective in these cases.

GPU-accelerated Monte Carlo simulation of MV-CBCT

Monte Carlo simulation (MCS) is one of the most accurate computation methods for dose calculation and image formation in radiation therapy. However, the high computational complexity and long execution time of MCS limits its broad use. In this paper, we present a novel strategy to accelerate MCS using a graphic processing unit (GPU), and we demonstrate the application in mega-voltage (MV) cone-beam computed tomography (CBCT) simulation. A new framework that generates a series of MV projections from a single simulation run is designed specifically for MV-CBCT acquisition. A Geant4-based GPU code for photon simulation is incorporated into the framework for the simulation of photon transport through a phantom volume. The FastEPID method, which accelerates the simulation of MV images, is modified and integrated into the framework. The proposed GPU-based simulation strategy was tested for its accuracy and efficiency in a Catphan 604 phantom and an anthropomorphic pelvis phantom with beam energies at 2.5 MV, 6 MV, and 6 MV FFF. In all cases, the proposed GPU-based simulation demonstrated great simulation accuracy and excellent agreement with measurement and CPU-based simulation in terms of reconstructed image qualities. The MV-CBCT simulation was accelerated by factors of roughly 900-2300 using an NVIDIA Tesla V100 GPU card against a 2.5 GHz AMD Opteron™ Processor 6380.

Dose rate determination for preclinical total body irradiation

The accuracy of delivered radiation dose and the reproducibility of employed radiotherapy methods are key factors for preclinical radiobiology applications and research studies. In this work, ionization chamber (IC) measurements and Monte Carlo (MC) simulations were used to accurately determine the dose rate for total body irradiation (TBI), a classic radiobiologic and immunologic experimental method. Several phantom configurations, including large solid water slab, small water box and rodentomorphic mouse and rat phantoms were simulated and measured for TBI setup utilizing a preclinical irradiator XRad320. The irradiator calibration and the phantom measurements were performed using an ADCL calibrated IC N31010 following the AAPM TG-61 protocol. The MC simulations were carried out using Geant4/GATE to compute absorbed dose distributions for all phantom configurations. All simulated and measured geometries had favorable agreement. On average, the relative dose rate difference was 2.3%. However, the study indicated large dose rate deviations, if calibration conditions are assumed for a given experimental setup as commonly done for a quick determination of irradiation times utilizing lookup tables and hand calculations. In a TBI setting, the reference calibration geometry at an extended source-to-surface distance and a large reference field size is likely to overestimate true photon scatter. Consequently, the measured and hand calculated dose rates, for TBI geometries in this study, had large discrepancies: 16% for a large solid water slab, 27% for a small water box, and 31%, 36%, and 30% for mouse phantom, rat phantom, and mouse phantom in a pie cage, respectively. Small changes in TBI experimental setup could result in large dose rate variations. MC simulations and the corresponding measurements specific to a designed experimental setup are vital for accurate preclinical dosimetry and reproducibility of radiobiological findings. This study supports the well-recognized need for physics consultation for all radiobiological investigations.

A general-purpose Monte Carlo particle transport code based on inverse transform sampling for radiotherapy dose calculation

The Monte Carlo (MC) method is widely used to solve various problems in radiotherapy. There has been an impetus to accelerate MC simulation on GPUs whereas thread divergence remains a major issue for MC codes based on acceptance-rejection sampling. Inverse transform sampling has the potential to eliminate thread divergence but it is only implemented for photon transport. Here, we report a MC package Particle Transport in Media (PTM) to demonstrate the implementation of coupled photon-electron transport simulation using inverse transform sampling. Rayleigh scattering, Compton scattering, photo-electric effect and pair production are considered in an analogous manner for photon transport. Electron transport is simulated in a class II condensed history scheme, i.e., catastrophic inelastic scattering and Bremsstrahlung events are simulated explicitly while subthreshold interactions are subject to grouping. A random-hinge electron step correction algorithm and a modified PRESTA boundary crossing algorithm are employed to improve simulation accuracy. Benchmark studies against both EGSnrc simulations and experimental measurements are performed for various beams, phantoms and geometries. Gamma indices of the dose distributions are better than 99.6% for all the tested scenarios under the 2%/2 mm criteria. These results demonstrate the successful implementation of inverse transform sampling in coupled photon-electron transport simulation.

A single neural network for cone-beam computed tomography-based radiotherapy of head-and-neck, lung and breast cancer

•

A deep learning network facilitated dose calculation from CBCT.

•

A single network achieved CBCT-based dose calculation generating synthetic CT for head-and-neck, lung, and breast cancer patients with similar performance to a network specifically trained for each anatomical site.

•

Generation of synthetic-CT can be achieved within 10 s, facilitating online adaptive radiotherapy scenarios.

Background and purpose Adaptive radiotherapy based on cone-beam computed tomography (CBCT) requires high CT number accuracy to ensure accurate dose calculations. Recently, deep learning has been proposed for fast CBCT artefact corrections on single anatomical sites. This study investigated the feasibility of applying a single convolutional network to facilitate dose calculation based on CBCT for head-and-neck, lung and breast cancer patients.

Materials and Methods Ninety-nine patients diagnosed with head-and-neck, lung or breast cancer undergoing radiotherapy with CBCT-based position verification were included in this study. The CBCTs were registered to planning CT according to clinical procedures. Three cycle-consistent generative adversarial networks (cycle-GANs) were trained in an unpaired manner on 15 patients per anatomical site generating synthetic-CTs (sCTs). Another network was trained with all the anatomical sites together. Performances of all four networks were compared and evaluated for image similarity against rescan CT (rCT). Clinical plans were recalculated on rCT and sCT and analysed through voxel-based dose differences and γ-analysis.

Results A sCT was generated in 10 s. Image similarity was comparable between models trained on different anatomical sites and a single model for all sites. Mean dose differences <0.5% were obtained in high-dose regions. Mean gamma (3%, 3 mm) pass-rates >95% were achieved for all sites.

Conclusion Cycle-GAN reduced CBCT artefacts and increased similarity to CT, enabling sCT-based dose calculations. A single network achieved CBCT-based dose calculation generating synthetic CT for head-and-neck, lung, and breast cancer patients with similar performance to a network specifically trained for each anatomical site.

A new open-source GPU-based microscopic Monte Carlo simulation tool for the calculations of DNA damages caused by ionizing radiation --- Part I: Core algorithm and validation

PURPOSE

Monte Carlo (MC) simulation of radiation interactions with water medium at physical, physicochemical, and chemical stages, as well as the computation of biologically relevant quantities such as DNA damages, are of critical importance for the understanding of microscopic basis of radiation effects. Due to the large problem size and many-body simulation problem in the chemical stage, existing CPU-based computational packages encounter the problem of low computational efficiency. This paper reports our development on a GPU-based microscopic Monte Carlo simulation tool gMicroMC using advanced GPU-acceleration techniques.

METHODS

gMicroMC simulated electron transport in the physical stage using an interaction-by-interaction scheme to calculate the initial events generating radicals in water. After the physicochemical stage, initial positions of all radicals were determined. Simulation of radicals' diffusion and reactions in the chemical stage was achieved using a step-by-step model using GPU-accelerated parallelization together with a GPU-enabled box-sorting algorithm to reduce the computations of searching for interaction pairs and therefore improve efficiency. A multi-scale DNA model of the whole lymphocyte cell nucleus containing ~6.2 Gbp of DNA was built.

RESULTS

Accuracy of physical stage simulation was demonstrated by computing stopping power and track length. The results agreed with published data and the data produced by GEANT4-DNA (version 10.3.3) simulations with 10 -20% difference in most cases. Difference of yield values of major radiolytic species from GEANT4-DNA results was within 10%. We computed DNA damages caused by monoenergetic 662 keV photons, approximately representing ¹³⁷ Cs decay. Single-strand break (SSB) and double-strand break (DSB) yields were 196 ± 8 SSB/Gy/Gbp and 7.3 ± 0.7 DSB/Gy/Gbp, respectively, which agreed with the result of 188 SSB/Gy/Gbp and 8.4 DSB/Gy/Gbp computed by Hsiao et al. Compared to computation using a single CPU, gMicroMC achieved a speedup factor of ~540x using an NVidia TITAN Xp GPU card.

CONCLUSIONS

The achieved accuracy and efficiency demonstrated that gMicroMC can facilitate research on microscopic radiation transport simulation and DNA damage calculation. gMicroMC is an open-source package available to the research community.

Technical advances in x-ray microbeam radiation therapy

In the last 25 years microbeam radiation therapy (MRT) has emerged as a promising alternative to conventional radiation therapy at large, third generation synchrotrons. In MRT, a multi-slit collimator modulates a kilovoltage x-ray beam on a micrometer scale, creating peak dose areas with unconventionally high doses of several hundred Grays separated by low dose valley regions, where the dose remains well below the tissue tolerance level. Pre-clinical evidence demonstrates that such beam geometries lead to substantially reduced damage to normal tissue at equal tumour control rates and hence drastically increase the therapeutic window. Although the mechanisms behind MRT are still to be elucidated, previous studies indicate that immune response, tumour microenvironment, and the microvasculature may play a crucial role. Beyond tumour therapy, MRT has also been suggested as a microsurgical tool in neurological disorders and as a primer for drug delivery. The physical properties of MRT demand innovative medical physics and engineering solutions for safe treatment delivery. This article reviews technical developments in MRT and discusses existing solutions for dosimetric validation, reliable treatment planning and safety. Instrumentation at synchrotron facilities, including beam production, collimators and patient positioning systems, is also discussed. Specific solutions reviewed in this article include: dosimetry techniques that can cope with high spatial resolution, low photon energies and extremely high dose rates of up to 15 000 Gy s-1, dose calculation algorithms-apart from pure Monte Carlo Simulations-to overcome the challenge of small voxel sizes and a wide dynamic dose-range, and the use of dose-enhancing nanoparticles to combat the limited penetrability of a kilovoltage energy spectrum. Finally, concepts for alternative compact microbeam sources are presented, such as inverse Compton scattering set-ups and carbon nanotube x-ray tubes, that may facilitate the transfer of MRT into a hospital-based clinical environment. Intensive research in recent years has resulted in practical solutions to most of the technical challenges in MRT. Treatment planning, dosimetry and patient safety systems at synchrotrons have matured to a point that first veterinary and clinical studies in MRT are within reach. Should these studies confirm the promising results of pre-clinical studies, the authors are confident that MRT will become an effective new radiotherapy option for certain patients.

Innovations in Computer Technologies Have Impacted Radiation Dosimetry Through Anatomically Realistic Phantoms and Fast Monte Carlo Simulations

Radiological physics principles have not changed in the past 60 y when computer technologies advanced exponentially. The research field of anatomical modeling for the purpose of radiation dose calculations has experienced an explosion in activity in the past two decades. Such an exciting advancement is due to the feasibility of creating three-dimensional geometric details of the human anatomy from tomographic imaging and of performing Monte Carlo radiation transport simulations on increasingly fast and cheap personal computers. The advent of a new type of high-performance computing hardware in recent years-graphics processing units-has made it feasible to carry out time-consuming Monte Carlo calculations at near real-time speeds. This paper introduces the history of three generations of computational human phantoms (the stylized medical internal radiation dosimetry-type phantoms, the voxelized tomographic phantoms, and the boundary representation deformable phantoms) and new development of the graphics processing unit-based Monte Carlo radiation dose calculations. Examples are given for research projects performed by my students in applying computational phantoms and a new Monte Carlo code, ARCHER, to problems in radiation protection, imaging, and radiotherapy. Finally, the paper discusses challenges and future opportunities for research.

Metropolis Monte Carlo simulation scheme for fast scattered X-ray photon calculation in CT

Monte Carlo (MC) method is commonly considered as the most accurate approach for particle transport simulation because of its capability to precisely model physics interactions and simulation geometry. Conventionally, MC simulation is performed in a particle-by-particle fashion. In certain problems such as computing scattered X-ray photon signal at a detector of CT, the conventional simulation scheme suffers from low efficiency mainly due to the fact that abundant photons are simulated but do not reach the detector. The computational resources spent on those photons are therefore wasted. To solve this problem, this study develops a novel GPU-based Metropolis MC (gMMC) with a novel path-by-path simulation scheme and demonstrates its effectiveness in an example problem of scattered X-ray photon calculation in CT. In contrast to the conventional MC approach, gMMC samples an entire photon path extending from the X-ray source to the detector using Metropolis-Hasting algorithm. The path-by-path simulation scheme ensures contribution of every sampled event to the signal of interest, improving overall efficiency. We benchmark gMMC against an in-house developed GPU-based MC tool, gMCDRR, which performs simulations in the conventional particle-by-particle fashion. gMMC reaches speed up factors of 37~48 times in simple phantom cases and 20-34 times in real patient cases. The results calculated by gMCDRR and gMMC agree well with average differences < 3%.

A fast, linear Boltzmann transport equation solver for computed tomography dose calculation (Acuros CTD)

PURPOSE

To improve dose reporting of CT scans, patient-specific organ doses are highly desired. However, estimating the dose distribution in a fast and accurate manner remains challenging, despite advances in Monte Carlo methods. In this work, we present an alternative method that deterministically solves the linear Boltzmann transport equation (LBTE), which governs the behavior of x-ray photon transport through an object.

METHODS

Our deterministic solver for CT dose (Acuros CTD) is based on the same approach used to estimate scatter in projection images of a CT scan (Acuros CTS). A deterministic method is used to compute photon fluence within the object, which is then converted to deposited energy by multiplying by known, material-specific conversion factors. To benchmark Acuros CTD, we used the AAPM Task Group 195 test for CT dose, which models an axial, fan beam scan (10 mm thick beam) and calculates energy deposited in each organ of an anthropomorphic phantom. We also validated our own Monte Carlo implementation of Geant4 to use as a reference to compare Acuros against for other common geometries like an axial, cone beam scan (160 mm thick beam) and a helical scan (40 mm thick beam with table motion for a pitch of 1).

RESULTS

For the fan beam scan, Acuros CTD accurately estimated organ dose, with a maximum error of 2.7% and RMSE of 1.4% when excluding organs with <0.1% of the total energy deposited. The cone beam and helical scans yielded similar levels of accuracy compared to Geant4. Increasing the number of source positions beyond 18 or decreasing the voxel size below 5 × 5 × 5 mm³ provided marginal improvement to the accuracy for the cone beam scan but came at the expense of increased run time. Across the different scan geometries, run time of Acuros CTD ranged from 8 to 23 s.

CONCLUSIONS

In this digital phantom study, a deterministic LBTE solver was capable of fast and accurate organ dose estimates.

Hybrid dose calculation: a dose calculation algorithm for microbeam radiation therapy

Microbeam radiation therapy (MRT) is still a preclinical approach in radiation oncology that uses planar micrometre wide beamlets with extremely high peak doses, separated by a few hundred micrometre wide low dose regions. Abundant preclinical evidence demonstrates that MRT spares normal tissue more effectively than conventional radiation therapy, at equivalent tumour control. In order to launch first clinical trials, accurate and efficient dose calculation methods are an inevitable prerequisite. In this work a hybrid dose calculation approach is presented that is based on a combination of Monte Carlo and kernel based dose calculation. In various examples the performance of the algorithm is compared to purely Monte Carlo and purely kernel based dose calculations. The accuracy of the developed algorithm is comparable to conventional pure Monte Carlo calculations. In particular for inhomogeneous materials the hybrid dose calculation algorithm out-performs purely convolution based dose calculation approaches. It is demonstrated that the hybrid algorithm can efficiently calculate even complicated pencil beam and cross firing beam geometries. The required calculation times are substantially lower than for pure Monte Carlo calculations.

A point kernel algorithm for microbeam radiation therapy

Microbeam radiation therapy (MRT) is a treatment approach in radiation therapy where the treatment field is spatially fractionated into arrays of a few tens of micrometre wide planar beams of unusually high peak doses separated by low dose regions of several hundred micrometre width. In preclinical studies, this treatment approach has proven to spare normal tissue more effectively than conventional radiation therapy, while being equally efficient in tumour control. So far dose calculations in MRT, a prerequisite for future clinical applications are based on Monte Carlo simulations. However, they are computationally expensive, since scoring volumes have to be small. In this article a kernel based dose calculation algorithm is presented that splits the calculation into photon and electron mediated energy transport, and performs the calculation of peak and valley doses in typical MRT treatment fields within a few minutes. Kernels are analytically calculated depending on the energy spectrum and material composition. In various homogeneous materials peak, valley doses and microbeam profiles are calculated and compared to Monte Carlo simulations. For a microbeam exposure of an anthropomorphic head phantom calculated dose values are compared to measurements and Monte Carlo calculations. Except for regions close to material interfaces calculated peak dose values match Monte Carlo results within 4% and valley dose values within 8% deviation. No significant differences are observed between profiles calculated by the kernel algorithm and Monte Carlo simulations. Measurements in the head phantom agree within 4% in the peak and within 10% in the valley region. The presented algorithm is attached to the treatment planning platform VIRTUOS. It was and is used for dose calculations in preclinical and pet-clinical trials at the biomedical beamline ID17 of the European synchrotron radiation facility in Grenoble, France.

GATE Monte Carlo simulation of dose distribution using MapReduce in a cloud computing environment

The GATE Monte Carlo simulation platform has good application prospects of treatment planning and quality assurance. However, accurate dose calculation using GATE is time consuming. The purpose of this study is to implement a novel cloud computing method for accurate GATE Monte Carlo simulation of dose distribution using MapReduce. An Amazon Machine Image installed with Hadoop and GATE is created to set up Hadoop clusters on Amazon Elastic Compute Cloud (EC2). Macros, the input files for GATE, are split into a number of self-contained sub-macros. Through Hadoop Streaming, the sub-macros are executed by GATE in Map tasks and the sub-results are aggregated into final outputs in Reduce tasks. As an evaluation, GATE simulations were performed in a cubical water phantom for X-ray photons of 6 and 18 MeV. The parallel simulation on the cloud computing platform is as accurate as the single-threaded simulation on a local server and the simulation correctness is not affected by the failure of some worker nodes. The cloud-based simulation time is approximately inversely proportional to the number of worker nodes. For the simulation of 10 million photons on a cluster with 64 worker nodes, time decreases of 41× and 32× were achieved compared to the single worker node case and the single-threaded case, respectively. The test of Hadoop's fault tolerance showed that the simulation correctness was not affected by the failure of some worker nodes. The results verify that the proposed method provides a feasible cloud computing solution for GATE.

Monte Carlo dose calculations for high-dose-rate brachytherapy using GPU-accelerated processing

PURPOSE

Current clinical brachytherapy dose calculations are typically based on the Association of American Physicists in Medicine Task Group report 43 (TG-43) guidelines, which approximate patient geometry as an infinitely large water phantom. This ignores patient and applicator geometries and heterogeneities, causing dosimetric errors. Although Monte Carlo (MC) dose calculation is commonly recognized as the most accurate method, its associated long computational time is a major bottleneck for routine clinical applications. This article presents our recent developments of a fast MC dose calculation package for high-dose-rate (HDR) brachytherapy, gBMC, built on a graphics processing unit (GPU) platform.

METHODS AND MATERIALS

gBMC-simulated photon transport in voxelized geometry with physics in (192)Ir HDR brachytherapy energy range considered. A phase-space file was used as a source model. GPU-based parallel computation was used to simultaneously transport multiple photons, one on a GPU thread. We validated gBMC by comparing the dose calculation results in water with that computed TG-43. We also studied heterogeneous phantom cases and a patient case and compared gBMC results with Acuros BV results.

RESULTS

Radial dose function in water calculated by gBMC showed <0.6% relative difference from that of the TG-43 data. Difference in anisotropy function was <1%. In two heterogeneous slab phantoms and one shielded cylinder applicator case, average dose discrepancy between gBMC and Acuros BV was <0.87%. For a tandem and ovoid patient case, good agreement between gBMC and Acruos BV results was observed in both isodose lines and dose-volume histograms. In terms of the efficiency, it took ∼47.5 seconds for gBMC to reach 0.15% statistical uncertainty within the 5% isodose line for the patient case.

CONCLUSIONS

The accuracy and efficiency of a new GPU-based MC dose calculation package, gBMC, for HDR brachytherapy make it attractive for clinical applications.

Comprehensive Monte Carlo study of patient doses from cone-beam CT imaging in radiotherapy

Accurate knowledge of ionizing radiation dose from cone-beam CT (CBCT) imaging in radiotherapy is important to allow concomitant risks to be estimated and for justification of imaging exposures. This study uses a Monte Carlo CBCT model to calculate imaging dose for a wide range of imaging protocols for male and female patients. The Elekta XVI CBCT system was modeled using GATE and simulated doses were validated against measurements in a water tank and thorax phantom. Imaging dose was simulated in the male and female ICRP voxel phantoms for a variety of anatomical sites and imager settings (different collimators, filters, full and partial rotation). The resulting dose distributions were used to calculate effective doses for each scan protocol. The Monte Carlo simulated doses agree with validation measurements within 5% and 10% for water tank and thorax phantom respectively. Effective dose for head CBCT scans was generally lower for scans centred on the pituitary than the larynx (0.03 mSv versus 0.06 mSv for male ICRP phantom). Pelvis CBCT scan effective dose was higher for the female than male phantom (5.11 mSv versus 2.80 mSv for M15 collimator scan), principally due to the higher dose received by gonads for the female scan. Medium field of view thorax scan effective doses ranged from 1.38-3.19 mSv depending on scan length and phantom sex. Effective dose for half rotation thorax scans with offset isocentre varied by almost a factor of three depending on laterality of the isocentre, patient sex and imaged field length. The CBCT imaging doses simulated here reveal large variations in dose depending on imaging isocentre location, patient sex and partial rotation angles. This information may be used to estimate risks from CBCT and to optimize CBCT imaging protocols.

Moving GPU-OpenCL-based Monte Carlo dose calculation toward clinical use: Automatic beam commissioning and source sampling for treatment plan dose calculation

We have previously developed a GPU‐based Monte Carlo (MC) dose engine on the OpenCL platform, named goMC, with a built‐in analytical linear accelerator (linac) beam model. In this paper, we report our recent improvement on goMC to move it toward clinical use. First, we have adapted a previously developed automatic beam commissioning approach to our beam model. The commissioning was conducted through an optimization process, minimizing the discrepancies between calculated dose and measurement. We successfully commissioned six beam models built for Varian TrueBeam linac photon beams, including four beams of different energies (6 MV, 10 MV, 15 MV, and 18 MV) and two flattening‐filter‐free (FFF) beams of 6 MV and 10 MV. Second, to facilitate the use of goMC for treatment plan dose calculations, we have developed an efficient source particle sampling strategy. It uses the pre‐generated fluence maps (FMs) to bias the sampling of the control point for source particles already sampled from our beam model. It could effectively reduce the number of source particles required to reach a statistical uncertainty level in the calculated dose, as compared to the conventional FM weighting method. For a head‐and‐neck patient treated with volumetric modulated arc therapy (VMAT), a reduction factor of ~2.8 was achieved, accelerating dose calculation from 150.9 s to 51.5 s. The overall accuracy of goMC was investigated on a VMAT prostate patient case treated with 10 MV FFF beam. 3D gamma index test was conducted to evaluate the discrepancy between our calculated dose and the dose calculated in Varian Eclipse treatment planning system. The passing rate was 99.82% for 2%/2 mm criterion and 95.71% for 1%/1 mm criterion. Our studies have demonstrated the effectiveness and feasibility of our auto‐commissioning approach and new source sampling strategy for fast and accurate MC dose calculations for treatment plans.

Prospective Randomized Trial for Image-Guided Biopsy Using Cone-Beam CT Navigation Compared with Conventional CT

PURPOSE

To compare cone-beam computed tomography (CT) navigation vs conventional CT image guidance during biopsies.

MATERIALS AND METHODS

Patients scheduled for image-guided biopsies were prospectively and randomly assigned to conventional CT guidance vs cone-beam CT navigation. Radiation dose, accuracy of final needle position, rate of histopathologic diagnosis, and number of needle repositions to reach the target (defined as pullback to adjust position) were compared.

RESULTS

A total of 58 patients (mean age, 57 y; 62.1% men) were randomized: 29 patients underwent 33 biopsies with CT guidance and 29 patients with 33 lesions underwent biopsy with cone-beam CT navigation. The average body mass index (BMI) was similar between groups, at 28.8 kg/m(2) ± 6.55 (P = .18). There was no difference between groups in terms of patient and lesion characteristics (eg, size, depth). The average lesion size was 29.1 ± 12.7mm for CT group vs 32.1mm ±16.8mm for cone-beam CT group (P < 0.59). Location of lesions was equally divided between the 2 groups, 20 lung lesions, 18 renal lesions and 20 other abdominal lesions. Mean number of needle repositions in the cone-beam CT group was 0.3 ± 0.5, compared with 1.9 ± 2.3 with conventional CT (P < .001). The average skin entry dose was 29% lower with cone-beam CT than with conventional CT (P < .04 accounting for BMI). The average estimated effective dose for the planning scan from phantom data was 49% lower with cone-beam CT vs conventional CT (P = .018). Accuracy, defined as the difference between planned and final needle positions, was 4.9 mm ± 4.1 for the cone-beam CT group, compared with 12.2 mm ± 8.1 for conventional CT (P < .001). Histopathologic diagnosis rates were similar between groups, at 90.9% for conventional CT and 93.9% for cone-beam CT (P = .67).

CONCLUSIONS

Cone-beam CT navigation for biopsies improved targeting accuracy with fewer needle repositions, lower skin entry dose, and lower effective dose for planning scan, and a comparable histopathologic diagnosis rate.

GPU-based RFA simulation for minimally invasive cancer treatment of liver tumours

PURPOSE

Radiofrequency ablation (RFA) is one of the most popular and well-standardized minimally invasive cancer treatments (MICT) for liver tumours, employed where surgical resection has been contraindicated. Less-experienced interventional radiologists (IRs) require an appropriate planning tool for the treatment to help avoid incomplete treatment and so reduce the tumour recurrence risk. Although a few tools are available to predict the ablation lesion geometry, the process is computationally expensive. Also, in our implementation, a few patient-specific parameters are used to improve the accuracy of the lesion prediction.

METHODS

Advanced heterogeneous computing using personal computers, incorporating the graphics processing unit (GPU) and the central processing unit (CPU), is proposed to predict the ablation lesion geometry. The most recent GPU technology is used to accelerate the finite element approximation of Penne's bioheat equation and a three state cell model. Patient-specific input parameters are used in the bioheat model to improve accuracy of the predicted lesion.

RESULTS

A fast GPU-based RFA solver is developed to predict the lesion by doing most of the computational tasks in the GPU, while reserving the CPU for concurrent tasks such as lesion extraction based on the heat deposition at each finite element node. The solver takes less than 3 min for a treatment duration of 26 min. When the model receives patient-specific input parameters, the deviation between real and predicted lesion is below 3 mm.

CONCLUSION

A multi-centre retrospective study indicates that the fast RFA solver is capable of providing the IR with the predicted lesion in the short time period before the intervention begins when the patient has been clinically prepared for the treatment.

A level set method for cupping artifact correction in cone-beam CT

PURPOSE

To reduce cupping artifacts and improve the contrast-to-noise ratio in cone-beam computed tomography (CBCT).

METHODS

A level set method is proposed to reduce cupping artifacts in the reconstructed image of CBCT. The authors derive a local intensity clustering property of the CBCT image and define a local clustering criterion function of the image intensities in a neighborhood of each point. This criterion function defines an energy in terms of the level set functions, which represent a segmentation result and the cupping artifacts. The cupping artifacts are estimated as a result of minimizing this energy.

RESULTS

The cupping artifacts in CBCT are reduced by an average of 90%. The results indicate that the level set-based algorithm is practical and effective for reducing the cupping artifacts and preserving the quality of the reconstructed image.

CONCLUSIONS

The proposed method focuses on the reconstructed image without requiring any additional physical equipment, is easily implemented, and provides cupping correction through a single-scan acquisition. The experimental results demonstrate that the proposed method successfully reduces the cupping artifacts.

A new scheme for real-time high-contrast imaging in lung cancer radiotherapy: a proof-of-concept study

Visualization of anatomy in real time is of critical importance for motion management in lung cancer radiotherapy. To achieve real-time, and high-contrast in-treatment imaging, we propose a novel scheme based on the measurement of Compton scatter photons. In our method, a slit x-ray beam along the superior-inferior direction is directed to the patient, (intersecting the lung region at a 2D plane) containing most of the tumor motion trajectory. X-ray photons are scattered off this plane primarily due to the Compton interaction. An imager with a pinhole or a slat collimator is placed at one side of the plane to capture the scattered photons. The resulting image, after correcting for incoming fluence inhomogeneity, x-ray attenuation, scatter angle variation, and outgoing beam geometry, represents the linear attenuation coefficient of Compton scattering. This allows the visualization of the anatomy on this plane. We performed Monte Carlo simulation studies both on a phantom and a patient for proof-of-principle purposes. In the phantom case, a small tumor-like structure could be clearly visualized. The contrast-resolution calculated using tumor/lung as foreground/background for kV fluoroscopy, cone beam computed tomography (CBCT), and scattering image were 0.037, 0.70, and 0.54, respectively. In the patient case, tumor motion could be clearly observed in the scatter images. Imaging dose to the voxels directly exposed by the slit beam was ~0.4 times of that under a single CBCT projection. These studies demonstrated the potential feasibility of the proposed imaging scheme to capture the instantaneous anatomy of a patient on a 2D plane with a high image contrast. Clear visualization of the tumor motion may facilitate marker-less tumor tracking.

GMctdospp: Description and validation of a CT dose calculation system

PURPOSE

To develop a Monte Carlo (MC)-based computed tomography (CT) dose estimation method with a graphical user interface with options to define almost arbitrary simulation scenarios, to make calculations sufficiently fast for comfortable handling, and to make the software free of charge for general availability to the scientific community.

METHODS

A framework called GMctdospp was developed to calculate phantom and patient doses with the MC method based on the EGSnrc system. A CT scanner was modeled for testing and was adapted to half-value layer, beam-shaping filter, z-profile, and tube-current modulation (TCM). To validate the implemented variance reduction techniques, depth-dose and cross-profile calculations of a static beam were compared against DOSXYZnrc/EGSnrc. Measurements for beam energies of 80 and 120 kVp at several positions of a CT dose-index (CTDI) standard phantom were compared against calculations of the created CT model. Finally, the efficiency of the adapted code was benchmarked against EGSnrc defaults.

RESULTS

The CT scanner could be modeled accurately. The developed TCM scheme was confirmed by the dose measurement. A comparison of calculations to DOSXYZnrc showed no systematic differences. Measurements in a CTDI phantom could be reproduced within 2% average, with a maximal difference of about 6%. Efficiency improvements of about six orders of magnitude were observed for larger organ structures of a chest-examination protocol in a voxelized phantom. In these cases, simulations took 25 s to achieve a statistical uncertainty of ∼0.5%.

CONCLUSIONS

A fast dose-calculation system for phantoms and patients in a CT examination was developed, successfully validated, and benchmarked. Influences of scan protocols, protection method, and other issues can be easily examined with the developed framework.

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

PURPOSE

Using the graphical processing units (GPU) hardware technology, an extremely fast Monte Carlo (MC) code ARCHERRT is developed for radiation dose calculations in radiation therapy. This paper describes the detailed software development and testing for three clinical TomoTherapy® cases: the prostate, lung, and head & neck.

METHODS

To obtain clinically relevant dose distributions, phase space files (PSFs) created from optimized radiation therapy treatment plan fluence maps were used as the input to ARCHERRT. Patient-specific phantoms were constructed from patient CT images. Batch simulations were employed to facilitate the time-consuming task of loading large PSFs, and to improve the estimation of statistical uncertainty. Furthermore, two different Woodcock tracking algorithms were implemented and their relative performance was compared. The dose curves of an Elekta accelerator PSF incident on a homogeneous water phantom were benchmarked against DOSXYZnrc. For each of the treatment cases, dose volume histograms and isodose maps were produced from ARCHERRT and the general-purpose code, GEANT4. The gamma index analysis was performed to evaluate the similarity of voxel doses obtained from these two codes. The hardware accelerators used in this study are one NVIDIA K20 GPU, one NVIDIA K40 GPU, and six NVIDIA M2090 GPUs. In addition, to make a fairer comparison of the CPU and GPU performance, a multithreaded CPU code was developed using OpenMP and tested on an Intel E5-2620 CPU.

RESULTS

For the water phantom, the depth dose curve and dose profiles from ARCHERRT agree well with DOSXYZnrc. For clinical cases, results from ARCHERRT are compared with those from GEANT4 and good agreement is observed. Gamma index test is performed for voxels whose dose is greater than 10% of maximum dose. For 2%/2mm criteria, the passing rates for the prostate, lung case, and head & neck cases are 99.7%, 98.5%, and 97.2%, respectively. Due to specific architecture of GPU, modified Woodcock tracking algorithm performed inferior to the original one. ARCHERRT achieves a fast speed for PSF-based dose calculations. With a single M2090 card, the simulations cost about 60, 50, 80 s for three cases, respectively, with the 1% statistical error in the PTV. Using the latest K40 card, the simulations are 1.7-1.8 times faster. More impressively, six M2090 cards could finish the simulations in 8.9-13.4 s. For comparison, the same simulations on Intel E5-2620 (12 hyperthreading) cost about 500-800 s.

CONCLUSIONS

ARCHERRT was developed successfully to perform fast and accurate MC dose calculation for radiotherapy using PSFs and patient CT phantoms.

Reconstructing cone-beam CT with spatially varying qualities for adaptive radiotherapy: a proof-of-principle study

With the aim of maximally reducing imaging dose while meeting requirements for adaptive radiation therapy (ART), we propose in this paper a new cone beam CT (CBCT) acquisition and reconstruction method that delivers images with a low noise level inside a region of interest (ROI) and a relatively high noise level outside the ROI. The acquired projection images include two groups: densely sampled projections at a low exposure with a large field of view (FOV) and sparsely sampled projections at a high exposure with a small FOV corresponding to the ROI. A new algorithm combining the conventional filtered back-projection algorithm and the tight-frame iterative reconstruction algorithm is also designed to reconstruct the CBCT based on these projection data. We have validated our method on a simulated head-and-neck (HN) patient case, a semi-real experiment conducted on a HN cancer patient under a full-fan scan mode, as well as a Catphan phantom under a half-fan scan mode. Relative root-mean-square errors (RRMSEs) of less than 3% for the entire image and ~1% within the ROI compared to the ground truth have been observed. These numbers demonstrate the ability of our proposed method to reconstruct high-quality images inside the ROI. As for the part outside ROI, although the images are relatively noisy, it can still provide sufficient information for radiation dose calculations in ART. Dose distributions calculated on our CBCT image and on a standard CBCT image are in agreement, with a mean relative difference of 0.082% inside the ROI and 0.038% outside the ROI. Compared with the standard clinical CBCT scheme, an imaging dose reduction of approximately 3-6 times inside the ROI was achieved, as well as an 8 times outside the ROI. Regarding computational efficiency, it takes 1-3 min to reconstruct a CBCT image depending on the number of projections used. These results indicate that the proposed method has the potential for application in ART.

Dose calculation and treatment plan optimization including imaging dose from kilovoltage cone beam computed tomography

BACKGROUND

With the increasing use of cone beam computed tomography (CBCT) for patient position verification and radiotherapy treatment adaptation, there is an increasing need to develop techniques that can take into account concomitant dose using a personalized approach.

MATERIAL AND METHODS

A total of 20 patients (10 pelvis and 10 head and neck) who had undergone radiation therapy using intensity modulated radiation therapy (IMRT) were selected and the dose from kV CBCT was retrospectively calculated using a treatment planning system previously commissioned for this purpose. The imaging dose was added to the CT images used for treatment planning and the difference in its addition prior to and after the planning was assessed.

RESULTS

The additional isocenter dose as a result of daily CBCT is in the order of 3-4 cGy for 35-fraction head and neck and 23-47 cGy for 25-fraction pelvis cases using the standard head and neck and pelvis image acquisition protocols. The pelvic dose is especially dependent on patient size and body mass index (BMI), being higher for patients with lower BMI. Due to the low energy of the kV CBCT beam, the maximum energy deposition is at or near the surface with the highest dose being on the patient's left side for the head and neck (∼7 cGy) and on the posterior for the pelvic cases (∼80 cGy). Addition of imaging dose prior to plan optimization resulted in an average reduction of 4% in the plan monitor units and 5% in the number of control points.

CONCLUSION

Dose from daily kV CBCT has been added to patient treatment plans using previously commissioned kV CBCT beams in a treatment planning system. Addition of imaging dose can be included in IMRT treatment plan optimization and would facilitate customization of imaging protocol based on patient anatomy and location of isocenter.

GPU-based high-performance computing for radiation therapy

Recent developments in radiotherapy therapy demand high computation powers to solve challenging problems in a timely fashion in a clinical environment. The graphics processing unit (GPU), as an emerging high-performance computing platform, has been introduced to radiotherapy. It is particularly attractive due to its high computational power, small size, and low cost for facility deployment and maintenance. Over the past few years, GPU-based high-performance computing in radiotherapy has experienced rapid developments. A tremendous amount of study has been conducted, in which large acceleration factors compared with the conventional CPU platform have been observed. In this paper, we will first give a brief introduction to the GPU hardware structure and programming model. We will then review the current applications of GPU in major imaging-related and therapy-related problems encountered in radiotherapy. A comparison of GPU with other platforms will also be presented.

Comprehensive evaluations of cone-beam CT dose in image-guided radiation therapy via GPU-based Monte Carlo simulations

Cone beam CT (CBCT) has been widely used for patient setup in image-guided radiation therapy (IGRT). Radiation dose from CBCT scans has become a clinical concern. The purposes of this study are (1) to commission a graphics processing unit (GPU)-based Monte Carlo (MC) dose calculation package gCTD for Varian On-Board Imaging (OBI) system and test the calculation accuracy, and (2) to quantitatively evaluate CBCT dose from the OBI system in typical IGRT scan protocols. We first conducted dose measurements in a water phantom. X-ray source model parameters used in gCTD are obtained through a commissioning process. gCTD accuracy is demonstrated by comparing calculations with measurements in water and in CTDI phantoms. Twenty-five brain cancer patients are used to study dose in a standard-dose head protocol, and 25 prostate cancer patients are used to study dose in pelvis protocol and pelvis spotlight protocol. Mean dose to each organ is calculated. Mean dose to 2% voxels that have the highest dose is also computed to quantify the maximum dose. It is found that the mean dose value to an organ varies largely among patients. Moreover, dose distribution is highly non-homogeneous inside an organ. The maximum dose is found to be 1-3 times higher than the mean dose depending on the organ, and is up to eight times higher for the entire body due to the very high dose region in bony structures. High computational efficiency has also been observed in our studies, such that MC dose calculation time is less than 5 min for a typical case.