Dosimetric Comparison of Integral Radiation Dose: Anisotropic Analytical Algorithm and Acuros XB in Breast Radiotherapy

The impact of the difference between Anisotropic Analytical Algorithm (AAA) and Acuros XB (AXB) in breast radiotherapy is not clearly due to different uses and further research is required to explain this effect. The aim of this study is to investigate the contribution of calculation differences between AAA and AXB to the integral radiation dose (ID) on critical organs. Seven field intensity modulated radiotherapy (IMRT) plans were generated using with AAA and AXB algorithms for twenty patients with early stage left breast cancer after breast conserving surgery. Volumetric and dosimetric differences, as well as, the Dmean, V5, V20 doses of the left and right-sided lung, the Dmean, V10, V20, V30 doses of heart and the Dmean, V5, V10 doses of the contralateral breast were investigated. The mean dose (Dmean), V5, V20 doses of the left-sided lung, the Dmean, V5, V10 doses of right-sided lung, the Dmean, V10, V20, V30 doses of heart and the Dmean, V5, V10 doses of the contralateral breast were found to be significantly higher with AAA. In this research integral dose was also higher in the AAA recalculated plan and the AXB plan with the average dose as follows left lung 2%, heart 2%, contralateral breast 8%, contralateral lung 4% respectively. Our study revealed that the calculation differences between Acuros XB (AXB) and Anisotropic Analytical Algorithm (AAA) in breast radiotherapy caused serious differences on the stored integral doses on critical organs. In addition, AXB plans showed significantly dosimetric improvements in multiple dosimetric parameters.


Introduction
Radiotherapy contributes to reducing the risk of local recurrence after surgery in the treatment of early stage breast cancer. Because of the proximity of surrounding tissues, the chest region has a very heterogeneous structure. Because of this region that contains heterogeneous tissues such as the treatment area, skin, lung, heart, and also anatomical difficulties; breast radiotherapy is difficult to apply. Depending on the technique used during breast radiotherapy, healthy tissues can take wide place in the treatment area; thus healthy tissues in and around the treatment area can be exposed to high integral doses.
The integral dose (ID) is the volume integral of the dose stored in a medium and it is equal to multiplication of the mean dose and volume in the medium in which the radiation is applied at any doses. The integral dose is also the area under the curve of a differential dose absolute-volume histogram. It is often stated that a large number of beam and monitor units (MU) used in IMRT cause an increase in ID and that the photon rays with high energy principally reduce the ID. D'Souza et al. reported that the change in ID with four or more beams is <1% as a function of the number of beams. As expected; high energy rays reduce the ID; they showed that the ID value was 1.5% -1.7% for nasopharynx, 0.9% -1.0% and 0.3% for the pancreas, and 0.4% for the prostate. These results showed that ID reduces with increasing tumor size for similar anatomic dimensions and the ID increases with increasing size of anatomical region for similar tumor sizes [1].
The human body is a structure that shows different density due to its chemical structural elements. Hounsfield Unit (HU) mean value is −1000 HU in air, +1000 HU in bone, −50 -100 HU in fat tissue and −500 HU in lungs. Classically, the Hounsfield Unit values obtained from computerized tomography are calculated with the decrease of radiation in tissue using tables specific to predefined density ranges, and CT data, CT calibration curves. Algorithms and dose calculation mechanisms that may include tissue composition may have different accuracy in determining the dose in each organ.
In the Eclipse 13.0 treatment planning system (Varian Medical Systems, Palo Alto, CA), the Anisotropic Analytical Algorithm (AAA) method is commonly used for dose calculation [2] [3] [4] [5]. In the literature, there are studies reporting that the calculated doses of AAA were significantly incorrectly calculated. Near the two mediums, especially in the transition from tissue to air, it was seen that the dose was incorrectly calculated [6]. Recently, a new dose calculation algorithm, Acuros XB (Varian Medical Systems, Palo Alto, CA), is started to be used in clinic to correct this situation. AXB uses a complex technique to solve Linear Boltzmann Transport Equation (LBTE) and provides an accurate approach to the calculation of patient dose in heterogeneities originating entirely from lung, bone, air and different density implants [7] [8] [9] [10].
LBTE defines the macroscopic behavior of the radiation beam in the medium through which they pass.

