Percentage depth dose evaluation in heterogeneous media using thermoluminescent dosimetry

GATE Monte Carlo approach to heterogeneity dose distribution in small fields used in radiation therapy

The Accurate dosage prediction in Radiation Therapy is challenging, prompting a need for precision beyond conventional clinical Treatment Planning Systems (TPS). Monte Carlo-based methods are sought for their superior accuracy. The aim of this study is to compare dose distributions between the ACUROS algorithm and the GATE platform in various tissue densities and field sizes, focusing on smaller fields. This study was initiated with a homogeneous validation of the TrueBeam STX system, using measurements obtained from the Centre Hospitalier Interregional Edith Cavell (CHIREC) in Brussels. The validation compared dosimetric functions (Percentage Depth Dose (PDD), Dose profile (DP) and Collimator scatter fraction (CSF)) employing the GAMMA index with a 2% / 2 mm criterion tolerance. Following this, heterogeneous studies examined dose distributions between the ACUROS algorithm and the GATE platform in various tissue densities and field sizes, with a specific focus on smaller fields. Simulations were conducted using both platforms on chest phantoms with heterogeneous slabs representing bone, lung, and heart, each housing a central tumor. The impact of electronic equilibrium on tumors for different small field sizes was evaluated. Results showed a remarkable 99% agreement between measurements and GATE calculations in the homogeneous validation of the TrueBeam STX system. However, in heterogeneous studies, ACUROS consistently overestimated lung doses by up to 8% compared to GATE simulation, especially evident with a flattening filter and smaller beam sizes at density interfaces. This highlights significant dose estimation discrepancies between ACUROS and GATE, emphasizing the need for precise calculations. The findings support exploring Monte Carlo-based methods for enhanced accuracy in Radiation Therapy treatment planning.

Dosimetric Analysis of Rhizophora-based Phantom Material in Radiation Therapy Applications Using Monte Carlo GATE Simulation

Purpose

This study aims to determine the percentage depth dose (PDD) of a phantom material made from soy-lignin bonded Rhizophora spp. particleboard coated with a gloss finish by using Monte Carlo Geant4 Application for Tomographic Emission (GATE) simulation.

Materials and Methods

The particleboard was fabricated using a hot pressing technique at target density of 1.0 g·cm-3 and the elemental fraction was recorded for the simulation. The PDD was simulated in the GATE simulation using the linear accelerator Elekta Synergy model for the water phantom and Rhizophora phantom, and the results were compared with the experimental PDD performed by several studies. Beam flatness and beam symmetry were also measured in this study.

Results

The simulated PDD for Rhizophora and water was in agreement with the experimental PDD of water with overall discrepancies of 0% to 8.7% at depth ranging from 1.0 to 15.0 cm. In the GATE simulation, all the points passed the clinical 3%/3 mm criterion in comparison with water, with the final percentage of 2.34% for Rhizophora phantom and 2.49% for the water phantom simulated in GATE. Both the symmetries are all within the range of an acceptable value of 2.0% according to the recommendation, with the beam symmetry of the water phantom and Rhizophora phantom at 0.58% and 0.28%, respectively.

Conclusions

The findings of this study provide the necessary foundation to confidently use the phantom for radiotherapy purposes, especially in treatment planning.

Evaluation of lung heterogeneity effects on dosimetric parameters in small photon fields using MAGIC polymer gel, Gafchromic film, and Monte Carlo simulation

In this work, the performance of MAGIC polymer gel in measuring dosimetric parameters beyond lung heterogeneity in small fields was investigated. All data were obtained using MAGIC, EBT2, and MC in four small field sizes. The maximum local differences between MAGIC and MC were less than 5.1, 3.9, 3.1, and 2.6% for PDD values behind lung heterogeneity at 5, 10, 20, and 30 mm field sizes, respectively. The findings showed that MAGIC is a suitable tool for dosimetry behind low-density heterogeneity.

A practical method for quantifying dose in bone and lung using TLDs when using 6 and 15 MV photon beams

