Impact of dose calculation algorithm on radiation therapy

Advanced External Beam Stereotactic Radiotherapy for Skull Base Reirradiation

BACKGROUND/OBJECTIVES

Stereotactic body radiation therapy (SBRT) for skull base reirradiation is particularly challenging, as patients have already received substantial radiation doses to the region, and nearby normal organs may have approached their tolerance limit from prior treatments. In this study, we reviewed the characteristics and capabilities of four advanced external beam radiation delivery systems and four modern treatment planning systems and evaluated the treatment plan quality of each technique using skull base reirradiation patient cases.

METHODS

SBRT plans were generated for sixteen skull base reirradiation patients using four modalities: the GK plan for the Elekta Leksell Gamma Knife Perfexion/ICON, the CyberKnife (CK) plan for the Accuray CyberKnife, the intensity-modulated proton therapy (IMPT) plan for the Hitachi ProBeat-FR proton therapy machine, and the volumetric-modulated arc therapy (VMAT) plan for the Varian TrueBeam STx. These plans were evaluated and compared using two novel gradient indices in addition to traditional dosimetry metrics for targets and organs at risk (OARs). The steepest border gradient quantified the percent prescription dose fall-off per millimeter at the boundary between the target and adjacent critical structures. This gradient index highlighted the system's ability to spare nearby critical OARs. The volume gradient assessed the extent of dose spread outside the target toward the patient's body.

RESULTS

All plans achieved comparable target coverage and conformity, while IMPT and VMAT demonstrated significantly better uniformity. The GK plans exhibited the highest border gradient, up to 20.9%/mm, followed by small-spot-size IMPT plans and CK plans. Additionally, IMPT plans showed the benefit of reduced dose spread in low-dose regions and the lowest maximum and mean doses to the brainstem and carotid artery.

CONCLUSIONS

The advanced external beam radiotherapy modalities evaluated in this study are well-suited for SBRT in skull base reirradiation, which demands precise targeting of tumors with highly conformal doses and steep dose gradients to protect nearby normal structures.

Evaluation of Dose Calculation Algorithms Accuracy for ISOgray Treatment Planning System in Motorized Wedged Treatment Fields

Background

Different dose calculation methods vary in accuracy and speed. While most methods sacrifice precision for efficiency Monte Carlo (MC) simulation offers high accuracy but slower calculation. ISOgray treatment planning system (TPS) uses Clarkson, collapsed cone convolution (CCC), and fast Fourier transform (FFT) algorithms for dose distribution. This study's primary goal is to evaluate the dose calculation accuracy for ISOgray TPS algorithms in the presence of a wedge.

Methods

This study evaluates the dose calculation algorithms using the ISOgray TPS in the context of radiation therapy. The authors compare ISOgray TPS algorithms on an Elekta Compact LINAC through MC simulations. The study compares MC simulations for open and wedge fields with ISOgray algorithms by using gamma index analysis for validation.

Results

The percentage depth dose results for all open and wedge fields showed a more than 98% pass rate for points. However, there were differences in the dose profile gamma index results. Open fields passed the gamma index analysis in the in-plane direction, but not all points passed in the cross-plane direction. Wedge fields passed in the cross-plane direction, but not all in the in-plane direction, except for the Clarkson algorithms.

Conclusion

In all investigated algorithms, error increases in the penumbra areas, outside the field, and at cross-plane of open fields and in-plane direction of wedged fields. By increasing the wedge angle, the discrepancy between the TPS algorithms and MC simulations becomes more pronounced. This discrepancy is attributed to the increased presence of scattered photons and the variation in the delivered dose within the wedge field, consequently impacts the beam quality. While the CCC and FFT algorithms had better accuracy, the Clarkson algorithm, particularly at larger effective wedge angles, exhibited greater effectiveness than the two mentioned algorithms.

Technical recommendations for implementation of Volumetric Modulated Arc Therapy and Helical Tomotherapy Total Body Irradiation

As a component of myeloablative conditioning before allogeneic hematopoietic stem cell transplantation (HSCT), Total Body Irradiation (TBI) is employed in radiotherapy centers all over the world. In recent and coming years, many centers are changing their technical setup from a conventional TBI technique to multi-isocenter conformal arc therapy techniques such as Volumetric Modulated Arc Therapy (VMAT) or Helical Tomotherapy (HT). These techniques allow better homogeneity and control of the target prescription dose, and provide more freedom for individualized organ-at-risk sparing. The technical design of multi-isocenter/multi-plan conformal TBI is complex and should be developed carefully. A group of early adopters with conformal TBI experience using different treatment machines and treatment planning systems came together to develop technical recommendations and share experiences, in order to assist departments wishing to implement conformal TBI, and to provide ideas for standardization of practices.

Are offline ART decisions for NSCLC impacted by the type of dose calculation algorithm?

Introduction

Decisions for plan-adaptations may be impacted by a transitioning from one dose-calculation algorithm to another. This study examines the impact on dosimetric-triggered offline adaptation in LA-NSCLC in the context of a transition from superposition/convolution dose calculation algorithm (Type-B) to linear Boltzmann equation solver dose calculation algorithms (Type-C).

