Dosimetric accuracy of tomotherapy dose calculation in thorax lesions

Implementation and Challenges of International Atomic Energy Agency/American Association of Physicists in Medicine TRS 483 Formalism for Field Output Factors and Involved Uncertainties Determination in Small Fields for TomoTherapy

PURPOSE

International Atomic Energy Agency published TRS-483 to address the issues of small field dosimetry. Our study calculates the output factor in the small fields of TomoTherapy using different detectors and dosimetric conditions. Furthermore, it estimates the various components of uncertainty and presents challenges faced during implementation.

MATERIALS AND METHODS

Beam quality TPR20,10(10) at the hypothetical field size of 10 cm × 10 cm was calculated from TPR20,10(S). Two ionization chambers based on the minimum field width required to satisfy the lateral charge particle equilibrium and one unshielded electron field diode (EFD) were selected. Output factor measurements were performed in various dosimetric conditions.

RESULTS

Beam quality TPR20,10(10) has a mean value of 0.627 ± 0.001. The maximum variation of output factor between CC01 chamber and EFD diode at the smallest field size was 11.80%. In source to surface setup, the difference between water and virtual water was up to 9.68% and 8.13%, respectively, for the CC01 chamber and EFD diode. The total uncertainty in the ionization chamber was 2.43 times higher compared to the unshielded EFD diode at the smallest field size.

CONCLUSIONS

Beam quality measurements, chamber selection procedure, and output factors were successfully carried out. A difference of up to 10% in output factor can occur if density scaling for electron density in virtual water is not considered. The uncertainty in output correction factors dominates, while positional and meter reading uncertainty makes a minor contribution to total uncertainty. An unshielded EFD diode is a preferred detector in small fields because of lower uncertainty in measurements compared to ionization chambers.

Skin dose calculation during radiotherapy of head and neck cancer using deformable image registration of planning and mega-voltage computed tomography scans

Background and Purpose

Head-Neck (HN) patients may experience severe acute skin complications that can cause treatment interruption and increase the risk of late fibrosis. This study assessed a method for accurately monitoring skin dose changes during helical tomotherapy for HN cancer based on deformable image registration of planning computed tomography (CT) and mega-voltage CT (MVCT).

Materials and Methods

Planning CTs of nine patients were deformably registered to mid-treatment MVCT (MV15) images resulting in CTdef images. The original plans were recalculated on both CTdef and mid-treatment kilo-voltage CT (CT15) taken as ground truth. Superficial layers (SL) of the body with thicknesses of 2, 3 and 5 mm (SL2, SL3, SL5) were considered as derma surrogates. SL V95%, V97%, V98%, V100%, V102%, V105% and V107% of the prescribed PTV dose were extracted for CT15/CTdef and compared (considering patients with skin dose > 95%). For comparison, doses were calculated directly on the calibrated MVCT and analyzed in the same way.

Results

Differences between SL2/SL3/SL5 V95%-V107% in CT15/CTdef were very small: for eight of nine patients the difference between the considered SL2 Vd% computed on CTdef and CT15 was less than 1.4 cm³ for all d%. A larger value was found when using MVCT for skin dose calculation (4.8 cm³ for SL2), although CTdef body contour matched CT15 body with accuracy similar to that of MV15.

Conclusions

Deforming the planning CT-to-MVCT was shown to be accurate considering external body contours and skin DVHs. The method was able to accurately identify superficial overdosing.

Voxel-by-voxel correlation between radiologically radiation induced lung injury and dose after image-guided, intensity modulated radiotherapy for lung tumors

PURPOSE

To correlate radiation dose to the risk of severe radiologically-evident radiation-induced lung injury (RRLI) using voxel-by-voxel analysis of the follow-up computed tomography (CT) of patients treated for lung cancer with hypofractionated helical Tomotherapy.

METHODS AND MATERIALS

The follow-up CT scans from 32 lung cancer patients treated with various regimens (5, 8, and 25 fractions) were registered to pre-treatment CT using deformable image registration (DIR). The change in density was calculated for each voxel within the combined lungs minus the planning target volume (PTV). Parameters of a Probit formula were derived by fitting the occurrences of changes of density in voxels greater than 0.361gcm^-3 to the radiation dose. The model's predictive capability was assessed using the area under receiver operating characteristic curve (AUC), the Kolmogorov-Smirnov test for goodness-of-fit, and the permutation test (P_test).

RESULTS

