Magnetic fields are causing small, but significant changes of the radiochromic EBT3 film response to 6 MV photons

Assessment of the deep learning-based gamma passing rate prediction system for 1.5 T magnetic resonance-guided linear accelerator

Measurement-based verification is impossible for the patient-specific quality assurance (QA) of online adaptive magnetic resonance imaging-guided radiotherapy (oMRgRT) because the patient remains on the couch throughout the session. We assessed a deep learning (DL) system for oMRgRT to predict the gamma passing rate (GPR). This study collected 125 verification plans [reference plan (RP), 100; adapted plan (AP), 25] from patients with prostate cancer treated using Elekta Unity. Based on our previous study, we employed a convolutional neural network that predicted the GPRs of nine pairs of gamma criteria from 1%/1 mm to 3%/3 mm. First, we trained and tested the DL model using RPs (n = 75 and n = 25 for training and testing, respectively) for its optimization. Second, we tested the GPR prediction accuracy using APs to determine whether the DL model could be applied to APs. The mean absolute error (MAE) and correlation coefficient (r) of the RPs were 1.22 ± 0.27% and 0.29 ± 0.10 in 3%/2 mm, 1.35 ± 0.16% and 0.37 ± 0.15 in 2%/2 mm, and 3.62 ± 0.55% and 0.32 ± 0.14 in 1%/1 mm, respectively. The MAE and r of the APs were 1.13 ± 0.33% and 0.35 ± 0.22 in 3%/2 mm, 1.68 ± 0.47% and 0.30 ± 0.11 in 2%/2 mm, and 5.08 ± 0.29% and 0.15 ± 0.10 in 1%/1 mm, respectively. The time cost was within 3 s for the prediction. The results suggest the DL-based model has the potential for rapid GPR prediction in Elekta Unity.

Proton dosimetry in a magnetic field: Measurement and calculation of magnetic field correction factors for a plane-parallel ionization chamber

BACKGROUND

The combination of magnetic resonance imaging and proton therapy offers the potential to improve cancer treatment. The magnetic field (MF)-dependent change in the dosage of ionization chambers in magnetic resonance imaging-integrated proton therapy (MRiPT) is considered by the correction factork B ⃗ , M , Q $k_{\vec{B},M,Q}$ , which needs to be determined experimentally or computed via Monte Carlo (MC) simulations.

PURPOSE

In this study,k B ⃗ , M , Q $k_{\vec{B},M,Q}$ was both measured and simulated with high accuracy for a plane-parallel ionization chamber at different clinical relevant proton energies and MF strengths.

MATERIAL AND METHODS

The dose-response of the Advanced Markus chamber (TM34045, PTW, Freiburg, Germany) irradiated with homogeneous 10 × $\times$ 10 cm2 $^2$ quasi mono-energetic fields, using 103.3, 128.4, 153.1, 223.1, and 252.7 MeV proton beams was measured in a water phantom placed in the MF of an electromagnet with MF strengths of 0.32, 0.5, and 1 T. The detector was positioned at a depth of 2 g/cm2 $^2$ in water, with chamber electrodes parallel to the MF lines and perpendicular to the proton beam incidence direction. The measurements were compared with TOPAS MC simulations utilizing COMSOL-calculated 0.32, 0.5, and 1 T MF maps of the electromagnet.k B ⃗ , M , Q $k_{\vec{B},M,Q}$ was calculated for the measurements for all energies and MF strengths based on the equation:k B ⃗ , M , Q = M Q M Q B ⃗ $k_{\vec{B},M,Q}=\frac{M_\mathrm{Q}}{M_\mathrm{Q}^{\vec{B}}}$ , whereM Q B ⃗ $M_\mathrm{Q}^{\vec{B}}$ andM Q $M_\mathrm{Q}$ were the temperature and air-pressure corrected detector readings with and without the MF, respectively. MC-based correction factors were determined ask B ⃗ , M , Q = D det D det B ⃗ $k_{\vec{B},M,Q}=\frac{D_\mathrm{det}}{D_\mathrm{det}^{\vec{B}}}$ , whereD det B ⃗ $D_\mathrm{det}^{\vec{B}}$ andD det $D_\mathrm{det}$ were the doses deposited in the air cavity of the ionization chamber model with and without the MF, respectively. Furthermore, MF effects on the chamber dosimetry are studied using MC simulations, examining the impact on the absorbed dose-to-water (D W $D_{W}$ ) and the shift in depth of the Bragg peak.

RESULTS

The detector showed a reduced dose-response for all measured energies and MF strengths, resulting in experimentally determinedk B ⃗ , M , Q $k_{\vec{B},M,Q}$ values larger than unity. For all energies and MF strengths examined,k B ⃗ , M , Q $k_{\vec{B},M,Q}$ ranged between 1.0065 and 1.0205. The dependence on the energy and the MF strength was found to be non-linear with a maximum at 1 T and 252.7 MeV. The MC simulatedk B ⃗ , M , Q $k_{\vec{B},M,Q}$ values agreed with the experimentally determined correction factors within their standard deviations with a maximum difference of 0.6%. The MC calculated impact onD W $D_{W}$ was smaller 0.2 %.

CONCLUSION

For the first time, measurements and simulations were compared for proton dosimetry within MFs using an Advanced Markus chamber. Good agreement ofk B ⃗ , M , Q $k_{\vec{B},M,Q}$ was found between experimentally determined and MC calculated values. The performed benchmarking of the MC code allows for calculatingk B ⃗ , M , Q $k_{\vec{B},M,Q}$ for various ionization chamber models, MF strengths and proton energies to generate the data needed for a proton dosimetry protocol within MFs and is, therefore, a step towards MRiPT.

Model-based machine learning for the recovery of lateral dose profiles of small photon fields in magnetic field

