Search

Background: There is increased interest in in vivo dosimetry for 192Ir brachytherapy (BT) treatments using high atomic number (Z) inorganic scintillators. Their high light output enables construction of small detectors with negligible stem effect and simple readout electronics. Experimental determination of absorbed‐dose energy dependence of detectors relative to water is prevalent, but it can be prone to high detector positioning uncertainties and does not allow for decoupling of absorbed‐dose energy dependence from other factors affecting detector response. Purpose: To investigate which measurement conditions and detector properties could affect their absorbed‐dose energy dependence in BT in vivo dosimetry. Methods: We used a general‐purpose Monte Carlo (MC) code PENELOPE for the characterization of high‐Z inorganic scintillators with the focus on ZnSe (Z¯=32$\bar{Z}=32$) Z. Two other promising media CsI (Z¯=54$\bar{Z}=54$) and Al2O3 (Z¯=11$\bar{Z}=11$) were included for comparison in selected scenarios. We determined absorbed‐dose energy dependence of crystals relative to water under different scatter conditions (calibration phantom 12 × 12 × 30 cm3, characterization phantoms 20 × 20 × 20 cm3, 30 × 30 × 30 cm3, 40 × 40 × 40 cm3, and patient‐like elliptic phantom 40 × 30 × 25 cm3). To mimic irradiation conditions during prostate treatments, we evaluated whether the presence of pelvic bones and calcifications affect ZnSe response. ZnSe detector design influence was also investigated. Results: In contrast to low‐Z organic and medium‐Z inorganic scintillators, ZnSe and CsI media have substantially greater absorbed‐dose energy dependence relative to water. The response was phantom‐size dependent and changed by 11% between limited‐ and full‐scatter conditions for ZnSe, but not for Al2O3. For a given phantom size, a part of the absorbed‐dose energy dependence of ZnSe is caused not due to in‐phantom scatter but due to source anisotropy. Thus, the absorbed‐dose energy dependence of high‐Z scintillators is a function of not only the radial distance but also the polar angle. Pelvic bones did not affect ZnSe response, whereas large and intermediate size calcifications reduced it by 9% and 5%, respectively, when placed midway between the source and the detector. Conclusions: Unlike currently prevalent low‐ and medium‐Z scintillators, high‐Z crystals are sensitive to characterization and in vivo measurement conditions. However, good agreement between MC data for ZnSe in the present study and experimental data for ZnSe:O by Jørgensen et al. (2021) suggests that detector signal is proportional to the average absorbed dose to the detector cavity. This enables an easy correction for non‐TG43‐like scenarios (e.g., patient sizes and calcifications) through MC simulations. Such information should be provided to the clinic by the detector vendors. [ABSTRACT FROM AUTHOR]

Purpose: An important characteristic of radiation dosimetry detectors is their energy response which consists of absorbed‐dose and intrinsic energy responses. The former can be characterized using Monte Carlo (MC) simulations, whereas the latter (i.e., detector signal per absorbed dose to detector) is extracted from experimental data. Such a characterization is especially relevant when detectors are used in nonrelative measurements at a beam quality that differs from the calibration beam quality. Having in mind the possible application of synthetic diamond detectors (microDiamond PTW 60019, Freiburg, Germany) for nonrelative dosimetry of low‐energy brachytherapy (BT) beams, we determined their intrinsic and absorbed‐dose energy responses in 25–250 kV beams relative to a 60Co beam, which is usually the reference beam quality for detector calibration in radiotherapy. Material and Methods: Three microDiamond detectors and, for comparison, two silicon diodes (PTW 60017) were calibrated in terms of air‐kerma free in air in six x‐ray beam qualities (from 25 to 250 kV) and in terms of absorbed dose to water in a 60Co beam at the national metrology laboratory in Sweden. The PENELOPE/penEasy MC radiation transport code was used to calculate the absorbed‐dose energy response of the detectors (modeled based on blueprints) relative to air and water depending on calibration conditions. The MC results were used to extract the relative intrinsic energy response of the detectors from the overall energy response. Measurements using an independent setup with a single ophthalmic BEBIG I25.S16 125I BT seed (effective photon energy of 28 keV) were used as a qualitative check of the extracted intrinsic energy response correction factors. Additionally, the impact of the thickness of the active volume as well as the presence of extra‐cameral components on the absorbed‐dose energy response of a microDiamond detector was studied using MC simulations. Results: The relative intrinsic energy response of the microDiamond detectors was higher by a factor of 2 in 25 and 50 kV beams compared to the 60Co beam. The variation in the relative intrinsic energy response of silicon diodes was within 10% over the investigated photon energy range. The use of relative intrinsic energy response correction factors improved the agreement among the absorbed dose to water values determined using microDiamond detectors and silicon diodes, as well as with the TG‐43 formalism‐based calculations for the 125I seed. MC study of microDiamond detector design features provided a possible explanation for inter‐detector response variation at low‐energy photon beams by differences in the effective thickness of the active volume. Conclusions: MicroDiamond detectors had a non‐negligible variation in the relative intrinsic energy response (factor of 2) which was comparable to that in the absorbed‐dose energy response relative to water at low‐energy photon beams. Silicon diodes, in contrast, had an absorbed‐dose energy dependence on photon energy that varied by a factor of 6, whereas the intrinsic energy dependence on beam quality was within 10%. It is important to decouple these two responses for a full characterization of detector energy response especially when the user and reference beam qualities differ significantly, and MC alone is not enough. [ABSTRACT FROM AUTHOR]

