Stopping-power ratios for clinical electron beams from a scatter-foil linear accelerator

Relationship between dose and stopping power values for electrons in skin and muscle tissues

Absorbed dose and stopping power of the target material are two important parameters for determining the radiation effects. In this work, the relationship between these two parameters has been investigated for electron beams incident on skin and muscle tissue. Absorbed dose was obtained by using the EGSnrc code and the stopping power values were calculated by considering the velocity-depended effective charge and mean excitation values. To obtain the relationship between absorbed dose and stopping power values, these parameters were graphed together and simple fitting functions have been obtained. The results obtained show that these parameters are linearly correlated with each other.

Dosimetry of a small air cavity for clinical electron beams: A Monte Carlo study

This study investigated dosimetric changes in a water phantom when a small air cavity was presented at the central axis of a clinical electron beam. We used 6-, 9-, and 16-MeV electron beams with a 10 x 10 cm(2) applicator and cutout produced by a Varian 21 EX linear accelerator. Percentage depth doses (PDDs) for different depths (0.5-7 cm), thicknesses (2-10 mm), and widths (1-5 cm) of air cavities were calculated using Monte Carlo simulations (EGSnrc code) validated by film measurements. By comparing PDDs of phantoms with and without the air cavity, it was found that when the depth or thickness of cavity was changed, the PDD curve below the cavity was shifted with a distance equal to the thickness of the cavity. However, when the width of the air cavity was changed, both the PDD curve and its slope within and below the cavity were changed. A larger width of the air cavity resulted in a shallower PDD curve within the cavity. The slope of the PDD curve below the cavity tended towards a value as the width of the air cavity was increased to 3-5 cm for the 6-, 9-, and 16-MeV electron beams. The dependence of the depth dose on the width of the air cavity is a result of the contribution of the electron side scattering in the water surrounding the cavity. The change in depth dose resulting from the presence of an air cavity can cause discrepancies between the calculated and actual dose during radiotherapy, unless the effects of the air cavity are properly characterized during treatment planning. From the dosimetry data in this study, neglecting an air cavity of 1-cm thickness in the build-up region of a 6-MeV electron beam resulted in a delivered dose 10-12% larger than the original prescription. Delivered doses 3% and 6% higher than the prescribed dose were observed when doses were prescribed at R(80) for a 16-MeV electron beam. These results were obtained by neglecting air cavities with thicknesses equal to 2 and 4 mm, respectively, at a depth of 5 cm.

Variations in energy spectra and water-to-material stopping-power ratios in three-dimensional conformal and intensity-modulated photon fields

Because of complex dose distributions and dose gradients that are created in three-dimensional conformal radiotherapy (3D-CRT) and intensity-modulated radiation therapy (IMRT), photon- and electron-energy spectra might change significantly with spatial locations and doses. This study examined variations in photon- and electron-energy spectra in 3D-CRT and IMRT photon fields. The effects of spectral variations on water-to-material stopping-power ratios used in Monte Carlo treatment planning systems and the responses of energy-dependent dosimeters, such as thermoluminescent dosimeters (TLDs) and radiographic films were further studied. The EGSnrc Monte Carlo code was used to simulate megavoltage 3D-CRT and IMRT photon fields. The photon- and electron-energy spectra were calculated in 3D water phantoms and anthropomorphic phantoms based on the fluence scored in voxel grids. We then obtained the water-to-material stopping-power ratios in the local voxels using the Spencer-Attix cavity theory. Changes in the responses of films and TLDs were estimated based on the calculated local energy spectra and published data on the dosimeter energy dependency. Results showed that the photon-energy spectra strongly depended on spatial positions and doses in both the 3D-CRT and IMRT fields. The relative fraction of low-energy photons (< 100 keV) increased inversely with the photon dose in low-dose regions of the fields. A similar but smaller effect was observed for electrons in the phantoms. The maximum variation of the water-to-material stopping-power ratio over the range of calculated dose for both 3D-CRT and IMRT was negligible (< 1.0%) for ICRU tissue, cortical bone, and soft bone and less than 3.6% for dry air and lung. Because of spectral softening at low doses, radiographic films in the phantoms could over-respond to dose by more than 30%, whereas the over-response of TLDs was less than 10%. Thus, spatial variations of the photon- and electron-energy spectra should be considered as important factors in 3D-CRT and IMRT dosimetry.

