Reference dosimetry in magnetic fields: formalism and ionization chamber correction factors

Clinical reference dosimetry for the 0.5 T inline rotating biplanar Linac-MR

BACKGROUND

The world's first clinical 0.5 T inline rotating biplanar Linac-MR system is commissioned for clinical use. For reference dosimetry, unique features to device, including an SAD = 120 cm, bore clearance of 60 cm × 110 cm, as well as 0.5 T inline magnetic field, provide some challenges to applying a standard dosimetry protocol (i.e., TG-51).

PURPOSE

In this work, we propose a simple and practical clinical reference dosimetry protocol for the 0.5T biplanar Linac-MR and validated its results.

METHODS

Our dosimetry protocol for this system is as follows: tissue phantom ratios at 20 and 10 cm are first measured and converted into %dd10x beam quality specifier using equations provided and Kalach and Rogers. The converted %dd10x is used to determine the ion chamber correction factor, using the equations in the TG-51 addendum for the Exradin A12 farmer chamber used, which is cross-calibrated with one calibrated at a standards laboratory. For a 0.5 T parallel field, magnetic field effect on chamber response is assumed to have no effect and is not explicitly corrected for. Once the ion chamber correction factor for a non-standard SAD (kQ,msr) is determined, TG-51 is performed to obtain dose at a depth of 10 cm at SAD = 120 cm. The dosimetry protocol is repeated with the magnetic field ramped down. To validate our dosimetry protocol, Monte Carlo (EGSnrc) simulations are performed to confirm the determined kQ,msr values. MC Simulations and magnetic Field On versus Field Off measurements are performed to confirm that the magnetic field has no effect. To validate our overall dosimetry protocol, external dose audits, based on optical simulated luminescent dosimeters, thermal luminescent dosimeters, and alanine dosimeters are performed on the 0.5 T Linac-MR system.

RESULTS

Our EGSnrc results confirm our protocol-determined kQ,msr values, as well as our assumptions about magnetic field effects (kB = 1) within statistical uncertainty for the A-12 chamber. Our external dosimetry procedures also validated our overall dosimetry protocol for the 0.5 T biplanar Linac-MR hybrid. Ramping down the magnetic field has resulted in a dosimetric difference of 0.1%, well within experimental uncertainty.

CONCLUSION

With the 0.5 T parallel magnetic field having minimal effect on the ion chamber response, a TPR20,10 approach to determine beam quality provides an accurate method to perform clinical dosimetry for the 0.5 T biplanar Linac-MR.

Experimental determination of magnetic field quality conversion factors for eleven ionization chambers in 1.5 T and 0.35 T MR-linac systems

BACKGROUND

The static magnetic field present in magnetic resonance (MR)-guided radiotherapy systems can influence dose deposition and charged particle collection in air-filled ionization chambers. Thus, accurately quantifying the effect of the magnetic field on ionization chamber response is critical for output calibration. Formalisms for reference dosimetry in a magnetic field have been proposed, whereby a magnetic field quality conversion factor kB,Q is defined to account for the combined effects of the magnetic field on the radiation detector. Determination of kB,Q in the literature has focused on Monte Carlo simulation studies, with experimental validation limited to only a few ionization chamber models.

PURPOSE

The purpose of this study is to experimentally measure kB,Q for 11 ionization chamber models in two commercially available MR-guided radiotherapy systems: Elekta Unity and ViewRay MRIdian.

METHODS

Eleven ionization chamber models were characterized in this study: Exradin A12, A12S, A28, and A26, PTW T31010, T31021, and T31022, and IBA FC23-C, CC25, CC13, and CC08. The experimental method to measure kB,Q utilized cross-calibration against a reference Exradin A1SL chamber. Absorbed dose to water was measured for the reference A1SL chamber positioned parallel to the magnetic field with its centroid placed at the machine isocenter at a depth of 10 cm in water for a 10 × 10 cm2 field size at that depth. Output was subsequently measured with the test chamber at the same point of measurement. kB,Q for the test chamber was computed as the ratio of reference dose to test chamber output, with this procedure repeated for each chamber in each MR-guided radiotherapy system. For the high-field 1.5 T Elekta Unity system, the dependence of kB,Q on the chamber orientation relative to the magnetic field was quantified by rotating the chamber about the machine isocenter.

RESULTS

Measured kB,Q values for our test dataset of ionization chamber models ranged from 0.991 to 1.002, and 0.995 to 1.004 for the Elekta Unity and ViewRay MRIdian, respectively, with kB,Q tending to increase as the chamber sensitive volume increased. Measured kB,Q values largely agreed within uncertainty to published Monte Carlo simulation data and available experimental data. kB,Q deviation from unity was minimized for ionization chamber orientation parallel or antiparallel to the magnetic field, with increased deviations observed at perpendicular orientations. Overall (k = 1) uncertainty in the experimental determination of the magnetic field quality conversion factor, kB,Q was 0.71% and 0.72% for the Elekta Unity and ViewRay MRIdian systems, respectively.

CONCLUSIONS

For a high-field MR-linac, the characterization of ionization chamber performance as angular orientation varied relative to the magnetic field confirmed that the ideal orientation for output calibration is parallel. For most of these chamber models, this study represents the first experimental characterization of chamber performance in clinical MR-linac beams. This is a critical step toward accurate output calibration for MR-guided radiotherapy systems and the measured kB,Q values will be an important reference data source for forthcoming MR-linac reference dosimetry protocols.

MR compatible detectors assessment for a 0.35 T MR-linac commissioning

PURPOSE

To assess a large panel of MR compatible detectors on the full range of measurements required for a 0.35 T MR-linac commissioning by using a specific statistical method represented as a continuum of comparison with the Monte Carlo (MC) TPS calculations. This study also describes the commissioning tests and the secondary MC dose calculation validation.

MATERIAL AND METHODS

Plans were created on the Viewray TPS to generate MC reference data. Absolute dose points, PDD, profiles and output factors were extracted and compared to measurements performed with ten different detectors: PTW 31010, 31021, 31022, Markus 34045 and Exradin A28 MR ionization chambers, SN Edge shielded diode, PTW 60019 microdiamond, PTW 60023 unshielded diode, EBT3 radiochromic films and LiF µcubes. Three commissioning steps consisted in comparison between calculated and measured dose: the beam model validation, the output calibration verification in four different phantoms and the commissioning tests recommended by the IAEA-TECDOC-1583.

MAIN RESULTS

The symmetry for the high resolution detectors was higher than the TPS data of about 1%. The angular responses of the PTW 60023 and the SN Edge were - 6.6 and - 11.9% compared to the PTW 31010 at 60°. The X/Y-left and the Y-right penumbras measured by the high resolution detectors were in good agreement with the TPS values except for the PTW 60023 for large field sizes. For the 0.84 × 0.83 cm2 field size, the mean deviation to the TPS of the uncorrected OF was - 1.7 ± 1.6% against - 4.0 ± 0.6% for the corrected OF whereas we found - 4.8 ± 0.8% for passive dosimeters. The mean absolute dose deviations to the TPS in different phantoms were 0 ± 0.4%, - 1.2 ± 0.6% and 0.5 ± 1.1% for the PTW 31010, PTW 31021 and Exradin A28 MR respectively.

CONCLUSIONS

The magnetic field effects on the measurements are considerably reduced at low magnetic field. The PTW 31010 ionization chamber can be used with confidence in different phantoms for commissioning and QA tests requiring absolute dose verifications. For relative measurements, the PTW 60019 presented the best agreement for the full range of field size. For the profile assessment, shielded diodes had a behaviour similar to the PTW 60019 and 60023 while the ionization chambers were the most suitable detectors for the symmetry. The output correction factors published by the IAEA TRS 483 seem to be applicable at low magnetic field pending the publication of new MR specific values.

Beam quality correction factors for ionization chambers in a 0.35 T magnetic resonance (MR)-linac - A Monte Carlo study

PURPOSE

The purpose of this study was to directly calculate [Formula: see text] correction factors for four cylindrical ICs for a 0.35 T MR-linac using the Monte Carlo (MC) method.

METHODS

A previously-validated TOPAS/GEANT4 MC head model of the 0.35 T MR-linac was employed. The MR-compatible Exradin A12, A1SL, A26, and A28 cylindrical ICs were modeled considering the dead volume in the air cavity. The [Formula: see text] correction factor was determined for initial electron energies of 5-7 MeV. The correction factor was calculated for all four angular orientations in the lateral plane. The impact of the 0.35 T magnetic field on the IC response was also investigated.

RESULTS

The maximum beam quality dependence in the [Formula: see text] exhibited by the A12, A1SL, A26, and A28 ICs was 1.10 %, 2.17 %, 0.81 %, and 1.75 %, respectively, considering all angular orientations. The magnetic field dependence was < 1 % and the maximum [Formula: see text] correction was < 2 % when the detector was aligned along the direction of the magnetic field at 0° and 180° angles. The A12 IC over-responded up to 5.40 % for the orthogonal orientation. An asymmetry in the response of up to 8.30 % was noted for the A28 IC aligned at 90° and 270° angles.

CONCLUSIONS

A parallel orientation for the IC, with respect to the magnetic field, is recommended for reference dosimetry in MRgRT. Both over and under-response in the IC signal was noted for the orthogonal orientations, which is highly dependent on the cavity diameter, cavity length, and the dead volume.

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.

Magnetic field influence on the light yield from fiber-coupled BCF-60 plastic scintillators of relevance for output factor dosimetry in MR-linacs

