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

Superficial Dosimetry Study of the Frequency of Bolus Using in Volumetric Modulated Arc Therapy after Modified Radical Mastectomy

OBJECTIVE

To investigate the effect of various frequencies of bolus use on the superficial dose of volumetric modulated arc therapy after modified radical mastectomy for breast cancer.

METHODS

Based on the computed tomography images of a female anthropomorphic breast phantom, a 0.5 cm silicone-based 3D-printed bolus was created. Nine points evenly distributed on the breast skin were selected for assessing the skin dose, and a volume of subcutaneous lymphatic drainage of the breast (noted as ROI2-3) was delineated for assessing the chest wall dose. The treatment plans with and without bolus (plan_wb and plan_nb) were separately designed using the prescription of 50 Gy in 25 fractions following the standard dose constraints of the adjacent organ at risk. To characterize the accuracy of treatment planning system (TPS) dose calculations, the doses of the nine points were measured five times by thermoluminescence dosimeters (TLDs) and then were compared with the TPS calculated dose.

RESULTS

Compared with Plan_nb (144.46 ± 10.32 cGy), the breast skin dose for plan_wb (208.75 ± 4.55 cGy) was significantly increased (t = -18.56, P < 0.001). The deviation of skin dose was smaller for Plan_wb, and the uniformity was significantly improved. The calculated value of TPS was in good agreement with the measured value of TLD, and the maximum deviation was within 5%. Skin and ROI2-3 doses were significantly increased with increasing frequencies of bolus applications. The mean dose of the breast skin and ROI2-3 for 15 and 23 times bolus applications were 45.33 Gy, 50.88 Gy and 50.36 Gy, 52.39 Gy, respectively.

CONCLUSION

3D printing bolus can improve the radiation dose and the accuracy of the planned dose. Setting Plan_wb to 15 times for T1-3N+ breast cancer patients and 23 times for T4N+ breast cancer patients can meet the clinical need. Quantitative analysis of the bolus application frequency for different tumor stages can provide a reference for clinical practice.

Sensitivity and specificity of secondary dose calculation for head and neck treatment plans

PURPOSE

Secondary dose calculation (SDC) with an independent algorithm is one option to perform plan-specific quality assurance (QA). While measurement-based QA can potentially detect errors in plan delivery, the dose values are typically only compared to calculations on homogeneous phantom geometries instead of patient CT data. We analyzed the sensitivity and specificity of an SDC software by purposely introducing different errors and determined thresholds for optimal decisions.

METHODS

Thirty head and neck VMAT plans and 30 modifications of those plans, including errors related to density and beam modelling, were recalculated using RadCalc with a Monte Carlo algorithm. Decision thresholds were obtained by receiver operating characteristics (ROC) analysis. For comparison, measurement-based QA using the ArcCHECK phantom was carried out and evaluated in the same way.

RESULTS

Despite optimized decision thresholds, none of the systems was able to reliably detect all errors. ArcCHECK analysis using a 2%/2 mm criterion with a threshold of 91.1% had an area under the curve (AUC) of 0.80. Evaluating differences in recalculated and planned DVH parameter of the target structures in RadCalc with a 2% threshold an AUC of 0.86 was achieved. Out-of-field deviations could be attributed to weaknesses in the beam model.

CONCLUSIONS

Secondary dose calculation with RadCalc is an alternative to established measurement-based phantom QA. Different tools catch different errors; therefore, a combination of approaches should be preferred.

A structure-based gamma evaluation method for identifying clinically relevant dose differences in organs at risk

Gamma evaluation is currently the most widely used dose comparison method for patient specific quality assurance (PSQA). However, existing methods for normalising the dose difference, using either the dose at the global maximum dose point or at each local point, can respectively lead to under- and over-sensitivity to dose differences in organ-at-risk structures. This may be of concern for plan evaluation from clinical perspectives. This study has explored and proposed a new method called structural gamma, which takes structural dose tolerances into consideration while performing gamma analysis for PSQA. As a demonstration of the structural gamma method, a total of 78 retrospective plans on four treatment sites were re-calculated on an in-house Monte Carlo system and compared with doses calculated from the treatment planning system. Structural gamma evaluations were performed using both QUANTEC dose tolerances and radiation oncologist specified dose tolerances, then compared with conventional global and local gamma evaluations. Results demonstrated that structural gamma evaluation is especially sensitive to errors in structures with restrictive dose constraints. The structural gamma map provides both geometric and dosimetric information on PSQA results, allowing straightforward clinical interpretation. The proposed structure-based gamma method accounts for dose tolerances for specific anatomical structures. This method can provide a clinically useful method to assess and communicate PSQA results, offering radiation oncologists a more intuitive way of examining agreement in surrounding critical normal structures.

Beam modeling and commissioning for Monte Carlo photon beam on an Elekta Versa HD LINAC

PURPOSE

This study aims at analyzing beam data commissioning along with verifying Monte Carlo-based treatment planning system on the basis of the manufacturer guidelines for Elekta Versa HD Linear Accelerator. Moreover, evaluating the beam match process in terms of quality index, photon profile and multi leaf collimator (MLC) offset is aimed as well.