Objective. To investigate the feasibility to train artificial neural networks (NN) to recover lateral dose profiles from detector measurements in a magnetic field. Approach. A novel framework based on a mathematical convolution model has been proposed to generate measurement-less training dataset. 2D dose deposition kernels and detector lateral fluence response functions of two air-filled ionization chambers and two diode-type detectors have been simulated without magnetic field and for magnetic field B = 0.35 and 1.5 T. Using these convolution kernels, training dataset consisting pairs of dose profiles D x , y and signal profiles M x , y were computed for a total of 108 2D photon fluence profiles ψ ( x , y ) (80% training/20% validation). The NN were tested using three independent datasets, where the second test dataset has been obtained from simulations using realistic phase space files of clinical linear accelerator and the third test dataset was measured at a conventional linac equipped with electromagnets. Main results. The convolution kernels show magnetic field dependence due to the influence of the Lorentz force on the electron transport in the water phantom and detectors. The NN show good performance during training and validation with mean square error reaching a value of 1e-6 or smaller. The corresponding correlation coefficients R reached the value of 1 for all models indicating an excellent agreement between expected D x , y and predicted D pred x , y . The comparisons between D x , y and D pred x , y using the three test datasets resulted in gamma indices (1 mm/1% global) <1 for all evaluated data points. Significance. Two verification approaches have been proposed to warrant the mathematical consistencies of the NN outputs. Besides offering a correction strategy not existed so far for relative dosimetry in a magnetic field, this work could help to raise awareness and to improve understanding on the distortion of detector’s signal profiles by a magnetic field.

Characterizing magnetically focused contamination electrons by off-axis irradiation on an inline MRI-Linac

Purpose

The aim of this study is to investigate off‐axis irradiation on the Australian MRI‐Linac using experiments and Monte Carlo simulations. Simulations are used to verify experimental measurements and to determine the minimum offset distance required to separate electron contamination from the photon field.

Methods

Dosimetric measurements were performed using a microDiamond detector, Gafchromic^® EBT3 film, and MOSkin^TM. Three field sizes were investigated including 1.9 × 1.9, 5.8 × 5.8, and 9.7 × 9.6 cm². Each field was offset a maximum distance, approximately 10 cm, from the central magnetic axis (isocenter). Percentage depth doses (PDDs) were collected at a source‐to‐surface distance (SSD) of 1.8 m for fields collimated centrally and off‐axis. PDD measurements were also acquired at isocenter for each off‐axis field to measure electron contamination. Monte Carlo simulations were used to verify experimental measurements, determine the minimum field offset distance, and demonstrate the use of a spoiler to absorb electron contamination.

Results

Off‐axis irradiation separates the majority of electron contamination from an x‐ray beam and was found to significantly reduce in‐field surface dose. For the 1.9 × 1.9, 5.8 × 5.8, and 9.7 × 9.6 cm² field, surface dose was reduced from 120.9% to 24.9%, 229.7% to 39.2%, and 355.3% to 47.3%, respectively. Monte Carlo simulations generally were within experimental error to MOSkin^TM and microDiamond, and used to determine the minimum offset distance, 2.1 cm, from the field edge to isocenter. A water spoiler 2 cm thick was shown to reduce electron contamination dose to near zero.

Conclusions

Experimental and simulation data were acquired for a range of field sizes to investigate off‐axis irradiation on an inline MRI‐Linac. The skin sparing effect was observed with off‐axis irradiation, a feature that cannot be achieved to the same extent with other methods, such as bolusing, for beams at isocenter.

Traceable reference dosimetry in MRI guided radiotherapy using alanine: calibration and magnetic field correction factors of ionisation chambers

Magnetic resonance imaging (MRI)-guided radiotherapy (RT) (MRIgRT) falls outside the scope of existing high energy photon therapy dosimetry protocols, because those protocols do not consider the effects of the magnetic field on detector response and on absorbed dose to water. The aim of this study is to evaluate and demonstrate the traceable measurement of absorbed dose in MRIgRT systems using alanine, made possible by the characterisation of alanine sensitivity to magnetic fields reported previously by Billaset al(2020Phys. Med. Biol.65115001), in a way which is compatible with existing standards and calibrations available for conventional RT. In this study, alanine is used to transfer absorbed dose to water to MRIgRT systems from a conventional linac. This offers an alternative route for the traceable measurement of absorbed dose to water, one which is independent of the transfer using ionisation chambers. The alanine dosimetry is analysed in combination with measurements with several Farmer-type chambers, PTW 30013 and IBA FC65-G, at six different centres and two different MRIgRT systems (Elekta Unity™ and ViewRay MRIdian™). The results are analysed in terms of the magnetic field correction factors, and in terms of the absorbed dose calibration coefficients for the chambers, determined at each centre. This approach to reference dosimetry in MRIgRT produces good consistency in the results, across the centres visited, at the level of 0.4% (standard deviation). Farmer-type ionisation chamber magnetic field correction factors were determined directly, by comparing calibrations in some MRIgRT systems with and without the magnetic field ramped up, and indirectly, by comparing calibrations in all the MRIgRT systems with calibrations in a conventional linac. Calibration coefficients in the MRIgRT systems were obtained with a standard uncertainty of 1.1% (Elekta Unity™) and 0.9% (ViewRay MRIdian™), for three different chamber orientations with respect to the magnetic field. The values obtained for the magnetic field correction factor in this investigation are consistent with those presented in the summary by de Pooteret al(2021Phys. Med. Biol.6605TR02), and would tend to support the adoption of a magnetic field correction factor which depends on the chamber type, PTW 30013 or IBA FC65-G.

The role of the construction and sensitive volume of compact ionization chambers on the magnetic field-dependent dose response

PURPOSE

The magnetic-field correction factors k B , Q of compact air-filled ionization chambers have been investigated experimentally and using Monte Carlo simulations up to 1.5 T. The role of the nonsensitive region within the air cavity and influence of the chamber construction on its dose response have been elucidated.

MATERIALS AND METHODS

The PTW Semiflex 3D 31021, PinPoint 3D 31022, and Sun Nuclear Cooperation SNC125c chambers were studied. The k B , Q factors were measured at the experimental facility of the German National Metrology Institute (PTB) up to 1.4 T using a 6 MV photon beam. The chambers were positioned with the chamber axis perpendicular to the beam axis (radial); and parallel to the beam axis (axial). In both cases, the magnetic field was directed perpendicular to both the beam axis and chamber axis. Additionally, the sensitive volumes of these chambers have been experimentally determined using a focused proton microbeam and finite element method. Beside the simulations of k B , Q factors, detailed Monte Carlo technique has been applied to analyse the secondary electron fluence within the air cavity, that is, the number of secondary electrons and the average path length as a function of the magnetic field strength.

RESULTS

