Dosimetric accuracy of delivering SBRT using dynamic arcs on Cyberknife

Purpose Several studies have demonstrated potential improvements in treatment time through the use of dynamic arcs for delivery of stereotactic body radiation therapy (SBRT) on Cyberknife. However, the delivery system has a finite accuracy, so that potential exists for dosimetric uncertainties. This study estimates the expected dosimetric accuracy of dynamic delivery of SBRT, based on realistic estimates of the uncertainties in delivery parameters. Methods Five SBRT patient cases (prostate A — conventional, prostate B — brachytherapy‐type, lung, liver, partial left breast) were retrospectively studied. Treatment plans were produced for a fixed arc trajectory using fluence optimization, segmentation, and direct aperture optimization. Dose rate uncertainty was modeled as a smoothly varying random fluctuation of ± 1.0%, ±2.0% or ± 5.0% over a time period of 10, 30 or 60 s. Multileaf collimator uncertainty was modeled as a lag in position of each leaf up to 0.25 or 0.5 mm. Robot pointing error was modeled as a shift of the target location, with the direction of the shift chosen as a random angle with respect to the multileaf collimator and with a random magnitude in the range 0.0–1.0 mm at the delivery nodes and with an additional random magnitude of 0.5–1.0 mm in between the delivery nodes. The impact of the errors was investigated using dose‐volume histograms. Results Uncertainty in dose rate has the effect of varying the total monitor units delivered, which in turn produces a variation in mean dose to the planning target volume. The random sampling of dose rate error produces a distribution of mean doses with a standard deviation proportional to the magnitude of the dose rate uncertainty. A lag in multileaf collimator position of 0.25 or 0.5 mm produces a small impact on the delivered dose. In general, an increase in the PTV mean dose of around 1% is observed. An error in robot pointing of the order of 1 mm produces a small increase in dose inhomogeneity to the planning target volume, sometimes accompanied by an increase in mean dose by around 1%. Conclusions Based upon the limited data available on the dose rate stability and geometric accuracy of the Cyberknife system, this study estimates that dynamic arc delivery can be accomplished with sufficient accuracy for clinical application. Dose rate variation produces a change in dose to the planning target volume according to the perturbation of total monitor units delivered, while multileaf collimator lag and robot pointing error typically increase the mean dose to the planning target volume by up to 1%.


INTRODUCTION
The Cyberknife system (Accuray Inc., Sunnyvale, CA) has shown itself to be a valuable device for treating patients with stereotactic body radiotherapy (SBRT). [1][2][3][4] The short-waveguide 6-MV flattening-filter-free linear accelerator is mounted on a robotic arm and is equipped with either a series of circular collimators, a variable circular diaphragm, or a multileaf collimator (MLC). [5][6][7] The MLC is widely used to allow treatment of larger tumors using fewer monitor units.
The Cyberknife currently delivers radiation from a number of static locations around the patient in a step-and-shoot arrangement. 8 However, a number of studies have demonstrated the potential for reduction in delivery time by the use of dynamic arc delivery, similar in nature to volumetric modulated arc therapy (VMAT), although from noncoplanar orientations as opposed to the more common coplanar arcs used for VMAT. For example, Kearney et al. 9 describe a noncoplanar arc optimization algorithm for Cyberknife with a circular collimator. They also describe an optimization method for producing dynamic arcs on the Cyberknife with MLC, using direct aperture optimization. 10 Simultaneously, data are beginning to emerge on the accuracy of the Cyberknife system. Yang et al. 11,12 describe a calibration procedure for the ArcCheck quality assurance device (Sun Nuclear, Melbourne, FL) that converts between a spatial error and a dosimetric error, so that the device can be used to measure the pointing accuracy of the robot. Wang et al. 13,14 use a scintillator and charge-coupled device camera to record a pair of spots located on the beam axis and calculate from these spots the position of the beam to an accuracy in the order of 0.1 mm. With these data, it is possible to estimate the geometric accuracy of a dynamic delivery, with a view to determining the dosimetric performance.
This study therefore investigates the performance of arc delivery using the Cyberknife with multileaf collimator, for the case of SBRT. The choice of trajectory is a key aspect of arc delivery, but as trajectory selection is an extensive subject, the reader is referred to previous studies for details. 8,10,15 In the present study, a previously investigated arc trajectory 15 is used to provide suitable dynamic baseline plans in the known range of collision-free operation of the Cyberknife robot. Based on published data and other estimates from Accuray, the expected variations in dose rate, leaf positioning, and robot target position are established and incorporated into models of these parameters (Fig. 1). The consequent impact on the dose distribution is calculated for several cases using these models. Beam delivery between control points is included in the model so as to give the most accurate possible estimation of the dose distribution.

