Monte carlo evaluation of the AAA treatment planning algorithm in a heterogeneous multilayer phantom and IMRT clinical treatments for an Elekta SL25 linear accelerator

The accuracy of dose calculations by anisotropic analytical algorithms for stereotactic radiotherapy in nasopharyngeal carcinoma

Nasopharyngeal tumors are commonly treated with intensity-modulated radiotherapy techniques. For photon dose calculations, problems related to loss of lateral electronic equilibrium exist when small fields are used. The anisotropic analytical algorithm (AAA) implemented in Varian Eclipse was developed to replace the pencil beam convolution (PBC) algorithm for more accurate dose prediction in an inhomogeneous medium. The purpose of this study was to investigate the accuracy of the AAA for predicting interface doses for intensity-modulated stereotactic radiotherapy boost of nasopharyngeal tumors. The central axis depth dose data and dose profiles of phantoms with rectangular air cavities for small fields were measured using a 6 MV beam. In addition, the air-tissue interface doses from six different intensity-modulated stereotactic radiotherapy plans were measured in an anthropomorphic phantom. The nasopharyngeal region of the phantom was especially modified to simulate the air cavities of a typical patient. The measured data were compared to the data calculated by both the AAA and the PBC algorithm. When using single small fields in rectangular air cavity phantoms, both AAA and PBC overestimated the central axis dose at and beyond the first few millimeters of the air-water interface. Although the AAA performs better than the PBC algorithm, its calculated interface dose could still be more than three times that of the measured dose when a 2 × 2 cm(2) field was used. Testing of the algorithms using the anthropomorphic phantom showed that the maximum overestimation by the PBC algorithm was 20.7%, while that by the AAA was 8.3%. When multiple fields were used in a patient geometry, the dose prediction errors of the AAA would be substantially reduced compared with those from a single field. However, overestimation of more than 3% could still be found at some points at the air-tissue interface.

Dosimetric verification of the anisotropic analytical algorithm in lung equivalent heterogeneities with and without bone equivalent heterogeneities

PURPOSE

In this study, the authors evaluated the accuracy of dose calculations performed by the convolution/superposition based anisotropic analytical algorithm (AAA) in lung equivalent heterogeneities with and without bone equivalent heterogeneities.

METHODS

Calculations of PDDs using the AAA and Monte Carlo simulations (MCNP4C) were compared to ionization chamber measurements with a heterogeneous phantom consisting of lung equivalent and bone equivalent materials. Both 6 and 10 MV photon beams of 4 x 4 and 10 x 10 cm(2) field sizes were used for the simulations. Furthermore, changes of energy spectrum with depth for the heterogeneous phantom using MCNP were calculated.

RESULTS

The ionization chamber measurements and MCNP calculations in a lung equivalent phantom were in good agreement, having an average deviation of only 0.64 +/- 0.45%. For both 6 and 10 MV beams, the average deviation was less than 2% for the 4 x 4 and 10 x 10 cm(2) fields in the water-lung equivalent phantom and the 4 x 4 cm(2) field in the water-lung-bone equivalent phantom. Maximum deviations for the 10 x 10 cm(2) field in the lung equivalent phantom before and after the bone slab were 5.0% and 4.1%, respectively. The Monte Carlo simulation demonstrated an increase of the low-energy photon component in these regions, more for the 10 X 10 cm(2) field compared to the 4 x 4 cm(2) field.

CONCLUSIONS

The low-energy photon by Monte Carlo simulation component increases sharply in larger fields when there is a significant presence of bone equivalent heterogeneities. This leads to great changes in the build-up and build-down at the interfaces of different density materials. The AAA calculation modeling of the effect is not deemed to be sufficiently accurate.

Dosimetry of oblique tangential photon beams calculated by superposition/convolution algorithms: a Monte Carlo evaluation