A nonsensitive volume within the air cavity adjacent to the chamber stem for the PTW chambers has been identified from the microbeam measurements and FEM calculations. The dose response of the three investigated ionization chambers does not deviate by more than 4% from the field-free case within the range of magnetic fields studied in this work for both the radial and axial orientations. The simulated k B , Q for the fully guarded PTW chambers deviate by up to 6% if their sensitive volumes are not correctly considered during the simulations. After the implementation of the sensitive volume derived from the microbeam measurements, an agreement of better than 1% between the experimental and Monte Carlo k B , Q factors for all three chambers can be achieved. Detailed analysis reveals that the stem of the PTW chambers could give rise to a shielding effect reducing the number of secondary electrons entering the air cavity in the presence of magnetic field. However, the magnetic field dependence of their path length within the air cavity is shown to be weaker than for the SNC125c chamber, where the length of the air cavity is larger than its diameter. For this chamber it is shown that the number of electrons and their path lengths in the cavity depend stronger on the magnetic field.

DISCUSSION AND CONCLUSION

For clinical measurements up to 1.5 T, the required k B , Q corrections of the three chambers could be kept within 3% in both the investigated chamber orientations. The results reiterate the importance of considering the sensitive volume of fully guarded chambers, even for the investigated compact chambers, in the Monte Carlo simulations of chamber response in magnetic field. The resulting magnetic field-dependent dose response has been demonstrated to depend on the chamber construction, such as the ratio between length and the diameter of the air cavity as well as the design of the chamber stem.

Clinical utility of Gafchromic film in an MRI-guided linear accelerator

Background

The purpose of this study is to comprehensively evaluate the suitability of Gafchromic EBT3 and EBT-XD film for dosimetric quality assurance in 0.35 T MR-guided radiotherapy.

Methods

A 0.35 T magnetic field strength was utilized to evaluate magnetic field effects on EBT3 and EBT-XD Gafchromic films by studying the effect of film exposure time within the magnetic field using two timing sequences and film not exposed to MR, the effect of magnetic field exposure on the crystalline structure of the film, and the effect of orientation of the film with respect to the bore within the magnetic field. The orientation of the monomer crystal was qualitatively evaluated using scanning electron microscopy (SEM) compared to unirradiated film. Additionally, dosimetric impact was evaluated through measurements of a series of open field irradiations (0.83 × 0.83-cm² to 19.92 × 19.92-cm²) and patient specific quality assurance measurements. Open fields were compared to planned dose and an independent dosimeter. Film dosimetry was applied to twenty conventional and twenty stereotactic body radiotherapy (SBRT) patient specific quality assurance cases.

Results

No visual changes in crystal orientation were observed in any evaluated SEM images nor were any optical density differences observed between films irradiated inside or outside the magnetic field for both EBT3 and EBT-XD film. At small field sizes, the average difference along dose profiles measured in film compared to the same points measured using an independent dosimeter and to predicted treatment planning system values was 1.23% and 1.56%, respectively. For large field sizes, the average differences were 1.91% and 1.21%, respectively. In open field tests, the average gamma pass rates were 99.8% and 97.2%, for 3%/3 mm and 3%/1 mm, respectively. The median (interquartile range) 3%/3 mm gamma pass rates in conventional QA cases were 98.4% (96.3 to 99.2%), and 3%/1 mm in SBRT QA cases were 95.8% (95.0 to 97.3%).

Conclusions

MR exposure at 0.35 T had negligible effects on EBT3 and EBT-XD Gafchromic film. Dosimetric film results were comparable to planned dose, ion chamber and diode measurements.

Initial clinical experience of patient-specific QA of treatment delivery in online adaptive radiotherapy using a 1.5 T MR-Linac

Purpose. This study aims to evaluate the performance of a commercial 1.5 T MR-Linac by analyzing its patient-specific quality assurance (QA) data collected during one full year of clinical operation.Methods and Materials. The patient-specific QA system consisted of offline delivery QA (DQA) and online calculation-based QA. Offline DQA was based on ArcCHECK-MR combined with an ionization chamber. Online QA was performed using RadCalc that calculated and compared the point dose calculation with the treatment planning system (TPS). A total of 24 patients with 189 treatment fractions were enrolled in this study. Gamma analysis was performed and the threshold that encompassed 95% of QA results (T95) was reported. The plan complexity metric was calculated for each plan and compared with the dose measurements to determine whether any correlation existed.Results. All point dose measurements were within 5% deviation. The mean gamma passing rates of the group data were found to be 96.8 ± 4.0% and 99.6 ± 0.7% with criteria of 2%/2mm and 3%/3mm, respectively. T95 of 87.4% and 98.2% was reported for the overall group with the two passing criteria, respectively. No statistically significant difference was found between adaptive treatments with adapt-to-position (ATP) and adapt-to-shape (ATS), whilst the category of pelvis data showed a better passing rate than other sites. Online QA gave a mean deviation of 0.2 ± 2.2%. The plan complexity metric was positively correlated with the mean dose difference whilst the complexity of the ATS cohort had larger variations than the ATP cohort.Conclusions. A patient-specific QA system based on ArcCHECK-MR, solid phantom and ionization chamber has been well established and implemented for validation of treatment delivery of a 1.5 T MR-Linac. Our QA data obtained over one year confirms that good agreement between TPS calculation and treatment delivery was achieved.

Effects on skin dose from unwanted air gaps under bolus in an MR-guided linear accelerator (MR-linac) system

Bolus is commonly used in MV photon radiotherapy to increase superficial dose and improve dose uniformity for treating shallow lesions. However, irregular patient body contours can cause unwanted air gaps between a bolus and patient skin. The resulting dosimetric errors could be exacerbated in MR-Linac treatments, as secondary electrons generated by photons are affected by the magnetic field. This study aimed to quantify the dosimetric effect of unwanted gaps between bolus and skin surface in an MR-Linac. A parallel-plate ionization chamber and EBT3 films were utilized to evaluate the surface dose under bolus with various gantry angles, field sizes, and different air gaps. The results of surface dose measurements were then compared to Monaco 5.40 Treatment Planning System (TPS) calculations. The suitability of using a parallel-plate chamber in MR-Linac measurement was validated by benchmarking the percentage depth dose and output factors with the microDiamond detector and air-filled ionization chamber measurements in water. A non-symmetric response of the parallel-plate chamber to oblique beams in the magnetic field was characterized. Unwanted air gaps significantly reduced the skin dose. For a frontal beam, skin dose was halved when there was a 5 mm gap, a much larger difference than in a conventional linac. Skin dose manifested a non-symmetric pattern in terms of gantry angle and gap size. The TPS overestimated skin dose in general, but shared the same trend with measurement when there was no air gap, or the gap size was larger than 5 mm. However, the calculated and measured results had a large discrepancy when the bolus-skin gap was below 5 mm. When treating superficial lesions, unwanted air gaps under the bolus will compromise the dosimetric goals. Our results highlight the importance of avoiding air gaps between bolus and skin when treating superficial lesions using an MR-Linac system.