Organic plastic scintillators are of interest for ionizing radiation dosimetry in megavoltage photon beams because plastic scintillators have a mass density very similar to that of water. This leads to insignificant perturbation of the electron fluence at the point of measurement in a water phantom. This feature is a benefit for dosimetry in strong magnetic fields (e.g., 1.5 T) as found in linacs with magnetic resonance imaging. The objective of this work was to quantify if the light yield per dose for the scintillating fiber BCF-60 material from Saint-Gobain Ceramics and Plastics Inc. is constant regardless of the magnetic flux density. This question is of importance for establishing traceable measurement in MR linacs using this detector type. Experiments were carried out using an accelerator combined with an electromagnet (max 0.7 T). Scintillator probes were read out using chromatic stem-removal techniques based on two optical channels or full spectral information. Reference dosimetry was carried out with PTW31010 and PTW31021 ionization chambers. TOPAS/GEANT4 was used for modelling. The light yield per dose for the BCF-60 was found to be strongly influenced by the magnitude of the magnetic field from about 1 mT to 0.7 T. The light yield per dose increased (1.3 ± 0.2)% (k = 1) from 1 mT to 10 mT and it increased (4.5 ± 0.9)% (k = 1) from 0 T to 0.7 T. Previous studies of the influence of magnetic fields on medical scintillator dosimetry have been unable to clearly identify if observed changes in scintillator response with magnetic field strength were related to changes in dose, stem signal removal, or scintillator light yield. In the current study of BCF-60, we see a clear change in light yield with magnetic field, and none of the other effects.

The influence of different versions of FLUKA and GEANT4 on the calculation of response functions of ionization chambers in clinical proton beams

Objective.To investigate the influence of different versions of the Monte Carlo codesgeant4 andflukaon the calculation of overall response functionsfQof air-filled ionization chambers in clinical proton beams.Approach. fQfactors were calculated for six plane-parallel and four cylindrical ionization chambers withgeant4 andfluka. These factors were compared to already published values that were derived using older versions of these codes.Main results.Differences infQfactors calculated with different versions of the same Monte Carlo code can be up to ∼1%. Especially forgeant4, the updated version leads to a more pronounced dependence offQon proton energy and to smallerfQfactors for high energies.Significance.Different versions of the same Monte Carlo code can lead to differences in the calculation offQfactors of up to ∼1% without changing the simulation setup, transport parameters, ionization chamber geometry modeling, or employed physics lists. These findings support the statement that the dominant contributor to the overall uncertainty of Monte Carlo calculatedfQfactors are type-B uncertainties.

Dose-response dependencies of OSL dosimeters in conventional linacs and 1.5T MR-linacs: an experimental and Monte Carlo study

Objective. This work focuses on the optically stimulated luminescence dosimetry (OSLD) dose-response characterization, with emphasis on 1.5T MR-Linacs.Approach. Throughout this study, the nanoDots OSLDs (Landauer, USA) were considered. In groups of three, the mean OSLD response was measured in a conventional linac and an MR-Linac under various irradiation conditions to investigate (i) dose-response linearity with and without the 1.5T magnetic field, (ii) signal fading rate and its dependencies, (iii) beam quality, detector orientation and dose rate dependencies in a conventional linac, (iii) potential MR imaging related effects on OSLD response and (iv) detector orientation dependence in an MR-Linac. Monte Carlo calculations were performed to further quantify angular dependence after rotating the detector around its central axis parallel to the magnetic field, and determine the magnetic field correction factors,kB,Q,for all cardinal detector orientations.Main results. OSLD dose-response supralinearity in an MR-Linac setting was found to agree within uncertainties with the corresponding one in a conventional linac, for the axial detector orientation investigated. Signal fading rate does not depend on irradiation conditions for the range of 3-30 d considered. OSLD angular (orientation) dependence is more pronounced under the presence of a magnetic field. OSLDs irradiated with and without real-time T2w MR imaging enabled during irradiation yielded the same response within uncertainties.kB,Qvalues were determined for all three cardinal orientations. Corrections needed reached up to 6.4%. However, if OSLDs are calibrated in the axial orientation and then irradiated in an MR-Linac placed again in the axial orientation (perpendicular to the magnetic field), then simulations suggest thatkB,Qcan be considered unity within uncertainties, irrespective of the incident beam angle.Significance. This work contributes towards OSLD dose-response characterization and relevant correction factors availability. OSLDs are suitable for QA checks in MR-based beam gating applications andin vivodosimetry in MR-Linacs.

Quantifying uncertainties associated with reference dosimetry in an MR-Linac

BACKGROUND

Magnetic resonance (MR)-guided radiation therapy provides capabilities to utilize high-resolution and real-time MR imaging before and during treatment, which is critical for adaptive radiotherapy. This emerging modality has been promptly adopted in the clinic settings in advance of adaptations to reference dosimetry formalism that are needed to account for the presence of strong magnetic fields. In particular, the influence of magnetic field on the uncertainty of parameters in the reference dosimetry equation needs to be determined in order to fully characterize the uncertainty budget for reference dosimetry in MR-guided radiation therapy systems.

PURPOSE

To identify and quantify key sources of uncertainty in the reference dosimetry of external high energy radiotherapy beams in the presence of a strong magnetic field.

METHODS

In the absence of a formalized Task Group report for reference dosimetry in MR-integrated linacs, the currently suggested formalism follows the TG-51 protocol with the addition of a quality conversion factor kBQ accounting for the effects of the magnetic field on ionization chamber response. In this work, we quantify various sources of uncertainty that impact each of the parameters in the formalism, and evaluate their overall contribution to the final dose. Measurements are done in a 1.5 T MR-Linac (Unity, Elekta AB, Stockholm, Sweden) which integrates a 1.5 T Philips MR scanner and a 7 MVFFF linac. The responses of several reference-class small volume ionization chambers (Exradin:A1SL, IBA:CC13, PTW:Semiflex-3D) and Farmer type ionization chambers (Exradin:A19, IBA:FC65-G) were evaluated throughout this process. Long-term reproducibility and stability of beam quality,TPR 10 20 ${\mathrm{TPR}}_{10}^{20}$ , was also measured with an in-house built phantom.

RESULTS

Relative to the conventional external high energy linacs, the uncertainty on overall reference dose in MR-linac is more significantly affected by the chamber setup: A translational displacement along y-axis of ± 3 mm results in dose variation of < |0.20| ± 0.02% (k = 1), while rotation of ± 5° in horizontal and vertical parallel planes relative to relative to the direction of magnetic field, did not exceed variation of < |0.44| ± 0.02% for all 5 ionization chambers. We measured a larger dose variation for xy-plane (horizontal) rotations (< |0.44| ± 0.02% (k = 1)) than for yz-plane (vertical) rotations (< ||0.28| ± 0.02% (k = 1)), which we associate with the gradient of kB,Q as a function of chamber orientation with respect to direction of the B0 -field. Uncertainty in Pion (for two depths), Ppol (with various sub-studies including effects of cable length, cable looping in the MRgRT bore, connector type in magnetic environment), and Prp were determined. Combined conversion factor kQ × kB,Q was provided for two reference depths at four cardinal angle orientations. Over a two-year period, beam quality was quite stable withTPR 10 20 ${\mathrm{TPR}}_{10}^{20}$ being 0.669 ± 0.01%. The actual magnitude ofTPR 10 20 ${\mathrm{TPR}}_{10}^{20}$ was measured using identical equipment and compared between two different Elekta Unity MR-Linacs with results agreeing to within 0.21%.

CONCLUSION

In this work, the uncertainty of a number of parameters influencing reference dosimetry was quantified. The results of this work can be used to identify best practice guidelines for reference dosimetry in the presence of magnetic fields, and to evaluate an uncertainty budget for future reference dosimetry protocols for MR-linac.

Investigation of Monte Carlo simulations of the electron transport in external magnetic fields using Fano cavity test

PURPOSE

Monte Carlo simulations are crucial for calculating magnetic field correction factors kB for the dosimetry in external magnetic fields. As in Monte Carlo codes the charged particle transport is performed in straight condensed history (CH) steps, the curved trajectories of these particles in the presence of external magnetic fields can only be approximated. In this study, the charged particle transport in presence of a strong magnetic field B→ was investigated using the Fano cavity test. The test was performed in an ionization chamber and a diode detector, showing how the step size restrictions must be adjusted to perform a consistent charged particle transport within all geometrical regions.

METHODS

Monte Carlo simulations of the charged particle transport in a magnetic field of 1.5 T were performed using the EGSnrc code system including an additional EMF-macro for the transport of charged particle in electro-magnetic fields. Detailed models of an ionization chamber and a diode detector were placed in a water phantom and irradiated with a so called Fano source, which is a monoenergetic, isotropic electron source, where the number of emitted particles is proportional to the local density.

RESULTS

The results of the Fano cavity test strongly depend on the energy of charged particles and the density within the given geometry. By adjusting the maximal length of the charged particle steps, it was possible to calculate the deposited dose in the investigated regions with high accuracy (<0.1%). The Fano cavity test was performed in all regions of the detailed detector models. Using the default value for the step size in the external magnetic field, the maximal deviation between Monte Carlo based and analytical dose value in the sensitive volume of the ion chamber and diode detector was 8% and 0.1%, respectively.

CONCLUSIONS

The Fano cavity test is a crucial validation method for the modeled detectors and the transport algorithms when performing Monte Carlo simulations in a strong external magnetic field. Special care should be given, when calculating dose in volumes of low density. This study has shown that the Fano cavity test is a useful method to adapt particle transport parameters for a given simulation geometry.

Monte Carlo calculated ionization chamber correction factors in clinical proton beams - deriving uncertainties from published data