Purpose: To investigate the applicability of output correction factors reported in TRS‐483 on 6‐MV small‐field detector‐reading ratios using four solid‐state detectors. Also, to investigate variations in 6‐MV small‐field output factors (OF) among nominally matched linear accelerators (linacs). Methods: The TRS‐483 Code of Practice (CoP) introduced and provided output correction factors to be applied to measured detector‐reading ratios to obtain OFs for several small‐field detectors. Detector readings for 0.5 cm × 0.5 cm to 8 cm × 8 cm fields were measured and normalized to that of 10 cm × 10 cm field giving the detector‐reading ratios. Three silicon diodes, IBA PFD, IBA EFD (IBA, Schwarzenbruck, Germany), PTW T60017, and one microdiamond, PTW T60019 (PTW, Freiburg, Germany), were used. Output correction factors from the CoP were applied to measured detector‐reading ratios. Measurements were performed on six Clinac and six TrueBeam linacs (Varian Medical Systems, Palo Alto, USA). An investigation of the relationship between the size of small fields and corresponding detector‐reading ratio among the linacs was performed by measuring lateral dose profiles for 0.5 cm × 0.5 cm fields to determine the full width half maximum (FWHM). The relationship between the linacs' focal spot size and the small‐field detector‐reading ratio was investigated by measuring 10 cm × 10 cm lateral dose profiles and determining the penumbra width reflecting the focal spot size. Measurement geometry was as follows: gantry angle = 0°, collimator angle = 0°, source‐to surface distance (SSD) = 90 cm, and depth in water = 10 cm. Results: For a given linac and 0.5 cm × 0.5 cm field, the deviations in detector‐reading ratios among the detectors were 9%–15% for the Clinacs and 4%–5% for the TrueBeams. Use of output correction factors reduced these deviations to 6%–12% and 3%–4%, respectively. For field sizes equal to or larger than 0.8 cm × 0.8 cm, the deviations were corrected to 1% using output correction factors for both Clinacs and TrueBeams. For a given detector and 0.5 cm × 0.5 cm field, the deviations in detector‐reading ratios among the linacs were 11%–17% for the Clinacs and 5–6% for the TrueBeams. For 1 cm × 1 cm the deviations were 1%–2% for Clinacs and 1% for TrueBeams. For field sizes larger than 1 cm × 1 cm the deviations were within 1% for both Clinacs and TrueBeams. No relationship between FWHMs and detector‐reading ratios for 0.5 cm × 0.5 cm was observed. For Clinacs, larger 10 cm × 10 cm penumbra width yielded lower 0.5 cm × 0.5 cm detector‐reading ratio indicating an effect of the focal spot size. For TrueBeams, the spread of penumbra widths was lower compared to Clinacs and no similar relationship was observed. Conclusions: Output correction factors from the TRS‐483 CoP are not sufficient for accurate determination of OF for 0.5 cm × 0.5 cm fields but are applicable for 0.8 cm × 0.8 cm to 8 cm × 8 cm fields. Nominally matched Clinacs and TrueBeams show large differences in detector‐reading ratios for fields smaller than 1 cm × 1 cm. [ABSTRACT FROM AUTHOR]

Purpose: High dose–rate brachytherapy is a method of radiotherapy for cancer treatment in which the radiation source is placed within the body. In addition to give a high enough dose to a tumor, it is also important to spare nearby healthy organs [organs at risk (OAR)]. Dose plans are commonly evaluated using the so‐called dosimetric indices; for the tumor, the portion of the structure that receives a sufficiently high dose is calculated, while for OAR it is instead the portion of the structure that receives a sufficiently low dose that is of interest. Models that include dosimetric indices are referred to as dose–volume models (DVMs) and have received much interest recently. Such models do not take the dose to the coldest (least irradiated) volume of the tumor into account, which is a distinct weakness since research indicates that the treatment effect can be largely impaired by tumor underdosage even to small volumes. Therefore, our aim is to extend a DVM to also consider the dose to the coldest volume. Methods: An improved DVM for dose planning is proposed. In addition to optimizing with respect to dosimetric indices, this model also takes mean dose to the coldest volume of the tumor into account. Results: Our extended model has been evaluated against a standard DVM in ten prostate geometries. Our results show that the dose to the coldest volume could be increased, while also computing times for the dose planning were improved. Conclusion: While the proposed model yields dose plans similar to other models in most aspects, it fulfils its purpose of increasing the dose to cold tumor volumes. An additional benefit is shorter solution times, and especially for clinically relevant times (of minutes) we show major improvements in tumour dosimetric indices. [ABSTRACT FROM AUTHOR]

Purpose: Experimental dosimetry of high‐dose‐rate (HDR) 192Ir brachytherapy (BT) sources is complicated due to high dose and dose‐rate gradients, and softening of photon energy spectrum with depth. A single crystal synthetic diamond detector microDiamond (PTW 60019, Freiburg, Germany) has a small active volume, high sensitivity, direct readout, and nearly water‐equivalent active volume. The purpose of this study was to evaluate the suitability of microDiamond detectors for the determination of absorbed dose to water around HDR 192Ir BT sources. Three microDiamond detectors were used, allowing for the comparison of their properties. Methods: In‐phantom measurements were performed using microSelectron and VariSource iX HDR 192Ir BT treatment units. Their treatment planning systems (TPSs), Oncentra (v. 4.3) and BrachyVision (v. 13.6), respectively, were used to create irradiation plans for a cubic PMMA phantom with the microDiamond positioned at one of three source‐to‐detector distances (SDDs) (1.5, 2.5, and 5.5 cm) at a time. The source was stepped in increments of 0.5 cm over a total length of 6 cm to yield absorbed dose of 2 Gy at the nominal reference‐point of the detector. Detectors were calibrated in 60Co beam in terms of absorbed dose to water, and Monte Carlo (MC) calculated beam quality correction factors were applied to account for absorbed‐dose energy dependence. Phantom correction factors were applied to account for differences in dimensions between the measurement phantom and a water phantom used for absorbed dose calculations made with a TPS. The same measurements were made with all three of the detectors. Additionally, dose‐rate dependence and stability of the detectors were evaluated in 60Co beam. Results: The percentage differences between experimentally determined and TPS‐calculated absorbed doses to water were from −1.3% to +2.9%. The values agreed to within experimental uncertainties, which were from 1.9% to 4.3% (k = 2) depending on the detector, SDD and treatment delivery unit. No dose‐rate or intrinsic energy dependence corrections were applied. All microDiamonds were comparable in terms of preirradiation dose, stability of the readings and energy response, and showed a good agreement. Conclusions: The results indicate that the microDiamond is potentially suitable for the determination of absorbed dose to water around HDR 192Ir BT sources and may be used for independent verification of TPS's calculations, as well as for QA measurements of HDR 192Ir BT treatment delivery units at clinical sites. [ABSTRACT FROM AUTHOR]