This paper presents a practical method for converting dose measured with thermoluminescent dosimeters (TLD) to dose in lung and bone for 6 MV and 15 MV photon beams. Monte Carlo (MC) simulations and Burlin cavity theory calculations were performed to calculate [Formula: see text], the dose-to-TLD to dose-to-medium conversion factor. A practical method was proposed for converting TLD-measured-dose to dose-in-medium using the TLD dose calibration in water and [Formula: see text] dose-to-medium to dose-to-water conversion factor. Theoretical calculations for [Formula: see text] were performed using photon spectrum weighted parameters and were compared with MC simulations. Verification of the proposed method was done using phantoms having either bone or lung equivalent slabs stacked in between solid water slabs. Percent depth dose (PDD) curves were measured using 0.089 cm thick LiF:Mg,Ti (TLD-100) dosemeters placed at various depths within these phantoms. They were then corrected with [Formula: see text] factors using the proposed dose conversion method, and were compared with the MC simulations. For 6 MV beam, the MC calculated [Formula: see text] factors were 0.942 and 1.002 for bone and lung, and for 15 MV it was 0.927 and 1.005 for bone and lung, respectively. The difference between the MC simulated and spectrum weighted theoretical [Formula: see text] factors were within 3% for both lung and bone. The PDD curves measured with TLD-100 chips that were corrected using the proposed method agreed well within 1.5% of the MC simulated PDD curves for both the water/lung/water and water/bone/water (WBW) phantoms. The dose-to-medium correction using MC simulated [Formula: see text] is convenient, easy, and accurate. Therefore, it can be used instead of Burlin cavity theory, especially in media with high atomic numbers such as bone for accurate dose quantification.

Validation of the RayStation Monte Carlo dose calculation algorithm using a realistic lung phantom

PURPOSE

Our purposes are to compare the accuracy of RaySearch's analytical pencil beam (APB) and Monte Carlo (MC) algorithms for clinical proton therapy and to present clinical validation data using a novel animal tissue lung phantom.

METHODS

We constructed a realistic lung phantom composed of a rack of lamb resting on a stack of rectangular natural cork slabs simulating lung tissue. The tumor was simulated using 70% lean ground lamb meat inserted in a spherical hole with diameter 40 ± 5 mm carved into the cork slabs. A single-field plan using an anterior beam and a two-field plan using two anterior-oblique beams were delivered to the phantom. Ion chamber array measurements were taken medial and distal to the tumor. Measured doses were compared with calculated RayStation APB and MC calculated doses.

RESULTS

Our lung phantom enabled measurements with the MatriXX PT at multiple depths in the phantom. Using the MC calculations, the 3%/3 mm gamma index pass rates, comparing measured with calculated doses, for the distal planes were 74.5% and 85.3% for the APB and 99.1% and 92% for the MC algorithms. The measured data revealed up to 46% and 30% underdosing within the distal regions of the target volume for the single and the two field plans when APB calculations are used. These discrepancies reduced to less than 18% and 7% respectively using the MC calculations.

CONCLUSIONS

RaySearch Laboratories' Monte Carlo dose calculation algorithm is superior to the pencil-beam algorithm for lung targets. Clinicians relying on the analytical pencil-beam algorithm should be aware of its pitfalls for this site and verify dose prior to delivery. We conclude that the RayStation MC algorithm is reliable and more accurate than the APB algorithm for lung targets and therefore should be used to plan proton therapy for patients with lung cancer.

Estimation of Dosimetric Parameters based on KNR and KNCSF Correction Factors for Small Field Radiation Therapy at 6 and 18 MV Linac Energies using Monte Carlo Simulation Methods

BACKGROUND

Estimating dosimetric parameters for small fields under non-reference conditions leads to significant errors if done based on conventional protocols used for large fields in reference conditions. Hence, further correction factors have been introduced to take into account the influence of spectral quality changes when various detectors are used in non-reference conditions at different depths and field sizes.

OBJECTIVE