2.A. Patient cases
Four patient cases were retrospectively investigated in this study, as follows: a Prostate A, planned for a homogeneous distribution of dose for a prescription of 36.25 Gy in five fractions according to RTOG 0938. 16 PTV volume 113 cm 3 . b Prostate B, the same case planned for a brachytherapystyle dose distribution using a prescription of 38.0 Gy in four fractions. 17 An SBRT technique was used in all cases, with at least 95% of the planning target volume (PTV) being required to receive the prescribed dose. The plans used a fixed isocenter, located at the center of the PTV.

2.B. Arc trajectory
The arc trajectory used for this study is shown in Fig. 2. It was a purpose-made path, distinct from the standard body path used by Cyberknife, and using different node positions. It consisted of 104 nodes, or control points, with a spacing of 5°in robot orientation. The trajectory was designed by an Accuray heuristic to give a uniform coverage of the space around the patient, while respecting hardware constraints on collisions and cabling. The performance of the arc in comparison with the standard body path was previously investigated. 15

2.C. Dose calculation and optimization scheme
Dose was calculated as: where d i was the dose at voxel i in the patient model, d ij was the dose contributed by beamlet j to voxel i, and w j was the Medical Physics, 47 (4), April 2020 beamlet weight. The dose calculation was based on an Accuray-supplied pencil-beam algorithm embedded into an inhouse computation and optimization framework. The commercial planning system for Cyberknife was not used in this study. The dose-influence matrices, d ij , of Eq. (1) were calculated by using this algorithm to calculate the dose for a series of bixel-sized fields. The dose grid was 2 9 CT pixel size in the transaxial direction and CT slice spacing in the longitudinal direction. Dose voxels which received less than 0.015% of the maximum dose of each d ij component were neglected so as to minimize the size of the dose matrices. Plans were produced for all of the cases using a three-step optimization scheme which produced fluence maps, then sequenced these into deliverable apertures, before performing direct aperture optimization. The resolution of the fluence map was 7.7 mm 9 5.0 mm at a nominal source-axis distance of 800 mm, the 7.7 mm being equal to two leaf widths so that the MLC leaves could be paired. 8 This was a practical feature designed to reduce the number of fluence bixels and hence reduce memory requirements and increase optimization speed. The optimization minimized an objective function, F, defined as: where the indices, k, referred to individual anatomical structures, each with objective value f k : In this equation, the variables a k and b k referred to the importance factors for structure k. Fluence was optimized using the L-BFGS optimization scheme, which was a standard gradient descent method using a recursion relation for calculation of search directions, so as to avoid having to store a large and memory-intensive inverse Hessian matrix. 22 A fluence map was produced by the optimizer for every third beam orientation. After fluence optimization using 40 iterations, the plan was sequenced using a well-established sequencer, 23 with three apertures being used to account for each fluence map. Two of these apertures were redistributed to the beam orientations either side of that for which fluence was optimized. The L-BFGS method was again used for the direct aperture optimization, with each aperture shape being used to define which fluence pixels were active, and then with the corresponding bixel doses being used to calculate dose and search directions. 22,24 The optimization itself was therefore based on dose calculated by summation of bixel doses.
During the direct aperture optimization, the method of Christiansen et al. 25 was used to model the motion of MLC leaves between control points. The motion of one leaf in a leaf pair produced a fluence ramp in one direction, and the motion of the opposing leaf produced a ramp of fluence in the other direction, and the total fluence was calculated as the difference between the two fluence patterns.
The direct aperture optimization also took into account the delivery constraints of the Cyberknife. 15 The approach taken was that the robot speed was given first priority, defining the delivery time for the arc. The robot speed was taken as 60 mms À1 , equivalent to the slowest speed observed in practice, so that the time to traverse 5°of arc was calculated to be 1.5 s. The MLC movement was then limited so that the time for MLC motion should not exceed the time taken by the robot to make its movement. The MLC leaf speed was taken to be 33 mms À1 , which was slightly faster than in current clinical practice, but a speed which was advised by Accuray to be achievable using the current generation of hardware. In the 1.5 s taken for the robot to traverse between nodes, the MLC leaves were therefore able to move 50 mm, and this limit was used by the optimizer. There was no limit to the minimum monitor units per segment, implying that the beam could be paused if required. There was also no limit on the maximum monitor units per segment, which meant that if the time for delivery of the monitor units was large, it could exceed the time taken for the robot motion. In practice, the robot would slow down in this situation. Further details of the optimization scheme are given elsewhere. 15 After the optimization, dose was recalculated using the Accuray dose calculation in the local computational framework, based on complete apertures rather than on the sums of bixel doses used during the optimization. The dose grid resolution used for this recalculation was the same as for the bixel-based dose calculation used during optimization. The recalculated aperture dose had a very similar relative dose distribution to the original bixel-based dose produced by the optimizer, but was seen to be offset by a case-specific factor of several percent. This resulted from the difference between calculating dose due to MLC-shaped apertures and due to collections of individual 7.7 mm 9 5 mm apertures. The aperture dose was actually the more accurate dose, but as the treatment plans had been created and clinical constraints had been met using bixel-based doses, the recalculated plan was renormalized to return the dose approximately to the bixelbased dose distribution, so that the clinical constraints were met as closely as possible. Since the difference between the aperture-based and bixel-based doses was a relative scaling, this renormalization was accomplished by multiplying the aperture dose by a simple scale factor. The same renormalization factor was used for all plans in each patient case, so that this did not confound any changes in dose resulting from the errors examined in the study. Throughout the study, the monitor units calculated by the optimization engine were used, without any application of a scale factor.
In the final dose recalculation, additional nodes were used to model the robot orientation and MLC leaf position as accurately as possible. Four intermediate nodes were added between each pair of nodes, such that the additional nodes, together with the second of the original pair of nodes, formed a set of five interpolated nodes. 15 A further 20 equally weighted interpolated apertures, allowing for MLC leaf motion, were calculated between each of the interpolated nodes, and dose was then calculated and summed for all of the interpolated apertures (Fig. 3). This detailed modeling of the MLC motion was found previously to be essential for Medical Physics, 47 (4), April 2020 accurate dose calculation with relatively large movements of small apertures. 26