Materials & Methods

Two dosimetric-triggered offline adaptive treatment workflows are compared in a retrospective planning study on 30 LA-NSCLC patients. One workflow uses a Type-B dose calculation algorithm and the other uses Type-C. Treatment plans were re-calculated on the anatomy of a mid-treatment synthetic-CT utilizing the same algorithm utilized for pre-treatment planning. Assessment for plan-adaptation was evaluated through a decision model based on target coverage and OAR constraint violation. The impact of algorithm during treatment planning was controlled for by recalculating the Type-B plan with Type-C.

Results

In the Type-B approach, 13 patients required adaptation due to OAR-constraint violations, while 15 patients required adaptation in the Type-C approach. For 8 out of 30 cases, the decision to adapt was opposite in both approaches. None of the patients in our dataset encountered CTV-target underdosage that necessitated plan-adaptation. Upon recalculating the Type-B approach with the Type-C algorithm, it was shown that 10 of the original Type-B plans revealed clinically relevant dose reductions (≥3%) on the CTV in their original plans. This re-calculation identified 21 plans in total that required ART.

Discussion

In our study, nearly one-third of the cases would have a different decision for plan-adaption when utilizing Type-C instead of Type-B. There was no substantial increase in the total number of plan-adaptations for LA-NSCLC. However, Type-C is more sensitive to altered anatomy during treatment compared to Type-B. Recalculating Type-B plans with the Type-C algorithm revealed an increase from 13 to 21 cases triggering ART.

Dosimetric comparison between Intensity Modulated Radiation Therapy (IMRT) vs dual arc Volumetric Arc Therapy (VMAT) for nasopharyngeal cancer (NPC): Systematic review and meta-analysis

BACKGROUND AND PURPOSE

Previous non-randomised studies comparing dosimetric outcomes between advanced techniques such as IMRT and VMAT reported conflicting findings. We thus sought to perform a systematic review and meta-analysis to consolidate the findings of these studies.

MATERIALS AND METHODS

We searched PUBMED and EMBASE for eligible studies from their time of inception to 10 March 2022. A random effects model was used to estimate the pooled mean differences (MDs) and their 95% confidence intervals(CIs) for target volume coverage, organ-at-risk(OAR) doses, monitor units(MUs) and treatment delivery times. We also performed a subgroup analysis to evaluate if different treatment planning systems (TPS) (Eclipse, Monaco and Pinnacle) used affected the pooled mean differences.

RESULTS

A total of 17 studies (383 patients) were eligible to be included. The pooled results showed that dual arc VMAT reduced D2% of PTV (MD=0.71Gy,95%CI=0.14-1.27,P=0.01), mean left cochlea dose (MD=2.6Gy,95%CI=0.03-5.16,P=0.05), mean right cochlea dose (MD=3.4Gy,95%CI=0.7-6.1,P=0.01), MUs (MD=554.9,95%CI=245.8-863.9,P=0.0004), treatment delivery times (MD=6.7mins,95%CI=4.5-8.9,P<0.0001) and integral dose (MD=0.97Gy,95%CI=0.28-1.67,P=0.006). None of the other indices were significantly better for the IMRT plans. The subgroup analysis showed that the integral dose was significantly lower only for Eclipse (MD=0.88Gy, 95%CI=0.14-1.63, P=0.02). The total MUs was significantly lower only for Eclipse (MD=1035.2, 95%CI=624.6-1445.9, P<0.0001) and Pinnacle (MD=293, 95%CI=15.6-570.5, P=0.04). Similarly, delivery time was also significantly lower only for Eclipse (MD=6.1mins, 95%CI=5.7-6.5, P<0.0001) and Pinnacle (MD=4.9mins, 95%CI=2.6-7.2, P<0.0001). The subgroup analysis however showed that target coverage was superior for the IMRT plans for both Pinnacle (MD=0.48Gy, 95%CI=0.31-0.66, P<0.0001) and Monaco (MD=0.12Gy, 95%CI=0.07-0.17, P<0.0001).

CONCLUSION

Dual-arc VMAT plans improved OAR doses, MUs and treatment times as compared to IMRT plans. The different TPS used may modify dosimetric outcomes.

An ultra-fast deep-learning-based dose engine for prostate VMAT via knowledge distillation framework with limited patient data