Although there are many works on evaluating dose calculations of the anisotropic analytical algorithm (AAA) using various homogeneous and heterogeneous phantoms, related work concerning dosimetry due to tangential photon beam is lacking. In this study, dosimetry predicted by the AAA and collapsed cone convolution (CCC) algorithm was evaluated using the tangential photon beam and phantom geometry. The photon beams of 6 and 15 MV with field sizes of 4 × 4 (or 7 × 7), 10 × 10 and 20 × 20 cm², produced by a Varian 21 EX linear accelerator, were used to test performances of the AAA and CCC using Monte Carlo (MC) simulation (EGSnrc-based code) as a benchmark. Horizontal dose profiles at different depths, phantom skin profiles (i.e., vertical dose profiles at a distance of 2 mm from the phantom lateral surface), gamma dose distributions, and dose-volume histograms (DVHs) of skin slab were determined. For dose profiles at different depths, the CCC agreed better with doses in the air-phantom region, while both the AAA and CCC agreed well with doses in the penumbra region, when compared to the MC. Gamma evaluations between the AAA/CCC and MC showed that deviations of 2D dose distribution occurred in both beam edges in the phantom and air-phantom interface. Moreover, the gamma dose deviation is less significant in the air-phantom interface than the penumbra. DVHs of skin slab showed that both the AAA and CCC underestimated the width of the dose drop-off region for both the 6 and 15 MV photon beams. When the gantry angle was 0°, it was found that both the AAA and CCC overestimated doses in the phantom skin profiles compared to the MC, with various photon beam energies and field sizes. The mean dose differences with doses normalized to the prescription point for the AAA and CCC were respectively: 7.6% ± 2.6% and 2.1% ± 1.3% for a 10 × 10 cm2 field, 6 MV; 16.3%± 2.1% and 6.7% ± 2.1% for a 20 × 20 cm2 field, 6 MV; 5.5% ± 1.2% and 1.7% ± 1.4% for a 10 × 10 cm2, 15 MV; 18.0% ± 1.3% and 8.3% ± 1.8% for a 20 × 20 cm², 15 MV. However, underestimations of doses in the phantom skin profile were found with small fields of 4 × 4 and 7 × 7 cm² for the 6 and 15 MV photon beams, respectively, when the gantry was turned 5° anticlockwise. As surface dose with tangential photon beam geometry is important in some radiation treatment sites such as breast, chest wall and sarcoma, it is found that neither of the treatment planning system algorithms can predict the dose well at depths shallower than 2 mm. The dosimetry data and beam and phantom geometry in this study provide a better knowledge of a dose calculation algorithm in tangential-like irradiation.

Intensity-modulated radiotherapy for locally advanced non-small-cell lung cancer: a dose-escalation planning study

PURPOSE

To evaluate the potential for dose escalation with intensity-modulated radiotherapy (IMRT) in positron emission tomography-based radiotherapy planning for locally advanced non-small-cell lung cancer (LA-NSCLC).

METHODS AND MATERIALS

For 35 LA-NSCLC patients, three-dimensional conformal radiotherapy and IMRT plans were made to a prescription dose (PD) of 66 Gy in 2-Gy fractions. Dose escalation was performed toward the maximal PD using secondary endpoint constraints for the lung, spinal cord, and heart, with de-escalation according to defined esophageal tolerance. Dose calculation was performed using the Eclipse pencil beam algorithm, and all plans were recalculated using a collapsed cone algorithm. The normal tissue complication probabilities were calculated for the lung (Grade 2 pneumonitis) and esophagus (acute toxicity, grade 2 or greater, and late toxicity).

RESULTS

IMRT resulted in statistically significant decreases in the mean lung (p <.0001) and maximal spinal cord (p = .002 and 0005) doses, allowing an average increase in the PD of 8.6-14.2 Gy (p ≤.0001). This advantage was lost after de-escalation within the defined esophageal dose limits. The lung normal tissue complication probabilities were significantly lower for IMRT (p <.0001), even after dose escalation. For esophageal toxicity, IMRT significantly decreased the acute NTCP values at the low dose levels (p = .0009 and p <.0001). After maximal dose escalation, late esophageal tolerance became critical (p <.0001), especially when using IMRT, owing to the parallel increases in the esophageal dose and PD.