A. Cakir, Z. Akgun
Up to now, the effect of the difference between AAA and AXB on breast radiotherapy is not known due to different uses and further research is required to explain this effect. The difference in calculation between the two algorithms, which is affected by parameters such as the energy of the incoming beam, the field size and the electron density of the medium, reveals dose responses.
The aim of this study is to investigate the effect of the calculation differences between the AXB and AAA algorithms on the integral doses of critical organs in the radiotherapy of patients who underwent breast-conserving surgery due to breast cancer.

Analytical Anisotropic Algorithm
The AAA dose calculation model is a 3D pencil beam and convolution superposition algorithm consisting of separate models for electrons emitted from primary photons, scattered photons, and beam modulators (primary collimator, beam straightening filter, and wedge filter). The functional forms for the basic physical quantities initiate a process by adding device properties to the account.
This often leads to a noticeable reduction in the calculation time required for such algorithms. Tissue heterogeneities are added to calculation as anisotropic using photon scattering kernels of the multiple lateral direction in the 3D neighborhood. The final dose distribution is formed by superimposing the contribution of photons and electron beams. The dose behind the air gap is a little over-calculated in the AAA algorithm with an error caused by the modeling of the scattered dose.

Acuros XB Algorithm
The Acuros XB (AXB) algorithm was developed for two strategic needs such as accuracy and speed in external photon beam treatment planning. AXB uses a complex technique to solve Linear Boltzmann Transport Equation (LBTE) and clarifies the patient dose calculation for heterogeneities originating entirely from lung, bone, air, and non-biological implants. LBTE, which is the linear form of the Boltzman Transport Equaiton (BTE) and defines the macroscopic behavior of the radiation particles, assumes that the interaction of the radiation particles in the medium is formed without interaction of each other and an external magnetic field. There are two solution approaches that try to explain LBTE. One of them is the Monte Carlo method, which does not clearly solve the commonly known LBTE and produces indirect solutions for LBTE. Second one is solving LBTE using numerical methods. The AXB algorithm uses the mass density information obtained from the CT images in each voxel for dose calculation.

Phantom Study
In order to measure the dosimetric difference between Analytical Anisotropic Algorithm and Acuros XB, measuring assembly was prepared as in Figure 1

Treatment Planning of Patients
Planning CT scans of 20 patients with early-stage left breast cancer, after breast-conserving surgery, were used for this study (Table 1)

Evaluation of Dose Volume Histograms
Dose volum histograms (DVH) were evaluated for RTOG target volumes and OARs of each patient plan. While the treatment of the patients was being planned, the primary dose limitation was 95.0% of the target and 95.0% of the dose prescribed for chest or chest wall was provided to be received as PTV breast .
From each patient's treatment planning; the V 5 Gy, V 20 Gy, mean (D mean ) doses were calculated for the ipsilateral lung, V 10 Gy, V 20 , V 30 Gy and D mean for the heart, V 5 , V 10 and D mean for the contralateral breast, V 5 , V 10 and D mean doses for the contralateral lung. Besides all these; Body-PTV doses which are not in the RTOG criteria were calculated. V 5 and V 20 doses were calculated separately for both algorithms, and small dose-volume responses of algorithms were investigated.

Evaluation of Phantom Measurements
In the measuring assembly shown in Figure 1 which we prepared to demonstrate the differences in calculation between Analytical Anisotropic Algorithm and Acuros XB, the difference in dose between AAA and TPS was found to be 4% in A, which is the first dose measurement point after 5 cm air cavity, 7.8% in B, which is the second dose measurement point, and 9.5% in C. Similarly; with AXB at the same measuring points, we achieved a difference of 1% at point A, 1.8% at point B and 3% at point C.

