A comparison of physical and dosimetric properties of lung substitute materials

Validation of an in vivo transit dosimetry algorithm using Monte Carlo simulations and ionization chamber measurements

PURPOSE

Transit dosimetry is a safety tool based on the transit images acquired during treatment. Forward-projection transit dosimetry software, as PerFRACTION, compares the transit images acquired with an expected image calculated from the DICOM plan, the CT, and the structure set. This work aims to validate PerFRACTION expected transit dose using PRIMO Monte Carlo simulations and ionization chamber measurements, and propose a methodology based on MPPG5a report.

METHODS

The validation process was divided into three groups of tests according to MPPG5a: basic dose validation, IMRT dose validation, and heterogeneity correction validation. For the basic dose validation, the fields used were the nine fields needed to calibrate PerFRACTION and three jaws-defined. For the IMRT dose validation, seven sweeping gaps fields, the MLC transmission and 29 IMRT fields from 10 breast treatment plans were measured. For the heterogeneity validation, the transit dose of these fields was studied using three phantoms: 10 , 30 , and a 3 cm cork slab placed between 10 cm of solid water. The PerFRACTION expected doses were compared with PRIMO Monte Carlo simulation results and ionization chamber measurements.

RESULTS

Using the 10 cm solid water phantom, for the basic validation fields, the root mean square (RMS) of the difference between PerFRACTION and PRIMO simulations was 0.6%. In the IMRT fields, the RMS of the difference was 1.2%. When comparing respect ionization chamber measurements, the RMS of the difference was 1.0% both for the basic and the IMRT validation. The average passing rate with a γ(2%/2 mm, TH = 20%) criterion between PRIMO dose distribution and PerFRACTION expected dose was 96.0% ± 5.8%.

CONCLUSION

We validated PerFRACTION calculated transit dose with PRIMO Monte Carlo and ionization chamber measurements adapting the methodology of the MMPG5a report. The methodology presented can be applied to validate other forward-projection transit dosimetry software.

Low-cost chest paediatric phantom for dose optimisation: construction and validation

INTRODUCTION AND OBJECTIVES

In order to perform chest dose optimisation studies, the imaging phantom should be adequate for image quality evaluation. Since high-end phantoms are cost prohibitive, there is a need for a low-cost construction method with fairly available tissue substitutes.

MATERIALS AND METHODS

Theoretical calculations of radiological characteristics were performed for each of lung, cortical bone and soft tissues in order to choose appropriate substitute, then, cork, P.V.C. (Polyvinyl chloride) and water were chosen, respectively. Validation included, firstly, measuring CT Hounsfield Units (HU) of a real patient's tissues then compared against their corresponding anatomies in the constructed phantom. Secondly, Signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) values were acquired in this study to evaluate the quality of images generated from the constructed phantom, then, compare their trends with a valid phantom under different exposure parameters (kVp and mAs).

RESULTS

From theoretical calculations, the percentage differences showed high accuracy of tissue substitutes when simulating real patient tissues; P.V.C. was ≥5.78%, cork was ≥4.46% and water ≥5%. The percentage difference (CT HU) between lung and cortical bone and their equivalent tissue substitutes were 10.44% and 0.53%-3.17%, respectively. Strong positive correlations were found for SNR when changing both kVp (0.79) and mAs (0.65). While the correlation strength of CNR values were found to be moderate when changing both kVp (0.58) and mAs (0.53).

CONCLUSIONS

Our low-cost phantom approved through CT HU that their materials replicate the radiological characteristics of real one-year-old child while SNR and SNR correlations confirmed its applicability in imaging and optimisation studies.

Evaluating impact of medium variation on dose calculated through planning system in a low cost in-house phantom

Purpose:In radiotherapy, accuracy in dose estimation of dose calculation methods is critical. The influence of deformity on radiation dose calculations derived by planning system is evaluated in present study. The goal of study was to create a low-cost inhomogeneous phantom for measuring absorbed dose using an Ionisation chamber and Gafchromic film, which was validated using treatment planning system (TPS) dose outcome.Methods and Materials:The central axis dose calculations were computed using Pencil Beam Convolution algorithm (PBC), Collapsed Cone Convolution (CCC) and Monte Carlo (MC) algorithm in the Monaco treatment planning system using an In-house phantom (20x20x20cm³) made up of acrylic sheet containing water and inhomogeneous material wooden powder equivalent to lung. Phantom was scanned in Computed Tomography (CT) scanner and image set was sent to the planning workstation. The depth dose evaluations were performed using ionization chamber and Gafchromic film with same beam settings and monitor units in every setup. Following that, the calculated doses obtained from TPS and measured depth doses were compared.Results:The results was reported for photon energies 6MV, 10MV, 15MV, 6FFF and 10FFF at varying field sizes of 4X4 cm², 5x5 cm², 10x10 cm², and 15x15 cm². MC maximum dose variation predicted was 2.06% in 15MV of measured chamber dose and -2.06% of measured gafchromic film dose in 6MVFFF. CCC maximum dose variation predicted was 2.68% of measured chamber dose in 6MV and 3.31% of measured gafchromic film dose in 6MV whereas PB maximum dose variation predicted was -5.94% in 15MV of measured chamber dose and -11.6% of measured gafchromic film dose in 6MVFFF.Conclusion:Low-cost in-house phantoms can be utilised to assess point and planar doses during patient-specific quality assurance in centres that don't have accessibility to phantoms due to the high cost of commercially available tools.