Objective. Deep-learning (DL)-based dose engines have been developed to alleviate the intrinsic compromise between the calculation accuracy and efficiency of the traditional dose calculation algorithms. However, current DL-based engines typically possess high computational complexity and require powerful computing devices. Therefore, to mitigate their computational burdens and broaden their applicability to a clinical setting where resource-limited devices are available, we proposed a compact dose engine via knowledge distillation (KD) framework that offers an ultra-fast calculation speed with high accuracy for prostate Volumetric Modulated Arc Therapy (VMAT).Approach. The KD framework contains two sub-models: a large pre-trained teacher and a small to-be-trained student. The student receives knowledge transferred from the teacher for better generalization. The trained student serves as the final engine for dose calculation. The model input is patient computed tomography and VMAT dose in water, and the output is DL-calculated patient dose. The ground-truth \dose was computed by the Monte Carlo module of the Monaco treatment planning system. Twenty and ten prostate cases were included for model training and assessment, respectively. The model's performance (teacher/student/student-only) was evaluated by Gamma analysis and inference efficiency.Main results. The dosimetric comparisons (input/DL-calculated/ground-truth doses) suggest that the proposed engine can effectively convert low-accuracy doses in water to high-accuracy patient doses. The Gamma passing rate (2%/2 mm, 10% threshold) between the DL-calculated and ground-truth doses was 98.64 ± 0.62% (teacher), 98.13 ± 0.76% (student), and 96.95 ± 1.02% (student-only). The inference time was 16 milliseconds (teacher) and 11 milliseconds (student/student-only) using a graphics processing unit device, while it was 936 milliseconds (teacher) and 374 milliseconds (student/student-only) using a central processing unit device.Significance. With the KD framework, a compact dose engine can achieve comparable accuracy to that of a larger one. Its compact size reduces the computational burdens and computing device requirements, and thus such an engine can be more clinically applicable.

Implementation and experimental evaluation of Mega-voltage fan-beam CT using a linear accelerator

BACKGROUND

Mega-voltage fan-beam Computed Tomography (MV-FBCT) holds potential in accurate determination of relative electron density (RED) and proton stopping power ratio (SPR) but is not widely available.

OBJECTIVE

To demonstrate the feasibility of MV-FBCT using a medical linear accelerator (LINAC) with a 2.5 MV imaging beam, an electronic portal imaging device (EPID) and multileaf collimators (MLCs).

METHODS

MLCs were used to collimate MV beam along z direction to enable a 1 cm width fan-beam. Projection data were acquired within one gantry rotation and preprocessed with in-house developed artifact correction algorithms before the reconstruction. MV-FBCT data were acquired at two dose levels: 30 and 60 monitor units (MUs). A Catphan 604 phantom was used to evaluate basic image quality. A head-sized CIRS phantom with three configurations of tissue-mimicking inserts was scanned and MV-FBCT Hounsfield unit (HU) to RED calibration was established for each insert configuration using linear regression. The determination coefficient ([Formula: see text]) was used to gauge the accuracy of HU-RED calibration. Results were compared with baseline single-energy kilo-voltage treatment planning CT (TP-CT) HU-RED calibration which represented the current standard clinical practice.

RESULTS

The in-house artifact correction algorithms effectively suppressed ring artifact, cupping artifact, and CT number bias in MV-FBCT. Compared to TP-CT, MV-FBCT was able to improve the prediction accuracy of the HU-RED calibration curve for all three configurations of insert materials, with [Formula: see text] > 0.9994 and [Formula: see text] < 0.9990 for MV-FBCT and TP-CT HU-RED calibration curves of soft-tissue inserts, respectively. The measured mean CT numbers of blood-iodine mixture inserts in TP-CT drastically deviated from the fitted values but not in MV-FBCT. Reducing the radiation level from 60 to 30 MU did not decrease the prediction accuracy of the MV-FBCT HU-RED calibration curve.

CONCLUSION

We demonstrated the feasibility of MV-FBCT and its potential in providing more accurate RED estimation.

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.

Investigation of tube voltage dependence on CT number and its effect on dose calculation algorithms using thorax phantom in Monaco treatment planning system for external beam radiation therapy

Introduction

The accuracy of dose calculation algorithms depends on the electron density and computed tomography (CT) number of medium scanned. Our study aimed to verify the impact of different CT scanning protocols on Hounsfield unit (HU) and effect on dose calculation algorithms.

Materials and Methods

CIRS thorax phantom with different density material plugs was scanned at varying tube voltages from CT scanner and HU values were measured in treatment planning system (TPS). Calibration curves of electron density at different tube voltages were plotted and used for dose calculation with different calculation algorithms at varying high energy megavoltage photon energies.

Results

Insignificant difference is obtained in electron density curves plotted at different tube voltages. The mean variation in HU values was found at different tube voltages for bone, lung, and water are 896.75 (standard deviation [SD] 122.88), -799.25 (SD 5.74), and -17.5 (SD 0.57), respectively. The estimated P values for change in HU values were 0.089, 0.258, and 0.121 for bone, lung, and water, respectively. Pencil beam (PB) convolution and collapsed cone algorithms show no significant dose difference, i.e., <1% variation and Monte Carlo (MC) shows maximum dose difference up to 1.4%.

Conclusion

Third-generation algorithms such as MC shows dependence on varying tube voltages in dose calculation. Calibration curves plotted at different kVp in TPS advised to be chosen wisely to avoid any dosimetric errors in different medium.

Impact of Treating Physician on Radiation Therapy Related Severe Toxicities in Men with Prostate Cancer

PURPOSE

The impact of treating physician on radiation therapy (RT) related toxicity is unclear. We carried out a secondary analysis of a randomized controlled study to determine whether the risk of RT-related late toxicities in patients with prostate cancer varies depending on the treating radiation oncologist.

METHODS AND MATERIALS

