Possibilities for tailoring dose distributions through the manipulation of electron beam characteristics

Delivery of intensity-modulated electron therapy by mechanical scanning: An algorithm study

Purpose

In principle, intensity-modulated electron therapy (IMET) can be delivered through mechanical scanning, with a robotic arm mounting a linac.

Materials and methods

Here is a scanning algorithm to identify the back-and-forth, top-to-bottom (zigzag) pattern scan sequence. The algorithm includes generating beam positions with a uniform resolution according to the applicator size; adopting discrete energies to achieve the depth of 90% dose by compositing energies; selecting energy by locating the target's distal edge; and employing the energy-by-energy scan strategy for step-and-shoot discrete scanning. After a zigzag scan sequence is obtained, the delivery order of the scan spots is optimized by fast simulated annealing (FSA) to minimize the path length. For algorithm evaluation, scan sequences were generated using the computed tomography data of 10 patients with pancreatic cancer undergoing intraoperative radiotherapy, and the results were compared between the zigzag path and an optimized path. A simple calculation of the treatment delivery time, which comprises the irradiation time, the total robotic arm moving time, the time for energy switch, and the time to stop and restart the beam, was also made.

Results

In these clinical cases, FSA optimization shortened the path lengths by 12%-43%. Assuming the prescribed dose was 15 Gy, machine dose rate was 15 Gy/s, energy switch time was 2 s, stop and restart beam time was 20 ms, and robotic arm move speed was 50 mm/s, the average delivery time was 124±38 s. The largest reduction in path length yielded an approximately 10% reduction in the delivery time, which can be further reduced by increasing the machine dose rate and the robotic arm speed, decreasing the time for energy switch, and/or developing more efficient algorithms.

Conclusion

Mechanically scanning IMET is potentially feasible and worthy of further exploration.

Introduction to passive electron intensity modulation

This work introduces a new technology for electron intensity modulation, which uses small area island blocks within the collimating aperture and small area island apertures in the collimating insert. Due to multiple Coulomb scattering, electrons contribute dose under island blocks and lateral to island apertures. By selecting appropriate lateral positions and diameters of a set of island blocks and island apertures, for example, a hexagonal grid with variable diameter circular island blocks, intensity modulated beams can be produced for appropriate air gaps between the intensity modulator (position of collimating insert) and the patient. Such a passive radiotherapy intensity modulator for electrons (PRIME) is analogous to using physical attenuators (metal compensators) for intensity modulated x‐ray therapy (IMXT). For hexagonal spacing, the relationship between block (aperture) separation (r) and diameter (d) and the local intensity reduction factor (IRF) is discussed. The PRIME principle is illustrated using pencil beam calculations for select beam geometries in water with half beams modulated by 70%–95% and for one head and neck field of a patient treated with bolus electron conformal therapy. Proof of principle is further illustrated by showing agreement between measurement and calculation for a prototype PRIME. Potential utilization of PRIME for bolus electron conformal therapy, segmented‐field electron conformal therapy, modulated electron radiation therapy, and variable surface geometries is discussed. Further research and development of technology for the various applications is discussed. In summary, this paper introduces a practical, new technology for electron intensity modulation in the clinic, demonstrates proof of principle, discusses potential clinical applications, and suggests areas of further research and development.

Dosimetric characteristics of an electron multileaf collimator for modulated electron radiation therapy

Modulated electron radiation therapy (MERT) has been proven as an effective way to deliver conformal dose distributions to shallow tumors while sparing distal critical structures and surrounding normal tissues. It had been shown that a dedicated electron multileaf collimator (eMLC) is necessary to reach the full potential of MERT. In this study, a manually-driven eMLC for MERT was investigated. Percentage depth dose (PDD) curves and profiles at different depths in a water tank were measured using ionization chamber and were also simulated using the Monte Carlo method. Comparisons have been performed between PDD curves and profiles collimated using the eMLC and conventional electron applicators with similar size of opening. Monte Carlo simulations were performed for all electron energies available (6, 9, 12, 15, 18 and 20 MeV) on a Varian 21EX accelerator. Monte Carlo simulation results were compared with measurements which showed good agreement (< 2%/1mm). The simulated dose distributions resulting from multiple static electron fields collimated by the eMLC agreed well with measurements. Further studies were carried out to investigate the properties of abutting electron beams using the eMLC, as it is an essential issue that needs to be addressed for optimizing the MERT outcome. A series of empirical formulas for abutting beams of different energies have been developed for obtaining the optimum gap sizes, which can highly improve the target dose uniformity.