Reference dosimetry in MRI-linacs: evaluation of available protocols and data to establish a Code of Practice

With the rapid increase in clinical treatments with MRI-linacs, a consistent, harmonized and sustainable ground for reference dosimetry in MRI-linacs is needed. Specific for reference dosimetry in MRI-linacs is the presence of a strong magnetic field. Therefore, existing Code of Practices (CoPs) are inadequate. In recent years, a vast amount of papers have been published in relation to this topic. The purpose of this review paper is twofold: to give an overview and evaluate the existing literature for reference dosimetry in MRI-linacs and to discuss whether the literature and datasets are adequate and complete to serve as a basis for the development of a new or to extend existing CoPs. This review is prefaced with an overview of existing MRI-linac facilities. Then an introduction on the physics of radiation transport in magnetic fields is given. The main part of the review is devoted to the evaluation of the literature with respect to the following subjects: • beam characteristics of MRI-linac facilities; • formalisms for reference dosimetry in MRI-linacs; • characteristics of ionization chambers in the presence of magnetic fields; • ionization chamber beam quality correction factors; and • ionization chamber magnetic field correction factors. The review is completed with a discussion as to whether the existing literature is adequate to serve as basis for a CoP. In addition, it highlights subjects for future research on this topic.

Experimental determination of magnetic field correction factors for ionization chambers in parallel and perpendicular orientations

Magnetic field correction factors are needed for absolute dosimetry in magnetic resonance (MR)-linacs. Currently experimental data for magnetic field correction factors, especially for small volume ionization chambers, are largely lacking. The purpose of this work is to establish, independent methods for the experimental determination of magnetic field correction factors [Formula: see text] in an orientation in which the ionization chamber is parallel to the magnetic field. The aim is to confirm previous experiments on the determination of Farmer type ionization chamber correction factors and to gather information about the usability of small-volume ionization chambers for absolute dosimetry in MR-linacs. The first approach to determine [Formula: see text] is based on a cross-calibration of measurements using a conventional linac with an electromagnet and an MR-linac. The absolute influence of the magnetic field in perpendicular orientation is quantified with the help of the conventional linac and the electromagnet. The correction factors for the parallel orientation are then derived by combining these measurements with relative measurements in the MR-linac. The second technique utilizes alanine electron paramagnetic resonance dosimetry. The alanine system as well as several ionization chambers were directly calibrated with the German primary standard for absorbed dose to water. Magnetic field correction factors for the ionization chambers were determined by a cross-calibration with the alanine in an MR-linac. Important quantities like [Formula: see text] for Farmer type ionization chambers in parallel orientation and the change of the dose to water due the magnetic field [Formula: see text] have been confirmed. In addition, magnetic field correction factors have been determined for small volume ionization chambers in parallel orientation. The electromagnet-based measurements of [Formula: see text] for [Formula: see text] MR-linacs and parallel ionization chamber orientations resulted in 0.9926(22), 0.9935(31) and 0.9841(27) for the PTW 30013, the PTW 31010 and the PTW 31021, respectively. The measurements based on the second technique resulted in values for [Formula: see text] of 0.9901(72), 0.9955(72), and 0.9885(71). Both methods show excellent accuracy and reproducibility and are therefore suitable for the determination of magnetic field correction factors. Small-volume ionization chambers showed a variation in the resulting values for [Formula: see text] and should be cross-calibrated instead of using tabulated values for correction factors.

Influence of 0.35 T magnetic field on the response of EBT3 and EBT-XD radiochromic films

PURPOSE

To investigate the inconsistency of recent literature on the effect of magnetic field on the response of radiochromic films, we studied the influence of 0.35 T magnetic field on dosimetric response of EBT3 and EBT-XD Gafchromic^TM films.

METHODS

Two different models of radiochromic films, EBT3 and EBT-XD, were investigated. Pieces of films samples from two different batches for each model were irradiated at different dose levels ranging from 1 to 20 Gy using 6 MV flattening filter free (FFF) x-rays generated by a clinical MR-guided radiotherapy system (B = 0.35 T). Film samples from the same batch were irradiated at corresponding dose levels using 6 MV FFF beam from a conventional linac (B = 0) for comparison. The net optical density was measured 48 h postirradiation using a flatbed scanner. The absorbance spectra were also measured over 500-700 nm wavelength range using a fiber-coupled spectrometer with 2.5 nm resolution. To study the effect of fractionated dose delivery to EBT3 (/EBT-XD) films, 8 (/16) Gy dose was delivered in four 2 (/4) Gy fractions with 24 h interval between fractions.

RESULTS

No significant difference was found in the net optical density and net absorbance of the films irradiated with or without the presence of magnetic field. No dependency on the orientation of the film during irradiation with respect to the magnetic field was observed. The fractionated dose delivery resulted in the same optical density as delivering the whole dose in a single fraction.

CONCLUSIONS

The 0.35 T magnetic field employed in the ViewRay^® MR-guided radiotherapy system did not show any significant influence on the response of EBT3 and EBT-XD Gafchromic^TM films.

Experimental verification the electron return effect around spherical air cavities for the MR-Linac using Monte Carlo calculation

PURPOSE

Dose deposition around unplanned air cavities during magnetic resonance-guided radiotherapy (MRgRT) is influenced by the electron return effect (ERE). This is clinically relevant for gas forming close to or inside organs at risk (OARs) that lie in the path of a single beam, for example, intestinal track during pelvic treatment. This work aims to verify Monte Carlo calculations that predict the dosimetric effects of ERE around air cavities. For this, we use GafChromic EBT3 film inside poly-methyl methacrylate (PMMA) -air phantoms.

METHOD

Four PMMA phantoms were produced. Three of the phantoms contained centrally located spherical air cavities (0.5, 3.5, 7.5 cm diameter), and one phantom contained no air. The phantoms were split to sandwich GafChromic EBT3 film in the center. The phantoms were irradiated on an Elekta Unity system using a single 10 × 10 cm2 7-MV photon beam under the influence of a 1.5-T transverse magnetic field. The measurements were replicated using the Elekta Monaco treatment planning system (TPS). Gamma analysis with pass criteria 3%/3 mm was used to compare the measured and calculated dose distributions. We also consider 3%/2 mm, 2%/3 mm, and 2%/2 mm pass criteria for interest.

RESULTS