Determining correction factors (K_NR and K_NCSF) recommended recently for small field dosimetry formalism by American Association of Physicists in Medicine (AAPM) for different detectors at 6 and 18 MV photon beams.

METHODS

EGSnrc Monte Carlo code was used to calculate the doses measured with different detectors located in a slab phantom and the recommended K_NR and K_NCSF correction factors for various circular small field sizes ranging from 5-30 mm diameters. K_NR and K_NCSF correction factors were determined for different active detectors (a pinpoint chamber, EDP-20 and EDP-10 diodes) in a homogeneous phantom irradiated to 6 and 18 MV photon beams of a Varian linac (2100C/D).

RESULTS

K_NR correction factor estimated for the highest small circular field size of 30 mm diameter for the pinpoint chamber, EDP-20 and EDP-10 diodes were 0.993, 1.020 and 1.054; and 0.992, 1.054 and 1.005 for the 6 and 18 MV beams, respectively. The K_NCSF correction factor estimated for the lowest circular field size of 5 mm for the pinpoint chamber, EDP-20 and EDP-10 diodes were 0.994, 1.023, and 1.040; and 1.000, 1.014, and 1.022 for the 6 and 18 MV photon beams, respectively.

CONCLUSION

Comparing the results obtained for the detectors used in this study reveals that the unshielded diodes (EDP-20 and EDP-10) can confidently be recommended for small field dosimetry as their correction factors (K_NR and K_NCSF) was close to 1.0 for all small field sizes investigated and are mainly independent from the electron beam spot size.

Development of a Geant4-based independent patient dose validation system with an elaborate multileaf collimator simulation model

Despite the improvements in the dose calculation models of the commercial treatment planning systems (TPS), their ability to accurately predict patient dose is still limited. One of the limitations is caused by the simplified model of the multileaf collimator (MLC). The aim of this study was to develop a Monte Carlo (MC) method‐based independent patient dose validation system with an elaborate MLC model for more accurate dose evaluation. Varian Clinac 2300 IX was simulated using Geant4 toolkits, after which MC commissioning with measurements was performed to validate the simulation model. A DICOM‐RT interface was developed to obtain the beam delivery conditions including the hundreds of MLC motions. Finally, the TPS dose distributions were compared with the MC dose distributions for water phantom cases and a patient case. Our results show that the TPS overestimated the absolute abutting leakage dose in the closed MLC field, with about 20% more of the maximum dose than that of the MC calculation. For water phantom cases, the dose distributions inside the target region were almost identical with the dose difference of less than 2%, while the dose near the edge of the target shows difference about 10% between Geant4 and TPS due to geometrical differences in MLC model. For the patient analysis, the Geant4 and TPS doses of all organs were matched well within 1.4% of the prescribed dose. However, for organs located in areas with high ratio of leaf pairs with distances less than 10 mm leaf pair (LP_(<10mm)), the maximum dose of TPS was overestimated by about 3% of the prescribed dose. These dose comparison results demonstrate that our system for calculating the patient dose is quite accurate. Furthermore, if the MLC sequences in treatment plan have a large ratio of LP_(short), more than 3% dose difference in normal tissue could be seen.

A depth dose study between AAA and AXB algorithm against Monte Carlo simulation using AIP CT of a 4D dataset from a moving phantom

Aim

To identifying depth dose differences between the two versions of the algorithms using AIP CT of a 4D dataset.

Background

Motion due to respiration may challenge dose prediction of dose calculation algorithms during treatment planning.

Materials and methods

The two versions of depth dose calculation algorithms, namely, Anisotropic Analytical Algorithm (AAA) version 10.0 (AAAv10.0), AAA version 13.6 (AAAv13.6) and Acuros XB dose calculation (AXB) algorithm version 10.0 (AXBv10.0), AXB version 13.6 (AXBv13.6), were compared against a full MC simulated 6X photon beam using QUASAR respiratory motion phantom with a moving chest wall. To simulate the moving chest wall, a 4 cm thick wax mould was attached to the lung insert of the phantom. Depth doses along the central axis were compared in the anterior and lateral beam direction for field sizes 2 × 2 cm², 4 × 4 cm² and 10 × 10 cm².

