Some considerations concerning volume-modulated arc therapy: a stepping stone towards a general theory

A polar-coordinate-based pencil beam algorithm for VMAT dose computation with high-resolution gantry angle sampling

PURPOSE

Most commercially available treatment planning systems (TPSs) approximate the continuous delivery of volumetric modulated arc therapy (VMAT) plans with a series of discretized static beams for treatment planning, which can make VMAT dose computation extremely inefficient. In this study, we developed a polar-coordinate-based pencil beam (PB) algorithm for efficient VMAT dose computation with high-resolution gantry angle sampling that can improve the computational efficiency and reduce the dose discrepancy due to the angular under-sampling effect.

METHODS AND MATERIALS

6 MV 1 × 1 mm² pencil beams were simulated on a uniform cylindrical phantom under an EGSnrc Monte Carlo (MC) environment. The MC-generated PB kernels were collected in the polar coordinate system for each bixel on a 40 × 40 cm² fluence map and subsequently fitted via a series of Gaussians. The fluence was calculated using a detectors' eye view with off-axis and MLC transmission factors corrected. Doses of VMAT arc on the phantom were computed by summing the convolution results between the corresponding PB kernels and fluence for each bixel in the polar coordinate system. The convolution was performed using Fast Fourier Transform to expedite the computing speed. The calculated doses were converted to the Cartesian coordinate system and compared with the reference dose computed by a collapsed cone convolution (CCC) algorithm of the TPS. A heterogeneous phantom was created to study the heterogeneity corrections using the proposed algorithm. Ten VMAT arcs were included to evaluate the algorithm performance. Gamma analysis and computation complexity theory were used to measure the dosimetric accuracy and computational efficiency, respectively.

RESULTS

The dosimetric comparisons on the homogeneous phantom between the proposed PB algorithm and the CCC algorithm for ten VMAT arcs demonstrate that the proposed algorithm can achieve a dosimetric accuracy comparable to that of the CCC algorithm with average gamma passing rates of 96% (2%/2mm) and 98% (3%/3mm). In addition, the proposed algorithm can provide better computational efficiency for VMAT dose computation using a PC equipped with a 4-core processor, compared to the CCC algorithm utilizing a dual 10-core server. Moreover, the computation complexity theory reveals that the proposed algorithm has a great advantage with regard to computational efficiency for VMAT dose computation on homogeneous medium, especially when a fine angular sampling rate is applied. This can support a reduction in dose errors from the angular under-sampling effect by using a finer angular sampling rate, while still preserving a practical computing speed. For dose calculation on the heterogeneous phantom, the proposed algorithm with heterogeneity corrections can still offer a reasonable dosimetric accuracy with comparable computational efficiency to that of the CCC algorithm.

CONCLUSIONS

We proposed a novel polar-coordinate-based pencil beam algorithm for VMAT dose computation that enables a better computational efficiency while maintaining clinically acceptable dosimetric accuracy and reducing dose error caused by the angular under-sampling effect. It also provides a flexible VMAT dose computation structure that allows adjustable sampling rates and direct dose computation in regions of interest, which makes the algorithm potentially useful for clinical applications such as independent dose verification for VMAT patient-specific QA. This article is protected by copyright. All rights reserved.

A sliding-window approach for improved VMAT dose calculation accuracy

The continuous delivery of volumetric modulated arc therapy (VMAT) plans is usually approximated by discrete apertures at evenly-spaced gantry angles for dose calculation purposes. This approximation can potentially lead to large dose calculation errors if the gantry angle spacings are large and/or there are large changes in the MLC apertures from one control point (CP) to the next. In this work, we developed a sliding-window (SW) method to improve VMAT dose calculation accuracy. For any 2 adjacent VMAT CPs ni and ni + 1, the dose distribution was approximated by a 2-CP SW IMRT beam with the starting MLC positions at CP ni and ending MLC positions at CP ni + 1, with the gantry angle fixed in the middle of the 2 VMAT CPs. Therefore, a VMAT beam with N CPs was approximated by a SW plan with N-1 SW beams. To validate the method, VMAT plans were generated for 10 patients in Pinnacle using 4° gantry spacing. Each plan was converted to a SW plan and dose was recalculated. Another VMAT plan, with 1° gantry spacing, was created by interpolating the original VMAT beam. The original plans were delivered on an Elekta Versa HD and measured with ArcCHECK. For both the isodose distribution and DVH, there were significant differences between the original VMAT plan and either the SW or the interpolated plan. However, they were indistinguishable between the SW and the interpolated plans. When compared with measurement, the average passing rates of the original VMAT plans were 87.3 ± 2.8% and 93.1 ± 1.0% for the 5 HN and 5 spine SBRT cases, respectively. On the other hand, the passing rates for both the VMAT1 and SW plans were above 95% for all the 10 cases studied. The dose calculation times of the original VMAT plans and the SW plans were very similar. We conclude that the proposed SW approach improves VMAT dose calculation accuracy without increase in dose calculation time.

EPID sensitivity to delivery errors for pre-treatment verification of lung SBRT VMAT plans

PURPOSE

To study the sensitivity of an Electronic Portal Imaging Device (EPID) in detecting delivery errors for VMAT lung stereotactic body radiotherapy (SBRT) using the Collapsed Arc method.

METHODS