The best-fit parameters for prediction of RRLI 6months post RT were D₅₀ of 73.0 (95% CI 59.2.4-85.3.7)Gy, and m of 0.41 (0.39-0.46) for hypofractionated (5 and 8 fractions) and D₅₀ of 96.8 (76.9-123.9)Gy, and m of 0.36 (0.34-0.39) for 25 fractions RT. According to the goodness-of-fit test the null hypothesis of modeled and observed occurrence of RRLI coming from the same distribution could not be rejected. The AUC was 0.581 (0.575-0.583) for fractionated and 0.579 (0.577-0.581) for hypofractionated patients. The predictive models had AUC>upper 95% band of the P_test.

CONCLUSIONS

The correlation of voxel-by-voxel density increase with dose can be used as a support tool for differential diagnosis of tumor from benign changes in the follow-up of lung IMRT patients.

Are simple IMRT beams more robust against MLC error? Exploring the impact of MLC errors on planar quality assurance and plan quality for different complexity beams

This study investigated the impact of beam complexities on planar quality assurance and plan quality robustness by introducing MLC errors in intensity‐modulate radiation therapy. Forty patients' planar quality assurance (QA) plans were enrolled in this study, including 20 dynamic MLC (DMLC) IMRT plans and 20 static MLC (SMLC) IMRT plans. The total beam numbers were 150 and 160 for DMLC and SMLC, respectively. Six different magnitudes of MLC errors were introduced to these beams. Gamma pass rates were calculated by comparing error‐free fluence and error‐induced fluence. The plan quality variation was acquired by comparing PTV coverage. Eight complexity scores were calculated based on the beam fluence and the MLC sequence. The complexity scores include fractal dimension, monitor unit, modulation index, fluence map complexity, weighted average of field area, weighted average of field perimeter, and small aperture ratio (<5cm2 and<50 cm2). The Spearman's rank correlation coefficient was calculated to analyze the correlation between these scores and gamma pass rate and plan quality variation. For planar QA, the most significant complexity index was fractal dimension for DMLC (p=−0.40) and weighted segment area for SMLC (p=0.27) at low magnitude MLC error. For plan quality, the most significant complexity index was weighted segment perimeter for DMLC (p=0.56) and weighted segment area for SMLC (p=0.497) at low magnitude MLC error. The sensitivity of planar QA was weakly associated with the field complexity with low magnitude MLC error, but the plan quality robustness was associated with beam complexity. Plans with simple beams were more robust to MLC error.

PACS number(s): 87.55

Dosimetric verification and quality assurance of running-start-stop (RSS) delivery in tomotherapy

The purpose of this study was to evaluate the dosimetric profiles and delivery accuracy of running-start-stop (RSS) delivery in tomotherapy and to present initial quality assurance (QA) results on the accuracy of the dynamic jaw motion, dosimetric penumbrae of the RSS dynamic jaw and the static jaw were measured by radiographic films. Delivery accuracy of the RSS was evaluated by gamma analysis on film measurements of 12 phantom plans. Consistency in the performance of RSS was evaluated by QA procedures over the first nine months after the installation of the feature. These QA were devised to check: 1) positional accuracy of moving jaws; 2) consistency of relative radiation output collimated by discrete and continuously sweeping jaws; 3) consistency of field widths and profiles. In the longitudinal direction, the dose penumbra in RSS delivery was reduced from 17.3mm to 10.2 mm for 2.5 cm jaw, and from 33.2 mm to 9.6 mm for 5 cm jaw. Gamma analysis on the twelve plans revealed that over 90% of the voxels in the proximity of the penumbra region satisfied the gamma criteria of 2% dose difference and 2 mm distance-to-agreement. The initial QA results during the first nine months after installation of the RSS are presented. Jaw motion was shown to be accurate with maximum encoder error less than 0.42 mm. The consistency of relative output for discrete and continuously sweeping jaws was within 1.2%. Longitudinal radiation profiles agreed to the reference profile with maximum gamma < 1 and field width error < 1.8%. With the same jaw width, RSS showed better dose penumbrae compared to those from static jaw delivery. The initial QA results on the accuracy of moving jaws, reproducibility of dosimetric output and profiles were satisfactory.

Monte Carlo evaluation of the dose calculation algorithm of TomoTherapy for clinical cases in dynamic jaws mode

PURPOSE