2.D. Uncertainty models
The uncertainty models are described below. The uncertainties were considered independently so as to determine their individual contributions to the delivery uncertainty.

2.D.1. Dose rate uncertainty
Dose rate variations were modeled as a series of pseudorandom errors defined at integer numbers of the error period, which was set to 10, 30 or 60 s in turn (Fig. 4). The first of these values represented a rapidly drifting error and the last represented a slowly drifting error. The errors were sampled from a uniform distribution spanning from ÀE to +E, where E was the error magnitude, which was either 0.01 (i.e., 1%), 0.025 or 0.05. Between the error sampling points, the error was determined by linear interpolation. Large and abrupt changes (e.g. due to arcing) were ignored.
Addition of the random errors to the treatment arcs was carried out using the 104-node optimized plans, before addition of interpolated nodes. The time, t n , at which segment n (n = 1. . .104) of the arc took place was calculated according to the formula: where M i was the number of monitor units per fraction at node i, D = 1000 monitor units per minute was the maximum dose rate of the linear accelerator, and T R = 1.5 s was the robot traversal time. In other words, the time taken to deliver each segment of the arc was equal to the fixed traversal time T R , unless the monitor units, delivered at maximum dose rate, required additional time for delivery, in which case the time was equal to the time to deliver the monitor units. According to the time that a given control point was delivered, the appropriate value of error was determined from the error model. The monitor units for that control point were changed by the error: where M i were the monitor units before addition of the error, M E were the monitor units after addition of the error, and E was the error. Strictly speaking, the delivery of M E monitor units took a different length of time to delivering M i monitor units, thereby influencing the time at which the error, E, should be sampled, but this secondary effect was neglected.
Using M E , five interpolated nodes, each with 20 interpolated apertures, were calculated as described above. The three error periods and three error magnitudes gave nine combinations of error. For each combination of error period and magnitude, the calculation was repeated 10 times to give an indication of the spread of outcomes resulting from the selection of random errors.