Characterisation of mega-voltage electron pencil beam dose distributions: viability of a measurement-based approach

The concept of electron pencil-beam dose distributions is central to pencil-beam algorithms used in electron beam radiotherapy treatment planning. The Hogstrom algorithm, which is a common algorithm for electron treatment planning, models large electron field dose distributions by the superposition of a series of pencil beam dose distributions. This means that the accurate characterisation of an electron pencil beam is essential for the accuracy of the dose algorithm. The aim of this study was to evaluate a measurement based approach for obtaining electron pencil-beam dose distributions. The primary incentive for the study was the accurate calculation of dose distributions for narrow fields as traditional electron algorithms are generally inaccurate for such geometries. Kodak X-Omat radiographic film was used in a solid water phantom to measure the dose distribution of circular 12 MeV beams from a Varian 21EX linear accelerator. Measurements were made for beams of diameter, 1.5, 2, 4, 8, 16 and 32 mm. A blocked-field technique was used to subtract photon contamination in the beam. The "error function" derived from Fermi-Eyges Multiple Coulomb Scattering (MCS) theory for corresponding square fields was used to fit resulting dose distributions so that extrapolation down to a pencil beam distribution could be made. The Monte Carlo codes, BEAM and EGSnrc were used to simulate the experimental arrangement. The 8 mm beam dose distribution was also measured with TLD-100 microcubes. Agreement between film, TLD and Monte Carlo simulation results were found to be consistent with the spatial resolution used. The study has shown that it is possible to extrapolate narrow electron beam dose distributions down to a pencil beam dose distribution using the error function. However, due to experimental uncertainties and measurement difficulties, Monte Carlo is recommended as the method of choice for characterising electron pencil-beam dose distributions.

Characterization of an add-on multileaf collimator for electron beam therapy

An add-on multileaf collimator for electrons (eMLC) has been developed that provides computer-controlled beam collimation and isocentric dose delivery. The design parameters result from the design study by Gauer et al (2006 Phys. Med. Biol. 51 5987-6003) and were configured such that a compact and light-weight eMLC with motorized leaves can be industrially manufactured and stably mounted on a conventional linear accelerator. In the present study, the efficiency of an initial computer-controlled prototype was examined according to the design goals and the performance of energy- and intensity-modulated treatment techniques. This study concentrates on the attachment and gantry stability as well as the dosimetric characteristics of central-axis and off-axis dose, field size dependence, collimator scatter, field abutment, radiation leakage and the setting of the accelerator jaws. To provide isocentric irradiation, the eMLC can be placed either 16 or 28 cm above the isocentre through interchangeable holders. The mechanical implementation of this feature results in a maximum field displacement of less than 0.6 mm at 90 degrees and 270 degrees gantry angles. Compared to a 10 x 10 cm applicator at 6-14 MeV, the beam penumbra of the eMLC at a 16 cm collimator-to-isocentre distance is 0.8-0.4 cm greater and the depth-dose curves show a larger build-up effect. Due to the loss in energy dependence of the therapeutic range and the much lower dose output at small beam sizes, a minimum beam size of 3 x 3 cm is necessary to avoid suboptimal dose delivery. Dose output and beam symmetry are not affected by collimator scatter when the central axis is blocked. As a consequence of the broader beam penumbra, uniform dose distributions were measured in the junction region of adjacent beams at perpendicular and oblique beam incidence. However, adjacent beams with a high difference in a beam energy of 6 to 14 MeV generate cold and hot spots of approximately 15% in the abutting region. In order to improve uniformity, the energy of adjacent beams must be limited to 6 to 10 MeV and 10 to 14 MeV respectively. At the maximum available beam energy of 14 MeV, radiation leakage results mainly from the intraleaf leakage of approximately 2.5% relative dose which could be effectively eliminated at off-axis distances remote from the field edge by adjusting the jaw field size to the respective opening of the eMLC. Additionally, the interleaf and leaf-end leakage could be reduced by using a tongue-and-groove leaf shape and adjoining the leaf-ends off-axis respectively.

