Biomedical subjects

T R Mackie

Publications and source records attributed to T R Mackie.

At least 73 records · Page 4Linked to original sources

PubMed1995

BEAM: a Monte Carlo code to simulate radiotherapy treatment units.

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.

Electrons↗

PubMed1995

Calculation of stopping-power ratios using realistic clinical electron beams.

The Spencer-Attix water/air restricted mass collision stopping-power ratio is calculated in realistic electron beams in the energy range from 5-50 MeV for a variety of clinical accelerators including the Varian Clinac 2100C, the Philips SL75-20, the Siemens KD2, the AECL Therac 20, and the Scanditronix Medical Microtron 50. The realistic clinical beams are obtained from full Monte Carlo simulations of the clinical linear accelerators using the code BEAM. The stopping-power ratios calculated using clinical beams are compared with those determined according to the AAPM and the IAEA protocols which were calculated by using monoenergetic parallel beams. Using the energy-range relationship of Rogers and Bielajew [Med. Phys. 13, 687-694 (1986)] leads to the most consistent picture in which the stopping-power ratios at dmax derived from mono-energetic calculations underestimate the stopping-power ratios calculated with the realistic beam by 0.3% at 5 MeV and up to 1.4% at 20 MeV. The stopping-power ratios at dmax determined according to the AAPM TG-21 protocol (1983) are shown to overestimate the realistic stopping-power ratios by up to 0.6% for a 5-MeV beam and underestimate them by up to 1.2% for a 20-MeV beam. Those determined according to the IAEA (1987) protocol overestimate the realistic stopping-power ratios by up to 0.3% for a 5-MeV beam and underestimate them by up to a 1.1% for a 20-MeV beam at reference depth. The causes of the differences in the stopping-power ratios between the realistic clinical mono-energetic beams are analyzed quantitatively. The changes in the stopping-power ratios at dmax are mainly due to the energy spread of the electron beam and the contaminant photons in the clinical beams. The effect of the angular spread of electrons is rather small except at the surface. Data are presented which give the corrected stopping-power ratios at dmax or reference depth starting from those determined according to protocols for any energy of clinical electron beams with scattering foils. For scanned clinical electron beams the correction to stopping-power ratios determined according to protocols is found to be less than 0.5% at dmax or reference depth for all beam energies studied. We quantify the differences in the stopping-power ratios determined using the depth of 50% ionization level and the depth of 50% dose level. The differences are very small except for very-high-energy beams (50 MeV) where they can be up to 0.8%.

Electrons↗

PubMed1995

Photon dose calculation incorporating explicit electron transport.

Significant advances have been made in recent years to improve photon dose calculation. However, accurate prediction of dose perturbation effects near the interfaces of different media, where charged particle equilibrium is not established, remain unsolved. Furthermore, changes in atomic number, which affect the multiple Coulomb scattering of the secondary electrons, are not accounted for by current photon dose calculation algorithms. As local interface effects are mainly due to the perturbation of secondary electrons, a photon-electron cascade model is proposed which incorporates explicit electron transport in the calculation of the primary photon dose component in heterogeneous media. The primary photon beam is treated as the source of many electron pencil beams. The latter are transported using the Fermi-Eyges theory. The scattered photon dose contribution is calculated with the dose spread array [T.R. Mackie, J.W. Scrimger, and J.J. Battista, Med. Phys. 12, 188-196 (1985)] approach. Comparisons of the calculation with Monte Carlo simulation and TLD measurements show good agreement for positions near the polystyrene-aluminum interfaces.

Algorithms↗

PubMed1994

Measurements of the electron dose distribution near inhomogeneities using a plastic scintillation detector.

PURPOSE: Accurate measurement of the electron dose distribution near an inhomogeneity is difficult with traditional dosimeters which themselves perturb the electron field. We tested the performance of a new high resolution, water-equivalent plastic scintillation detector which has ideal properties for this application. METHODS AND MATERIALS: A plastic scintillation detector with a 1 mm diameter, 3 mm long cylindrical sensitive volume was used to measure the dose distributions behind standard benchmark inhomogeneities in water phantoms. The plastic scintillator material is more water equivalent than polystyrene in terms of its mass collision stopping power and mass scattering power. Measurements were performed for beams of electrons having initial energies of 6 and 18 MeV at depths from 0.2-4.2 cm behind the inhomogeneities. RESULTS: The detector reveals hot and cold spots behind heterogeneities at resolutions equivalent to typical film digitizer spot sizes. Plots of the dose distributions behind air, aluminum, lead, and formulations for cortical and inner bone-equivalent materials are presented. CONCLUSION: The plastic scintillation detector is suited for measuring the electron dose distribution near an inhomogeneity.