For the update of the IAEA TRS-398 Code of Practice (CoP), global ionization chamber factors (fQ) and beam quality correction factors (kQ) for air-filled ionization chambers in clinical proton beams have been calculated with different Monte Carlo codes. In this study, average Monte Carlo calculated fQ and kQ factors are provided and the uncertainty of these factors is estimated. Average fQ factors in monoenergetic proton beams with energies between 60 MeV and 250 MeV were derived from Monte Carlo calculated fQ factors published in the literature. Altogether, 195 fQ factors for six plane-parallel and three cylindrical ionization chambers calculated with penh, fluka and geant4 were incorporated. Additionally, a weighted standard deviation of fQ factors was calculated, where the same weight was assigned to each Monte Carlo code. From average fQ factors, kQ factors were derived and compared to the values from the IAEA TRS-398 CoP published in 2000 as well as to the values of the upcoming version. Average Monte Carlo calculated fQ factors are constant within 0.6% over the energy range investigated. In general, the different Monte Carlo codes agree within 1% for low energies and show larger differences up to 2% for high energies. As a result, the standard deviation of fQ factors increases with energy and is ∼0.3% for low energies and ∼0.8% for high energies. kQ factors derived from average Monte Carlo calculated fQ factors differ from the values presented in the IAEA TRS-398 CoP by up to 2.4%. The overall estimated uncertainty of Monte Carlo calculated kQ factors is ∼0.5%-1% smaller than the uncertainties estimated in IAEA TRS-398 CoP since the individual ionization chamber characteristics (e.g. fluence perturbations) are considered in detail in Monte Carlo calculations. The agreement between Monte Carlo calculated kQ factors and the values of the upcoming version of IAEA TRS-398 CoP is better with deviations smaller than 1%.

Monte Carlo simulation for proton dosimetry in magnetic fields: Fano test and magnetic field correction factorskBfor Farmer-type ionization chambers

Objective. In this contribution we present a special Fano test for charged particles in presence of magnetic fields in the MC code TOol for PArticle Simulation (TOPAS), as well as the determination of magnetic field correction factorskBfor Farmer-type ionization chambers using proton beams.Approach. Customized C++ extensions for TOPAS were implemented to model the special Fano tests in presence of magnetic fields for electrons and protons. The Geant4-specific transport parameters,DRoverRandfinalRange,were investigated to optimize passing rate and computation time. ThekBwas determined for the Farmer-type PTW 30013 ionization chamber, and 5 custom built ionization chambers with same geometry but varying inner radius, testing magnetic flux density ranging from 0 to 1.0 T and two proton beam energies of 157.43 and 221.05 MeV.Main results. Using the investigated parameters, TOPAS passed the Fano test within 0.39 ± 0.15% and 0.82 ± 0.42%, respectively for electrons and protons. The chamber response (kB,M,Q) gives a maximum at different magnetic flux densities depending of the chamber size, 1.0043 at 1.0 T for the smallest chamber and 1.0051 at 0.2 T for the largest chamber. The local dose differencecBremained ≤ 0.1% for both tested energies. The magnetic field correction factorkB, for the chamber PTW 30013, varied from 0.9946 to 1.0036 for both tested energies.Significance. The developed extension for the special Fano test in TOPAS MC code with the adjusted transport parameters, can accurately transport electron and proton particles in magnetic field. This makes TOPAS a valuable tool for the determination ofkB. The ionization chambers we tested showed thatkBremains small (≤0.72%). To the best of our knowledge, this is the first calculations ofkBfor proton beams. This work represents a significant step forward in the development of MRgPT and protocols for proton dosimetry in presence of magnetic field.

Experimental and Monte Carlo-based determination of magnetic field correction factors k B , Q $k_{B,Q}$ in high-energy photon fields for two ionization chambers

BACKGROUND

The integration of magnetic resonance tomography into clinical linear accelerators provides high-contrast, real-time imaging during treatment and facilitates online-adaptive workflows in radiation therapy treatments. The associated magnetic field also bends the trajectories of charged particles via the Lorentz force, which may alter the dose distribution in a patient or a phantom and affects the dose response of dosimetry detectors.

PURPOSE

To perform an experimental and Monte Carlo-based determination of correction factors k B , Q $k_{B,Q}$ , which correct the response of ion chambers in the presence of external magnetic fields in high-energy photon fields.

METHODS

The response variation of two different types of ion chambers (Sun Nuclear SNC125c and SNC600c) in strong external magnetic fields was investigated experimentally and by Monte Carlo simulations. The experimental data were acquired at the German National Metrology Institute, PTB, using a clinical linear accelerator with a nominal photon energy of 6 MV and an external electromagnet capable of generating magnetic flux densities of up to 1.5 T in opposite directions. The Monte Carlo simulation geometries corresponded to the experimental setup and additionally to the reference conditions of IAEA TRS-398. For the latter, the Monte Carlo simulations were performed with two different photon spectra: the 6 MV spectrum of the linear accelerator used for the experimental data acquisition and a 7 MV spectrum of a commercial MRI-linear accelerator. In each simulation geometry, three different orientations of the external magnetic field, the beam direction and the chamber orientation were investigated.

RESULTS

Good agreement was achieved between Monte Carlo simulations and measurements with the SNC125c and SNC600c ionization chambers, with a mean deviation of 0.3% and 0.6%, respectively. The magnitude of the correction factor k B , Q $k_{B,Q}$ strongly depends on the chamber volume and on the orientation of the chamber axis relative to the external magnetic field and the beam directions. It is greater for the SNC600c chamber with a volume of 0.6 cm³ than for the SNC125c chamber with a volume of 0.1 cm³ . When the magnetic field direction and the chamber axis coincide, and they are perpendicular to the beam direction, the ion chambers exhibit a calculated overresponse of less than 0.7(6)% (SNC600c) and 0.3(4)% (SNC125c) at 1.5 T and less than 0.3(0)% (SNC600c) and 0.1(3)% (SNC125c) for 0.35 T for nominal beam energies of 6 MV and 7 MV. This chamber orientation should be preferred, as k B , Q $k_{B,Q}$ may increase significantly in other chamber orientations. Due to the special geometry of the guard ring, no dead-volume effects have been observed in any orientation studied. The results show an intra-type variation of 0.17% and 0.07% standard uncertainty (k=1) for the SNC125c and SNC600c, respectively.

CONCLUSION

Magnetic field correction factors k B , Q $k_{B,Q}$ for two different ion chambers and for typical clinical photon beam qualities were presented and compared with the few data existing in the literature. The correction factors may be applied in clinical reference dosimetry for existing MRI-linear accelerators.

ACPSEM position paper: dosimetry for magnetic resonance imaging linear accelerators

Consistency and clear guidelines on dosimetry are essential for accurate and precise dosimetry, to ensure the best patient outcomes and to allow direct dose comparison across different centres. Magnetic Resonance Imaging Linac (MRI-linac) systems have recently been introduced to Australasian clinics. This report provides recommendations on reference dosimetry measurements for MRI-linacs on behalf of the Australiasian College of Physical Scientists and Engineers in Medicine (ACPSEM) MRI-linac working group. There are two configurations considered for MRI-linacs, perpendicular and parallel, referring to the relative direction of the magnetic field and radiation beam, with different impacts on dose deposition in a medium. These recommendations focus on ion chambers which are most commonly used in the clinic for reference dosimetry. Water phantoms must be MR safe or conditional and practical limitations on phantom set-up must be considered. Solid phantoms are not advised for reference dosimetry. For reference dosimetry, IAEA TRS-398 recommendations cannot be followed completely due to physical differences between conventional linac and MRI-linac systems. Manufacturers' advice on reference conditions should be followed. Beam quality specification of TPR_20,10 is recommended. The configuration of the central axis of the ion chamber relative to the magnetic field and radiation beam impacts the chamber response and must be considered carefully. Recommended corrections to delivered dose are [Formula: see text], a correction for beam quality and [Formula: see text], for the impact of the magnetic field on dosimeter response in the magnetic field. Literature based values for [Formula: see text] are given. It is important to note that this is a developing field and these recommendations should be used together with a review of current literature.

Dosimetry in 1.5 T MR-Linacs: Monte Carlo determination of magnetic field correction factors and investigation of the air gap effect

BACKGROUND

In Magnetic Resonance-Linac (MR-Linac) dosimetry formalisms, a new correction factor, k_B,Q , has been introduced to account for corresponding changes to detector readings under the beam quality, Q, and the presence of magnetic field, B.

PURPOSE

This study aims to develop and implement a Monte Carlo (MC)-based framework for the determination of k_B,Q correction factors for a series of ionization chambers utilized for dosimetry protocols and dosimetric quality assurance checks in clinical 1.5 T MR-Linacs. Their dependencies on irradiation setup conditions are also investigated. Moreover, to evaluate the suitability of solid phantoms for dosimetry checks and end-to-end tests, changes to the detector readings due to the presence of small asymmetrical air gaps around the detector's tip are quantified.

METHODS

Phase space files for three irradiation fields of the ELEKTA Unity 1.5 T/7 MV flattening-filter-free MR-Linac were provided by the manufacturer and used as source models throughout this study. Twelve ionization chambers (three farmer-type and nine small-cavity detectors, from three manufacturers) were modeled (including their dead volume) using the EGSnrc MC code package. k_B,Q values were calculated for the 10 × 10 cm² irradiation field and for four cardinal orientations of the detectors' axes with respect to the 1.5 T magnetic field. Potential dependencies of k_B,Q values with respect to field size, depth, and phantom material were investigated by performing additional simulations. Changes to the detectors' readings due to the presence of small asymmetrical air gaps (0.1 up to 1 mm) around the chambers' sensitive volume in an RW3 solid phantom were quantified for three small-cavity chambers and two orientations.

RESULTS

For both parallel (to the magnetic field) orientations, k_B,Q values were found close to unity. The maximum correction needed was 1.1%. For each detector studied, the k_B,Q values calculated for the two parallel orientations agreed within uncertainties. Larger corrections (up to 5%) were calculated when the detectors were oriented perpendicularly to the magnetic field. Results were compared with corresponding ones found in the literature, wherever available. No considerable dependence of k_B,Q with respect to field size (down to 3 × 3 cm² ), depth, or phantom material was noticed, for the detectors investigated. As compared to the perpendicular one, in the parallel to the magnetic field orientation, the air gap effect is minimized but is still considerable even for the smallest air gap considered (0.1 mm).

