Decomposition of the dose conversion factor based on fluence spectra of secondary charged particles: Application to lateral dose profiles in photon fields

Investigations on the beam quality correction factor for ionization chambers in high-energy brachytherapy dosimetry

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.

Enhancement of the EGSnrc code egs_chamber for fast fluence calculations of charged particles

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.

Cema-based formalism for the determination of absorbed dose for high-energy photon beams

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 k_Q,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, k_Q,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.

Suitability of the microDiamond detector for experimental determination of the anisotropy function of High Dose Rate 192 Ir brachytherapy sources

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) ¹⁹² Ir brachytherapy sources.

METHODS

The HDR ¹⁹² 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 ¹⁹² 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 ¹⁹² Ir brachytherapy sources and support its employment for the determination and validation of TG-43 parameters of such sources.

Density compensated diodes for small field dosimetry: comprehensive testing and implications for design

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 cm². 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 cm² 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 cm² 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.

[Detector Based Determination of Water Absorbed Dose According to DIN 6800 Teil 1: Suggestion for an Extension of the Fundamental Formalism]

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.