Electrodes↗

PubMed1994

Fine tuning of linear accelerator accessories for stereotactic radiotherapy.

PURPOSE: Experience with the University of Wisconsin's stereotactic radiotherapy (SRT) accessory system was applied to build a new system, facilitate alignment of linac photon beams with a Brown-Roberts-Wells (BRW) stereotaxy, and increase the versatility and stability of the stereotaxy. METHODS AND MATERIALS: High tensile strength stainless steel was used in the floor stand to increase the range of gantry rotation relative to ranges allowed by truss-mounted stands. The collimator assembly and floor stand were each fitted with two-axis gimbal and translation adjustments in addition to the floor stand's three-axis adjustments. The head ring positioning assembly was fitted with two braces to prevent the head ring from deforming with patient motion. Six MV linac photon beam characteristics were measured with a computer-controlled scanning system and a diode in water, at source to surface distances (SSD) of 80 and 100 cm, and for 13 divergent collimators ranging in diameter from 1-4 cm at 100 cm SSD. Quality assurance software was applied to screen data for questionable consistency or symmetry. Integrity of the stereotaxy was evaluated with target simulation films and repeated measurements which were part of the quality assurance of clinical treatments. A method was developed using a glass etched contact reticle to obtain average simulated target to beam center distances (delta av) from target simulation films. RESULTS AND CONCLUSION: New aspects of the current system have improved the ability to fine tune and analyze stereotactic alignment. Beam characteristics met stringent output criteria and penumbral widths were the same or narrower than penumbral widths reported elsewhere. The precision of measuring delta av was 0.1 mm, and delta av averaged over 50 target simulation films was 0.7 +/- 0.1 mm. Results suggest that it may be useful to determine delta av from target simulation films with the method described here.

Humans↗

PubMed1994

A comparison of three inverse treatment planning algorithms.

Three published inverse treatment planning algorithms for physical optimization of external beam radiotherapy are compared. All three algorithms attempt to minimize a quadratic objective function of the dose distribution. It is shown that the algorithms are based on the common framework of Newton's method of multi-dimensional function minimization. The approximations used within this framework to obtain the different algorithms are described. The use of these algorithms requires that the number of weights of elemental dose distributions be equal to the number of sample points taken in the dose volume. The primary factor in determining how the algorithms are implemented is the dose computation model. Two of the algorithms use pencil beam dose models and therefore directly optimize individual pencil beam weights, whereas the third algorithm is implemented to optimize groups of pencil beams, each group converging upon a common point. All dose computation models assume that the irradiated medium is homogeneous. It is shown that the two different implementations produce similar results for the simple optimization problem of conforming dose to a convex target shape. Complex optimization problems consisting of non-convex target shapes and dose limiting structures are shown to require a pencil beam optimization method.

Algorithms↗

PubMed1994

A filtered backprojection dose calculation method for inverse treatment planning.

An efficient method of calculating dose distributions in homogeneous media for megavoltage photons is described. The method is similar to filtered backprojection image reconstruction and is based on the analogy between external beam radiotherapy and SPECT image reconstruction. The filtered backprojection dose calculation significantly reduces the computation time for a large number of x-ray beams compared to a conventional convolution dosimetry method. A factor of 20 reduction in computation time is demonstrated for a 2D implementation of the model. The method has proved useful for speeding up an inverse treatment planning algorithm for conformal radiotherapy, and has the potential to be implemented in the reconstruction hardware of a radiotherapy CT simulator. Results of computer simulations based on the model are presented.

Humans↗

PubMed1993

The application of correlated sampling to the computation of electron beam dose distributions in heterogeneous phantoms using the Monte Carlo method.

Although the Monte Carlo method is capable of computing the dose distribution in heterogeneous phantoms directly, there are some advantages to computing a heterogeneity correction factor. If this approach is adopted there are savings in time using correlated sampling. This technique forces histories to have the same energy, position, direction and random number seed as incident on both the heterogeneous and homogeneous water phantom. This ensures that a history that has, by chance, travelled through only water in the heterogeneous phantom will have the same path as it would have through the homogeneous phantom, resulting in a reduced variance when a ratio of heterogeneous dose to homogeneous dose is formed. Metrics to describe the distributions of uncertainty, efficiency, and degree of correlation are defined. EGS4 Monte Carlo calculation of the dose distribution from a 20 MeV electron beam on water phantoms containing aluminum or air slab heterogeneities illustrate that this technique is the most efficient when the heterogeneity is deep within the phantom, but that improved efficiency can be realized even when the heterogeneity is at or near the surface. This is because some correlation between the two histories is retained despite passage through the heterogeneity.

