1
|
Failing T, Hensley FW, Keil B, Zink K. Investigations on the beam quality correction factor for ionization chambers in high-energy brachytherapy dosimetry. Phys Med Biol 2024; 69:165002. [PMID: 39009012 DOI: 10.1088/1361-6560/ad638b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Accepted: 07/15/2024] [Indexed: 07/17/2024]
Abstract
Objective. To enhance the investigations on MC calculated beam quality correction factors of thimble ionization chambers from high-energy brachytherapy sources and to develop reliable reference conditions in source and detector setups in water.Approach. The response of five different ionization chambers from PTW-Freiburg and Standard Imaging was investigated for irradiation by a high dose rate Ir-192 Flexisource in water. For a setup in a Beamscan water phantom, Monte Carlo simulations were performed to calculate correction factors for the chamber readings. After exact positioning of source and detector the absorbed dose rate at the TG-43 reference point at one centimeter nominal distance from the source was measured using these factors and compared to the specification of the calibration certificate. The Monte Carlo calculations were performed using the restricted cema formalism to gain further insight into the chamber response. Calculations were performed for the sensitive volume of the chambers, determined by the methods currently used in investigations of dosimetry in magnetic fields.Main results. Measured dose rates and values from the calibration certificate agreed within the combined uncertainty (k= 2) for all chambers except for one case in which the full air cavity was simulated. The chambers showed a distinct directional dependence. With the restricted cema formalism calculations it was possible to examine volume averaging and energy dependence of the perturbation factors contributing to the beam quality correction factor also differential in energy.Significance. This work determined beam quality correction factors to measure the absorbed dose rate from a brachytherapy source in terms of absorbed dose to water for a variety of ionization chambers. For the accurate dosimetry of brachytherapy sources with ionization chambers it is advisable to use correction factors based on the sensitive volume of the chambers and to take account for the directional dependence of chamber response.
Collapse
Affiliation(s)
- T Failing
- Institute of Medical Physics and Radiation Protection (IMPS), University of Applied Sciences, Gießen, Germany
| | - F W Hensley
- Department for Radiotherapy and Radiooncology, University Medical Center Heidelberg, Heidelberg, Germany
| | - B Keil
- Institute of Medical Physics and Radiation Protection (IMPS), University of Applied Sciences, Gießen, Germany
- Department for Diagnostic and Interventional Radiology, Philipps-University Marburg, Marburg, Germany
- LOEWE Research Cluster for Advanced Medical Physics in Imaging and Therapy (ADMIT), TH Mittelhessen University of Applied Sciences, Giessen, Germany
| | - K Zink
- Institute of Medical Physics and Radiation Protection (IMPS), University of Applied Sciences, Gießen, Germany
- LOEWE Research Cluster for Advanced Medical Physics in Imaging and Therapy (ADMIT), TH Mittelhessen University of Applied Sciences, Giessen, Germany
- Department for Radiotherapy and Radiooncology, University Medical Center Giessen-Marburg, Marburg, Germany
- Marburg Iontherapy Center (MIT), Marburg, Germany
| |
Collapse
|
2
|
Failing T, Hartmann GH, Hensley FW, Keil B, Zink K. Enhancement of the EGSnrc code egs_chamber for fast fluence calculations of charged particles. Z Med Phys 2022; 32:417-427. [PMID: 35643800 PMCID: PMC9948836 DOI: 10.1016/j.zemedi.2022.04.003] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2022] [Revised: 03/29/2022] [Accepted: 04/14/2022] [Indexed: 11/19/2022]
Abstract
PURPOSE Simulation of absorbed dose deposition in a detector is one of the key tasks of Monte Carlo (MC) dosimetry methodology. Recent publications (Hartmann and Zink, 2018; Hartmann and Zink, 2019; Hartmann et al., 2021) have shown that knowledge of the charged particle fluence differential in energy contributing to absorbed dose is useful to provide enhanced insight on how response depends on detector properties. While some EGSnrc MC codes provide output of charged particle spectra, they are often restricted in setup options or limited in calculation efficiency. For detector simulations, a promising approach is to upgrade the EGSnrc code egs_chamber which so far does not offer charged particle calculations. METHODS Since the user code cavity offers charged particle fluence calculation, the underlying algorithm was embedded in egs_chamber. The modified code was tested against two EGSnrc applications and DOSXYZnrc which was modified accordingly by one of the authors. Furthermore, the gain in efficiency achieved by photon cross section enhancement was determined quantitatively. RESULTS Electron and positron fluence spectra and restricted cema calculated by egs_chamber agreed well with the compared applications thus demonstrating the feasibility of the new code. Additionally, variance reduction techniques are now applicable also for fluence calculations. Depending on the simulation setup, considerable gains in efficiency were obtained by photon cross section enhancement. CONCLUSION The enhanced egs_chamber code represents a valuable tool to investigate the response of detectors with respect to absorbed dose and fluence distribution and the perturbation caused by the detector in a reasonable computation time. By using intermediate phase space scoring, egs_chamber offers parallel calculation of charged particle fluence spectra for different detector configurations in one single run.
Collapse
Affiliation(s)
- Thomas Failing
- Department for Radiotherapy and Radiooncology, University Medical Center Göttingen, Göttingen 37075, Germany; Institute of Medical Physics and Radiation Protection (IMPS), University of Applied Sciences, Gießen 35390, Germany.
| | | | - Frank W Hensley
- Department for Radiotherapy and Radiooncology, University Medical Center Heidelberg, Heidelberg 69120, Germany
| | - Boris Keil
- Institute of Medical Physics and Radiation Protection (IMPS), University of Applied Sciences, Gießen 35390, Germany; Diagnostic and Interventional Radiology, Philipps-University Marburg, Marburg 35043, Germany
| | - Klemens Zink
- Institute of Medical Physics and Radiation Protection (IMPS), University of Applied Sciences, Gießen 35390, Germany; Department for Radiotherapy and Radiooncology, University Medical Center Giessen-Marburg, Marburg 35043, Germany; Marburg Iontherapy Center (MIT), Marburg 35043, Germany
| |
Collapse
|
3
|
Hartmann GH, Andreo P, Kapsch RP, Zink K. Cema-based formalism for the determination of absorbed dose for high-energy photon beams. Med Phys 2021; 48:7461-7475. [PMID: 34613620 DOI: 10.1002/mp.15266] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Revised: 07/26/2021] [Accepted: 09/23/2021] [Indexed: 11/08/2022] Open
Abstract
PURPOSE Determination of absorbed dose is well established in many dosimetry protocols and considered to be highly reliable using ionization chambers under reference conditions. If dosimetry is performed under other conditions or using other detectors, however, open questions still remain. Such questions frequently refer to appropriate correction factors. A converted energy per mass (cema)-based approach to formulate such correction factors offers a good understanding of the specific response of a detector for dosimetry under various measuring conditions and thus an estimate of pros and cons of its application. METHODS Determination of absorbed dose requires the knowledge of the beam quality correction factor kQ,Qo , where Q denotes the quality of a user beam and Qo is the quality of the radiation used for calibration. In modern Monte Carlo (MC)-based methods, kQ,Qo is directly derived from the MC-calculated dose conversion factor, which is the ratio between the absorbed dose at a point of interest in water and the mean absorbed dose in the sensitive volume of an ion chamber. In this work, absorbed dose is approximated by the fundamental quantity cema. This approximation allows the dose conversion factor to be substituted by the cema conversion factor. Subsequently, this factor is decomposed into a product of cema ratios. They are identified as the stopping power ratio water to the material in the sensitive detector volume, and as the correction factor for the fluence perturbation of the secondary charged particles in the detector cavity caused by the presence of the detector. This correction factor is further decomposed with respect to the perturbation caused by the detector cavity and that caused by external detector properties. The cema-based formalism was subsequently tested by MC calculations of the spectral fluence of the secondary charged particles (electrons and positrons) under various conditions. RESULTS MC calculations demonstrate that considerable fluence perturbation may occur particularly under non-reference conditions. Cema-based correction factors to be applied in a 6-MV beam were obtained for a number of ionization chambers and for three solid-state detectors. Feasibility was shown at field sizes of 4 × 4 and 0.5 cm × 0.5 cm. Values of the cema ratios resulting from the decomposition of the dose conversion factor can be well correlated with detector response. Under the small field conditions, the internal fluence correction factor of ionization chambers is considerably dependent on volume averaging and thus on the shape and size of the cavity volume. CONCLUSIONS The cema approach is particularly useful at non-reference conditions including when solid-state detectors are used. Perturbation correction factors can be expressed and evaluated by cema ratios in a comprehensive manner. The cema approach can serve to understand the specific response of a detector for dosimetry to be dependent on (a) radiation quality, (b) detector properties, and (c) electron fluence changes caused by the detector. This understanding may also help to decide which detector is best suited for a specific measurement situation.
Collapse
Affiliation(s)
- Günther H Hartmann
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Pedro Andreo
- Department of Medical Radiation Physics and Nuclear Medicine, Karolinska University Hospital, Stockholm, Sweden.,Department of Oncology-Pathology, Karolinska Institutet, Stockholm, Sweden
| | | | - Klemens Zink
- Institut fuer Medizinische Physik und Strahlenschutz (IMPS), University of Applied Sciences, Giessen, Giessen, Germany.,Department for Radiotherapy and Radiooncology, University Medical Center Giessen-Marburg, Marburg, Germany.,Marburg Iontherapy Center (MIT), Marburg, Germany
| |
Collapse
|
4
|
Rossi G, Gainey M, Kollefrath M, Hofmann E, Baltas D. Suitability of the microDiamond detector for experimental determination of the anisotropy function of High Dose Rate 192 Ir brachytherapy sources. Med Phys 2020; 47:5838-5851. [PMID: 32970875 DOI: 10.1002/mp.14488] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2019] [Revised: 09/07/2020] [Accepted: 09/13/2020] [Indexed: 11/06/2022] Open
Abstract
PURPOSE To investigate the suitability of the microDiamond detector (mDD) type 60019 (PTW-Freiburg, Germany) to measure the anisotropy function F(r,θ) of High Dose Rate (HDR) 192 Ir brachytherapy sources. METHODS The HDR 192 Ir brachytherapy source, model mHDR-v2r (Elekta AB, Sweden), was placed inside a water tank within a 4F plastic needle. Four mDDs (mDD1, mDD2, mDD3, and mDD4) were investigated. Each mDD was placed laterally with respect to the source, and measurements were performed at radial distances r = 1 cm, 3 and 5 cm, and polar angles θ from 0° to 168°. The Monte Carlo (MC) system EGSnrc was used to simulate the measurements and to calculate phantom effect, energy dependence and volume-averaging correction factors. F(r,θ) was determined according to TG-43 formalism from the detector reading corrected with the MC-based factors and compared to the consensus anisotropy function CON F(r,θ). RESULTS At 1 cm, the differences between measurements and MC simulations ranged from -0.8% to +0.8% for θ = 0° and from -2.1% to + 2.3% for θ ≠ 0°. At 3 and 5 cm, the differences ranged from +1.4% to +3.9% for θ = 0°, and from -0.4% to +2.9% for θ ≠ 0°. All differences were within the uncertainties (k = 2). At small angles, the phantom effect correction was up to -1.9%. This effect was mainly caused by the air between source and needle tip. The energy correction was angle-independent everywhere. For small angles at 1 cm, the volume-averaging correction was up to -2.9% and became less important for larger angles and distances. The differences of the measured F(r,θ) corrected with the MC-based factors to CON F(r,θ) ranged from -1.0% to +3.4% for mDD1, -2.2% to +4.2% for mDD2, -2.5% to +4.0% for mDD3, and -2.6% to +3.4% for mDD4. All differences were within the uncertainties (k = 2) except one at (3 cm, 0°). For all the mDDs, F(r,0°) was always higher than CON F(r,0°), with average differences of +3.1% (1 cm), +3.6% (3 cm), and +1.9% (5 cm). The inter-detector variability was within 2.9% (1 cm), 1.8% (3 cm), and 3.4% (5 cm). CONCLUSIONS A reproducible method and experimental setup were presented for measuring and validating F(r,θ) of an HDR 192 Ir brachytherapy source in a water phantom using the mDD. The phantom effect and the volume-averaging need to be taken into account, especially for the smaller distances and angles. Good agreement to CON F(r,θ) was obtained. The discrepancies at (1 cm, 0°), accurately predicted by the MC results, may suggest a reconsideration of CON F(r,θ), at least for this position. The slight overestimations at (3 cm,0°) and (5 cm,0°), both in comparison to CON F(r,θ) and MC results, may be due to an underestimation of the air volume between source and needle tip, dark current, intrinsic over-response of the mDDs, or radiation-induced charge imbalance in the detector's components. The results indicate that the mDD is a valuable tool for measurements with HDR 192 Ir brachytherapy sources and support its employment for the determination and validation of TG-43 parameters of such sources.
Collapse
Affiliation(s)
- Giulio Rossi
- Division of Medical Physics, Department of Radiation Oncology, Medical Center, Faculty of Medicine, University of Freiburg, German Cancer Consortium (DKTK) Partner Site Freiburg, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Mark Gainey
- Division of Medical Physics, Department of Radiation Oncology, Medical Center, Faculty of Medicine, University of Freiburg, German Cancer Consortium (DKTK) Partner Site Freiburg, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Michael Kollefrath
- Division of Medical Physics, Department of Radiation Oncology, Medical Center, Faculty of Medicine, University of Freiburg, German Cancer Consortium (DKTK) Partner Site Freiburg, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Elena Hofmann
- Division of Medical Physics, Department of Radiation Oncology, Medical Center, Faculty of Medicine, University of Freiburg, German Cancer Consortium (DKTK) Partner Site Freiburg, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Dimos Baltas
- Division of Medical Physics, Department of Radiation Oncology, Medical Center, Faculty of Medicine, University of Freiburg, German Cancer Consortium (DKTK) Partner Site Freiburg, German Cancer Research Center (DKFZ), Heidelberg, Germany
| |
Collapse
|
5
|
Georgiou G, Kumar S, Würfel JU, Underwood TSA, Thompson JM, Hill MA, Rowbottom CG, Fenwick JD. Density compensated diodes for small field dosimetry: comprehensive testing and implications for design. Phys Med Biol 2020; 65:155011. [PMID: 32392539 DOI: 10.1088/1361-6560/ab91d9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
PURPOSE In small megavoltage photon fields, the accuracies of an unmodified PTW 60017-type diode dosimeter and six diodes modified by adding airgaps of thickness 0.6-1.6 mm and diameter 3.6 mm have been comprehensively characterized experimentally and computationally. The optimally thick airgap for density compensation was determined, and detectors were micro-CT imaged to investigate differences between experimentally measured radiation responses and those predicted computationally. METHODS Detectors were tested on- and off-axis, at 5 and 15 cm depths in 6 and 15 MV fields ≥ 0.5 × 0.5 cm2. Computational studies were carried out using the EGSnrc/BEAMnrc Monte Carlo radiation transport code. Experimentally, radiation was delivered using a Varian TrueBeam linac and doses absorbed by water were measured using Gafchromic EBT3 film and ionization chambers, and compared with diode readings. Detector response was characterized via the [Formula: see text] formalism, choosing a 4 × 4 cm2 reference field. RESULTS For the unmodified 60017 diode, the maximum error in small field doses obtained from diode readings uncorrected by [Formula: see text] factors was determined as 11.9% computationally at +0.25 mm off-axis and 5 cm depth in a 15 MV 0.5 × 0.5 cm2 field, and 11.7% experimentally at -0.30 mm off-axis and 5 cm depth in the same field. A detector modified to include a 1.6 mm thick airgap performed best, with maximum computationally and experimentally determined errors of 2.2% and 4.1%. The 1.6 mm airgap deepened the modified dosimeter's effective point of measurement by 0.5 mm. For some detectors significant differences existed between responses in small fields determined computationally and experimentally, micro-CT imaging indicating that these differences were due to within-tolerance variations in the thickness of an epoxy resin layer. CONCLUSIONS The dosimetric performance of a 60017 diode detector was comprehensively improved throughout 6 and 15 MV small photon fields via density compensation. For this approach to work well with good detector-to-detector reproducibility, tolerances on dense component dimensions should be reduced to limit associated variations of response in small fields, or these components should be modified to have more water-like densities.
Collapse
Affiliation(s)
- Georgios Georgiou
- Department of Molecular and Clinical Cancer Medicine, Institute of Translational Medicine, University of Liverpool, The Sherrington Building, Ashton Street, Liverpool L69 3BX, United Kingdom. Department of Physics, Clatterbridge Cancer Centre, Clatterbridge Road, Wirral CH63 4JY, United Kingdom. Department of Physics, University of Liverpool, Oliver Lodge Laboratory, Oxford Street, Liverpool L69 7ZE, United Kingdom
| | | | | | | | | | | | | | | |
Collapse
|
6
|
Hartmann GH, Hensley F, Kapsch RP, Poppe B, Sauer O, Würfel J, Zink K. [Detector Based Determination of Water Absorbed Dose According to DIN 6800 Teil 1: Suggestion for an Extension of the Fundamental Formalism]. Z Med Phys 2019; 30:24-39. [PMID: 31585786 DOI: 10.1016/j.zemedi.2019.05.001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2019] [Revised: 05/13/2019] [Accepted: 05/20/2019] [Indexed: 10/25/2022]
Abstract
For any detector to be used for the determination of absorbed dose at the point of measurement in water a basic equation is required to convert the reading of the detector into absorbed dose in water. The German DIN 6800 part 1 provides a general formalism for that. A further differentiated formalism applicable to photon dosimetry is suggested in this work. This modified formalism presents the two following still general and at the same time fundamental properties of any dosimetry detector: a) a clear distinction of correction factors with respect to the physical processes involved during the measurement, and b) the fact that the process of energy absorption in the detector is determined by the spectral distribution of the fluence of the secondary charged particles. It is concluded that based on the modified formalism and knowing this spectral distribution within the detector a general method is available with benefits for ionization chambers as well as for any other dosimetry detector and which is applicable at reference as well as non-reference conditions without any preconditions.
Collapse
Affiliation(s)
| | | | - Ralf-Peter Kapsch
- Physikalisch-Technische Bundesanstalt (PTB), Braunschweig, Deutschland
| | - Bjoern Poppe
- Carl von Ossietzky Universität, Oldenburg, Deutschland
| | - Otto Sauer
- Universitätsklinikum, Würzburg, Deutschland
| | - Jan Würfel
- Physikalisch-Technische Werkstätten (PTW), Freiburg, Deutschland
| | - Klemens Zink
- Technische Hochschule Mittelhessen, Gießen, Deutschland
| |
Collapse
|
7
|
Hartmann GH, Zink K. A Monte Carlo study on the PTW 60019 microDiamond detector. Med Phys 2019; 46:5159-5172. [DOI: 10.1002/mp.13721] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2018] [Revised: 07/02/2019] [Accepted: 07/03/2019] [Indexed: 11/12/2022] Open
Affiliation(s)
| | - Klemens Zink
- Institute of Medical Physics and Radiation Protection (IMPS) University of Applied Sciences Giessen 35390Giessen Germany
- Department for Radiotherapy and Radiooncology University Medical Center Giessen‐Marburg 35043Marburg Germany
- Frankfurt Institute for Advanced Studies (FIAS), Goethe‐University 60438Frankfurt Germany
| |
Collapse
|
8
|
Hartmann GH, Andreo P. Fluence calculation methods in Monte Carlo dosimetry simulations. Z Med Phys 2019; 29:239-248. [DOI: 10.1016/j.zemedi.2018.08.003] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2018] [Revised: 07/04/2018] [Accepted: 08/26/2018] [Indexed: 11/25/2022]
|