Baseline VMAT plans and plans with errors intentionally introduced were generated for 15 lung SBRT patients. Three types of errors were introduced by modifying collimator angles and multi-leaf collimator (MLC) field sizes (MLCFS) and MLC shifts by ±5, ±2, and ±1° or millimeters. A total of 103 plans were measured with EPID on an Elekta Synergy Linear Accelerator (Agility MLC) and compared to both the original treatment planning system (TPS) Collapsed Arc dose matrix and the no-error plan baseline EPID measurements. Gamma analysis was performed using the OmniPro-I'mRT (IBA Dosimetry) software and gamma criteria of 1%/1 mm, 2%/1 mm, 2%/2 mm, and 3%/3.

RESULTS

When the error-introduced EPID measured dose matrices were compared to the TPS matrices, the majority of simulated errors were detected with gamma tolerance of 2%/1 mm and 1%/1 mm. When the error-introduced EPID measured dose matrices were compared to the baseline EPID measurements, all the MLCFS and MLC shift errors, and ±5°collimator errors were detected using 2%/1 mm and 1%/1 mm gamma criteria.

CONCLUSION

This work demonstrates the feasibility and effectiveness of the collapsed arc technique and EPID for pre-treatment verification of lung SBRT VMAT plans. The EPID was able to detect the majority of MLC and the larger collimator errors with sensitivity to errors depending on the gamma tolerances.

Angular under-sampling effect on VMAT dose calculation: An analysis and a solution strategy

PURPOSE

Most VMAT algorithms compute the dose on discretized apertures with small angular separations for practical reasons. However, machines deliver the VMAT dose with a continuously moving MLC and gantry and a continuously changing dose rate. The computed dose can deviate from the delivered dose, especially if no, or loose, MLC movement constraints are applied for the VMAT optimization. The goal of this paper is to establish a simplified mathematical model to analyze the discrepancy between the VMAT plan calculation dose and the delivered dose and to provide a reasonable solution for clinical implementation.

METHODS

A simplified metric is first introduced to describe the discrepancy between doses computed with discretized apertures and a continuous delivery model. The delivery fluences were formed separately for six different leaf movement scenarios. The formula was then rewritten in a more general form. The correlation between discretized and continuous fluence is summarized using this general form. The Fourier analysis for the impacts from three separate factors - dose kernel width, aperture width, aperture distance - to the dose discrepancy is also presented in order to provide insight into the dose discrepancy caused by under-sampling in the frequency domain. Finally, a weighting-based interpolation (WBI) algorithm, which can improve the aperture interpolation efficiency, is proposed. The associated evaluation methods and criteria for the proposed algorithm are also given.

RESULTS

The comparisons between the WBI algorithm and the equal angular interpolation (EAI) method suggested that the proposed algorithm has a great advantage with regard to aperture number efficiency. To achieve a 90% gamma passing rate using the dose computed with apertures generated with 0.5° EAI, with the initial optimization apertures as the standard for the comparison, the WBI needs only 66% and 54% of the aperture numbers that the EAI method needs for a 2° and a 4° angular separation of the VMAT optimization, respectively. The results also suggested that the weighted dose error index value, Θ, can be used as a stopping criterion for an interpolation algorithm, e.g., WBI or EAI, or as an indicator for sampling level evaluations. The phantom results indicate that the gamma passing rate decreases with increasing depth, from the phantom surface to the iso center, for the plans computed with under-sampled apertures. No obvious variation trends were observed for the plans computed with well-sampled apertures.

CONCLUSIONS

The mathematical analysis suggests that the dose discrepancies due to under-sampling are strongly correlated with the aperture width, the distance between apertures, and the width of the dose kernel. The WBI algorithm proves to be an efficient aperture interpolation strategy and is useful for dose computation of VMAT plans.

Analysis of dose distribution in organs at risk in patients with prostate cancer treated with the intensity-modulated radiation therapy and arc technique

This study describes a comparative analysis of treatment plans in 48 patients with prostate cancer treated with ionizing radiation. Each patient was subjected to the intensity-modulated radiation therapy (IMRT) and arc technique. In each treatment plan, the organs at risk were assessed: the urinary bladder, rectum and heads of the femur, as well as the volume of normal tissue. The following features were compared: treatment time, conformity indices for the planning target volume, mean doses and standard deviation in organs at risk, and organ volumes for each particular dose. The treatment period in the arc technique is 13.7% shorter than in the IMRT technique. Comparing the results of the IMRT and arc techniques (arc vs. IMRT), the mean values were 29.21 ± 12.91 Gy versus 28.36 ± 13.79 Gy for the bladder, 20.36 ± 3.16 Gy versus 18.17 ± 5.11 Gy for the right femoral head, and 18.98 ± 3.28 Gy versus 16.67 ± 5.15 Gy for the left femoral head. For the rectum, lower values were obtained after application of the arc technique, not the IMRT technique: 35.84 ± 12.28 Gy versus 35.90 ± 13.05 Gy. The results indicate that the applied therapy has a statistically significant influence on the volume for a particular dose with regard to the urinary bladder. It is advisable to apply the IMRT technique to patients who need the femur heads and urinary bladder protected by exposing them to low irradiation doses.

Dosimetric and radiobiological characterizations of prostate intensity-modulated radiotherapy and volumetric-modulated arc therapy: A single-institution review of ninety cases