Improvement of field matching in segmented-field electron conformal therapy using a variable-SCD applicator

The purpose of the present study is to demonstrate that the use of an electron applicator with energy-dependent source-to-collimator distances (SCDs) will significantly improve the dose homogeneity for abutted electron fields in segmented-field electron conformal therapy (ECT). Multiple Coulomb scattering theory was used to calculate and study the P(80-20) penumbra width of off-axis dose profiles as a function of air gap and depth. Collimating insert locations with air gaps (collimator-to-isocenter distance) of 5.0, 7.5, 11.5, 17.5 and 19.5 cm were selected to provide equal P(80-20) at a depth of 1.5 cm in water for energies of 6, 9, 12, 16 and 20 MeV, respectively, for a Varian 2100EX radiation therapy accelerator. A 15 x 15 cm(2) applicator was modified accordingly, and collimating inserts used in the variable-SCD applicator for segmented-field ECT were constructed with diverging edges using a computer-controlled hot-wire cutter, which resulted in 0.27 mm accuracy in the abutted edges. The resulting electron beams were commissioned for the pencil-beam algorithm (PBA) on the Pinnacle(3) treatment planning system. Four hypothetical planning target volumes (PTVs) and one patient were planned for segmented-field ECT using the new variable-SCD applicator, and the resulting dose distributions were compared with those calculated for the identical plans using the conventional 95 cm SCD applicator. Also, a method for quality assurance of segmented-field ECT dose plans using the variable-SCD applicator was evaluated by irradiating a polystyrene phantom using the treatment plans for the hypothetical PTVs. Treatment plans for all four of the hypothetical PTVs using the variable-SCD applicator showed significantly improved dose homogeneity in the abutment regions of the segmented-field ECT plans. This resulted in the dose spread (maximum dose-minimum dose), sigma, and D(90-10) in the PTV being reduced by an average of 32%, 29% and 32%, respectively. Reductions were most significant for abutted fields of nonadjacent energies. Planning segmented-field ECT using the variable-SCD applicator for a patient with recurrent squamous cell carcinoma of the left ear showed the dose spread, sigma, and D(90-10) of the dose distribution in the PTV being reduced by an average of 38%, 22% and 22%, respectively. The measured and calculated dose in a polystyrene phantom resulting from the variable-SCD, segmented-field ECT plans for the hypothetical PTVs showed good agreement; however, isolated differences between dose calculation and measurement indicated the need for a more accurate dose algorithm than the PBA for segmented-field ECT. These results confirmed our hypothesis that using the variable-SCD applicator for segmented-field ECT results in the PTV dose distribution becoming more homogenous and being within the range of 85-105% of the 'given dose'. Clinical implementation of this method requires variable-SCD applicators, and the design used in the present work should be acceptable, as should our methods for construction of the inserts. Dose verification measurements in a polystyrene phantom and the recommended improvements in dose calculation should be appropriate for quality assurance of segmented-field ECT.

Design of a computer-controlled multileaf collimator for advanced electron radiotherapy

