- Open Access
Dose calculation of Acuros XB and Anisotropic Analytical Algorithm in lung stereotactic body radiotherapy treatment with flattening filter free beams and the potential role of calculation grid size
Radiation Oncology volume 10, Article number: 53 (2015)
The study aimed to appraise the dose differences between Acuros XB (AXB) and Anisotropic Analytical Algorithm (AAA) in stereotactic body radiotherapy (SBRT) treatment for lung cancer with flattening filter free (FFF) beams. Additionally, the potential role of the calculation grid size (CGS) on the dose differences between the two algorithms was also investigated.
SBRT plans with 6X and 10X FFF beams produced from the CT scan data of 10 patients suffering from stage I lung cancer were enrolled in this study. Clinically acceptable treatment plans with AAA were recalculated using AXB with the same monitor units (MU) and identical multileaf collimator (MLC) settings. Furthermore, different CGS (2.5 mm and 1 mm) in the two algorithms was also employed to investigate their dosimetric impact. Dose to planning target volumes (PTV) and organs at risk (OARs) between the two algorithms were compared. PTV was separated into PTV_soft (density in soft-tissue range) and PTV_lung (density in lung range) for comparison.
The dose to PTV_lung predicted by AXB was found to be 1.33 ± 1.12% (6XFFF beam with 2.5 mm CGS), 2.33 ± 1.37% (6XFFF beam with 1 mm CGS), 2.81 ± 2.33% (10XFFF beam with 2.5 mm CGS) and 3.34 ± 1.76% (10XFFF beam with 1 mm CGS) lower compared with that by AAA, respectively. However, the dose directed to PTV_soft was comparable. For OARs, AXB predicted a slightly lower dose to the aorta, chest wall, spinal cord and esophagus, regardless of whether the 6XFFF or 10XFFF beam was utilized. Exceptionally, dose to the ipsilateral lung was significantly higher with AXB.
AXB principally predicts lower dose to PTV_lung compared to AAA and the CGS contributes to the relative dose difference between the two algorithms.
The Anisotropic Analytical Algorithm (AAA), a convolution/superposition method, is widely utilized for dose calculation in Eclipse treatment planning system [1,2]. However, AAA was reported to significantly overestimate the dose near air-tissue interfaces . Recently, a new dose calculation algorithm named Acuros XB (AXB) has been implemented for clinical use. This new algorithm explicitly solves the linear Boltzmann transport equation that describes the macroscopic behavior of radiation particles as they travel through and interact with matters . AXB was found to provide good agreement with Monte Carlo (MC) algorithm [5-7] and achieved more accurate dose prediction than AAA in heterogeneous phantom [2,6-10].
The dose difference between the two algorithms has became an issue of concern in the past few years and several publications have investigated the dosimetric impact of AXB in breast cancer, non-small cell lung cancer (NSCLC) and nasopharyngeal carcinoma with AAA as a benchmark [4,11,12]. The dose difference between the two algorithms was reported to be closely related to the beam energy, field size and the density of the materials . However, they failed to consider the potential impact of the calculation grid size (CGS), which was reported to be associated with dose variation [13,14]. So far, the dose difference between AAA and AXB due to the use of different CGS in lung stereotactic body radiotherapy (SBRT) treatment remains unknown and the effect of it requires further investigation.
On the other hand, although two studies have evaluated the dosimetric impact of AXB in SBRT plan for lung cancer [15,16], little information has been focused on the dose difference generated from flattening filter free (FFF) beams, which was widely used in lung SBRT treatment for its improved treatment efficiency [17,18]. Because FFF beams provide softer spectrum than flattened beams , investigation of their impact is also essential.
Therefore, the aim of this study is to quantitatively determine the dose differences between AXB and AAA in SBRT treatment for lung cancer with FFF beams. Furthermore, the potential role of the CGS on the dose differences between the two algorithms was also discussed.
The same set of beam data (including percentage depth dose curve, profiles and output factors) used by AAA and measured in a three-dimensional water scanning system (PTW, Germany) for field sizes from 3 × 3 to 40 × 40 cm2 were imported in Eclipse treatment planning system (Version 10.0, Varian Medical Systems, Palo Alto, CA) for the configuration of AXB. All data presented in this study were collected from a commissioned Varian Truebeam™ accelerator equipped with a Millennium 120 multileaf collimator (MLC, with spatial resolution of 5 and 10 mm for the central and outer 20 cm, respectively.
CT scanning and contouring of organs at risk (OARs)
This retrospective study was approved by the ethics committee at the Cancer Hospital of Shantou University Medical College. Informed consent was obtained from all subjects prior to the study. During March 2012 and May 2013, 10 patients suffering from stage I NSCLC at our hospital were enrolled in this study. The CT datasets were acquired using a 16-slice CT scanner (The Philips Brilliance CT Big Bore Oncology Configuration, Cleveland, OH, USA) for all patients. Four dimensional computed tomography (4DCT) scanning and conventional enhanced 3D scanning with intravenous (IV) contrast was performed for 3 and 7 patients, respectively. Scanning was acquired at a 3 mm slice thickness for both 3D and 4DCT. CT images were then transferred to Eclipse treatment planning system. For patients with conventional enhanced scanning, gross tumor volume (GTV) was contoured by an experienced radiation oncologist under the CT pulmonary windows, and the planning target volume (PTV) was acquired according to the tumor motion under fluoroscopic examination with the aid of a simulator. For patients with 4DCT scanning, GTV accounting for tumor motion on all 10 phases of the 4DCT were contoured in the same way. These 10 phases of the GTV were then combined to form the internal target volume (ITV). To account for set-up uncertainties and potential baseline tumor shift, PTV was expanded with a uniform 5 mm margin from ITV. The PTV was then split into two substructures: PTV_soft, with a density in the soft-tissue range (0.590-0.985 g/cm3), and PTV_lung with a density in the lung range (0.011-0.590 g/cm3) . All plans were carried out on the enhanced CT. OARs, including the aorta, bronchial tree, chest wall, spinal cord, esophagus, heart, ipsilateral and contralateral lung were contoured according to the RTOG 0915 report .
The prescription and dose constraints for the target and OARs in this study followed RTOG 0915 criteria . Prescription was set to 4 fractions of 12 Gy. Treatment planning was designed with two volumetric modulated arc therapy (VMAT)-based partial arcs. Rotation angles were chosen to avoid the beams from entering the contralateral lung. Collimator angles for all plans were set to 30° in one arc and the complementary angle for another. The same beam arrangement and optimization constraints were employed for the 6XFFF and 10XFFF beams. Plans of different beam energy were optimized, selecting a maximum dose rate of 1400 MU/min and 2400 MU/min for 6XFFF and 10XFFF beams, respectively. Plan optimization was performed with the progressive resolution optimizer (PRO) algorithm implemented in Eclipse 10.0. Minimum and maximum dose constraints to the PTV were added during the planning optimization. Additionally, a dose-constraining structure (2 cm away from the PTV) was created to ensure a steep dose fall-off outside the PTV. The final dose calculation was normalized to guarantee that 95% of PTV received the prescribed dose. Dose calculations were carried out using AAA and AXB with a CGS of 2.5 mm and 1 mm, taking into account heterogeneity correction. For AXB calculation, we used the same monitor units (MU) and identical MLC settings as AAA algorithm and the “dose-to-medium” dose reporting mode was employed [12,21]. Accordingly, treatment plans of AAA with 2.5 mm and 1.0 mm CGS and AXB with 2.5 mm and 1.0 mm CGS were performed for each patient, respectively.
Evaluation of the plans and statistical analysis
For the PTV, the evaluation parameters included the minimum, maximum and the mean dose. For the OARs, the analysis included the maximum, mean dose and a set of appropriate Vx values. Vx was the volume of the organ receiving a dose of x or more. For example, V40 was the volume of organ receiving a dose of 40 Gy or more. The absolute doses in this study were presented in Gy and all the data were rounded to the nearest tenth. To determine the difference between plans, Wilcoxon signed-rank tests were performed using SPSS 17.0 (Chicago, IL). Differences were considered statistically significant when p < 0.05.
The dose difference between the AAA and AXB algorithm for SBRT plans of 10 lung cancer patients were summarized in Table 1 (6XFFF beam) and Table 2 (10XFFF beam). A considerable difference between the two algorithms was observed in PTV_lung, particularly when the 10XFFF beam plan with 1 mm CGS was used. A lower minimum dose and mean dose to PTV_lung were estimated by AXB than by AAA, regardless of whether 6XFFF or 10XFFF beam was used. Specifically, the dose directed to PTV_lung predicted by AXB was 1.33 ± 1.12% (6XFFF beam with 2.5 mm CGS) and 2.81 ± 2.33% (10XFFF beam with 2.5 mm CGS) lower compared with AAA. Furthermore, the dose difference between the two algorithms became larger by using 1 mm CGS (Figure 1). It was 2.33 ± 1.37% (6XFFF beam with 1 mm CGS) and 3.34 ± 1.76% (10XFFF beam with 1 mm CGS) lower compared with AAA. However, the minimum, maximum and mean dose directed to PTV_soft was comparable for the two algorithms, regardless of whether the 6XFFF or 10XFFF beam was used. When OARs were concerned, AXB always predicted a slightly lower dose than AAA. For example, AXB predicted a slightly lower dose to the aorta, chest wall, spinal cord and esophagus for the 6XFFF and 10XFFF beams. All of the dose change was statistically significant with p < 0.05. However, V5 and V10 of the ipsilateral lung were statistically higher with the AXB algorithm compared to the AAA. Figure 2 displayed the mean dose volume histograms (DVHs) of PTV and OARs at 6XFFF beam with 2.5 mm CGS, 6XFFF beam with 1 mm CGS, 10XFFF beam with 2.5 mm CGS and 10XFFF beam with 1 mm CGS, respectively. The dose distribution between the two algorithms from one representative case was illustrated in Figure 3.
The International Commission on Radiation Units (ICRU) recommended an overall accuracy of 5% in the delivery of absorbed dose , and the American Association of Physicists in Medicine (AAPM) advised that the uncertainty in the computed dose distribution should be less than 2% . Both of the AXB and AAA algorithms exhibit similar dose prediction in homogeneous media. However, AXB was reported to be more accurate in heterogeneous medium [2,5-7]. Accordingly, it is essential for us to accurately determine the dose delivered to lung tumor in which the water density target is always surrounded by low density lung tissue. For lung radiotherapy treatment, PTV in the low density range is originated either from the tumor motion uncertainty  or uncertainty in patient positioning and alignment of the therapeutic beams during the treatment . We determined that the mean dose of PTV_lung predicted by the AXB algorithm was up to 3.34% lower compared to the AAA by using the 10XFFF beam energy. Although the dose to PTV_lung was reduced, the PTV_soft dose was not affected. We acknowledged this as it was an interesting finding. For lung SBRT treatment, the PTV_lung dose changes were of positive significance because they could be associated with significant changes in tumor control probability (TCP) and normal tissue complication probabilities (NTCP) according to the AAPM Report No. 85 . Therefore, determination of the dosimetric effect with AXB on TCP will be an interesting work for our future study.
The dosimetric impact of the AXB algorithm has been investigated in several studies. Fogliata et al.  compared the AXB to AAA in patients with stage III NSCLC and concluded that a lower mean dose to soft tissue and a slightly higher mean dose in low-density lung tissue were obtained by the 6 MV beam with the AXB algorithm. Kan et al.  addressed this issue using IMRT and RapidArc techniques in nasopharyngeal carcinoma patients and found that AXB presented a 1% higher dose in air compared with AAA. The same mean target dose between the two algorithms was demonstrated by Liu et al.  in SBRT treatment for NSCLC patients. Kroon et al.  reported that D98% of PTV was lower using AXB, whereas the mean dose of PTV was comparable for the lung SBRT treatment plan. Our finding indicating that the mean dose of PTV_lung was obviously lower with the AXB algorithm differs from the results of these studies. These contradictions are reasonable because the differences in lung doses calculated by AAA and AXB were influenced by many factors, such as beam energy, field size, field number and the density of lung tissue [4,6]. Another possibility that may lead to the inconsistency was the use of FFF beam in our research. Unlike flattened beams, FFF ones can provide softer energy and a cone-like profile. Moreover, some differences in the electron contamination parameters were acquired by removal of the flattening filter, which was not observed in the flattened beams . These differences may influence the dose differences between the two algorithms and further investigation is necessary.
Additionally, we also determined that AXB estimated a lower dose for most of the OARs, whereas V5 and V10 of the ipsilateral lung were significantly higher with AXB. The higher dose prediction in ipsilateral lung could be explained by the shape of the profiles obtained by the AXB algorithm. Kan et al.  reported that the dose profile in air indicated a lower dose inside the field but a higher out-of-field dose predicted by AXB. When the gantry rotated around the patient, all OARs were inside the treatment field at some degrees of gantry rotation during the entire treatment process except for the ipsilateral lung. This helped to explain why V5 and V10 of the ipsilateral lung increased with the AXB algorithm. V5 was reported to be a dosimetric factor that was predictive of radiation-induced pneumonitis (RIP) , however, as the dose difference was minor (0.6-1.6% higher) with the AXB and its clinical effect could be neglected.
Our finding that the CGS also contributed to the dose difference between AAA and AXB is interesting. CGS was reported to be associated with dose estimation and calculation accuracy. Chung et al.  and Mittauer et al.  found that CGS influenced dose estimation in head and neck treatment. Ong et al.  found that 1.0 mm CGS improved the accuracy of dose calculation compared with 2.5 mm for the AAA algorithm. Kan et al.  also reported that AXB with 1.0 mm CGS resulted in a significant improvement in the dose accuracy (within 3%). Smaller grid resolution reduces the averaging effect and results in better sampling of the structure voxels to the calculation grid . Interestingly, our result demonstrated that CGS also influenced the relative dose difference between the AAA and AXB algorithm. By using 1 mm CGS, dose difference in PTV_lung between the two algorithms was greater compared with 2.5 mm CGS, regardless of whether the 6XFFF beam or 10XFFF beam was used. That’s due to the ability of AAA and AXB algorithms to adapt the CGS to ensure the dose accuracy in the axis perpendicular to the image slices. If the slice thickness is larger than the defined grid size, AAA and AXB may calculate the dose on dose planes between the image slices. And if the slice spacing is smaller than the defined grid size, AAA and AXB may skip calculating the dose on some slices. Based on the fact that the AXB algorithm predicts lower dose in each image slice, the dose difference between the two algorithms will undoubtedly become larger with 1 mm CGS because it involves more image slices in dose prediction. Compared with previous research, our result can more truly reflect the dosimetric difference between the two algorithms with 1 mm CGS due to its better calculation accuracy than 2.5 mm CGS. This was another factor that led to the inconsistency with other previous studies. To our knowledge, the potential impact of CGS on the dose difference between AAA and AXB has not been reported and analyzed and this study is the first to address this issue.
Our study has certain limitations. Although some recent publications have demonstrated the accuracy of AXB algorithm, few studies have verified the dosimetric performance of AXB algorithm in a SBRT plan produced from FFF beams. Our further work will focus on lung SBRT plans where small fields and air cavity are involved. Additionally, Kan et al.  recommended the use of 1 mm CGS instead of 2.5 mm whenever possible for stereotactic plans using a slice spacing of 1.25 mm, particularly when PTV contains relatively large areas of low density tissue. Whether different slice spacing will influence the relative dose difference between the two algorithms will be another area of interest for us.
The use of AXB mainly presents a lower dose in PTV_lung compared to the dose that is calculated from AAA in SBRT treatment for lung cancer, regardless of whether the 6XFFF or 10XFFF beam is used. The dose differences become larger by using higher beam quality and smaller CGS.
Anisotropic Analytical Algorithm
Stereotactic body radiotherapy
Calculation grid size
Flattening filter free
Planning target volume
Gross tumor volume
Internal target volume
Organs at risk
Non-small cell lung cancer
4 dimensional computed tomography
Radiation Therapy Oncology Group
Progressive resolution optimizer
Dose volume histogram
International Commission on Radiation Units & Measurements
American Association of Physicists in Medicine
Tumor control probability
Normal tissue complication probabilities
Ono K, Endo S, Tanaka K, Hoshi M, Hirokawa Y. Dosimetric verification of the anisotropicanalytical algorithm in lung equivalent heterogeneities with and without bone equivalent heterogeneities. Med Phys. 2010;37:4456–63.
Han T, Followill D, Mikell J, Repchak R, Molineu A, Howell R, et al. Dosimetric impact of Acuros XB deterministic radiation transport algorithm for heterogeneous dose calculation in lung cancer. Med Phys. 2013;40:051710.
Kan MW, Cheung JY, Leung LH, Lau BM, Yu PK. The accuracy of dose calculations by anisotropic analytical algorithms for stereotactic radiotherapy in nasopharyngeal carcinoma. Phys Med Biol. 2011;56:397–413.
Kan MW, Leung LH, Yu PK. Dosimetric impact of using the Acuros XB algorithm for intensity modulated radiation therapy and RapidArc planning in nasopharyngeal carcinomas. Int J Radiat Oncol Biol Phys. 2013;85:e73–80.
Bush K, Gagne IM, Zavgorodni S, Ansbacher W, Beckham W. Dosimetric validation ofAcuros XB with Monte Carlo methods for photon dose calculations. Med Phys. 2011;38:2208–21.
Fogliata A, Nicolini G, Clivio A, Vanetti E, Cozzi L. Dosimetric evaluation of Acuros XB Advanced Dose Calculation algorithm in heterogeneous media. Radiat Oncol. 2011;6:82.
Han T, Mikell JK, Salehpour M, Mourtada F. Dosimetric comparison of Acuros XB deterministic radiation transport method with Monte Carlo and model-based convolution methods in heterogeneous media. Med Phys. 2011;38:2651–64.
Han T, Mourtada F, Kisling K, Mikell J, Followill D, Howell R. Experimental validation of deterministic Acuros XB algorithm for IMRT and VMAT dose calculations with the Radiological Physics Center's head and neck phantom. Med Phys. 2012;39:2193–202.
Kan MW, Leung LH, Yu PK. Verification and dosimetric impact of Acuros XB algorithm on intensity modulated stereotactic radiotherapy for locally persistent nasopharyngeal carcinoma. Med Phys. 2012;39:4705–14.
Hoffmann L, Jorgensen MB, Muren LP, Petersen JB. Clinical validation of the Acuros XB photon dose calculation algorithm, a grid-based Boltzmann equation solver. Acta Oncol. 2012;51:376–85.
Fogliata A, Nicolini G, Clivio A, Vanetti E, Cozzi L. On the dosimetric impact of inhomogeneity management in the Acuros XB algorithm for breast treatment. Radiat Oncol. 2011;6:103.
Fogliata A, Nicolini G, Clivio A, Vanetti E, Cozzi L. Critical appraisal of Acuros XB and Anisotropic Analytic Algorithm dose calculation in advanced non-small-cell lung cancer treatments. Int J Radiat Oncol Biol Phys. 2012;83:1587–95.
Chung H, Jin H, Palta J, Suh TS, Kim S. Dose variations with varying calculation grid size in head and neck IMRT. Phys Med Biol. 2006;51:4841–56.
Mittauer K, Lu B, Yan G, Kahler D, Gopal A, Amdur R, et al. A study of IMRT planning parameters on planning efficiency, delivery efficiency, and plan quality. Med Phys. 2013;40:061704.
HW L i, Nugent Z, Clayton R, Dunscombe P, Lau H, Khan R. Clinical impact of using the deterministic patient dose calculation algorithm Acuros XB for lung stereotactic body radiation therapy. Acta Onco. 2013;9:1–6.
Kroon PS, Hol S, Essers M. Dosimetric accuracy and clinical quality of Acuros XB and AAA dose calculation algorithm for stereotactic and conventional lung volumetric modulated arc therapy plans. Radiat Oncol. 2013;8:149.
Ong CL, Verbakel WF, Dahele M, Cuijpers JP, Slotman BJ, Senan S. Fast arc delivery for stereotactic body radiotherapy of vertebral and lung tumors. Int J Radiat Oncol Biol Phys. 2012;83:e137–43.
Lang S, Shrestha B, Graydon S, Cavelaars F, Linsenmeier C, Hrbacek J, et al. Clinical application of flattening filter free beams for extracranial stereotactic radiotherapy. Radiother Oncol. 2013;106:255–9.
Vassiliev ON, Titt U, Ponisch F, Kry SF, Mohan R, Gillin MT. Dosimetric properties of photon beams from a flattening filter free clinical accelerator. Phys Med Biol. 2006;51:1907–17.
Radiation Therapy Oncology Group. A Randomized Phase II Study Comparing 2 Stereotactic Body Radiation Therapy (SBRT) Schedules for Medically Inoperable Patients with Stage I Peripheral Non-Small Cell Lung Cancer. 2014. [http://www.rtog.org/ClinicalTrials/ProtocolTable/StudyDetails.aspx?action=openFile&FileID=4673]. Accessed 12 Jan 2015.
Ma CM, Li J. Dose specification for radiation therapy: dose to water or dose to medium? Phys Med Biol. 2011;56:3073–89.
International Commission on Radiation Units & Measurements. Determination of absorbed dose in a patient irradiated by beams of X or gamma rays in radiotherapy procedures. ICRU Report 24; 1976.
American Asociation of Physicists in Medicine. Tissue inhomogeneity corrections for megavoltage photon beams. 2004. https://www.aapm.org/pubs/reports/rpt_85.pdf. Accessed 12 Jan 2015.
Ezhil M, Vedam S, Balter P, Choi B, Mirkovic D, Starkschall G, et al. Determination of patient-specific internal gross tumor volumes for lung cancer using four-dimensional computed tomography. Radiat Oncol. 2009;4:4.
International Commission on Radiation Units & Measurements. Prescribing, recording, and reporting photon-beam therapy. J ICRU. 2007;7:NP.
Fogliata A, Nicolini G, Clivio A, Vanetti E, Mancosu P, Cozzi L. Dosimetric validation of the Acuros XB Advanced Dose Calculation algorithm: fundamental characterization in water. Phys Med Biol. 2011;56:1879–904.
Baker R, Han G, Sarangkasiri S, DeMarco M, Turke C, Stevens CW, et al. Clinical and dosimetric predictors of radiation pneumonitis in a large series of patients treated with stereotactic body radiation therapy to the lung. Int J Radiat Oncol Biol Phys. 2013;85:190–5.
Ong CL, Cuijpers JP, Senan S, Slotman BJ, Verbakel WF. Impact of the calculation resolution of AAA for small fields and RapidArc treatment plans. Med Phys. 2011;38:4471–9.
This work was sponsored by Shantou University Medical College Clinical Research Enhancement Initiative No. 201424, and partly supported by research and travel grants from Science and Technology Key Project of Shantou City, China, 2012, No. 2012-165.
All authors certify that his manuscript has not been published in whole or in part nor is it being considered for publication elsewhere. The authors declare that they have no competing interests.
Each author had participated sufficiently in the work to take public responsibility for appropriate portions of the content. BH conceived, designed the experiments, collected the data and wrote the manuscript draft. LW configured the beam data of AXB algorithm; PL analyzed the data. CC revised the manuscript. All authors read and approved the final manuscript.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Huang, B., Wu, L., Lin, P. et al. Dose calculation of Acuros XB and Anisotropic Analytical Algorithm in lung stereotactic body radiotherapy treatment with flattening filter free beams and the potential role of calculation grid size. Radiat Oncol 10, 53 (2015). https://doi.org/10.1186/s13014-015-0357-0
- Dose difference
- Acuros XB
- Anisotropic analytical algorithm
- Stereotactic body radiotherapy
- Calculation grid size