This study reviewed prostate volumetric-modulated arc therapy (VMAT) plans with intensity-modulated radiotherapy (IMRT) plans after prostate IMRT technique was replaced by VMAT in an institution. Characterizations of dosimetry and radiobiological variation in prostate were determined based on treatment plans of 40 prostate IMRT patients (planning target volume = 77.8–335 cm³) and 50 VMAT patients (planning target volume = 120–351 cm³) treated before and after 2013, respectively. Both IMRT and VMAT plans used the same dose-volume criteria in the inverse planning optimization. Dose-volume histogram, mean doses of target and normal tissues (rectum, bladder and femoral heads), dose-volume points (D_99% of planning target volume; D_30%, D_50%, V_{30 Gy} and V_{35 Gy} of rectum and bladder; D_5%, V_{14 Gy}, V_{22 Gy} of femoral heads), conformity index (CI), homogeneity index (HI), gradient index (GI), prostate tumor control probability (TCP), and rectal normal tissue complication probability (NTCP) based on the Lyman-Burman-Kutcher algorithm were calculated for each IMRT and VMAT plan. From our results, VMAT plan was found better due to its higher (1.05%) CI, lower (0.83%) HI and (0.75%) GI than IMRT. Comparing doses in normal tissues between IMRT and VMAT, it was found that IMRT mostly delivered higher doses of about 1.05% to the normal tissues than VMAT. Prostate TCP and rectal NTCP were found increased (1%) for VMAT than IMRT. It is seen that VMAT technique can decrease the dose-volume evaluation criteria for the normal tissues. Based on our dosimetric and radiobiological results in treatment plans, it is concluded that our VMAT implementation could produce comparable or slightly better target coverage and normal tissue sparing with a faster treatment time in prostate radiotherapy.

Dose calculation for hypofractionated volumetric-modulated arc therapy: approximating continuous arc delivery and tongue-and-groove modeling

Hypofractionated treatments generally increase the complexity of a treatment plan due to the more stringent constraints of normal tissues and target coverage. As a result, treatment plans contain more modulated MLC motions that may require extra efforts for accurate dose calculation. This study explores methods to minimize the differences between in-house dose calculation and actual delivery of hypofractionated volumetric-modulated arc therapy (VMAT), by focusing on arc approximation and tongue-and-groove (TG) modeling. For dose calculation, the continuous delivery arc is typically approximated by a series of static beams with an angular spacing of 2°. This causes significant error when there is large MLC movement from one beam to the next. While increasing the number of beams will minimize the dose error, calculation time will increase significantly. We propose a solution by inserting two additional apertures at each of the beam angle for dose calculation. These additional apertures were interpolated at two-thirds' degree before and after each beam. Effectively, there were a total of three MLC apertures at each beam angle, and the weighted average fluence from the three apertures was used for calculation. Because the number of beams was kept the same, calculation time was only increased by about 6%-8%. For a lung plan, areas of high local dose differences (> 4%) between film measurement and calculation with one aperture were significantly reduced in calculation with three apertures. Ion chamber measurement also showed similar results, where improvements were seen with calculations using additional apertures. Dose calculation accuracy was further improved for TG modeling by developing a sampling method for beam fluence matrix. Single element point sampling for fluence transmitted through MLC was used for our fluence matrix with 1 mm resolution. For Varian HDMLC, grid alignment can cause fluence sampling error. To correct this, transmission volume averaging was applied. For three paraspinal HDMLC cases, the average dose difference was greatly reduced in film and calculation comparisons with our new approach. The gamma (3%, 3 mm) pass rates have improved significantly from 74.1%, 90.0%, and 90.4% to 99.2%, 97.9%, and 97.3% for three cases, for calculation without volume averaging and calculation with volume averaging, respectively. Our results indicate that more accurate MLC leaf position and transmission sampling can improve accuracy and agreement between calculation and measurement, and are particularly important for hypofractionated VMAT that consists of large MLC movement.

Total-variation regularization based inverse planning for intensity modulated arc therapy

Intensity modulated arc therapy (IMAT) delivers conformal dose distributions through continuous gantry rotation with constant or variable speed while modulating the field aperture shape and weight. The enlarged angular space and machine delivery constraints make inverse planning of IMAT more intractable as compared to its counterpart of fixed gantry IMRT. Currently, IMAT inverse planning is being done using two extreme methods: the first one computes in beamlet domain with a subsequent arc leaf sequencing, and the second proceeds in machine parameter domain with entire emphasis placed on a pre-determined delivery method without exploring potentially better alternative delivery schemes. Towards truly optimizing the IMAT treatment on a patient specific basis, in this work we propose a total-variation based inverse planning framework for IMAT, which takes advantage of the useful features of the above two existing approaches while avoiding their shortcomings. A quadratic optimization algorithm has been implemented to demonstrate the performance and advantage of the proposed approach. Applications of the technique to a prostate case and a head and neck case indicate that the algorithm is capable of generating IMAT plans with patient specific numbers of arcs efficiently. Superior dose distributions and delivery time are achieved with a maximum number of apertures of three for each field. As compared to conventional beamlet-based algorithms, our method regularizes the field modulation complexity during optimization, and permits us to obtain the best possible plan with a pre-set modulation complexity of fluences. As illustrated in both prostate and head-and-neck case studies, the proposed method produces more favorable dose distributions than the segment-based algorithms, by optimally accommodating the clinical need of intensity modulation levels for each individual field. On a more fundamental level, our formulation preserves the convexity of optimization and makes the search of the global optimal solution possible with a deterministic method.

Commissioning and comprehensive evaluation of the ArcCHECK cylindrical diode array for VMAT pretreatment delivery QA