Electrons↗

PubMed1992

Defining the role of radiosurgery in the management of brain metastases.

The role of stereotactic radiosurgery in the management of recurrent and newly diagnosed brain metastases was evaluated prospectively. From December 1988 to March 1991, 58 lesions in 40 patients were treated with accelerator-based stereotactic radiosurgery. All patients were followed for a minimum of 6 months or to death. The primary purpose was to determine the impact of radiosurgery on local control and its subsequent effects on quality of life. An overall tumor control rate of 82% with a complete response rate of 43% were achieved. As anticipated, the response rate for smaller tumors was substantially better than that for larger tumors (78% for lesions < 2 cm3; 50% for lesions > or = 10 cm3). Although the overall in-field progression rate was 18.5%, only 1/23 (4%) complete responders subsequently recurred. The in-field failure rate is highly comparable with recently published surgical data. Progression outside the brain was noted in two-thirds of patients. One quarter of the deaths were neurologic. The median survival for this minimally selected patient population was 6.5 months. Stereotactic radiosurgery was also associated with improved quality of life as measured by Karnofsky score, neurologic function, and steroid dependence. Long-term steroid dependence was encountered in only four patients. We conclude that stereotactic radiosurgery can be used effectively in patients with brain metastases. In this series, a high tumor response rate was achieved which was associated with improved quality of life.

Adult↗

PubMed1992

Water-equivalent plastic scintillation detectors for high-energy beam dosimetry: I. Physical characteristics and theoretical consideration.

A minimally perturbing plastic scintillation detector has been developed for the dosimetry of high-energy beams in radiotherapy. The detector system consists of two identical parallel sets of radiation-resistant optical fibre bundles, each connected to independent photomultiplier tubes (PMTs). One fibre bundle is connected to a miniature water equivalent plastic scintillator and so scintillation as well as Cerenkov light generated in the fibres is detected at its PMT. The other 'background' bundle is not connected to the scintillator and so only Cerenkov light is detected by its PMT. The background signal is subtracted to yield only the signal from the scintillator. The water-equivalence of plastic scintillation detectors is studied for photon and electron beams in the radiotherapy range. Application of Burlin cavity theory shows that the energy dependence of such detectors is expected to be better than the commonly used systems (ionization chambers, LiF thermoluminescent dosimeters, film and Si diodes). It is also shown that they are not affected by temperature variations and exhibit much less radiation damage than either photon or electron diode detectors.

Fiber Optic Technology↗

PubMed1992

Water-equivalent plastic scintillation detectors for high-energy beam dosimetry: II. Properties and measurements.

The properties of a new scintillation detector system for use in dosimetry of high-energy beams in radiotherapy have been measured. The most important properties of these detectors are their hgh spatial resolution and their nearly water-equivalence. Measurements have shown that they have excellent reproducibility and stability, and a linear response versus dose-rate. It is also shown that they have better spatial resolution than ionization chambers and have much less energy or depth dependence in electron fields due to the removal of the influence of the polarization effect. Dose distributions in water, using miniature plastic scintillation detectors, have been measured for different high-energy photon and electron beams.

Fiber Optic Technology↗

PubMed1991

A unified approach to the optimization of brachytherapy and external beam dosimetry.

Semi-automated optimization of dose distributions is possible using techniques borrowed from imaging science. The ideal distribution of dose is first deconvolved by a convolution kernel yielding an ideal weighting distribution in the patient. The weighting distribution describes the total energy released per unit mass of the irradiated medium. For internal and external radiation sources, this is directly related to the amount and distribution of radioactivity and energy fluence in the medium, respectively. For external sources, the exponential Radon transform is used to obtain ideal fluence projections incident on the patient. In both instances negative values are produced, which when set to zero result in perturbed dose distributions. This may necessitate iterative techniques to reduce the 'residual dose' produced by the zeroing process. Application of the approach is presented for the optimization of 2-dimensional dose distributions in external beam therapy, radioimmunotherapy, and brachytherapy sources.

Brachytherapy↗

PubMed1991

A three-dimensional volume visualization package applied to stereotactic radiosurgery treatment planning.