CONCLUSION

For the 10 × 10 cm² field, magnetic field correction factors for 12 ionization chambers and four orientations were determined. For each detector, the k_B,Q value may be also applied for dosimetry procedures under different irradiation parameters provided that the orientation is taken into account. Moreover, if solid phantoms are used, even the smallest asymmetrical air gap may still bias small-cavity chamber response. This work substantially expands the availability and applicability of k_B,Q correction factors that are detector- and orientation-specific, enabling more options in MR-Linac dosimetry checks, end-to-end tests, and quality assurance protocols.

Dosimetric characterization of a novel commercial plastic scintillation detector with an MR-Linac

BACKGROUND

Plastic scintillators have been used as radiation detectors for the past few years, as they are water-equivalent and independent of the dose, dose rate, and angle of incidence. In addition, they are also independent of the presence of a magnetic field and could be used for in vivo dosimetry in an MR-Linac. With the advent of a new commercial scintillation detector, Blue Physics Model 10, its characterization has been performed on an MR-Linac with a view to future applications.

PURPOSE

To perform the dosimetric characterization and study potential applications of a novel commercial plastic scintillation detector in a MR-Linac.

METHODS

Scintillation detector description, calibration procedure, short-term repeatability, dose-response linearity, dose-rate dependence, angular dependence, and temperature dependence have been studied. Percent-depth-dose (PDD) and beam profiles were measured for small fields and a standard field, as well as output factors, for comparison with other PTW detectors: a diamond diode and PinPoint and Semiflex 3D ionization chambers. The suitability of the plastic scintillator for in vivo dosimetry in a magnetic field has also been studied measuring the dose to a point in an anthropomorphic phantom while acquiring MR imaging. This measured dose was compared with that calculated with Monaco planning system and with that measured with a PTW Semiflex 3D chamber, the latter without acquiring MR images.

RESULTS

Short-term repeatability presented negligible variations (<0.4%) for 100 and 20 MU. Similar results were obtained for dose-response linearity and dose-rate dependence. A small angular dependence was determined, while the scintillator resulted practically independent of the temperature. PDDs showed excellent agreement except in the build-up region, and calculated penumbras with the profiles given by the scintillator were between the ones obtained with the diamond detector and the PinPoint ionization chamber. Measured OF with the scintillator were the highest between all detectors, 1.26% higher than the value obtained with the microdiamond for the smallest field measured, 0.5 × 0.5 cm² . Finally, the total dose to a point measured with the scintillator was 0.51% higher compared to that calculated by the planning system.

CONCLUSION

The Blue Physics model 10 scintillation system showed excellent dosimetric characteristics. Its response independent of the temperature and the presence of a magnetic field make it suitable for in vivo dosimetry in an MR-Linac while acquiring MR images, which could solve the impossibility of performing a dosimetric QA for each adapted plan. Furthermore, its temporal resolution allows independent radiation pulses to be measured and visualized, which could be used in future applications.

Elekta Unity MR-linac commissioning: mechanical and dosimetry tests

We report the commissioning results of Elekta Unity for the dosimetric performance and mechanical quality assurance (QA), and propose additional commissioning procedures. Mechanical tests included multi-leaf collimator (MLC) positional accuracy, radiation isocenter diameter at the center and off-center position, and coincidence between the magnetic resonance (MR) image center and radiation isocenter. Comparisons between the measurements and calculations of the simple irradiated field, intensity modulated radiation therapy (IMRT) commissioning, MLC output factor ratio, validation of independent dose calculation software and end-to-end testing were performed to evaluate dosimetric performance. The average values of the MLC positional accuracy for film- and imaging device-based analysis were −0.1 and 0.3 mm, respectively. The measured radiation isocenter size was 0.41 mm, and the off-center results were within 1 mm. The coincidence was −0.21, −1.19 and 0.49 mm along the x-, y- and z-axes, respectively. The calculated percent depth doses (PDD) and profiles agreed with the measurements. The results of independent dose calculation were within the action level recommended by American Associations of Physicist in Medicine. The gamma passing rate (GPR) for IMRT commissioning was 98.6 ± 0.9%, and end-to-end testing of adapted plans showed agreement within 2% between the measurement and calculation. We reported the results of mechanical and dosimetric performances of Elekta Unity, and proposed novel commissioning procedures. Our results should provide knowledge to the physics community for enhancing the QA programs.

Practical guidelines of online MR-guided adaptive radiotherapy

The first magnetic resonance (MR)-guided radiotherapy system in Japan was installed in May 2017. Implementation of online MR-guided adaptive radiotherapy (MRgART) began in February 2018. Online MRgART offers greater treatment accuracy owing to the high soft-tissue contrast in MR-images (MRI), compared to that in X-ray imaging. The Japanese Society for Magnetic Resonance in Medicine (JSMRM), Japan Society of Medical Physics (JSMP), Japan Radiological Society (JRS), Japanese Society of Radiological Technology (JSRT), and Japanese Society for Radiation Oncology (JASTRO) jointly established the comprehensive practical guidelines for online MRgART. These guidelines propose the essential requirements for clinical implementation of online MRgART with respect to equipment, personnel, institutional environment, practice guidance, and quality assurance/quality control (QA/QC). The minimum requirements for related equipment and QA/QC tools, recommendations for safe operation of MRI system, and the implementation system are described. The accuracy of monitor chamber and detector in dose measurements should be confirmed because of the presence of magnetic field. The ionization chamber should be MR-compatible. Non-MR-compatible devices should be used in an area that is not affected by the static magnetic field (outside the five Gauss line), and their operation should be checked to ensure that they do not affect the MR image quality. Dose verification should be performed using an independent dose verification system that has been confirmed to be reliable through commissioning. This guideline proposes the checklists to ensure the safety of online MRgART. Successful clinical implementation of online MRgART requires close collaboration between physician, radiological technologist, nurse, and medical physicist.

MRI-guided Radiotherapy (MRgRT) for treatment of Oligometastases: Review of clinical applications and challenges

PURPOSE

Early clinical results on the application of magnetic resonance imaging (MRI) coupled with a linear accelerator to deliver MR-guided radiation therapy (MRgRT) have demonstrated feasibility for safe delivery of stereotactic body radiotherapy (SBRT) in treatment of oligometastatic disease. Here we set out to review the clinical evidence and challenges associated with MRgRT in this setting.

METHODS AND MATERIALS

We performed a systematic review of the literature pertaining to clinical experiences and trials on the use of MRgRT primarily for the treatment of oligometastatic cancers. We reviewed the opportunities and challenges associated with the use of MRgRT.

RESULTS

Benefits of MRgRT pertaining to superior soft-tissue contrast, real-time imaging and gating, and online adaptive radiotherapy facilitate safe and effective dose escalation to oligometastatic tumors while simultaneously sparing surrounding healthy tissues. Challenges concerning further need for clinical evidence and technical considerations related to planning, delivery, quality assurance (QA) of hypofractionated doses, and safety in the MRI environment must be considered.

CONCLUSIONS

The promising early indications of safety and effectiveness of MRgRT for SBRT-based treatment of oligometastatic disease in multiple treatment locations should lead to further clinical evidence to demonstrate the benefit of this technology.

The impact of ion chamber components onkB,Qfor reference dosimetry in MRgRT

For reference dosimetry in MRgRT,k_B,Qis used to correct for the impact of the magnetic field on the chamber calibration coefficient. It has been demonstrated that for accurate simulation ofk_B,Qthe dead volume (DV) must be considered. This work goes one step further by analysing the contribution of secondary electrons generated in the various chamber components tok_B,Q. The Farmer-type chamber PTW 30013 geometry was modelled for two different DVs. Monte Carlo simulations were performed for a⁶⁰Co source and a 7 MV MRI-linac and the model was validated against measurements. Both parallel (α = 0° or 180°) and perpendicular (α = 90° or 270°) orientations of the chamber and the magnetic (B) field were considered, and severalB-field strengths between 0 T and 1.5 T. To study the dose contribution to the reduced volume (RV = cavity - DV) from the secondary electrons produced in certain components of the chamber the labelling of the particles was implemented in the PENELOPE user code PENMAIN. A separate model with each solid component of the chamber modelled as liquid water was used to investigate the impact of material choice onk_B,Q. Results show that simulatedk_B,Qvalues agree better with the measuredk_B,Qwhen the DV is considered. It is demonstrated that small components of the chamber impactk_B,Qconsiderably, since the contribution to the RV-dose from the bodies closer to the RV is higher than withoutB. Moreover, it is seen that the impact to the dose in the RV is reduced when the material of each component is modelled as liquid water. Therefore, chamber design and, to a lesser extent, choice of material affectk_B,Q, and an accurate geometrical model of the chamber components and its further validation are important for correct calculations ofk_B,Q.

Monte Carlo optimization and experimental validation of a prototype ionization chamber for accurate magnetic resonance image guided radiation therapy (MRgRT) daily output constancy measurements in solid phantoms

PURPOSE

To optimize the design, develop and test a prototype ionization chamber for accurate daily output constancy measurements in solid phantoms in clinical MRgRT radiotherapy beams. Up to 4 % variations in response using commercial ionization chambers have been previously reported; the prototype ionization chamber developed here aims to minimize these variations.

METHODS

