Converting absorbed dose to medium to absorbed dose to water for Monte Carlo based photon beam dose calculations

Experimental validation of absorbed dose-to-medium calculation algorithms in heterogeneous media

Objective.The aim of this work was to determine heterogeneous correction factorshQclin,Qreffclin,frefdetm,wto validate absorbed dose-to-mediumDm,Qclinm,fclincalculation algorithms from detector readings. The impact of detector orientation perpendicular and parallel to the beam central axis on the correction factors was also investigated.Approach.ThehQclin,Qreffclin,frefdetm,wfactors were calculated for four types of detectors (PTW PinPoint T31016, PTW microDiamond T60019, PTW microSilicon T60023 and EBT3 film) placed in different media (cortical bone, lung, adipose tissue, Teflon and RW3) for the 6 MV energy beam with a 10 × 10 cm2field size. These corrections were then applied to the detector measurements performed at different depths in heterogeneous phantoms.Main results.ThehQclin,Qreffclin,frefdetm,wfactors mainly depended on the media and slightly on the type of detector. Considering all detectors, the largest corrections were found in high-density media with values ranging from 0.911 to 0.934 in cortical bone. For comparison, the corrections in other media were closer to unity with values from 0.966 (lung and RW3) to 0.991 (adipose tissue). Except for the PinPoint T31016, detector orientation-dependence was observed especially in high-density media. A good agreement (≤1.5%) was found betweenDm,Qclinm,fclincalculations and the detector readings corrected with thehQclin,Qreffclin,frefdetm,wfactor for all studied heterogeneous phantoms.Significance.This paper could serve as an initial guideline for medical physicists involved in the validation of the advanced type-b dose calculation algorithms reportingDm,Qclinm,fclin. To our knowledge, this is the first study to assess the impact of the orientation of different detectors in heterogeneous media. The orientation dependence of the detector response observed in water may not reflect what is observed in heterogeneous media, especially in high-density media. The knowledge of thehQclin,Qreffclin,frefdetm,wfactors becomes mandatory for accurate interpretation of detector readings and comparisons withDm,Qclinm,fclincalculations.

Adjuvant Radiation Treatment of Breast Cancer After Mastectomy: Advanced Algorithms and Partial Bolus Improve the Dose Calculation Accuracy in the Case of Thin-Chest-Wall Irradiation

Purpose

The aim of this study was to examine measured and calculated dose distributions in a thin-chest-wall phantom and estimate the variations in the dose-volume histogram (DVH) parameters used in plan evaluation for patient geometries with chest-wall thicknesses <15 mm with and without bolus implementation.

Methods and Materials

Measurements were made using thermoluminescent dosimeters in a chest-wall phantom. The Monte Carlo method, anisotropic analytical algorithm, and Acuros XB Eclipse algorithms were used to calculate dose distributions for clinical plans. DVH parameters for clinical target volume tumor (CTVT) and planning target volume (PTV) and mean doses were evaluated for 15 patients with a chest-wall thickness of 8 to 15 mm with and without partial bolus and for 10 patients with a chest-wall thickness of 20 to 25 mm without bolus.

Results

Measurements showed that the dose at a depth of 2 to 12 mm at the beam entrance and laterally was within 90% of the dose at 8 mm depth. Monte Carlo and Acuros XB calculations were well aligned with the experimental data, whereas the anisotropic analytical algorithm underestimated the beam entrance and lateral doses. The DVH parameters for the patients with a thin chest wall were sensitive to calculation algorithm, resolution, body structure definition, and patient geometry. The parameters CTVT_V95%, CTVT_D98%, and PTV_D98% were much lower than the tolerance criteria. Partial bolus improved the values for all algorithms and decreased the variations due to patient geometry. Dose calculations for patients with a chest-wall thickness of 20 to 25 mm resulted in sufficient target coverage and low dependence on patient geometry and calculation algorithm without the use of bolus.

Conclusions

Dose calculations using advanced algorithms and resolution <2 mm are recommended for patients with a thin chest wall. Specific DVH criteria or the implementation of partial bolus was needed to facilitate plan development and evaluation for this patient group.

Investigation of a water equivalent depth method for dosimetric accuracy evaluation of synthetic CT

PURPOSE

To provide a metric that reflects the dosimetric utility of the synthetic CT (sCT) and can be rapidly determined.

METHODS

Retrospective CT and atlas-based sCT of 62 (53 IMRT and 9 VMAT) prostate cancer patients were used. For image similarity measurements, the sCT and reference CT (rCT) were aligned using clinical registration parameters. Conventional image similarity metrics including the mean absolute error (MAE) and mean error (ME) were calculated. The water equivalent depth (WED) was automatically determined for each patient on the rCT and sCT as the distance from the skin surface to the treatment plan isocentre at 36 equidistant gantry angles, and the mean WED difference (ΔWED¯) between the two scans was calculated. Doses were calculated on each scan pair for the clinical plan in the treatment planning system. The image similarity measurements and ΔWED¯ were then compared to the isocentre dose difference (ΔD_iso) between the two scans.

RESULTS

While no particular relationship to dose was observed for the other image similarity metrics, the ME results showed a linear trend against ΔD_iso with R² = 0.6, and the 95 % prediction interval for ΔD_iso between -1.2 and 1 %. The ΔWED¯ results showed an improved linear trend (R² = 0.8) with a narrower 95 % prediction interval from -0.8 % to 0.8 %.

CONCLUSION

ΔWED¯ highly correlates with ΔD_iso for the reference and synthetic CT scans. This is easy to calculate automatically and does not require time-consuming dose calculations. Therefore, it can facilitate the process of developing and evaluating new sCT generation algorithms.

Calorimeter measurements of absolute dose in aluminum, a surrogate of bone, to validate dose-to-medium in Acuros XB

Objective. While the accuracy of dose calculations in water with Acuros XB is well established, experimental validation of dose in bone is limited. Acuros XB reports both dose-to-medium and dose-to-water, and these values differ in bone, but there are no reports of measurements of validation in bone. This work compares Acuros XB calculations to measurements of absolute dose in aluminum (medium similar to bone). The validity of using selected relative dosimeters in aluminum is also investigated.Approach. A calorimeter with an aluminum core embedded in an aluminum phantom was selected as bone surrogate for the measurement of absolute dose. Matching the medium of the core to the medium of the phantom allowed eliminating the calculation of the conversion between media. The dose was measured at the fixed depth of 3.3 cm in aluminum (∼9 g·cm^-2) with 6X, 10X, 6FFF and 10FFF photon beams from a TrueBeam Varian linac. In addition, experimental cross-calibration between water and aluminum was performed for an IBA CC13 ionization chamber, a PTW microDiamond and EBT3 Gafchromic film.Main results. Calculations with Acuros XB dose-to-medium in aluminum differed from the calorimetry data by -2.8% to -3.5%, depending on the beam. Use of dose-to-water would have resulted in about 39% discrepancy. The cross calibration coefficient between water and aluminum yielded values of about 0.87 for the CC13 chamber, 0.91 for the microDiamond, and 0.88 for the film, and independent of the beam within about ±1%.Significance. It was demonstrated the value of the dose-to-medium in aluminum (surrogate of bone) computed with Acuros XB is close to the value of the absolute dose measured with a calorimeter, and there is a significant discrepancy when dose-to-water is used instead. The use of an ionization chamber, a microDiamond and Gafchromic film in aluminum required a considerable correction from calibration in water.