Results

For the lateral beam direction, the moving chest wall highlighted differences of up to 105% for AAAv10.0 and 40% for AXBv10.0 from MC calculations in the surface and buildup doses. AAAv13.6 and AXBv13.6 agrees with MC predictions to within 10% at similar depth. For anterior beam doses, dose differences predicted for both versions of AAA and AXB algorithm were within 7% and results were consistent with static heterogeneous studies.

Conclusions

The presence of the moving chest wall was capable of identifying depth dose differences between the two versions of the algorithms. These differences could not be identified in the static chest wall as shown in the anterior beam depth dose calculations.

Evaluation of AAA and XVMC Algorithms for Dose Calculation in Lung Equivalent Heterogeneity in Photon Fields: A Comparison of Calculated Results with Measurements

AIMS

The aims of the present work are (1) to evaluate dose calculation accuracy of two commonly used algorithms for 15 MV small photon fields in a medium encompassing heterogeneity and (2) to compare them with measured results obtained from gafchromic film EBT2.

MATERIALS AND METHODS

Authors employed kailwood (Pinus Wallichiana) to mimic lung. Briefly, seven Kailwood plates, each measuring 25x25 cm² of varying thicknesses totaling 13 cm equivalent to the mean thickness of an adult human lung, were sandwiched between 5 cm tissue equivalent material from top and 10 cm below. Physical measurements were performed using Radiochromic film EBT2. The field sizes of 1x1, 2x2, 5x5 and 10x10 cm² were selected at 100 cm SSD. Simulations were performed using EGSnrc/DOSRZnrc Monte Carlo code. The dose variation inside the inhomogeneity and near the interface was calculated using AAA & XVMC algorithm.

RESULTS

Preliminary results show that there is large variation of dose inside inhomogeneity. The maximum variation of dose inside the inhomogeneity for 1x1 cm² was found 40% by AAA and 4.5% by XVMC compared to measured/simulated results. For the field size of 2x2 cm², these figures were 27% by AAA & 3.5% by XVMC. For 5x5 cm² field size, the variation is small which becomes insignificant for larger fields.

CONCLUSION

The results presented in this work indicate that for smaller fields, XVMC algorithm gives more realistic prediction, while there is the need for caution on using AAA algorithm for dose calculations involving small area irradiation encompassing heterogeneities and low-density media.

Independent dose validation system for Gamma Knife radiosurgery, using a DICOM-RT interface and Geant4

Leksell GammaPlan was specifically designed for Gamma Knife (GK) radiosurgery planning, but it has limited accuracy for estimating the dose distribution in inhomogeneous areas, such as the embolization of arteriovenous malformations. We aimed to develop an independent patient dose validation system based on a patient-specific model, constructed using a DICOM-RT interface and the Geant4 toolkit. Leksell Gamma Knife Perfexion was designed in Geant4.10.00 and includes a DICOM-RT interface. Output factors for each collimator in a sector and dose distributions in a spherical water phantom calculated using a Monte Carlo (MC) algorithm were compared with the output factors calculated by the tissue maximum ratio (TMR) 10 algorithm and dose distributions measured using film, respectively. Studies using two types of water phantom and two patient simulation cases were evaluated by comparing the dose distributions calculated by the MC, the TMR and the convolution algorithms. The water phantom studies showed that if the beam size is small and the target is located in heterogeneous media, the dose difference could be up to 11%. In the two patient simulations, the TMR algorithm overestimated the dose by about 4% of the maximum dose if a complex and large bony structure was located on the beam path, whereas the convolution algorithm showed similar results to those of the MC algorithm. This study demonstrated that the in-house system could accurately verify the patient dose based on full MC simulation and so would be useful for patient cases where the dose differences are suspected.