The gamma analysis showed that >95% of the points agreed between the TPS-calculated and measured dose distributions, using 3%/3 mm criteria. The phantom containing the largest air cavity had the lowest agreement, with most of the disagreeing points lying inside the air cavity (dose to air region).

CONCLUSIONS

The dose effects due to ERE around spherical air cavities are being calculated in the TPS with sufficient accuracy for clinical use.

MRI-LINAC beam profile measurements using a plastic scintillation dosimeter

Plastic scintillation dosimeters (PSDs) possess many desirable qualities for dosimetry with LINACs. These qualities are expected to make PSDs effective for MRI-LINAC dosimetry, however little research has been conducted investigating their dosimetric performance with MRI-LINACs. In this work, an in-house PSD was used to measure 8 beam profiles with an in-line MRI-LINAC, compared with film measurements. One dimensional global gamma indices (γ) and corresponding γ pass rates were calculated to compare PSD and film profiles for the 1%/1 mm, 2%/2 mm and 3%/3 mm criterion. The mean global pass rates were 85.8%, 97.5% and 99.4% for the 1%/1 mm, 2%/2 mm and 3%/3 mm criteria, respectively. The majority of the γ failures occurred in the penumbral regions. Penumbra widths were measured to be slightly narrower with the PSD compared to film, however, the uncertainties in the measured penumbra widths brought the PSD and film penumbra widths into agreement. Differences in dose were calculated between the PSD and film, and remained within 2.2% global agreement for the central regions and 1.5% global agreement for out of field regions. These values for range of agreement were similar to the those reported in the literature for other dosimeters which are trusted for relative MRI-LINAC dosimetry.

Multichannel Film Dosimetry for Quality Assurance of Intensity Modulated Radiotherapy Treatment Plans Under 0.35 T Magnetic Field

Purpose

To evaluate the intensity modulated radiotherapy (IMRT) quality assurance (QA) results of the multichannel film dosimetry analysis with single scan method by using Gafchromic™ EBT3 (Ashland Inc., Covington, KY, USA) film under 0.35 T magnetic field.

Methods

Between September 2018 and June 2019, 70 patients were treated with ViewRay MRIdian^® (ViewRay Inc., Mountain View, CA) linear accelerator (Linac). Film dosimetry QA plans were generated for all IMRT treatments. Multichannel film dosimetry for red, green and blue (RGB) channels were compared with treatment planning system (TPS) dose maps by gamma evaluation analysis.

Results

The mean gamma passing rates of RGB channels are 97.3% ± 2.26%, 96.0% ± 3.27% and 96.2% ± 3.14% for gamma evaluation with 2% DD/2 mm distance to agreement (DTA), respectively. Moreover, the mean gamma passing rates of RGB channels are 99.7% ± 0.41%, 99.6% ± 0.59% and 99.5% ± 0.67% for gamma evaluation with 3% DD/3 mm DTA, respectively.

Conclusion

The patient specific QA using Gafchromic™ EBT3 film with multichannel film dosimetry seems to be a suitable tool to implement for MR-guided IMRT treatments under 0.35 T magnetic field. Multichannel film dosimetry with Gafchromic™ EBT3 is a consistent QA tool for gamma evaluation of the treatment plans even with 2% DD/2 mm DTA under 0.35 T magnetic field presence.

Optical imaging method to quantify spatial dose variation due to the electron return effect in an MR-linac

PURPOSE

Treatment planning systems (TPSs) for MR-linacs must employ Monte Carlo-based simulations of dose deposition to model the effects of the primary magnetic field on dose. However, the accuracy of these simulations, especially for areas of tissue-air interfaces where the electron return effect (ERE) is expected, is difficult to validate due to physical constraints and magnetic field compatibility of available detectors. This study employs a novel dosimetric method based on remotely captured, real-time optical Cherenkov and scintillation imaging to visualize and quantify the ERE.

METHODS

An intensified CMOS camera was used to image two phantoms with designed ERE cavities. Phantom A was a 40 cm × 10 cm × 10 cm clear acrylic block drilled with five holes of increasing diameters (0.5, 1, 2, 3, 4 cm). Phantom B was a clear acrylic block (25 cm × 20 cm × 5 cm) with three cavities of increasing diameter (3, 2, 1 cm) split into two halves in the transverse plane to accommodate radiochromic film. Both phantoms were imaged while being irradiated by 6 MV flattening filter free (FFF) beams within a MRIdian Viewray (Viewray, Cleveland, OH) MR-linac (0.34 T primary field). Phantom A was imaged while being irradiated by 6 MV FFF beams on a conventional linac (TrueBeam, Varian Medical Systems, San Jose, CA) to serve as a control. Images were post processed in Matlab (Mathworks Inc., Natick, MA) and compared to TPS dose volumes.

RESULTS

Control imaging of Phantom A without the presence of a magnetic field supports the validity of the optical image data to a depth of 6 cm. In the presence of the magnetic field, the optical data shows deviations from the commissioned TPS dose in both intensity and localization. The largest air cavity examined (3 cm) indicated the largest dose differences, which were above 20% at some locations. Experiments with Phantom B illustrated similar agreement between optical and film dosimetry comparisons with TPS data in areas not affected by ERE.

CONCLUSION

There are some appreciable differences in dose intensity and spatial dose distribution observed between the novel experimental data set and the dose models produced by the current clinically implemented MR-IGRT TPS.

Dosimetric Optimization and Commissioning of a High Field Inline MRI-Linac

Purpose: Unique characteristics of MRI-linac systems and mutual interactions between their components pose specific challenges for their commissioning and quality assurance. The Australian MRI-linac is a prototype system which explores the inline orientation, with radiation beam parallel to the main magnetic field. The aim of this work was to commission the radiation-related aspects of this system for its application in clinical treatments.

Methods: Physical alignment of the radiation beam to the magnetic field was fine-tuned and magnetic shielding of the radiation head was designed to achieve optimal beam characteristics. These steps were guided by investigative measurements of the beam properties. Subsequently, machine performance was benchmarked against the requirements of the IEC60976/77 standards. Finally, the geometric and dosimetric data was acquired, following the AAPM Task Group 106 recommendations, to characterize the beam for modeling in the treatment planning system and with Monte Carlo simulations. The magnetic field effects on the dose deposition and on the detector response have been taken into account and issues specific to the inline design have been highlighted.