We present commissioning and comprehensive evaluation for ArcCHECK as a QA equipment for volumetric-modulated arc therapy (VMAT), using the 6 MV photon beam with and without the flattening filter, and the SNC patient software (version 6.2). In addition to commissioning involving absolute dose calibration, array calibration, and PMMA density verification, ArcCHECK was evaluated for its response dependency on linac dose rate, instantaneous dose rate, radiation field size, beam angle, and couch insertion. Scatter dose characterization, consistency and symmetry of response, and dosimetry accuracy evaluation for fixed aperture arcs and clinical VMAT patient plans were also investigated. All the evaluation tests were performed with the central plug inserted and the homogeneous PMMA density value. Results of gamma analysis demonstrated an overall agreement between ArcCHECK-measured and TPS-calculated reference doses. The diode based field size dependency was found to be within 0.5% of the reference. The dose rate-based dependency was well within 1% of the TPS reference, and the angular dependency was found to be ± 3% of the reference, as tested for BEV angles, for both beams. Dosimetry of fixed arcs, using both narrow and wide field widths, resulted in clinically acceptable global gamma passing rates on the 3%/3mm level and 10% threshold. Dosimetry of narrow arcs showed an improvement over published literature. The clinical VMAT cases demonstrated high level of dosimetry accuracy in gamma passing rates.

Optimal set of grid size and angular increment for practical dose calculation using the dynamic conformal arc technique: a systematic evaluation of the dosimetric effects in lung stereotactic body radiation therapy

PURPOSE

To recommend the optimal plan parameter set of grid size and angular increment for dose calculations in treatment planning for lung stereotactic body radiation therapy (SBRT) using dynamic conformal arc therapy (DCAT) considering both accuracy and computational efficiency.

MATERIALS AND METHODS

Dose variations with varying grid sizes (2, 3, and 4 mm) and angular increments (2°, 4°, 6°, and 10°) were analyzed in a thorax phantom for 3 spherical target volumes and in 9 patient cases. A 2-mm grid size and 2° angular increment are assumed sufficient to serve as reference values. The dosimetric effect was evaluated using dose-volume histograms, monitor units (MUs), and dose to organs at risk (OARs) for a definite volume corresponding to the dose-volume constraint in lung SBRT. The times required for dose calculations using each parameter set were compared for clinical practicality.

RESULTS

Larger grid sizes caused a dose increase to the structures and required higher MUs to achieve the target coverage. The discrete beam arrangements at each angular increment led to over- and under-estimated OARs doses due to the undulating dose distribution. When a 2° angular increment was used in both studies, a 4-mm grid size changed the dose variation by up to 3-4% (50 cGy) for the heart and the spinal cord, while a 3-mm grid size produced a dose difference of <1% (12 cGy) in all tested OARs. When a 3-mm grid size was employed, angular increments of 6° and 10° caused maximum dose variations of 3% (23 cGy) and 10% (61 cGy) in the spinal cord, respectively, while a 4° increment resulted in a dose difference of <1% (8 cGy) in all cases except for that of one patient. The 3-mm grid size and 4° angular increment enabled a 78% savings in computation time without making any critical sacrifices to dose accuracy.

CONCLUSIONS

A parameter set with a 3-mm grid size and a 4° angular increment is found to be appropriate for predicting patient dose distributions with a dose difference below 1% while reducing the computation time by more than half for lung SBRT using DCAT.

Geometrical splitting technique to improve the computational efficiency in Monte Carlo calculations for proton therapy

PURPOSE

To present the implementation and validation of a geometrical based variance reduction technique for the calculation of phase space data for proton therapy dose calculation.

METHODS

The treatment heads at the Francis H Burr Proton Therapy Center were modeled with a new Monte Carlo tool (TOPAS based on Geant4). For variance reduction purposes, two particle-splitting planes were implemented. First, the particles were split upstream of the second scatterer or at the second ionization chamber. Then, particles reaching another plane immediately upstream of the field specific aperture were split again. In each case, particles were split by a factor of 8. At the second ionization chamber and at the latter plane, the cylindrical symmetry of the proton beam was exploited to position the split particles at randomly spaced locations rotated around the beam axis. Phase space data in IAEA format were recorded at the treatment head exit and the computational efficiency was calculated. Depth-dose curves and beam profiles were analyzed. Dose distributions were compared for a voxelized water phantom for different treatment fields for both the reference and optimized simulations. In addition, dose in two patients was simulated with and without particle splitting to compare the efficiency and accuracy of the technique.

RESULTS

A normalized computational efficiency gain of a factor of 10-20.3 was reached for phase space calculations for the different treatment head options simulated. Depth-dose curves and beam profiles were in reasonable agreement with the simulation done without splitting: within 1% for depth-dose with an average difference of (0.2 ± 0.4)%, 1 standard deviation, and a 0.3% statistical uncertainty of the simulations in the high dose region; 1.6% for planar fluence with an average difference of (0.4 ± 0.5)% and a statistical uncertainty of 0.3% in the high fluence region. The percentage differences between dose distributions in water for simulations done with and without particle splitting were within the accepted clinical tolerance of 2%, with a 0.4% statistical uncertainty. For the two patient geometries considered, head and prostate, the efficiency gain was 20.9 and 14.7, respectively, with the percentages of voxels with gamma indices lower than unity 98.9% and 99.7%, respectively, using 2% and 2 mm criteria.

CONCLUSIONS

The authors have implemented an efficient variance reduction technique with significant speed improvements for proton Monte Carlo simulations. The method can be transferred to other codes and other treatment heads.

Comparison of dosimetric variation between prostate IMRT and VMAT due to patient's weight loss: Patient and phantom study

AIM