Dosimetric comparison of analytical anisotropic algorithm and the two dose reporting modes of Acuros XB dose calculation algorithm in volumetric modulated arc therapy of carcinoma lung and carcinoma prostate

Volumetric Modulated Arc Therapy (VMAT) is an important modality for radical radiotherapy of all major treatment sites. This study aims to compare Analytical Anisotropic Algorithm (AAA) and the two dose-reporting modes of Acuros XB (AXB) algorithm -the dose to medium option (D_m) and the dose to water option (D_w) in Volumetric Modulated Arc Therapy (VMAT) of carcinoma lung and carcinoma prostate. We also compared the measured dose with Treatment Planning System calculated dose for AAA and the two dose reporting options of Acuros XB using Electronic Portal Imaging Device (EPID) and ArcCHECK phantom. Treatment plans of twenty patients each who have already undergone radiotherapy for cancer of lung and cancer of prostate were selected for the study. Three sets of VMAT plans were generated in Eclipse Treatment Planning System (TPS), one with AAA and two plans with Acuros-D_m and Acuros-D_w options. The Dose Volume Histograms (DVHs) were compared and analyzed for Planning Target Volume (PTV) and critical structures for all the plans. Verification plans were created for each plan and measured doses were compared with TPS calculated doses using EPID and ArcCHECK phantom for all the three algorithms. For lung plans, the mean dose to PTV in the AXB-D_w plans was higher by 1.7% and in the AXB-D_m plans by 0.66% when compared to AAA plans. For prostate plans, the mean dose to PTV in the AXB-D_w plans was higher by 3.0% and in the AXB-D_m plans by 1.6% when compared to AAA plans. There was no difference in the Conformity Index (CI) between AAA and AXB-D_m and between AAA and AXB-D_w plans for both sites. But the homogeneity worsened in AXB-D_w and AXB-D_m plans when compared to AAA plans for both sites. AXB-D_w calculated higher dose values for PTV and all the critical structures with significant differences with one or two exceptions. Point dose measurements in ArcCHECK phantom showed that AXB-D_m and AXB-D_w options showed very small deviations with measured dose distributions than AAA for both sites. Results of EPID QA also showed better pass rates for AXB-D_w and AXB-D_m than AAA for both sites when gamma analysis was done for 3%/3 mm and 2%/2 mm criteria. With reference to the results, it is always better to choose Acuros algorithm for dose calculations if it is available in the TPS. AXB-D_w plans showed very high dose values in the PTV when compared to AAA and AXB-D_m in both sites studied. Also, the volume of PTV receiving 107% dose was significantly high in AXB-D_w plans compared to AXB-D_m plans in sites involving high density bones. Considering the results of dosimetric comparison and QA measurements, it is always better to choose AXB-D_m algorithm for dose calculations for all treatment sites especially when high density bony structures and complex treatment techniques are involved. For patient specific QA purposes, choosing AXB-D_m or AXB-D_w does not make any significant difference between calculated and measured dose distributions.

Treatment Planning Study for Microbeam Radiotherapy Using Clinical Patient Data

Microbeam radiotherapy (MRT) is a novel, still preclinical dose delivery technique. MRT has shown reduced normal tissue effects at equal tumor control rates compared to conventional radiotherapy. Treatment planning studies are required to permit clinical application. The aim of this study was to establish a dose comparison between MRT and conventional radiotherapy and to identify suitable clinical scenarios for future applications of MRT. We simulated MRT treatment scenarios for clinical patient data using an inhouse developed planning algorithm based on a hybrid Monte Carlo dose calculation and implemented the concept of equivalent uniform dose (EUD) for MRT dose evaluation. The investigated clinical scenarios comprised fractionated radiotherapy of a glioblastoma resection cavity, a lung stereotactic body radiotherapy (SBRT), palliative bone metastasis irradiation, brain metastasis radiosurgery and hypofractionated breast cancer radiotherapy. Clinically acceptable treatment plans were achieved for most analyzed parameters. Lung SBRT seemed the most challenging treatment scenario. Major limitations comprised treatment plan optimization and dose calculation considering the tissue microstructure. This study presents an important step of the development towards clinical MRT. For clinical treatment scenarios using a sophisticated dose comparison concept based on EUD and EQD2, we demonstrated the capability of MRT to achieve clinically acceptable dose distributions.

The effect on tumour control probability of using AXB algorithm in replacement of AAA for SBRT of hepatocellular carcinoma located at lung-liver boundary region

Objective:

To retrospectively analyze the clinical impact on stereotactic body radiation therapy (SBRT) for hepatocellular carcinoma (HCC) located at lung–liver boundary due to the use of Acuros XB algorithm (AXB) in replacement of anisotropic analytical algorithm (AAA).

Methods:

23 SBRT volumetric modulated arc therapy (VMAT) plans for HCC located at lung–liver boundary were calculated using AAA and AXB respectively with the same treatment parameters. The dose–volume data of the planned target volumes (PTVs) were compared. A published tumour control probability (TCP) model was used to calculate the effect of dosimetric difference between AAA and AXB on tumour control probability.

Results:

For dose calculated by AXB (Dose to medium), the D95% and D98% of the PTV were on average 2.4 and 3.1% less than that calculated by AAA. For dose calculated by AXB (dose to water), the D95% and D98% of the PTV were on average 1.8%, and 2.7% less than that calculated by AAA. Up to 5% difference in D95% and 8% difference in D98% were observed in the worst cases. The significant decrease in D95% calculated by AXB compared to AAA could result in a % decrease in 2 year TCP up to 8% in the worst case (from 46.8 to 42.9%).

Conclusion:

The difference in dose calculated by AAA and AXB could lead to significant difference in TCP for HCC SBRT located at lung–liver boundary region.

Advances in knowledge:

The difference in calculated dose and tumour control probability for HCC SBRT between AAA and AXB algorithm at lung–liver boundary region was compared.

Impact of the dose quantity used in MV photon optimization on dose distribution, robustness, and complexity

PURPOSE

Convolution/superposition algorithms used in megavoltage (MV) photon radiotherapy model radiation transport in water, yielding dose to water-in-water (D_w,w ). Advanced algorithms constitute a step forward, but their dose distributions in terms of dose to medium-in-medium (D_m,m ) or dose to water-in-medium (D_w,m ) can be problematic when used in plan optimization due to their different dose responses to some atomic composition heterogeneities. Failure to take this into account can lead to undesired overcorrections and thus to unnoticed suboptimal and unrobust plans. Dose to reference-like medium (D_ref,m* ) was recently introduced to overcome these limitations while ensuring accurate transport. This work evaluates and compares the performance of these four dose quantities in planning target volume (PTV)-based optimization.

