Letter in response to "Optimization of radiation therapy and the development of multileaf collimation" by Anders Brahme

A practical Monte Carlo MU verification tool for IMRT quality assurance

Quality assurance (QA) for intensity-modulated radiation therapy (IMRT) treatment planning and beam delivery, using ionization chamber measurements and film dosimetry in a phantom, is time consuming. The Monte Carlo method is the most accurate method for radiotherapy dose calculation. However, a major drawback of Monte Carlo dose calculation as currently implemented is its slow speed. The goal of this work is to bring the efficiency of Monte Carlo into a practical range by developing a fast Monte Carlo monitor unit (MU) verification tool for IMRT. A special estimator for dose at a point called the point detector has been used in this research. The point detector uses the next event estimation (NEE) method to calculate the photon energy fluence at a point of interest and then converts it to collision kerma by the mass energy absorption coefficient assuming the presence of transient charged particle equilibrium. The MU verification tool has been validated by comparing the calculation results with measurements. It can be used for both patient dose verification and phantom QA calculation. The dynamic leaf-sequence log file is used to rebuild the actual MLC leaf sequence in order to predict the dose actually received by the patient. Dose calculations for 20 patient plans have been performed using the point detector method. Results were compared with direct Monte Carlo simulations using EGS4/MCSIM, which is a well-benchmarked Monte Carlo code. The results between the point detector and MCSIM agreed to within 2%. A factor of 20 speedup can be achieved with the point detector method compared with direct Monte Carlo simulations.

Intensity modulation in radiotherapy: photons versus protons in the paranasal sinus

PURPOSE

The purpose of this study is to investigate whether successive tightening of normal tissue constraints on an intensity modulated X-ray therapy plan might be able to improve it to the point of clinical comparability with the corresponding intensity modulated proton therapy plan.

MATERIALS AND METHODS

Photon and proton intensity modulated plans were calculated for a paranasal sinus case using nominal dose constraints. Additional photon plans were then calculated in an effort to match the dose-volume histograms of the critical structures to those of the proton plan.

RESULTS

On reducing the low dose contribution to both orbits in the photon plan by tightening the constraints on these structures, an increased dose heterogeneity across the target resulted. When all critical structures were more strictly constrained, target dose homogeneity and conformity was further compromised. An increased integral dose to the non-critical normal tissues was observed for the photon plans as dose was progressively removed from the critical structures.

CONCLUSIONS

Both modalities were found to provide comparable target volume conformation and sparing of critical structures, when the nominal dose constraints were applied. However, the use of intensity modulated protons provided the only method by which critical structures could be spared at all dose levels, whilst simultaneously providing acceptable dose homogeneity within the target volume.

Clinical implementation of intensity-modulated arc therapy

PURPOSE

Intensity-modulated arc therapy (IMAT) is a method for delivering intensity-modulated radiation therapy (IMRT) using rotational beams. During delivery, the field shape, formed by a multileaf collimator (MLC), changes constantly. The objectives of this study were to (1) clinically implement the IMAT technique, and (2) evaluate the dosimetry in comparison with conventional three-dimensional (3D) conformal techniques.

METHODS AND MATERIALS

Forward planning with a commercial system (RenderPlan 3D, Precision Therapy International, Inc., Norcross, GA) was used for IMAT planning. Arcs were approximated as multiple shaped fields spaced every 5-10 degrees around the patient. The number and ranges of the arcs were chosen manually. Multiple coplanar, superimposing arcs or noncoplanar arcs with or without a wedge were allowed. For comparison, conventional 3D conformal treatment plans were generated with the same commercial forward planning system as for IMAT. Intensity-modulated treatment plans were also created with a commercial inverse planning system (CORVUS, Nomos Corporation). A leaf-sequencing program was developed to generate the dynamic MLC prescriptions. IMAT treatment delivery was accomplished by programming the linear accelerator (linac) to deliver an arc and the MLC to step through a sequence of fields. Both gantry rotation and leaf motion were enslaved to the delivered MUs. Dosimetric accuracy of the entire process was verified with phantoms before IMAT was used clinically. For each IMAT treatment, a dry run was performed to assess the geometric and dosimetric accuracy. Both the central axis dose and dose distributions were measured and compared with predictions by the planning system.