Dosimetric evaluation of two phases of respiratory movement using a lung equivalent material for radiotherapy treatment planning

Background/aim:

Radiation dosimetry requires special phantoms which are comparable with organs and tissues of a human body. The lung is one of the organs with a low density. Therefore, it is important to create and use lung equivalent phantoms in dosimetric controls. The aim of this study was to investigate the importance of using lung equivalent phantoms for different respiratory phases during measurements with both computed tomography (CT) and linear accelerator.

Materials and methods:

The maximum lung inhalation phantom (LIP) and lung exhalation phantom (LEP) were created for two respiratory phases. The Hounsfield Unit (HU) values based on the selected slice thickness and CT tube voltages were investigated, as well as the difference between energy and algorithms used in the treatment planning system.

Results:

It was found that the change in HU values according to slice thickness were more significant in measurements for respiratory phases. The dose difference between LEP and LIP at a point which is located 1 cm below the surface of the phantoms was found as 1·0% for 6 megavolt (MV) and 2·8% for 18 MV. The highest difference between the two algorithms was found to be 7·22% for 6 MV and 10·93% for 18 MV for LIP phantom.

Conclusion:

It can be said that the LIP and LEP phantoms prepared in accordance with respiratory phases can be a simple and inexpensive method to investigate any difference in dosimetry during respiratory phases. Also, measured and calculated dose values are in good agreement when thinner slice thickness was chosen.

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.

Fabrication and characterization of polymeric cellular foams for low-density computed tomography phantom applications

Computed tomography imaging phantom devices have proven to be beneficial in improving computed tomography diagnostic techniques. Though commercial phantoms are available with tissue mimicking properties, there is a lack of low-density tissue specificity and variety. This study proposes a method for the fabrication of various low-density tissue mimicking computed tomography imaging phantoms. By illustrating the fabrication technique, material properties can be shown to be controlled and assessed against characteristic computed tomography imaging properties, most particularly, the computed tomography number in Hounsfield Units. A batch cellular foaming technique was utilized on thermoplastic polyurethane with ranging heated water bath foaming times from 0.5 to 10 min to fabricate polymeric computed tomography phantoms of controlled foam material properties. Computed tomography number values were experimentally measured. Additionally, separate experimental measurements were made on the foam characteristic properties of fabricated thermoplastic polyurethane foams. A relative decreasing trend was exhibited between the foam characteristic properties of cell density, average cell size, and material density to computed tomography number.

The dosimetric effect of electron density overrides in 3DCRT Lung SBRT for a range of lung tumor dimensions

The combined effects of lung tumor motion and limitations of treatment planning system dose calculations in lung regions increases uncertainty in dose delivered to the tumor and surrounding normal tissues in lung stereotactic body radiotherapy (SBRT). This study investigated the effect on plan quality and accuracy when overriding treatment volume electron density values. The QUASAR phantom with modified cork cylindrical inserts, each containing a simulated spherical tumor of 15‐mm, 22‐mm, or 30‐mm diameter, was used to simulate lung tumor motion. Using Monaco 5.1 treatment planning software, two standard plans (50% central phase (50%) and average intensity projection (AIP)) were compared to eight electron density overridden plans that focused on different target volumes (internal target volume (ITV), planning target volume (PTV), and a hybrid plan (HPTV)). The target volumes were set to a variety of electron densities between lung and water equivalence. Minimal differences were seen in the 30‐mm tumor in terms of target coverage, plan conformity, and improved dosimetric accuracy. For the smaller tumors, a PTV override showed improved target coverage as well as better plan conformity compared to the baseline plans. The ITV plans showed the highest gamma pass rate agreement between treatment planning system (TPS) and measured dose (P < 0.040). However, the low electron density PTV and HPTV plans also showed improved gamma pass rates (P < 0.035, P < 0.011). Low‐density PTV overrides improved the plan quality and accuracy for tumor diameters less than 22 mm only. Although an ITV override generated the most significant increase in accuracy, the low‐density PTV plans had the additional benefit of plan quality improvement. Although this study and others agreed that density overrides improve the treatment of SBRT, the optimal density override and the conditions under which it should be applied were found to be department specific, due to variations in commissioning and calculation methods.

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.

