Evaluation of an MRI linac magnetic isocenter walkout with gantry rotation in the presence of angle-specific corrections.

Objective. To reduce the magnetic isocenter position variation with gantry rotation on an 0.35 T MRI-guided linac to a practically negligible level.Approach. Central fRequency (CF) offset, eddy current calibration, cross-term calibration, gradient delay, and gradient offsets are tuned for each MR linac installation at every 30° of gantry rotation and stored in a look-up table (LUT). During treatment, the CF is tuned only once in the beginning at an arbitrary gantry angle. After that, imaging paramters are offset based on the stored LUT values for any given gantry angle.Main results. For the same hardware configuration, the implementation of the gantry-angle-specific parameter corrections reduced the total isocenters range of travel in the transverse plane from 1.1 to 0.3 mm and from 0.8 to 0.2 mm in horizontal and vertical directions, respectively. With the longitudinal shift always being negligible (≤0.2 mm), the radius of the sphere encompassing the isocenter locations was reduced from 0.6 to 0.2 mm. Geometric distortion improved as well; in particular, the gantry-angle-averaged maximum longitudinal distortion within a 35 cm diameter sphere was reduced from 1.4 to 0.8 mm. Since the CF is tuned only once during treatment, imaging may resume promptly after the gantry reaches the next target position.Significance. The MRI-guided linear accelerator was conceived primarily as an instrument for precision image-guided therapy. Thus, it is important to keep the treatment and imaging isocentres as close as possible while minimizing the geometric distortion. The described solution reduces the walkout of the imaging isocenter to a fraction of 1 mm, while keeping geometric distortion in a substantial volume below 1 mm. The approach is robust and does not increase the overall procedure time.

An independent Monte Carlo-based IMRT QA tool for a 0.35 T MRI-guided linear accelerator.

PURPOSE: To develop an independent log file-based intensity-modulated radiation therapy (IMRT) quality assurance (QA) tool for the 0.35 T magnetic resonance-linac (MR-linac) and investigate the ability of various IMRT plan complexity metrics to predict the QA results. Complexity metrics related to tissue heterogeneity were also introduced. METHODS: The tool for particle simulation (TOPAS) Monte Carlo code was utilized with a previously validated linac head model. A cohort of 29 treatment plans was selected for IMRT QA using the developed QA tool and the vendor-supplied adaptive QA (AQA) tool. For 27 independent patient cases, various IMRT plan complexity metrics were calculated to assess the deliverability of these plans. A correlation between the gamma pass rates (GPRs) from the AQA results and calculated IMRT complexity metrics was determined using the Pearson correlation coefficients. Tissue heterogeneity complexity metrics were calculated based on the gradient of the Hounsfield units. RESULTS: The median and interquartile range for the TOPAS GPRs (3%/3 mm criteria) were 97.24% and 3.75%, respectively, and were 99.54% and 0.36% for the AQA tool, respectively. The computational time for TOPAS ranged from 4 to 8 h to achieve a statistical uncertainty of <1.5%, whereas the AQA tool had an average calculation time of a few minutes. Of the 23 calculated IMRT plan complexity metrics, the AQA GPRs had correlations with 7 out of 23 of the calculated metrics. Strong correlations (|r| > 0.7) were found between the GPRs and the heterogeneity complexity metrics introduced in this work. CONCLUSIONS: An independent MC and log file-based IMRT QA tool was successfully developed and can be clinically deployed for offline QA. The complexity metrics will supplement QA reports and provide information regarding plan complexity.

A multi-institutional comparison of dosimetric data for a 0.35 T MR-linac.

Objective.A comparison of percent depth dose (PDD) curves, lateral beam profiles, output factors (OFs), multileaf collimator (MLC) leakage, and couch transmission factors was performed between ten institutes for a commercial 0.35 T MR-linac.Approach.The measured data was collected during acceptance testing of the MR-linac. The PDD curves were measured for the 3.32 × 3.32 cm2, 9.96 × 9.96 cm2, and 27.20 × 24.07 cm2field sizes. The lateral beam profiles were acquired for a 27.20 × 24.07 cm2field size using an ion chamber array and penumbra was defined as the distance between 80% of the maximum dose and 20% of the maximum dose after normalizing the profiles to the dose at the inflection points. The OFs were measured using solid-state dosimeters, whereas radiochromic films were utilized to measure radiation leakage through the MLC stacks. The relative couch transmission factors were measured for various gantry angles. The variation in the multi-institutional data was quantified using the percent standard deviation metric.Main results.Minimal variations (<1%) were found between the PDD data, except for the build-up region and the deeper regions of the PDD curve. The in-field region of the lateral beam profiles varied <1.5% between different institutions and a small variation (<0.7 mm) in penumbra was observed. A variation of <1% was observed in the OF data for field sizes above 1.66 × 1.66 cm2, whereas large variations were shown for small-field sizes. The average and maximum MLC leakage was calculated to be <0.3% and <0.6%, which was well below the international electrotechnical commission (IEC) leakage thresholds. The couch transmission was smallest for oblique beams and ranged from 0.83 to 0.87.Significance.The variation in the data was found to be relatively small and the different 0.35 T MR-linacs were concluded to have similar dosimetric characteristics.