MATERIALS AND METHODS

The process of collecting beam data for Monaco 5.51 Treatment Planning System commissioning was done based on the instructions provided by the manufacturer as well as AAPM TG-106. Monte Carlo analysis was done for output factors in water, percent depth dose and beam profiles. A set of eight static and intensity modulated radiation therapy fields were used to verify the MLC parameters.

RESULTS

The difference between the measured and modeled penetrative quality (D₁₀) was achieved to be 0.54%. The output factors for 6 MV photon energy were measured and the difference between the measured and Monte Carlo output results was smaller than 1% for all the fields. The average percentage of passing the gamma criteria for commissioning test fields was (95+-4)%, however, the minimum agreement was 87.5% belonging to "7SEGA". Additionally, the agreement between both LINAC is 96%, however, the second LINAC reveals a positive offset in the point of approximately -4 cm on the x-axis at the crossplane.

CONCLUSION

Test commissioning was successfully verified using a homogeneous phantom for point dose measurements, post modelling MLC parameters and patient QA plans. All plan parameters pass the gamma criteria. 6 MV photon beam was successfully commissioned for Elekta VersaHD LINAC and is ready for clinical implementation.

Recommendations for simulating and measuring with biofabricated lung equivalent materials based on atomic composition analysis

Monte Carlo simulations of lung equivalent materials often involve the density being artificially lowered rather than a true lung tissue (or equivalent plastic) and air composition being simulated. This study used atomic composition analysis to test the suitability of this method. Atomic composition analysis was also used to test the suitability of 3D printing PLA or ABS with air to simulate lung tissue. It was found that there was minimal atomic composition difference when using an artificially lowered density, with a 0.8 % difference in Nitrogen the largest observed. Therefore, excluding infill pattern effects, lowering the density of the lung tissue (or plastic) in simulations should be sufficiently accurate to simulate an inhaled lung, without the need to explicitly include the air component. The average electron density of 3D printed PLA and air, and ABS and air were just 0.3 % and 1.3 % different to inhaled lung, confirming their adequacy for MV photon dosimetry. However large average atomic number differences (5.6 % and 20.4 % respectively) mean that they are unlikely to be suitable for kV photon dosimetry.

Monte Carlo calculations of radiotherapy dose in "homogeneous" anatomy

Given the substantial literature on the use of Monte Carlo (MC) simulations to verify treatment planning system (TPS) calculations of radiotherapy dose in heterogeneous regions, such as head and neck and lung, this study investigated the potential value of running MC simulations of radiotherapy treatments of nominally homogeneous pelvic anatomy. A pre-existing in-house MC job submission and analysis system, built around BEAMnrc and DOSXYZnrc, was used to evaluate the dosimetric accuracy of a sample of 12 pelvic volumetric arc therapy (VMAT) treatments, planned using the Varian Eclipse TPS, where dose was calculated with both the Analytical Anisotropic Algorithm (AAA) and the Acuros (AXB) algorithm. In-house TADA (Treatment And Dose Assessor) software was used to evaluate treatment plan complexity, in terms of the small aperture score (SAS), modulation index (MI) and a novel exposed leaf score (ELS/ELA). Results showed that the TPS generally achieved closer agreement with the MC dose distribution when treatments were planned for smaller (single-organ) targets rather than larger targets that included nodes or metastases. Analysis of these MC results with reference to the complexity metrics indicated that while AXB was useful for reducing dosimetric uncertainties associated with density heterogeneity, the residual TPS dose calculation uncertainties resulted from treatment plan complexity and TPS model simplicity. The results of this study demonstrate the value of using MC methods to recalculate and check the dose calculations provided by commercial radiotherapy TPSs, even when the treated anatomy is assumed to be comparatively homogeneous, such as in the pelvic region.

Exploring the gamma surface: A new method for visualising modulated radiotherapy quality assurance results

PURPOSE

This work presents a novel method of visualising the results of patient-specific quality assurance (QA) for modulated radiotherapy treatment plans, using a three-dimensional distribution of gamma pass rates, referred to as the "gamma surface". The method was developed to aid in comparing borderline and failing QA plans, and to better compare patient-specific QA results between departments.

METHODS

Gamma surface plots were created for a representative sample of situations encountered during patient-specific QA. To produce a gamma surface plot, for each QA result, gamma pass rates were plotted as a heat map, with dose difference on one axis and distance-to-agreement on the other. This involved the calculation of 100 × 100 gamma pass rates over a dose difference and distance-to-agreement grid. As examples, five 220 × 680 arrays of dose points from radiotherapy treatment plans were compared against measurement data consisting of 21 × 66 arrays of dose points spaced 10 mm apart.

RESULTS

The gamma surface plots facilitated the rapid evaluation of criteria combinations for each plan, clearly highlighting the difference between plans that are modelled and delivered well, and those that are not. Large scale features were also evident in each surface, hinting at potential over-modulation, systematic dose errors, and small or large scale areas of disagreement in the distributions.

CONCLUSIONS

Gamma surface plots are a useful tool for investigating QA failures and borderline results, and have the capacity to grant insights into treatment plan QA performance that may otherwise be missed.