Introduction In photon brachytherapy (BT), experimental dosimetry is needed to verify treatment plans if planning algorithms neglect varying attenuation, absorption or scattering conditions. The detector’s response is energy dependent, including the detector material to water dose ratio and the intrinsic mechanisms. The local mean photon energy E ¯ ( r ) must be known or another equivalent energy quality parameter used. We propose the brachytherapy photon radiation quality index Q BT ( E ¯ ) , to characterize the photon radiation quality in view of measurements of distributions of the absorbed dose to water, D w , around BT sources. Materials and methods While the external photon beam radiotherapy (EBRT) radiation quality index Q EBRT ( E ¯ ) = TPR 10 20 ( E ¯ ) is not applicable to BT, the authors have applied a novel energy dependent parameter, called brachytherapy photon radiation quality index , defined as Q BT ( E ¯ ) = D prim ( r = 2 cm, θ 0 = 90 ° ) / D prim ( r 0 = 1 cm, θ 0 = 90 ° ) , utilizing precise primary absorbed dose data, D prim , from source reference databases, without additional MC-calculations. Results and discussion For BT photon sources used clinically, Q BT ( E ¯ ) enables to determine the effective mean linear attenuation coefficient μ ¯ ( E ) and thus the effective energy of the primary photons E prim eff ( r 0 , θ 0 ) at the TG-43 reference position P ref ( r 0 = 1 cm, θ 0 = 90 ° ) , being close to the mean total photon energy E ¯ tot ( r 0 , θ 0 ) . If one has calibrated detectors, published E ¯ tot ( r ) and the BT radiation quality correction factor k Q , Q 0 BT ( E ¯ , r , θ ) for different BT radiation qualities Q and Q 0 , the detector’s response can be determined and D w ( r , θ ) measured in the vicinity of BT photon sources. Conclusions This novel brachytherapy photon radiation quality index Q BT characterizes sufficiently accurate and precise the primary photon‘s penetration probability and scattering potential. [ABSTRACT FROM AUTHOR]

Purpose: In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) 192Ir source and a virtual water phantom were designed, which can be imported into a TPS. Methods: A hypothetical, generic HDR 192Ir source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic 192Ir source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra® Brachy with advanced collapsed-cone engine (ACE) and BrachyVision ACUROS™ ]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and PENELOPE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)³ voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR 192Ir source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods. Results: TG-43 and PSS datasets were generated for the generic source, the PSS data for use with theACE algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109±0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ACE at clinically relevant distances. Conclusions: A hypothetical, generic HDR 192Ir source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs. [ABSTRACT FROM AUTHOR]

Abstract: Purpose: According to the American Association of Physicists in Medicine Task Group No. 43 (TG-43) formalism used for dose calculation in brachytherapy treatment planning systems, the absolute level of absorbed dose is determined through coupling with the measurable quantity air-kerma strength or the numerically equal reference air-kerma rate (RAKR). Traceability to established standards is important for accurate dosimetry in laying the ground for reliable comparisons of results and safety in adoption of new treatment protocols. The purpose of this work was to compare the source strength for a high-dose rate (HDR) 192Ir source as measured using equipment traceable to different standard laboratories in Europe and the United States. Methods and Materials: Source strength was determined for one HDR 192Ir source using four independent systems, all with traceability to different primary or interim standards in the United States and Europe. Results: The measured HDR 192Ir source strengths varied by 0.8% and differed on average from the vendor value by 0.3%. Measurements with the well chambers were 0.5% ± 0.1% higher than the vendor-provided source strength. Measurements with the Farmer chamber were 0.7% lower than the average well chamber results and 0.2% lower than the vendor-provided source strength. All of these results were less than the reported source calibration uncertainties (k=2) of each measurement system. Conclusions: In view of the uncertainties in ion chamber calibration factors, the maximum difference in source strength found in this study is small and confirms the consistency between calibration standards in use for HDR 192Ir brachytherapy. [Copyright &y& Elsevier]