Comparison of PENELOPE Monte Carlo dose calculations with Fricke dosimeter and ionization chamber measurements in heterogeneous phantoms (18 MeV electron and 12 MV photon beams)

Different measurements of depth-dose curves and dose profiles were performed in heterogeneous phantoms and compared to dose distributions calculated by a Monte Carlo code. These heterogeneous phantoms consisted of lung and/or bone heterogeneities. Irradiations and simulations were carried out for an 18 MeV electron beam and a 12 MV photon beam. Depth-dose curves were measured with Fricke dosimeters and with plane and cylindrical ionization chambers. Dose profiles were measured with a small cylindrical ionization chamber at different depths. The LINAC was modelled using the PENELOPE code and phase space files were used as input data for the calculations of the dose distributions in every simulation. The detectors (Fricke dosimeters and ionization chambers) were not modelled in the geometry. There is generally a good agreement between the measurements and PENELOPE. Some discrepancies exist, near interfaces, between the ionization chamber and PENELOPE due to the attenuation of the lower energy electrons by the wall of the ionization chamber.

A virtual-accelerator-based verification of a Monte Carlo dose calculation algorithm for electron beam treatment planning in homogeneous phantoms

By introducing Monte Carlo (MC) techniques to the verification procedure of dose calculation algorithms in treatment planning systems (TPSs), problems associated with conventional measurements can be avoided and properties that are considered unmeasurable can be studied. The aim of the study is to implement a virtual accelerator, based on MC simulations, to evaluate the performance of a dose calculation algorithm for electron beams in a commercial TPS. The TPS algorithm is MC based and the virtual accelerator is used to study the accuracy of the algorithm in water phantoms. The basic test of the implementation of the virtual accelerator is successful for 6 and 12 MeV (gamma < 1.0, 0.02 Gy/2 mm). For 18 MeV, there are problems in the profile data for some of the applicators, where the TPS underestimates the dose. For fields equipped with patient-specific inserts, the agreement is generally good. The exception is 6 MeV where there are slightly larger deviations. The concept of the virtual accelerator is shown to be feasible and has the potential to be a powerful tool for vendors and users.

Measurements of output factors with different detector types and Monte Carlo calculations of stopping-power ratios for degraded electron beams

The aim of the present study was to investigate three different detector types (a parallel-plate ionization chamber, a p-type silicon diode and a diamond detector) with regard to output factor measurements in degraded electron beams, such as those encountered in small-electron-field radiotherapy and intraoperative radiation therapy (IORT). The Monte Carlo method was used to calculate mass collision stopping-power ratios between water and the different detector materials for these complex electron beams (nominal energies of 6, 12 and 20 MeV). The diamond detector was shown to exhibit excellent properties for output factor measurements in degraded beams and was therefore used as a reference. The diode detector was found to be well suited for practical measurements of output factors, although the water-to-silicon stopping-power ratio was shown to vary slightly with treatment set-up and irradiation depth (especially for lower electron energies). Application of ionization-chamber-based dosimetry, according to international dosimetry protocols, will introduce uncertainties smaller than 0.3% into the output factor determination for conventional IORT beams if the variation of the water-to-air stopping-power ratio is not taken into account. The IORT system at our department includes a 0.3 cm thin plastic scatterer inside the therapeutic beam, which furthermore increases the energy degradation of the electrons. By ignoring the change in the water-to-air stopping-power ratio due to this scatterer, the output factor could be underestimated by up to 1.3%. This was verified by the measurements. In small-electron-beam dosimetry, the water-to-air stopping-power ratio variation with field size could mostly be ignored. For fields with flat lateral dose profiles (>3 x 3 cm2), output factors determined with the ionization chamber were found to be in close agreement with the results of the diamond detector. For smaller field sizes the lateral extension of the ionization chamber hampers its use. We therefore recommend that the readily available silicon diode detector should be used for output factor measurements in complex electron fields.