2.D.2. MLC position uncertainty
To model MLC leaf position uncertainty, at each control point of the arc (before interpolation), the position of each MLC leaf was moved 0.25 mm towards its position at the previous control point. The MLC leaf positions and the movement of 0.25 mm were defined at the nominal sourceaxis distance of 800 mm. If the position of the MLC leaf at the previous control point was <0.25 mm away from its position at the current control point, the leaf was only moved as far as the position at the previous control point. This limitation occurred for between 4% and 13% of leaf movements, depending on the patient case, when considering all MLC leaves at all control points and neglecting closed leaf pairs. No change in position occurred for the first control point in the arc. The direction of MLC perturbation (i.e., use of a lag in position) and the representative magnitudes of uncertainty were based on worst-case observations of the physical Cyberknife device.
After the change in position had been effected, five interpolated nodes, each with 20 intermediate apertures, were introduced and the dose was calculated, based on the interpolated plan. The experiment was repeated with a lag of 0.5 mm.

2.D.3. Robot pointing uncertainty
Each robot position was specified by a source point and target point, the target point corresponding to the isocenter on a conventional linear accelerator, but in this case varying in source-target distance. To model the uncertainty in robot pointing, at each control point of the arc (before interpolation), the target point of the node was moved according to a uniformly sampled random direction with respect to the multileaf collimator. The error directions of the interpolated control points were then computed by linear interpolation between these directions. The interpolated directions were chosen to take the shortest route around the error circle (Fig. 5).
The magnitude of the error at the nodes before interpolation was uniformly sampled from 0. . .N, where N was either 0 or 1 mm. It was expected that the positioning of the robot would be less accurate between nodes than at the actual node positions, since between nodes, the control system would be carrying out some form of interpolation. An additional uniformly sampled error of 0. . .M was therefore added midway between the nodes, to reflect the additional uncertainty between nodes, where M was either 0.5 mm or 1.0 mm. The total error was thus (n 1 + n 2 )/2 + m at the mid-point between the nodes, where n 1 and n 2 were randomly sampled from a uniform distribution of width N and m was randomly sampled from a uniform distribution of width M. The magnitude of the error between the original nodes and the midpoint between them was calculated by linear interpolation, so that the total error changed smoothly from n 1 to (n 1 + n 2 )/2 + m and then smoothly to n 2 .

2.D.4. Combination of uncertainties
The above uncertainties were applied individually to the treatment plans, but their impact if they were all present simultaneously was also investigated. A dose rate uncertainty of magnitude 0.025 with a time period of 30 s, an MLC lag of 0.25 mm and a robot pointing uncertainty of 1.0 mm at nodes with an additional uncertainty of 0.5 mm between nodes was used for this investigation. All of the uncertainties were applied to the same treatment plan and the effect on the dose distribution was observed. Three runs were carried out to assess the impact of the random variations.
A summary of the numerical values of uncertainty used in the model, together with their sources, is given in Table I.

3.A. Optimized treatment plans
The calculated monitor units and estimated delivery times for the five arcing treatment plans are given in Table II. The estimated delivery times are calculated using Eq. (4) for the final node, and the corresponding times for static delivery are also shown. Representative dose distributions for the treatment plans without delivery uncertainties are shown in a previous study, together with more details of comparison with static delivery. 15

3.B. Dose rate uncertainty
The variation in monitor units per fraction with control point for the prostate A case, and the random uncertainty from the uncertainty model, are shown in Fig. 6. The error distribution [ Fig. 6(b)] varies from run to run, due to the selection of a different set of random numbers each time. The   results of 10 runs of the calculation are shown in Fig. 7 for 5% uncertainty at time periods of 10, 30, and 60 s. The graphs show a range of doses to the PTV, and this also occurs for the other patients. The statistics for the runs giving the lowest and highest doses are shown for all patients in Table III with 5% uncertainty and a time period of 60 s, that is, for the situation in Fig. 7(c). The variation in PTV dose is accompanied by a variation in total monitor units per fraction. The predominant effect of the random variation of the monitor units is to change the total number of monitor units delivered, and this correlates closely with the variation in dose to the PTV (Fig. 8). The longer time period results in fewer random samples of error over the length of the delivery, so that the errors are less likely to cancel than with many samples. Consequently, the longer time period gives rise to a larger fluctuation in PTV dose than with the shorter time period.
The other error magnitudes of 0.01 and 0.025 give similar results to an error magnitude of 0.05, but with a correspondingly scaled down effect. The standard deviation of the error is expected to be proportional to the error magnitude, so a curve of the form r ¼ M ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi aT 2 þ bT p is fitted to the graphs of standard deviation vs error magnitude, where r is the standard deviation, in Gy, of the PTV mean dose, M is the magnitude of the dose rate error, as a dimensionless proportion of the correct dose rate, and T is the time period, in seconds, of the dose rate variation. The values of a and b are constant and to be determined. By first plotting r against M for several values of T and then fitting a polynomial of the form aT 2 þ bT to a plot of the squared gradient of these curves against T, the behavior is found to be modeled by the equation: As a percentage of the PTV mean dose for the plan without errors, the relationship between standard deviation and error magnitude and period is: where r is now the standard deviation of the PTV mean dose, as a percentage of the correct PTV mean dose. The overall variation in PTV mean dose is shown in Fig. 9 for the prostate A case, with the modeled results using this equation included. The results of applying Eq. (7) to all of the patient cases are shown in Table IV for the largest error magnitude.  For context, the observed data are based on a sample size of 10, so the standard error of the actual standard deviation, r, shown in Fig. 9 and Table IV is r ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2 n À 1 ð Þ p % r=4. In other words, the inherent variability of the data is approximately one quarter of the standard deviations given.

