An automation of Monte Carlo workflow for dosimetry study of an Elekta LINAC delivery system in radiotherapy.

Purpose: This study aims to automate the Monte Carlo (MC) workflow utilized for radiotherapy dosimetry, focusing on an Elekta LINAC delivery system. It addresses the challenge of integrating MC simulations into routine clinical practice, making this accurate yet complex method more accessible and efficient for radiotherapy dosimetry. Methods and Materials: We developed a user-friendly software featuring a graphical user interface (GUI) that integrates EGSnrc for MC simulations. The software streamlines the process from retrieving Digital Imaging and Communications in Medicine (DICOM) data to executing dose calculations and comparing dose distributions. To validate our proposed tool, we compared its computed doses for IMRT and VMAT plans from the Pinnacle TPS for an Elekta Versa HD linear accelerator against MC simulation results. This comparison utilized our in-house software and GUI as the tool, covering various treatment sites and prescriptions. Results: The automated MC workflow demonstrated high accuracy in dose calculations and streamlined integration with clinical workflows. The comparison between the MC-simulated and TPS-calculated doses revealed excellent agreement, highlighting the reliability of MC for independent dose verification in complex treatment scenarios. Conclusions: The automated MC workflow developed represents a substantial improvement in the practicality and efficiency of MC simulations in radiotherapy. This advancement not only simplifies the dosimetry process but also ensures high accuracy, establishing it as a valuable tool for routine patient-specific quality assurance and the development of specialized treatment procedures.

MR-linac MLC positioning QA by digitally stitching dual double-exposed films.

PURPOSE: The picket fence (PF) test is highly recommended for multi-leaf collimator (MLC) quality assurance. However, since the electronic portal imaging device (EPID) on the Elekta Unity only covers a small area, it is not feasible to perform the PF test for the entire MLC. Here, we propose a technique for the PF test by stitching two double-exposed films. METHODS: Two EBT3 films were used to encompass the entire MLC, with each one covering one half of the area. Two fields were employed to apply double exposure: a PF pattern consisting of 11 2 mm wide pickets and a 2.84 cm x 22 cm open field. The edges of the open field defined by the diaphragms were used to correct film rotation as well as align them horizontally. The PF pattern was also measured with the EPID where the pickets were used to align the films vertically. Individual leaf positions were detected on the merged film for quantitative analysis. Various MLC positioning errors were introduced to evaluate the technique's sensitivity. RESULTS: The merged films covered 72 leaf pairs properly (four leaf pairs on both sides were outside the treatment couch). With the EPID, the leaf positioning accuracy was -0.02 ± 0.07 mm (maximum: 0.29 mm) and the picket width variation was 0.00 ± 0.03 mm (maximum: 0.11 mm); with the films, the position accuracy and width variation were -0.03 ± 0.13 mm (maximum: 0.80 mm) and 0.00 ± 0.13 mm (maximum: 0.74 mm), respectively. The EPID was able to detect errors of 0.5 mm or above with submillimeter accuracy; the films were only able to detect errors > 1.0 mm. CONCLUSION: We developed a quantitative technique for the PF test on the Elekta Unity. The merged films covered nearly the entire MLC leaf banks. The technique exhibited clinically acceptable accuracy and sensitivity to MLC positioning errors.

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.

An ultra-fast deep-learning-based dose engine for prostate VMAT via knowledge distillation framework with limited patient data.

