Errors in off-axis treatment planning for a 4-MeV machine

Computer treatment planning systems allow dose computation in planes parallel to the central one (off-axis plans). The beam data may consist of, e.g., percentage depth doses along the central axis plus off-axis ratios (OAR) at several depths. In some systems, the calculation of an off-axis plan is based on the assumption that the OAR can be represented by a separable function: OAR(x,y) = f(x).g(y), where x and y are the symmetry axes perpendicular to the beam axis and the functions f and g are equal for a square open field. The errors of this assumption for a 4-MeV machine were measured for open fields and wedged fields at five different depths. The measured dose was compared with that predicted by the above equation for 50%, 75%, and 88% of the half field width from the beam axis. Maximum deviation of more than 10% was observed with the probable error of the measurement being 1%.     

Evaluation of a commercial three-dimensional electron beam treatment planning system

We evaluated a commercial three-dimensional (3D) electron beam treatment planning system (CADPLAN V.2.7.9) using both experimentally measured and Monte Carlo calculated dose distributions to compare with those predicted by CADPLAN calculations. Tests were carried out at various field sizes and electron beam energies from 6 to 20 MeV. For a homogeneous water phantom the agreement between measured and CADPLAN calculated dose distributions is very good except at the phantom surface. CADPLAN is able to predict hot and cold spots caused by a simple 3D inhomogeneity but unable to predict dose distributions for a more complex geometry where CADPLAN underestimates dose changes caused by inhomogeneity. We discussed possible causes for the inaccuracy in the CADPLAN dose calculations. In addition, we have tested CADPLAN treatment monitor unit and electron cut-out factor calculations and found that CADPLAN predictions generally agree with manual calculations.     

Margins for treatment planning of proton therapy

For protons and other charged particles, the effect of set-up errors on the position of isodoses is considerably less in the direction of the incident beam than it is laterally. Therefore, the margins required between the clinical target volume (CTV) and planning target volume (PTV) can be less in the direction of the incident beam than laterally. Margins have been calculated for a typical head plan and a typical prostate plan, for a single field, a parallel opposed and a four-field arrangement of protons, and compared with margins calculated for photons, assuming identical geometrical uncertainties for each modality. In the head plan, where internal motion was assumed negligible, the CTV-PTV margin reduced from approximately 10 mm to 3 mm in the axial direction for the single field and parallel opposed plans. For a prostate plan, where internal motion cannot be ignored, the corresponding reduction in margin was from 11 mm to 7 mm. The planning organ at risk (PRV) margin in the axial direction reduced from 6 mm to 2 mm for the head plan, and from 7 mm to 4 mm for the prostate plan. No reduction was seen on the other axes, or for any axis of the four-field plans. Owing to the shape of proton dose distributions, there are many clinical cases in which good dose distributions can be obtained with one or two fields. When this is done, it is possible to use smaller PTV and PRV margins. This has the potential to convert untreatable cases, in which the PTV and PRV overlap, into cases with a gap between PTV and PRV of adequate size for treatment planning.     

Monte Carlo evaluation of a photon pencil kernel algorithm applied to fast neutron therapy treatment planning

When dedicated software is lacking, treatment planning for fast neutron therapy is sometimes performed using dose calculation algorithms designed for photon beam therapy. In this work Monte Carlo derived neutron pencil kernels in water were parametrized using the photon dose algorithm implemented in the Nucletron TMS (treatment management system) treatment planning system. A rectangular fast-neutron fluence spectrum with energies 0-40 MeV (resembling a polyethylene filtered p(41)+Be spectrum) was used. Central axis depth doses and lateral dose distributions were calculated and compared with the corresponding dose distributions from Monte Carlo calculations for homogeneous water and heterogeneous slab phantoms. All absorbed doses were normalized to the reference dose at 10 cm depth for a field of radius 5.6 cm in a 30 x 40 x 20 cm3 water test phantom. Agreement to within 7% was found in both the lateral and the depth dose distributions. The deviations could be explained as due to differences in size between the test phantom and that used in deriving the pencil kernel (radius 200 cm, thickness 50 cm). In the heterogeneous phantom, the TMS, with a directly applied neutron pencil kernel, and Monte Carlo calculated absorbed doses agree approximately for muscle but show large deviations for media such as adipose or bone. For the latter media, agreement was substantially improved by correcting the absorbed doses calculated in TMS with the neutron kerma factor ratio and the stopping power ratio between tissue and water. The multipurpose Monte Carlo code FLUKA was used both in calculating the pencil kernel and in direct calculations of absorbed dose in the phantom.     

Clinical use of a commercial Monte Carlo treatment planning system for electron beams