This is a secondary analysis of a phase 3 randomized controlled study in which patients with prostate cancer with Gleason score ≤7, clinical stage T1b-T3a, and prostate-specific antigen <30 ng/mL were randomized to receive androgen suppression for 6 months, starting either 4 months before or concurrently with definitive prostate radiation therapy. Incidence of late RT-related toxicity was estimated using Kaplan-Meier methods. We applied multivariable semiparametric shared frailty models with gamma distribution to determine the between-physician variation in the hazard of late RT-related grade ≥3 gastrointestinal, genitourinary, or overall toxicity. Patient level covariables included age, risk group, year of enrollment, and treatment regimen. Frailty variance, a measure of unexplained heterogeneity, was estimated with 95% confidence intervals (CIs). Statistical significance was suggested when the lower limit of the 95% CI for the frailty variance was >0. The Commenges-Andersen test was used for P value estimation.

RESULTS

Overall, 426 patients were treated by 9 radiation oncologists. On log-rank test, there was a significant difference in the cumulative incidence of overall grade ≥3 toxicities (P = .001) and grade ≥3 gastrointestinal toxicity (P = .01) among the physician-based clusters. The frailty variance for overall late grade ≥3 toxicity was 0.31 (95% CI, 0.02-1.39; P = .01). The frailty variance for the grade ≥3 gastrointestinal and genitourinary toxicity was 0.84 (95% CI, 0.00-4.20; P = .11) and 0.11 (95% CI, 0.00-1.13; P = .31), respectively.

CONCLUSIONS

In our study, the hazard of overall RT-related late grade ≥3 toxicity varied significantly depending on treating radiation oncologist. Further studies are required to explore the underlying processes that lead to such variations in clinical trials involving radiation therapy in prostate cancer.

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.

Dosimetric comparison of 3-dimensional conformal radiotherapy (3D-CRT) and volumetric-modulated arc therapy (VMAT) in locally advanced cancer cervix

Introduction:

Dosimetric advantages of volumetric-modulated arc therapy (VMAT) over three-dimensional conformal radiotherapy (3D-CRT) are not established in a head-on comparison of a uniform group of locally advanced carcinoma of the cervix (LACC). Therefore, we conducted a dosimetric comparison of these two techniques in LACC patients.

Materials and methods:

Computed tomography (CT) data of histologically proven de novo LACC, including Stage IIB–IIIB and earlier stages deemed inoperable, were included in this prospective observational dosimetric study. Planning was initially done by 3D-CRT technique (dose of 45–50·4 Gy @ 1·8–2 Gy/# was used in the actual treatment), followed by VMAT planning and appropriate dosimetric comparisons were done in 39 cases.

Results:

For planning target volume coverage, D95, D98 and D100 (p < 0·0001 for all parameters) and V95 and V100 (p = 0·002 and <0·0001, respectively) were significantly improved with VMAT. The conformity index (CI) was significantly better with VMAT (p = 0·03), while 3D-CRT had a significantly better homogeneity index (HI)(p = 0·003). Dose to the urinary bladder was significantly reduced with VMAT compared to 3D-CRT for V20–V50, except V10. The doses to the rectum and abdominal cavity were significantly reduced with VMAT compared to 3D-CRT plans for all parameters (V10–V50). The number of organs at risks (OARs) for which constraints were met was higher with VMAT plans than with 3D-CRT plans, with at least four out of the five OARs protected in 46·1 versus 5·1% and all constraints achieved in 15·4% versus none.

Conclusion:

We conclude that in dosimetric terms, VMAT is superior to 3D-CRT for LACC.

Boosting radiotherapy dose calculation accuracy with deep learning

In radiotherapy, a trade‐off exists between computational workload/speed and dose calculation accuracy. Calculation methods like pencil‐beam convolution can be much faster than Monte‐Carlo methods, but less accurate. The dose difference, mostly caused by inhomogeneities and electronic disequilibrium, is highly correlated with the dose distribution and the underlying anatomical tissue density. We hypothesize that a conversion scheme can be established to boost low‐accuracy doses to high‐accuracy, using intensity information obtained from computed tomography (CT) images. A deep learning‐driven framework was developed to test the hypothesis by converting between two commercially available dose calculation methods: Anisotropic analytic algorithm (AAA) and Acuros XB (AXB). A hierarchically dense U‐Net model was developed to boost the accuracy of AAA dose toward the AXB level. The network contained multiple layers of varying feature sizes to learn their dose differences, in relationship to CT, both locally and globally. Anisotropic analytic algorithm and AXB doses were calculated in pairs for 120 lung radiotherapy plans covering various treatment techniques, beam energies, tumor locations, and dose levels. For each case, the CT and the AAA dose were used as the input and the AXB dose as the “ground‐truth” output, to train and test the model. The mean squared errors (MSEs) and gamma passing rates (2 mm/2% & 1 mm/1%) were calculated between the boosted AAA doses and the “ground‐truth” AXB doses. The boosted AAA doses demonstrated substantially improved match to the “ground‐truth” AXB doses, with average (± s.d.) gamma passing rate (1 mm/1%) 97.6% (±2.4%) compared to 87.8% (±9.0%) of the original AAA doses. The corresponding average MSE was 0.11(±0.05) vs 0.31(±0.21). Deep learning is able to capture the differences between dose calculation algorithms to boost the low‐accuracy algorithms. By combining a less accurate dose calculation algorithm with a trained deep learning model, dose calculation can potentially achieve both high accuracy and efficiency.

Investigation of central electrode artefacts of ionisation chamber effect on dose calculation using advanced calculation algorithms AAA and Acuros XB

Aim:

To investigate the central electrode artefact effect of different ion chambers in the verification phantom using the dose calculation algorithms Analytical Anisotropic Algorithm (AAA) and Acuros XB.

Materials and methods:

The dosimetric study was conducted using an in-house fabricated polymethyl methacrylate head phantom. The treatment planning system (TPS)-calculated doses in the phantom with detectors were compared against the dummy detector fillets using AAA and Acuros XB algorithm. The planned and measured doses were compared for the study.

Results:

The mean percentage variation in volumetric-modulated arc therapy plans using Acuros XB and the measurement in the head phantom are statistically significant (p-value = 0.001) for FC65 and CC13 chambers. In small volume chambers (A14SL and CC01), the measured and TPS-calculated dose shows a good agreement.

Findings:

The study confirmed the CT set of the phantom with detectors (FC65 and CC13) give more artefacts/heterogeneity caused a significant variation in dose calculation using Acuros XB. Therefore, the study suggests a method of using phantom CT set with the dummy detector for mean dose calculation for the Acuros XB algorithm.

Mid-position treatment strategy for locally advanced lung cancer: a dosimetric study

OBJECTIVE

The internal target volume (ITV) strategy generates larger planning target volumes (PTVs) in locally advanced non-small cell lung cancer (LA-NSCLC) than the Mid-position (Mid-p) strategy. We investigated the benefit of the Mid-p strategy regarding PTV reduction and dose to the organs at risk (OARs).

METHODS

44 patients with LA-NSCLC were included in a randomized clinical study to compare ITV and Mid-p strategies. GTV were delineated by a physician on maximum intensity projection images and on Mid-p images from four-dimensional CTs. CTVs were obtained by adding 6 mm uniform margin for microscopic extension. CTV to PTV margins were calculated using the van Herk's recipe for setup and delineation errors. For the Mid-p strategy, the mean target motion amplitude was added as a random error. For both strategies, three-dimensional conformal plans delivering 60-66 Gy to PTV were performed. PTVs, dose-volume parameters for OARs (lung, esophagus, heart, spinal cord) were reported and compared.

RESULTS

With the Mid-p strategy, the median of volume reduction was 23.5 cm3 (p = 0.012) and 8.8 cm3 (p = 0.0083) for PTVT and PTVN respectively; the median mean lung dose reduction was 0.51 Gy (p = 0.0057). For 37.1% of the patients, delineation errors led to smaller PTV with the ITV strategy than with the Mid-p strategy.

CONCLUSION

PTV and mean lung dose were significantly reduced using the Mid-p strategy. Delineation uncertainty can unfavorably impact the advantage.

ADVANCES IN KNOWLEDGE

To the best of our knowledge, this is the first dosimetric comparison study between ITV and Mid-p strategies for LA-NSCLC.

Migration of treatment planning system using existing commissioned planning system

Introduction:

Commissioning of a new planning system involves extensive data acquisition which can be onerous involving significant clinic downtime. This could be circumvented by extracting data from existing treatment planning system (TPS) to speed up the process.

Material and methods:

In this study, commissioning beam data was obtained from a clinically commissioned TPS (Pinnacle™) using Matlab™ generated Pinnacle™ executable scripts to commission an independent 3D dose verification TPS (Eclipse™). Profiles and output factors for commissioning as required by Eclipse™ were computed on a 50 × 50 × 50 cm3 water phantom at a dose grid resolution of 2 mm3. Verification doses were computed and compared to clinical TPS dose profiles based on TG-106 guidelines. Standard patient plans from Pinnacle™ including intensity modulated radiation therapy and volumetric modulated arc therapy were re-computed on Eclipse™ TPS while maintaining the same monitor units. Computed dose was exported back to Pinnacle for comparison with the original plans. This methodology enabled us to alleviate all ambiguities that arise in such studies.

Results:

Profile analysis using in-house software showed that for all field sizes including small multi-leaf collimator-generated fields, >95% of infield and penumbra data points of Eclipse™ match Pinnacle™ generated and measured profiles with 2%/2 mm gamma criteria. Excellent agreement was observed in the penumbra regions, with >95% of the data points passing distance to agreement criteria for complex C-shaped and S-shaped profiles. Dose volume histograms and isodose lines of patient plans agreed well to within a 0·5% for target coverage.

Findings:

Migration of TPS is possible without compromising accuracy or enduring the cumbersome measurement of commissioning data. Economising time for commissioning such a verification system or for migration of TPS can add great QA value and minimise downtime.

Comparative photon and proton dosimetry for patients with mediastinal lymphoma in the era of Monte Carlo treatment planning and variable relative biological effectiveness

Background

Existing pencil beam analytical (PBA) algorithms for proton therapy treatment planning are not ideal for sites with heterogeneous tissue density and do not account for the spatial variations in proton relative biological effectiveness (vRBE). Using a commercially available Monte Carlo (MC) treatment planning system, we compared various dosimetric endpoints between proton PBA, proton MC, and photon treatment plans among patients with mediastinal lymphoma.

Methods

Eight mediastinal lymphoma patients with both free breathing (FB) and deep inspiration breath hold (DIBH) CT simulation scans were analyzed. The original PBA plans were re-calculated with MC. New proton plans that used MC for both optimization and dose calculation with equivalent CTV/ITV coverage were also created. A vRBE model, which uses a published model for DNA double strand break (DSB) induction, was applied on MC plans to study the potential impact of vRBE on cardiac doses. Comparative photon plans were generated on the DIBH scan.

Results

Re-calculation of FB PBA plans with MC demonstrated significant under coverage of the ITV V99 and V95. Target coverage was recovered by re-optimizing the PT plan with MC with minimal change to OAR doses. Compared to photons with DIBH, MC-optimized FB and DIBH proton plans had significantly lower dose to the mean lung, lung V5, breast tissue, and spinal cord for similar target coverage. Even with application of vRBE in the proton plans, the putative increase in RBE at the end of range did not decrease the dosimetric advantages of proton therapy in cardiac substructures.

Conclusions

MC should be used for PT treatment planning of mediastinal lymphoma to ensure adequate coverage of target volumes. Our preliminary data suggests that MC-optimized PT plans have better sparing of the lung and breast tissue compared to photons. Also, the potential for end of range RBE effects are unlikely to be large enough to offset the dosimetric advantages of proton therapy in cardiac substructures for mediastinal targets, although these dosimetric findings require validation with late toxicity data.

Dosimetric comparison of advanced radiotherapy approaches using photon techniques and particle therapy in the postoperative management of thymoma

BACKGROUND

The purpose of this study was to compare dosimetric differences related to target volume and organs-at-risk (OAR) using 3D-conformal radiotherapy (3DCRT), volumetric modulated arc therapy (VMAT), TomoTherapy (Tomo), proton radiotherapy (PRT), and carbon ion radiotherapy (CIRT) as part of postoperative thymoma irradiation.

MATERIAL AND METHODS

This single-institutional analysis included 10 consecutive patients treated with adjuvant radiotherapy between December 2013 and September 2016. CT-datasets and respective RT-structures were anonymized and plans for all investigated RT modalities (3DCRT, VMAT, Tomo, PRT, CIRT) were optimized for a total dose of 50 Gy in 25 fractions. Comparisons between target volume and OAR dosimetric parameters were performed using the Wilcoxon rank-sum test.

RESULTS

The best target volume coverage (mean PTV V_95% for all patients) was observed for Tomo (97.9%), PRT (97.6%), and CIRT (96.6%) followed by VMAT (85.4%) and 3DCRT (74.7%). PRT and CIRT both significantly reduced mean doses to the lungs, breasts, heart, and esophagus, as well as the spinal cord maximum dose compared with photon modalities. Among photon-based techniques, VMAT showed improved OAR sparing over 3DCRT. Tomo was associated with considerable low-dose exposure to the lungs, breasts, and heart.

CONCLUSIONS

Particle radiotherapy (PRT, CIRT) showed superior OAR sparing and optimal target volume coverage. The observed dosimetric advantages are expected to reduce toxicity rates. However, their clinical impact must be investigated prospectively.

Technical and dosimetric implications of respiratory induced density variations in a heterogeneous lung phantom

BACKGROUND

Stereotactic Body Radiotherapy (SBRT) is an ablative dose delivery technique which requires the highest levels of precision and accuracy. Modeling dose to a lung treatment volume has remained a complex and challenging endeavor due to target motion and the low density of the surrounding media. When coupled together, these factors give rise to pulmonary induced tissue heterogeneities which can lead to inaccuracies in dose computation. This investigation aims to determine which combination of imaging techniques and computational algorithms best compensates for time dependent lung target displacements.

METHODS

A Quasar phantom was employed to simulate respiratory motion for target ranges up to 3 cm. 4DCT imaging was used to generate Average Intensity Projection (AIP), Free Breathing (FB), and Maximum Intensity Projection (MIP) image sets. In addition, we introduce and compare a fourth dataset for dose computation based on a novel phase weighted density (PWD) technique. All plans were created using Eclipse version 13.6 treatment planning system and calculated using the Analytical Anisotropic Algorithm and Acuros XB. Dose delivery was performed using Truebeam STx linear accelerator where radiochromic film measurements were accessed using gamma analysis to compare planned versus delivered dose.

RESULTS

In the most extreme case scenario, the mean CT difference between FB and MIP datasets was found to be greater than 200 HU. The near maximum dose discrepancies between AAA and AXB algorithms were determined to be marginal (< 2.2%), with a greater variability occurring within the near minimum dose regime (< 7%). Radiochromatic film verification demonstrated all AIP and FB based computations exceeded 98% passing rates under conventional radiotherapy tolerances (gamma 3%, 3 mm). Under more stringent SBRT tolerances (gamma 3%, 1 mm), the AIP and FB based treatment plans exhibited higher pass rates (> 85%) when compared to MIP and PWD (< 85%) for AAA computations. For AXB, however, the delivery accuracy for all datasets were greater than 85% (gamma 3%,1 mm), with a corresponding reduction in overall lung irradiation.

CONCLUSIONS

Despite the substantial density variations between computational datasets over an extensive range of target movement, the dose difference between CT datasets is small and could not be quantified with ion chamber. Radiochromatic film analysis suggests the optimal CT dataset is dependent on the dose algorithm used for evaluation. With AAA, AIP and FB resulted in the best conformance between measured versus calculated dose for target motion ranging up to 3 cm under both conventional and SBRT tolerance criteria. With AXB, pass rates improved for all datasets with the PWD technique demonstrating slightly better conformity over AIP and FB based computations (gamma 3%, 1 mm). As verified in previous studies, our results confirm a clear advantage in delivery accuracy along with a relative decrease in calculated dose to the lung when using Acuros XB over AAA.

Advanced radiation techniques for locally advanced non-small cell lung cancer: intensity-modulated radiation therapy and proton therapy

Radiation therapy (RT) represents an integral part of a multimodality treatment plan in the definitive, preoperative and postoperative management of non-small cell lung cancer (NSCLC). Technological advances in RT have enabled a shift from two-dimensional radiotherapy to more conformal techniques. Three-dimensional conformal radiotherapy (3DCRT), the current minimum technological standard for treating NSCLC, allows for more accurate delineation of tumor burden by using computed tomography-based treatment planning instead of two-dimensional radiographs. Intensity-modulated RT (IMRT) and proton therapy represent advancements over 3DCRT that aim to improve the conformity of RT and provide the possibility for dose escalation to the tumor by minimizing radiation dose to organs at risk. Both techniques likely confer benefits to certain anatomic subgroups of NSCLC requiring RT. This article reviews pertinent studies evaluating the use of IMRT and proton therapy in locally advanced NSCLC, and outlines challenges, indications for use, and areas for future research.

G4DARI: Geant4/GATE based Monte Carlo simulation interface for dosimetry calculation in radiotherapy

Monte Carlo (MC) simulation is widely recognized as an important technique to study the physics of particle interactions in nuclear medicine and radiation therapy. There are different codes dedicated to dosimetry applications and widely used today in research or in clinical application, such as MCNP, EGSnrc and Geant4. However, such codes made the physics easier but the programming remains a tedious task even for physicists familiar with computer programming. In this paper we report the development of a new interface GEANT4 Dose And Radiation Interactions (G4DARI) based on GEANT4 for absorbed dose calculation and for particle tracking in humans, small animals and complex phantoms. The calculation of the absorbed dose is performed based on 3D CT human or animal images in DICOM format, from images of phantoms or from solid volumes which can be made from any pure or composite material to be specified by its molecular formula. G4DARI offers menus to the user and tabs to be filled with values or chemical formulas. The interface is described and as application, we show results obtained in a lung tumor in a digital mouse irradiated with seven energy beams, and in a patient with glioblastoma irradiated with five photon beams. In conclusion, G4DARI can be easily used by any researcher without the need to be familiar with computer programming, and it will be freely available as an application package.

Advanced Proton Beam Dosimetry Part I: review and performance evaluation of dose calculation algorithms

The accuracy of dose calculation is vital to the quality of care for patients undergoing proton beam therapy (PBT). Currently, the dose calculation algorithms available in commercial treatment planning systems (TPS) in PBT are classified into two classes: pencil beam (PB) and Monte-Carlo (MC) algorithms. PB algorithms are still regarded as the standard of practice in PBT, but they are analytical approximations whereas MC algorithms use random sampling of interaction cross-sections that represent the underlying physics to simulate individual particles trajectories. This article provides a brief review of PB and MC dose calculation algorithms employed in commercial treatment planning systems and their performance comparison in phantoms through simulations and measurements. Deficiencies of PB algorithms are first highlighted by a simplified simulation demonstrating the transport of a single sub-spot of proton beam that is incident at an oblique angle in a water phantom. Next, more typical cases of clinical beams in water phantom are presented and compared to measurements. The inability of PB to correctly predict the range and subsequently distal fall-off is emphasized. Through the presented examples, it is shown how dose errors as high as 30% can result with use of a PB algorithm. These dose errors can be minimized to clinically acceptable levels of less than 5%, if MC algorithm is employed in TPS. As a final illustration, comparison between PB and MC algorithm is made for a clinical beam that is use to deliver uniform dose to a target in a lung section of an anthropomorphic phantom. It is shown that MC algorithm is able to correctly predict the dose at all depths and matched with measurements. For PB algorithm, there is an increasing mismatch with the measured doses with increasing tissue heterogeneity. The findings of this article provide a foundation for the second article of this series to compare MC vs. PB based lung cancer treatment planning.

Advanced proton beam dosimetry part II: Monte Carlo vs. pencil beam-based planning for lung cancer

BACKGROUND

Proton pencil beam (PB) dose calculation algorithms have limited accuracy within heterogeneous tissues of lung cancer patients, which may be addressed by modern commercial Monte Carlo (MC) algorithms. We investigated clinical pencil beam scanning (PBS) dose differences between PB and MC-based treatment planning for lung cancer patients.

METHODS

With IRB approval, a comparative dosimetric analysis between RayStation MC and PB dose engines was performed on ten patient plans. PBS gantry plans were generated using single-field optimization technique to maintain target coverage under range and setup uncertainties. Dose differences between PB-optimized (PBopt), MC-recalculated (MCrecalc), and MC-optimized (MCopt) plans were recorded for the following region-of-interest metrics: clinical target volume (CTV) V95, CTV homogeneity index (HI), total lung V20, total lung VRX (relative lung volume receiving prescribed dose or higher), and global maximum dose. The impact of PB-based and MC-based planning on robustness to systematic perturbation of range (±3% density) and setup (±3 mm isotropic) was assessed. Pairwise differences in dose parameters were evaluated through non-parametric Friedman and Wilcoxon sign-rank testing.

RESULTS

In this ten-patient sample, CTV V95 decreased significantly from 99-100% for PBopt to 77-94% for MCrecalc and recovered to 99-100% for MCopt (P<10-5). The median CTV HI (D95/D5) decreased from 0.98 for PBopt to 0.91 for MCrecalc and increased to 0.95 for MCopt (P<10-3). CTV D95 robustness to range and setup errors improved under MCopt (ΔD95 =-1%) compared to MCrecalc (ΔD95 =-6%, P=0.006). No changes in lung dosimetry were observed for large volumes receiving low to intermediate doses (e.g., V20), while differences between PB-based and MC-based planning were noted for small volumes receiving high doses (e.g., VRX). Global maximum patient dose increased from 106% for PBopt to 109% for MCrecalc and 112% for MCopt (P<10-3).

CONCLUSIONS

MC dosimetry revealed a reduction in target dose coverage under PB-based planning that was regained under MC-based planning along with improved plan robustness. MC-based optimization and dose calculation should be integrated into clinical planning workflows of lung cancer patients receiving actively scanned proton therapy.

Assessment of radiation-induced secondary cancer risk in the Brazilian population from left-sided breast-3D-CRT using MCNPX

OBJECTIVE

The aim of this work was to simulate a 6MV conventional breast 3D conformational radiation therapy (3D-CRT) with physical wedges (50 Gy/25#) in the left breast, calculate the mean absorbed dose in the body organs using robust models and computational tools and estimate the secondary cancer-incidence risk to the Brazilian population.

METHODS

The VW female phantom was used in the simulations. Planning target volume (PTV) was defined in the left breast. The 6MV parallel-opposed fields breast-radiotherapy (RT) protocol was simulated with MCNPx code. The absorbed doses were evaluated in all the organs. The secondary cancer-incidence risk induced by radiotherapy was calculated for different age groups according to the BEIR VII methodology.

RESULTS

RT quality indexes indicated that the protocol was properly simulated. Significant absorbed dose values in red bone marrow, RBM (0.8 Gy) and stomach (0.6 Gy) were observed. The contralateral breast presented the highest risk of incidence of a secondary cancer followed by leukaemia, lung and stomach. The risk of a secondary cancer-incidence by breast-RT, for the Brazilian population, ranged between 2.2-1.7% and 0.6-0.4%.

CONCLUSION

RBM and stomach, usually not considered as OAR, presented high second cancer incidence risks of 0.5-0.3% and 0.4-0.1%, respectively. This study may be helpful for breast-RT risk/benefit assessment. Advances in knowledge: MCNPX-dosimetry was able to provide the scatter radiation and dose for all body organs in conventional breast-RT. It was found a relevant risk up to 2.2% of induced-cancer from breast-RT, considering the whole thorax organs and Brazilian cancer-incidence.

Dosimetric characterization of 3D printed bolus at different infill percentage for external photon beam radiotherapy

BACKGROUND AND PURPOSE

3D printing is rapidly evolving and further assessment of materials and technique is required for clinical applications. We evaluated 3D printed boluses with acrylonitrile butadiene styrene (ABS) and polylactide (PLA) at different infill percentage.

MATERIAL AND METHODS

A low-cost 3D printer was used. The influence of the air inclusion within the 3D printed boluses was assessed thoroughly both with treatment planning system (TPS) and with physical measurements. For each bolus, two treatment plans were calculated with Monte Carlo algorithm, considering the computed tomography (CT) scan of the 3D printed bolus or modelling the 3D printed bolus as a virtual bolus structure with a homogeneous density. Depth dose measurements were performed with Gafchromic films.

RESULTS

High infill percentage corresponds to high density and high homogeneity within bolus material. The approximation of the bolus in the TPS as a homogeneous material is satisfying for infill percentages greater than 20%. Measurements performed with PLA boluses are more comparable to the TPS calculated profiles. For boluses printed at 40% and 60% infill, the discrepancies between calculated and measured dose distribution are within 5%.

CONCLUSIONS

3D printing technology allows modulating the shift of the build-up region by tuning the infill percentage of the 3D printed bolus in order to improve superficial target coverage.