Monte Carlo simulations with the EGSnrc code system are used to optimize an ionization chamber design by increasing the thickness of a brass (high-density, non-ferromagnetic, easy-to-machine) wall until results consistent with no air gap are produced for simulations with a 1.5 T and 0.35 T magnetic field, with a 0.2 mm air gap and varying the placement of the chamber model within the air gap. Based on the results of these simulations, prototype ionization chambers are manufactured and tested in conventional linac beams and in a 7 MV Elekta Unity MR-linac. The chambers are rotated about their axes, both parallel and perpendicular to the 1.5 T magnetic field, through 360 degrees in a plastic phantom with measurements made at each cardinal angle. This reveals any variation in chamber response by varying the thickness of the air gap between the chamber and the phantom.

RESULTS

Monte Carlo simulations demonstrate that the optimal thickness of the chamber wall to mitigate the effect of an asymmetric air gap between the chamber and the plastic phantom is 1.1 mm of brass. With this thickness, the differences between simulations with and without an air gap and with asymmetric placement of the chamber within the air gap are less than 0.2 %. A prototype chamber constructed with a 1.1 mm brass wall thickness exhibits less than 0.3 % variation in response when rotated about its axis in the plastic phantom in a beam from an MR-linac, independent of whether its axis is parallel or perpendicular to the magnetic field.

CONCLUSION

The optimized ionization chamber design and validated prototype for accurate MR-linac daily output constancy measurements allows utilization of conventional phantoms and procedures in MRgRT systems. This can minimize disruption to clinical workflow for MR-linac QA measurements. This article is protected by copyright. All rights reserved.

Experimental characterisation of the magnetic field correction factor,kB⃗,for Roos chambers in a parallel MRI-linac

Objective.Reference dosimetry on an MRI-linac requires a chamber specific magnetic field correction factor,kB⃗.This work aims to measure the correction factor for a parallel plate chamber on a parallel MRI-linac.Approach.kB⃗is defined as the ratio of the absorbed dose to water calibration coefficient in the presence of the magnetic field,ND,wB⃗relative to that under 0 T conditions,ND,w0T.kB⃗was measured via aND,wtransfer to a field chamber at each magnetic field strength from a chamber with knownND,wandkB⃗.This was achieved on the parallel MRI-linac by moving the measurement set-up between a high magnetic field strength region at the MRI-isocentre and a low magnetic field strength region at the end of the bore whilst maintaining consistent set-up and scatter conditions. Three PTW 34001 Roos chambers were investigated as well as a PTW 30013 Farmer used to validate methodology.Main Results.The beam quality used for the measurements ofkB⃗wasTPR_20/10 = 0.632. ThekB⃗for the PTW Farmer chamber at 1 T on a parallel MRI-linac was 0.993 ± 0.013 (k = 1). The averagekB⃗factor measured for the three Roos chambers on a 1 T parallel MRI-linac was 0.999 ± 0.014 (k = 1).Significance.The results presented are the first measurements ofkB⃗for a Roos chamber on a parallel MRI-linac. The Roos chamber results demonstrate the potential for the chamber as a reference dosimeter in parallel MRI-linacs.

Compact bunker shielding assessment for 1.5 T MR-Linac

This study evaluated the effect of the 1.5 T magnetic field of the magnetic resonance-guided linear accelerator (MR-Linac) on the radiation leakage doses penetrating the bunker radiation shielding wall. The evaluated 1.5 T MR-Linac Unity system has a bunker of the minimum recommended size. Unlike a conventional Linac, both primary beam transmission and secondary beam leakage were considered independently in the design and defined at the machine boundary away from the isocenter. Moreover, additional shielding was designed considering the numerous ducts between the treatment room and other rooms. The Linac shielding was evaluated by measuring the leakage doses at several locations. The intrinsic vibration and magnetic field were inspected at the proposed isocenter of the system. For verification, leakage doses were measured before and after applying the magnetic field. The intrinsic vibration and magnetic field readings were below the permitted limit. The leakage dose (0.05–12.2 µSv/week) also complied with internationally stipulated limits. The special shielding achieved a five-fold reduction in leakage dose. Applying the magnetic field increased the leakage dose by 0.12 to 4.56 µSv/week in several measurement points, although these values fall within experimental uncertainty. Thus, the effect of the magnetic field on the leakage dose could not be ascertained.

Quality management in radiotherapy treatment delivery

Radiation Oncology continues to rely on accurate delivery of radiation, in particular where patients can benefit from more modulated and hypofractioned treatments that can deliver higher dose to the target while optimising dose to normal structures. These deliveries are more complex, and the treatment units are more computerised, leading to a re-evaluation of quality assurance (QA) to test a larger range of options with more stringent criteria without becoming too time and resource consuming. This review explores how modern approaches of risk management and automation can be used to develop and maintain an effective and efficient QA programme. It considers various tools to control and guide radiation delivery including image guidance and motion management. Links with typical maintenance and repair activities are discussed, as well as patient-specific quality control activities. It is demonstrated that a quality management programme applied to treatment delivery can have an impact on individual patients but also on the quality of treatment techniques and future planning. Developing and customising a QA programme for treatment delivery is an important part of radiotherapy. Using modern multidisciplinary approaches can make this also a useful tool for department management.

Commissioning measurements on an Elekta Unity MR-Linac

Magnetic resonance-guided radiotherapy technology is relatively new and commissioning publications, quality assurance (QA) protocols and commercial products are limited. This work provides guidance for implementation measurements that may be performed on the Elekta Unity MR-Linac (Elekta, Stockholm, Sweden). Adaptations of vendor supplied phantoms facilitated determination of gantry angle accuracy and linac isocentre, whereas in-house developed phantoms were used for end-to-end testing and anterior coil attenuation measurements. Third-party devices were used for measuring beam quality, reference dosimetry and during treatment plan commissioning; however, due to several challenges, variations on standard techniques were required. Gantry angle accuracy was within 0.1°, confirmed with pixel intensity profiles, and MV isocentre diameter was < 0.5 mm. Anterior coil attenuation was approximately 0.6%. Beam quality as determined by TPR_20,10 was 0.705 ± 0.001, in agreement with treatment planning system (TPS) calculations, and gamma comparison against the TPS for a 22.0 × 22.0 cm² field was above 95.0% (2.0%, 2.0 mm). Machine output was 1.000 ± 0.002 Gy per 100 MU, depth 5.0 cm. During treatment plan commissioning, sub-standard results indicated issues with machine behaviour. Once rectified, gamma comparisons were above 95.0% (2.0%, 2.0 mm). Centres which may not have access to specialized equipment can use in-house developed phantoms, or adapt those supplied by the vendor, to perform commissioning work and confirm operation of the MRL within published tolerances. The plan QA techniques used in this work can highlight issues with machine behaviour when appropriate gamma criteria are set.

Monte Carlo investigation of electron fluence perturbation in MRI-guided radiotherapy beams using six commercial radiation detectors

With the integration of treatments with MRI-linacs to the clinical workflow, the understanding and characterization of detector response in reference dosimetry in magnetic fields are required. The external magnetic field perturbs the electron fluence. The degree of perturbation depends on the irradiation conditions and on the detector type. The purpose of this study is to evaluate the magnetic field impact on the electron fluence spectra in several detectors to provide a deeper understanding of detector response in these conditions. Monte Carlo calculations of the electron fluence are performed in six detectors (solid-state: PTW60012 and PTW60019, ionization chambers: PTW30013, PTW31010, PTW31021, and PTW31022) in water and irradiated by a 7 MV FFF photon beam with a small and a reference field, at 0 and 1.5 T. Three chamber axis orientations are investigated: parallel or perpendicular (either the Lorentz force pointing towards the stem or the tip) to the magnetic field and always perpendicular to the photon beam. One orientation for the solid-state detector is studied: parallel to the photon beam and perpendicular to the magnetic field. Additionally, electron fluence spectra are calculated in modified detector geometries to identify the underlying physical mechanisms behind the fluence perturbations. The total electron fluence in the Farmer chamber varies up to 1.24% and 5.12% at 1.5 T, in the parallel and perpendicular orientation, respectively. The interplay between the gyration radius and the Farmer chamber cavity length significantly affects the electron fluence in the perpendicular orientation. For the small-cavity chambers, the maximal variation in total electron fluence is 0.19% in the parallel orientation for the reference field. Significant small-field effects occur in these chambers; the magnetic field reduces the total electron fluence (with respect to the no field case) between 9.86% and 14.50%, depending on the orientation. The magnetic field strongly impacted the solid-state detectors in both field sizes, probably due to the high-Z components and cavity density. The maximal reductions of total electron fluence are 15.06 ± 0.09% (silicon) and 16.00 ± 0.07% (microDiamond). This work provides insights into detector response in magnetic fields by illustrating the interplay between several factors causing dosimetric perturbation effects: (1) chamber and magnetic field orientation, (2) cavity size and shape, (3) extracameral components, (4) air gaps and their asymmetry, (5) electron energy. Low-energy electron trajectories are more susceptible to change in magnetic fields, and are associated with detector response perturbation. Detectors with higher density and high-Z extracameral components exhibit more significant perturbations in the presence of a magnetic field, regardless of field size.

A portable magnet for radiation biology and dosimetry studies in magnetic fields

BACKGROUND AND PURPOSE

In the current and rapidly evolving era of real-time MRI-guided radiotherapy, our radiation biology and dosimetry knowledge is being tested in a novel way. This paper presents the successful design and implementation of a portable device used to generate strong localized magnetic fields. These are ideally suited for small scale experiments that mimic the magnetic field environment inside an MRI-linac system, or more broadly MRI-guided particle therapy as well.

MATERIALS AND METHODS

A portable permanent magnet based device employing an adjustable steel yoke and magnetic field focusing cones has been designed, constructed and tested. The apparatus utilises two banks of Nd₂ Fe₁₄ B permanent magnets totalling around 50 kg in mass to generate a strong magnetic field throughout a small volume between two pole tips. The yoke design allows adjustment of the pole tip gap and exchanging of the focusing cones. Further to this, beam portal holes are present in the yoke and focusing cones, allowing for radiation beams of up to 5 x 5 cm² to pass through the region of high magnetic field between the focusing cone tips. Finite element magnetic modelling was performed to design and characterise the performance of the device. Automated physical measurements of the magnetic field components at various locations were measured to confirm the performance. The adjustable pole gap and interchangeable cones allows rapid changing of the experimental set-up to allow different styles of measurements to be performed.

