PubMed Health⌕ Search

Biomedical subjects

A F Bielajew

Publications and source records attributed to A F Bielajew.

At least 19 recordsLinked to original sources
PubMed2006

A proposed alternative to phase-space recycling using the adaptive kernel density estimator method.

We have implemented a nonparametric density estimation technique, the adaptive kernel density estimator (AKDE), to generate additional phase space (PS) variables in the vicinity of simulated PS points in Monte Carlo linear accelerator simulation. The method involves the placement of kernels at simulated PS points that have a "window width" that depends on the density of simulated PS points. This method has been tested on known one-dimensional (1-D) and two-dimensional (2-D) probability density functions (PDFs) and has been used to sample (photons only) from PS files generated from accelerator simulations. The original simulated PS vector (x, y, u, v, E) was reduced to a rotationally invariant PS vector (r, theta, alpha, E) that takes advantage of the azimuthal symmetry (phi) above the collimating jaws. The new PS vector (r', theta', alpha', E') is sampled in the vicinity of the sampled PS vector (r, theta, alpha, E). The first step in assessing the accuracy of the method was a correlation analysis among the AKDE generated PS variables compared with correlations among the original PS variables. "In-air" particle fluence distributions between AKDE samples and the original PS distribution showed agreement within 2% (-8.8% to 6.8%) across the entire phase space plane. Central axis energy distributions and angular distributions agreed on average to within 1.5% (range = -1.5% to 6.6%) and 0.1% (range = 0 to 3.0%), respectively. Dose profiles were calculated for field sizes 3 x 3 cm2, 10 x 10 cm2, and 30 x 30 cm2 for AKDE and compared against calculations performed with PS recycling. AKDE calculated depth doses and profiles were within 2% and 2%/1 mm, respectively, of those computed using PS recycling.

Algorithms↗
PubMed2001

An apparatus for applying strong longitudinal magnetic fields to clinical photon and electron beams.

Monte Carlo studies have recently renewed interest in the use of the effect of strong transverse and longitudinal magnetic fields to manipulate the dose characteristics of clinical photon and electron beams. A 3.5 T superconducting solenoidal magnet was used to evaluate the effect of a longitudinal field on both photon and electron beams. This note describes the apparatus and demonstrates some of the effects on the beam trajectory and dose distributions for measurements in a homogeneous phantom. The effects were studied using film in air and in phantoms which fit in the magnet bore. The magnetic field focused and collimated the electron beams. The converging, non-uniform field confined the beam and caused it to converge with increasing depth in the phantom. Due to the field's collecting and focusing effect, the beam flux density increased, leading to increased dose deposition near the magnetic axis, especially near the surface of the phantom. This study illustrates some benefits and challenges associated with the use of non-uniform longitudinal magnetic fields in conjunction with clinical electron and photon beams.

Electrons↗
PubMed2001

Description and dosimetric verification of the PEREGRINE Monte Carlo dose calculation system for photon beams incident on a water phantom.

PEREGRINE is a three-dimensional Monte Carlo dose calculation system written specifically for radiotherapy. This paper describes the implementation and overall dosimetric accuracy of PEREGRINE physics algorithms, beam model, and beam commissioning procedure. Particle-interaction data, tracking geometries, scoring, variance reduction, and statistical analysis are described. The BEAM code system is used to model the treatment-independent accelerator head, resulting in the identification of primary and scattered photon sources and an electron contaminant source. The magnitude of the electron source is increased to improve agreement with measurements in the buildup region in the largest fields. Published measurements provide an estimate of backscatter on monitor chamber response. Commissioning consists of selecting the electron beam energy, determining the scale factor that defines dose per monitor unit, and describing treatment-dependent beam modifiers. We compare calculations with measurements in a water phantom for open fields, wedges, blocks, and a multileaf collimator for 6 and 18 MV Varian Clinac 2100C photon beams. All calculations are reported as dose per monitor unit. Aside from backscatter estimates, no additional, field-specific normalization is included in comparisons with measurements. Maximum discrepancies were less than either 2% of the maximum dose or 1.2 mm in isodose position for all field sizes and beam modifiers.

Algorithms↗
PubMed2000

Towards the elimination of Monte Carlo statistical fluctuation from dose volume histograms for radiotherapy treatment planning.

The Monte Carlo calculation of dose for radiotherapy treatment planning purposes introduces unavoidable statistical noise into the prediction of dose in a given volume element (voxel). When the doses in these voxels are summed to produce dose volume histograms (DVHs), this noise translates into a broadening of differential DVHs and correspondingly flatter DVHs. A brute force approach would entail calculating dose for long periods of time--enough to ensure that the DVHs had converged. In this paper we introduce an approach for deconvolving the statistical noise from DVHs, thereby obtaining estimates for converged DVHs obtained about 100 times faster than the brute force approach described above. There are two important implications of this work: (a) decisions based upon DVHs may be made much more economically using the new approach and (b) inverse treatment planning or optimization methods may employ Monte Carlo dose calculations at all stages of the iterative procedure since the prohibitive cost of Monte Carlo calculations at the intermediate calculation steps can be practically eliminated.