A multileaf collimator for electrons (eMLC) has been designed that fulfils the technical requirements for providing advanced irradiation techniques with electrons. In the present work, the basic design parameters of leaf material, leaf height, leaf width and number of leaves as well as leaf overtravel and leaf shape were determined such that an eMLC with motorized leaves can be manufactured by a company specialized in MLC technology. For this purpose, a manually driven eMLC with variable source-to-collimator distance (SCD) was used to evaluate the chosen leaf specification and investigate the impact of the SCD on the off-axis dose distribution. In order to select the final SCD of the eMLC, a compromise had to be found between maximum field size, minimum beam penumbra and necessary distance between eMLC and isocentre to eliminate patient realignments during gantry rotation. As a result, the eMLC is placed according to the target position at 72 and 84 cm SCD, respectively. This feature will be achieved by interchangeable distance holders. At these SCDs, the corresponding maximum field sizes at 100 cm source-to-isocentre distance are 20 x 20 cm and 17 x 17 cm, respectively. Finally, the off-axis dose distribution at the maximum opening of the eMLC was improved by fine-tuning the settings of the accelerator jaws and introducing trimmer bars above the eMLC. Following this optimization, a prototype eMLC consisting of 2 x 24 computer-controlled brass leaves is manufactured by 3D Line Medical Systems.

Intensity modulated radiation therapy with electrons using algorithm based energy/range selection methods

BACKGROUND AND PURPOSE

In recent years photon intensity modulated radiation therapy (IMRT) has gained attention due to its ability to improve conformity of dose distributions. A potential advantage of electron-IMRT is that the dose fall off in the depth dose curve makes it possible to modulate the dose distribution in the direction of the beam by selecting different electron energies. This paper examines the use of a computer based energy selection in combination with the IMRT technique to optimise the electron dose distribution.

MATERIALS AND METHODS

One centimetre square electron beamlets ranging from 2.5 to 50 MeV were pre-calculated in water using Monte Carlo methods. A modified IMRT optimisation tool was then used to find an optimum mix of electron energies and intensities. The main principles used are illustrated in some simple geometries and tested on two clinical cases of post-operated ca. mam.

RESULTS

It is clearly illustrated that the energy optimisation procedure lowers the dose to lung and heart and makes the dose in the target more homogeneous. Increasing the energy at steep gradients compensates for lack of target coverage at beam edges and steep gradients. Comparison with a clinically acceptable four segment plan indicates the advantage of the used electron IMRT technique.

CONCLUSIONS

Using an intensity optimised mix of computer selected electron energies has the potential to improve electron treatments for mastectomy patients with good target coverage and reduced dose to normal tissue such as lung and heart.

A comparative dosimetric study on tangential photon beams, intensity-modulated radiation therapy (IMRT) and modulated electron radiotherapy (MERT) for breast cancer treatment

Recently, energy- and intensity-modulated electron radiotherapy (MERT) has garnered a growing interest for the treatment of superficial targets. In this work. we carried out a comparative dosimetry study to evaluate MERT, photon beam intensity-modulated radiation therapy (IMRT) and conventional tangential photon beams for the treatment of breast cancer. A Monte Carlo based treatment planning system has been investigated, which consists of a set of software tools to perform accurate dose calculation, treatment optimization, leaf sequencing and plan analysis. We have compared breast treatment plans generated using this home-grown treatment optimization and dose calculation software forthese treatment techniques. The MERT plans were planned with up to two gantry angles and four nominal energies (6, 9, 12 and 16 MeV). The tangential photon treatment plans were planned with 6 MV wedged photon beams. The IMRT plans were planned using both multiple-gantry 6 MV photon beams or two 6 MV tangential beams. Our results show that tangential IMRT can reduce the dose to the lung, heart and contralateral breast compared to conventional tangential wedged beams (up to 50% reduction in high dose volume or 5 Gy in the maximum dose). MERT can reduce the maximum dose to the lung by up to 20 Gy and to the heart by up to 35 Gy compared to conventional tangential wedged beams. Multiple beam angle IMRT can significantly reduce the maximum dose to the lung and heart (up to 20 Gy) but it induces low and medium doses to a large volume of normal tissues including lung, heart and contralateral breast. It is concluded that MERT has superior capabilities to achieve dose conformity both laterally and in the depth direction, which will be well suited for treating superficial targets such as breast cancer.