Results: Alignment of the radiation beam axis and the imaging isocentre within 2 mm tolerance was obtained. The system was commissioned at two source-to-isocentre distances (SIDs): 2.4 and 1.8 m. Reproducibility and proportionality of the dose monitoring system met IEC criteria at the larger SID but slightly exceeded it at the shorter SID. Profile symmetry remained under 103% for the fields up to ~34 × 34 and 21 × 21 cm² at the larger and shorter SID, respectively. No penumbra asymmetry, characteristic for transverse systems, was observed. The electron focusing effect, which results in high entrance doses on central axis, was quantified and methods to minimize it have been investigated.

Conclusion: Methods were developed and employed to investigate and quantify the dosimetric properties of an inline MRI-Linac system. The Australian MRI-linac system has been fine-tuned in terms of beam properties and commissioned, constituting a key step toward the application of inline MRI-linacs for patient treatments.

High resolution silicon array detector implementation in an inline MRI-linac

PURPOSE

Dynamic dosimaging is a concept whereby a detector in motion is tracked with magnetic resonance imaging (MRI) to validate the amount and position of dose in a radiation therapy treatment on an MRI-linac. This work takes steps toward the realization of dynamic dosimaging with the novel high resolution silicon array detector: MagicPlate-512 (M512). The performance of the M512 was assessed in a 1.0 T inline MRI-linac, without simultaneous imaging and then during an imaging sequence, both during dosimetry. MR images were acquired to determine the effect of the detector and its components on image quality.

METHODS

Beam profiles were measured using the M512 on the Australian MRI-Linac and a comparison made with Gafchromic EBT3 film to investigate any intrinsic magnetic field effects in the silicon. The M512 has 512 sensitive volumes, each 0.5 × 0.5 × 0.037 mm³ in dimension, organized in a two-dimensional array. Small field sizes up to 4.2 × 3.8 cm² were investigated in both solid water and then solid lung phantoms. Beam profiles taken at 1.0 T were compared to 0 T conditions, and also to profiles taken during a gradient echo (GRE) imaging sequence. Differences in 80%-20% penumbral width and full width at half maximum (FWHM) were investigated. Localizer MR images were acquired of the detector adjacent to a water phantom.

RESULTS

Good agreement was observed between the M512 and film, with average differences in penumbral width and FWHM of <1 mm in the absence of the imaging sequence. Concurrent imaging widened the penumbra by up to 1.2 mm due to RF noise affecting the detector; film profiles were unchanged. Magnetic resonance images were affected by noise, in particular, due to the large amount of aluminum present, as well as from the USB cable, which acted as an antenna. Unfortunately, due to these issues, suitable dynamic dose imaging was not achieved with the current M512/phantom configuration and the MRI-linac. However, progress was made toward achieving this goal for future work.

CONCLUSIONS

The M512 silicon array detector successfully measured high-resolution beam profiles in agreement with Gafchromic film to within an average of <1 mm on the first MRI-linac in Australia. More effective noise reduction will be required for the achievement of dynamic dosimaging in the future.

Experimental characterization of magnetically focused electron contamination at the surface of a high-field inline MRI-linac

PURPOSE

The fringe field of the Australian MRI-linac causes contaminant electrons to be focused along the central axis resulting in a high surface dose. This work aims to characterize this effect using Gafchromic film and high-resolution detectors, MOSkin^TM and microDiamond. The secondary aim is to investigate the influence of the inline magnetic field on the relative dose response of these detectors.

METHODS

The Australian MRI-linac has the unique feature that the linac is mounted on rails allowing for measurements to be performed at different magnetic field strengths while maintaining a constant source-to-surface distance (SSD). Percentage depth doses (PDD) were collected at SSD 1.82 m in a solid water phantom positioned in a low magnetic field region and then at isocenter of the MRI where the magnetic field is 1 T. Measurements for a range of field sizes were taken with the MOSkin^TM , microDiamond, and Gafchromic^® EBT3 film. The detectors' relative responses at 1 T were compared to the near 0 T PDD beyond the region of electron contamination, that is, 20 mm depth. The near surface measurements inside the MRI bore were compared among the different detectors.

RESULTS

Skin dose in the MRI, as measured with the MOSkin^TM , was 104.5% for 2.1 × 1.9 cm² , 185.6% for 6.1 × 5.8 cm² , 369.1% for 11.8 × 11.5 cm² , and 711.1% for 23.5 × 23 cm² . The detector measurements beyond the electron contamination region showed agreement between the relative response at 1 T and near 0 T. Film was in agreement with both detectors in this region further demonstrating their relative response is unaffected by the magnetic field.

CONCLUSIONS

Experimental characterization of the high electron contamination at the surface was performed for a range of field sizes. The relative response of MOSkin^TM and microDiamond detectors, beyond the electron contamination region, were confirmed to be unaffected by the 1-T inline magnetic field.

A feasibility study for high-resolution silicon array detector performance in the magnetic field of a permanent magnet system

PURPOSE

Magnetic field effects on dose distribution and detector functionality must be well understood. The detector utilized to investigate these magnetic field effects was the DUO silicon array detector; the performance of this high spatial resolution detector was assessed under these conditions. The results were compared to Gafchromic EBT3 film to highlight any intrinsic magnetic field effects in the silicon. The results were also compared to previously published MagicPlate-512 (M512) data. The DUO has an improved spatial resolution (200 µm) over the M512 (2 mm).

METHODS

A permanent magnet named Magnetic Apparatus for RaDiation Oncology Studies (MARDOS) paired with a standard linear accelerator (linac) enables either transverse (1.2 T) or inline (0.95 T) orientations of the magnetic field with respect to the radiation beam. A 6 MV Varian 2100C Linac provided the radiation component for the measurements. The DUO detector has 505 sensitive volumes (each volume measuring 800 × 40 × 100 µm³ ) organized in two orthogonal, linear arrays. The DUO was embedded in a solid water phantom in the first set-up and then a solid lung phantom in the second set-up and placed between the magnet cones. Beam profiles were compared under the magnetic field conditions and 0 T. Small field sizes from 0.8 × 0.8 cm² up to 2.3 × 2.3 cm² were investigated. The size of the air gap above the sensitive volumes of the DUO was investigated in the transverse orientation to assess the anticipated magnetic field effects. Full width at half maximum (FWHM), 80-20% penumbral widths and maximum dose differences between detectors and between the presence/absence of a magnetic field were investigated. Symmetry was also assessed for investigation of profile skewness under the transverse field.

RESULTS