The charge of Task Group 186 (TG-186) is to provide guidance for early adopters of model-based dose calculation algorithms (MBDCAs) for brachytherapy (BT) dose calculations to ensure practice uniformity. Contrary to external beam radiotherapy, heterogeneity correction algorithms have only recently been made available to the BT community. Yet, BT dose calculation accuracy is highly dependent on scatter conditions and photoelectric effect cross-sections relative to water. In specific situations, differences between the current water-based BT dose calculation formalism (TG-43) and MBDCAs can lead to differences in calculated doses exceeding a factor of 10. MBDCAs raise three major issues that are not addressed by current guidance documents: (1) MBDCA calculated doses are sensitive to the dose specification medium, resulting in energy-dependent differences between dose calculated to water in a homogeneous water geometry (TG-43), dose calculated to the local medium in the heterogeneous medium, and the intermediate scenario of dose calculated to a small volume of water in the heterogeneous medium. (2) MBDCA doses are sensitive to voxel-by-voxel interaction cross sections. Neither conventional single-energy CT nor ICRU/ICRP tissue composition compilations provide useful guidance for the task of assigning interaction cross sections to each voxel. (3) Since each patient-source-applicator combination is unique, having reference data for each possible combination to benchmark MBDCAs is an impractical strategy. Hence, a new commissioning process is required. TG-186 addresses in detail the above issues through the literature review and provides explicit recommendations based on the current state of knowledge. TG-43-based dose prescription and dose calculation remain in effect, with MBDCA dose reporting performed in parallel when available. In using MBDCAs, it is recommended that the radiation transport should be performed in the heterogeneous medium and, at minimum, the dose to the local medium be reported along with the TG-43 calculated doses. Assignments of voxel-by-voxel cross sections represent a particular challenge. Electron density information is readily extracted from CT imaging, but cannot be used to distinguish between different materials having the same density. Therefore, a recommendation is made to use a number of standardized materials to maintain uniformity across institutions. Sensitivity analysis shows that this recommendation offers increased accuracy over TG-43. MBDCA commissioning will share commonalities with current TG-43-based systems, but in addition there will be algorithm-specific tasks. Two levels of commissioning are recommended: reproducing TG-43 dose parameters and testing the advanced capabilities of MBDCAs. For validation of heterogeneity and scatter conditions, MBDCAs should mimic the 3D dose distributions from reference virtual geometries. Potential changes in BT dose prescriptions and MBDCA limitations are discussed. When data required for full MBDCA implementation are insufficient, interim recommendations are made and potential areas of research are identified. Application of TG-186 guidance should retain practice uniformity in transitioning from the TG-43 to the MBDCA approach. [ABSTRACT FROM AUTHOR]

Purpose: Experimental radiation dosimetry with thermoluminescent dosimeters (TLDs), calibrated in a 60Co or megavoltage (MV) photon beam, is recommended by AAPM TG-43U1for verification of Monte Carlo calculated absorbed doses around brachytherapy sources. However, it has been shown by Carlsson Tedgren et al. [Med. Phys. 38, 5539-5550 (2011)] that for TLDs of LiF:Mg,Ti, detector response was 4% higher in a 137Cs beam than in a 60Co one. The aim of this work was to investigate if similar over-response exists when measuring absorbed dose to water around 192Ir sources, using LiF:Mg,Ti dosimeters calibrated in a 6 MV photon beam. Methods: LiF dosimeters were calibrated to measure absorbed dose to water in a 6 MV photon beam and used to measure absorbed dose to water at distances of 3, 5, and 7 cm from a clinical high dose rate (HDR) 192Ir source in a polymethylmethacrylate (PMMA) phantom. Measured values were compared to values of absorbed dose to water calculated using a treatment planning system (TPS) including corrections for the difference in energy absorption properties between calibration quality and the quality in the users' 192Ir beam and for the use of a PMMA phantom instead of the water phantom underlying dose calculations in the TPS. Results: Measured absorbed doses to water around the 192Ir source were overestimated by 5% compared to those calculated by the TPS. Corresponding absorbed doses to water measured in a previous work with lithium formate electron paramagnetic resonance (EPR) dosimeters by Antonovic et al. [Med. Phys. 36, 2236-2247 (2009)], using the same irradiation setup and calibration procedure as in this work, were 2% lower than those calculated by the TPS. The results obtained in the measurements in this work and those obtained using the EPR lithium formate dosimeters were, within the expanded (k = 2) uncertainty, in agreement with the values derived by the TPS. The discrepancy between the results using LiF:Mg,Ti TLDs and the EPR lithium formate dosimeters was, however, statistically significant and in agreement with the difference in relative detector responses found for the two detector systems by Carlsson Tedgren et al. [Med. Phys. 38, 5539-5550 (2011)] and by Adolfsson et al. [Med. Phys. 37, 4946-4959 (2010)]. Conclusions: When calibrated in 60Co or MV photon beams, correction for the linear energy transfer (LET) dependence of LiF:Mg,Ti detector response will be needed as to measure absorbed doses to water in a 192Ir beam with highest accuracy. Such corrections will depend on the manufacturing process (MTS-N Poland or Harshaw TLD-100) and details of the annealing and read-out schemes used. [ABSTRACT FROM AUTHOR]

Purpose: Dose plans generated with optimization models hitherto used in high-dose-rate (HDR) brachytherapy have shown a tendency to yield longer dwell times than manually optimized plans. Concern has been raised for the corresponding undesired hot spots, and various methods to mitigate these have been developed. The hypotheses upon this work is based are (a) that one cause for the long dwell times is the use of objective functions comprising simple linear penalties and (b) that alternative penalties, as these are piecewise linear, would lead to reduced length of individual dwell times. Methods: The characteristics of the linear penalties and the piecewise linear penalties are analyzed mathematically. Experimental comparisons between the two types of penalties are carried out retrospectively for a set of prostate cancer patients. Results: When the two types of penalties are compared, significant changes can be seen in the dwell times, while most dose-volume parameters do not differ significantly. On average, total dwell times were reduced by 4.2%, with a reduction of maximum dwell times by 25%, when the alternative penalties were used. Conclusions: The use of linear penalties in optimization models for HDR brachytherapy is one cause for the undesired long dwell times that arise in mathematically optimized plans. By introducing alternative penalties, a significant reduction in dwell times can be achieved for HDR brachytherapy dose plans. Although various measures for mitigating the long dwell times are already available, the observation that linear penalties contribute to their appearance is of fundamental interest. [ABSTRACT FROM AUTHOR]