On the MLC leaves alignment in the direction orthogonal to movement.

The main focus of the recommended spatial accuracy tests for the multi-leaf collimators (MLC) is calibration of the leaf position along the movement direction and overall alignment to the radiation isocenter. No explicit attention was typically paid to the alignment of the leaves from the opposing banks in the direction orthogonal to movement. This paper is a case study demonstrating that verification of such alignment at the time of acceptance testing is prudent. The original standard MLC (SMLC) on an MRIdian MRI-guided linac (ViewRay Inc., Mountain View, CA, USA) was upgraded to a high-speed MLC (HSMLC), which is supposed to be mechanically identical to the SMLC except for the higher drive screw pitch. The results of the end-to-end IMRT tests demonstrated unacceptable dosimetric results exemplified by an average and maximum ion chamber (IC) point dose error in the high-dose low-gradient region of 2.5 ± 1.4% and 4.6%, respectively. Before the upgrade, those values were 0.3 ± 0.7% and 0.9%, respectively. An exhaustive analysis of possible failure modes eventually zeroed in on the average misalignment of about 1 mm in the Y (along the couch) direction between the right and left upper MLC banks. The MLC was replaced, reducing the Y-direction misalignment to 0.4 mm. As a result, the average and maximum IC dose-errors became acceptable 1.0 ± 0.7% and 1.6%, respectively. Simple film and/or chamber array tests during acceptance testing can easily detect Y-direction misalignments between opposing leaves banks measuring a fraction of a mm at isocenter. Left undetected, such misalignment can cause nontrivial dosimetric consequences.

Development and evaluation of a GEANT4-based Monte Carlo Model of a 0.35 T MR-guided radiation therapy (MRgRT) linear accelerator.

PURPOSE: The aim of this work was to develop and benchmark a magnetic resonance (MR)-guided linear accelerator head model using the GEANT4 Monte Carlo (MC) code. The validated model was compared to the treatment planning system (TPS) and was also used to quantify the electron return effect (ERE) at a lung-water interface. METHODS: The average energy, including the spread in the energy distribution, and the radial intensity distribution of the incident electron beam were iteratively optimized in order to match the simulated beam profiles and percent depth dose (PDD) data to measured data. The GEANT4 MC model was then compared to the TPS model using several photon beam tests including oblique beams, an off-axis aperture, and heterogeneous phantoms. The benchmarked MC model was utilized to compute output factors (OFs) with the 0.35 T magnetic field turned on and off. The ERE was quantified at a lung-water interface by simulating PDD curves with and without the magnetic field for 6.6 × 6.6  cm 2 and 2.5 × 2.5  cm 2 field sizes. A 2%/2 mm gamma criterion was used to compare the MC model with the TPS data throughout this study. RESULTS: The final incident electron beam parameters were 6.0 MeV average energy with a 1.5 MeV full width at half maximum (FWHM) Gaussian energy spread and a 1.0 mm FWHM Gaussian radial intensity distribution. The MC-simulated OFs were found to be in agreement with the TPS-calculated and measured OFs, and no statistical difference was observed between the 0.35 T and 0.0 T OFs. Good agreement was observed between the TPS-calculated and MC-simulated data for the photon beam tests with gamma pass rates ranging from 96% to 100%. An increase of 4.3% in the ERE was observed for the 6.6 × 6.6  cm 2 field size relative to the 2.5 × 2.5  cm 2 field size. The ratio of the 0.35 T PDD to the 0.0 T PDD was found to be up to 1.098 near lung-water interfaces for the 6.6 × 6.6  cm 2 field size using the MC model. CONCLUSIONS: A vendor-independent Monte Carlo model has been developed and benchmarked for a 0.35 T/6 MV MR-linac. Good agreement was obtained between the GEANT4 and TPS models except near heterogeneity interfaces.

Technical Note: Validation of an automatic ACR phantom quality assurance tool for an MR-guided radiotherapy system.