Comprehensive Investigation on Controlling for CT Imaging Variabilities in Radiomics Studies

Radiomics has shown promise in improving models for predicting patient outcomes. However, to maximize the information gain of the radiomics features, especially in larger patient cohorts, the variability in radiomics features owing to differences between scanners and scanning protocols must be accounted for. To this aim, the imaging variability of radiomics feature values was evaluated on 100 computed tomography scanners at 35 clinics by imaging a radiomics phantom using a controlled protocol and the commonly used chest and head protocols of the local clinic. We used a linear mixed-effects model to determine the degree to which the manufacturer and individual scanners contribute to the overall variability. Using a controlled protocol reduced the overall variability by 57% and 52% compared to the local chest and head protocols respectively. The controlled protocol also reduced the relative contribution of the manufacturer to the total variability. For almost all variabilities (manufacturer, scanner, and residual with different preprocesssing), the controlled protocol scans had a significantly smaller variability than the local protocol scans did. For most radiomics features, the imaging variability was small relative to the inter-patient feature variability in non-small cell lung cancer and head and neck squamous cell carcinoma patient cohorts. From this study, we conclude that using controlled scans can reduce the variability in radiomics features, and our results demonstrate the importance of using controlled protocols in prospective radiomics studies.

Simultaneous characterization of electron density and effective atomic number for radiotherapy planning using stoichiometric calibration method and dual energy algorithms

Relative electron densities of body tissues (ρ_e) for radiotherapy treatment planning are normally obtained by CT scanning of tissue substitute materials (TSMs) and producing a Hounsfield Unit-ρ_e calibration curve. Aiming for more accurate, simultaneous characterization of ρ_e and effective atomic number (Z_eff) of real tissues, an in-house phantom (including 10 water solutions plus composite cork as TSMs) was constructed and scanned at 4 kVps. Dual-energy algorithms were applied to 80-140 and 100-140 kVp combination scans, for better differentiation of tissues with same attenuation coefficient at 120 kVp but different ρ_e and Z_eff. Stoichiometric calibration and closeness of the ρ_e of the 11 TSMs to real tissues (≤ 0.5%) resulted in smaller ρ_e calculation discrepancies, compared to studies with commercial phantoms (p < 0.024). Applying an energy subtraction algorithm further mitigated errors by spectral separation and reduction of beam hardening artifacts and noise, reducing the mean and standard deviation of the absolute difference of ρ_e at 80-140 kVp (p < 0.003) and 100-140 kVp (p < 0.0001) scans, compared to 120 kVp scan, respectively. Moreover, a parametrization algorithm decreased the Z_eff discrepancy from real tissues at 80-140 kVp scans; for thyroid, the residual error was ≤ 0.18 units of Z_eff (vs. 0.2 with the Gammex 467 phantom from a previous study). These results further suggest that a dual-energy algorithm in combination with stoichiometry can decrease errors in calculation of the ρ_e of real tissues to ameliorate inhomogeneity for dose calculation in radiotherapy treatment planning, especially when the energy spectrum of the X-ray tube of the CT machine is not available.

Ethylene-vinyl acetate foam as a new lung substitute in radiotherapy

PURPOSE

The purpose of this study was to evaluate ethylene-vinyl acetate (EVA) foam as a new lung substitute in radiotherapy and to study its physical and dosimetric characteristics.

METHODS

We calculated the ideal vinyl acetate (VA) content of EVA foam sheets to mimic the physical and dosimetric characteristics of the ICRU lung tissue. We also computed the water-to-medium mass collision stopping power ratios, mass attenuation coefficients, CT numbers, effective atomic numbers and electron densities for: ICRU lung tissue, the RANDO commercial phantom, scaled WATER and EVA foam sheets with varying VA contents in a range between the minimum and maximum values supplied by the manufacturer. For all these substitutes, we simulated percent depth-dose curves with EGSnrc Monte Carlo (MC PDDs) in a water-lung substitute-water slab phantom expressed as dose-to-medium and dose-to-water for 3 × 3- and 10 × 10-cm² field sizes. PDD for the 10 × 10-cm² field size was also calculated with the MultiGrid Superposition algorithm (MGS PDD) for a relative electron density to water ratio of 0.26. The latter was compared with the MC PDDs in dose-to-water for scaled WATER and EVA foam sheets with the VA content that was most similar to the calculated ideal content that is physically achievable in practice.

RESULTS

