The importance of accurate linear accelerator head modelling for IMRT Monte Carlo calculations.

Two Monte Carlo dose engines for radiotherapy treatment planning, namely a beta release of Peregrine and MCDE (Monte Carlo dose engine), were compared with Helax-TMS (collapsed cone superposition convolution) for a head and neck patient for the Elekta SLi plus linear accelerator. Deviations between the beta release of Peregrine and MCDE up to 10% were obtained in the dose volume histogram of the optical chiasm. It was illustrated that the differences are not caused by the particle transport in the patient, but by the modelling of the Elekta SLi plus accelerator head and more specifically the multileaf collimator (MLC). In MCDE two MLC modules (MLCQ and MLCE) were introduced to study the influence of the tongue-and-groove geometry, leaf bank tilt and leakage on the actual dose volume histograms. Differences in integral dose in the optical chiasm up to 3% between the two modules have been obtained. For single small offset beams though the FWHM of lateral profiles obtained with MLCE can differ by more than 1.5 mm from profiles obtained with MLCQ. Therefore, and because the recent version of MLCE is as fast as MLCQ, we advise to use MLCE for modelling the Elekta MLC. Nevertheless there still remains a large difference (up to 10%) between Peregrine and MCDE. By studying small offset beams we have shown that the profiles obtained with Peregrine are shifted, too wide and too flat compared with MCDE and phantom measurements. The overestimated integral doses for small beam segments explain the deviations observed in the dose volume histograms. The Helax-TMS results are in better agreement with MCDE, although deviations exceeding 5% have been observed in the optical chiasm. Monte Carlo dose deviations of more than 10% as found with Peregrine are unacceptable as an influence on the clinical outcome is possible and as the purpose of Monte Carlo treatment planning is to obtain an accuracy of 2%. We would like to emphasize that only the Elekta MLC has been tested in this work, so it is certainly possible that alpha releases of Peregrine provide more accurate results for other accelerators.

MCDE: a new Monte Carlo dose engine for IMRT.

A new accurate Monte Carlo code for IMRT dose computations, MCDE (Monte Carlo dose engine), is introduced. MCDE is based on BEAMnrc/DOSXYZnrc and consequently the accurate EGSnrc electron transport. DOSXYZnrc is reprogrammed as a component module for BEAMnrc. In this way both codes are interconnected elegantly, while maintaining the BEAM structure and only minimal changes to BEAMnrc.mortran are necessary. The treatment head of the Elekta SLiplus linear accelerator is modelled in detail. CT grids consisting of up to 200 slices of 512 x 512 voxels can be introduced and up to 100 beams can be handled simultaneously. The beams and CT data are imported from the treatment planning system GRATIS via a DICOM interface. To enable the handling of up to 50 x 10(6) voxels the system was programmed in Fortran95 to enable dynamic memory management. All region-dependent arrays (dose, statistics, transport arrays) were redefined. A scoring grid was introduced and superimposed on the geometry grid, to be able to limit the number of scoring voxels. The whole system uses approximately 200 MB of RAM and runs on a PC cluster consisting of 38 1.0 GHz processors. A set of in-house made scripts handle the parallellization and the centralization of the Monte Carlo calculations on a server. As an illustration of MCDE, a clinical example is discussed and compared with collapsed cone convolution calculations. At present, the system is still rather slow and is intended to be a tool for reliable verification of IMRT treatment planning in the case of the presence of tissue inhomogeneities such as air cavities.

Leaf position optimization for step-and-shoot IMRT.