A Monte Carlo investigation of fluence profiles collimated by an electron specific MLC during beam delivery for modulated electron radiation therapy

Modulated electron radiation therapy (MERT) is able to deliver conformal dose to shallow tumors while significantly reducing dose to distal structures and surrounding tissues. An electron specific multileaf collimator (eMLC) has been proposed and constructed as an effective means of delivering electron beams for MERT. The aim of this work is to apply the Monte Carlo method to investigate the fluence profiles collimated by the eMLC in order to achieve accurate beam delivery for MERT. In this work, the EGS4/BEAM code was used to simulate the eMLC collimated electron beams of 6-20 MeV generated from a Varian Clinac 2100C linear accelerator. An attempt was made to describe the fluence profiles with an analytic Sigmoid function. The function parameters were determined by the fittings of the Monte Carlo simulated fluence profiles. How the function parameters depend on the eMLC aperture size, the off-axis location, and the electron beam energy has been investigated. It has been found that the eMLC collimated fluence profiles are dependent on beam energy, while almost independent of leaf location or dimension of MLC aperture. There is little difference in the fluence profiles collimated by the leaf side and the leaf end for the straight-edged leaves. It is possible that these energy-dependent Sigmoid functions can serve as operators to account for the energy dependence of the eMLC collimated fluence profiles. These operators can be incorporated into the inverse planning algorithm to derive desired dose distributions using a set of electron beams of variable energy and field size suitable for delivery by the eMLC.

[IMCO(++)--a Monte Carlo based IMRT system]

The application of intensity modulated radiotherapy (IMRT) to dose escalation in the target volume sets particular demands in terms of accuracy of dose calculation. Dose calculation errors due to approximations are compensated by the optimization algorithm, a procedure that ultimately leads to incorrect fluence modulation. Such inaccuracies affect particularly the dose distribution in areas with secondary electron disequilibrium. In case tissues heterogeneity predominates, conventional dose calculation methods (such as Pencil Beam) can produce relative errors up to more than 10%. The accuracy can be significantly improved by the application of a Monte-Carlo (MC) algorithm. This paper describes a MC-based inverse treatment planning algorithm (IMCO++), based on a non-iterative approach with a feedback-controlling process. The convergence behavior of IMCO++ was investigated and the used MC dose-calculation codes MMms and XVMC were compared by means of a heterogeneous phantom. IMCO++ plans were optimized in various phantoms. All plans showed conformity in terms of dose distribution of the target volume and dose reduction in risk organs (according to the requirements of the target parameter), as well as a very fast convergence of the algorithm (in less than 10 optimization steps).

Electron conformal radiotherapy using bolus and intensity modulation

PURPOSE

Conformal electron beam therapy can be delivered using shaped bolus, which varies the penetration of the electrons across the incident beam so that the 90% isodose surface conforms to the distal surface of the planning target volume (PTV). Previous use of this modality has shown that the irregular proximal surface of the bolus causes the dose heterogeneity in the PTV to increase from 10%, the typical dose spread of a flat-water surface to approximately 20%. The present work evaluates the ability to restore dose homogeneity by varying the incident electron intensity.

METHODS AND MATERIALS

Three patients, one each with chest wall, thorax, and head-and-neck cancer, were planned using electron conformal therapy with bolus, with and without intensity modulation. Resulting dose distributions and dose-volume histograms were compared with non-intensity-modulated bolus plans.

RESULTS

In all cases, the DeltaD(90%-10%) for the PTV was reduced; for example, for the head-and-neck case, the DeltaD(90%-10%) for the PTV was reduced from 14.9% to 9.2%. This reduction in dose spread is a direct result of intensity modulation.

CONCLUSIONS

The results showed that intensity-modulated electron beams could significantly improve the dose homogeneity in the PTV for patients treated with electron conformal therapy using shaped bolus.