Objective. Deep-learning (DL)-based dose engines have been developed to alleviate the intrinsic compromise between the calculation accuracy and efficiency of the traditional dose calculation algorithms. However, current DL-based engines typically possess high computational complexity and require powerful computing devices. Therefore, to mitigate their computational burdens and broaden their applicability to a clinical setting where resource-limited devices are available, we proposed a compact dose engine via knowledge distillation (KD) framework that offers an ultra-fast calculation speed with high accuracy for prostate Volumetric Modulated Arc Therapy (VMAT).Approach. The KD framework contains two sub-models: a large pre-trained teacher and a small to-be-trained student. The student receives knowledge transferred from the teacher for better generalization. The trained student serves as the final engine for dose calculation. The model input is patient computed tomography and VMAT dose in water, and the output is DL-calculated patient dose. The ground-truth \dose was computed by the Monte Carlo module of the Monaco treatment planning system. Twenty and ten prostate cases were included for model training and assessment, respectively. The model's performance (teacher/student/student-only) was evaluated by Gamma analysis and inference efficiency.Main results. The dosimetric comparisons (input/DL-calculated/ground-truth doses) suggest that the proposed engine can effectively convert low-accuracy doses in water to high-accuracy patient doses. The Gamma passing rate (2%/2 mm, 10% threshold) between the DL-calculated and ground-truth doses was 98.64 ± 0.62% (teacher), 98.13 ± 0.76% (student), and 96.95 ± 1.02% (student-only). The inference time was 16 milliseconds (teacher) and 11 milliseconds (student/student-only) using a graphics processing unit device, while it was 936 milliseconds (teacher) and 374 milliseconds (student/student-only) using a central processing unit device.Significance. With the KD framework, a compact dose engine can achieve comparable accuracy to that of a larger one. Its compact size reduces the computational burdens and computing device requirements, and thus such an engine can be more clinically applicable.

A polar-coordinate-based pencil beam algorithm for VMAT dose computation with high-resolution gantry angle sampling.

PURPOSE: Most commercially available treatment planning systems (TPSs) approximate the continuous delivery of volumetric modulated arc therapy (VMAT) plans with a series of discretized static beams for treatment planning, which can make VMAT dose computation extremely inefficient. In this study, we developed a polar-coordinate-based pencil beam (PB) algorithm for efficient VMAT dose computation with high-resolution gantry angle sampling that can improve the computational efficiency and reduce the dose discrepancy due to the angular under-sampling effect. METHODS AND MATERIALS: 6 MV 1 × 1 m m 2 $1 \times 1{\rm{\;m}}{{\rm{m}}^2}$ pencil beams were simulated on a uniform cylindrical phantom under an EGSnrc Monte Carlo (MC) environment. The MC-generated PB kernels were collected in the polar coordinate system for each bixel on a 40 × 40 c m 2 $40 \times 40{\rm{\;c}}{{\rm{m}}^2}$ fluence map and subsequently fitted via a series of Gaussians. The fluence was calculated using a detectors' eye view with off-axis and MLC transmission factors corrected. Doses of VMAT arc on the phantom were computed by summing the convolution results between the corresponding PB kernels and fluence for each bixel in the polar coordinate system. The convolution was performed using fast Fourier transform to expedite the computing speed. The calculated doses were converted to the Cartesian coordinate system and compared with the reference dose computed by a collapsed cone convolution (CCC) algorithm of the TPS. A heterogeneous phantom was created to study the heterogeneity corrections using the proposed algorithm. Ten VMAT arcs were included to evaluate the algorithm performance. Gamma analysis and computation complexity theory were used to measure the dosimetric accuracy and computational efficiency, respectively. RESULTS: The dosimetric comparisons on the homogeneous phantom between the proposed PB algorithm and the CCC algorithm for 10 VMAT arcs demonstrate that the proposed algorithm can achieve a dosimetric accuracy comparable to that of the CCC algorithm with average gamma passing rates of 96% (2%/2mm) and 98% (3%/3mm). In addition, the proposed algorithm can provide better computational efficiency for VMAT dose computation using a PC equipped with a 4-core processor, compared to the CCC algorithm utilizing a dual 10-core server. Moreover, the computation complexity theory reveals that the proposed algorithm has a great advantage with regard to computational efficiency for VMAT dose computation on homogeneous medium, especially when a fine angular sampling rate is applied. This can support a reduction in dose errors from the angular under-sampling effect by using a finer angular sampling rate, while still preserving a practical computing speed. For dose calculation on the heterogeneous phantom, the proposed algorithm with heterogeneity corrections can still offer a reasonable dosimetric accuracy with comparable computational efficiency to that of the CCC algorithm. CONCLUSIONS: We proposed a novel polar-coordinate-based pencil beam algorithm for VMAT dose computation that enables a better computational efficiency while maintaining clinically acceptable dosimetric accuracy and reducing dose error caused by the angular under-sampling effect. It also provides a flexible VMAT dose computation structure that allows adjustable sampling rates and direct dose computation in regions of interest, which makes the algorithm potentially useful for clinical applications such as independent dose verification for VMAT patient-specific QA.