In 2002 we fully implemented clinically a commercial Monte Carlo based treatment planning system for electron beams. The software, developed by MDS Nordion (presently Nucletron), is based on Kawrakow's VMC++ algorithm. The Monte Carlo module is integrated with our Theraplan Plustrade mark treatment planning system. An extensive commissioning process preceded clinical implementation of this software. Using a single virtual 'machine' for each electron beam energy, we can now calculate very accurately the dose distributions and the number of MU for any arbitrary field shape and SSD. This new treatment planning capability has significantly impacted our clinical practice. Since we are more confident of the actual dose delivered to a patient, we now calculate accurate three-dimensional (3D) dose distributions for a greater variety of techniques and anatomical sites than we have in the past. We use the Monte Carlo module to calculate dose for head and neck, breast, chest wall and abdominal treatments with electron beams applied either solo or in conjunction with photons. In some cases patient treatment decisions have been changed, as compared to how such patients would have been treated in the past. In this paper, we present the planning procedure and some clinical examples.     

Dosimetric evaluation in heterogeneous tissue of anterior electron beam irradiation for treatment of retinoblastoma

A dosimetric study of anterior electron beam irradiation for treatment of retinoblastoma was performed to evaluate the influence of tissue heterogeneities on the dose distribution within the eye and the accuracy of the dose calculated by a pencil beam algorithm. Film measurements were made in a variety of polystyrene phantoms and in a removable polystyrene eye incorporated into a tissue substitute phantom constructed from a human skull. Measurements in polystyrene phantoms were used to demonstrate the algorithm's ability to predict the effect of a lens block placed in the beam, as well as the eye's irregular surface shape. The eye phantom was used to measure dose distributions within the eye in both the sagittal and transverse planes in order to test the algorithm's ability to predict the dose distribution when bony heterogeneities are present. Results show (1) that previous treatment planning conclusions based on flat, uniform phantoms for central-axis depth dose are adequate; (2) that a three-dimensional heterogeneity correction is required for accurate dose calculations; and (3) that if only a two-dimensional heterogeneity correction is used in calculating the dose, it is more accurate for the sagittal than the transverse plane.     

Comparison of anatomy-based, fluence-based and aperture-based treatment planning approaches for VMAT

Volumetric modulated arc therapy (VMAT) has the potential to reduce treatment times while producing comparable or improved dose distributions relative to fixed-field intensity-modulated radiation therapy. In order to take full advantage of the VMAT delivery technique, one must select a robust inverse planning tool. The purpose of this study was to evaluate the effectiveness and efficiency of VMAT planning techniques of three categories: anatomy-based, fluence-based and aperture-based inverse planning. We have compared these techniques in terms of the plan quality, planning efficiency and delivery efficiency. Fourteen patients were selected for this study including six head-and-neck (HN) cases, and two cases each of prostate, pancreas, lung and partial brain. For each case, three VMAT plans were created. The first VMAT plan was generated based on the anatomical geometry. In the Elekta ERGO++ treatment planning system (TPS), segments were generated based on the beam's eye view (BEV) of the target and the organs at risk. The segment shapes were then exported to Pinnacle TPS followed by segment weight optimization and final dose calculation. The second VMAT plan was generated by converting optimized fluence maps (calculated by the Pinnacle TPS) into deliverable arcs using an in-house arc sequencer. The third VMAT plan was generated using the Pinnacle SmartArc IMRT module which is an aperture-based optimization method. All VMAT plans were delivered using an Elekta Synergy linear accelerator and the plan comparisons were made in terms of plan quality and delivery efficiency. The results show that for cases of little or modest complexity such as prostate, pancreas, lung and brain, the anatomy-based approach provides similar target coverage and critical structure sparing, but less conformal dose distributions as compared to the other two approaches. For more complex HN cases, the anatomy-based approach is not able to provide clinically acceptable VMAT plans while highly conformal dose distributions were obtained using both aperture-based and fluence-based inverse planning techniques. The aperture-based approach provides improved dose conformity than the fluence-based technique in complex cases.     

Monte Carlo based treatment planning for modulated electron beam radiation therapy