3.C. MLC position uncertainty
Dose-volume histograms showing the impact of a 0.5 mm MLC lag are shown in Fig. 10 for the prostate B and lung cases. (Due to the large overlap of the PTV with the proximal bronchus in the lung case, the plan is unable to meet the constraint that the dose to 4 cm 3 of proximal bronchus should not exceed 18 Gy.) The MLC lag has the effect of modifying the dose to the PTV by around 1%, with a slightly larger effect on the urethra in the prostate B case. Table V shows the statistics for all of the cases, for both 0.25 and 0.5 mm lag.

3.D. Robot pointing uncertainty
The effect of a randomly sampled robot pointing error of 1.0 mm at the delivery nodes, with an additional randomly

Case
Magnitude Period (s) Predicted % SD Actual % SD Medical Physics, 47 (4), April 2020 sampled error 1.0 mm between the nodes is shown in Fig. 11 for the liver and partial breast cases. The effect of the error is to modify the PTV dose slightly, with an increase in dose usually being observed. Table VI shows the statistics for all cases and robot pointing errors.

3.E. Combination of uncertainties
The effect of a dose rate uncertainty of magnitude 0.025, an MLC lag of 0.25 mm and a robot pointing uncertainty of 1.0 mm with an additional 0.5 mm between nodes, is to adjust the PTV mean dose by around 1%, mostly with the minimum and maximum doses following the trend of the mean dose. The results of three runs for each patient are shown in Table VII.

