A Monte Carlo-based dosimetric characterization of Esteya® , an electronic surface brachytherapy unit

Dosimetric impact of physics libraries for electronic brachytherapy Monte Carlo studies

BACKGROUND

Low-energy x-ray beams used in electronic brachytherapy (eBT) present significant dosimetric challenges due to their high depth-dose gradients, the dependence of detector response on materials, and the lack of standardized dose-to-water references. These challenges have driven the need for Monte Carlo (MC) simulations to ensure accurate dosimetry. However, discrepancies in the physics models used by different MC systems have raised concerns about their dosimetric consistency, particularly in modeling bremsstrahlung interactions.

PURPOSE

To assess the dosimetric impact of using different physics approaches in three state-of-the-art MC systems for eBT, focusing on the disagreements observed when different MC methods are used to evaluate bremsstrahlung interactions.

METHODS

The MC studies of the Axxent S700, the Esteya, and the INTRABEAM eBT systems were performed using two EGSnrc applications (egs_brachy and egs_kerma), TOPAS, and PENELOPE-2018 (PEN18). The fluence spectra and depth doses were compared for simplified x-ray tube models, which maintain the target mode (transmission or reflection), the target material and thickness, and the surface applicators' source-to-surface distance. An extra simulation was made to evaluate the utility of the simplified models as proxies in predicting the most important characteristics of an accurate applicator's simulation (detailed model of INTRABEAM's 30 mm surface applicator). The EGSnrc applications and PEN18 utilized their default bremsstrahlung angular emission approaches. TOPAS used two physics lists: g4em-livermore (TOPASliv) and g4em-penelope (TOPASpen).

RESULTS

The most significant differences between MC codes were observed for the transmission target mode. The bremsstrahlung component of the fluence spectra differed by about 15% on average, comparing PEN18, EGSnrc applications, and TOPASliv, with PEN18's fluences consistently lower. EGSnrc and PEN18 agreed within 3% for their characteristic spectrum components. However, PEN18's characteristic lines overreached TOPASliv's by 40%. Those spectral characteristics generated depth dose differences, where PEN18, on average, scored 9% lower than EGSnrc and TOPASliv. Considering TOPASpen in the transmission mode, PEN18's fluence spectrum presented a lower bremsstrahlung (5%) but a higher characteristic component (10%); these spectral differences compensated, generating depth dose differences within 1% average. In the reflection target mode, EGSnrc and PEN18 agreed within 4% for the bremsstrahlung and characteristic components of the fluence spectra. With TOPASpen in the reflection mode, PEN18 presents 12% lower fluences in the bremsstrahlung component but 6% higher characteristic lines. This spectral behavior diminished the depth dose differences up to 3%.

CONCLUSION

This work found considerable disagreements between three state-of-the-art MC systems commonly used in medical applications when simulating bremsstrahlung in eBT. The differences arose when the bremsstrahlung angular distribution and the atomic relaxation processes in the target became relevant. More theoretical and experimental studies are necessary to evaluate the impact of these differences on related calculations.

A formalism for traceable dosimetry in superficial electronic brachytherapy (eBT)

Objective. Despite the number of treatments performed with electronic brachytherapy (eBT) there is no uniform methodology for reference dosimetry for international traceability to primary dosimetry standards in different eBT systems. The objective of this study is to propose a formalism for traceability reference dosimetry in superficial eBT, that is easy to apply in the clinic. This method was investigated for an Elekta Esteya with one applicator.Approach. The calibration x-ray spectrum at the primary standards dosimetry laboratory was matched to the measured eBT photon spectrum. Subsequently, two ionization chambers of different types were calibrated at the primary standard dosimetry laboratory (PSDL) in terms of air kerma against a primary standard. The chambers were used to measure ionization chamber reading ratios in-air at different distances from the applicator. Monte Carlo based air kerma ratios were calculated at different positions from the eBT applicator as well as backscatter factors in water and average mass energy absorption ratios in water and in air. Relative measurements with radiochromic films were performed in a water phantom to determine the ratio of absorbed dose to water,Dw, at the surface toDwat 1 cm depth in water. These were compared with Monte Carlo calculations.Main results. Calculations and measurements were combined to estimate theDwat the surface and at 1 cm depth in water. Ionization chamber agreement of the surface dose was 1.7%, within an uncertainty of 6.8% (k= 2). They agreed with the manufacturer dosimetry within 1.8%, with an uncertainty of 5.0% (k= 2). The feasibility of the formalism and methodology for the Esteya system was demonstrated.Significance. This study proposes a method for harmonization of traceable reference dosimetry for eBT contact treatments which does not involve a detailed simulation of the ionization chamber. The method demonstrated feasibility for one eBT system using one surface applicator. In the future the method could be applied for different eBT systems.

A study of Type B uncertainties associated with the photoelectric effect in low-energy Monte Carlo simulations

Purpose.To estimate Type B uncertainties in absorbed-dose calculations arising from the different implementations in current state-of-the-art Monte Carlo (MC) codes of low-energy photon cross-sections (<200 keV).Methods.MC simulations are carried out using three codes widely used in the low-energy domain: PENELOPE-2018, EGSnrc, and MCNP. Three dosimetry-relevant quantities are considered: mass energy-absorption coefficients for water, air, graphite, and their respective ratios; absorbed dose; and photon-fluence spectra. The absorbed dose and the photon-fluence spectra are scored in a spherical water phantom of 15 cm radius. Benchmark simulations using similar cross-sections have been performed. The differences observed between these quantities when different cross-sections are considered are taken to be a good estimator for the corresponding Type B uncertainties.Results.A conservative Type B uncertainty for the absorbed dose (k = 2) of 1.2%-1.7% (<50 keV), 0.6%-1.2% (50-100 keV), and 0.3% (100-200 keV) is estimated. The photon-fluence spectrum does not present clinically relevant differences that merit considering additional Type B uncertainties except for energies below 25 keV, where a Type B uncertainty of 0.5% is obtained. Below 30 keV, mass energy-absorption coefficients show Type B uncertainties (k = 2) of about 1.5% (water and air), and 2% (graphite), diminishing in all materials for larger energies and reaching values about 1% (40-50 keV) and 0.5% (50-75 keV). With respect to their ratios, the only significant Type B uncertainties are observed in the case of the water-to-graphite ratio for energies below 30 keV, being about 0.7% (k = 2).Conclusions.In contrast with the intermediate (about 500 keV) or high (about 1 MeV) energy domains, Type B uncertainties due to the different cross-sections implementation cannot be considered subdominant with respect to Type A uncertainties or even to other sources of Type B uncertainties (tally volume averaging, manufacturing tolerances, etc). Therefore, the values reported here should be accommodated within the uncertainty budget in low-energy photon dosimetry studies.

Depth-dose measurement corrections for the surface electronic brachytherapy beams of an Esteya® unit: a Monte Carlo study

Three different correction factors for measurements with the parallel-plate ionization chamber PTW T34013 on the Esteya electronic brachytherapy unit have been investigated. This chamber type is recommended by AAPM TG-253 for depth-dose measurements in the 69.5 kV x-ray beam generated by the Esteya unit.

Monte Carlo simulations using the PENELOPE-2018 system were performed to determine the absorbed dose deposited in water and in the chamber sensitive volume at different depths with a Type A uncertainty smaller than 0.1%. Chamber-to-chamber differences have been explored performing measurements using three different chambers. The range of conical applicators available, from 10 to 30 mm in diameter, has been explored.

Using a depth-independent global chamber perturbation correction factor without a shift of the effective point of measurement yielded differences between the absorbed dose to water and the corrected absorbed dose in the sensitive volume of the chamber of up to 1% and 0.6% for the 10 mm and 30 mm applicators, respectively. Calculations using a depth-dependent perturbation factor, including or excluding a shift of the effective point of measurement, resulted in depth-dose differences of about ± 0.5% or less for both applicators. The smallest depth-dose differences were obtained when a shift of the effective point of measurement was implemented, being displaced 0.4 mm towards the center of the sensitive volume of the chamber. The correction factors were obtained with combined uncertainties of 0.4% (k = 2). Uncertainties due to chamber-to-chamber differences are found to be lower than 2%.

The results emphasize the relevance of carrying out detailed Monte Carlo studies for each electronic brachytherapy device and ionization chamber used for its dosimetry.

On the use of the absorbed depth-dose measurements in the beam calibration of a surface electronic high-dose-rate brachytherapy unit, a Monte Carlo-based study

PURPOSE

To evaluate the use of the absorbed depth-dose as a surrogate of the half-value layer in the calibration of a high-dose-rate electronic brachytherapy (eBT) equipment. The effect of the manufacturing tolerances and the absorbed depth-dose measurement uncertainties in the calibration process are also addressed.

METHODS

The eBT system Esteya^® (Elekta Brachytherapy, Veenendaal, The Netherlands) has been chosen as a proof-of-concept to illustrate the feasibility of the proposed method, using its 10 mm diameter applicator. Two calibration protocols recommended by the AAPM (TG-61) and the IAEA (TRS-398) for low-energy photon beams were evaluated. The required Monte Carlo (MC) simulations were carried out using PENELOPE2014. Several MC simulations were performed modifying the flattening filter thickness and the x-ray tube potential, generating one absorbed depth-dose curve and a complete set of parameters required in the beam calibration (i.e., HVL, backscatter factor (B_w ), and mass energy-absorption coefficient ratios (µ_en /ρ)_water,air ), for each configuration. Fits between each parameter and some absorbed dose-ratios calculated from the absorbed depth-dose curves were established. The effect of the manufacturing tolerances and the absorbed dose-ratio uncertainties over the calibration process were evaluated by propagating their values over the fitting function, comparing the overall calibration uncertainties against reference values. We proposed four scenarios of uncertainty (from 0% to 10%) in the dose-ratio determination to evaluate its effect in the calibration process.

RESULTS

The manufacturing tolerance of the flattening filter (±0.035 mm) produces a change of 1.4% in the calculated HVL and a negligible effect over the B_w , (µ_en /ρ)_water,air , and the overall calibration uncertainty. A potential variation of 14% of the electron energies due to manufacturing tolerances in the x-ray tube (69.5 ± ~10 keV) generates a variation of 10% in the HVL. However, this change has a negligible effect over the B_w and (µ_en /ρ)_water,air , adding 0.1% to the overall calibration uncertainty. The fitting functions reproduce the data with an uncertainty (k = 2) below 1%, 0.5%, and 0.4% for the HVL, B_w , and (µ_en /ρ)_water,air , respectively. The four studied absorbed dose-ratio uncertainty scenarios add, in the worst-case scenario, 0.2% to the overall uncertainty of the calibration process.

CONCLUSIONS

This work shows the feasibility of using the absorbed depth-dose curve in the calibration of an eBT system with minimal loss of precision.