A dosimetry system using lithium formate monohydrate (HCO2Li·H2O) as detector material and electron paramagnetic resonance (EPR) spectroscopy for readout has been used to measure absorbed dose distributions around clinical 192Ir sources. Cylindrical tablets with diameter of 4.5 mm, height of 4.8 mm, and density of 1.26 g/cm3 were manufactured. Homogeneity test and calibration of the dosimeters were performed in a 6 MV photon beam. 192Ir irradiations were performed in a PMMA phantom using two different source models, the GammaMed Plus HDR and the microSelectron PDR-v1 model. Measured absorbed doses to water in the PMMA phantom were converted to the corresponding absorbed doses to water in water phantoms of dimensions used by the treatment planning systems (TPSs) using correction factors explicitly derived for this experiment. Experimentally determined absorbed doses agreed with the absorbed doses to water calculated by the TPS to within ±2.9%. Relative standard uncertainties in the experimentally determined absorbed doses were estimated to be within the range of 1.7%–1.3% depending on the radial distance from the source, the type of source (HDR or PDR), and the particular absorbed doses used. This work shows that a lithium formate dosimetry system is well suited for measurements of absorbed dose to water around clinical HDR and PDR 192Ir sources. Being less energy dependent than the commonly used thermoluminescent lithium fluoride (LiF) dosimeters, lithium formate monohydrate dosimeters are well suited to measure absorbed doses in situations where the energy dependence cannot easily be accounted for such as in multiple-source irradiations to verify treatment plans. Their wide dynamic range and linear dose response over the dose interval of 0.2–1000 Gy make them suitable for measurements on sources of the strengths used in clinical applications. The dosimeter size needs, however, to be reduced for application to single-source dosimetry. [ABSTRACT FROM AUTHOR]

Abstract: Background and purpose: To investigate the status of source strength determination in terms of reference air kerma rate (RAKR) for HDR and PDR 192Ir brachytherapy in Sweden. Materials and methods: RAKR was determined in each of the 14 Swedish afterloaders, using calibrated equipment from the Swedish Secondary Standard Dosimetry Laboratory. Results: Values of RAKR from the external audit, the hospitals and vendors agreed within the uncertainty limits guaranteed by the vendors. Conclusions: The accuracy in RAKR determination has increased over the last years as a result of increased availability of interpolation standards for HDR 192Ir and the increased use of robust well-type ion chambers designed for brachytherapy. It is recommended to establish a ratio between the RAKR value from own measurements at the hospital and that of the vendor since such a ratio embeds constant systematic differences due to e.g. varying traceability and therefore has the potential of being less uncertain than the RAKR alone. Traceability to primary standards for HDR 192Ir sources will in the future significantly decrease the uncertainty in RAKR of 192Ir brachytherapy. [Copyright &y& Elsevier]

Background Purpose Methods Results Conclusions Model‐based dose calculation algorithms (MBDCA), such as the Advanced Collapsed cone Engine (ACE) in Oncentra Brachy® can be used to overcome the limitations of the TG‐43 formalism. ACE is a point kernel superposition algorithm that calculates the total dose separated into primary, first‐scatter, and multiple‐scatter dose. Albeit ACE yields accurate results under most circumstances, several studies have reported underestimations of the dose to cortical bone. These underestimations are likely caused by approximations in the handling of multiple‐scatter dose for non‐water media. Such would result in noticeable deviations where the multiple‐scatter is a considerable fraction of the total dose, that is, at greater distances from the source.To improve and test the accuracy of the multiple‐scatter dose component in the ACE algorithm to remedy its inaccuracy for non‐water geometries.A careful analysis of the transport and absorption of the multiple‐scatter energy fluence revealed an inconsistency in the scaling of energy absorption ratios for non‐water media of the multiple‐scatter kernel. We implemented an updated algorithm version, ACEcorr, and tested it for three different geometries. All had a single 192Ir‐source at the center of a cubic water phantom with a box‐shaped heterogeneity of either cortical bone or air, positioned at different distances from the source. Dose distributions for the three cases were calculated with ACE and ACEcorr and compared to Monte Carlo simulations, using the percentage dose difference ratio as figure‐of‐merit. All dose calculation methods scored separately the dose deposited by primary, first‐scattered, and multiple‐scattered photons.The accuracy of the updated algorithm ACEcorr was superior to ACE. In the cortical bone heterogeneity, the mean percentage dose difference ratio for the total dose improved from −11.7%$ - 11.7{\mathrm{\% }}$ to −2.2%$ - 2.2{\mathrm{\% }}$ (in the worst case) by our update. Less impact was seen in the air heterogeneity, where both ACE and ACEcorr deviated less than 2% from the Monte Carlo results. The algorithm update mainly concerns the multiple‐scattered dose component, but an accompanying data processing update also had a small effect (≤$ \le $0.5% difference) on the primary and first‐scattered dose. The calculation times were not affected.The updates to ACE improved the accuracy of multiple‐scatter dose calculation for non‐water media, without increasing calculation times. [ABSTRACT FROM AUTHOR]

• Most Radiotherapy Departments in Europe have HDR-PDR brachytherapy equipment with Ir-192 or Co-60. • The present GEC-ESTRO ACROP Recommendations provide guidance on the calibration of such sources. • Practical aspects and issues not specifically accounted for in well-accepted societal recommendations are included. • The aim is to provide a European-wide standard in HDR-PDR BT source calibration at the hospital level. The vast majority of radiotherapy departments in Europe using brachytherapy (BT) perform temporary implants of high- or pulsed-dose rate (HDR-PDR) sources with photon energies higher than 50 keV. Such techniques are successfully applied to diverse pathologies and clinical scenarios. These recommendations are the result of Working Package 21 (WP-21) initiated within the BRAchytherapy PHYsics Quality Assurance System (BRAPHYQS) GEC-ESTRO working group with a focus on HDR-PDR source calibration. They provide guidance on the calibration of such sources, including practical aspects and issues not specifically accounted for in well-accepted societal recommendations, complementing the BRAPHYQS WP-18 Report dedicated to low energy BT photon emitting sources (seeds). The aim of this report is to provide a European-wide standard in HDR-PDR BT source calibration at the hospital level to maintain high quality patient treatments. [ABSTRACT FROM AUTHOR]