DISCUSSION
The ability of the Cyberknife system to treat from a variety of noncoplanar beam directions 27 is useful for high-quality treatment of SBRT. [1][2][3][4] Treating from such a range of orientations is time-consuming, so that the total treatment time is long. However, several recent studies have shown that use of a dynamic arc delivery can reduce the treatment time by a factor of approximately two. For example, Kearney et al. 10 report for prostate and brain patients a speedup of 1.5 AE 0.3, depending on the parameters used by the optimizer, and Bedford et al. 15 report a median speedup of 1.90 (range 1.53 to 2.36). Similarly, using a C-arm linear accelerator, Wild et al. 24 predict a delivery time of 6.5 min on average for noncoplanar VMAT. The choice of arc trajectory is likely to have an impact on the delivery time. 28 Since the Cyberknife does not currently deliver radiation in dynamic mode, it is not possible to measure the dosimetric accuracy of dynamic delivery in the same way that is possible for a C-arm linear accelerator. 28 However, several recent studies have reported on the geometric accuracy of the device, which have been used in the present study to give estimates of the expected dosimetric performance. [11][12][13][14] Variations in dose rate are shown to have negligible effect on the shape of the dose distribution, but rather to impact on the total monitor units delivered, which in turn affect the scaling of the plan, most noticeably in the mean dose to the PTV. There is a corresponding effect on the organs at risk in proportion to the dose received by these organs at risk. For example, proximal structures such as the urethra and the proximal bronchus are impacted similarly to the PTV, whereas distal structures are hardly affected at all. The magnitude of the dosimetric error determined in this study is likely to be a worst-case estimate, as the variation in the monitor units produced by the optimizer is quite high (Fig. 6), with some nodes not delivering any dose at all.
The impact of an MLC lag is shown to cause a small change in PTV dose. This is almost always in the form of an increase in PTV dose of around 1% for a 0.5 mm lag in MLC leaf position. This is in accord with the work of Christiansen et al., 25   work uses a dynamic beam model throughout, in order to model the dose delivery as accurately as possible, but the comparison of Christiansen et al. 25 with a static situation is very similar to a situation with a severe MLC lag.
In terms of practical measurements on the current generation of Cyberknife, perhaps the most that is known is for the accuracy of robot pointing at static nodes. 11,12 The effect of a random error in the order of 1.0 mm at control points, with an additional error between control points, is to modify the PTV dose by approximately 1%. An uncertainty in robot position is similar in nature to an uncertainty in MLC leaf position, so the similarities between the results for MLC lag and robot position are rational. The uncertainty in robot pointing is modeled in this study as a shift in the target point, with the robot itself positioned correctly (but not orientated correctly). The results might be affected by using a model in which the robot itself is positioned incorrectly. However, Wang and Nelson 14 show that the overall result of translation and rotation errors is a shift in the target point of less than 1 mm, so the model used in the present study is representative.
Taking these results together, the study shows that the Cyberknife is expected to be able to deliver dynamic arc beams with a dosimetric accuracy of the order of 1-2%. This is comparable to accuracy studies for Cyberknife with static beam delivery, on the one hand, and conventional linear accelerators delivering VMAT, on the other hand. For Cyberknife with static beam delivery, Dieterich et al. 29 and Moore et al. 30 recommend a dosimetric accuracy of 90% within 2% and 2 mm in comparison with the treatment planning system. Some of this tolerance is used by the treatment planning system, but the accuracy predicted by the present study is within this order of magnitude. Similarly, for coplanar VMAT delivery, [31][32][33] Ling et al. 34 report MLC leaf positions generally within 0.5 mm of their expected position, and relative dose accuracy of 0.7%. Bedford and Warrington 35 report a dynamically delivered dose generally within 2% of the same plan delivered statically, and a gamma agreement of 3% and 3 mm for delivery of a complete treatment plan. Mans et al. 36 recommend a dose rate dependence with gantry angle of better than 0.5% and MLC leaf positioning of better than 1.0 mm and preferably better than 0.5 mm. For comparison with a treatment planning system, 90% of the delivered dose should be within 3% and 3 mm.
The results have been obtained for PTV volumes between 14 cm 3 (lung case) and 113 cm 3 (prostate A case). A larger relative impact is observed with MLC position uncertainty for the lung case than the prostate A case, which may be partly due to the difference in volume. It is very likely that different results may occur for very small stereotactic  volumes below 10 cm 3 , in which case the present results should be treated with caution. The choice of arc trajectory may have some impact on the dosimetric accuracy. The trajectory used in this study is the result of a heuristic which aims to provide a uniform coverage of the beam orientation space, while taking into account cabling requirements and collision avoidance. This trajectory is asymmetric, and has been shown in a previous study 15 to compare well with the standard (symmetric) body path for the Cyberknife in terms of dose distribution and conformity. The trajectory affects the number of monitor units, the MLC shapes and the robot positions, and is therefore likely to affect the delivery uncertainties and hence the dosimetric results of the study to some degree. However, it is thought that the principal conclusions of this study are likely to be similar to whichever trajectory is chosen.
Further improvement in the delivery accuracy could be obtained by the use of robust treatment planning. 37 Although this is usually applied to overcome uncertainties in target position, setup uncertainties, or range uncertainties, it could be used to overcome the limitations of linear accelerator performance. By informing the optimization engine of the possible scenarios that can occur during treatment delivery, the treatment plan can be designed to give an acceptable dose distribution in most or all situations. For this reason, models of treatment delivery accuracy are important.

CONCLUSIONS
Estimated uncertainties in dose rate, MLC leaf position, and robot target position are shown to have a dosimetric impact of around 1-2% during the arc delivery of SBRT using Cyberknife. The dose rate uncertainty affects the total number of monitor units delivered, which in turn affects the delivered dose. A lag of around 0.5 mm in MLC leaf position affects the delivered dose by approximately 1%, usually with an increase in dose, and an uncertainty in robot pointing has a similar effect. The dose to the PTV is affected most, with a lesser impact on organs at risk that receive less dose in the treatment plan.