CONCLUSION

In LA-NSCLC, IMRT offers the potential to significantly escalate the PD, dependent on the lung and spinal cord tolerance. However, parallel increases in the esophageal dose abolished the advantage, even when using collapsed cone algorithms. This is important to consider in the context of concomitant chemoradiotherapy schedules using IMRT.

Breast in vivo dosimetry by EPID

An electronic portal imaging device (EPID) is an effective detector for in vivo transit dosimetry. In fact, it supplies two‐dimensional information, does not require special efforts to be used during patient treatment, and can supply data in real time. In the present paper, a new procedure has been proposed to improve the EPID in vivo dosimetry accuracy by taking into account the patient setup variations. The procedure was applied to the breast tangential irradiation for the reconstruction of the dose at the breast midpoint, Dm. In particular, the patient setup variations were accounted for by comparing EPID images versus digitally reconstructed radiographies. In this manner, EPID transit signals were obtained corresponding to the geometrical projections of the breast midpoint on the EPID for each therapy session. At the end, the ratios R between Dm and the doses computed by the treatment planning system (TPS) at breast midpoints, Dm,TPS, were determined for 800 therapy sessions of 20 patients. Taking into account the method uncertainty, tolerance levels equal to ±5% have been determined for the ratio R.

The improvement of in vivo dosimetry results obtained (taking into account patient misalignment) has been pointed out comparing the R values obtained with and without considering patient setup variations. In particular, when patient misalignments were taken into account, the R values were within ± 5% for 93% of the checks; when patient setup variations were not taken into account, the R values were within ± 5% in 72% of the checks. This last result points out that the transit dosimetry method overestimates the dose discrepancies if patient setup variations are not taken into account for dose reconstruction. In this case, larger tolerance levels have to be adopted as a trade‐off between workload and ability to detect errors, with the drawback being that some errors (such as the ones in TPS implementation or in beam calibration) cannot be detected, limiting the in vivo dosimetry efficacy.

The paper also reports preliminary results about the possibility of reconstructing a dose profile perpendicular to the beam central axis reaching from the apex to the lung and passing through the middle point of the breast by an algorithm, similar to the one used for dose reconstruction at breast midpoint. In particular, the results have shown an accuracy within ± 3% for the dose profile reconstructed in the breast (excluding the interface regions) and an underestimation of the lung dose.

PACS numbers: 87.55.Qr, 87.55.km, 87.53.Bn

Verification of IMRT dose calculations using AAA and PBC algorithms in dose buildup regions

The purpose of this comparative study was to test the accuracy of anisotropic analytical algorithm (AAA) and pencil beam convolution (PBC) algorithms of Eclipse treatment planning system (TPS) for dose calculations in the low- and high-dose buildup regions. AAA and PBC algorithms were used to create two intensity-modulated radiotherapy (IMRT) plans of the same optimal fluence generated from a clinically simulated oropharynx case in an in-house fabricated head and neck phantom. The TPS computed buildup doses were compared with the corresponding measured doses in the phantom using thermoluminescence dosimeters (TLD 100). Analysis of dose distribution calculated using PBC and AAA shows an increase in gamma value in the dose buildup region indicating large dose deviation. For the surface areas of 1, 50 and 100 cm2, PBC overestimates doses as compared to AAA calculated value in the range of 1.34%-3.62% at 0.6 cm depth, 1.74%-2.96% at 0.4 cm depth, and 1.96%-4.06% at 0.2 cm depth, respectively. In high-dose buildup region, AAA calculated doses were lower by an average of -7.56% (SD = 4.73%), while PBC was overestimated by 3.75% (SD = 5.70%) as compared to TLD measured doses at 0.2 cm depth. However, at 0.4 and 0.6 cm depth, PBC overestimated TLD measured doses by 5.84% (SD = 4.38%) and 2.40% (SD = 4.63%), respectively, while AAA underestimated the TLD measured doses by -0.82% (SD = 4.24%) and -1.10% (SD = 4.14%) at the same respective depth. In low-dose buildup region, both AAA and PBC overestimated the TLD measured doses at all depths except -2.05% (SD = 10.21%) by AAA at 0.2 cm depth. The differences between AAA and PBC at all depths were statistically significant (p < 0.05) in high-dose buildup region, whereas it is not statistically significant in low-dose buildup region. In conclusion, AAA calculated the dose more accurately than PBC in clinically important high-dose buildup region at 0.4 cm and 0.6 cm depths. The use of an orfit cast increases the dose buildup effect, and this buildup effect decreases with depth.