Dosimetry characteristics of degraded electron beams investigated by Monte Carlo calculations in a setup for intraoperative radiation therapy

Degraded electron beams, as used for intraoperative radiation therapy (IORT) or similar complicated dosimetric situations, have different characteristics compared to conventional electron therapy beams. If international dosimetry protocols are applied in a direct manner to such degraded beams, uncertainties will be introduced in the absorbed dose determination. The Monte Carlo method has been used to verify experimentally determined relative absorbed dose distributions and output factors in an IORT geometry. Monte Carlo generated dose distributions are mostly within +/-2% or +/-2 mm of measured data. The simulated output variation between the IORT cones (relative output factors) are mostly within 2% of measured values. By comparing IORT and conventional electron beam characteristics (e.g. energy spectra, angular distributions and the contributions of different system components to these quantities) limitations and uncertainties of commonly used dosimetric techniques in IORT electron fields are quantified. The intraoperative treatment field contains a larger amount of scattered electrons, which leads to a broader energy spectrum as well as a wider angular distribution of electrons at the phantom surface. The dose from the scattered electrons can contribute up to 40% of the total dose at a depth of dose maximum, compared to approximately 10% for standard beams. A study of the energy spectra at the reference depth reveals that an uncertainty of the order of 1% can be introduced if ionization chamber based dosimetry is used to determine output factors for the investigated IORT system. We recommend that relative absorbed dose distributions and output factors in IORT electron beams and for similar complicated dosimetric situations should be determined with detectors having a small energy and angular dependence (e.g. diamond detectors or p-Si diodes).

Monte Carlo calculation of output factors for circular, rectangular, and square fields of electron accelerators (6-20 MeV)

Monte Carlo (MC) techniques can be used to build a simulation model of an electron accelerator to calculate output factors for electron fields. This can be useful during commissioning of electron beams from a linac and in clinical practice where irregular fields are also encountered. The Monte Carlo code BEAM/EGS4 was used to model electron beams (6-20 MeV) from a Varian 2100C linear accelerator. After optimization of the Monte Carlo simulation model, agreement within 1% to 2% was obtained between calculated and measured (with a Si diode) lateral and depth dose distributions or within 1 mm in the penumbral regions. Output factors for square, rectangular, and circular fields were measured using two different plane-parallel ion chambers (Markus and NACP) and compared to MC simulations. The agreement was usually within 1% to 2%. This study was not primarily concerned with minimizing the simulation time required to obtain output factors but some considerations with respect to this are presented. It would be particularly useful if the MC model could also be used to calculate output factors for other, similar linacs. To see if this was possible, the primary electron energies in the MC model were retuned to model a recently commissioned similar linac. Good agreement between calculated and measured output factors was obtained for most field sizes for this second accelerator.

Quantifying effects of lead shielding in electron beams: a Monte Carlo study

Lead shielding in contact with the patient's skin is often encountered in radiotherapy with electron beams. The influence of the lead shielding on dose distributions in the patient cannot fully be assessed using modern treatment planning systems. In this work the problem of quantifying the effect of lead shielding on dose distributions is addressed. Monte Carlo dose calculations were performed in a half-blocked water phantom shielded by lead, using a realistic model for the fluence of an electron linear accelerator. Electron beam energies of 6-20 MeV and lead thicknesses of 1-7 mm are used for 10 x 10 cm2 and 5 x 5 cm2 fields. The perturbation of the particle fluence and dose distributions in water introduced by the lead shielding is quantified. The effect of oblique electron beams on the dose perturbation is shown. A fictitious clinical example, the shielding of an eye in electron beam treatment, is used to demonstrate the usefulness of Monte Carlo based treatment planning algorithms that can incorporate the effects of lead shielding.