Investigation of effect of filter on the stand-up technique for total skin irradiation by Monte Carlo simulation.

PURPOSE: The aim of this study was to investigate dosimetric effects of scattering filter on the stand-up technique for total skin irradiation (TSI) with a single electron field by Monte Carlo (MC) simulation. METHODS: MC simulations were performed with BEAMnrc and DOSXYZnrc packages under EGSnrc environment. Scattering filter of a metal disc was mounted in the accessory slot. The filter materials (Cu, Fe, Au, Zn, Ag) were investigated, with thickness ranging from 0.05 to 0.55 mm, depending on material. The extended source to skin distance (SSD) ranging from 250 to 350 cm was studied. The following dosimetric quantities were evaluated: percent depth dose (PDD), profiles and output factor at depth of maximum, and composite dose distribution on a 30-cm diameter cylindrical phantom. They were compared with the standard dual beam technique used at our clinic. The effects on different patient sizes were also studied. RESULTS: No filter produced acceptable profile flatness (±10% within the central 160 cm) at 250 cm SSD. At 300 cm SSD, Au (0.1 mm), Ag (0.25 mm), Cu (0.5 mm) produced acceptable flatness while Zn (0.45 mm) required 325 cm SSD. For these four configurations, the dmax was 0.90-0.99 cm, similar to dual beam (0.97 cm); R50 was 1.85-1.91 cm, compared with dual beam of 2.06 cm; the output factor ranged from 0.025 to 0.029, lower than the dual beam (0.080). With the composite fields for four configurations, the dmax was 0.10 cm, compared with dual beam (0.16 cm). The surface dose was 97%, similar to dual beam (96%). B-factor was 3.3-3.4, compared with dual beam of 3.1. The maximum X-ray contamination was 3%, higher than dual beam (1%). CONCLUSIONS: The investigation suggests the TSI stand-up technique can be implemented using a single electron beam if a customized filter is used. More dosimetric measurements are needed to validate the MC results and clinical implementation.

Validation of the dosimetry of total skin irradiation techniques by Monte Carlo simulation.

PURPOSE: To validate the dose measurements for two total skin irradiation techniques with Monte Carlo simulation, providing more information on dose distributions, and guidance on further technique optimization. METHODS: Two total skin irradiation techniques (stand-up and lay-down) with different setup were simulated and validated. The Monte Carlo simulation was primarily performed within the EGSnrc environment. Parameters of jaws, MLCs, and a customized copper (Cu) filter were first tuned to match the profiles and output measured at source-to-skin distance (SSD) of 100 cm where the secondary source is defined. The secondary source was rotated to simulate gantry rotation. VirtuaLinac, a cloud-based Monte Carlo package, was used for Linac head simulation as a secondary validation. The following quantities were compared with measurements: for each field/direction at the treatment SSDs, the percent depth dose (PDD), the profiles at the depth of maximum, and the absolute dosimetric output; the composite dose distribution on cylindrical phantoms of 20 to 40 cm diameters. RESULTS: Cu filter broadened the FWHM of the electron beam by 44% and degraded the mean energy by 0.7 MeV. At SSD = 100 cm, MC calculated PDDs agreed with measured data within 2%/2 mm (except for the surface voxel) and lateral profiles agreed within 3%. At the treatment SSD, profiles and output factors of individual field matched within 4%; dmax and R80 of the simulated PDDs also matched with measurement within 2 mm. When all fields were combined on the cylindrical phantom, the dmax shifted toward the surface. For lay-down technique, the maximum x-ray contamination at the central axis was (MC: 2.2; Measurement: 2.1)% and reduced to 0.2% at 40 cm off the central axis. CONCLUSIONS: The Monte Carlo results in general agree well with the measurement, which provides support in our commissioning procedure, as well as the full three-dimensional dose distribution of the patient phantom.