Comparison of CCC and ETAR dose calculation algorithms in pituitary adenoma radiation treatment planning; Monte Carlo evaluation

Aims

To verify the accuracy of two common absorbed dose calculation algorithms in comparison to Monte Carlo (MC) simulation for the planning of the pituitary adenoma radiation treatment.

Materials and methods

After validation of Linac's head modelling by MC in water phantom, it was verified in Rando phantom as a heterogeneous medium for pituitary gland irradiation. Then, equivalent tissue-air ratio (ETAR) and collapsed cone convolution (CCC) algorithms were compared for a conventional three small non-coplanar field technique. This technique uses 30 degree physical wedge and 18 MV photon beams.

Results

Dose distribution findings showed significant difference between ETAR and CCC of delivered dose in pituitary irradiation. The differences between MC and dose calculation algorithms were 6.40 ± 3.44% for CCC and 10.36 ± 4.37% for ETAR. None of the algorithms could predict actual dose in air cavity areas in comparison to the MC method.

Conclusions

Difference between calculation and true dose value affects radiation treatment outcome and normal tissue complication probability. It is of prime concern to select appropriate treatment planning system according to our clinical situation. It is further emphasised that MC can be the method of choice for clinical dose calculation algorithms verification.

Impact of the calculation resolution of AAA for small fields and RapidArc treatment plans

PURPOSE

To investigate the impact of the calculation resolution of the anisotropic analytical algorithms (AAA) for a variety of small fields in homogeneous and heterogeneous media and for RapidArc plans.

METHODS

Dose distributions calculated using AAA version 8.6.15 (AAA8) and 10.0.25 (AAA10) were compared to measurements performed with GafChromic EBT film, using phantoms made of polystyrene or a combination of polystyrene and cork. The accuracy of the algorithms calculated using grid resolutions of 2.5 and 1.0 mm was investigated for different field sizes, and for a limited selection of RapidArc plans (head and neck, small meningioma, and lung). Additional plans were optimized to create excessive multileaf collimator modulation and measured on a homogenous phantom. Gamma evaluation criterion of 3% dose difference and 2- or 1-mm distance to agreement (DTA) were applied to evaluate the accuracy of the algorithms.

RESULTS

For fields < or = 3 x 3 cm2, both versions of AAA predicted lower peak doses and broader penumbra widths than the measurements. However, AAA10 and a finer calculation grid improved the agreement. For RapidArc plans with many small multileaf collimator (MLC) segments and relatively high number of monitor units (MU), AAA8 failed to identify small dose peaks within the target. Both versions performed better in polystyrene than in cork. In homogeneous cork layers, AAA8 underestimated the average target dose for a clinical lung plan. This was improved with AAA10 calculated using a 1 mm grid.

CONCLUSIONS

AAA10 improves the accuracy of dose calculations, and calculation grid of 1.0 mm is superior to using 2.5 mm, although calculation times increased by factor of 5. A suitable upper MU constraint should be assigned during optimization to avoid plans with high modulation. For plans with a relative high number of monitor units, calculations using 1 mm grid resolution are recommended. For planning target volume (PTV) which contains relatively large area of low density tissue, users should be aware of possible dose underestimation in the low density region and recalculation with AAA10 grid 1.0 mm is recommended.

Effect of low-density heterogeneities in telecobalt therapy and validation of dose calculation algorithm of a treatment planning system