A Monte Carlo based treatment planning system for modulated electron radiation therapy (MERT) is presented. This new variation of intensity modulated radiation therapy (IMRT) utilizes an electron multileaf collimator (eMLC) to deliver non-uniform intensity maps at several electron energies. In this way, conformal dose distributions are delivered to irregular targets located a few centimetres below the surface while sparing deeper-lying normal anatomy. Planning for MERT begins with Monte Carlo generation of electron beamlets. Electrons are transported with proper in-air scattering and the dose is tallied in the phantom for each beamlet. An optimized beamlet plan may be calculated using inverse-planning methods. Step-and-shoot leaf sequences are generated for the intensity maps and dose distributions recalculated using Monte Carlo simulations. Here, scatter and leakage from the leaves are properly accounted for by transporting electrons through the eMLC geometry. The weights for the segments of the plan are re-optimized with the leaf positions fixed and bremsstrahlung leakage and electron scatter doses included. This optimization gives the final optimized plan. It is shown that a significant portion of the calculation time is spent transporting particles in the leaves. However, this is necessary since optimizing segment weights based on a model in which leaf transport is ignored results in an improperly optimized plan with overdosing of target and critical structures. A method of rapidly calculating the bremsstrahlung contribution is presented and shown to be an efficient solution to this problem. A homogeneous model target and a 2D breast plan are presented. The potential use of this tool in clinical planning is discussed.     

Comprehensive evaluation of a commercial macro Monte Carlo electron dose calculation implementation using a standard verification data set

A commercial electron dose calculation software implementation based on the macro Monte Carlo algorithm has recently been introduced. We have evaluated the performance of the system using a standard verification data set comprised of two-dimensional (2D) dose distributions in the transverse plane of a 15 X 15 cm2 field. The standard data set was comprised of measurements performed for combinations of 9-MeV and 20-MeV beam energies and five phantom geometries. The phantom geometries included bone and air heterogeneities, and irregular surface contours. The standard verification data included a subset of the data needed to commission the dose calculation. Additional required data were obtained from a dosimetrically equivalent machine. In addition, we performed 2D dose measurements in a water phantom for the standard field sizes, a 4 cm X 4 cm field, a 3 cm diameter circle, and a 5 cm X 13 cm triangle for the 6-, 9-, 12-, 15-, and 18-MeV energies of a Clinac 21EX. Output factors were also measured. Synthetic CT images and structure contours duplicating the measurement configurations were generated and transferred to the treatment planning system. Calculations for the standard verification data set were performed over the range of each of the algorithm parameters: statistical precision, grid-spacing, and smoothing. Dose difference and distance-to-agreement were computed for the calculation points. We found that the best results were obtained for the highest statistical precision, for the smallest grid spacing, and for smoothed dose distributions. Calculations for the 21EX data were performed using parameters that the evaluation of the standard verification data suggested would produce clinically acceptable results. The dose difference and distance-to-agreement were similar to that observed for the standard verification data set except for the portion of the triangle field narrower than 3 cm for the 6- and 9-MeV electron beams. The output agreed with measurements to within 2%, with the exception of the 3-cm diameter circle and the triangle for 6 MeV, which were within 5%. We conclude that clinically acceptable results may be obtained using a grid spacing that is no larger than approximately one-tenth of the distal falloff distance of the electron depth dose curve (depth from 80% to 20% of the maximum dose) and small relative to the size of heterogeneities. For judicious choices of parameters, dose calculations agree with measurements to better than 3% dose difference and 3-mm distance-to-agreement for fields with dimensions no less than about 3 cm.     

The influence of lateral electronic disequilibrium on the radiation treatment planning for lung cancer irradiation

Using higher energy photons can obtain better target dose uniformity and skin sparing for treating deep lesions, but the effect of lacking lateral scattering in the low-density lung may degrade the target coverage. To analyze the influence of lateral electronic disequilibrium on the radiation treatment planning for lung cancer, three dimension conformal treatment (3D-CRT) plans of using 6 MV and 18 MV X-ray respectively for a lung cancer case have been worked out by using pencil beam algorithm and collapsed cone algorithm provided by Helax-TMS treatment planning system for the same radiation field arrangement for both energies. Dose volume histogram (DVH) in target and organs at risk (OARs) are used for comparison of different plans. The study shows that using pencil beam algorithm, the target DVH are similar for 6 MV and 18 MV plan. However, using collapsed cone algorithm that can make account of lateral electron scattering, the target is underdosed. The change is even more pronounced for 18 MV plan. The doses for lung and spinal cord are similar for these two energies and two algorithms. Therefore, for lung cancer, dose calculation algorithm should have the ability of handling accurately the effect of the tissue density heterogeneity. It is better to use the lower-energy photons (6 MV) than to use the higher-energy photons (18 MV).     

Intraoperative cone-beam CT for guidance of head and neck surgery: Assessment of dose and image quality using a C-arm prototype