RESULTS

A mostly uniform magnetic field of 1.2 T can be achieved over a volume of at least 3 x 3 x 3 cm³ . This can be reduced in strength to 0.3 T but increased in volume to 10 x 10 x 10 cm³ via removal of the cone tips and/or adjustment of the steel yoke. Although small, these volumes are sufficient to house radiation detectors, cell culture dishes and various phantom arrangements targeted at examining small radiation field dosimetry inside magnetic field strengths that can be changed with ease. Most important is the ability to align the magnetic field both perpendicular to, or inline with the radiation beam. To date, the system has been successfully used to conduct published research in the areas of radiation detector performance, lung phantom dosimetry, and how small clinical electron beams behave in these strong magnetic fields.

CONCLUSIONS

A portable, relatively inexpensive, and simple to operate device has successfully been constructed and used for performing radiation oncology studies around the theme of MRI-guided radiotherapy. This can be in either inline and perpendicular magnetic fields of up to 1.2 T with x-ray and particle beams.

A Geant4 Fano test for novel very high energy electron beams

Objective.The boundary crossing algorithm available in Geant4 10.07-p01 general purpose Monte Carlo code has been investigated for a 12 and 200 MeV electron source by the application of a Fano cavity test.Approach.Fano conditions were enforced through all simulations whilst varying individual charged particle transport parameters which control particle step size, ionisation and single scattering.Main Results.At 12 MeV, Geant4 was found to return excellent dose consistency within 0.1% even with the default parameter configurations. The 200 MeV case, however, showed significant consistency issues when default physics parameters were employed with deviations from unity of more than 6%. The effect of the inclusion of nuclear interactions was also investigated for the 200 MeV beam and was found to return good consistency for a number of parameter configurations.Significance.The Fano test is a necessary investigation to ensure the consistency of charged particle transport available in Geant4 before detailed detector simulations can be conducted.

Dosimetric evaluation of irradiation geometry and potential air gaps in an acrylic miniphantom used for external audit of absolute dose calibration for a hybrid 1.5 T MR-linac system

Introduction

To investigate the impact of partial lateral scatter (LS), backscatter (BS) and presence of air gaps on optically stimulated luminescence dosimeter (OSLD) measurements in an acrylic miniphantom used for dosimetry audit on the 1.5 T magnetic resonance‐linear accelerator (MR‐linac) system.

Methods

The following irradiation geometries were investigated using OSLDs, A26 MR/A12 MR ion chamber (IC), and Monaco Monte Carlo system: (a) IC/OSLD in an acrylic miniphantom (partial LS, partial BS), (b) IC/OSLD in a miniphantom placed on a solid water (SW) stack at a depth of 1.5 cm (partial LS, full BS), (c) IC/OSLD placed at a depth of 1.5 cm inside a 3 cm slab of SW/buildup (full LS, partial BS), and (d) IC/OSLD centered inside a 3 cm slab of SW/buildup at a depth of 1.5 cm placed on top of a SW stack (full LS, full BS). Average of two irradiated OSLDs with and without water was used at each setup. An air gap of 1 and 2 mm, mimicking presence of potential air gap around the OSLDs in the miniphantom geometry was also simulated. The calibration condition of the machine was 1 cGy/MU at SAD = 143.5 cm, d = 5 cm, G90, and 10 × 10 cm².

Results

The Monaco calculation (0.5% uncertainty and 1.0 mm voxel size) for the four setups at the measurement point were 108.2, 108.1, 109.4, and 110.0 cGy. The corresponding IC measurements were 109.0 ± 0.03, 109.5 ± 0.06, 110.2 ± 0.02, and 109.8 ± 0.03 cGy. Without water, OSLDs measurements were ∼10% higher than the expected. With added water to minimize air gaps, the measurements were significantly improved to within 2.2%. The dosimetric impacts of 1 and 2 mm air gaps were also verified with Monaco to be 13.3% and 27.9% higher, respectively, due to the electron return effect.

Conclusions

A minimal amount of air around or within the OSLDs can cause measurement discrepancies of 10% or higher when placed in a high b‐field MR‐linac system. Care must be taken to eliminate the air from within and around the OSLD.

Evaluation of the electron transport algorithm in magnetic field in EGS5 Monte Carlo code

PURPOSE

To evaluate the accuracy of electron transport in the magnetic field of Electron Gamma Shower version 5 (EGS5) by using the special Fano cavity test.

METHODS

To simulate electron transport in the magnetic field, the trajectory of the electron was reconstructed with a short step length to restrict fractional energy loss, and the maximum user step length (mxustep) was set at 0.01 cm or 0.001 cm. For the special Fano cavity test, three-layer slab Fano test geometry was used, and uniform and isotropic per unit mass mono-energetic electrons with 0.01, 0.1, 1.0, and 10 MeV were permitted from the central axis of geometry in 0.35 T and 1.5 T. Furthermore, the magnetic field strength was scaled based on the mass density of the material. The relative difference between the calculated dose to gap and the theoretical value was evaluated. Furthermore, the special Fano cavity test was also performed using EGSnrc with the electron-enhanced electric and magnetic field macros under the same conditions, and the results were compared with those of EGS5.

RESULTS

Deviations in 0.35 T were within 0.3% regardless of the parameter settings. In 1.5 T, stable results within 0.3% were obtained using 0.001 cm as the mxustep, except for one at 10 MeV. Further, the accuracy of EGSnrc was within 0.2%, except for 10 MeV for a 0.2-cm gap in 1.5 T.

CONCLUSIONS

EGS5 with the appropriate parameter settings enable electron transport in magnetic fields similar with the accuracy of EGSnrc.

Monte Carlo calculation of detector perturbation and quality correction factors in a 1.5 T magnetic resonance guided radiation therapy small photon beams

Objective.With future advances in magnetic resonance imaging-guided radiation therapy, small photon beams are expected to be included regularly in clinical treatments. This study provides physical insights on detector dose-response to multiple megavoltage photon beam sizes coupled to magnetic fields and determines optimal orientations for measurements.Approach.Monte Carlo simulations determine small-cavity detector (solid-state: PTW60012 and PTW60019, ionization chambers: PTW31010, PTW31021, and PTW31022) dose-responses in water to an Elekta Unity 7 MV FFF photon beam. Investigations are performed for field widths between 0.25 and 10 cm in four detector axis orientations with respect to the 1.5 T magnetic field and the photon beam. The magnetic field effect on the overall perturbation factor (P_MC) accounting for the extracameral components, atomic composition, and density is quantified in each orientation. The density (P_ρ) and volume averaging (P_vol) perturbation factors and quality correction factors (kQB,QfB,f) accounting for the magnetic field are also calculated in each orientation.Main results.Results show thatP_volremains the most significant perturbation both with and without magnetic fields. In most cases, the magnetic field effect onP_volis 1% or less. The magnetic field effect onP_ρis more significant on ionization chambers than on solid-state detectors. This effect increases up to 1.564 ± 0.001 with decreasing field size for chambers. On the contrary, the magnetic field effect on the extracameral perturbation factor is higher on solid-state detectors than on ionization chambers. For chambers, the magnetic field effect onP_MCis only significant for field widths <1 cm, while, for solid-state detectors, this effect exhibits different trends with orientation, indicating that the beam incident angle and geometry play a crucial role.Significance.Solid-state detectors' dose-response is strongly affected by the magnetic field in all orientations. The magnetic field impact on ionization chamber response increases with decreasing field size. In general, ionization chambers yieldkQB,QfB,fcloser to unity, especially in orientations where the chamber axis is parallel to the magnetic field.

Longitudinal assessment of quality assurance measurements in a 1.5T MR-linac: Part I-Linear accelerator

Purpose

To describe and report longitudinal quality assurance (QA) measurements for the mechanical and dosimetric performance of an Elekta Unity MR‐linac during the first year of clinical use in our institution.

Materials and methods

The mechanical and dosimetric performance of the MR‐linac was evaluated with daily, weekly, monthly, and annual QA testing. The measurements monitor the size of the radiation isocenter, the MR‐to‐MV isocenter concordance, MLC and jaw position, the accuracy and reproducibility of step‐and‐shoot delivery, radiation output and beam profile constancy, and patient‐specific QA for the first 50 treatments in our institution. Results from end‐to‐end QA using anthropomorphic phantoms are also included as a reference for baseline comparisons. Measurements were performed in water or water‐equivalent plastic using ion chambers of various sizes, an ion chamber array, MR‐compatible 2D/3D diode array, portal imager, MRI, and radiochromic film.

Results

The diameter of the radiation isocenter and the distance between the MR/MV isocenters was (μ ± σ) 0.39 ± 0.01 mm and 0.89 ± 0.05 mm, respectively. Trend analysis shows both measurements to be well within the tolerance of 1.0 mm. MLC and jaw positional accuracy was within 1.0 mm while the dosimetric performance of step‐and‐shoot delivery was within 2.0%, irrespective of gantry angle. Radiation output and beam profile constancy were within 2.0% and 1.0%, respectively. End‐to‐end testing performed with ion‐chamber and radiochromic film showed excellent agreement with treatment plan. Patient‐specific QA using a 3D diode array identified gantry angles with low‐pass rates allowing for improvements in plan quality after necessary adjustments.

Conclusion

The MR‐linac operates within the guidelines of current recommendations for linear accelerator performance, stability, and safety. The analysis of the data supports the recently published guidance in establishing clinically acceptable tolerance levels for relative and absolute measurements.

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.