Algorithms↗
PubMed2000

150-250 meV electron beams in radiation therapy.

High-energy electron beams in the range 150-250 MeV are studied to evaluate the feasibility for radiotherapy. Monte Carlo simulation results from the PENELOPE code are presented and used to determine lateral spread and penetration of these beams. It is shown that the penumbra is comparable to photon beams at depths less than 10 cm and the practical range (Rp) of these beams is greater than 40 cm. The depth dose distribution of electron beams compares favourably with photon beams. Effects caused by nuclear reactions are evaluated, including increased dose due to neutron production and induced radioactivity resulting in an increased relative biological effectiveness (RBE) factor of < 1.03.

Algorithms↗
PubMed2000

DPM, a fast, accurate Monte Carlo code optimized for photon and electron radiotherapy treatment planning dose calculations.

A new Monte Carlo (MC) algorithm, the 'dose planning method' (DPM), and its associated computer program for simulating the transport of electrons and photons in radiotherapy class problems employing primary electron beams, is presented. DPM is intended to be a high accuracy MC alternative to the current generation of treatment planning codes which rely on analytical algorithms based on an approximate solution of the photon/electron Boltzmann transport equation. For primary electron beams, DPM is capable of computing 3D dose distributions (in 1 mm3 voxels) which agree to within 1% in dose maximum with widely used and exhaustively benchmarked general-purpose public-domain MC codes in only a fraction of the CPU time. A representative problem, the simulation of 1 million 10 MeV electrons impinging upon a water phantom of 128(3) voxels of 1 mm on a side, can be performed by DPM in roughly 3 min on a modern desktop workstation. DPM achieves this performance by employing transport mechanics and electron multiple scattering distribution functions which have been derived to permit long transport steps (of the order of 5 mm) which can cross heterogeneity boundaries. The underlying algorithm is a 'mixed' class simulation scheme, with differential cross sections for hard inelastic collisions and bremsstrahlung events described in an approximate manner to simplify their sampling. The continuous energy loss approximation is employed for energy losses below some predefined thresholds, and photon transport (including Compton, photoelectric absorption and pair production) is simulated in an analogue manner. The delta-scattering method (Woodcock tracking) is adopted to minimize the computational costs of transporting photons across voxels.

Algorithms↗
PubMed1998

A thin target approach for portal imaging in medical accelerators.

A new thin-target method (patent pending) is described for portal imaging with low-energy (tens of keV) photons from a medical linear accelerator operating in a special mode. Low-energy photons are usually produced in the accelerator target, but are absorbed by the target and flattening filter, both made of medium- or high-Z materials such as Cu or W. Since the main contributor to absorption of the low-energy photons is self-absorption by the thick target through the photoelectric effect, it is proposed to lower the thickness of the portal imaging target to the minimum required to get the maximum low-energy photon fluence on the exit side of the target, and to lower the atomic number of the target so that predominantly photoelectric absorption is reduced. To determine the minimum thickness of the target, EGS4 Monte Carlo calculations were performed. As a result of these calculations, it was concluded that the maximum photon fluence for a 4 MeV electron beam is obtained with a 1.5 mm Cu target. This value is approximately five times less than the thickness of the Cu target routinely used for bremsstrahlung production in radiotherapeutic practice. Two sets of experiments were performed: the first with a 1.5 mm Cu target and the second with a 5 mm Al target (Cu mass equivalent) installed in the linear accelerator. Portal films were taken with a Rando anthropomorphic phantom. To emphasize the low-energy response of the new thin target we used a Kodak Min-R mammographic film and cassette combination, with a strong low-energy response. Because of its high sensitivity, only 1 cGy is required. The new portal images show a remarkable improvement in sharpness and contrast in anatomical detail compared with existing ones. It is also shown that further lowering of the target's atomic number (for example to C or Be) produces no significant improvement.

Algorithms↗
PubMed1997

Fluence non-uniformity effects in air kerma determination around brachytherapy sources.