PURPOSE: To modify and evaluate an automatic American College of Radiology (ACR) phantom analysis toolbox for ACR quality assurance (QA) on a low-field MR-guided radiotherapy system (ViewRay). METHODS: An open-source toolbox was modified for ACR QA of a 0.35T MRI system (ViewRay MRIdian). A total of 17 ACR datasets were evaluated, including 10 datasets acquired from different systems across the world, and seven datasets acquired at our center between 2014 and 2020. All required ACR tests, geometric accuracy (GA), high-contrast spatial resolution (HCSR), slice thickness accuracy (ST), slice position accuracy (SP), percent integral uniformity (PIU), percentage signal ghosting (PSG), and low-contrast object detectability (LCOD), were assessed manually and using the toolbox automatically. Measurements between manual and automatic analysis were compared. Precision, recall, and accuracy were calculated, where the manual results were used as the ground truth. RESULTS: The software took less than 2 min to complete all seven tests, which usually requires 40 min or more if analyzed manually. Overall, the automatic measurement was consistent with the manual result. The absolute differences between the two measurements were 0.72 ± 0.66 (mm), 0.01 ± 0.03, 0.50 ± 0.60 (mm), 0.39 ± 0.41 (mm), 1.01 ± 1.00 (%), 0.0016 ± 0.0019, and 2.79 ± 2.29 for the seven tests. The precision of the automatic toolbox was 100% for all tests, indicating that a test would 100% pass a manual analysis if it had passed the automatic analysis. Recall and accuracy were ≥ 96% for GA, HCSR, SP, PIU, PSG, and LCOD tests, and 91% for the ST test. CONCLUSION: An automatic ACR QA tool was adopted and evaluated for the low-field MR-guided radiotherapy (MRgRT) system. Overall, the toolbox provided comparable results as manual analysis, and reduced the processing time from over 40 min to <2 min. This toolbox holds the potential to be widely adopted either as a second check tool or partially replace human measurement for MRgRT programs using the same system.

Clinical experience of MRI4D QUASAR motion phantom for latency measurements in 0.35T MR-LINAC.

PURPOSE: In MRgRT, accuracy of treatment depends on the gating latency, when real-time targeting and gating is enabled. Gating latency is dependent on image acquisition, processing time, accuracy, efficacy of target tracking algorithms, and radiation beam delivery latency. In this report, clinical experience of the MRI4D QUASAR motion phantom for latency measurements on a 0.35-T magnetic resonance-linear accelerator (MR-LINAC) with two imaging speeds and four tracking algorithms was studied. MATERIALS/METHODS: Beam-control latency was measured on a 0.35-T MR-LINAC system with four target tracking algorithms and two real-time cine imaging sequences [four and eight frames per second (FPS)]. Using an MR-compatible motion phantom, the delays between phantom beam triggering signal and linac radiation beam control signal were evaluated for three motion periods with a rigid target. The gating point was set to be 8 mm above the full exhalation position. The beam-off latency was measured for a total of 24 combinations of tracking algorithm, imaging FPS, and motion periods. The corresponding gating target margins were determined using the target motion speed multiplied by the beam-off latency. RESULTS: The largest measured beam-off latency was 302 ± 20 ms with the Large Deforming Targets (LDT) algorithm and 4 s motion period imaged with 8-FPS cine MRI. The corresponding gating uncertainty based on target motion speed was 3.0 mm. The range of the average beam-off latency was 128-243 ms in 4-FPS imaging and 47-302 ms in 8-FPS imaging. CONCLUSIONS: The gating latency was measured using an MRI4D QUASAR motion phantom in a 0.35-T MR-LINAC. The latency measurements include time delay related to MR imaging method, target tracking algorithm and system delay. The gating uncertainty was estimated based on the beam-off latency measurements and the target motion.

Characterization of radiotherapy component impact on MR imaging quality for an MRgRT system.

Radiotherapy components of an magnetic resonnace-guided radiotherapy (MRgRT) system can alter the magnetic fields, causing spatial distortion and image deformation, altering imaging and radiation isocenter coincidence and the accuracy of dose calculations. This work presents a characterization of radiotherapy component impact on MR imaging quality in terms of imaging isocenter variation and spatial integrity changes on a 0.35T MRgRT system, pre- and postupgrade of the system. The impact of gantry position, MLC field size, and treatment table power state on imaging isocenter and spatial integrity were investigated. A spatial integrity phantom was used for all tests. Images were acquired for gantry angles 0-330° at 30° increments to assess the impact of gantry position. For MLC and table power state tests all images were acquired at the home gantry position (330°). MLC field sizes ranged from 1.66 to 27.4 cm edge length square fields. Imaging isocenter shift caused by gantry position was reduced from 1.7 mm at gantry 150° preupgrade to 0.9 mm at gantry 120° postupgrade. Maximum spatial integrity errors were 0.5 mm or less pre- and postupgrade for all gantry angles, MLC field sizes, and treatment table power states. However, when the treatment table was powered on, there was significant reduction in SNR. This study showed that gantry position can impact imaging isocenter, but spatial integrity errors were not dependent on gantry position, MLC field size, or treatment table power state. Significant isocenter variation, while reduced postupgrade, is cause for further investigation.