Optimized treatment planning for prostate cancer comparing IMPT, VHEET and 15 MV IMXT

The merits of intensity-modulated very-high energy electron therapy (VHEET) and intensity-modulated proton therapy (IMPT) in relation to intensity-modulated x-ray therapy (IMXT) with respect to the treatment of the prostate have been quantified. Optimized dose distributions were designed for 5-11 beams of 250 MeV VHEET and 15 MV IMXT as well as 1-9 beam ports of IMPT. In the case of the comparison between 250 MeV VHEET and 15 MV IMXT, it was found that the quality of target coverage achievable with VHEET was comparable to or sometimes better than that provided by IMXT. However, VHEET provided an improvement over IMXT in the dose sparing of the sensitive structures and normal tissues. Compared to IMXT, VHEET decreased the mean rectal dose and bladder dose by up to 10% of the prescribed target dose, while reducing by up to 12% of the prescribed target dose the integral dose to normal tissues. In quantifying the merits of IMPT relative to IMXT, it was found that using intensity-modulated proton beams for inverse planning instead of intensity-modulated photon beams improved target dose homogeneity by up to 1.3% of the prescribed target dose, while reducing the mean rectal dose, bladder dose, and normal tissue integral dose by up to 27%, 30% and 28% of the prescribed target dose respectively. The comparison of optimized planning for IMPT and VHEET showed that the quality of target coverage achievable with IMPT is comparable to or better (by up to 1.3% of the prescribed target dose) than that provided by VHEET. Compared to VHEET, IMPT delivered a mean rectal dose and a bladder dose that was lower by up to 17% and 23% of prescribed target dose respectively, and also reduced the integral dose to normal tissues by up to 17% of the prescribed target dose. These results indicate that of the three modalities the greatest dose escalation will be possible with IMPT, then VHEET, and then IMXT. It follows that IMPT will result in the highest probability of complication-free tumour control, while IMXT will provide the lowest probability.

A comparison of physically and radiobiologically based optimization for IMRT

Many optimization techniques for intensity modulated radiotherapy have now been developed. The majority of these techniques including all the commercial systems that are available are based on physical dose methods of assessment. Some techniques have also been based on radiobiological models. None of the radiobiological optimization techniques however have assessed the clinically realistic situation of considering both tumor and normal cells within the target volume. This study considers a ratio-based fluence optimizing technique to compare a dose-based optimization method described previously and two biologically based models. The biologically based methods use the values of equivalent uniform dose calculated for the tumor cells and integral biological effective dose for normal cells. The first biologically based method includes only tumor cells in the target volume while the second considers both tumor and normal cells in the target volume. All three methods achieve good conformation to the target volume. The biologically based optimization without the normal tissue in the target volume shows a high dose region in the center of the target volume while this is reduced when the normal tissues are also considered in the target volume. This effect occurs because the normal tissues in the target volume require the optimization to reduce the dose and therefore limit the maximum dose to that volume.

Optimization of intensity-modulated very high energy (50-250 MeV) electron therapy

This work evaluates the potential of very high energy (50-250 MeV) electron beams for dose conformation and identifies those variables that influence optimized dose distributions for this modality. Intensity-modulated plans for a prostate cancer model were optimized as a function of the importance factors, beam energy and number of energy bins, number of beams, and the beam orientations. A trial-and-error-derived constellation of importance factors for target and sensitive structures to achieve good conformal dose distributions was 500, 50, 10 and I for the target, rectum, bladder and normal tissues respectively. Electron energies greater than 100 MeV were found to be desirable for intensity-modulated very high energy electron therapy (VHEET) of prostate cancer. Plans generated for lower energy beams had relatively poor conformal dose distributions about the target region and delivered high doses to sensitive structures. Fixed angle beam treatments utilizing a large number of fields in the range 9-21 provided acceptable plans. Using more than 21 beams at fixed gantry angles had an insignificant effect on target coverage, but resulted in an increased dose to sensitive structures and an increased normal tissue integral dose. Minor improvements in VHEET plans utilizing a 'small' number (< or =9) of beams may be achieved if, in addition to intensity modulation, energy modulation is implemented using a small number (< or =3) of beam energies separated by 50 to 100 MeV. Rotation therapy provided better target dose homogeneity but unfortunately resulted in increased rectal dose, bladder dose and normal tissue integral dose relative to the 21-field fixed angle treatment plan. Modulation of the beam energy for rotation therapy had no beneficial consequences on the optimized dose distributions. Lastly, selection of beam orientations influenced the optimized treatment plan even when a large number of beams (approximately 15) were employed.