In ionization chamber dosimetry close to brachytherapy sources, the non-uniform photon fluence in the vicinity of the air cavity of the ionization chamber must be corrected for. In this paper, experimentally determined relative non-uniformity correction factors are determined for a PTW23331 chamber and a number of Farmer-type chambers. It is concluded that the 'anisotropic' theory of Bielajew agrees better with the experimental results than the 'isotropic' theory of Kondo and Randolph. The experiments show that neither the material choice nor the radius of the central electrode have any significant effect on the non-uniformity correction factors. Experiments on the wall-material dependence in the non-uniformity correction factor, which is predicted by the anisotropic theory, were inconclusive. Parameters for the Farmer-type chambers for calculation of theoretical non-uniformity correction factors are given. Data to derive anisotropic fluence non-uniformity correction factors are given in tabular form. An error in the Kondo and Randolph equations is reported and the corrected equations presented.

Biophysical Phenomena↗
PubMed1995

An improved energy-range relationship for high-energy electron beams based on multiple accurate experimental and Monte Carlo data sets.

A theoretically based analytical energy-range relationship has been developed and calibrated against well established experimental and Monte Carlo calculated energy-range data. Only published experimental data with a clear statement of accuracy and method of evaluation have been used. Besides published experimental range data for different uniform media, new accurate experimental data on the practical range of high-energy electron beams in water for the energy range 10-50 MeV from accurately calibrated racetrack microtrons have been used. Largely due to the simultaneous pooling of accurate experimental and Monte Carlo data for different materials, the fit has resulted in an increased accuracy of the resultant energy-range relationship, particularly at high energies. Up to date Monte Carlo data from the latest versions of the codes ITS3 and EGS4 for absorbers of atomic numbers between four and 92 (Be, C, H2O, PMMA, Al, Cu, Ag, Pb and U) and incident electron energies between 1 and 100 MeV have been used as a complement where experimental data are sparse or missing. The standard deviation of the experimental data relative to the new relation is slightly larger than that of the Monte Carlo data. This is partly due to the fact that theoretically based stopping and scattering cross-sections are used both to account for the material dependence of the analytical energy-range formula and to calculate ranges with the Monte Carlo programs. For water the deviation from the traditional energy-range relation of ICRU Report 35 is only 0.5% at 20 MeV but as high as -2.2% at 50 MeV. An improved method for divergence and ionization correction in high-energy electron beams has also been developed to enable use of a wider range of experimental results.

Electrons↗
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↗
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↗
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↗
PubMedDate not supplied

Theoretical and experimental investigation of dose enhancement due to charge storage in electron-irradiated phantoms.

Recent measurements have shown that significant errors in radiation dosimetry can arise by the use of insulating plastic phantoms which have been exposed to electron beams. The effect has been attributed to the generation of large electric fields in the phantom by charge storage causing alteration of electron trajectories and an increase in the measured dose. In this report, we examine this hypothesis theoretically by calculating the change in response to radiation of an ion chamber in a cylindrical cavity in an electron-irradiated polymethylmethacrylate phantom. The electric field distribution is determined using a model which allows for charge leakage by radiation-induced conductivity, and the dose in the cavity is determined by a Monte Carlo simulation using the EGS (electron gamma shower) code modified to account for electron trajectories in the electric field. The theoretical results are shown to agree well with new and previously published experimental dose enhancement data. The agreement is taken as confirmation of the reported explanation of the effect. The use of conducting phantoms in radiation dosimetry is advocated.

Electric Conductivity↗
PubMedDate not supplied

Calculated buildup curves for photons with energies up to 60Co.

This report investigates the buildup region of the depth-dose curve for broad parallel beams of 60Co and lower energy photons incident on slabs of tissue as defined by the International Commission on Radiation Units and Measurements (ICRU) and other materials. The EGS Monte Carlo system was used to simulate the transport of photons and electrons in the materials. Calculated buildup regions of the depth-dose curves are presented for photons incident on ICRU tissue. For 60Co beams, the primary and scatter dose components are presented, both for the buildup region and at depth in the phantom. Effects of beam radius were investigated and found to have only a small effect in the buildup region of the central axis depth-dose curve. The calculated effects of beam radius are larger at depth and in good agreement with experiment. Effects of scatter and attenuation in a small miniphantom of water were investigated as a function of photon energy and were found to cancel within 0.8%. The mean distance of electron transport was calculated for 60Co beams incident on carbon, water, polystyrene, aluminum, and lead.

Cobalt Radioisotopes↗
PubMedDate not supplied

The effect of free electrons on ionization chamber saturation curves.

The effect of the collection of free electrons--those which have not formed negative ions in the gas of a parallel-plate ionization chamber--is considered for the case of continuous radiation. A theory and its solution, valid for the entire voltage range, is given. The solution is based upon the saturation theory of Townsend, and free electrons are included in a perturbative fashion. The collection of free electrons becomes important at exposure rates above 430 R/min. Since the electron contribution involves an extra dependence on exposure rate, some of the wide variation in the cavity gas parameter (m value) reported in the literature may be due to the presence of free electrons.

Electrons↗