METHODS

We considered three cases with heterogeneities inside the PTV: virtual phantom with water surrounded by bone; head and neck; and lung. These cases were planned with volumetric modulated arc therapy (VMAT) technique, optimizing with the same setup and objectives for each dose quantity. We used different algorithms of the Varian Eclipse treatment planning system (TPS): Acuros XB (AXB) for D_m,m and D_w,m , and Analytical Anisotropic Algorithm (AAA) for D_w,w . D_ref,m* was obtained from D_m,m distributions using an in-house software considering water as the reference medium (D_w,m* ). The optimization process consisted of: (1) common first optimization, (2) dose distribution computed for each quantity, (3) re-optimization, and (4) final calculation for each dose quantity. The dose distribution, robustness to patient setup errors, and complexity of the plans were analyzed and compared.

RESULTS

The quantities showed similar dose distributions after the optimization but differed in terms of plan robustness. The cases with soft tissue and high-density heterogeneities followed the same pattern. For AXB D_m,m , cold regions appeared in the heterogeneities after the first optimization. They were compensated in the second optimization through local fluence increases, but any positional mismatch impacted robustness, with clinical target volume (CTV) variations from the nominal scenario around +3% for bone and up to +7% for metal. For AXB D_w,m the pattern was inverse (hot regions compensated by fluence decreases) and more pronounced, with CTV dose variations around -7% for bone and up to -17% for metal. Neither AXB D_w,m* nor AAA D_w,w presented these dose inhomogeneities, which resulted in more robust plans. However, D_w,w differed markedly from the other quantities in the lung case because of its lower radiation transport accuracy. AXB D_m,m was the most complex of the four dose quantities and AXB D_w,m* the least complex, though we observed no major differences in this regard.

CONCLUSIONS

The dose quantity used in MV photon optimization can affect plan robustness. D_w,w distributions from convolution/superposition algorithms are robust but may not provide sufficient radiation transport accuracy in some cases. D_m,m and D_w,m from advanced algorithms can compromise robustness because their different responses to some composition heterogeneities introduce additional fluence compensations. D_ref,m* offers advantages in plan optimization and evaluation, producing accurate and robust plans without increasing complexity. D_ref,m* can be easily implemented as a built-in feature of the TPS and can facilitate and simplify the treatment planning process when using advanced algorithms. Final reporting can be kept in D_m,m or D_w,m for clinical correlations.

Effects of lung tissue characterization in radiotherapy of breast cancer under deep inspiration breath hold when using Monte Carlo dosimetry

PURPOSE

To investigate the sensitivity of Monte Carlo (MC) calculated lung dose distributions to lung tissue characterization in external beam radiotherapy of breast cancer under Deep Inspiration Breath Hold (DIBH).

METHODS

EGSnrc based MC software was employed. Mean lung densities for one hundred patients were analysed. CT number frequency and clinical dose distributions were calculated for 15 patients with mean lung density below 0.14 g/cm³. Lung volume with a pre-defined CT numbers was also considered. Lung tissue was characterized by applying different CT calibrations in the low-density region and air-lung tissue thresholds. Dose impact was estimated by Dose Volume Histogram (DVH) parameters.

RESULTS

Mean lung densities below 0.14 g/cm³ were found in 10% of the patients. CT numbers below -960 HU dominated the CT frequency distributions with a high rate of CT numbers at -990 HU. Mass density conversion approach influenced the DVH shape. V_4Gy and V_8Gy varied by 7% and 5% for the selected patients and by 9% and 3.5% for the pre-defined lung volume. V_16Gy and V_20Gy, were within 2.5%. Regions above 20 Gy were affected. Variations in air- lung tissue differentiation resulted in DVH parameters within 1%. Threshold at -990 HU was confirmed by the CT number frequency distributions.

CONCLUSIONS

Lung dose distributions were more sensitive to variations in the CT calibration curve below lung (inhale) density than to air-lung tissue differentiation. Low dose regions were mostly affected. The dosimetry effects were found to be potentially important to 10% of the patients treated under DIBH.

Report dose-to-medium in clinical trials where available; a consensus from the Global Harmonisation Group to maximize consistency

PURPOSE

To promote consistency in clinical trials by recommending a uniform framework as it relates to radiation transport and dose calculation in water versus in medium.

METHODS

The Global Quality Assurance of Radiation Therapy Clinical Trials Harmonisation Group (GHG; www.rtqaharmonization.org) compared the differences between dose to water in water (D_w,w), dose to water in medium (D_w,m), and dose to medium in medium (D_m,m). This was done based on a review of historical frameworks, existing literature and standards, clinical issues in the context of clinical trials, and the trajectory of radiation dose calculations. Based on these factors, recommendations were developed.

RESULTS

No framework was found to be ideal or perfect given the history, complexity, and current status of radiation therapy. Nevertheless, based on the evidence available, the GHG established a recommendation preferring dose to medium in medium (D_m,m).

CONCLUSIONS

Dose to medium in medium (D_m,m) is the preferred dose calculation and reporting framework. If an institution's planning system can only calculate dose to water in water (D_w,w), this is acceptable.

Evaluation of clinical implications in the use of dose to water versus dose to medium by using NTCP and TCP models for urinary bladder tumours

Purpose: To analyze the dosimetric and radiobiological differences between dose to water versus dose to medium for patients with carcinoma of the urinary bladder.

Materials and Methods: 15 patients with cancer of urinary bladder were selected for the study. VMAT plans were generated for each patient. The dose distributions were calculated in the modes dose to water and to medium with the Monaco treatment planning system. A dosimetric comparative analysis has been made between the two modes of planning in this study. Subsequently, NTCP and TCP were determined for OARs and targets respectively.

Results: The mean dose to 2 cc of the rectum, small bowel, left and right femoral heads respectively was higher by 0.8, 1.2, 2.7, and 2.2% for the dose to water calculation. Similarly, the mean dose to D2, D50, and D98 for PTV was higher by 0.4, 0.3, and 0.3% for dose to water calculation. Such small dose differences had little effect on the values of TCP and NTCP.

Conclusion: For patients with the urinary bladder there were very small differences between results between calculations carried out in dose to medium and dose to water modes.

Towards the standardization of the absorbed dose report mode in high energy photon beams