For the TomoTherapy(®) system, longitudinal conformation can be improved by selecting a smaller field width but at the expense of longer treatment time. Recently, the TomoEdge(®) feature has been released with the possibility to move dynamically the jaws at the edges of the target volume, improving longitudinal penumbra and enabling faster treatments. Such delivery scheme requires additional modeling of treatment delivery. Using a previously validated Monte Carlo model (TomoPen), we evaluated the accuracy of the implementation of TomoEdge in the new dose engine of TomoTherapy for 15 clinical cases.

METHODS

TomoPen is based on PENELOPE. Particle tracking in the treatment head is performed almost instantaneously by 1) reading a particle from a phase-space file corresponding to the largest field and 2) correcting the weight of the particle depending on the actual jaw and MLC configurations using Monte Carlo pre-generated data. 15 clinical plans (5 head-and-neck, 5 lung and 5 prostate tumors) planned with TomoEdge and with the last release of the treatment planning system (VoLO(®)) were re-computed with TomoPen. The resulting dose-volume histograms were compared.

RESULTS

Good agreement was achieved overall, with deviations for the target volumes typically within 2% (D95), excepted for small lung tumors (17 cm(3)) where a maximum deviation of 4.4% was observed for D95. The results were consistent with previously reported values for static field widths.

CONCLUSIONS

For the clinical cases considered in the present study, the introduction of TomoEdge did not impact significantly the accuracy of the computed dose distributions.

Fast, simple, and informative patient-specific dose verification method for intensity modulated total marrow irradiation with helical tomotherapy

Background

Patient-specific dose verification for treatment planning in helical tomotherapy is routinely performed using a homogeneous virtual water cylindrical phantom of 30 cm diameter and 18 cm length (Cheese phantom). Because of this small length, treatment with total marrow irradiation (TMI) requires multiple deliveries of the dose verification procedures to cover a wide range of the target volumes, which significantly prolongs the dose verification process. We propose a fast, simple, and informative patient-specific dose verification method which reduce dose verification time for TMI with helical tomotherapy.

Methods

We constructed a two-step solid water slab phantom (length 110 cm, height 8 cm, and two-step width of 30 cm and 15 cm), termed the Whole Body Phantom (WB phantom). Three ionization chambers and three EDR-2 films can be inserted to cover extended field TMI treatment delivery. Three TMI treatment plans were conducted with a TomoTherapy HiArt Planning Station and verified using the WB phantom with ion chambers and films. Three regions simulating the head and neck, thorax, and pelvis were covered in a single treatment delivery. The results were compared to those with the cheese phantom supplied by Accuray, Inc. following three treatment deliveries to cover the body from head to pelvis.

Results

Use of the WB phantom provided point doses or dose distributions from head and neck to femur in a single treatment delivery of TMI. Patient-specific dose verification with the WB phantom was 62% faster than with the cheese phantom. The average pass rate in gamma analysis with the criteria of a 3-mm distance-to-agreement and 3% dose differences was 94% ± 2% for the three TMI treatment plans. The differences in pass rates between the WB and cheese phantoms at the upper thorax to abdomen regions were within 2%. The calculated dose agreed with the measured dose within 3% for all points in all five cases in both the WB and cheese phantoms.

Conclusions

Our dose verification method with the WB phantom provides simple and rapid quality assurance without limiting dose verification information in total marrow irradiation with helical tomotherapy.

Monte Carlo and ray tracing algorithms in the cyberknife treatment planning for lung tumours- comparison and validation

Multiplan treatment planning system, used with Cyberknife system, provides the option of using either the ray tracing algorithm or the Monte Carlo algorithm for the final dose calculation. In order to compare and validate the dose calculations of these algorithms, especially in a heterogeneous medium, a lung phantom study was carried out. Validation has been done with thermoluminiscent dosimetry (TLD) using lithium fluoride rods for the point doses and film dosimetry using EBT2 films for the dose distribution. In the point dose measurements, an agreement of 100.1+2.6 % (1 SD) is observed with the Monte Carlo dose calculation, whereas it is only 91.2+ 3.2% (1 SD) with the ray tracing calculation. On subjecting the dose distributions from irradiated EBT2 films for validation of Monte Carlo calculation MC , over 96% of the pixels pass the gamma criteria of 3mm and 3cGy.On analyzing the dose profiles from EBT2 films and the corresponding profiles from the plan calculated using the Monte Carlo algorithm, it is seen that the maximum distance-to-agreement values are within the 3mm criteria set, whereas the maximum values are as high as 8 mm when compared with plan calculated using ray tracing algorithm. The results of the actual measurements are more consistent with the dose calculation by the Monte Carlo algorithm.