PURPOSE: To describe the theoretical basis, the algorithm, and implementation of a tool that optimizes segment shapes and weights for step-and-shoot intensity-modulated radiation therapy delivered by multileaf collimators. METHODS AND MATERIALS: The tool, called SOWAT (Segment Outline and Weight Adapting Tool) is applied to a set of segments, segment weights, and corresponding dose distribution, computed by an external dose computation engine. SOWAT evaluates the effects of changing the position of each collimating leaf of each segment on an objective function, as follows. Changing a leaf position causes a change in the segment-specific dose matrix, which is calculated by a fast dose computation algorithm. A weighted sum of all segment-specific dose matrices provides the dose distribution and allows computation of the value of the objective function. Only leaf position changes that comply with the multileaf collimator constraints are evaluated. Leaf position changes that tend to decrease the value of the objective function are retained. After several possible positions have been evaluated for all collimating leaves of all segments, an external dose engine recomputes the dose distribution, based on the adapted leaf positions and weights. The plan is evaluated. If the plan is accepted, a segment sequencer is used to make the prescription files for the treatment machine. Otherwise, the user can restart SOWAT using the new set of segments, segment weights, and corresponding dose distribution. The implementation was illustrated using two example cases. The first example is a T1N0M0 supraglottic cancer case that was distributed as a multicenter planning exercise by investigators from Rotterdam, The Netherlands. The exercise involved a two-phase plan. Phase 1 involved the delivery of 46 Gy to a concave-shaped planning target volume (PTV) consisting of the primary tumor volume and the elective lymph nodal regions II-IV on both sides of the neck. Phase 2 involved a boost of 24 Gy to the primary tumor region only. SOWAT was applied to the Phase 1 plan. Parotid sparing was a planning goal. The second implementation example is an ethmoid sinus cancer case, planned with the intent of bilateral visus sparing. The median PTV prescription dose was 70 Gy with a maximum dose constraint to the optic pathway structures of 60 Gy. RESULTS: The initial set of segments, segment weights, and corresponding dose distribution were obtained, respectively, by an anatomy-based segmentation tool, a segment weight optimization tool, and a differential scatter-air ratio dose computation algorithm as external dose engine. For the supraglottic case, this resulted in a plan that proved to be comparable to the plans obtained at the other institutes by forward or inverse planning techniques. After using SOWAT, the minimum PTV dose and PTV dose homogeneity increased; the maximum dose to the spinal cord decreased from 38 Gy to 32 Gy. The left parotid mean dose decreased from 22 Gy to 19 Gy and the right parotid mean dose from 20 to 18 Gy. For the ethmoid sinus case, the target homogeneity increased by leaf position optimization, together with a better sparing of the optical tracts. CONCLUSIONS: By using SOWAT, the plans improved with respect to all plan evaluation end points. Compliance with the multileaf collimator constraints is guaranteed. The treatment delivery time remains almost unchanged, because no additional segments are created.

An implementation strategy for IMRT of ethmoid sinus cancer with bilateral sparing of the optic pathways.

PURPOSE: To develop a protocol for the irradiation of ethmoid sinus cancer, with the aim of sparing binocular vision; of developing a strategy of intensity-modulated radiation therapy (IMRT) planning that produces dose distributions that (1) are consistent with the protocol prescriptions and (2) are deliverable by static segmental IMRT techniques within a 15-minute time slot; of fine tuning the implementation strategy to a class solution approach that is sufficiently automated and efficient, allowing routine clinical application; of reporting on the early clinical implementation involving 11 patients between February 1999 and July 2000. patients and methods: Eleven consecutive T1-4N0M0 ethmoid sinus cancer patients were enrolled in the study. For Patients 1-8, a first protocol was implemented, defining a planning target volume prescription dose of 60 to 66 Gy in 30-33 fractions and a maximum dose (Dmax) of 50 Gy to optic pathway structures and spinal cord and limit of 60 Gy to brainstem. For Patients 9-11, an adapted (now considered mature) protocol was implemented, defining a (planning target volume) prescription dose of 70 Gy in 35 fractions and a Dmax to optic pathway structures and brainstem of 60 Gy and to spinal cord of 50 Gy. RESULTS: The class solution-directed strategy developed during this study reduced the protocol translation process from a few days to about 2 hours of planner time. The mature class solution involved the use of 7 beam incidences (20-37 segments), which could be delivered within a 15-minute time slot. Acute side effects were limited and mild. None of the patients developed dry eye syndrome or other visual disturbances. The follow-up period is too short for detection of retinopathy or optic nerve and chiasm toxicity. CONCLUSION: Conventional radiotherapy of ethmoid sinus tumors is associated with serious morbidity, including blindness. We hypothesize that IMRT has the potential to save binocular vision. The dose to the optic pathway structures can be reduced selectively by IMRT. Further enrollment of patients and longer follow-up will show whether the level of reduction tested by the clinical protocol is sufficient to save binocular vision. An adaptive strategy of IMRT planning was too inefficient for routine clinical practice. A class solution-directed strategy improved efficiency by eliminating human trial and error during the IMRT planning process.