The benefits of using an algorithm that reports absorbed dose-to-medium have been jeopardized by the clinical experience and the experimental protocols that have mainly relied on absorbed dose-to-water. The aim of the present work was to investigate the physical aspects that govern the dosimetry in heterogeneous media using Monte Carlo method and to introduce a formalism for the experimental validation of absorbed dose-to-medium reporting algorithms. Particle fluence spectra computed within the sensitive volume of two simulated detectors (T31016 Pinpoint 3D ionization chamber and EBT3 radiochromic film) placed in different media (water, RW3, lung and bone) were compared to those in the undisturbed media for 6 MV photon beams. A heterogeneity correction factor that takes into account the difference between the detector perturbation in medium and under reference conditions as well as the stopping-power ratios was then derived for all media using cema calculations. Furthermore, the different conversion approaches and Eclipse treatment planning system algorithms were compared against the Monte Carlo absorbed dose reports. The detectors electron fluence perturbation in RW3 and lung media were close to that in water (≤1.5%). However, the perturbation was greater in bone (∼4%) and impacted the spectral shape. It was emphasized that detectors readings should be corrected by the heterogeneity correction factor that ranged from 0.932 in bone to 0.985 in lung. Significant discrepancies were observed between all the absorbed dose reports and conversions, especially in bone (exceeding 10%) and to a lesser extent in RW3. Given the ongoing advances in dose calculation algorithms, it is essential to standardize the absorbed dose report mode with absorbed dose-to-medium as a favoured choice. It was concluded that a retrospective conversion should be avoided and switching from absorbed dose-to-water to absorbed dose-to-medium reporting algorithm should be carried out by a direct comparison of both algorithms.

A new dose quantity for evaluation and optimisation of MV photon dose distributions when using advanced algorithms: proof of concept and potential applications

Advanced algorithms used in MV photon radiotherapy model radiation transport in any media. They represent a step forward but introduce new uncertainties and questions, including whether to report the doses to water (D_w,m) or medium (D_m,m) voxels, and the impact of fluence changes introduced by surrounding media. These aspects can compromise consistency between both reporting modes and with previous algorithms in which clinical experience is based. This study introduces a new dose quantity, the dose-to-reference-like medium, to overcome the aforementioned shortcomings. It is linked to a reference medium, water in this study (D_w,m*), and defined as the absorbed dose in a voxel of this reference medium surrounded by a reference-like medium with the same radiation transport characteristics as the original one. We propose to derive D_w,m* distributions by post-processing D_w,m or D_m,m applying a correction factor (CF) to each voxel which depends on its composition. We present and justify a simple and straightforward method to obtain CFs that only involves two phantoms with the same density: one with the considered composition and the other with that of the reference medium. A proof of concept was performed in a clinical environment for Acuros XB (AXB) advanced algorithm and 6 MV photon beams. The CFs were derived for the tissues characterised in Acuros. D_w,m* was compared to D_w,m, D_m,m, and D_w,w from AAA analytical algorithm for some virtual and clinical cases. All the previous quantities presented limitations that can be solved by D_w,m*. This new quantity allows the applicability of evaluation parameters, traceability to clinical experience, and isolation of heterogeneity effects to identify optimum plans, offering useful characteristics for plan evaluation and optimisation in clinical practice. Additionally, it also has potential applications in automated treatment planning and multi-centre activities such as clinical trials, audits, benchmarking, and shared models for automation.

Feasibility of using a commercial collapsed cone dose engine for 1.5T MR-LINAC online independent dose verification

PURPOSE

To validate the feasibility and accuracy of commonly used collapsed cone (CC) dose engine for Elekta Unity 1.5T MR-LINAC online independent dose verification.

MATERIALS AND METHODS

The Unity beam model was built and commissioned in RayStation treatment planning system with CC dose engine. Four AAPM TG-119 test plans were created and measured with ArcCHECK phantom for comparison, another thirty patient plans from six tumor sites were also included. The dosimetric criteria for various ROIs and 3D gamma passing rates were quantitatively evaluated, and the effects of magnetic field and dose deposition type on the dose difference between two systems were further analyzed.

RESULTS

ArcCHECK based measurement showed a clear magnetic field induced profile shift between CC with both measurement and GPUMCD. For clinical plans, gamma passing rates with criteria (3%, 3 mm) between GPUMCD and CC large than 90% can be achieved for most tumor sites except esophagus and lung cases, the mean dose difference of 3% can be satisfied for most ROIs from all tumor sites. The magnetic field caused a large dose impact on low density areas, the average gamma passing rates were improved from 85.54% to 96.43% and 87.40% to 99.54% for esophagus and lung cases when the magnetic field effect was excluded.

CONCLUSIONS

It is feasible to use CC dose engine as a secondary dose calculation tool for Elekta Unity system for most tumor sites, while the accuracy is limited and should be used carefully for low density areas, such as esophagus and lung cases.

Monte Carlo simulation of 6-MV dynamic wave VMAT deliveries by Vero4DRT linear accelerator using EGSnrc moving sources

The commissioning and benchmark of a Monte Carlo (MC) model of the 6‐MV Brainlab‐Mitsubishi Vero4DRT linear accelerator for the purpose of quality assurance of clinical dynamic wave arc (DWA) treatment plans is reported. Open‐source MC applications based on EGSnrc particle transport codes are used to simulate the medical linear accelerator head components. Complex radiotherapy irradiations can be simulated in a single MC run using a shared library format combined with BEAMnrc “source20.” Electron energy tuning is achieved by comparing measured vs simulated percentage depth doses (PDDs) for MLC‐defined field sizes in a water phantom. Electron spot size tuning is achieved by comparing measured and simulated inplane and crossplane beam profiles. DWA treatment plans generated from RayStation (RaySearch) treatment planning system (TPS) are simulated on voxelized (2.5 mm³) patient CT datasets. Planning target volume (PTV) and organs at risk (OAR) dose–volume histograms (DVHs) are compared to TPS‐calculated doses for clinically deliverable dynamic volumetric modulated arc therapy (VMAT) trajectories. MC simulations with an electron beam energy of 5.9 MeV and spot size FWHM of 1.9 mm had the closest agreement with measurement. DWA beam deliveries simulated on patient CT datasets results in DVH agreement with TPS‐calculated doses. PTV coverage agreed within 0.1% and OAR max doses (to 0.035 cc volume) agreed within 1 Gy. This MC model can be used as an independent dose calculation from the TPS and as a quality assurance tool for complex, dynamic radiotherapy treatment deliveries. Full patient CT treatment simulations are performed in a single Monte Carlo run in 23 min. Simulations are run in parallel using the Condor High‐Throughput Computing software¹ on a cluster of eight servers. Each server has two physical processors (Intel Xeon CPU E5‐2650 0 @2.00 GHz), with 8 cores per CPU and two threads per core for 256 calculation nodes.

Dosimetric impact of switching from AAA to Acuros dose-to-water and dose-to-medium for RapidArc plans of nasopharyngeal carcinomas

PURPOSE

This work aims to evaluate the dosimetric consequences of replacing the Anisotropic Analytical Algorithm (AAA) by Acuros XB (AXB), dose-to-water (Dw) or dose-to-medium (Dm), for RapidArc plans of nasopharyngeal carcinomas (NPC).

MATERIALS AND METHODS

