Implementation of a new EGSnrc particle source class for computed tomography: validation and uncertainty quantification.

Objective. Personalized dose monitoring and risk management are of increasing significance with the growing number of computer tomography (CT) examinations. These require high-quality Monte Carlo (MC) simulations that are of the utmost importance for the new developments in personalized CT dosimetry. This work aims to extend the MC framework EGSnrc source code with a new particle source. This, in turn, allows CT-scanner-specific dose and image calculations for any CT scanner. The novel method can be used with all modern EGSnrc user codes, particularly for the simulation of the effective dose based on DICOM images and the calculation of CT images.Approach. The new particle source can be used with input data derived by the user. The input data can be generated by the user based on a previously developed method for the experimental characterization of any CT scanner (doi.org/10.1016/j.ejmp.2015.09.006). Furthermore, the new particle source was benchmarked by air kerma measurements in an ionization chamber at a clinical CT scanner. For this, the simulated angular distribution and attenuation characteristics were compared to measurements to verify the source output free in air. In a second validation step, simulations of air kerma in a homogenous cylindrical and an anthropomorphic thorax phantom were performed and validated against experimentally determined results. A detailed uncertainty evaluation of the simulated air kerma values was developed.Main results. We successfully implemented a new particle source class for the simulation of realistic CT scans. This method can be adapted to any CT scanner. For the attenuation characteristics, there was a maximal deviation of 6.86% between the measurement and the simulation. The mean deviation for all tube voltages was 2.36% (σ= 1.6%). For the phantom measurements and simulations, all the values agreed within 5.0%. The uncertainty evaluation resulted in an uncertainty of 5.5% (k=1).

On the robustness of detective quantum efficiency within the limits of IEC 61267 RQA standard radiation qualities.

IEC 61267 allows a certain leeway regarding the establishment of radiation qualities in order to enable the use of X-ray tubes having different anode angles and inherent filtrations. This allowance has a direct impact on the calculation of the detective quantum efficiency and may potentially complicate any comparison of different imaging detectors based on this quantity. This work investigates this effect by applying computational methods. To this end, an algorithm was implemented to calculate the variation of the squared signal-to-noise ratio per air kerma for RQA standard radiation qualities and to deduce corresponding uncertainties based on GUM Supplement 2. For RQA standard radiation qualities, the results show standard uncertainties for the squared signal-to-noise ratio per air kerma of between 0.05 and 2.1%. Comparing imaging detectors based on detective quantum efficiency is associated with substantial uncertainty for some radiation qualities. This is due to the different photon fluences with respect to energy that are allowed by IEC 61267 for identical standard radiation qualities.

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

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

Chemical radiation dosimetry in magnetic fields: characterization of a Fricke-type chemical detector in 6 MV photon beams and magnetic fields up to 1.42 T.

In magnetic resonance guided radiotherapy (MRgRT) radiation dose measurements needs to be performed in the presence of a magnetic field. In this study, the influence of magnetic fields on the readings of a Fricke detector, a chemical dosimeter, have been investigated in 6 MV photon beams. This type of detector has been chosen, as the Federal Office of Metrology (METAS, Switzerland) has great experience with Fricke dosimetry and since it is not expected that this detector is greatly affected by the presence of a magnetic field. Magnetic fields with field strengths between 0 T and 1.42 T were applied during the detector irradiation. In a 5 × 10 cm2 irradiation field, the Fricke readings are affected less than 0.9% by the applied magnetic fields. Taking the altered dose distribution due to the magnetic field ([Formula: see text]) into account, the magnetic field correction factors ([Formula: see text]) for the Fricke detector at 0.35 T and 1.42 T are determined to be 0.9948 and 0.9980, respectively. These small corrections hardly exceed the measurement uncertainties. Hence, we could proof that the Fricke detector is not significantly influenced by the presence of a magnetic field. The Fricke detector was also tested for the feasibility of measuring output factors in the presence of magnetic fields. For irradiation field sizes larger than the detector (>2 × 2 cm2), comparable results were obtained as for other detectors. The output factors decrease when a magnetic field is applied. This effect is more pronounce for larger magnetic field strengths and smaller irradiation fields due to shifts of the depth dose curves and asymmetry of lateral dose profiles.

Influence of beam quality on reference dosimetry correction factors in magnetic resonance guided radiation therapy.

