Introduction {#S0001}
============

Stereotactic body radiation therapy (SBRT) has become a very important treatment technique for early stage non-small cell lung cancer,[@CIT0001],[@CIT0002] which delivers high radiation doses to a tumor target in a hypo-fractionated schedule with the minimum dose to the tissue around the target area. Therefore, it can protect the normal tissues and organs around the target. However, for non-homogeneous tissues such as lung, the usage of SBRT technology was challenged with the problem of accurate dose calculation.[@CIT0003] Therefore, SBRT of non-small cell lung cancer needs a more precise algorithm, which takes into account tissue non-uniformity correction to ensure the accuracy of dose calculation. Varian Eclipse Treatment Planning System (Varian Medical Systems, Palo Alto, USA) is a widely used treatment planning system in the current radiotherapy physics program with different versions. The latest Eclipse Treatment Planning System (Version 13.5) is equipped with Acuros XB (AXB) algorithm. However, a more precise AXB algorithm is not widely used due to many factors.[@CIT0004]--[@CIT0010] Most of the studies only focus on the differences of algorithms in motifs. In order to intuitively understand the difference between different algorithms in the practice of patient planning, here a systematic comparison between AXB, anisotropic analytic algorithm (AAA), and pencil beam convolution (PBC) on dosimetric differences of SBRT in non-small cell lung cancer was studied with the aim to provide references for the treatment planning of clinical radiotherapy.

Patients and methods {#S0002}
====================

Patients selection and target volume delineation {#S0002-S2001}
------------------------------------------------

The study was approved by the Ethics Committee of Zhongnan Hospital of Wuhan University. All patients signed informed consent prior to SBRT treatment, which included allowing their data to be used for further research. Owing to the retrospective nature of the present study, it was granted an exemption in writing by the Ethics Committee of Zhongnan Hospital of Wuhan. A total of 38 patients with non-small cell lung cancer who received SBRT from February 2014 to June 2018 in Zhongnan Hospital of Wuhan University were included. Thirty-eight patients consisted of 27 males and 11 females, with a median age of 65 years (range, 30--81 years). Thirty-five patients were stage IVB, and three patients were stage IA (According to the AJCC 2017 staging). These patients with stage IV had already received systemic chemotherapy or targeted therapy before receiving the treatment of SBRT. All of them experienced a progress of the disease or cannot continue to afford the other systemic therapy such as the third- generation TKI or immunotherapy due to the bad social economy status. The location of tumor target was as follows: 16 cases in the right upper lung, 4 cases in the right lower lung, 13 cases in the left upper lung, and 5 cases in left lower lung. The mean GTV volume is 24.08 cm^3^ (range, 0.5--114.1 cm^3^). Patient simulation was positioned using the whole-body immobilization system (R624-SCF, Klarity, Guangzhou, China). CT simulation images with 3 mm thickness were gained with a large aperture simulation CT scanner (Siemens Somatom Sensation, Munich, German) for each patients. Tumor motion was assessed with four-dimensional CT imaging with Varian real-time position management (Varian Medical Systems, Palo Alto, CA). The abdominal compression was employed for all patients. CT images were transferred to Eclipse version 13.5 planning system for target volume delineation and plan design. Gross tumor volume and internal target volume were segmented on a four-dimensional CT scan, respectively. A 5-mm uniform margin from ITV was typically used for creation of the PTV.[@CIT0005] Organs at risk (OARs) were contoured including the heart, the spinal cord, the chest wall, the lung and so on. The same attending radiation oncologists delineated the target volumes and OARs. SBRT plans were generated by the same senior medical physicist.

SBRT planning {#S0002-S2002}
-------------

Varian Eclipse TPS (Version 13.5) was introduced to calculate all plans. Linear accelerator Clinac IX (Varian Medical Systems, Palo Alto, USA) was used for the radiotherapy treatment. The prescription dose of 38 patients for all tumors was 50 Gy in 5 fraction at the 90% isodose line of the maximum dose delivered every other day. All the dose limits to the critical organs are according to the protocol guidelines of the two trials (RTOG0813 and RTOG0618) and the recommendations of the American Association of Physicists in Medicine Task Group 101 report. The recommended maximal dose to the trachea and proximal bronchial tree was 30, 35, and 40 Gy for 3-, 4-, and 5-fraction regimens, respectively. The recommended esophageal maximal dose ranged from 25.2 to 35 Gy, and maximal heart dose ranged from 30 to 38 Gy, increasing with the number of fractions used. For 5-fraction regimens, mediastinal structures were permitted to have a maximal point dose of 105% of the prescription dose. All plans were originally calculated with AAA algorithm, the computational grid was 1.5 mm ⨰ 1.5 mm, and all plans are corrected for non-uniformity, then with the AXB and PBC algorithms with the same prescription dose and field parameters. Each plan met with the requirements of the RTOG 0813 Report.[@CIT0005] One hundred percent of the prescription dose includes at least 95% of PTV, 90% of the prescription dose included at least 99% of PTV.

Evaluation of dosimetry {#S0002-S2003}
-----------------------

In this study, the target coverage (V100% and V90% of the PTV), Conformal index (CI), Homogeneity index (HI), minimum, maximum, and average dose of PTV (D~min~, D~max~, D~mean~) were compared. V100% was defined as the volume of the prescription dose line for 100% includes the volume of PTV, V90% was defined as the volume of PTV contained in the prescription dose line for 90%, and the volume was normalized to the volume of the target area PTV. CI was calculated with Paddick Conformal index：$\documentclass[12pt]{minimal}
\usepackage{wasysym}
\usepackage[substack]{amsmath}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage[mathscr]{eucal}
\usepackage{mathrsfs}
\DeclareFontFamily{T1}{linotext}{}
\DeclareFontShape{T1}{linotext}{m}{n} {linotext }{}
\DeclareSymbolFont{linotext}{T1}{linotext}{m}{n}
\DeclareSymbolFontAlphabet{\mathLINOTEXT}{linotext}
\begin{document}
$${\rm{C}}{{\rm{I}}_{{\rm{paddick}}}}{\rm{ = }}{{{{{\rm{(T}}{{\rm{V}}_{{\rm{PI}}}}{\rm{)}}}^{\rm{2}}}} \over {{\rm{PI*TV}}}}$$
\end{document}$, PI was the volume included in the prescription dose line, TV was the volume of PTV, TV~PI~ was the volume of the prescription dose line including the volume of the PTV. A perfect plan was TV=PI=TV~PI~ and CI~Paddick~=1. HI was calculated according to the ICRU 83 Report：$\documentclass[12pt]{minimal}
\usepackage{wasysym}
\usepackage[substack]{amsmath}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage[mathscr]{eucal}
\usepackage{mathrsfs}
\DeclareFontFamily{T1}{linotext}{}
\DeclareFontShape{T1}{linotext}{m}{n} {linotext }{}
\DeclareSymbolFont{linotext}{T1}{linotext}{m}{n}
\DeclareSymbolFontAlphabet{\mathLINOTEXT}{linotext}
\begin{document}
$${\rm{HI = }}{{{{\rm{D}}_{{\rm{2\% }}}}{\rm{ - }}{{\rm{D}}_{{\rm{98\% }}}}} \over {{{\rm{D}}_{{\rm{50\% }}}}}}$$
\end{document}$, D~2%~ was the maximum dose of volume accepted by 2% of PTV, D~98%~ was the maximum dose of volume accepted by 98% of PTV, D~50%~ was the maximum dose of volume accepted by 50% of PTV. The value of the HI indicates the uniform dose distribution. When HI=0 indicates that there is no drop in the target area, the uniformity is the best.

Statistical analysis {#S0002-S2004}
--------------------

Data were analyzed using the Statistical Package for Social Sciences version 22.0 (IBM Corporation, Armonk, NY, USA). The significance of differences between mean values was determined by one-way ANOVA (Bonferroni test). *P*-values \<0.05 were considered statistically significant.

Results {#S0003}
=======

Clinical characteristics {#S0003-S2001}
------------------------

Patient and tumor characteristics of the 38 patients are described in [Table 1](#T0001){ref-type="table"}.Table 1Patient and tumor characteristicsParameterTotalPatients (n=38)Female =11, male =27Median age in years (range)65 (30--81)Median GTV in cm^3^ (range)24.08 (0.5--114.1)Tumor location (n=38)13 LUL, 16 RUL, 5 LLL, 4 RLL[^2]

V100% and V90% of the PTV {#S0003-S2002}
-------------------------

For 38 patients, the average values and standard viations are summarized in [Table 2](#T0002){ref-type="table"} and [Figure 1](#F0001){ref-type="fig"} for V100% and V90% of the PTV. AXB plan resulted in decreased V100% with a mean difference 6.14% compared with PBC plan (For V100%, AXB vs AAA vs PBC=93.44% vs 95.54 vs 99.58%, *P*\<0.05). Three plans showed no significant difference as to the parameter V90% (For V90%, AXB vs AAA vs PBC=98.66% vs 99.75% vs 99.83%, *P*\>0.05).Table 2The V100% and V90% of PTV of the three algorithmsParameterMean ± SD (%)*P*-valueAXBAAAPBCAXB vs AAAAXB vs PBCAAAvs PBCV100%93.44±2.0395.54±0.4799.58±1.010.210.000.02V90%99.66±0.7399.75±0.6399.83±0.490.320.260.99[^3][^4] Figure 1The comparison of the V100% (**A**) and V90% (**B**) of the PTV in the three algorithms. \*means *P*\<0.05.**Notes:** V100% was defined as the volume of the prescription dose line for 100% includes the volume of PTV. V90% was defined as the volume of PTV contained in the prescription dose line for 90%.**Abbreviation:** SD, standard deviation.

D~min~, D~max~, and D~mean~ of PTV {#S0003-S2003}
----------------------------------

D~min~, D~max~, and D~mean~ of PTV in 38 patients are indicated in [Table 3](#T0003){ref-type="table"} and [Figure 2](#F0002){ref-type="fig"}. AXB plan leaded to reduced D~min~ of PTV compared with other two algorithms (For D~min~ of PTV, AXB vs AAA vs PBC=4048cGy vs 4365Gy vs 4873Gy, *P*\<0.05). PBC induced the enhanced trend of D~max~ of PTV compared with other two algorithms (For D~max~ of PTV, AXB vs AAA vs PBC=5897cGy vs 5904Gy vs 6122Gy; AXB vs AAA *P*=0.41; AXB vs PBC *P*=0.06; AAA vs PBC: *P*=0.01), and increased the D~mean~ of PTV in three algorithms with significant difference (For D~mean~ of PTV, AXB vs AAA vs PBC=5332cGy vs 5330Gy vs 5785Gy, *P*\<0.05).Table 3D~min~, D~max~, and D~mean~ of PTVParameterMean ± SD (cGy)*P-*valueAXBAAAPBCAXB vs AAAAXB vs PBCAAA vs PBCD~min~4048±3274365±2314873±3060.040.000.00D~max~5897±2455804±2406122±1750.410.060.01D~mean~5332±1165330±1325785±1660.970.000.00[^5][^6] Figure 2The comparison of D~min~(**A**), D~max~(**B**), and D~mean~ (**C**) of PTV in the three algorithms. \*means *P*\<0.05.

CI and HI of the PTV {#S0003-S2004}
--------------------

CI and HI of the PTV in 38 patients are shown in [Table 4](#T0004){ref-type="table"} and [Figure 3](#F0003){ref-type="fig"}. AXB algorithm achieved a similar plan conformity with other two algorithms (For CI, AXB vs AAA vs PBC=0.80 vs 0.85 vs 0.71, *P*\>0.05) but slightly higher tend of target dose heterogeneity (For HI, AXB vs AAA vs PBC=0.20 vs 0.14 vs 0.31, *P*\>0.05).Table 4CI and HI of PTVParameterMean±SD*P*-valueAXBAAAPBCAXB vs AAAAXB vs PBCAAA vs PBCCI0.80±0.080.85±0.060.71±0.050.140.010.00HI0.20±0.070.14±0.080.31±0.070.130.080.02[^7][^8] Figure 3The comparison of CI (**A**) and HI (**B**) of the PTV in the three algorithms. \*means *P*\<0.05.**Abbreviations:** SD, standard deviation; CI, conformal index; HI, Homogeneity index.

Dosimetric distribution for the plan target volume {#S0003-S2005}
--------------------------------------------------

The isodose distributions in the axial view at the isocenter and dose volume histogram for the same patient calculated from the three algorithms are presented in [Figure 4](#F0004){ref-type="fig"}. The PTV coverage of AXB plan was lower than that of the other two plans. It can be seen from [Figure 4](#F0004){ref-type="fig"} that the 100% prescription dose line of the PBC algorithm includes the largest target area.Figure 4The isodose distributions in the axial view at the isocenter and dose volume histogram (DVH) of a patient. The red line was PTV, the green line is the isodose line of prescription dose 5000cGy (100%), and the light blue line is an isodose line of 2500cGy (50%).

Discussion {#S0004}
==========

The application of different algorithms in the radiotherapy physics plan has always been the research direction of scholars.[@CIT0003],[@CIT0004],[@CIT0007]--[@CIT0013] The AXB algorithm is similar to the Monte Carlo algorithm, the physical action of the ray in the medium was clearly simulated, and the linear Boltzmann linear transmission problem was solved in the discrete and multidimensional space.[@CIT0006] Many scholars have confirmed that the AXB algorithm was the closest to the Monte Carlo algorithm,[@CIT0010],[@CIT0014] the accuracy was higher than the AAA and PBC algorithms.[@CIT0015]--[@CIT0017] The AAA algorithm is based on the superposition convolution technique, and the uniformity correction is considered in the calculation.[@CIT0009] The PBC algorithm is a pencil beam algorithm, which calculates dose distribution through a very narrow beam, each pencil takes into account the attenuation factor of different beams, even the method of equivalent tissue air ratio is used for non-uniformity correction. But the lateral non-uniformity correction is not considered in the PBC algorithm.[@CIT0007]

There is a strong correlation between tumor control probability and radiation dose in SBRT radiotherapy, accurate radiation dose to the target volume is very important.[@CIT0008] Therefore, it is necessary to adopt a more accurate dose calculation algorithm in radiotherapy physics plan. The present research showed that V100% in group of the PBC algorithm was more than that in the group of the AXB and AAA algorithms, which indicated that the PBC algorithm significantly overestimates the dose of PTV in dose calculation. Meanwhile, the V100% in the group of AAA algorithm was 3% larger than that in the group of AXB without statistical significance (*P*=0.21). The results here were consistent with the previous results of Rana.[@CIT0013] Rana and his colleagues found that the V100% in the group of AAA is 8.2% higher than that in the group of AXB algorithm in 16 cases of non-small cell lung cancer. It showed that the AAA algorithm would have an effect on target dose. For V90%, the difference between AAA and PBC algorithm was not statistically significant compared with AXB algorithm. The value of V90% in three algorithms was more than 98%, which showed that 90% of the prescription dose lines obtained by the three algorithms can effectively encircle the target area.

The present study also found that the trend of the D~min~, the D~max~ and the D~mean~ for PTV was the same as that of V100%. Compared with AXB, 20% increase of D~min~ in PBC algorithm (*P*=0.00) was found, 4% increase for D~max~ (*P*=0.06), about 8% increase for D~mean~ (*P*=0.00). All of these data also indicated that the PBC algorithm significantly could overestimate the target dose. It was also found that the D~min~ calculated by the AAA algorithm was about 8% (*P*=0.04) more that gained by the AXB algorithm, but there was no significant difference between the D~max~ and the D~mean~ between the two algorithms. These implied that the AAA algorithm does not significantly overestimate the target dose.

The CI for PTV in the group of the AAA was about 6% more than that in the AXB algorithm, and in the PBC algorithm was the smallest. The probable reason was overestimating the dose of the target area by the PBC algorithm which results in a larger volume size of the prescription dose. HI of PTV gained from PBC algorithm was maximum, which showed that the dose of the PTV of the PBC algorithm was the most heterogeneous. The foregoing results indicated that the PBC algorithm can significantly overestimate the target dose because the lung is considered as a homogeneous tissue, the AAA algorithm would also overestimate the dose of PTV. AXB can more authentically simulate the physical process of radiation in heterogeneous tissue.

Conclusion {#S0005}
==========

The present study demonstrated that the PBC and the AAA algorithm could overestimate the dose of PTV for SBRT physics plan for non-small cell lung cancer compared to the AXB algorithm. Therefore, the AAA and PBC algorithms were inaccurate for the PTV calculation of the planned target area. In the clinical plan of SBRT for non-small cell lung cancer, the AXB algorithm should be suggested to improve the accuracy of the dose plan.

This research was partially supported by Zhongnan Hospital of Wuhan University Science, Technology and Innovation Seed Fund (Grant No. ZNPY2016043), Natural Science Foundation of Hubei Province (Grant No. 2018CFB561) and Scientific Fund of Health Commission of Hubei Province (Grant No. WJ2019M264).

Ethic statement {#S0006}
===============

This study has been approved by the Ethics Committee of Zhongnan Hospital of Wuhan University.

Disclosure {#S0007}
==========

The authors declare no conflicts of interest in this work.

[^1]: These authors contributed equally to this work

[^2]: **Abbreviations:** LUL, left upper lobe; RUL, right upper lobe; LLL, left lower lobe; RLL, right lower Lobe.

[^3]: **Note:** *P\<*0.05 considered statistically significant.

[^4]: **Abbreviation:** SD, standard deviation.

[^5]: **Note:** *P\<*0.05 considered statistically significant.

[^6]: **Abbreviation:** SD, standard deviation.

[^7]: **Note:** *P\<*0.05 considered statistically significant.

[^8]: **Abbreviation:** SD, standard deviation.