Seventeen NPC plans created with AAA (v15.6) were recalculated with AXB (v15.6) Dw and Dm. The dose-volume parameters to the planning target volumes (PTV) and relevant organs at risk (OAR) were compared. The high dose PTV was divided into bone, air and tissue components and the comparison was performed for each of them.

RESULTS

AXB Dw revealed no significant differences in the PTVs compared to AAA. Lower values were observed to spinal cord, brainstem, oral cavity and parotids (0.5% to 2.3%), and higher values to cochleas (up to 5.4%) and mandible (up to 6.7%). AXB Dm predicted lower values than AAA for all PTVs and OARs (2.0% to 6.1%). For the bone PTV subvolume, AXB Dw and Dm predicted respectively higher (2.4%) and lower (2.2% to 3.4%) values. No significant differences were noted in air. AXB predicted lower values than AAA in soft tissues (0.4% to 1.6%). The largest difference was found to the mandible V_60Gy parameter, with median differences of 6.7% for AXB Dw and -6.0% for AXB Dm.

CONCLUSION

Significant dose differences are expected when switching from AAA to AXB in NPC cases. The dose prescriptions and the tolerance limits for some OARs, especially those of high density, may need to be adjusted depending on the selected dose calculation algorithm and reporting mode.

Dose Calculation Comparisons between Three Modern Treatment Planning Systems

Purpose:

Monaco treatment planning system (TPS) version 5.1 uses a Monte-Carlo (MC)-based dose calculation engine. The aim of this study is to verify and compare the Monaco-based dose calculations with both Pinnacle³ collapsed cone convolution superposition (CCCS) and Eclipse anisotropic analytical algorithm (AAA) calculations.

Materials and Methods:

For this study, 18 previously treated lung and head-and-neck (HN) cancer patients were chosen to compare the dose calculations between Pinnacle, Monaco, and Eclipse. Plans were chosen from those that had been treated using the Elekta VersaHD or a Novalis Tx linac. All of the treated volumetric-modulated arc therapy plans used 6 MV or 10 MV photon beams. The original plans calculated with CCCS or AAA along with the recalculated ones using MC from the three TPS were exported into Velocity software for intercomparison.

Results:

To compare the dose calculations, Planning target volume (PTV) heterogeneity indexes and conformity indexes were calculated from the dose volume histograms (DVH) of all plans. While mean lung dose (MLD), lung V5 and V20 values were recorded for lung plans, the computed dose to parotids, brainstem, and mandible were documented for HN plans. In plan evaluation, percent differences of the above dosimetric values in Monaco computation were compared against each of the other TPS computations.

Conclusion:

It could be concluded through this research that there can be differences in the calculation of dose across different TPSs. Although relatively small, these differences could become apparent when compared using DVH. These differences most likely arise from the different dose calculation algorithms used in each TPS. Monaco employs the MC allowing it to have much more detailed calculations that result in it being seen as the most accurate and the gold standard.

Experimental validation of Monte Carlo based treatment planning system in bone density equivalent media

Introduction Advanced, Monte Carlo (MC) based dose calculation algorithms, determine absorbed dose as dose to medium-in-medium (Dm,m) or dose to water-in-medium (Dw,m). Some earlier studies identified the differences in the absorbed doses related to the calculation mode, especially in the bone density equivalent (BDE) media. Since the calculation algorithms built in the treatment planning systems (TPS) should be dosimetrically verified before their use, we analyzed dose differences between two calculation modes for the Elekta Monaco TPS. We compared them with experimentally determined values, aiming to define a supplement to the existing TPS verification methodology. Materials and methods In our study, we used a 6 MV photon beam from a linear accelerator. To evaluate the accuracy of the TPS calculation approaches, measurements with a Farmer type chamber in a semi-anthropomorphic phantom were compared to those obtained by two calculation options. The comparison was made for three parts of the phantom having different densities, with a focus on the BDE part. Results Measured and calculated doses were in agreement for water and lung equivalent density materials, regardless of the calculation mode. However, in the BDE part of the phantom, mean dose differences between the calculation options ranged from 5.7 to 8.3%, depending on the method used. In the BDE part of the phantom, neither of the two calculation options were consistent with experimentally determined absorbed doses. Conclusions Based on our findings, we proposed a supplement to the current methodology for the verification of commercial MC based TPS by performing additional measurements in BDE material.

Clinical implementation of a Monte Carlo based independent TPS dose checking system

The increase in complexity of treatment plans over time through modalities such as intensity modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT) has often not been met with an increase in capability of the secondary dose calculation checking systems typically used to verify the treatment planning system. Monte Carlo (MC) codes such as EGSnrc have become easily available and are capable of performing calculations of highly complex radiotherapy treatments. This educational note demonstrates a method for implementing and using a fully automated system for performing and analysing full MC calculations of conformal, IMRT and VMAT radiotherapy plans. Example calculations were based on BEAMnrc/DOSXYZnrc and are performed automatically after either uploading exported plan DICOM data through a Python-based web interface, or exporting DICOM data to a monitored network location. This note demonstrates how completed MC calculations can then be analysed using an automatically generated dose point comparison report, or easily re-imported back into the treatment planning system. Agreement between the TPS and MC calculation was an improvement on agreement between RadCalc and the TPS, with differences ranging from 1.2 to 5.5% between RadCalc and the treatment planning system (TPS), and 0.1-1.7% between MC and TPS. Comparison of the dose-volume histogram (DVH) parameters [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text] for the example VMAT plans showed agreement for the mean planning target volume dose within [Formula: see text], [Formula: see text] and [Formula: see text] generally within [Formula: see text] with the exception of a brain case, and [Formula: see text] within [Formula: see text]. Overall, this note provides a demonstration of a system that has been integrated well into existing clinical workflow, and has been shown to be a valuable additional tool in the secondary checking of treatment plan calculations.

Clinical implication in the use of the AAA algorithm versus the AXB in nasopharyngeal carcinomas by comparison of TCP and NTCP values

Purpose

Retrospective analysis of volumetric modulated arc therapy treatment plans to investigate qualitative, possible, clinical consequences of the use of AAA versus AXB in nasopharyngeal cancer (NPC) cases.

Methods

The dose distribution of 26 treatment plans, produced using RapidArc technique and AAA algorithm, were recalculated using AXB and the same number of monitor units provided by AAA and clinically delivered to each patient. The potential clinical effect of dosimetric differences in the planning target volume (PTV) and in organs at risk (OAR) were evaluated by comparing TCP and NTCP values. The Wilcoxon Signed Rank test was used for statistical comparison of all results obtained from the use of the two algorithms.

Results

The poorer coverage of the PTV, with higher prescribed dose, was reflected in the TCP, which was significantly lower when AXB was used, the median value was 81.55% (range: 74.90, 88.60%) and 84.10% (range: 77.70, 89.90%) for AAA (p < 0.001). OAR mean dose was lower in the AXB recalculated plan than the AAA plan and the difference was statistically significant for all the structures. The NTCP for developing mandible necrosis showed the largest median percentage difference between AAA and AXB (56.6%), the NTCP of risk for larynx edema of Grade ≥ 2 followed with 12.2%.