Background The Stockholm Hemangioma Cohort is important for evaluation of late effects after exposure to ionizing radiation during childhood. Dose estimates in this cohort were based on both measurements and calculations using an old treatment planning system. Methods We compare previously published and calculated dose estimates with new ones, obtained by Monte Carlo simulations, which mimic the hemangioma treatments with 226 Ra needles and tubes. The distances between the 226 Ra sources and the thyroid and breasts, respectively, were reassessed. Result The Monte Carlo calculations showed significantly lower dose values than those obtained earlier. The differences depended both on the modeling of the sources and on further individualized distances from the sources. The mean value of the new calculated doses was 25% of the old breast doses and 46% of the old thyroid doses. Conclusion New dosimetry for hemangioma treatments gives significantly lower organ doses for the few cases receiving the highest absorbed dose values. This implies that radiation risk estimates will increase and have to be recalculated. For retrospective studies it is now possible to calculate organ doses from radium treatments using modern treatment planning systems by modeling the source geometry carefully and apply the TG-43 formalism. It is important to be aware of the large uncertainties in calculated absorbed dose values. [ABSTRACT FROM AUTHOR]

Purpose: To develop a prognostic score that correlates to a low, medium, and high incidence of treatment failure after plaque brachytherapy of uveal melanoma (UM)., Methods and Materials: All patients who have received plaque brachytherapy for posterior UM at St. Erik Eye Hospital in Stockholm, Sweden from 1995 through 2019 were included (n = 1636). Treatment failure was defined as tumor recurrence, lack of tumor regression, or any other condition requiring a secondary transpupillary thermotherapy (TTT), plaque brachytherapy, or enucleation. The total sample was randomized into 1 training and 1 validation cohort, and a prognostic score for the risk for treatment failure was developed., Results: In multivariate Cox regression, low visual acuity, tumor distance to the optic disc ≤2 mm, American Joint Committee on Cancer (AJCC) stage, and a tumor apical thickness of >4 (for Ruthenium-106) or >9 mm (for Iodine-125) were independent predictors of treatment failure. No reliable threshold could be identified for tumor diameter or cancer stage. In competing risk analyses of the validation cohort, the cumulative incidence of treatment failure, as well as of secondary enucleation, increased with the prognostic score: In the low, intermediate, and high-risk classes, the 10-year incidence of treatment failure was 19, 28, and 35% and of secondary enucleation 7, 19, and 25 %, respectively., Conclusions: Low visual acuity, American Joint Committee on Cancer stage, tumor thickness, and tumor distance to the optic disc are independent predictors of treatment failure after plaque brachytherapy for UM. A prognostic score was devised that identifies low, medium, and high risk for treatment failure., (© 2022 The Author(s).)

Background: There is increased interest in in vivo dosimetry for 192 Ir brachytherapy (BT) treatments using high atomic number (Z) inorganic scintillators. Their high light output enables construction of small detectors with negligible stem effect and simple readout electronics. Experimental determination of absorbed-dose energy dependence of detectors relative to water is prevalent, but it can be prone to high detector positioning uncertainties and does not allow for decoupling of absorbed-dose energy dependence from other factors affecting detector response ., Purpose: To investigate which measurement conditions and detector properties could affect their absorbed-dose energy dependence in BT in vivo dosimetry., Methods: We used a general-purpose Monte Carlo (MC) code PENELOPE for the characterization of high-Z inorganic scintillators with the focus on ZnSe ( Z ¯ = 32 $\bar{Z}=32$ ) Z. Two other promising media CsI ( Z ¯ = 54 $\bar{Z}=54$ ) and Al 2 O 3 ( Z ¯ = 11 $\bar{Z}=11$ ) were included for comparison in selected scenarios. We determined absorbed-dose energy dependence of crystals relative to water under different scatter conditions (calibration phantom 12 × 12 × 30 cm 3 , characterization phantoms 20 × 20 × 20 cm 3 , 30 × 30 × 30 cm 3 , 40 × 40 × 40 cm 3 , and patient-like elliptic phantom 40 × 30 × 25 cm 3 ). To mimic irradiation conditions during prostate treatments, we evaluated whether the presence of pelvic bones and calcifications affect ZnSe response. ZnSe detector design influence was also investigated., Results: In contrast to low-Z organic and medium-Z inorganic scintillators, ZnSe and CsI media have substantially greater absorbed-dose energy dependence relative to water. The response was phantom-size dependent and changed by 11% between limited- and full-scatter conditions for ZnSe, but not for Al 2 O 3 . For a given phantom size, a part of the absorbed-dose energy dependence of ZnSe is caused not due to in-phantom scatter but due to source anisotropy. Thus, the absorbed-dose energy dependence of high-Z scintillators is a function of not only the radial distance but also the polar angle. Pelvic bones did not affect ZnSe response, whereas large and intermediate size calcifications reduced it by 9% and 5%, respectively, when placed midway between the source and the detector., Conclusions: Unlike currently prevalent low- and medium-Z scintillators, high-Z crystals are sensitive to characterization and in vivo measurement conditions. However, good agreement between MC data for ZnSe in the present study and experimental data for ZnSe:O by Jørgensen et al. (2021) suggests that detector signal is proportional to the average absorbed dose to the detector cavity. This enables an easy correction for non-TG43-like scenarios (e.g., patient sizes and calcifications) through MC simulations. Such information should be provided to the clinic by the detector vendors., (© 2022 American Association of Physicists in Medicine.)