Cone-beam computed tomography (CBCT) with a flat-panel detector represents a promising modality for intraoperative imaging in interventional procedures, demonstrating sub-mm three-dimensional (3D) spatial resolution and soft-tissue visibility. Measurements of patient dose and in-room exposure for CBCT-guided head and neck surgery are reported, and the 3D imaging performance as a function of dose and other acquisition/reconstruction parameters is investigated. Measurements were performed on a mobile isocentric C-arm (Siemens PowerMobil) modified in collaboration with Siemens Medical Solutions (Erlangen, Germany) to provide flat-panel CBCT. Imaging dose was measured in a custom-built 16 cm cylindrical head phantom at four positions (isocenter, anterior, posterior, and lateral) as a function of kVp (80-120 kVp) and C-arm trajectory ("tube-under" and "tube-over" half-rotation orbits). At 100 kVp, for example ("tube-under" orbit), the imaging dose was 0.059 (isocenter), 0.022 (anterior), 0.10 (posterior), and 0.056 (lateral) mGy/ mAs, with scans at approximately 50 and approximately 170 mAs typical for visualization of bony and soft-tissue structures, respectively. Dose to radiosensitive structures (viz., the eyes and thyroid) were considered in particular: significant dose sparing to the eyes (a factor of 5) was achieved using a "tube-under" (rather than "tube-over") half-rotation orbit; a thyroid shield (0.5 mm Pb-equivalent) gave moderate reduction in thyroid dose due to x-ray scatter outside the primary field of view. In-room exposure was measured at positions around the operating table and up to 2 m from isocenter. A typical CBCT scan (10 mGy to isocenter) gave in-air exposure ranging from 29 mR (0.26 mSv) at 35 cm from isocenter, to <0.5 mR (<0.005 mSv) at 2 m from isocenter. Three-dimensional (3D) image quality was assessed in CBCT reconstructions of an anthropomorphic head phantom containing contrast-detail spheres (11-103 HU; 1.6-12.7 mm) and a natural human skeleton. The contrast-to-noise ratio (CNR) was evaluated across a broad range of dose (0.6-23.3 mGy). CNR increased as the square root of dose, with excellent visualization of bony and soft-tissue structures achieved at approximately 3 mGy (0.10 mSv) and approximately 10 mGy (0.35 mSv), respectively. The prototype C-arm demonstrates CBCT image quality sufficient for guidance of head and neck procedures based on soft-tissue and bony anatomy at dose levels low enough for repeat intraoperative imaging, with total dose over the course of the procedure comparable to or less than the effective dose of a typical (2 mSv) diagnostic CT of the head.     

Measurement of the dose components of fast and thermal neutrons and photons from a 0.1 mg 252Cf source in water for brachytherapy treatment planning

The use of 252Cf in brachytherapy is expected to be more effective with the therapy of bulky tumors than the conventional therapy with photons. For treatment planning a code developed for calculation of gamma dose was used to generate the dose distributions of fast and 10B enhanced thermal neutrons and photons. Dose distributions of these components measured with ionization chambers and a GM counter were fitted to analytical functions as required by the modified treatment planning program. A comparison of these experimental results and the treatment planning output indicate good agreement. Therefore, the program may be used to optimize the brachytherapy procedure considering all three dose components. A realistic case of a patient being treated with conventional brachytherapy has been used to calculate the dose distribution that would be obtained by use of the 252Cf source.     

Dosimetric verification of a commercial inverse treatment planning system

A commercial three-dimensional (3D) inverse treatment planning system, Corvus (Nomos Corporation, Sewickley, PA), was recently made available. This paper reports our preliminary results and experience with commissioning this system for clinical implementation. This system uses a simulated annealing inverse planning algorithm to calculate intensity-modulated fields. The intensity-modulated fields are divided into beam profiles that can be delivered by means of a sequence of leaf settings by a multileaf collimator (MLC). The treatments are delivered using a computer-controlled MLC. To test the dose calculation algorithm used by the Corvus software, the dose distributions for single rectangularly shaped fields were compared with water phantom scan data. The dose distributions predicted to be delivered by multiple fields were measured using an ion chamber that could be positioned in a rotatable cylindrical water phantom. Integrated charge collected by the ion chamber was used to check the absolute dose of single- and multifield intensity modulated treatments at various spatial points. The measured and predicted doses were found to agree to within 4% at all measurement points. Another set of measurements used a cubic polystyrene phantom with radiographic film to record the radiation dose distribution. The films were calibrated and scanned to yield two-dimensional isodose distributions. Finally, a beam imaging system (BIS) was used to measure the intensity-modulated x-ray beam patterns in the beam's-eye view. The BIS-measured images were then compared with a theoretical calculation based on the MLC leaf sequence files to verify that the treatment would be executed accurately and without machine faults. Excellent correlation (correlation coefficients > or = 0.96) was found for all cases. Treatment plans generated using intensity-modulated beams appear to be suitable for treatment of irregularly shaped tumours adjacent to critical structures. The results indicated that the system has potential for clinical radiation treatment planning and delivery and may in the future reduce treatment complexity.     