Conclusions

Differences in dose distribution of NPC treatment plans recalculated with AXB are of clinical significance in those situations where the PTV and OAR involve air or bone, media in which AXB has been shown to more accurately represent the true dose distribution. The availability of AXB algorithm could improve patient dose estimation, increasing the data consistency of clinical trials.

An Integrated Framework Based on Full Monte Carlo Simulations for Double-Scattering Proton Therapy

PURPOSE

We developed an integrated framework that employs a full Monte Carlo (MC) model for treatment-plan simulations of a passive double-scattering proton system.

MATERIALS AND METHODS

We have previously validated a virtual machine source model for full MC proton-dose calculations by comparing the percentage of depth-dose curves, spread-out Bragg peaks, and lateral profiles against measured commissioning data. This study further expanded our previous work by developing an integrate framework that facilitates its clinical use. Specifically, we have (1) constructed patient-specific applicator and compensator numerically from the plan data and incorporated them into the beamline, (2) created the patient anatomy from the computed tomography image and established the transformation between patient and machine coordinate systems, and (3) developed a graphical user interface to ease the whole process from importing the treatment plan in the Digital Imaging and Communications in Medicine format to parallelization of the MC calculations. End-to-end tests were performed to validate the functionality, and 3 clinical cases were used to demonstrate clinical utility of the framework.

RESULTS

The end-to-end tests demonstrated that the framework functioned correctly for all tested functionality. Comparisons between the treatment planning system calculations and MC results in 3 clinical cases revealed large dose difference up to 17%, especially in the beam penumbra and near the end of beam range. The discrepancy likely originates from a variety of sources, such as the dose algorithms, modeling of the beamline, and the dose metric. The agreement for other regions was acceptable.

CONCLUSION

An integrated framework was developed for full MC simulations of double-scattering proton therapy. It can be a valuable tool for dose verification and plan evaluation.

Commissioning and performance evaluation of RadCalc for the Elekta unity MRI-linac

Recent availability of MRI‐guided linear accelerators has introduced a number of clinical challenges, particularly in the context of online plan adaptation. Paramount among these is verification of plan quality prior to patient treatment. Currently, there are no commercial products available for monitor unit verification that fully support the newly FDA cleared Elekta Unity 1.5 T MRI‐linac. In this work, we investigate the accuracy and precision of RadCalc for this purpose, which is a software package that uses a Clarkson integration algorithm for point dose calculation. To this end, 18 IMRT patient plans (186 individual beams) were created and used for RadCalc point dose calculations. In comparison with the primary treatment planning system (Monaco), mean point dose deviations of 0.0 ± 1.0% (n = 18) and 1.7 ± 12.4% (n = 186) were obtained on a per‐plan and per‐beam basis, respectively. The dose plane comparison functionality within RadCalc was found to be highly inaccurate, however, modest improvements could be made by artificially shifting jaws and multi leaf collimator positions to account for the dosimetric shift due to the magnetic field (67.3% vs 96.5% mean 5%/5 mm gamma pass rate).

Dose to water versus dose to medium from cavity theory applied to small animal irradiation with kilovolt x-rays

Dose reporting is a matter of concern in the preclinical field as the different dose descriptors dose-to-water-in-medium [Formula: see text] and dose-to-medium-in-medium [Formula: see text] coexist. For kV photons differences between both quantities are expected to be amplified due to photon energy absorption coefficients differences for different media, and could represent a limiting factor for accurate translation of pre-clinical research into clinical trials. The main goal of this study was to analyse the relationship between [Formula: see text] and [Formula: see text] for kV irradiation of small animals, using different flavours of the intermediate cavity theory (ICT). Irradiations of mathematical phantoms and a mouse CT scan, both with different voxel sizes and materials, were investigated. A modified version of the Monte Carlo code DOSXYZnrc was used to derive [Formula: see text] and convert to [Formula: see text] using ICT. Local photon spectra were generated in different regions of the mouse. Depending on energy and cavity size, which we equate to the voxel size, [Formula: see text] ranged from 0.68 to 4.37 times [Formula: see text]. Higher kV energy combined with very small cavity sizes yielded decreased [Formula: see text] in comparison to [Formula: see text]; this behaviour was reversed for larger cavities combined with lower kV energies. Hence, the impact of the cavity dimensions on estimated [Formula: see text] is significant on pre-clinical kV beams. [Formula: see text] and [Formula: see text] in the ex vivo male mouse were found to differ by  -29% to 286%. Caution is advised when using the ICT due to a lack of consensus on weighting factor (d-parameter) deriving methods; for the same irradiation conditions, different d-values affected [Formula: see text] up to 20%. Pre-clinically, such divergence between dose descriptors could enable biological damage. The abiding debate over which quantity to favour is foreseen to linger while it is unclear which quantity correlates better with the biological effects of ionizing irradiation: preclinical radiotherapy might represent an ideal platform for measurement-based studies to settle this fundamental question. Finally, dose distribution comparisons require caution and should use the same reporting quantity.

Local control rates in stereotactic body radiotherapy (SBRT) of lung metastases associated with the biologically effective dose

Aim

To evaluate dose differences in lung metastases treated with stereotactic body radiotherapy (SBRT), and the correlation with local control, regarding the dose algorithm, target volume and tissue density.

Background

Several studies showed excellent local control rates in SBRT for lung metastases, with different fractionation schemes depending on the tumour location or size. These results depend on the dose distributions received by the lesions in terms of the tissue heterogeneity corrections performed by the dose algorithms.

Materials and methods

Forty-seven lung metastases treated with SBRT, using intrafraction control and respiratory gating with internal fiducial markers as surrogates (ExacTrac, BrainLAB AG), were calculated using Pencil Beam (PB) and Monte Carlo (MC) (iPlan, BrainLAB AG).Dose differences between both algorithms were obtained for the dose received by 99% (D _99%) and 50% (D _50%) of the planning treatment volume (PTV). The biologically effective dose delivered to 99% (BED_99%) and 50% (BED_50%) of the PTV were estimated from the MC results. Local control was evaluated after 24 months of median follow-up (range: 3-52 months).

Results

The greatest variations (40.0% in ΔD _99% and 38.4% in ΔD _50%) were found for the lower volume and density cases. The BED_99% and BED_50% were strongly correlated with observed local control rates: 100% and 61.5% for BED_99% > 85 Gy and <85 Gy (p < 0.0001), respectively, and 100% and 58.3% for BED_50% > 100 Gy and <100 Gy (p < 0.0001), respectively.

Conclusions

Lung metastases treated with SBRT, with delivered BED_99% > 85 Gy and BED_50% > 100 Gy, present better local control rates than those treated with lower BED values (p = 0.001).

Comparative study between Acuros XB algorithm and Anisotropic Analytical Algorithm in the case of heterogeneity for the treatment of lung cancer