This study compared the dosimetric impact between prostate IMRT and VMAT due to patient's weight loss.

BACKGROUND

Dosimetric variation due to change of patient's body contour is difficult to predict in prostate IMRT and VMAT, since a large number of small and irregular segmental fields is used in the delivery.

MATERIALS AND METHODS

Five patients with prostate volumes ranging from 32.0 to 86.5 cm(3) and a heterogeneous pelvis phantom were used for prostate IMRT and VMAT plans using the same set of dose-volume constraints. Doses in IMRT and VMAT plans were recalculated with the patient's and phantom's body contour reduced by 0.5-2 cm to mimic size reduction. Dose coverage/criteria of the PTV and CTV and critical organs (rectum, bladder and femoral heads) were compared between IMRT and VMAT.

RESULTS

In IMRT plans, increases of the D99% for the PTV and CTV were equal to 4.0 ± 0.1% per cm of reduced depth, which were higher than those in VMAT plans (2.7 ± 0.24% per cm). Moreover, increases of the D30% of the rectum and bladder per reduced depth in IMRT plans (4.0 ± 0.2% per cm and 3.5 ± 0.5% per cm) were higher than those of VMAT (2.2 ± 0.2% per cm and 2.0 ± 0.6% per cm). This was also true for the increase of the D5% for the right femoral head in a patient or phantom with size reduction due to weight loss.

CONCLUSIONS

VMAT would be preferred to IMRT in prostate radiotherapy, when a patient has potential to suffer from weight loss during the treatment.

Variation of PTV dose distribution on patient size in prostate VMAT and IMRT: a dosimetric evaluation using the PTV dose–volume factor

Background

We propose to use the PTV dose–volume factor (PDVF) to evaluate treatment plans of prostate volumetric modulated arc therapy (VMAT) and intensity modulated radiotherapy (IMRT).

Purpose

PDVF was used to compare the variation of planning target volume (PTV) coverage between VMAT and IMRT because of weight loss of patient.

Materials and methods

VMAT and IMRT plans of five patients (prostate volume = 32–86·5 cm3) using the 6 MV photon beams were created with the external contour reduced by depths of 0·5–2 cm to reflect the weight loss. Moreover, integral doses (volume integral of the patient dose) and prostate tumour control probability (TCP) were calculated.

Results

We found that reduced depth resulted in PDVF decreasing 0·03 ± 4·7 × 10−4 (VMAT) and 0·04 ± 9·7 × 10−3 (IMRT) per cm for patients. The decrease of PDVF or degradation of PTV coverage was found more significant in IMRT plans than VMAT with patient size reduction. The integral dose did not change significantly between VMAT and IMRT, while the prostate TCP increased with an increase of reduced depth.

Conclusion

We concluded that PDVF can be successfully used to evaluate the variation of PTV coverage because of weight loss of patient in prostate VMAT and IMRT. Degradation of PTV coverage in prostate VMAT is less significant than IMRT regarding patient size reduction.

Prostate volumetric-modulated arc therapy: dosimetry and radiobiological model variation between the single-arc and double-arc technique

This study investigates the dosimetry and radiobiological model variation when a second photon arc was added to prostate volumetric‐modulated arc therapy (VMAT) using the single‐arc technique. Dosimetry and radiobiological model comparison between the single‐arc and double‐arc prostate VMAT plans were performed on five patients with prostate volumes ranging from 29−68.1 cm3. The prescription dose was 78 Gy/39 fractions and the photon beam energy was 6 MV. Dose‐volume histogram, mean and maximum dose of targets (planning and clinical target volume) and normal tissues (rectum, bladder and femoral heads), dose‐volume criteria in the treatment plan (D99% of PTV; D30%,D50%,V17Gy and V35Gy of rectum and bladder; D5% of femoral heads), and dose profiles along the vertical and horizontal axis crossing the isocenter were determined using the single‐arc and double‐arc VMAT technique. For comparison, the monitor unit based on the RapidArc delivery method, prostate tumor control probability (TCP), and rectal normal tissue complication probability (NTCP) based on the Lyman‐Burman‐Kutcher algorithm were calculated. It was found that though the double‐arc technique required almost double the treatment time than the single‐arc, the double‐arc plan provided a better rectal and bladder dose‐volume criteria by shifting the delivered dose in the patient from the anterior–posterior direction to the lateral. As the femoral head was less radiosensitive than the rectum and bladder, the double‐arc technique resulted in a prostate VMAT plan with better prostate coverage and rectal dose‐volume criteria compared to the single‐arc. The prostate TCP of the double‐arc plan was found slightly increased (0.16%) compared to the single‐arc. Therefore, when the rectal dose‐volume criteria are very difficult to achieve in a single‐arc prostate VMAT plan, it is worthwhile to consider the double‐arc technique.

PACS number: 87.55.D‐, 87.55.dk, 87.55.K‐, 87.55.Qr

Dosimetry estimation on variations of patient size in prostate volumetric-modulated arc therapy