An isocenter position verification device for electronic portal imaging: physical and dosimetrical characteristics.

The physical and dosimetrical characteristics of a device, designed to visualize the isocenter position on electronic portal images, were examined. The device, to be mounted on the gantry head of the accelerator, containing five spheric lead markers, was designed in order to visualize the isocenter position on portal images. A quality control device was designed to check the reliability of this technique. The disturbance of the dose distribution by the markers was studied with gel dosimetry. The use of markers resulted in a precise and accurate method to visualize the isocenter on portal images. A maximum underdosage of 11%, due to attenuation by the markers, was observed. The use of markers to visualize the isocenter position on portal images, is a fast and reliable method when analyzing patient setup errors with online electronic portal imaging.

Validation of MR-based polymer gel dosimetry as a preclinical three-dimensional verification tool in conformal radiotherapy.

The aim of this work was to investigate MR-based polymer gel dosimetry as a three-dimensional (3D) dosimetry technique in conformal radiotherapy. A cylindrical container filled with polymer gel was placed in a water-filled torso phantom to verify a treatment plan for the conformal irradiation of a mediastinal tumor located near the esophagus. Magnetic resonance spin-spin relaxation rate images were acquired and, after calibration, converted to absorbed dose distributions. The dose maps were compared with dose distributions measured using radiographic film. The average root-mean-square structural deviation, for the complete dose distribution, amounted to less than 3% between gel and film dose maps. It may be expected that MR gel dosimetry will become a valuable tool in the verification of 3D dose distributions. The influence of imaging artifacts arising from eddy currents, temperature drift during scanning, and B1 field inhomogeneity on the dose maps was taken into account and minimized.

Clinical delivery of intensity modulated conformal radiotherapy for relapsed or second-primary head and neck cancer using a multileaf collimator with dynamic control.

BACKGROUND AND PURPOSE: Concave dose distributions generated by intensity modulated radiotherapy (IMRT) were applied to re-irradiate three patients with pharyngeal cancer. PATIENTS, MATERIALS AND METHODS: Conventional radiotherapy for oropharyngeal (patients 1 and 3) or nasopharyngeal (patient 2) cancers was followed by relapsing or new tumors in the nasopharynx (patients 1 and 2) and hypopharynx (patient 3). Six non-opposed coplanar intensity modulated beams were generated by combining non-modulated beamparts with intensities (weights) obtained by minimizing a biophysical objective function. Beamparts were delivered by a dynamic MLC (Elekta Oncology Systems, Crawley, UK) forced in step and shoot mode. RESULTS AND CONCLUSIONS: Median PTV-doses (and ranges) for the three patients were 73 (65-78), 67 (59-72) and 63 (48-68) Gy. Maximum point doses to brain stem and spinal cord were, respectively, 67 Gy (60% of volume below 30 Gy) and 32 Gy (97% below 10 Gy) for patient 1; 60 Gy (69% below 30 Gy) and 34 Gy (92% below 10 Gy) for patient 2 and 21 Gy (96% below 10 Gy) at spinal cord for patient 3. Maximum point doses to the mandible were 69 Gy for patient 1 and 64 Gy for patient 2 with, respectively, 66 and 92% of the volume below 20 Gy. A treatment session, using the dynamic MLC, was finished within a 15-min time slot.

Postoperative radiotherapy of paranasal sinus tumours: a challenge for intensity modulated radiotherapy.