Brachytherapy is a mature treatment modality. The literature is abundant in terms of review articles and comprehensive books on the latest established as well as evolving clinical practices. The intent of this article is to part ways and look beyond the current state-of-the-art and review emerging technologies that are noteworthy and perhaps may drive the future innovations in the field. There are plenty of candidate topics that deserve a deeper look, of course, but with practical limits in this communicative platform, we explore four topics that perhaps is worthwhile to review in detail at this time. First, intensity modulated brachytherapy (IMBT) is reviewed. The IMBT takes advantage of anisotropic radiation profile generated through intelligent high-density shielding designs incorporated onto sources and applicators such to achieve high quality plans. Second, emerging applications of 3D printing (i.e. additive manufacturing) in brachytherapy are reviewed. With the advent of 3D printing, interest in this technology in brachytherapy has been immense and translation swift due to their potential to tailor applicators and treatments customizable to each individual patient. This is followed by, in third, innovations in treatment planning concerning catheter placement and dwell times where new modelling approaches, solution algorithms, and technological advances are reviewed. And, fourth and lastly, applications of a new machine learning technique, called deep learning, which has the potential to improve and automate all aspects of brachytherapy workflow, are reviewed. We do not expect that all ideas and innovations reviewed in this article will ultimately reach clinic but, nonetheless, this review provides a decent glimpse of what is to come. It would be exciting to monitor as IMBT, 3D printing, novel optimization algorithms, and deep learning technologies evolve over time and translate into pilot testing and sensibly phased clinical trials, and ultimately make a difference for cancer patients. Today's fancy is tomorrow's reality. The future is bright for brachytherapy., (© 2021 Institute of Physics and Engineering in Medicine.)

Cerium oxide nanoparticles with integrated gadolinium have been proved to be useful as contrast agents in magnetic resonance imaging. Of question is their performance in dual-energy computed tomography. The aims of this work are to determine (1) the relation between the computed tomography number and the concentration of the I, Gd or Ce contrast agent and (2) under what conditions it is possible to resolve the type of contrast agent. Hounsfield values of iodoacetic acid, gadolinium acetate and cerium acetate dissolved in water at molar concentrations of 10, 50 and 100 mM were measured in a water phantom using the Siemens SOMATOM Definition Force scanner; gadolinium- and cerium acetate were used as substitutes for the gadolinium-integrated cerium oxide nanoparticles. The relation between the molar concentration of the I, Gd or Ce contrast agent and the Hounsfield value was linear. Concentrations had to be sufficiently high to resolve the contrast agents., (© The Author(s) 2021. Published by Oxford University Press.)

Dual-energy computed tomography (CT) can be used in radiotherapy treatment planning for the calculation of absorbed dose distributions. The aim of this work is to evaluate whether there is room for improvement in the accuracy of the Monoenergetic Plus algorithm by Siemens Healthineers. A Siemens SOMATOM Force scanner was used to scan a cylindrical polymethyl methacrylate phantom with four rod-inserts made of different materials. Images were reconstructed using ADMIRE and processed with Monoenergetic Plus. The resulting CT numbers were compared with tabulated values and values simulated by the proof-of-a-concept algorithm DIRA developed by the authors. Both the Monoenergetic Plus and DIRA algorithms performed well; the accuracy of attenuation coefficients was better than about ±1% at the energy of 70 keV. Compared with DIRA, the worse performance of Monoenergetic Plus was caused by its (i) two-material decomposition to iodine and water and (ii) imperfect suppression of the beam hardening artifact in ADMIRE., (© The Author(s) 2021. Published by Oxford University Press.)

The choice of the material base to which the material decomposition is performed in dual-energy computed tomography may affect the quality of reconstructed images. The aim of this work is to investigate how the commonly used bases (water, bone), (water, iodine) and (photoelectric effect, Compton scattering) affect the reconstructed linear attenuation coefficient in the case of the Alvarez-Macovski method. The performance of this method is also compared with the performance of the Dual-energy Iterative Reconstruction Algorithm (DIRA). In both cases, the study is performed using simulations. The results show that the Alvarez-Macovski method produced artefacts when iodine was present in the phantom together with human tissues since this method can only work with one doublet. It was shown that these artefacts could be avoided with DIRA using the (water, bone) doublet for tissues and the (water, iodine) doublet for the iodine solution., (© The Author(s) 2021. Published by Oxford University Press.)

Diamond detectors are increasingly employed in dosimetry. Their response has been investigated by means of Monte Carlo (MC) methods, but there is no consensus on what mass density ρ , mean excitation energy I and number of conduction electrons per atom n ce to use in the simulations. The ambiguity occurs due to its seeming similarity with graphite (both are carbon allotropes). Except for the difference in ρ between crystalline graphite (2.265 g cm -3 ) and diamond (3.515 g cm -3 ), their dielectric properties are assumed to be identical. This is incorrect, and the two materials should be distinguished: ( ρ = 2.265 g cm -3 , I = 81.0 eV, n ce = 1) for graphite and ( ρ = 3.515 g cm -3 , I = 88.5 eV, n ce = 0) for diamond. Simulations done with the MC codepenelopeshow that the energy imparted in diamond decreases by up to 1% with respect to 'pseudo-diamond' ( ρ = 3.515 g cm -3 , I = 81.0 eV, n ce = 0) depending on the beam quality and cavity thickness. The energy imparted changed the most in cavities that are small compared with the range of electrons. The difference in the density-effect term relative to graphite was the smallest for diamond owing to an interplay effect that ρ , I and n ce have on this term, in contrast to pseudo-diamond media when either ρ or I alone were adjusted. The study also presents a parameterized density-effect correction function for diamond that may be used by MC codes like EGSnrc. Theestarprogram assumes that n ce = 2 for all carbon-based materials, hence it delivers an erroneous density-effect correction term for graphite and diamond. Despite the small changes of the energy imparted in diamond simulated with two different I values and expected close-to-negligible deviation from the published small-field output correction data, it is important to pay attention to material properties and model the medium faithfully., (Creative Commons Attribution license.)