This study investigated the dosimetric variations of the target and critical organs of patients who had weight loss associated with prostate volumetric-modulated arc therapy (VMAT). Five patients with prostate volumes ranging from 32-86.5 cm³ were selected from a group of 30 patients. Prostate VMAT plans were carried out on each patient using the 6-MV photon beam with a single 360° arc. Decrease of patient size as a result of weight loss was mimicked by contracting the patient's external contour in the anterior, left, and right directions with depths from 0.5-2 cm. Soft tissue excluded by the contracted external contour was replaced by air and the dose distribution was recalculated using the same beam geometry and dose prescription. Dose-volume histograms and dose-volume points such as D99% and D5% for the planning target volume (PTV), clinical target volume (CTV), rectum, bladder, and femoral heads were calculated with variations of reduced depth. In addition, the minimum, maximum, and mean doses for the target and critical organs were determined. PTV and CTV D99% were found to have increased 2.86 ± 0.30% per cm and 2.75 ± 0.38% per cm of reduced depth ranging from 0.5-2 cm. Moreover, the rectal and bladder D30% increased 2.20 ± 0.20% per cm and 2.31 ± 0.83% per cm, and the femoral head D5% increased 3.30 ± 0.11% per cm of reduced depth. Results from variations of the minimum, maximum, and mean doses of the PTV, CTV, rectum, bladder, and femoral heads showed that there was a >5% increase of dose when the reduced depth reached 2 cm. This study provided dosimetry estimation for radiation oncology staff to justify dose variations of the target and critical organs when patients' weight loss occurred in prostate VMAT. Dose variations >5% were seen when the patients' reduced depth was equal to 2 cm.

Searching standard parameters for volumetric modulated arc therapy (VMAT) of prostate cancer

Background

Since December 2009 a new VMAT planning system tool is available in Oncentra® MasterPlan v3.3 (Nucletron B.V.). The purpose of this study was to work out standard parameters for the optimization of prostate cancer.

Methods

For ten patients with localized prostate cancer plans for simultaneous integrated boost were optimized, varying systematically the number of arcs, collimator angle, the maximum delivery time, and the gantry spacing. Homogeneity in clinical target volume, minimum dose in planning target volume, median dose in the organs at risk, maximum dose in the posterior part of the rectum, and number of monitor units were evaluated using student’s test for statistical analysis. Measurements were performed with a 2D-array, taking the delivery time, and compared to the calculation by the gamma method.

Results

Plans with collimator 45° were superior to plans with collimator 0°. Single arc resulted in higher minimum dose in the planning target volume, but also higher dose values to the organs at risk, requiring less monitor units per fraction dose than dual arc. Single arc needs a higher value (per arc) for the maximum delivery time parameter than dual arc, but as only one arc is needed, the measured delivery time was shorter and stayed below 2.5 min versus 3 to 5 min. Balancing plan quality, dosimetric results and calculation time, a gantry spacing of 4° led to optimal results.

Conclusion

A set of parameters has been found which can be used as standard for volumetric modulated arc therapy planning of prostate cancer.

Implementing RapidArc into clinical routine: a comprehensive program from machine QA to TPS validation and patient QA

PURPOSE

With the increased commercial availability of intensity modulated arc therapy (IMAT) comes the need for comprehensive QA programs, covering the different aspects of this newly available technology. This manuscript proposes such a program for the RapidArc (RA) (Varian Medical Systems, Palo Alto) IMAT solution.

METHODS

The program was developed and tested out for a Millennium120 MLC on iX Clinacs and a HighDefinition MLC on a Novalis TX, using a variety of measurement equipment including Gafchromic film, 2D ion chamber arrays (Seven29 and StarCheck, PTW, Freiburg, Germany) with inclinometer and Octavius phantom, the Delta4 systam (ScandiDos, Uppsala, Sweden) and the portal imager (EPID). First, a number of complementary machine QA tests were developed to monitor the correct interplay between the accelerating/decelerating gantry, the variable dose rate and the MLC position, straining the delivery to the maximum allowed limits. Second, a systematic approach to the validation of the dose calculation for RA was adopted, starting with static gantry and RA specific static MLC shapes and gradually moving to dynamic gantry, dynamic MLC shapes. RA plans were then optimized on a series of artificial structures created within the homogeneous Octavius phantom and within a heterogeneous lung phantom. These served the double purpose of testing the behavior of the optimization algorithm (PRO) as well as the precision of the forward dose calculation. Finally, patient QA on a series of clinical cases was performed with different methods. In addition to the well established in-phantom QA, we evaluated the portal dosimetry solution within the Varian approach.

RESULTS

For routine machine QA, the "Snooker Cue" test on the EPID proved to be the most sensitive to overall problem detection. It is also the most practical one. The "Twinkle" and "Sunrise" tests were useful to obtain well differentiated information on the individual treatment delivery components. The AAA8.9 dose calculations showed excellent agreement with all corresponding measurements, except in areas where the 2.5 mm fixed fluence resolution was insufficient to accurately model the tongue and groove effect or the dose through nearly closed opposing leafs. Such cases benefited from the increased fluence resolution in AAA10.0. In the clinical RA fields, these effects were smeared out spatially and the impact of the fluence resolution was considerably less pronounced. The RA plans on the artificial structure sets demonstrated some interesting characteristics of the PRO8.9 optimizer, such as a sometimes unexpected dependence on the collimator rotation and a suboptimal coverage of targets within lung tissue. Although the portal dosimetry was successfully validated, we are reluctant to use it as a sole means of patient QA as long as no gantry angle information is embedded.

CONCLUSIONS

The all-in validation program allows a systematic approach in monitoring the different levels of RA treatments. With the systematic approach comes a better understanding of both the capabilities and the limits of the used solution. The program can be useful for implementation, but also for the validation of major upgrades.

Dynamic tomotherapy delivery

PURPOSE

Several dynamic techniques are introduced to speed up TomoTherapy delivery and improve longitudinal conformity. These techniques include dynamic jaws, dynamic couch, and their combinations.

METHODS