BACKGROUND AND PURPOSE: Intensity modulated radiotherapy (IMRT) is used in our department for treatment of paranasal sinuses. We describe the methodology that was developed together with the clinical implementation, illustrated by a case report. MATERIAL AND METHODS: Patient history, treatment and short follow-up are described. An IMRT, obtained by superposition of static beam segments was implemented. Electronic portal images, compared to digitally reconstructed radiographs (DRR) were used to evaluate and adjust patient positioning. RESULTS, DISCUSSION AND CONCLUSION: IMRT is an appropriate and feasible treatment technique for head and neck cancer in anatomical regions that are difficult to treat. A high tumour dose can be combined with a good sparing of the surrounding organs at risk (OAR's).

Conformal radiotherapy of Stage III non-small cell lung cancer: a class solution involving non-coplanar intensity-modulated beams.

PURPOSE: We developed a semiautomatic class solution to irradiate centrally located Stage III non-small cell lung cancer (NSCLC), involving a beam intensity modulation technique and optimization using a biophysical cost function. METHODS AND MATERIALS: Treatment for 10 patients with Stage III NSCLC was planned, using a conventional three- or four-beam three-dimensional (3D) technique and two techniques involving, respectively, seven (BIM1) and five (BIM2) noncoplanar beam incidences with intensity modulation. Two planning target volumes were defined: PTV1 included macroscopic tumor volume and PTV2 included macroscopic and microscopic disease. Beams were divided into beam parts (segments) and their outlines were defined during virtual simulation. Optimization using a biophysical cost function determined beam weights, segment weights, and wedge angles. Biological end points included tumor control probability of both target volumes (TCP1 and TCP2) and normal tissue complication probability (NTCP) of heart, lung, and spinal cord. The resulting uncomplicated local control probability (UCLP) was calculated. Physical end points included dose at PTV1 expressed as a dose minimum and dose maximum. Target-dose inhomogeneity was constrained in all plans. RESULTS: Concerning tumor evaluation, TCP1 was 74% (range 54-89%) for the 3D plan, 78.0% (range 62-94%) for BIM1, and 86.0% (range 59-93%) for BIM2. TCP1*TCP2 was, respectively, 67.0% (range 39-81%), 73.0% (range 56-94%), and 81.0% (range 54-93%). Minimum doses to PTV1 were 85, 80, and 88 Gy with the three respective techniques, while dose maxima were 89, 101, and 100 Gy. NTCPs of lung were 45.0% (range 11-75%) for 3D, 19.5% (range 8-59%) for BIM1, and 24.5% (range 3-61%) for BIM2. NTCPs of heart and spinal cord were comparable for all techniques. ULCPs were 37.0% (range 9-73%), 52.5% (range 22-86%), and 60.0% (range 20-85%), respectively. Applying physical limits to ensure clinical safety, minimum doses at PTV1 were recalculated. These were 72, 71, and 74 Gy for 3D, BIM1, and BIM2, respectively. CONCLUSION: The BIM2 plan is a candidate class solution for dose escalation studies in centrally located Stage III NSCLC.

3D conformal intensity-modulated radiotherapy planning: interactive optimization by constrained matrix inversion.

BACKGROUND AND PURPOSE: This paper presents a method for interactive optimization of 3D conformal intensity-modulated radiotherapy plans employing a quadratic objective that also contains dose limitations in the organs at risk. This objective function is minimized by constrained matrix inversion (CMI) that follows the same approach as the gradient technique using matrix notation. MATERIALS AND METHODS: Sherouse's GRATIS radiotherapy design system is used to determine the outlines of the target volume and the organs at risk and to input beam segments which are given by the beam segmentation technique. This technique defines the beam incidences and the beam segmentation. The weights of the segments are then calculated using a quadratic objective function and CMI. The objective function to be minimized consists of two components based on the planning target volume (PTV) and the organ at risk (OAR) with an importance factor w associated with the OAR. RESULTS: Optimization is tested for concave targets in the head and neck region wrapping around the spinal cord. For a predefined w-value, segment weights are optimized within a few seconds on a DEC Alpha 3000. In practice, 5-10 w-values have to be tested, making optimization a less than 5 min procedure. This optimization procedure predicts the possibility of target dose escalation for a tumour in the lower neck to 120-150 Gy without exceeding the spinal cord tolerance, whereas human planners could not increase the dose above 65-80 Gy. CONCLUSIONS: Treatment plans optimized using a quadratic objective function and the CMI algorithm are superior to those which are generated by human planners. The optimization algorithm is very fast and allows interactive use. Quadratic optimization by CMI is routinely used by clinicians at the Division of Radiotherapy, U.Z.-Gent.

Three-dimensional dosimetry using polymer gel and magnetic resonance imaging applied to the verification of conformal radiation therapy in head-and-neck cancer.

BACKGROUND AND PURPOSE: It was our aim to investigate NMR-based BANG gel dosimetry as a three-dimensional dosimetry technique in conformal radiotherapy. MATERIALS AND METHODS: The BANG gel consisting of gelatin, water and co-monomers was first validated in a cylindrical glass flask for a single standard beam. Next, the gel contained in a human neck-shaped cast was used to verify a treatment plan for the conformal irradiation of a concave tumour in the lower neck. Magnetic resonance relaxation rate images were acquired and, based on an appropriate calibration of the gel, converted to absorbed dose distributions. The resulting maps were compared with dose distributions measured using radiographic film. RESULTS: The gel-measured dose profiles of standard beams agreed within 3% (root mean square difference) with the profiles measured with high spatial resolution by a diamond detector. For the multi-beam conformal treatment, the difference map between gel-measured and film-measured dose distributions revealed a noise component and a more systematic deviation including structural or space-coherent patterns. The mean absolute value of the difference amounted to 8%. A number of possible causes for this deviation are designated. CONCLUSIONS: Polymer gel dosimetry in combination with magnetic resonance imaging is a promising method for dosimetric verification of conformal radiotherapy.

Non-coplanar beam intensity modulation allows large dose escalation in stage III lung cancer.

PURPOSE: To evaluate the feasibility of dose escalation in stage III non-small cell lung cancer, we compared standard coplanar (2D) with non-coplanar beam arrangements, without (3D) and with beam intensity modulation (3D-BIM). MATERIALS AND METHODS: This study was a planning effort performed on a non-selected group of 10 patients. Starting from a serial CT scan, treatment planning was performed using Sherouse's GRATIS 3D planning system. Two target volumes were defined; gross tumor volume (GTV) defined a high-dose target volume that had to receive a dose of at least 80 Gy and GTV plus the lymph node regions with >10% probability of invasion defined an intermediate-dose target volume (GTV + N). It was our intention to irradiate GTV + N up to 56 Gy or more. If the prescribed doses on GTV and GTV + N could not be reached with either the 2D or 3D technique, a 3D-BIM plan was performed. The 3D-BIM plan was a class solution involving identical gantry angles, segment arrangements and relative segment weights for all patients. Dose volume histograms for GTV, GTV + N, lung and spinal cord were calculated. Criteria for tolerance were met if no points inside the spinal cord exceeded 50 Gy and if at least 50% of the lung volume received less than 20 Gy. Under these constraints, maximal achievable doses to GTV and GTV + N were calculated. RESULTS: In all 2D plans, spinal cord was the limiting factor and the prescribed doses for GTV and GTV + N could not be reached in any patient. The non-coplanar 3D plan resulted in a satisfying solution in 4 out of 10 patients under the same constraints. In comparison with 2D, the minimum dose in GTV + N was increased. Six patients had to be planned with the 3D-BIM technique. The theoretical minimum dose to GTV + N ranged between 56 and 98 Gy. The delivery of 80 Gy or more to GTV was possible in all patients. For a minimal dose of 80 Gy to GTV, the maximal dose to any point of the spinal cord varied between 27 and 46 Gy. The lung volume receiving more than 20 Gy ranged from 26 to 46%. CONCLUSION: The potential of 3D-BIM for dose escalation is explained as follows: (i) compared to other planning techniques, a larger amount of lung tissue can be spared by using beam directions that are well-aligned with the mediastinal structures. Such beam directions have narrow angles with the sagittal plane; (ii) dividing all beams into segments with well-specified geometrical restrictions in relation to the spinal cord and well-defined relative weights results in a lower dose to the spinal cord.