This paper describes BEAM, a general purpose Monte Carlo code to simulate the radiation beams from radiotherapy units including high-energy electron and photon beams, 60Co beams and orthovoltage units. The code handles a variety of elementary geometric entities which the user puts together as needed (jaws, applicators, stacked cones, mirrors, etc.), thus allowing simulation of a wide variety of accelerators. The code is not restricted to cylindrical symmetry. It incorporates a variety of powerful variance reduction techniques such as range rejection, bremsstrahlung splitting and forcing photon interactions. The code allows direct calculation of charge in the monitor ion chamber. It has the capability of keeping track of each particle's history and using this information to score separate dose components (e.g., to determine the dose from electrons scattering off the applicator). The paper presents a variety of calculated results to demonstrate the code's capabilities. The calculated dose distributions in a water phantom irradiated by electron beams from the NRC 35 MeV research accelerator, a Varian Clinac 2100C, a Philips SL75-20, an AECL Therac 20 and a Scanditronix MM50 are all shown to be in good agreement with measurements at the 2 to 3% level. Eighteen electron spectra from four different commercial accelerators are presented and various aspects of the electron beams from a Clinac 2100C are discussed. Timing requirements and selection of parameters for the Monte Carlo calculations are discussed.

Tomotherapy, literally "slice therapy," is a proposal for the delivery of radiation therapy with intensity-modulated strips of radiation. The proposed method employs a linear accelerator, or another radiation-emitting device, which would be mounted on a ring gantry like a CT scanner. The patient would move through the bore of the gantry simultaneously with gantry rotation. The intensity modulation would be performed by temporally modulated multiple independent leaves that open and close across the slit opening. At any given time, any leaf would be (1) closed, covering a portion of the slit, (2) open, allowing radiation through, or (3) changing between these states. This method would result in the delivery of highly conformal radiation. Overall treatment times should be comparable with contemporary treatment delivery times. The ring gantry would make it convenient to mount a narrow multisegmented megavoltage detector system for beam verification and a CT scanner on the treatment unit. Such a treatment unit could become a powerful tool for treatment planning, conformal treatment, and verification using tomographic images. The physical properties of this treatment delivery are evaluated and the fundamental design specifications are justified.

Monte Carlo studies of dose distributions in patients treated with radiotherapy electron beams would benefit from generalized models of clinical beams if such models introduce little error into the dose calculations. Methodology is presented for the design of beam models, including their evaluation in terms of how well they preserve the character of the clinical beam, and the effect of the beam models on the accuracy of dose distributions calculated with Monte Carlo. This methodology has been used to design beam models for electron beams from two linear accelerators, with either a scanned beam or a scattered beam. Monte Carlo simulations of the accelerator heads are done in which a record is kept of the particle phase-space, including the charge, energy, direction, and position of every particle that emerges from the treatment head, along with a tag regarding the details of the particle history. The character of the simulated beams are studied in detail and used to design various beam models from a simple point source to a sophisticated multiple-source model which treats particles from different parts of a linear accelerator as from different sub-sources. Dose distributions calculated using both the phase-space data and the multiple-source model agree within 2%, demonstrating that the model is adequate for the purpose of Monte Carlo treatment planning for the beams studied. Benefits of the beam models over phase-space data for dose calculation are shown to include shorter computation time in the treatment head simulation and a smaller disk space requirement, both of which impact on the clinical utility of Monte Carlo treatment planning.

A realistic model of photon beams generated by clinical linear accelerators has been incorporated in a convolution/superposition method to compute dose distributions in photon treatment fields. In this beam model, a primary photon source represents photons directly from the target, and an extra-focal photon source represents scattered photons from the primary collimator and the flattening filter. Monte Carlo simulation was used to study clinical linear accelerators producing photon beams. From the output of the Monte Carlo simulation, the fluence and spectral distributions of each photon component, as well as the geometrical characteristics of each photon source with respect to its distance to the isocenter and its source distribution, were analyzed. These quantities were used to reproduce realistic photon distributions in treatment fields, and thus to compute dose distributions using the convolution method. Our results showed that compared to the primary photon fluence, the extra-focal photon fluence from the primary collimator and the flattening filter was 11%-16% at the isocenter, among which 70% was contributed by the flattening filter. The variation of extra-focal photons in different treatment fields was predicted accurately by accounting for the finite size of the extra-focal source. Compared to measurements, dose distributions in photon treatment fields, including those of asymmetric jaw settings and at different SSDs were calculated accurately, particularly in the penumbral region, by using the convolution method with the new dual source photon beam model.

No abstract