In previous RENEB interlaboratory comparisons based on the manual scoring of dicentric chromosomes, a tendency for systematic overestimation for doses > 2.5 Gy was found. However, these exercises included only very few doses in the high dose range, and they were heterogeneous in terms of radiation quality and evaluation mode, and comparable only to a limited extent. Here, this presumed deviation was explored by investigating three doses > 2.5 Gy. Blood samples were irradiated (2.56, 3.41 and 4.54 Gy) using a 60Co source and sent to 14 member laboratories of the RENEB network, which performed the dicentric chromosome assay (manual and/or semi-automatic scoring) and reported dose estimates. Most participants provided estimates that agreed very well with the physical reference doses and all provided dose estimates were in the correct clinical category (> 2 Gy). The previously observed tendency for a systematic bias across all laboratories was not confirmed. However, tendencies for systematic underestimation were detected for dose estimations for reference doses given in terms of absorbed dose to blood and for some participants, a laboratory-specific trend of systematic under- or overestimation was observed. The importance of regularly performed quality checks for a broad dose range became obvious to avoid misinterpretation of results.
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).
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.
Objective After a nuclear detonation, people will be exposed to varying mixtures of neutrons and gamma radiation and the biological effectiveness of mixed beams is not well known. Additionally, it is not known how far both radiation types interact, a question that is relevant for generating calibration curves for the purpose of biological dosimetry. The objective of this study was to investigate the potential impact of two different combinations of neutron and gamma radiation on gene expression and dicentric chromosomes in peripheral blood mononuclear cells (PBMC). Methods Whole blood from 3 human donors was exposed to neutrons with an energy spectrum similar to that of the Hiroshima uranium bomb, to gamma radiation from a 60Co source and to a 50:50 combination of both radiations, given in two orders of sequence. In all cases the total doses were 0.5, 0.75 and 1.0 Gy. Dicentric chromosomes were analyzed by light microscopy and the expression of six known radiation-responsive genes BBC3, CDKN1A, FDXR, GADD45A, MDM2, and XPC were analyzed by RT-qPCR. Results Per unit dose, exposure to neutrons lead to a higher level of dicentrics and gene expression as compared to gamma radiation. Dose-response relationships for both endpoints were linear, allowing calculating the expected outcome of combined exposure by arithmetic. For dicentric chromosomes, the RBE values for 60Co→neutrons, neutrons→60Co and neutrons were 4.05, 3.62 and 7.30 respectively. For gene expression the RBE values were gene-specific, but showed values in the range of 1.14–3.01 for 60Co→neutrons, 1.33–2.68 for neutrons→60Co and 1.39–3.91 for neutrons. Conclusions The results demonstrate that combined exposure to neutrons and gamma radiation, regardless of the order of sequence, leads to an additive response at both endpoints. This indicates that calibration curves for mixed beams can be constructed from dose response relationships of the single beam components.
Introduction. To be better suited for measurements at air kerma rates which are common for the calibration of detectors for dosimetry in brachytherapy and radiation protection, two graphite-wall, air-filled cavity ionization chambers (PS-10 and PS-50 built by PTW of Freiburg, Germany) were characterized for use as primary standards at the Physikalisch-Technische Bundesanstalt (PTB) for 60Co, 137Cs and 192Ir gamma ray sources. Methods. Applying experimental or Monte Carlo methods, corrections were determined for wall effects, stem scattering effects, point-source non-uniformity, saturation effects, and deviations from Spencer-Attix cavity theory. The geometrical cavity volume of the ionization chambers was determined by a special experimental method. Non-charge-collecting parts of the volume were calculated by finite-element simulations. In addition, mass electronic stopping powers and mass-energy absorption coefficients were reevaluated in close agreement with ICRU 90. To compare the results to an already established air kerma standard, measurements of air kerma rates were taken for gamma rays from all three sources. Results. In contrast to the standards currently in use, the newly characterized ionization chambers show significant wall effects up to 4.17% (192Ir). For the PS-10 ionization chamber, the geometric volume differs by 0.42% from the charge collecting volume. For the ^192Ir source, Spencer-Attix corrections were applied exceeding 0.30% for both ionization chambers. Applying the derived corrections factors, air kerma rates measured with the two ionization chambers at the three sources did not differ by more than 0.13%. Comparing the results to an already established air kerma standard, deviations ranged up to 0.61%. The uncertainty budget for the determination of air kerma rates resulted in a relative combined standard uncertainty of 0.35%. Conclusion. Both ionization chambers can be used as primary air kerma standards for 60Co, 137Cs and 192Ir gamma ray sources using the values of corresponding physical quantities and correction factors derived in this work.
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.
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 x 10 cm(2) 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 (C-B) into account, the magnetic field correction factors (k(B,Q)) 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 x 2 cm(2)), 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.
Correction factors for reference dosimetry in magnetic resonance (MR) imaging-guided radiation therapy (kB→,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, kB→,M,Q was assessed for different tissue-phantom ratios (TPR20,10). kB→,M,Q differed by 0.0029(43) between TPR20,10=0.6790(23)(6 MV linac) and TPR20,10=0.7028(14) (7 MV MR-linac) at 1.5T. The agreement was best in an orientation in which the secondary electrons were deflected to the stem of the ionization chamber.
The Metrological Electron Accelerator Facility (MELAF) of the Physikalisch-Technische Bundesanstalt (PTB) offers access to well characterized high-energy (0.5–50 MeV) electron and photon radiation fields also for external researchers. This work outlines the capabilities of the facility to foster new collaborations. As example, the experimental determination of ionization chamber typical correction factors for Magnetic-Resonance guided Radiotherapy (MRgRT) is presented.
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 k((B) over right arrow ,M,Q). This factor can be determined experimentally or calculated by means of Monte Carlo (MC) simulations. To date, a small number of experimental values for k((B) over right arrow ,M,Q) 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, k((B) over right arrow ,M,Q) 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, k((B) over right arrow ,M,Q) 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 k((B) over right arrow ,M,Q) and the experimental values for k((B) over right arrow ,M,Q) was 0.31(30)% and the minimum difference was 0.02(24)%. For the PTW 30013 ionization chambers, the experimental values for k((B) over right arrow ,M,Q) 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 k((B) over right arrow ,M,Q) may be neglectable for ionization chambers of the PTW 30013 type.
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 kB. The purpose of this study is to systematically investigate kB 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.
Abstract Today, hybrid systems of linear accelerator and MRI scanner are clinically available. Therefore it is important to investigate the feasibility of reference dosimetry with ionization chambers in the presence of a magnetic field and determine correction factors. In this work, correction factors under various conditions that influence the chamber response were experimentally investigated, using a conventional 6 MV linear accelerator together with a stand-alone magnet. We found that the correction factor for a PTW31010 ionization chamber ranges from 0.9873 to 1.009 depending on the magnetic field strength, magnetic field orientation and magnetic field size. The phantom material also does have an influence on the measured signal. Therefore, reference dosimetry with ionization chambers in the presence of a magnetic field is feasible, but requires dedicated correction factors, which depend on the experimental setup.
The interest in hybrid systems combining magnetic resonance imaging and medical linear accelerator (MR-Linac) is rapidly increasing due to the clinical availability of different systems. Reference dosimetry is a critical issue for integrating these devices into clinical practice. However, the response of ionization chambers changes according to the distinct orientation of the chamber with respect to the magnetic field. In this study, we have carried out Monte Carlo simulations to identify an optimal orientation for thimble type chambers in MRgRT reference dosimetry. Our findings suggest that an orientation where the chamber axis is parallel to the magnetic field axis should be preferred.
Abstract An integrated system of a linear accelerator and a magnetic resonance imaging (MRI) device may provide real-time imaging during radiotherapy treatments. This work investigated changes affecting the dose deposition caused by a magnetic field (B-field) transverse to the beam direction by means of Monte Carlo simulations. Two different phantoms were used: A water phantom (Ph1) and a water-air phantom (Ph2) with a 4-2-4 cm water-air-water cross section. Dose depositions were scored for B-field values of 0 T, 0.35 T, 0.5 T, 1.5 T, 3 T and 5 T. Beams were based on a precalculated photon spectrum taken from an earlier simulated Elekta 6 MV FFF accelerator. All lateral profiles in Ph1 showed a Lorentz force driven shift w.r.t. the B-field strength, presenting a steeper penumbra in the shift's direction. Depositions were shifted up to 0.3 cm for 5 T, showing a constant central axis plateau-dose or an increase by 2.3 % for small fields. Depth-dose curves in Ph1 showed a shift of the dose maximum towards the beam entrance direction for increasing B-field of up to 1.1 cm; the maximum dose was increased by 6.9 %. In Ph2, an asymmetric dose increase by up to 36.9 % was observed for 1.5 T at the water-air boundary, resulting from the electron return effect (ERE). In our scenario, B-field dependent dose shifts and local build-ups were observed, which consequently affect the resulting dose distribution and need to be considered in magnetic resonance guided radiotherapy treatment planning.