We calculated an ideal VA content of 55%; however, the maximum physically achievable content with current manufacturing techniques is 40%. The physical characteristics of the EVA foam sheets with a VA content of 40% (EVA40) are very close to those of the ICRU lung reference. The physical densities of the EVA40 foam sheets ranged from 0.030 to 0.965 g/cm³ , almost covering the entire physical density range of the inflated/deflated lung (0.260-1.050 g/cm³ ). Its mass attenuation coefficient at the effective energy of a 6-MV photon beam agrees within 0.8% of the ICRU reference value, and its CT number agrees within 6 HU. The effective atomic number for EVA40 varies by less than 0.42 of the ICRU value, and its effective electron density is within 0.9%. PDDs expressed in dose-to-medium and dose-to-water agree with the ICRU curve within 2% in all regions. PDDs calculated with both MC and MGS were within 1.5%.

CONCLUSIONS

The EVA40 is an excellent cork-like lung substitute for radiotherapy applications. From a sole material used in footwear, it is possible to obtain a lung substitute that mimics the physical and dosimetric characteristics of ICRU lung tissue even better than the RANDO commercial phantom.

A study on a dental device for the prevention of mucosal dose enhancement caused by backscatter radiation from dental alloy during external beam radiotherapy

The changes in dose distribution caused by backscatter radiation from a common commercial dental alloy (Au-Ag-Pd dental alloy; DA) were investigated to identify the optimal material and thicknesses of a dental device (DD) for effective prevention of mucositis. To this end, 1 cm³ of DA was irradiated with a 6-MV X-ray beam (100 MU) in a field size of 10 × 10 cm² using a Novalis TX linear accelerator. Ethylene vinyl acetate copolymer, polyolefin elastomer, and polyethylene terephthalate (PET) were selected as DD materials. The depth dose along the central axis was determined with respect to the presence/absence of DA and DDs at thicknesses of 1-10 mm using a parallel-plate ionization chamber. The dose in the absence of DDs showed the lowest value at a distance of 5 mm from the DA surface and gradually increased with distance between the measurement point and the DA surface for distances of ≥5 mm. Except for PET, no significant difference between the DA dose curves for the presence and absence of DDs was observed. In the dose curve, PET showed a slightly higher dose for DA with DD than for DA without DD for thicknesses of ≥4 mm. The findings herein suggest that the optimal DD material for preventing local dose enhancement of the mucosa caused by DA backscatter radiation should have a relatively low atomic number and physical density and that optimal DD thickness should be chosen considering backscatter radiation and percentage depth dose.

Evaluation of various boluses in dose distribution for electron therapy of the chest wall with an inward defect

Delivering radiotherapy to the postmastectomy chest wall can be achieved using matched electron fields. Surgical defects of the chest wall change the dose distribution of electrons. In this study, the improvement of dose homogeneity using simple, nonconformal techniques of thermoplastic bolus application on a defect is evaluated. The proposed phantom design improves the capability of film dosimetry for obtaining dose profiles of a patient's anatomical condition. A modeled electron field of a patient with a postmastectomy inward surgical defect was planned. High energy electrons were delivered to the phantom in various settings, including no bolus, a bolus that filled the inward defect (PB0), a uniform thickness bolus of 5 mm (PB1), and two 5 mm boluses (PB2). A reduction of mean doses at the base of the defect was observed by any bolus application. PB0 increased the dose at central parts of the defect, reduced hot areas at the base of steep edges, and reduced dose to the lung and heart. Thermoplastic boluses that compensate a defect (PB0) increased the homogeneity of dose in a fixed depth from the surface; adversely, PB2 increased the dose heterogeneity. This study shows that it is practical to investigate dose homogeneity profiles inside a target volume for various techniques of electron therapy.

Dose correction in lung for HDR breast brachytherapy

Purpose

To evaluate the dosimetric impact of lung tissue in Ir-192 APBI.

Material and methods

In a 40 × 40 × 40 cm³ water tank, an Accelerated Partial Breast Irradiation (APBI) brachytherapy balloon inflated to 4 cm diameter was situated directly below the center of a 30 × 30 × 1 cm³ solid water slab. Nine cm of solid water was stacked above the 1 cm base. A parallel plate ion chamber was centered above the base and ionization current measurements were taken from the central HDR source dwell position for channels 1, 2, 3 and 5 of the balloon. Additional ionization data was acquired in the 9 cm stack at 1 cm increments. A comparable data set was also measured after replacing the 9 cm solid water stack with cork slabs. The ratios of measurements in the two phantoms were calculated and compared to predicted results of a commercial treatment planning system.

Results

Lower dose was measured in the cork within 1 cm of the cork/solid water interface possibly due to backscatter effects. Higher dose was measured beyond 1 cm from the cork/solid water interface, increasing with path length up to 15% at 9 cm depth in cork. The treatment planning system did not predict either dose effect.

Conclusions

This study investigates the dosimetry of low density material when the breast is treated with Ir-192 brachytherapy. HDR dose from Ir-192 in a cork media is shown to be significantly different than in unit density media. These dose differences are not predicted in most commercial brachytherapy planning systems. Empirical models based on measurements could be used to estimate lung dose associated with HDR breast brachytherapy.