The aim of this study was to investigate the impact of heterogeneity on the dose calculation for two algorithms implemented in the TPS “Analytical Anisotropic Algorithm (AAA) and Acuros XB” and validated the use of Acuros XB algorithm in clinical routine. First, we compare the dose calculated by these algorithms and the dose measured at the given point P, which is found after heterogeneity insert. Second, we extend our work on clinical cases that present a complex heterogeneity. By evaluating the impact of the choice of the algorithm on the dose coverage of the tumor, and the dose received by the organs at risk for 20 patients affected by lung cancer.

The result of our phantom study showed a good agreement with several studies that showed the superiority of the Acuros XB over the AAA in predicting dose when it concerns heterogeneous media. The treatment plans for 20 lung cancers were calculated by two algorithms AAA and Acuros XB, the results show a statistical significant difference between algorithms for Homogeneity Index and the maximum dose of planning target volume (HI: 0.11±0.01 vs 0.05±0.01 p = 0.04; Dmax: 69.30±3.12 vs 68.51±2.64 p = 0.02). Instead, no statistically significant difference was observed for conformity index CI and mean dose (CI: 0.98±0.18 vs 0.99±0.14 p = 0.33; Dmean: 66.3±0.65 vs 66.10 ±0.61 p = 0.54). For organs at risk, the maximum dose for spinal cord, mean dose and D37 % of lung minus GTV (dose receiving 37% of lung volume) were found to be lower for AAA plans than Acuros XB and the differences were statistically significant (p<0.05). For the heart D33% and D67% were found to be higher for AAA plans than Acuros XB and the differences were statistically significant (p<0.05), but No difference was observed for D100% of the heart.

The use of the AXB algorithm is suitable in the case of presence of heterogeneity, because it allows to have a better accuracy close to the Monte Carlo calculation.

Clinical evaluation for the difference of absorbed doses calculated to medium and calculated to water by Monte Carlo method

BACKGROUND

To evaluate the difference of absorbed doses calculated to medium and to water by a Monte Carlo (MC) algorithm based treatment planning system (TPS), and to assess the potential clinical impact to dose prescription.

METHODS

Thirty patients, 10 nasopharyngeal cancer (NPC), 10 lung cancer and 10 bone metastases cases, were selected for this study. For each case, the treatment plan was generated using a commercial MC based TPS and dose was calculated to medium (D_m). The plan was recalculated for dose to water (D_w) using the same Monitor Units (MU) and control points. The differences between D_m and D_w were qualitatively evaluated by dose-volume parameters and by the plan subtraction method. All plans were measured using the MapCheck2, and gamma passing rates were calculated.

RESULTS

For NPC and Lung cases, the mean differences between D_w and D_m for the targets were less than 2% and the maximum difference was 3.9%. The maximum difference of D_2% for the organs at risk (OARs) was 6.7%. The maximum differences between D_w and D_m were as high as 10% in certain high density regions. For bone metastases cases, the mean differences between D_w and D_m for the targets were more than 2.2% and the maximum difference was 7.1%. The differences between D_w and D_m for the OARs were basically negligible. At 3%&3 mm criterion, the gamma passing rate of D_w plan and D_m plan were close (> 94%).

CONCLUSION

The differences between D_w and D_m has little clinical impact for most clinical cases. In bony structures the differences may become clinically significant if the target/OAR is receiving doses close to its tolerance limit which can potentially influence the selection or rejection of a particular plan.

Monte Carlo simulations in radiotherapy dosimetry

BACKGROUND

The use of the Monte Carlo (MC) method in radiotherapy dosimetry has increased almost exponentially in the last decades. Its widespread use in the field has converted this computer simulation technique in a common tool for reference and treatment planning dosimetry calculations.

METHODS

This work reviews the different MC calculations made on dosimetric quantities, like stopping-power ratios and perturbation correction factors required for reference ionization chamber dosimetry, as well as the fully realistic MC simulations currently available on clinical accelerators, detectors and patient treatment planning.

CONCLUSIONS

Issues are raised that include the necessity for consistency in the data throughout the entire dosimetry chain in reference dosimetry, and how Bragg-Gray theory breaks down for small photon fields. Both aspects are less critical for MC treatment planning applications, but there are important constraints like tissue characterization and its patient-to-patient variability, which together with the conversion between dose-to-water and dose-to-tissue, are analysed in detail. Although these constraints are common to all methods and algorithms used in different types of treatment planning systems, they make uncertainties involved in MC treatment planning to still remain "uncertain".

Verification of Acuros XB dose algorithm using 3D printed low-density phantoms for clinical photon beams

The transport-based dose calculation algorithm Acuros XB (AXB) has been shown to accurately account for heterogeneities primarily through comparisons with Monte Carlo simulations. This study aims to provide additional experimental verification of AXB for clinically relevant flattened and unflattened beam energies in low density phantoms of the same material. Polystyrene slabs were created using a bench-top 3D printer. Six slabs were printed at varying densities from 0.23 to 0.68 g/cm3 , corresponding to different density humanoid tissues. The slabs were used to form different single and multilayer geometries. Dose was calculated with Eclipse™ AXB 11.0.31 for 6MV, 15MV flattened and 6FFF (flattening filter free) energies for field sizes of 2 × 2 and 5 × 5 cm2 . EBT3 film was inserted into the phantoms, which were irradiated. Absolute dose profiles and 2D Gamma analyses were performed for 96 dose planes. For all single slab configurations and energies, absolute dose differences between the AXB calculation and film measurements remained <3% for both fields in the high-dose region, however, larger disagreement was seen within the penumbra. For the multilayered phantom, percentage depth dose with AXB was within 5% of discrete film measurements. The Gamma index at 2%/2 mm averaged 98% in all combinations of fields, phantoms and photon energies. The transport-based dose algorithm AXB is in good agreement with the experimental measurements for small field sizes using 6MV, 6FFF and 15MV beams adjacent to various low-density heterogeneous media. This work provides preliminary experimental grounds to support the use of AXB for heterogeneous dose calculation purposes.

Dosimetrical and radiobiological approach to manage the dosimetric shift in the transition of dose calculation algorithm in radiation oncology: how to improve high quality treatment and avoid unexpected outcomes?

BACKGROUND

For a given prescribed dose of radiotherapy, with the successive generations of dose calculation algorithms, more monitor units (MUs) are generally needed. This is due to the implementation of successive improvements in dose calculation: better heterogeneity correction and more accurate estimation of secondary electron transport contribution. More recently, there is the possibility to report the dose-to-medium, physically more accurate compared to the dose-to-water as the reference one. This last point is a recent concern and the main focus of this study.

METHODS

In this paper, we propose steps for a general analysis procedure to estimate the dosimetric alterations, and the potential clinical changes, between a reference algorithm and a new one. This includes dosimetric parameters, gamma index, radiobiology indices based on equivalent uniform dose concept and statistics with bootstrap simulation. Finally, we provide a general recommendation on the clinical use of new algorithms regarding the dose prescription or dose limits to the organs at risks.