Quality Assurance for Small-Field VMAT SRS and Conventional-Field IMRT Using the Exradin W1 Scintillator

BACKGROUND

Plastic scintillator detector (PSD) Exradin W1 has shown promising performance in small field dosimetry due to its water equivalence and small sensitive volume. However, few studies reported its capability in measuring fields of conventional sizes. Therefore, the purpose of this study is to assess the performance of W1 in measuring point dose of both conventional IMRT plans and VMAT SRS plans.

METHODS

Forty-seven clinical plans (including 29 IMRT plans and 18 VMAT SRS plans with PTV volume less than 8 cm³) from our hospital were included in this study. W1 and Farmer-Type ionization chamber Exradin A19 were used in measuring IMRT plans, and W1 and microchamber Exradin A16 were used in measuring SRS plans. The agreement between the results of different types of detectors and TPS was evaluated.

RESULTS

For IMRT plans, the average differences between measurements and TPS in high-dose regions were 0.27% ± 1.66% and 0.90% ± 1.78% (P = 0.056), and were -0.76% ± 1.47% and 0.37% ± 1.34% in low-dose regions (P = 0.000), for W1 and A19, respectively. For VMAT SRS plans, the average differences between measurements and TPS were -0.19% ± 0.96% and -0.59% ± 1.49% for W1 and A16 with no statistical difference (P = 0.231).

CONCLUSION

W1 showed comparable performance with application-dedicated detectors in point dose measurements for both conventional IMRT and VMAT SRS techniques. It is a potential one-stop solution for general radiotherapy platforms that deliver both IMRT and SRS plans.

The dose response of PTW microDiamond and microSilicon in transverse magnetic field under small field conditions

The aim of the present work is to investigate the behavior of two diode-type detectors (PTW microDiamond 60019 and PTW microSilicon 60023) in transverse magnetic field under small field conditions. A formalism based on TRS 483 has been proposed serving as the framework for the application of these high-resolution detectors under these conditions. Measurements were performed at the National Metrology Institute of Germany (PTB, Braunschweig) using a research clinical linear accelerator facility. Quadratic fields corresponding to equivalent square field sizesSbetween 0.63 and 4.27 cm at the depth of measurement were used. The magnetic field strength was varied up to 1.4 T. Experimental results have been complemented with Monte Carlo simulations up to 1.5 T. Detailed simulations were performed to quantify the small field perturbation effects and the influence of detector components on the dose response. The does response of both detectors decreases by up to 10% at 1.5 T in the largest field size investigated. InS = 0.63 cm, this reduction at 1.5 T is only about half of that observed in field sizesS > 2 cm for both detectors. The results of the Monte Carlo simulations show agreement better than 1% for all investigated conditions. Due to normalization at the machine specific reference field, the resulting small field output correction factors for both detectors in magnetic fieldkQclin,QmsrBare smaller than those in the magnetic field-free case, where correction up to 6.2% at 1.5 T is required for the microSilicon in the smallest field size investigated. The volume-averaging effect of both detectors was shown to be nearly independent of the magnetic field. The influence of the enhanced-density components within the detectors has been identified as the major contributors to their behaviors in magnetic field. Nevertheless, the effect becomes weaker with decreasing field size that may be partially attributed to the deficiency of low energy secondary electrons originated from distant locations in small fields.

Acceptance procedure for the linear accelerator component of the 1.5 T MRI-linac

Purpose

To develop and implement an acceptance procedure for the new Elekta Unity 1.5 T MRI‐linac.

Methods

Tests were adopted and, where necessary adapted, from AAPM TG106 and TG142, IEC 60976 and NCS 9 and NCS 22 guidelines. Adaptations were necessary because of the atypical maximum field size (57.4 × 22 cm), FFF beam, the non‐rotating collimator, the absence of a light field, the presence of the 1.5 T magnetic field, restricted access to equipment within the bore, fixed vertical and lateral table position, and the need for MR image to MV treatment alignment. The performance specifications were set for stereotactic body radiotherapy (SBRT).

Results

The new procedure was performed similarly to that of a conventional kilovoltage x‐ray (kV) image guided radiation therapy (IGRT) linac. Results were acquired for the first Unity system.

Conclusions

A comprehensive set of tests was developed, described and implemented for the MRI‐linac. The MRI‐linac met safety requirements for patients and operators. The system delivered radiation very accurately with, for example a gantry rotation locus of isocenter of radius 0.38 mm and an average MLC absolute positional error of 0.29 mm, consistent with use for SBRT. Specifications for clinical introduction were met.

Technical Note: End-to-end verification of an MR-Linac using a dynamic motion phantom

PURPOSE

MR-Linac integrates an MRI scanner and a linear accelerator to provide adaptive radiation treatment. Superior tissue contrast and real-time imaging can give the clinicians confidence to reduce the margins of the planning target volume (PTV). The purpose of this study was to verify the dosimetric accuracy of an MR-Linac system in treating a moving target and assess the error with different motion patterns and adaptation methods.

METHODS

We performed an end-to-end test for Elekta Unity (Elekta) using the 4D Dynamic Thorax Phantom (CIRS MRgRT 008Z), comparing the measured and planned dose. The moving phantom had four measurement locations in the tumor, liver, kidney, and spinal cord regions with a PTW30013 ion chamber. For seven different motion patterns, we first acquired simulation CT using a slow-scanning protocol, based on which we generated reference plans. The treatment technique was the standard intensity-modulated radiation therapy (IMRT). We tested both adaptation workflows: the Adapt-to-Position (ATP) and the Adapt-to-Shape (ATS). The three-dimensional (3D) distribution was measured using a diode array phantom (Sun Nuclear Inc.) to check the dose distribution accuracy as part of the routine QA process. We also performed end-to-end tests on a conventional Linac. Finally, we used SPSS Statistics 22.0 (Inc., Chicago, IL, USA) for data analysis.

RESULTS

All pretreatment reference plans and delivered plans had excellent QA results with a better than 95% passing rate of relative gamma analysis (2%/2 mm criteria). The adaptive planning for MR-Linac produced quality plans. The measured dose in the target agreed with the calculated dose.

CONCLUSIONS

The adaptive treatment on the MR-Linac system investigated met the expected performance with tumor motions. The outline of the target could be visualized and accurately contoured on the 3D MR for online planning. Under different motion patterns, the difference between the measured and calculated dose was acceptable clinically.

Automatic 3D Monte-Carlo-based secondary dose calculation for online verification of 1.5 T magnetic resonance imaging guided radiotherapy

•

First implementation of an independent 3D-secondary dose calculation (3D-SDC).

•

Validation of the 3D-SDC solution using patient plans and experimental plan QA.

•

Online SDC of central targets is feasible with a median calculation time of 1:23 min.

•

Peripheral targets with small beam numbers need alternative validation strategies.

Background and purpose

Hybrid magnetic resonance linear accelerator (MR-Linac) systems represent a novel technology for online adaptive radiotherapy. 3D secondary dose calculation (SDC) of online adapted plans is required to assure patient safety. Currently, no 3D-SDC solution is available for 1.5T MR-Linac systems. Therefore, the aim of this project was to develop and validate a method for online automatic 3D-SDC for adaptive MR-Linac treatments.

Materials and methods

An accelerator head model was designed for an 1.5T MR-Linac system, neglecting the magnetic field. The use of this model for online 3D-SDC of MR-Linac plans was validated in a three-step process: (1) comparison to measured beam data, (2) investigation of performance and limitations in a planning phantom and (3) clinical validation using n = 100 patient plans from different tumor entities, comparing the developed 3D-SDC with experimental plan QA.

Results

The developed model showed median gamma passing rates compared to MR-Linac base data of 84.7%, 100% and 99.1% for crossplane, inplane and depth-dose-profiles, respectively. Comparison of 3D-SDC and full dose calculation in a planning phantom revealed that with ⩾5 beams gamma passing rates >95% can be achieved for central target locations. With a median calculation time of 1:23 min, 3D-SDC of online adapted clinical MR-Linac plans demonstrated a median gamma passing rate of 98.9% compared to full dose calculation, whereas experimental plan QA reached 99.5%.

Conclusion

Here, we describe the first technical 3D-SDC solution for online adaptive MR-guided radiotherapy. For clinical situations with peripheral targets and a small number of beams additional verification appears necessary. Further improvement may include 3D-SDC with consideration of the magnetic field.

MR-guided proton therapy: Impact of magnetic fields on the detector response

PURPOSE

To investigate the response of detectors for proton dosimetry in the presence of magnetic fields.

MATERIAL AND METHODS

Four ionization chambers (ICs), two thimble-type and two plane-parallel-type, and a diamond detector were investigated. All detectors were irradiated with homogeneous single-energy-layer fields, using 252.7 MeV proton beams. A Farmer IC was additionally irradiated in the same geometrical configuration, but with a lower nominal energy of 97.4 MeV. The beams were subjected to magnetic field strengths of 0, 0.25, 0.5, 0.75, and 1 T produced by a research dipole magnet placed at the room's isocenter. Detectors were positioned at 2 cm water equivalent depth, with their stem perpendicular to both the magnetic field lines and the proton beam's central axis, in the direction of the Lorentz force. Normality and two sample statistical Student's t tests were performed to assess the influence of the magnetic field on the detectors' responses.

RESULTS

For all detectors, a small but significant magnetic field-dependent change of their response was found. Observed differences compared to the no magnetic field case ranged from +0.5% to -0.7%. The magnetic field dependence was found to be nonlinear and highest between 0.25 and 0.5 T for 252.7 MeV proton beams. A different variation of the Farmer chamber response with magnetic field strength was observed for irradiations using lower energy (97.4 MeV) protons. The largest magnetic field effects were observed for plane-parallel ionization chambers.