RESULTS

By the end of May 2001, 50 patients had completed their treatments with the IMAT technique. Two to five arcs were needed to achieve highly conformal dose distributions. The IMAT plans provided better dose uniformity in the target and lower doses to normal structures than 3D conformal plans. The results varied when the comparison was made with fixed gantry IMRT. In general, IMAT plans provided more uniform dose distributions in the target, whereas the inverse-planned fixed gantry treatments had greater flexibility in controlling dose to the critical structures. Because the field sizes and shapes used in the IMAT were similar to those used in conventional treatments, the dosimetric uncertainty was very small. Of the first 32 patients treated, the average difference between the measured and predicted doses was -0.54 +/- 1.72% at isocenter. The 80%-95% isodose contours measured with film dosimetry matched those predicted by the planning system to within 2 mm. The planning time for IMAT was slightly longer than for generating conventional 3D conformal plans. However, because of the need to create phantom plans for the dry run, the overall planning time was doubled. The average time a patient spent on the table for IMAT treatment was similar to conventional treatments.

CONCLUSION

Initial results demonstrated the feasibility and accuracy of IMAT for achieving highly conformal dose distributions for different sites. If treatment plans can be optimized for IMAT cone beam delivery, we expect IMAT to achieve dose distributions that rival both slice-based and fixed-field IMRT techniques. The efficient delivery with existing linac and MLC makes IMAT a practical choice.

A simple model for examining issues in radiotherapy optimization

Convolution/superposition software has been used to produce a library of photon pencil beam dose matrices. This library of pencil beams is designed to serve as a tool for both education and investigation in the field of radiotherapy optimization. The elegance of this pencil beam model stems from its cylindrical symmetry. Because of the symmetry, the dose distribution for a pencil beam from any arbitrary angle can be determined through a simple rotation of a pre-computed dose matrix. Rapid dose calculations can thus be performed while maintaining the accuracy of a convolution/superposition based dose computation. The pencil beam data sets have been made publicly available. It is hoped that the data sets will facilitate a comparison of a variety of optimization and delivery approaches. This paper will present a number of studies designed to demonstrate the usefulness of the pencil beam data sets. These studies include an examination of the extent to which a treatment plan can be improved through either an increase in the number of beam angles and/or a decrease in the collimator size. A few insights into the significance of heterogeneity corrections for treatment planning for intensity modulated radiotherapy will also be presented.

Optimizing the planning of intensity-modulated radiotherapy