The penumbral widths measured by the DUO detector demonstrated good agreement with film and the M512 to within an average of 0.5 mm (within uncertainty: ±1 mm). The static inline magnetic field had minimal effect on the profiles in solid water. As expected, the lower density of solid lung meant that this material was more susceptible to demonstrating magnetic field effects in the dose deposited. The greatest penumbral narrowing due to the inline field (0.7 mm) occurred in lung. Central axis dose increase was greatest in lung (maximum: 9%). The transverse field widened penumbra, most notably in the solid lung phantom, by a maximum of 2.3 mm. The largest asymmetry due to the transverse field (4.6%) was also in solid lung. When the air gap above the DUO was filled with bolus, the dose maximum measured by the DUO was within 1.4% of film.

CONCLUSIONS

The DUO detector has been shown to be successful in accurately describing the dose changes for small field sizes to within a 200-µm resolution in an environment resembling that of an MRI-linac. The DUO measurements were in agreement with both film and the M512 measurements, and therefore the DUO was found to be an appropriate alternative to the M512, with improvement in terms of its higher spatial resolution. MARDOS provided a suitable environment for these preliminary tests before progressing to the MRI-linac.

Characterization of EBT3 radiochromic films for dosimetry of proton beams in the presence of magnetic fields

PURPOSE

Radiochromic film dosimetry is extensively used for quality assurance in photon and proton beam therapy. So far, GafchromicTM EBT3 film appears as a strong candidate to be used in future magnetic resonance (MR) based therapy systems. The response of Gafchromic EBT3 films in the presence of magnetic fields has already been addressed for different MR-linacs systems. However, a detailed evaluation of the influence of external magnetic fields on the film response and calibration curves for proton therapy has not yet been reported. This study aims to determine the dose responses of EBT3 films for clinical proton beams exposed to magnetic field strengths up to 1 T in order to investigate the feasibility of EBT3 film as an accurate dosimetric tool for a future MR particle therapy system (MRPT).

METHODS

The dosimetric characteristics of EBT3 films were studied for a proton beam passing through magnetic field strengths of B = 0, 0.5, and 1 T. Absorbed dose calibration and measurements were performed using clinical proton beams in the nominal energy range of 62.4-252.6 MeV. Irradiations were done using an in-house developed PMMA slab phantom placed in the center of a dipole research magnet. Monte Carlo (MC) simulations using the GATE/Geant4 toolkit were performed to predict the effect of magnetic fields on the energy deposited by proton beams in the phantom. Planned and measured doses from 3D box cube irradiations were compared to assess the accuracy of the dosimetric method using EBT3 films with/without the external magnetic field.

RESULTS

Neither for the mean pixel value nor for the net optical density, any significant deviations were observed due to the presence of an external magnetic field (B ≤ 1T) for doses up to 10 Gy. Dose-response curves for the red channel were fitted by a three-parameter function for the field-free case and for B = 1T, showing for both cases an R-square coefficient of unity and almost identical fitting parameters. Independently of the magnetic field, EBT3 films showed an under-response as high as 8% in the Bragg peak region, similarly to previously reported effects for particle therapy. No noticeable influence of the magnetic field strength was observed on the quenching effect of the EBT3 films.

CONCLUSIONS

For the first time detailed absorbed dose calibrations of EBT3 films for proton beams in magnetic field regions were performed. Results showed that EBT3 films represent an attractive solution for the dosimetry of a future MRPT system. As film response functions for protons are not affected by the magnetic field strenght, they can be used for further investigations to evaluate the dosimetric effects induced due to particle beams bending in magnetic fields regions.

Investigation of TLD and EBT3 performance under the presence of 1.5T, 0.35T, and 0T magnetic field strengths in MR/CT visible materials

PURPOSE

The aim of this study was to investigate thermoluminescent dosimeters (TLD) and radiochromic EBT3 film inside MR/CT visible geometric head and thorax phantoms in the presence of: 0, 0.35, and 1.5 T magnetic fields.

METHODS

Thermoluminescent Dosimeters reproducibility studies were examined by irradiating IROC-Houston's TLD acrylic block five times under 0 and 1.5 T configurations of Elekta's Unity system and three times under 0 and 0.35 T configurations of ViewRay's MRIdian Cobalt-60 (⁶⁰ Co) system. Both systems were irradiated with an equivalent 10 × 10 cm² field size, and a prescribed dose of 3 Gy to the maximum depth deposition (dmax). EBT3 film and TLDs were investigated using two geometrical Magnetic Resonance (MR)-guided Radiation Therapy (MRgRT) head and thorax phantoms. Each geometrical phantom had eight quadrants that combined to create a centrally located rectangular tumor (3 × 3 × 5 cm³ ) surrounded by tissue to form a 15 × 15 × 15 cm³ cubic phantom. Liquid polyvinyl chloride plastic and Superflab were used to simulate the tumor and surrounding tissue in the head phantom, respectively. Synthetic ballistic gel and a heterogeneous in-house mixture were used to construct the tumor and surrounding tissue in the thorax phantom, respectively. EBT3 and double-loaded TLDs were used in the phantoms to compare beam profiles and point dose measurements with and without magnetic fields. GEANT4 Monte Carlo simulations were performed to validate the detectors for both Unity 0 T/1.5 T and MRIdian 0 T/0.35 T configurations.

RESULTS

Average TLD block measurements which, compared the magnetic field effects (magnetic field vs 0 T) on the Unity and MRIdian systems, were 0.5% and 0.6%, respectively. The average ratios between magnetic field effects for the geometric thorax and head phantoms under the Unity system were -0.2% and 1.6% and for the MRIdian system were 0.2% and -0.3%, respectively. Beam profiles generated with both systems agreed with Monte Carlo measurements and previous literature findings.

CONCLUSIONS

TLDs and EBT3 film dosimeters could potentially be used in MR/CT visible tissue equivalent phantoms that will experience a magnetic field environment.

Monte Carlo simulations of out-of-field skin dose due to spiralling contaminant electrons in a perpendicular magnetic field

Purpose

The purpose of this study was to evaluate the potential skin dose toxicity contribution of spiralling contaminant electrons (SCE) generated in the air in an MR‐linac with a 0.35 or 1.5 T magnetic field using the EGSnrc Monte Carlo (MC) code. Comparisons to experimental results at 1.5 T are also performed.

Methods

An Elekta generated phase space file for the Unity MR‐linac is used in conjunction with the EGSnrc enhanced electric and magnetic field transport macros to simulate surface dose profiles and depth‐dose curves in panels located 5 cm away from the beam edge and positioned either parallel or perpendicular to the magnetic field. Electrons generated in the air will spiral along the magnetic field lines, and though surface doses within the field will be reduced, the electrons can contribute to out‐of‐field surface doses.