CONCLUSION

Small magnetic field-dependent changes in the detector response were identified, which should be corrected for dosimetric applications.

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.

Machine QA for the Elekta Unity system: A Report from the Elekta MR-linac consortium

Over the last few years, magnetic resonance image‐guided radiotherapy systems have been introduced into the clinic, allowing for daily online plan adaption. While quality assurance (QA) is similar to conventional radiotherapy systems, there is a need to introduce or modify measurement techniques. As yet, there is no consensus guidance on the QA equipment and test requirements for such systems. Therefore, this report provides an overview of QA equipment and techniques for mechanical, dosimetric, and imaging performance of such systems and recommendation of the QA procedures, particularly for a 1.5T MR‐linac device. An overview of the system design and considerations for QA measurements, particularly the effect of the machine geometry and magnetic field on the radiation beam measurements is given. The effect of the magnetic field on measurement equipment and methods is reviewed to provide a foundation for interpreting measurement results and devising appropriate methods. And lastly, a consensus overview of recommended QA, appropriate methods, and tolerances is provided based on conventional QA protocols. The aim of this consensus work was to provide a foundation for QA protocols, comparative studies of system performance, and for future development of QA protocols and measurement methods.

Monte Carlo calculation of perturbation correction factors for air-filled ionization chambers in clinical proton beams using TOPAS/GEANT

INTRODUCTION

Current dosimetry protocols for clinical protons using air-filled ionization chambers assume that the perturbation correction factor is equal to unity for all ionization chambers and proton energies. Since previous Monte Carlo based studies suggest that perturbation correction factors might be significantly different from unity this study aims to determine perturbation correction factors for six plane-parallel and four cylindrical ionization chambers in proton beams at clinical energies.

MATERIALS AND METHODS

The dose deposited in the air cavity of the ionization chambers was calculated with the help of the Monte Carlo code TOPAS/Geant4 while specific constructive details of the chambers were removed step by step. By comparing these dose values the individual perturbation correction factors p_cel, p_stem, p_sleeve, p_wall, p_cav⋅p_dis as well as the total perturbation correction factor p_Q were derived for typical clinical proton energies between 80 and 250MeV.

RESULTS

The total perturbation correction factor p_Q was smaller than unity for almost every ionization chamber and proton energy and in some cases significantly different from unity (deviation larger than 1%). The maximum deviation from unity was 2.0% for cylindrical and 1.5% for plane-parallel ionization chambers. Especially the factor p_wall was found to differ significantly from unity. It was shown that this is due to the fact that secondary particles, especially alpha particles and fragments, are scattered from the chamber wall into the air cavity resulting in an overresponse of the chamber.

CONCLUSION

Perturbation correction factors for ionization chambers in proton beams were calculated using Monte Carlo simulations. In contrast to the assumption of current dosimetry protocols the total perturbation correction factor p_Q can be significantly different from unity. Hence, beam quality correction factors [Formula: see text] that are calculated with the help of perturbation correction factors that are assumed to be unity come with a corresponding additional uncertainty.

Technical Challenges of Real-Time Adaptive MR-Guided Radiotherapy

In the past few years, radiotherapy (RT) has experienced a major technological innovation with the development of hybrid machines combining magnetic resonance (MR) imaging and linear accelerators. This new technology for MR-guided cancer treatment has the potential to revolutionize the field of adaptive RT due to the opportunity to provide high-resolution, real-time MR imaging before and during treatment application. However, from a technical point of view, several challenges remain which need to be tackled to ensure safe and robust real-time adaptive MR-guided RT delivery. In this manuscript, several technical challenges to MR-guided RT are discussed. Starting with magnetic field strength tradeoffs, the potential and limitations for purely MR-based RT workflows are discussed. Furthermore, the current status of real-time 3D MR imaging and its potential for real-time RT are summarized. Finally, the potential of quantitative MR imaging for future biological RT adaptation is highlighted.

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.

End-to-end validation of the geometric dose delivery performance of MR linac adaptive radiotherapy

The clinical introduction of hybrid magnetic resonance (MR) guided radiotherapy (RT) delivery systems has led to the need to validate the end-to-end dose delivery performance on such machines. In the current study, an MR visible phantom was developed and used to test the spatial deviation between planned and delivered dose at two 1.5 T MR linear accelerator (MR linac) systems, including pre-treatment imaging, dose planning, online imaging, image registration, plan adaptation, and dose delivery. The phantom consisted of 3D printed plastic and MR visible silicone rubber. It was designed to minimise air gaps close to the radiochromic film used as a dosimeter. Furthermore, the phantom was designed to allow submillimetre, reproducible positioning of the film in the phantom. At both MR linac systems, 54 complete adaptive, MR guided RT workflow sessions were performed. To test the dose delivery performance of the MR linac systems in various adaptive RT (ART) scenarios, the sessions comprised a range of systematic positional shifts of the phantom and imaging or plan adaptation conditions. In each workflow session, the positional translation between the film and the adaptive planned dose was determined. The results showed that the accuracy of the MR linac systems was between 0.1 and 0.9 mm depending on direction. The highest mean deviance observed was in the posterior-anterior direction, and the direction of the error was consistent between centres. The precision of the systems was related to whether the workflow utilized the internal image registration algorithm of the MR linac. Workflows using the internal registration algorithm led to a worse precision (0.2-0.7 mm) compared to workflows where the algorithm was decoupled (0.2 mm). In summary, the spatial deviation between planned and delivered dose of MR-guided ART at the two MR linac systems was well below 1 mm and thus acceptable for clinical use.

Towards the standardization of the absorbed dose report mode in high energy photon beams

The benefits of using an algorithm that reports absorbed dose-to-medium have been jeopardized by the clinical experience and the experimental protocols that have mainly relied on absorbed dose-to-water. The aim of the present work was to investigate the physical aspects that govern the dosimetry in heterogeneous media using Monte Carlo method and to introduce a formalism for the experimental validation of absorbed dose-to-medium reporting algorithms. Particle fluence spectra computed within the sensitive volume of two simulated detectors (T31016 Pinpoint 3D ionization chamber and EBT3 radiochromic film) placed in different media (water, RW3, lung and bone) were compared to those in the undisturbed media for 6 MV photon beams. A heterogeneity correction factor that takes into account the difference between the detector perturbation in medium and under reference conditions as well as the stopping-power ratios was then derived for all media using cema calculations. Furthermore, the different conversion approaches and Eclipse treatment planning system algorithms were compared against the Monte Carlo absorbed dose reports. The detectors electron fluence perturbation in RW3 and lung media were close to that in water (≤1.5%). However, the perturbation was greater in bone (∼4%) and impacted the spectral shape. It was emphasized that detectors readings should be corrected by the heterogeneity correction factor that ranged from 0.932 in bone to 0.985 in lung. Significant discrepancies were observed between all the absorbed dose reports and conversions, especially in bone (exceeding 10%) and to a lesser extent in RW3. Given the ongoing advances in dose calculation algorithms, it is essential to standardize the absorbed dose report mode with absorbed dose-to-medium as a favoured choice. It was concluded that a retrospective conversion should be avoided and switching from absorbed dose-to-water to absorbed dose-to-medium reporting algorithm should be carried out by a direct comparison of both algorithms.

Monte Carlo modeling of the influence of strong magnetic fields on the stem-effect in plastic scintillation detectors used in radiotherapy dosimetry

PURPOSE

To investigate the impact of strong magnetic fields on the stem-effect in plastic scintillation detectors (PSDs) using Monte Carlo methods.

METHODS

Prior to building the light guide model, the properties of the Čerenkov process in GEANT4 were investigated by simulating depth-dose and depth-Čerenkov emission profiles in water as functions of Čerenkov process input parameters. In addition, profile simulations were performed for magnetic field strengths ranging from 0 T to 1.5 T. A PMMA light guide was constructed in GEANT4 using data from the manufacturer and literature. Simulations were performed with the model as functions of depth and fiber-beam angle where the simulated stem-effect spectrum and the Čerenkov light ratio (CLR) were scored and compared to measured data in the literature. The light guide optical properties were iteratively adjusted until agreement between the simulated and measured data was achieved. Simulations were performed with the validated model as functions of depth and magnetic field strength and the simulated data were compared to measured data in the literature. The model was also used to evaluate the sensitivity of the CLR to the various optical properties of the light guide in different irradiation conditions.

RESULTS

No significant changes in the depth-dose or depth-Čerenkov emission profiles were observed with step-size restrictions imposed by the Čerenkov process input parameters, which was attributed to the condensed history algorithm and transport parameters used in this work. Similar changes in the depth-dose and depth-Čerenkov emission profiles were observed with increasing magnetic field strength, which indicates the Čerenkov process is not adversely impacted by the presence of the magnetic field. Following optimization of the light guide optical properties, agreement within two standard deviations was observed between the simulated and measured optical data for all validation geometries considered. Agreement within one standard deviation was observed between the simulated and measured data for all depths and field strengths ≥0 T whereas discrepancies were observed for magnetic field strengths <-0.35 T. These significant differences were attributed to insufficient measurement data for this irradiation configuration during model validation. Of the light guide optical properties investigated, the fluorescence signal had the greatest impact on the CLR sensitivity to the magnetic field.

CONCLUSIONS

No significant change in the Čerenkov emission per dose in water was observed for magnetic field strengths up to 1.5 T. The nominal fiber fluorescence signal was found to have a significant impact on the CLR sensitivity to varying irradiation conditions where changes up to 11.7% were observed whereas the mirror reflectivity and fiber attenuation had a modest impact with maximum CLR changes of 2.6% and 1.2% relative to 0 T, respectively. The results of this work suggest light guides with low fiber fluorescence should be used with PSDs for dosimetry measurements in magnetic fields to minimize the impact of the magnetic field on the CLR correction.