In general, dynamic jaws techniques allow jaws to move during a treatment. On the one hand, the jaws open wide to increase efficiency and thus reduce beam-on time. On the other hand, the jaws can close and follow the target border when sharp penumbra is required near the superior and inferior borders of tumor sites, which results in improved longitudinal dose conformity. The main purpose of the dynamic couch technique is to move the couch as fast as possible at variable speed to reduce beam-on time. Delivering most conformal dose as fast as possible requires a combination of dynamic jaws and dynamic couch techniques (DJDC). Motions of the jaws and couch are determined from the longitudinal fluence profile, which is calculated from an optimized leaf sinogram of small jaw width regular delivery or running start and stop delivery (RSS). We focused on RSS and DJDC in this study and also discussed other delivery techniques.

RESULTS

Several conceptual cases are simulated to compare different delivery techniques. The results show that beam-on time can be reduced by about 60% compared to regular delivery with a 2.5 cm jaw width (REG 2.5 cm) for these cases and arbitrary longitudinal fluence profiles can be delivered. Two clinical cases, a prostate and a head-and-neck case, with different delivery techniques are calculated. The results show that plan quality yielded by DJDC with a maximum 5.0 cm jaw width is overall comparable to or better than that of the existing REG 2.5 cm.

CONCLUSIONS

The DVH comparisons show better critical structure avoidance with the dynamic techniques. At the same time, beam-on time is reduced by about one half compared to REG 2.5 cm. Dynamic delivery techniques provide users more tools to speed up delivery and/or improve plan quality.

Curcumin decreases survival of Hep3B liver and MCF-7 breast cancer cells: the role of HIF

BACKGROUND

Curcumin, a commonly used spice, affects the activities of cytokines, enzymes, and transcription factors that are linked to inflammation. Furthermore, curcumin has been assigned tumor growth inhibiting effects, possibly mediated by promoting hypoxia-inducible factor (HIF) degradation. HIFs are transcription factors that play a central role in the adaptation and response to low oxygen levels in metazoan cells. However, curcumin also exhibits properties of an iron chelator indicating its potential of inhibiting HIF-α prolyl hydroxylase (PHD) activity.

METHODS

We were interested in clarifying these divergent actions of curcumin in due consideration of the effects on radio-therapy. Thus, concentration- and time-dependent effects of curcumin on HIF-α and -β protein levels and activity in hepatoma and breast carcinoma cell cultures under normoxic and hypoxic conditions were studied.

RESULTS

It was shown that HIF-1α accumulated in normoxia after the application of higher doses of the drug. Curcumin proved to lower HIF-1α and HIF-2α protein levels in hypoxia. HIF-1β (ARNT; arylhydrocarbon nuclear translocator) protein levels and HIF transcriptional activity were reduced in normoxia and hypoxia after 4 h and 24 h incubation periods. Furthermore, curcumin treatment negatively impacted on clonogenic cell survival of Hep3B hepatoma and MCF-7 breast carcinoma cells.

CONCLUSION

Effects of curcumin on cell growth and survival factor expression suggest its potential benefit in the treatment of cancer without a direct radiosensitizing influence of curcumin on these cells.

Intensity-modulated arc therapy: principles, technologies and clinical implementation

Intensity-modulated arc therapy (IMAT) was proposed by Yu (1995 Phys. Med. Biol. 40 1435-49) as an alternative to tomotherapy. Over more than a decade, much progress has been made. The advantages and limitations of the IMAT technique have also been better understood. In recent years, single-arc forms of IMAT have emerged and become commercially adopted. The leading example is the volumetric-modulated arc therapy (VMAT), a single-arc form of IMAT that delivers apertures of varying weights with a single-arc rotation that uses dose-rate variation of the treatment machine. With commercial implementation of VMAT, wide clinical adoption has quickly taken root. However, there remains a lack of general understanding for the planning of such arc treatments, as well as what delivery limitations and compromises are made. Commercial promotion and competition add further confusion for the end users. It is therefore necessary to provide a summary of this technology and some guidelines on its clinical implementation. The purpose of this review is to provide a summary of the works from the radiotherapy community that led to wide clinical adoption, and point out the issues that still remain, providing some perspective on its further developments. Because there has been vast experience in IMRT using multiple intensity-modulated fields, comparisons between IMAT and IMRT are also made in the review within the areas of planning, delivery and quality assurance.

Evaluation of a new VMAT QA device, or the "X" and "O" array geometries

We introduce a logical process of three distinct phases to begin the evaluation of a new 3D dosimetry array. The array under investigation is a hollow cylinder phantom with diode detectors fixed in a helical shell forming an “O” axial detector cross section (ArcCHECK), with comparisons drawn to a previously studied 3D array with diodes fixed in two crossing planes forming an “X” axial cross section (Delta⁴). Phase I testing of the ArcCHECK establishes: robust relative calibration (response equalization) of the individual detectors, minor field size dependency of response not present in a 2D predecessor, and uncorrected angular response dependence in the axial plane. Phase II testing reveals vast differences between the two devices when studying fixed‐width full circle arcs. These differences are primarily due to arc discretization by the TPS that produces low passing rates for the peripheral detectors of the ArcCHECK, but high passing rates for the Delta⁴. Similar, although less pronounced, effects are seen for the test VMAT plans modeled after the AAPM TG119 report. The very different 3D detector locations of the two devices, along with the knock‐on effect of different percent normalization strategies, prove that the analysis results from the devices are distinct and noninterchangeable; they are truly measuring different things. The value of what each device measures, namely their correlation with – or ability to predict – clinically relevant errors in calculation and/or delivery of dose is the subject of future Phase III work.