Results

Surface dose profiles showed good agreement with experimental findings and the maximum simulated doses at surfaces perpendicular to the magnetic field were 3.77 ± 0.01% and 3.55 ± 0.01% for 1.5 and 0.35 T. These results are expressed as a percentage of the maximum dose to water delivered by the photon beam. The surface dose variations in the out‐of‐field region converge to the 0 T doses within the first 0.5 cm of material. An asymmetry in the dose distribution in surfaces positioned on either side of the photon beam and aligned parallel to the magnetic field is determined to be due to the magnetic field directing electrons deeper into, or localizing them to the surface of, the measurement panel.

Conclusions

These results confirm the SCE dose contribution in surfaces perpendicular to the magnetic field and show these doses to be of the order of a few percentage of the maximum dose to water of the beam. Good agreement in the dose profiles is seen in comparisons between the MC simulations and experimental work. The effect is apparent in 0.35 and 1.5 T magnetic fields and dissipates within the first few millimeters of material. It should be noted that only SCEs from beam anteriorly incident on the patient will influence the patient surface dose, and the use of beams incident over different angles will reduce the dose to any particular patient surface.

Monte Carlo modeling of a 60Co MRI-guided radiotherapy system on Geant4 and experimental verification of dose calculation under a magnetic field of 0.35 T

Our purpose was to establish the commissioning procedure of Monte Carlo modeling on a magnetic resonance imaging-guided radiotherapy system (MRIdian, Viewray Inc.) under a magnetic field of 0.345 T through experimental measurements. To do this, we sought (i) to assess the depth-dose and lateral profiles generated by the Geant4 using either EBT3 film or the BJR-25 data; (ii) to assess the calculation accuracy under a magnetic field of 0.345 T. The radius of the electron trajectory caused by the electron return effect (ERE) in a vacuum was obtained both by the Geant4 and the theoretical methods. The surface dose on the phantom was calculated and compared with that obtained from the film measurements. The dose distribution in a phantom having two air gaps was calculated and measured with EBT 3 film. (i) The difference of depth-dose profile generated by the Geant4 from the BJR-25 data was 0.0 ± 0.8% and 0.3 ± 1.5% for field sizes of 4.5 and 27.3 cm2, respectively. Lateral dose profiles generated by Geant4 agreed well with those generated from the EBT3 film data. (ii) The radius of the electron trajectory generated by Geant4 agreed well with the theoretical values. A maximum of ~50% reduction of the surface dose under a magnetic field of 0.345 T was observed due to elimination of the electron contamination caused by the magnetic field, as determined by both the film measurements and the Geant4. Changes in the dose distributions in the air gaps caused by the ERE were observed on the Geant4 and in the film measurements. Gamma analysis (3%/3 mm) showed a pass rate of 95.1%. Commissioning procedures for the MRI-guided radiotherapy system on the Geant4 were established, and we concluded that the Geant4 had provided high calculation accuracy under a magnetic field of 0.345 T.

Microstructure changes in radiochromic films due to magnetic field and radiation

PURPOSE

To correlate the dose response and changes in microscopic structures of the radiochromic films exposed to the clinical magnetic field in the range 1.5-3 T with standard and flattening filter-free (FFF) photon beams.

METHODS

The radiochromic film was cut into 5 × 5 cm² sheets/samples from one batch. These samples were exposed to a 1.5-T and/or 3-T B-fields from an MRI scanner using an abdominal sequence for 7 min before and after irradiation with 6 MV and/or 6 MV FFF beams. Films were placed in a reference condition at 5 cm depth in a solid water phantom and exposed up to 20 Gy. The sample orientation was maintained the same during exposure, readout, and scanning electron microscopic (SEM) analysis. The samples were scanned with an Epson Expression 11000XL in a 48-bit RGB color mode at 300 dpi with red channel. Scanned images were processed in Image J and red channel mean intensity values were recorded. The samples were then coated with 6 nm gold and imaged by SEM Teneo (5 kV, 13 pA) under 2000, 2500, and 3000 magnifications for texture analysis.

RESULTS

The changes in the microstructure of the films in magnetic fields (1.5- and 3.0-T) were dose dependent. The orientation and granular size of samples at higher doses were altered compared to the controls. Needle-shaped structures of the active layer were longer and aligned for samples exposed to higher doses and magnetic field. However, no significant changes in optical density due to the presence of a magnetic field pre/postirradiation up to 20 Gy were observed.

CONCLUSION

Fine structures of the film represent the polymerization characteristics that are affected by the radiation dose in the magnetic field. Upon exposure to radiation, diacetylene monomers undergo polymerization that forms longer chains with a temporal response. Even though this study did not notice any significant changes in optical density due to the presence of magnetic field, this should be studied in simultaneous application of the magnetic field during treatment in a dedicated MR-linac unit.

Technical Note: Experimental verification of magnetic field-induced beam deflection and Bragg peak displacement for MR-integrated proton therapy

PURPOSE

Given its sensitivity to anatomical variations, proton therapy is expected to benefit greatly from integration with magnetic resonance imaging for online anatomy monitoring during irradiation. Such an integration raises several challenges, as both systems mutually interact. The proton beam will experience quasi-continuous energy loss and energy-dependent electromagnetic deflection at the same time, giving rise to a deflected beam trajectory and an altered dose distribution with a displaced Bragg peak. So far, these effects have only been predicted using Monte Carlo and analytical models, but no clear consensus has been reached and experimental benchmark data are lacking. We measured proton beam trajectories and Bragg peak displacement in a homogeneous phantom placed inside a magnetic field and compared them to simulations.

METHODS

Planar dose distributions of proton pencil beams (80-180 MeV) traversing the field of a 0.95 T NdFeB permanent magnet while depositing energy in a PMMA slab phantom were measured using EBT3 radiochromic films and simulated using the Geant4 toolkit. Deflected beam trajectories and the Bragg peak displacement were extracted from the measured planar dose distributions and compared against the simulations.

RESULTS

The lateral beam deflection was clearly visible on the EBT3 films and ranged from 1 to 10 mm for 80 to 180 MeV, respectively. Simulated and measured beam trajectories and Bragg peak displacement agreed within 0.8 mm for all studied proton energies.

CONCLUSIONS

These results prove that the magnetic field-induced Bragg peak displacement is both measurable and accurately predictable in a homogeneous phantom at 0.95 T, and allows Monte Carlo simulations to be used as gold standard for proton beam trajectory prediction in similar frameworks for MR-integrated proton therapy.