Target volume dose considerations in proton beam treatment planning for lung tumors

We performed a treatment planning study in order to gather basic insight in the effect of setup errors and breathing motion on the cumulative proton dose to a lung tumor. We used a simplified geometry that simulates a 50 mm diameter gross tumor volume (GTV) located centrally inside lung tissue. The GTV was expanded with a uniform 5 mm margin into a clinical target volume (CTV) and into a variety of planning target volume (PTV's). Proton beam apertures were designed to conform the prescribed dose laterally to the PTV while the range compensator was designed to provide distal coverage of the CTV. Different smearing distances were applied to the range compensators, and the cumulative dose in the CTV was evaluated for different combinations of breathing motion and systematic setup errors. Evaluation parameters were the dose to 99% of the CTV (D99) and the equivalent uniform dose (EUD), with a surviving fraction at 2 Gy of SF2 = 0.5. For a single proton field designed to a 15 mm expansion of the CTV and without smearing applied to the range compensator, D99 of the CTV reduced from 96% for no tumor displacement to 41% and 13% for systematic setup errors of 5 and 10 mm, respectively. For a representative clinical combination, of 5 mm systematic error and 10 mm breathing amplitude, the EUD of the CTV was about 40 Gy (prescribed dose 70 Gy) regardless the CTV to PTV margin, and without smearing. Smearing the range compensator increases the dose to the CTV substantially with a lateral margin and smearing distance of 7.5 mm providing ample tumor coverage. In this latter case, D99 of the target volume increased to 87% for a single field treatment plan. Smearing does, however, lead to an increase in dose to normal tissues distal to the clinical target volume. Next to countering geometric mismatches due to patient setup, smearing can also be used to counter the detrimental effects of breathing motion on the dose to the clinical target volume. We show that the lateral margin and smearing distance can be substantially smaller than the maximum tumor displacement due to setup errors and patient breathing, as measured by the D99 and the EUD.     

Monte Carlo based modulated electron beam treatment planning using a few-leaf electron collimator--feasibility study

Energy modulated electron beam therapy with conventional clinical accelerators has lagged behind photon IMRT despite its potential to achieve highly conformal dose distributions in superficial targets. One of the reasons for this is the absence of an automated collimating device that allows for the flexible delivery of a series of variable field openings. Electron-specific multileaf collimators attached to the bottom of the applicator require the use of a large number of motors and suffer from being relatively bulky and impractical for head and neck sites. In this work, we investigate the treatment planning aspects of a proposed 'few-leaf' electron collimator (FLEC) that consists of four motor-driven trimmer bars at the end of the applicator. The device is designed to serve as an accessory to standard equipment and allows for the shaping of any irregular field by combination of rectangular fieldlets. Using a Monte Carlo model of the FLEC, dose distributions are optimized using a simulated annealing (SA) inverse planning algorithm based on a limited number of Monte Carlo pre-generated, realistic phantom-specific dose kernels and user-specified dose-volume constraints. Using a phantom setup with an artificial target enclosed by organs at risk (OAR) as well as using a realistic patient case, we demonstrate that highly conformal distributions can be generated. Estimates of delivery times are made and show that a full treatment fraction can be kept to 15 min or less.     

Treatment planning for a small animal using Monte Carlo simulation

The development of a small animal model for radiotherapy research requires a complete setup of customized imaging equipment, irradiators, and planning software that matches the sizes of the subjects. The purpose of this study is to develop and demonstrate the use of a flexible in-house research environment for treatment planning on small animals. The software package, called DOSCTP, provides a user-friendly platform for DICOM computed tomography-based Monte Carlo dose calculation using the EGSnrcMP-based DOSXYZnrc code. Validation of the treatment planning was performed by comparing the dose distributions for simple photon beam geometries calculated through the Pinnacle3 treatment planning system and measurements. A treatment plan for a mouse based on a CT image set by a 360-deg photon arc is demonstrated. It is shown that it is possible to create 3D conformal treatment plans for small animals with consideration of inhomogeneities using small photon beam field sizes in the diameter range of 0.5-5 cm, with conformal dose covering the target volume while sparing the surrounding critical tissue. It is also found that Monte Carlo simulation is suitable to carry out treatment planning dose calculation for small animal anatomy with voxel size about one order of magnitude smaller than that of the human.