Electron spectra derived from depth dose distributions

The technique of extracting electron energy spectra from measured distributions of dose along the central axis of clinical electron beams is explored in detail. Clinical spectra measured with this simple spectroscopy tool are shown to be sufficient in accuracy and resolution for use in Monte Carlo treatment planning. A set of monoenergetic depth dose curves of appropriate energy spacing, precalculated with Monte Carlo for a simple beam model, are unfolded from the measured depth dose curve. The beam model is comprised of a point electron and photon source placed in vacuum with a source-to-surface distance of 100 cm. Systematic error introduced by this model affects the calculated depth dose curve by no more than 2%/2 mm. The component of the dose due to treatment head bremsstrahlung, subtracted prior to unfolding, is estimated from the thin-target Schiff spectrum within 0.3% of the maximum total dose (from electrons and photons) on the beam axis. Optimal unfolding parameters are chosen, based on physical principles. Unfolding is done with the public-domain code FERDO. Comparisons were made to previously published spectra measured with magnetic spectroscopy and to spectra we calculated with Monte Carlo treatment head simulation. The approach gives smooth spectra with an average resolution for the 27 beams studied of 16+/-3% of the mean peak energy. The mean peak energy of the magnetic spectrometer spectra was calculated within 2% for the AECL T20 scanning beam accelerators, 3% for the Philips SL25 scattering foil based machine. The number of low energy electrons in Monte Carlo spectra is estimated by unfolding with an accuracy of 2%, relative to the total number of electrons in the beam. Central axis depth dose curves calculated from unfolded spectra are within 0.5%/0.5 mm of measured and simulated depth dose curves, except near the practical range, where 1%/1 mm errors are evident.

Photon beam characterization and modelling for Monte Carlo treatment planning

Photon beams of 4, 6 and 15 MV from Varian Clinac 2100C and 2300C/D accelerators were simulated using the EGS4/BEAM code system. The accelerators were modelled as a combination of component modules (CMs) consisting of a target, primary collimator, exit window, flattening filter, monitor chamber, secondary collimator, ring collimator, photon jaws and protection window. A full phase space file was scored directly above the upper photon jaws and analysed using beam data processing software, BEAMDP, to derive the beam characteristics, such as planar fluence, angular distribution, energy spectrum and the fractional contributions of each individual CM. A multiple-source model has been further developed to reconstruct the original phase space. Separate sources were created with accurate source intensity, energy, fluence and angular distributions for the target, primary collimator and flattening filter. Good agreement (within 2%) between the Monte Carlo calculations with the source model and those with the original phase space was achieved in the dose distributions for field sizes of 4 cm x 4 cm to 40 cm x 40 cm at source surface distances (SSDs) of 80-120 cm. The dose distributions in lung and bone heterogeneous phantoms have also been found to be in good agreement (within 2%) for 4, 6 and 15 MV photon beams for various field sizes between the Monte Carlo calculations with the source model and those with the original phase space.

Monte Carlo modelling of electron beams from medical accelerators

Monte Carlo simulation of radiation transport is considered to be one of the most accurate methods of radiation therapy dose calculation. With the rapid development of computer technology, Monte Carlo based treatment planning for radiation therapy is becoming practical. A basic requirement for Monte Carlo treatment planning is a detailed knowledge of the radiation beams from medical accelerators. A practical approach to obtain the above is to perform Monte Carlo simulation of radiation transport in the medical accelerator. Additionally, Monte Carlo modelling of the treatment machine head can also improve our understanding of clinical beam characteristics, help accelerator design and improve the accuracy of clinical dosimetry by providing more realistic beam data. This paper summarizes work over the past two decades on Monte Carlo simulation of clinical electron beams from medical accelerators.