PACS number: 87.55Qr

Validation of Pinnacle treatment planning system for use with Novalis delivery unit

For an institution that already owns the licenses, it is economically advantageous and technically feasible to use Pinnacle TPS (Philips Radiation Oncology Systems, Fitchburg, WI) with the BrainLab Novalis delivery system (BrainLAB A.G., Heimstetten, Germany). This takes advantage of the improved accuracy of the convolution algorithm in the presence of heterogeneities compared with the pencil beam calculation, which is particularly significant for lung SBRT treatments. The reference patient positioning DRRs still have to be generated by the BrainLab software from the CT images and isocenter coordinates transferred from Pinnacle. We validated this process with the end‐to‐end hidden target test, which showed an isocenter positioning error within one standard deviation from the previously established mean value. The Novalis treatment table attenuation is substantial (up to 6.2% for a beam directed straight up and up to 8.4% for oblique incidence) and has to be accounted for in calculations. A simple single‐contour treatment table model was developed, resulting in mean differences between the measured and calculated attenuation factors of 0.0%–0.2%, depending on the field size. The maximum difference for a single incidence angle is 1.1%. The BrainLab micro‐MLC (mMLC) leaf tip, although not geometrically round, can be represented in Pinnacle by an arch with satisfactory dosimetric accuracy. Subsequently, step‐and‐shoot (direct machine parameter optimization) IMRT dosimetric agreement is excellent. VMAT (called “SmartArc” in Pinnacle) treatments with constant gantry speed and dose rate are feasible without any modifications to the accelerator. Due to the 3 mm‐wide mMLC leaves, the use of a 2 mm calculation grid is recommended. When dual arcs are used for the more complex cases, the overall dosimetric agreement for the SmartArc plans compares favorably with the previously reported results for other implementations of VMAT: γ(3%,3mm) for absolute dose obtained with the biplanar diode array passing rates above 97% with the mean of 98.6%. However, a larger than expected dose error with the single‐arc plans, confined predominantly to the isocenter region, requires further investigation

PACS numbers: 87.55Qr, 87.56Nk

The number of beams in IMRT--theoretical investigations and implications for single-arc IMRT

The first purpose of this paper is to shed some new light on the old question of selecting the number of beams in intensity-modulated radiation therapy (IMRT). The second purpose is to illuminate the related issue of discrete static beam angles versus rotational techniques, which has recently re-surfaced due to the advancement of volumetric modulated arc therapy (VMAT). A specific objective is to find analytical expressions that allow one to address the points raised above. To make the problem mathematically tractable, it is assumed that the depth dose is flat and that the lateral dose profile can be approximated by polynomials, specifically Chebyshev polynomials of the first kind, of finite degree. The application of methods known from image reconstruction then allows one to answer the first question above as follows: the required number of beams is determined by the maximum degree of the polynomials used in the approximation of the beam profiles, which is a measure of the dose variability. There is nothing to be gained by using more beams. In realistic cases, in which the variability of the lateral dose profile is restricted in several ways, the required number of beams is of the order of 10-20. The consequence of delivering the beams with a 'leaf sweep' technique during continuous rotation of the gantry, as in VMAT, is also derived in an analytical form. The main effect is that the beams fan out, but the effect near the axis of rotation is small. This result can serve as a theoretical justification of VMAT. Overall the analytical derivations in this paper, albeit based on strong simplifications, provide new insights into, and a deeper understanding of, the beam angle problem in IMRT. The decomposition of the beam profiles into well-behaved and easily deliverable smooth functions, such as Chebyshev polynomials, could be of general interest in IMRT treatment planning.

Initial dosimetric evaluation of SmartArc - a novel VMAT treatment planning module implemented in a multi-vendor delivery chain

We performed an initial dosimetric evaluation of SmartArc – a novel VMAT planning module for the Philips Pinnacle treatment planning system. It was implemented in a multi‐vendor environment, with the other two major components of the delivery chain being MOSAIQ record and verify system (IMPAC Medical Systems, Sunnyvale, CA) and a Trilogy linac (Varian Medical Systems, Palo Alto, CA). A test suite of structure sets and dose objectives provided by the AAPM for multi‐institutional comparison of IMRT dosimetry was used. A total of fifty plans were successfully delivered. The effect of control point spacing on dosimetric accuracy was investigated. When calculated with the 4° spacing, the overall mean point dose errors measured with an ion chamber were 0.5±1.4and −0.3±1.4% for the PTV and OAR, respectively. The γ(3%, 3 mm) passing rate, measured for absolute dose with a biplanar diode array, was 98.2±1.6% (range 94.5–99.9%). Ninety percent of the passing rate values were above 97.7%. With the 6° control point spacing, the highly modulated plans exhibited large dosimetric errors (e.g. γ(3%, 3 mm) passing rates below 90% and ion chamber point dose errors of 6–12%), while the results were still acceptable for the simpler cases. The data show that the practical accuracy of the small‐arc approximation, which is at the heart of VMAT dose calculations, depends not only on the control point spacing, but also on the size and relative position of the MLC openings corresponding to the consecutive control points. The effect of the minimum allowed separation between the opposing leaves was found to be minimal. It appears that 4° control point spacing may be a good compromise between calculation speed and accuracy. However each institution is encouraged to establish its own treatment planning guidelines based on the case complexity and acceptable error level.

PACS number: 87.55Qr