Purpose: Experimental dosimetry of high-dose-rate (HDR) 192 Ir brachytherapy (BT) sources is complicated due to high dose and dose-rate gradients, and softening of photon energy spectrum with depth. A single crystal synthetic diamond detector microDiamond (PTW 60019, Freiburg, Germany) has a small active volume, high sensitivity, direct readout, and nearly water-equivalent active volume. The purpose of this study was to evaluate the suitability of microDiamond detectors for the determination of absorbed dose to water around HDR 192 Ir BT sources. Three microDiamond detectors were used, allowing for the comparison of their properties., Methods: In-phantom measurements were performed using microSelectron and VariSource iX HDR 192 Ir BT treatment units. Their treatment planning systems (TPSs), Oncentra (v. 4.3) and BrachyVision (v. 13.6), respectively, were used to create irradiation plans for a cubic PMMA phantom with the microDiamond positioned at one of three source-to-detector distances (SDDs) (1.5, 2.5, and 5.5 cm) at a time. The source was stepped in increments of 0.5 cm over a total length of 6 cm to yield absorbed dose of 2 Gy at the nominal reference-point of the detector. Detectors were calibrated in 60 Co beam in terms of absorbed dose to water, and Monte Carlo (MC) calculated beam quality correction factors were applied to account for absorbed-dose energy dependence. Phantom correction factors were applied to account for differences in dimensions between the measurement phantom and a water phantom used for absorbed dose calculations made with a TPS. The same measurements were made with all three of the detectors. Additionally, dose-rate dependence and stability of the detectors were evaluated in 60 Co beam., Results: The percentage differences between experimentally determined and TPS-calculated absorbed doses to water were from -1.3% to +2.9%. The values agreed to within experimental uncertainties, which were from 1.9% to 4.3% (k = 2) depending on the detector, SDD and treatment delivery unit. No dose-rate or intrinsic energy dependence corrections were applied. All microDiamonds were comparable in terms of preirradiation dose, stability of the readings and energy response, and showed a good agreement., Conclusions: The results indicate that the microDiamond is potentially suitable for the determination of absorbed dose to water around HDR 192 Ir BT sources and may be used for independent verification of TPS's calculations, as well as for QA measurements of HDR 192 Ir BT treatment delivery units at clinical sites., (© 2017 American Association of Physicists in Medicine.)

Monte Carlo (MC) calculated detector-specific output correction factors for small photon beam dosimetry are commonly used in clinical practice. The technique, with a geometry description based on manufacturer blueprints, offers certain advantages over experimentally determined values but is not free of weaknesses. Independent MC calculations of output correction factors for a PTW-60019 micro-diamond detector were made using the EGSnrc and PENELOPE systems. Compared with published experimental data the MC results showed substantial disagreement for the smallest field size simulated ([Formula: see text] mm). To explain the difference between the two datasets, a detector was imaged with x rays searching for possible anomalies in the detector construction or details not included in the blueprints. A discrepancy between the dimension stated in the blueprints for the active detector area and that estimated from the electrical contact seen in the x-ray image was observed. Calculations were repeated using the estimate of a smaller volume, leading to results in excellent agreement with the experimental data. MC users should become aware of the potential differences between the design blueprints of a detector and its manufacturer production, as they may differ substantially. The constraint is applicable to the simulation of any detector type. Comparison with experimental data should be used to reveal geometrical inconsistencies and details not included in technical drawings, in addition to the well-known QA procedure of detector x-ray imaging.

Purpose: In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) (192)Ir source and a virtual water phantom were designed, which can be imported into a TPS., Methods: A hypothetical, generic HDR (192)Ir source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic (192)Ir source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra(®) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision ACUROS™ ]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and PENELOPE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)(3) voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR (192)Ir source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods., Results: TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ace algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 ± 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ace at clinically relevant distances., Conclusions: A hypothetical, generic HDR (192)Ir source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.

Interest in high dose rate (HDR) electronic brachytherapy operating at 50 kV is increasing. For quality assurance it is important to identify dosimetry systems that can measure the absorbed doses in absolute terms which is difficult in this energy region. In this work a comparison is made between two dosimetry systems, EPR lithium formate dosimeters and radiochromic EBT2 film. Both types of dosimeters were irradiated simultaneously in a PMMA phantom using the Axxent EBS. Absorbed dose to water was determined at distances of 10 mm, 30 mm and 50 mm from the EBS. Results were traceable to different primary standards as regards to absorbed dose to water (EPR) and air kerma (EBT2). Monte Carlo simulations were used in absolute terms as a third estimate of absorbed dose to water. Agreement within the estimated expanded (k = 2) uncertainties (5% (EPR), 7% (EBT2)) was found between the results at 30 mm and 50 mm from the x-ray source. The same result was obtained in 4 repetitions of irradiation, indicating high precision in the measurements with both systems. At all distances, agreement between EPR and Monte Carlo simulations was shown as was also the case for the film measurements at 30mm and 50mm. At 10mm the geometry for the film measurements caused too large uncertainty in measured values depending on the exact position (within sub-mm distances) of the EBS and the 10 mm film results were exculded from comparison. This work has demonstrated good performance of the lithium formate EPR dosimetry system in accordance with earlier experiments at higher photon energies ((192)Ir HDR brachytherapy). It was also highlighted that there might be issues regarding the energy dependence and intrinsic efficiency of the EBT2 film that need to be considered for measurements using low energy sources.