Intensity and energy modulated radiotherapy with proton beams: variables affecting optimal prostate plan

Inverse planning for intensity- and energy-modulated radiotherapy (IEMRT) with proton beams involves the selection of (i) the relative importance factors to control the relative importance of the target and sensitive structures, (ii) an appropriate energy resolution to achieve an acceptable depth modulation, (iii) an appropriate beamlet width to modulate the beam laterally, and (iv) a sufficient number of beams and their orientations. In this article we investigate the influence of these variables on the optimized dose distribution of a simulated prostate cancer IEMRT treatment. Good dose conformation for this prostate case was achieved using a constellation of I factors for the target, rectum, bladder, and normal tissues of 500, 50, 15, and 1, respectively. It was found that for an active beam delivery system, the energy resolution should be selected on the basis of the incident beams' energy spread (sigmaE) and the appropriate energy resolution varied from 1 MeV at sigmaE = 0.0 to 5 MeV at sigmaE= 2.0 MeV. For a passive beam delivery system the value of the appropriate depth resolution for inverse planning may not be critical as long as the value chosen is at least equal to one-half the FWHM of the primary beam Bragg peak. Results indicate that the dose grid element dimension should be equal to or no less than 70% of the beamlet width. For this prostate case, we found that a maximum of three to four beam ports is required since there was no significant advantage to using a larger number of beams. However for a small number (< or = 4) of beams the selection of beam orientations, while having only a minor effect on target coverage, strongly influenced the sensitive structure sparing and normal tissue integral dose.

A method for physically based radiotherapy optimization with intelligent tissue weight determination

An important aspect of the intensity modulated radiotherapy (IMRT) process is that of optimizing the beam intensity profiles. Most such methods use a comparison between the defined dose distribution (including prescription dose and dose limits for critical structures) and the current distribution. The comparison may involve dose differences or dose ratios. This paper investigates four different ratio-based methods, with a hypothetical U-shaped target and cylindrical phantom as the test case. All methods are shown to give satisfactory results. There are differences in the dose distributions produced and these can be related to the formulation of the respective fluence-update methods. In common with other IMRT beam optimization methods the techniques used here include tissue-weighting factors which express the relative importance of achieving the goal dose level or limit. Altering the values of these weights can improve the ability of the optimization algorithm to adhere to the goal dose levels. Typically the weights are adjusted on a "trial and error" basis until the planner is satisfied with the results. In this paper a technique for automatic updating of the weights according to the optimization results is applied to each of the four fluence update methods. With this technique the difference between the current dose level in each tissue and the specified dose constant for that tissue is used to increase the weight should the constraint not be met. It is found that the four methods yield similar results if the weights are updated in this way.

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.

The energy-dependent electron loss model for pencil beam dose kernels

The 'monoenergetic' electron loss model was derived in a previous work to account for pathlength straggling in the Fermi-Eyges pencil beam problem. In this paper, we extend this model to account for energy-loss straggling and secondary knock-on electron transport in order to adequately predict a depth dose curve. To model energy-loss straggling, we use a weighted superposition of a discrete number of monoenergetic pencil beams with different initial energies where electrons travel along the depth-energy characteristics in the continuous slowing down approximation (CSDA). The energy straggling spectrum at depth determines the weighting assigned to each monoenergetic pencil beam. Supplemented by a simple transport model for the secondary knock-on electrons, the 'energy-dependent' electron loss model predicts both lateral and depth dose distributions from the electron pencil beams in good agreement with Monte Carlo calculations and measurements. The calculation of dose distribution from a pencil beam takes 0.2 s on a Pentium III 500 MHz computer. Being computationally fast, the 'energy-dependent' electron loss model can be used for the calculation of 3D energy deposition kernels in dose optimization schemes without using precalculated or measured data.