Evaluation of Dose-Volume Histograms
The dosimetric results for the treatment plan of 20 breast cancer patients formed using two different algorithms, are shown in Table 2. PTV breast minimum, PTV breast D mean , PTV breast D 95% doses; V 5 Gy, V 20 Gy and D mean doses for the contralateral lung; V 10 Gy, V 20 , V 30 Gy and D mean doses for the heart; V 5 , V 10 and D mean doses for the contralateral breast; V 5 , V 10 and D mean doses for the contralateral lung are shown in Table 2.
In Table 2; V 5 and V 20 doses are shown separately for both algorithms by calculating the Body-PTV doses, which are not in the RTOG criteria.
Dose-Volume Histogram (DVH) for a patient calculated using the AAA and AXB dose algorithms is shown in Figure 2. There was a difference of 13% between the PTV breast Dmin doses calculated by the AAA and AXB algorithms (p < 0.05). This may be related to the fact that the PTV breast volume was calculated as high by the AAA algorithm although the PTV breast volume was removed by 0.5 cm from the skin. There was no significant difference between the PTV D95% dosescalculated by the AAA and AXB algorithms (p > 0.05).
There was a significant difference between the V 5 , V 20 and D mean doses calculated by AXB and AAA algorithms for the left side lung (p < 0.05). V 5 , V 20 and D mean doses of the left side lung are higher in all AAA plans than the dosimetric results of the AXB plans.
There was a significant difference between the V 10 , V 20 , V 30 and D mean doses of the heart in all the plans of the AAA and AXB algorithms (p < 0.05). Figure 3 shows comparison of dose distribution calculated using AAA and AXB dose calculation algorithms for a patient early-stage left breast cancer. The same 7 field IMRT plan recalculated using the AXB dose calculation algorithm (for the same patient, with the same dose scale).
A. Cakir, Z. Akgun Table 2. PTV breast D min : The minimum dose of PTV breast volume, PTV breast D mean : The mean dose of PTV breast volume, PTV breast D 95% : The dose prescribed to the 95% isodose line for the PTV breast , V 5 : % Volume receiving dose of 5 Gy, V 10 : % Volume receiving dose of 10 Gy, V 20 : % Volume receiving dose of 20 Gy, V 30 : % Volume receiving dose of 30 Gy, Body-PTV: The volume remaining after the removal of the PTV volume from the total body volume.   When D mean doses of the contralateral lung doses were evaluated, there was a significant difference between AAA plans and AXB plans (p < 0.05). There was also a significant difference between the doses of V 5 and V 10 in the contralateral lung (p < 0.05).
There was a significant difference between the D mean , V 5 and V 10 doses of the contralateral breast (p < 0.05).
Finally; in our study comparing the V 5 and V 20 doses for Body-PTV, it was found that there was no significant relationship between AAA plans and AXB planning [11]. To reduce this effect in our study, all plans for both algorithms were calculated at 2.5 mm CGS.
In a planning study performed for nasopharyngeal cancers by Kan et al., they reported that a lower dose between 3% and 6% was obtained with AXB on critical organs compared to AAA [12]. Chung [15]. In our dosimetric phantom study, we detected that there was a 3% compliance was between the depth dose measurement and Acuros AXB algorithm. However, we found there was a dose difference of 9.5% between the depth dose measurement and AAA algorithm.
In the phantom study with 6 MV photon rays by Bush et al., they showed that there was difference of 4.5% between AXB and Monte Carlo algorithms in the transition from air to tissue and this difference increased to 13% between AAA and Monte Carlo algorithms [16]. In parallel with this study, Kan et al. found A. Cakir, Z. Akgun that there was a difference of 3% between the measurement and AXB algorithm in the transition from air to tissue and this difference reached 10% with AAA algorithm (12). In the rando phantom study study developed for head and neck treatments by Han et al., they found the maximum difference between measurement and calculation as 4.6% for AAA and 3.6% for AXB [17].
In a study comparing the treatment plans of deep inspriation condition and free-breathing breast patients by Fogliata

Conclusions
The most remarkable side of our study is that the dose of breast tissue that starts immediately after the skin was calculated as higher with the AAA algorithm.
There was difference of 13% between AAA plans and AXB plans in determining PTV min dose. The AAA algorithm calculated a higher dose than should be in the build-up area between air and tissue. The higher dose should be in PTV, and will increase the maximum dose effect in the hot dose regions as a result of the normalization of the plan to the treatment dose.
The integral dose (ID) is the volume integral of the dose stored in a medium and it is equal to multiplication of the mean dose (D mean ) and volume in the medium. Often, the large number of beams and high irradiation times used in the IMRT cause an increase in the ID and this increase in the ID is prevented by using high energy photon rays. In the literature, there is no study comparing integral doses of healthy organs between AAA and AXB algorithms.
Our study revealed that the calculation differences between Acuros XB (AXB) and Anisotropic Analytical Algorithm (AAA) in breast radiotherapy caused serious differences on the stored integral doses on critical organs.
As a result; in the region where we are pushing critical dose limits such as skin, lung, heart, and contralateral breast, the administering of breast radiotherapy and the accuracy of dose delivery are very difficult due to both anatomical and different density tissues. In addition to very different tissue densities within the treatment area, many devices that increase the dosimetric uncertainty associated with the patient stabilizing devices found that the outside treatment area affects the dose in the patient. It should be kept in mind that treatment planning algorithms do not yet have the ability to accurately calculate the dose in air-tissue transition zones.

Conflicts of Interest
None.