Correction factors for reference dosimetry in magnetic resonance (MR) imaging-guided radiation therapy ( k B → , M , Q ) are often determined in setups that combine a conventional 6 MV linac with an electromagnet. This study investigated whether results based on these measurements were applicable for a 7 MV MR-linac using Monte Carlo simulations. For a Farmer-type ionization chamber, k B → , M , Q was assessed for different tissue-phantom ratios ( TPR 20 , 10 ). k B → , M , Q differed by 0.0029 ( 43 ) between TPR 20 , 10 = 0.6790 ( 23 ) (6 MV linac) and TPR 20 , 10 = 0.7028 ( 14 ) (7 MV MR-linac) at 1.5 T . The agreement was best in an orientation in which the secondary electrons were deflected to the stem of the ionization chamber.

A finite element method for the determination of the relative response of ionization chambers in MR-linacs: simulation and experimental validation up to 1.5 T.

In magnetic resonance (MR) guided radiotherapy, the magnetic field-dependent change in the dose response of ionization chambers is typically included by means of a correction factor [Formula: see text]. This factor can be determined experimentally or calculated by means of Monte Carlo (MC) simulations. To date, a small number of experimental values for [Formula: see text] at magnetic flux densities above 1.2 T have been available to benchmark these simulations. Furthermore, MC simulations of the dose response of ionization chambers in magnetic fields (where such simulations are based on manufacturer blueprints) have been shown to converge with results that deviate considerably from experimental values for orientations where the magnetic field is perpendicular to the axis of the ionization chamber and the influence of the magnetic field is largest. In this work, [Formula: see text] was simulated for a PTW 30013 Farmer ionization chamber using an approach based on finite element simulations. First, the electrical field inside the ionization chamber was simulated using finite element methods. The collecting volume of the ionization was not defined in terms of the physical dimensions of the detector but in terms of the simulated electrical field lines inside the chamber. Then, an MC simulation of the dose response of a Farmer type chamber (PTW 30013) was performed using EGSnrc with a dedicated package to consider the effect of the magnetic field. In the second part, [Formula: see text] was determined experimentally for two different PTW 30013 ionization chambers for a range of magnetic flux densities between B = 0 and 1.5 T, covering the range of commercially available MR-linacs. In the perpendicular orientation, the maximum difference between the simulated values for [Formula: see text] and the experimental values for [Formula: see text] was 0.31(30)% and the minimum difference was 0.02(24)%. For the PTW 30013 ionization chambers, the experimental values for [Formula: see text] were 0.9679(1) and 0.9681(1) for a magnetic flux density of 1.5 T. The value resulting from the simulation was 0.967(3). The comparison of the correction factors simulated using this new approach with the experimental values determined in this study shows excellent agreement for all magnetic flux densities up to 1.5 T. Integrating the explicit simulation of the collection volume inside the ionization chambers into the MC simulation model significantly improves simulations of the chamber response in magnetic fields. The results presented suggest that intra-type variations for [Formula: see text] may be neglectable for ionization chambers of the PTW 30013 type.

Ionization chamber correction factors for MR-linacs.

Previously, readings of air-filled ionization chambers have been described as being influenced by magnetic fields. To use these chambers for dosimetry in magnetic resonance guided radiotherapy (MRgRT), this effect must be taken into account by introducing a correction factor k B. The purpose of this study is to systematically investigate k B for a typical reference setup for commercially available ionization chambers with different magnetic field strengths. The Monte Carlo simulation tool EGSnrc was used to simulate eight commercially available ionization chambers in magnetic fields whose magnetic flux density was in the range of 0-2.5 T. To validate the simulation, the influence of the magnetic field was experimentally determined for a PTW30013 Farmer-type chamber for magnetic flux densities between 0 and 1.425 T. Changes in the detector response of up to 8% depending on the magnetic flux density, on the chamber geometry and on the chamber orientation were obtained. In the experimental setup, a maximum deviation of less than 2% was observed when comparing measured values with simulated values. Dedicated values for two MR-linac systems (ViewRay MRIdian, ViewRay Inc, Cleveland, United States, 0.35 T/ 6 MV and Elekta Unity, Elekta AB, Stockholm, Sweden, 1.5 T/7 MV) were determined for future use in reference dosimetry. Simulated values for thimble-type chambers are in good agreement with experiments as well as with the results of previous publications. After further experimental validation, the results can be considered for definition of standard protocols for purposes of reference dosimetry in MRgRT.