A comprehensive software package has been developed for visualization and analysis of 3-dimensional data sets. The system offers a variety of 2- and 3-dimensional display facilities including highly realistic volume rendered images generated directly from the data set. The package has been specifically modified and successfully used for stereotactic radiosurgery treatment planning. The stereotactic coordinate transformation is determined by finding the localization frame automatically in the CT volume. Treatment arcs are specified interactively and displayed as paths on 3-dimensional anatomical surfaces. The resulting dose distribution is displayed using traditional 2-dimensional displays or as an isodose surface composited with underlying anatomy and the target volume. Dose volume histogram analysis is an integral part of the system. This paper gives an overview of volume rendering methods and describes the application of these tools to stereotactic radiosurgery treatment planning.

Brain Diseases↗

PubMed1991

The dose distribution surrounding 192Ir and 137Cs seed sources.

Dose distributions in water have been measured using LiF thermoluminescent dosimeters for 192Ir seed sources with stainless steel and with platinum encapsulation to determine the effect of differing encapsulation. The dose distribution has also been measured for a 137Cs seed source. In addition, dose distributions surrounding these sources have been calculated using the EGS4 Monte Carlo code and have been compared to the measured data. The two methods are in good agreement for all three sources. Tables are given which describe the dose distribution surrounding each source as a function of distance and angle. In addition, specific dose constants have been determined from results of Monte Carlo simulation. This work has confirmed the utility of the EGS4 Monte Carlo code in modelling 192Ir and 137Cs seed sources to obtain brachytherapy dose distributions.

Brachytherapy↗

PubMed1991

Effect of source encapsulation on the energy spectra of 192Ir and 137Cs seed sources.

The effect of source encapsulation on the energy spectra of 192Ir and 137Cs seed sources, both with stainless steel and with platinum encapsulation, was determined from results of Monte Carlo simulation. The fractional scatter dose around these sources has also been determined from Monte Carlo simulation. The platinum-encapsulated 192Ir source exhibited greater attenuation of the primary spectrum, as expected, and, consistent with this greater attenuation, exhibited more scattered radiation. Significantly less scatter was seen with the 137Cs source than with either 192Ir source, as is consistent with the higher-energy photons from 137Cs.

Brachytherapy↗

PubMed1990

Monte Carlo and convolution dosimetry for stereotactic radiosurgery.

The dosimetry of small photon beams used for stereotactic radiosurgery was investigated using Monte Carlo simulation, convolution calculations, and measurements. A Monte Carlo code was used to simulate radiation transport through a linear accelerator to produce and score energy spectrum and angular distribution of 6 MV bremsstrahlung photons exiting from the accelerator treatment head. These photons were then transported through a stereotactic collimator system and into a water phantom placed at isocenter. The energy spectrum was also used as input for the convolution method of photon dose calculation. Monte Carlo and convolution results were compared with the measured data obtained using an ionization chamber, a diode, and film.

Brain Neoplasms↗

PubMed1988

Generation of photon energy deposition kernels using the EGS Monte Carlo code.

The EGS Monte Carlo code was used to generate photon energy deposition kernels which describe the energy deposited by charged particles set in motion by primary, first scattered, second scattered, multiple scattered and bremsstrahlung plus annihilation photons. These were calculated for a water medium irradiated with monoenergetic photons with energies in the range 0.1-50 MeV. In addition to the primary energy deposition kernels, primary charged particle transport was further characterised by computing the effective centre of the voxels, and the effective penetration depth, effective radius and effective lateral distance travelled by these particles. The dose per unit collision kerma for parallel monoenergetic primary photons beta' was calculated. Additional applications of the energy deposition kernels are discussed.

Computer Simulation↗

PubMed1982

Contamination of a 15-MV photon beam by electrons and scattered photons.

The 15-MV photon beam of a linear accelerator (Siemens Mevatron 20) was studied for electron and scattered photon contamination. The surface dose, attributable almost entirely to contamination electrons, has a Gaussian lateral distribution, a linear dependence on field width for square fields, and an inverse square dependence on distance from the bottom of the fixed head assembly. This geometrical dependence is consistent with the proposal that the field flattening filter is the main source of electron contamination when accessories are absent. A tissue-maximum-ratio curve in the build-up region for the electron and photon contamination was produced utilizing the linearity of dose with respect to field width. The derived contamination curve inside was similar to the measured build-up curve outside the field. The primary photon component, obtained by subtracting the contaminant contribution, showed no dependence on field size, source-to-probe distance, or presence of accessories.

Electrons↗