[Comparison of RTPS and Monte Carlo dose distributions in heterogeneous phantoms for photon beams]

The purpose of this study was to compare dose distributions from three different RTPS with those from Monte Carlo (MC) calculations and measurements, in heterogeneous phantoms for photon beams. This study used four algorithms for RTPS: AAA (analytical anisotropic algorithm) implemented in the Eclipse (Varian Medical Systems) treatment planning system, CC (collapsed cone) superposition from the Pinnacle (Philips), and MGS (multigrid superposition) and FFT (fast Fourier transform) convolution from XiO (CMS). The dose distributions from these algorithms were compared with those from MC and measurements in a set of heterogeneous phantoms. Eclipse/AAA underestimated the dose inside the lung region for low energies of 4 and 6 MV. This is because Eclipse/AAA do not adequately account for a scaling of the spread of the pencil (lateral electron transport) based on changes in the electron density at low photon energies. The dose distributions from Pinnacle/CC and XiO/MGS almost agree with those of MC and measurements at low photon energies, but increase errors at high energy of 15 MV, especially for a small field of 3x3 cm(2). The FFT convolution extremely overestimated the dose inside the lung slab compared to MC. The dose distributions from the superposition algorithms almost agree with those from MC as well as measured values at 4 and 6 MV. The dose errors for Eclipse/AAA are lager in lung model phantoms for 4 and 6 MV. It is necessary to use the algorithms comparable to superposition for accuracy of dose calculations in heterogeneous regions.

Variations of lung density and geometry on inhomogeneity correction algorithms: a Monte Carlo dosimetric evaluation

This work contributed the following new information to the study of inhomogeneity correction algorithm: (1) Evaluation of lung dose calculation methods as a function of lung relative electron density (rhoe,lung) and treatment geometry and (2) comparison of doses calculated using the collapsed cone convolution (CCC) and adaptive convolution (AC) in lung using the Monte Carlo (MC) simulation with the EGSnrc-based code. The variations of rhoe,lung and geometry such as the position and dimension of the lung were studied with different photon beam energies and field sizes. Three groups of inhomogeneous lung phantoms, namely, "slab," "column," and "cube," with different positions, volumes, and shapes of lung in water as well as clinical computed tomography lung images were used. The rhoe,lung in each group of phantoms vary from 0.05 to 0.7. 6 and 18 MV photon beams with small (4 x 4 cm2) and medium (10 x 10 cm2) field sizes produced by a Varian 21 EX linear accelerator were used. This study reveals that doses in the inhomogeneous lung calculated by the CCC match well with those by AC within +/- 1%, indicating that the AC, with an advantage of shorter computing times (three to four times shorter than CCC), is a good substitute for CCC. Comparing the CCC and AC to MC in general, significant dose deviations are found when the rhoe,lung is < or =0.3. The degree of deviation depends on the photon beam energy and field size and is relatively large when high-energy photon beams with small fields are used. For penumbra widths (20%-80%), the CCC and AC agree well with MC for the slab and cube phantoms with the lung volumes at the central beam axis (CAX). However, deviations (>2 mm) occur in the column phantoms, with two lung volumes separated by a unit density column along the CAX in the middle using the 18 MV beam with 4 x 4 cm2 field for rhoe,lung < or =0.1. This study provides new dosimetric data to evaluate the impact of the variations of rhoe,lung and geometry on dose calculations in inhomogeneous media using CCC and AC.