Energy- and intensity-modulated electron beams for radiotherapy

This work investigates the feasibility of optimizing energy- and intensity-modulated electron beams for radiation therapy. A multileaf collimator (MLC) specially designed for modulated electron radiotherapy (MERT) was investigated both experimentally and by Monte Carlo simulations. An inverse-planning system based on Monte Carlo dose calculations was developed to optimize electron beam energy and intensity to achieve dose conformity for target volumes near the surface. The results showed that an MLC with 5 mm leaf widths could produce complex field shapes for MERT. Electron intra- and inter-leaf leakage had negligible effects on the dose distributions delivered with the MLC, even at shallow depths. Focused leaf ends reduced the electron scattering contributions to the dose compared with straight leaf ends. As anticipated, moving the MLC position toward the patient surface reduced the penumbra significantly. There were significant differences in the beamlet distributions calculated by an analytic 3-D pencil beam algorithm and the Monte Carlo method. The Monte Carlo calculated beamlet distributions were essential to the accuracy of the MERT dose distribution in cases involving large air gaps, oblique incidence and heterogeneous treatment targets (at the tissue-bone and bone-lung interfaces). To demonstrate the potential of MERT for target dose coverage and normal tissue sparing for treatment of superficial targets, treatment plans for a hypothetical treatment were compared using photon beams and MERT.

Electron beam modeling and commissioning for Monte Carlo treatment planning

A hybrid approach for commissioning electron beam Monte Carlo treatment planning systems has been studied. The approach is based on the assumption that accelerators of the same type have very similar electron beam characteristics and the major difference comes from the on-site tuning of the electron incident energy at the exit window. For one type of accelerator, a reference machine can be selected and simulated with the Monte Carlo method. A multiple source model can be built on the full Monte Carlo simulation of the reference beam. When commissioning electron beams from other accelerators of the same type, the energy spectra in the source model are tuned to match the measured dose distributions. A Varian Clinac 2100C accelerator was chosen as the reference machine and a four-source beam model was established based on the Monte Carlo simulations. This simplified beam model can be used to generate Monte Carlo dose distributions accurately (within 2%/2 mm compared to those calculated with full phase space data) for electron beams from the reference machine with various nominal energies, applicator sizes, and SSDs. Three electron beams were commissioned by adjusting the energy spectra in the source model. The dose distributions calculated with the adjusted source model were compared with the dose distributions calculated using the phase space data for these beams. The agreement is within 1% in most of cases and 2% in all situations. This preliminary study has shown the capability of the commissioning approach for handling large variation in the electron incident energy. The possibility of making the approach more versatile is also discussed.

Monte Carlo-based inverse treatment planning

A Monte Carlo based inverse treatment planning system (MCI) has been developed which combines arguably the most accurate dose calculation method (Monte Carlo particle transport) with a 'guaranteed' optimization method (simulated annealing). A distribution of photons is specified in the tumour volume; they are transported using an adjoint calculation method to outside the patient surface to build up an intensity distribution. This intensity distribution is used as the initial input into an optimization algorithm. The dose distribution from each beam element from a number of fields is pre-calculated using Monte Carlo transport. Simulated annealing optimization is then used to find the weighting of each beam element, to yield the optimal dose distribution for the given criteria and constraints. MCI plans have been generated in various theoretical phantoms and patient geometries. These plans show conformation of the dose to the target volume and avoidance of critical structures. To verify the code, an experiment was performed on an anthropomorphic phantom.