Telecobalt machines are still prominently used for the treatment of a variety of cancer cases in developing countries. The human body is a heterogeneous composition of variety of tissues and cavities which vary widely in their physical and radiological properties. The presence of heterogeneities in the path of telecobalt beam presents an altered dose distribution in the region of clinical interests. A computerized treatment planning system (TPS) is generally used for calculating the dose distribution in the patient. Experimental measurements were carried out in a telecobalt beam with the objectives to study the effects of low-density heterogeneities and to verify the ability of the ASHA radiotherapy TPS in predicting the altered dose distribution along the central axis and off-axis of the beam. Locally available kailwood was tested for its lung equivalence and measurements were carried out in a polymethyl methacrylate phantom by introducing lung equivalent and air gap heterogeneities. A comparison of experimentally measured and TPS calculated dose values indicates that the TPS overestimates the dose by 11.6% in lung equivalent (kailwood) heterogeneity along the central axis. Similarly, it was found that the TPS overestimates the dose by 3.9% and 5.9%, respectively, with air heterogeneity of 1.0 and 2.0 cm. While testing the adequacy of TPS in off-axis region, it was found that the TPS calculation does not indicate the widening of the beam profile in the low-density heterogeneity region. This study suggests that the effective path length based algorithm of the ASHA radiotherapy TPS is unable to achieve the recommended 3% accuracy of clinical dose calculation in heterogeneous media.

Dosimetric evaluation of Acuros XB Advanced Dose Calculation algorithm in heterogeneous media

BACKGROUND

A study was realised to evaluate and determine relative figures of merit of a new algorithm for photon dose calculation when applied to inhomogeneous media.

METHODS

The new Acuros XB algorithm implemented in the Varian Eclipse treatment planning system was compared against a Monte Carlo method (VMC++), and the Analytical Anisotropic Algorithm (AAA). The study was carried out in virtual phantoms characterized by simple geometrical structures. An insert of different material and density was included in a phantom built of skeletal-muscle and HU = 0 (setting "A"): Normal Lung (lung, 0.198 g/cm3); Light Lung (lung, 0.035 g/cm3); Bone (bone, 1.798 g/cm3); another phantom (setting "B") was built of adipose material and including thin layers of bone (1.85 g/cm3), adipose (0.92 g/cm3), cartilage (1.4745 g/cm3), air (0.0012 g/cm3). Investigations were performed for 6 and 15 MV photon beams, and for a large (13 × 13 cm2) and a small (2.8 × 13 cm2) field.

RESULTS

Results are provided in terms of depth dose curves, transverse profiles and Gamma analysis (3 mm/3% and 2 mm/2% distance to agreement/dose difference criteria) in planes parallel to the beam central axis; Monte Carlo simulations were assumed as reference. Acuros XB gave an average gamma agreement, with a 3 mm/3% criteria, of 100%, 86% and 100% for Normal Lung, Light Lung and Bone settings, respectively, and dose to medium calculations. The same figures were 86%, 11% and 100% for AAA, where only dose rescaled to water calculations are possible.

CONCLUSIONS

In conclusion, Acuros XB algorithm provides a valid and accurate alternative to Monte Carlo calculations for heterogeneity management.

Four-dimensional dosimetry validation and study in lung radiotherapy using deformable image registration and Monte Carlo techniques

Thoracic cancer treatment presents dosimetric difficulties due to respiratory motion and lung inhomogeneity. Monte Carlo and deformable image registration techniques have been proposed to be used in four-dimensional (4D) dose calculations to overcome the difficulties. This study validates the 4D Monte Carlo dosimetry with measurement, compares 4D dosimetry of different tumor sizes and tumor motion ranges, and demonstrates differences of dose-volume histograms (DVH) with the number of respiratory phases that are included in 4D dosimetry. BEAMnrc was used in dose calculations while an optical flow algorithm was used in deformable image registration and dose mapping. Calculated and measured doses of a moving phantom agreed within 3% at the center of the moving gross tumor volumes (GTV). 4D CT image sets of lung cancer cases were used in the analysis of 4D dosimetry. For a small tumor (12.5 cm³) with motion range of 1.5 cm, reduced tumor volume coverage was observed in the 4D dose with a beam margin of 1 cm. For large tumors and tumors with small motion range (around 1 cm), the 4D dosimetry did not differ appreciably from the static plans. The dose-volume histogram (DVH) analysis shows that the inclusion of only extreme respiratory phases in 4D dosimetry is a reasonable approximation of all-phase inclusion for lung cancer cases similar to the ones studied, which reduces the calculation in 4D dosimetry.