Technical note: Heterogeneity dose calculation accuracy in IMRT: study of five commercial treatment planning systems using an anthropomorphic thorax phantom

The purpose of this study was to determine the accuracy of five commonly used intensity-modulated radiation therapy (IMRT) treatment planning systems (TPSs), 3 using convolution superposition algorithms or the analytical anisotropic algorithm (CSA/AAAs) and 2 using pencil beam algorithms (PBAs), in calculating the absorbed dose within a low-density, heterogeneous region when compared with measurements made in an anthropomorphic thorax phantom. The dose predicted in the target center met the test criteria (5% of the dose normalization point or 3 mm distance to agreement) for all TPSs tested; however, at the tumor-lung interface and at the peripheral lung in the vicinity of the tumor, the CSA/AAAs performed better than the PBAs (85% and 50%, respectively, of pixels meeting the 5%/3-mm test criteria), and thus should be used to determine dose in heterogeneous regions.

Treatment planning comparison between conformal radiotherapy and helical tomotherapy in the case of locally advanced-stage NSCLC

BACKGROUND AND PURPOSE

To investigate the impact of Helical Tomotherapy (HT) upon the dose distribution when compared to our routinely delivered 3D conformal radiotherapy (CRT) in the case of patients affected by stage III non-small-cell lung cancer (NSCLC).

MATERIAL AND METHODS

Thirteen stage III inoperable NSCLC patients were scheduled to receive 61.2-70.2Gy, 1.8Gy/fraction. Two treatment techniques (HT and CRT) were considered, and in the case of CRT the dose calculation was performed using both the pencil beam (PB) and Anisotropic Analytical Algorithm (AAA) available on the Varian Eclipse planning system. Dose volume constraints for PTV coverage and OAR sparing were assessed for the HT inverse planning with the highest priority upon PTV coverage and spinal cord sparing. The three plans were compared in terms of dose-volume histograms (DVHs) and normal tissue complication probability (NTCP). A statistical analysis was performed using non-parametric Wilcoxon matched pairs tests.

RESULTS

In CRT the use of a less accurate algorithm (PB) decreased the monitor unit number by 2.4%. HT significantly improved dose homogeneity within PTV compared with CRT_AAA. For lung parenchyma V20-V40 were lower with HT, corresponding to a decrease of 7% in the risk of radiation pneumonitis. The volume of the heart and esophagus irradiated to >45-60Gy were reduced using HT plans. For eight PTs with an esophagus-PTV overlap >5%, HT significantly reduced both late and acute esophageal complication probability.

CONCLUSIONS

Our findings obtained in stage III NSCLC patients underline that HT guarantees an important sparing of lungs and esophagus, thus HT has the potential to improve therapeutic ratio, when compared with CRT, by means of dose escalation and/or combined treatment strategy. In CRT of locally advanced lung cancers, the use of a more advanced algorithm would give significantly better modeling of target dose and coverage.

A 3D pencil-beam-based superposition algorithm for photon dose calculation in heterogeneous media

In this work, a novel three-dimensional superposition algorithm for photon dose calculation is presented. The dose calculation is performed as a superposition of pencil beams, which are modified based on tissue electron densities. The pencil beams have been derived from Monte Carlo simulations, and are separated into lateral and depth-directed components. The lateral component is modeled using exponential functions, which allows accurate modeling of lateral scatter in heterogeneous tissues. The depth-directed component represents the total energy deposited on each plane, which is spread out using the lateral scatter functions. Finally, convolution in the depth direction is applied to account for tissue interface effects. The method can be used with the previously introduced multiple-source model for clinical settings. The method was compared against Monte Carlo simulations in several phantoms including lung- and bone-type heterogeneities. Comparisons were made for several field sizes for 6 and 18 MV energies. The deviations were generally within (2%, 2 mm) of the field central axis d(max). Significantly larger deviations (up to 8%) were found only for the smallest field in the lung slab phantom for 18 MV. The presented method was found to be accurate in a wide range of conditions making it suitable for clinical planning purposes.