RESULTS

The dosimetrical and radiobiological data showed a significant effect, which might exceed 5-10%, of the calculation method on the dose the distribution and clinical outcomes for lung cancer patients. Wilcoxon signed rank paired comparisons indicated that the delivered dose in MUs was significantly increased (> 2%) using more advanced dose calculation methods as compared to the reference one.

CONCLUSION

This paper illustrates and explains the use of dosimetrical, radiobiologcal and statistical tests for dosimetric comparisons in radiotherapy. The change of dose calculation algorithm may induce a dosimetric shift, which has to be evaluated by the physicists and the oncologists. This includes the impact on tumor control and on the risk of toxicity based on normal tissue dose constraints. In fact, the alteration in dose distribution makes it hard to keep exactly the same tumor control probability along with the same normal tissue complication probability.

An automated Monte Carlo QC system for volumetric modulated arc therapy: Possibilities and challenges

PURPOSE

To develop and implement an automated Monte Carlo (MC) system for patient specific VMAT quality control in a patient geometry that generates treatment planning system (TPS) compliant DICOM objects and includes a module for 3D analysis of dose deviations. Also, the aims were to recommend diagnose specific tolerance criteria and an evaluation procedure.

METHODS

The EGSnrc code package formed the basis for development of the MC system. The workflow consists of a number of modules connected to a TPS by means of manual DICOM exports and imports which were executed sequentially without user interaction. DVH comparison was performed in the TPS. In addition, MC- and TPS dose distributions were analysed by applying the normalized dose difference (NDD) formalism. NDD failure maps and a pass rate for a certain threshold were obtained. 170 clinical plans (prostate, thorax, head-and-neck and gynecological) were selected for analysis.

RESULTS

Agreement within 1.5% was found between clinical- and MC data for the mean dose to the target volumes and within 3% for parameters more sensitive to the shape of the DVH e.g. D_98% PTV. Regarding the NDD analysis, tolerance criteria 2%/3 mm were established for prostate plans and 3%/3 mm for the rest of the cases.

CONCLUSIONS

An automated MC system was developed and implemented. Evaluation procedure is recommended with NDD-analysis as a first step. For pass rate < 95%, the evaluation continues with comparison of DVH parameters. For deviations larger than 2%, a visual inspection of the clinical- and MC dose distributions is performed.

On the conversion of dose to bone to dose to water in radiotherapy treatment planning systems

Background and purpose

Conversion factors between dose to medium (D_m,m) and dose to water (D_w,w) provided by treatment planning systems that model the patient as water with variable electron density are currently based on stopping power ratios. In the current paper it will be illustrated that this conversion method is not correct.

Materials and methods

Monte Carlo calculations were performed in a phantom consisting of a 2 cm bone layer surrounded by water. D_w,w was obtained by modelling the bone layer as water with the electron density of bone. Conversion factors between D_w,w and D_m,m were obtained and compared to stopping power ratios and ratios of mass-energy absorption coefficients in regions of electronic equilibrium and interfaces. Calculations were performed for 6 MV and 20 MV photon beams.

Results

In the region of electronic equilibrium the stopping power ratio of water to bone (1.11) largely overestimates the conversion obtained using the Monte Carlo calculations (1.06). In that region the MC dose conversion corresponds to the ratio of mass energy absorption coefficients. Near the water to bone interface, the MC ratio cannot be determined from stopping powers or mass energy absorption coefficients.

Conclusion

Stopping power ratios cannot be used for conversion from D_m,m to D_w,w provided by treatment planning systems that model the patient as water with variable electron density, either in regions of electronic equilibrium or near interfaces. In regions of electronic equilibrium mass energy absorption coefficient ratios should be used. Conversions at interfaces require detailed MC calculations.

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.

Comparison of calculated dose distributions reported as dose-to-water and dose-to-medium for intensity-modulated radiotherapy of nasopharyngeal cancer patients

Advanced dose calculation algorithms for radiation therapy treatment planning can report external beam photon dose 2-sided, in terms of dose-to-medium (D_m) and dose-to-water (D_w). The purpose of our study was to determinate the effect of D_w and D_m reporting modes built in Elekta Monaco treatment planning system on intensity-modulated radiotherapy dose distributions for patients with nasopharyngeal cancer. For 13 patients involved in this retrospective study, 2 plans were created: 1 using D_w and another according to D_m reporting mode. Treatment plans were normalized such that 100% planning target volume should be covered by 95% of prescribed dose. Dose-volume constraints were assigned according to international standards. The comparison between dose distributions was performed evaluating quantities important for respective volumes of interest. For target volumes, heterogeneity index and conformity index methodology were used along with the maximum dose concept. Also, for the comparisons over particular organ at risk, maximum dose or mean dose as well as dose-volume concepts were used. For all target volumes and majority of organs at risk, the differences between 2 reporting modes are statistically insignificant, but this is not the case for bony structured organs at risks: mandible and cochlea. It was observed that D_w is higher than D_m with mean difference of 9.91% (p = 0.000009) of the mandible volume covered with 70 Gy. The same trend was observed for left and right cochlea with difference in mean dose of 8.74% (p = 0.037) and 6.87% (p = 0.029), respectively. The comparative analysis of dosimetric parameters in this study shows that the selection of reporting modes in Monaco treatment planning system can produce dose differences up to 15% in high-density volumes such as mandible and cochlea, which might have clinical consequences.

Application of dual-energy CT to suppression of metal artefact caused by pedicle screw fixation in radiotherapy: a feasibility study using original phantom

The objective of the present study was the determination of the potential dosimetric benefits of using metal-artefact-suppressed dual-energy computed tomography (DECT) images for cases involving pedicle screw implants in spinal sites. A heterogeneous spinal phantom was designed for the investigation of the dosimetric effect of the pedicle-screw-related artefacts. The dosimetric comparisons were first performed using a conventional two-directional opposed (AP-PA) plan, and then a volumetric modulated arc therapy (VMAT) plan, which are both used for the treatment of spinal metastases in our institution. The results of Acuros^® XB dose-to-medium (Dm) and dose-to-water (Dw) calculations using different imaging options were compared with experimental measurements including the chamber and film dosimetries in the spinal phantom. A dual-energy composition image with a weight factor of  -0.2 and a dual-energy monochromatic image (DEMI) with an energy level of 180 keV were found to have superior abilities for artefact suppression. The Dm calculations revealed greater dosimetric effects of the pedicle screw-related artefacts compared to the Dw calculations. The results of conventional single-energy computed tomography showed that, although the pedicle screws were made from low-Z titanium alloy, the metal artefacts still have dosimetric effects, namely, an average (maximum) Dm error of 4.4% (5.6%) inside the spinal cord for a complex VMAT treatment plan. Our findings indicate that metal-artefact suppression using the proposed DECT (DEMI) approach is promising for improving the dosimetric accuracy near the implants and inside the spinal cord (average (maximum) Dm error of 1.1% (2.0%)).