A method of computing optimized intensity-modulated beam profiles has been further developed and used to generate highly conformal radiotherapy dose distributions. The features of these distributions are shown to be strongly dependent on the tuning built into the algorithm. The optimization aims to achieve a specified dose prescription with a posteriori computation of probabilistic biological response. A method has been developed to show the effect of stratifying the intensity-modulated beam profiles into a number of finite intensity increments. It is shown that, provided the number of intensity strata is not too small, highly conformal dose distributions can be achieved with a number of fields (e.g. 7 or 9) which is not excessively large. This number, however, depends on the exact shape of the planning target volume (PTV and its disposition with respect to juxtaposed organs at risk (OARs). These intensity-modulated profiles can therefore be delivered either by apparatus for 'tomotherapy' or by using the multileaf collimator at each gantry orientation to deliver a sequence of fixed fields with different field sizes, constructing the beam profile via finite increments of beam intensity. When the PTV and OARS overlap, due to including a finite margin on the clinical target volume to account for tissue movement, it is shown that the dose delivered to the overlap region provides a limit on what can be achieved with conformal therapy. This problem is encountered, for example, when treating the prostate which lies next to part of the rectum and bladder. Some comment is provided on, but not a solution for, the problem of optimizing field orientation.

Optimized planning using physical objectives and constraints

Intensity-modulated radiation therapy (IMRT) allows one to achieve a better conformation of the high-dose region to the prescribed tumor target volume than uniform beam therapy, especially in complex treatment situations. Still, perfect conformation is impossible. Hence the goal of optimized IMRT planning or inverse planning is to find the beam profiles that yield the optimum among the physically achievable treatment plans. The principal physical advantage of IMRT is best exploited if the optimization is driven by physical criteria. This article presents an overview of such physical, yet clinically relevant, criteria along with optimization algorithms that take these criteria into account. Practical computer implementations are described, which allow one to perform the optimization in an interactive manner within a few minutes. The application of these methods to some complex clinical example cases is presented, and the results are compared with uniform beam treatment plans and with biologically optimized plans.

Beam profiles for x-ray rotation therapy

Determining the beam configuration necessary to deliver a desired dose distribution with rotation therapy is equivalent to solving an integral equation. The equation has been solved analytically for a handful of dose distributions having specific radial variation and either rotational or reflective angular symmetry. In this work a numerical method for calculating beam profiles appropriate for producing distributions having arbitrary radial variation and angular symmetry of order l > or = 2 is presented. The accuracy of the technique is demonstrated by comparison with one of the few dose distributions for which an analytic solution exists, and the ability to produce both more general and conformal distributions is also shown. The problems of negative intensity and scatter are discussed.

Number and orientations of beams in intensity-modulated radiation treatments

The fundamental question of how many equispaced coplanar intensity-modulated photon beams are required to obtain an optimum treatment plan is investigated in a dose escalation study for a typical prostate tumor. Furthermore, optimization of beam orientations to improve dose distributions is explored. A dose-based objective function and a fast gradient technique are employed for optimizing the intensity profiles (inverse planning). An exhaustive search and fast simulated annealing techniques (FSA) are used to optimize beam orientations. However, to keep computation times reasonable, the intensity profiles for each beam arrangement are still optimized using inverse planning. A pencil beam convolution algorithm is employed for dose calculation. All calculations are performed in three-dimensional (3D) geometry for 15 MV photons. DVHs, dose displays, TCP, NTCP, and biological score functions are used for evaluation of treatment plans. It is shown that for the prostate case presented here, the minimum required number of equiangular beams depends on the prescription dose level and ranges from three beams for 70 Gy plans to seven to nine beams for 81 Gy plans. For the highest dose level (81 Gy), beam orientations are optimized and compared to equiangular spaced arrangements. It is shown that (1) optimizing beam orientations is most valuable for a small numbers of beams (< or = 5) and the gain diminishes rapidly for higher numbers of beams; (2) if sensitive structures (for example rectum) are partially enclosed by the target volume, beams coming from their direction tend to be preferable, since they allow greater control over dose distributions; (3) while FSA and an exhaustive search lead to the same results, computation times using FSA are reduced by two orders of magnitude to clinically acceptable values. Moreover, characteristics of and demands on biology-based and dose-based objective functions for optimization of intensity-modulated treatments are discussed.

Evaluation of dose calculation algorithm of the peacock system for multileaf intensity modulation collimator

PURPOSE

To evaluate the dose calculation algorithm used in the inverse treatment planning computer system for the intensity modulation multileaf collimator.

METHODS AND MATERIALS

The inverse treatment-planning computer system calculates the intensities of multiple pencil beams to achieve an optimal distribution and modulates the beam intensity through the special multileaf collimator. The system's dose calculation algorithm made the two basic assumptions: (a) The tissue-maximum ratios (TMRs) of a single pencil beam have the same values as TMRs for raylines through each pencil beam that are determined from percentage depth dose isodose curves along the long axis of the 2 x 20 cm2 field with all leaves open; and (b) the relative output factors (ROF) of each pencil beam also have the same values as the rayline TMR at d(max) of the 2 x 20 cm2 field. To verify these two assumptions, a special multileaf collimator was installed to our linear accelerator which produces 4 MV x-rays. The TMRs and ROFs for the single leaves 1 through 10 were measured using an ion chamber and TLD dosimeter in either a water or a polystyrene phantom. The values of rayline TMRs were calculated from the measured crossplane isodose curves of the 2 x 20 cm2 field. Comparisons were made between these two sets of data.

RESULTS

Based on our measurements, we found that the ROFs of a pencil beam obtained from the rayline TMRs at d(max) are as much as 7.6% greater than that of single pencil beams. The ROF of the 1 x 1 cm2 pencil beam is 4 and 6.5% less than that of a cluster of four neighboring pencil beams forming a 2 x 2 cm2, and a 2 x 20 cm2 field respectively. However, the rayline TMRs are generally larger than the TMRs of a single pencil beam. At a depth of 8 cm, the average depth in the middle of intracranial space, the rayline TMRs of the pencil beams of leaves 1 and 10 are 5.4 and 9% higher than a single pencil beam TMR at the same depth, respectively. Also interesting is to note that the TMRs of each of the single pencil beams were found to be equal.

CONCLUSIONS

In our article, evaluations and comparisons of TMRs and ROFs were made for two extreme conditions. The measured values of TMRs and ROFs of a single beam have been shown to be significantly different from those used in the calculations. Because both the TMR and ROF are influenced by the scattering radiation in the same direction, the deviations for these two factors would be expected to be magnified. Thus, for the two extreme situations we have investigated, dose deviations would be on the order of 15%. In real patient treatment; of course, these deviations may be somewhat less, but still significant. Our results, however, show that further investigations are warranted.

Optimization of intensity-modulated 3D conformal treatment plans based on biological indices

To overcome the limitations of the intensity modulation optimization techniques based on dose criteria, we introduce a method for optimizing intensity distributions in which we employ an objective function based on biological indices. The objective function also includes constraints on dose and dose-volume combinations to ensure that the results are consistent with the physician's judgement. We apply a variant of the steepest-descent method to optimize the objective function. The method is three-dimensional and incorporates scattered radiation in the optimization process using an iterative scheme employing the pencil beam convolution method. Previously we had shown that the inverse technique of obtaining optimum intensity distributions, for which the objectives are defined in terms of a desired uniform dose to the target volume and desired upper limits of dose to normal organs, produces satisfactory approximations of the desired dose distributions for prostate plans. However, for lung, the performance of this technique was considerably inferior. Our conclusion was that, in general, it is not sufficient to specify the objectives of optimization purely in terms of a desired pattern of dose and that the objectives should also incorporate biology, perhaps in the form of biological indices. We demonstrate that the biology-based approach produces lung plans that are superior to those produced when only dose-based objectives are used. For the treatment of prostate, the two methods produce comparable dose distributions.

Intensity-modulated arc therapy with dynamic multileaf collimation: an alternative to tomotherapy

The desire to improve local tumour control and cure more cancer patients, coupled with advances in computer technology and linear accelerator design, has spurred the developments of three-dimensional conformal radiotherapy techniques. Optimized treatment plans, aiming to deliver high dose to the target while minimizing dose to the surrounding tissues, can be delivered with multiple fields each with spatially modulated beam intensities or with multiple-slice treatments. This paper introduces a new method, intensity-modulated arc therapy (IMAT), for delivering optimized treatment plans to improve the therapeutic ratio. It utilizes continuous gantry motion as in conventional arc therapy. Unlike conventional arc therapy, the field shape, which is conformed with the multileaf collimator, changes during gantry rotation. Arbitrary two-dimensional beam intensify distributions at different beam angles are delivered with multiple superimposing arcs. A system capable of delivering the IMAT has been implemented. An example is given that illustrates the feasibility of this new method. Advantages of this new technique over tomotherapy and other slice-based delivery schemes are also discussed.

External beam radiotherapy of localized prostatic adenocarcinoma. Evaluation of conformal therapy, field number and target margins

The purpose of the present study was to identify factors of importance in the planning of external beam radiotherapy of prostatic adenocarcinoma. Seven patients with urogenital cancers were planned for external radiotherapy of the prostate. Four different techniques were used, viz. a 4-field box technique and four-, five- or six-field conformal therapy set-ups combined with three different margins (1-3 cm). The evaluations were based on the doses delivered to the rectum and the urinary bladder. A normal tissue complication probability (NTCP) was calculated for each plan using Lyman's dose volume reduction method. The most important factors that resulted in a decrease of the dose delivered to the rectum and the bladder were the use of conformal therapy and smaller margins. Conformal therapy seemed more important for the dose distribution in the urinary bladder. Five- and six-field set-ups were not significantly better than those with four fields. NTCP calculations were in accordance with the evaluation of the dose volume histograms. To conclude, four-field conformal therapy utilizing reduced margins improves the dose distribution to the rectum and the urinary bladder in the radiotherapy of prostatic adenocarcinoma.