Lateral Dose Profiles Research Articles

Extended collimator model for pencil-beam dose calculation in proton radiotherapy

We have developed a simple collimator model to improve the accuracy of penumbra behaviour in pencil-beam dose calculation for proton radiotherapy. In this model, transmission of particles through a three-dimensionally extended opening of a collimator is calculated in conjunction with phase-space distribution of the particles. Comparison of the dose distributions calculated using the new three-dimensional collimator model and the conventional two-dimensional model to lateral dose profiles experimentally measured with collimated proton beams showed the superiority of the new model over the conventional one.

Determination of dosimetrical quantities used in microbeam radiation therapy (MRT) with Monte Carlo simulations

Microbeam radiation therapy (MRT) is being performed by using an array of narrow rectangular x-ray beams (typical beam sizes 25 microm X 1 cm), positioned close to each other (typically 200 microm separation), to irradiate a target tissue. The ratio of peak-to-valley doses (PVDR's) in the composite dose distribution has been found to be strongly correlated with the normal tissue tolerance and the therapeutic effect of MRT. In this work a Monte Carlo (MC) study of the depth- and lateral-dose profiles in water for single x-ray microbeams of different shapes and energies has been performed with the MC code PENELOPE. The contributions to the dose deposition from different interaction types have been determined at different distances from the center of the microbeam. The dependence of the peak dose, in a water phantom, on the microbeam field size used in the preclinical trials, has been demonstrated. Composite dose distributions for an array of microbeams were obtained using superposition algorithms and PVDR's were determined and compared with literature results obtained with other Monte Carlo codes. The dependence of the PVDR's on microbeam width, x-ray energy used, and on the separation between adjacent microbeams has been studied in detail.

Monte Carlo study of backscatter in a flattening filter free clinical accelerator

In conventional linear accelerators, the flattening filter provides a uniform lateral dose profile. In intensity modulated radiation therapy applications, however, the flatness of the photon field and hence the presence of a flattening filter, is not necessary. Removing the filter may provide some advantages, such as faster treatments and smaller out-of-field doses to the patients. In clinical accelerators the backscattered radiation dose from the collimators must be taken into account when the dose to the target volume in the patient is being determined. In the case of a conventional machine, this backscatter is known to great precision. In a flattening filter free accelerator, however, the amount of backscatter may be different. In this study we determined the backscatter contribution to the monitor chamber signal in a flattening filter free clinical accelerator (Varian Clinac 21EX) with Monte Carlo simulations. We found that with the exception of very small fields in the 18-MV photon mode, the contribution of backscattered radiation to the monitor signal did not differ from that of conventional machines with a flattening filter. Hence, a flattening filter free clinical accelerator would not necessitate a different backscatter correction.

Dosimetric properties of photon beams from a flattening filter free clinical accelerator

Basic dosimetric properties of 6 MV and 18 MV photon beams from a Varian Clinac 21EX accelerator operating without the flattening filter have been measured. These include dose rate data, depth dose dependencies and lateral profiles in a water phantom, total scatter factors and transmission factors of a multileaf collimator. The data are reviewed and compared with measurements for the flattened beams. The unflattened beams have the following: a higher dose rate by factors of 2.3 (6 MV) and 5.5 (18 MV) on the central axis; lower out-of-field dose due to reduced head scatter and softer spectra; less variation of the total scatter factor with field size; and less variation of the shape of lateral dose profiles with depth. The findings suggest that with a flattening filter free accelerator better radiation treatments can be developed, with shorter delivery times and lower doses to normal tissues and organs.

The dose-area product, a new parameter for the dosimetry of narrow photon beams

In the dosimetry of narrow photon fields with side lengths of the order of 1 cm, the traditional parametrisation via the absolute dose on the beam axis and the relative lateral dose distribution has to deal with the difficulty to find sufficiently small detectors and to adjust them accurately on the narrow-beam axis. This can be avoided by reconsidering the parametrisation, using as normalization factor the surface integral of the dose in the plane perpendicular to the beam axis, abbreviated as the “dose-area product” (DAP). We investigated and confirmed the ability of a large-area parallel-plate ionisation chamber, with a sensitive volume shaped as a flat cylinder of 81.6 mm diameter and 2 mm thickness, to perform the integration over the full lateral dose profile of narrow photon beams with side lengths up to 5 cm. The lateral adjustment of this large-area detector relative to a narrow photon beam is not critical. The large-area ionisation chamber was calibrated in terms of the DAP by reference to a 0.3 cm3 ionisation chamber. A field-size dependent “modified output factor” was defined as the ratio of the DAP, measured at 5 cm phantom depth for 100 cm SSD, and the monitor reading. A prominent phenomenon of narrow photon fields is the field-size and source-distance independence of the relative axial profile of the DAP as function of the thickness of a pre-absorber or of the depth in a phantom. For narrow-beam treatment planning in IMRT, the DAP is combined with the energy- and field size-dependent relative lateral dose distribution which is represented, for example, by a Gaussian convolution kernel. Another useful feature of the DAP is the possibility of its direct control during patient irradiation by means of an on-line monitor with spatial resolution, arranged in the accessory holder.

Benchmarking analytical calculations of proton doses in heterogeneous matter

A proton dose computational algorithm, performing an analytical superposition of infinitely narrow proton beamlets (ASPB) is introduced. The algorithm uses the standard pencil beam technique of laterally distributing the central axis broad beam doses according to the Moliere scattering theory extended to slablike varying density media. The purpose of this study was to determine the accuracy of our computational tool by comparing it with experimental and Monte Carlo (MC) simulation data as benchmarks. In the tests, parallel wide beams of protons were scattered in water phantoms containing embedded air and bone materials with simple geometrical forms and spatial dimensions of a few centimeters. For homogeneous water and bone phantoms, the proton doses we calculated with the ASPB algorithm were found very comparable to experimental and MC data. For layered bone slab inhomogeneity in water, the comparison between our analytical calculation and the MC simulation showed reasonable agreement, even when the inhomogeneity was placed at the Bragg peak depth. There also was reasonable agreement for the parallelepiped bone block inhomogeneity placed at various depths, except for cases in which the bone was located in the region of the Bragg peak, when discrepancies were as large as more than 10%. When the inhomogeneity was in the form of abutting air-bone slabs, discrepancies of as much as 8% occurred in the lateral dose profiles on the air cavity side of the phantom. Additionally, the analytical depth-dose calculations disagreed with the MC calculations within 3% of the Bragg peak dose, at the entry and midway depths in the phantom. The distal depth-dose 20%-80% fall-off widths and ranges calculated with our algorithm and the MC simulation were generally within 0.1 cm of agreement. The analytical lateral-dose profile calculations showed smaller (by less than 0.1 cm) 20%-80% penumbra widths and shorter fall-off tails than did those calculated by the MC simulations. Overall, this work validates the usefulness of our ASPB algorithm as a reasonably fast and accurate tool for quality assurance in planning wide beam proton therapy treatment of clinical sites either composed of homogeneous materials or containing laterally extended inhomogeneities that are comparable in density and located away from the Bragg peak depths.

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 × 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.

Influence of initial electron beam characteristics on Monte Carlo calculated absorbed dose distributions for linear accelerator electron beams

The least known parameters in a Monte Carlo simulation of a linear accelerator treatment head are often the properties of the initial electron beam directed onto the exit vacuum window. Several initial beams with different spatial fluence distributions, angular divergences and energy spectra have been transported through the geometry of a scattering foil accelerator. The electron beam characteristics (energy spectrum and angular distribution) at the phantom surface and the subsequent relative absorbed dose distribution in a water phantom were calculated. The dose distribution was found to be insensitive to the geometrical properties of the initial beam. Furthermore, the lateral dose profiles are unaffected by the energy spectrum of the initial beam. The effect on the depth–dose curve is negligible if the initial energy spectrum is symmetric (e.g., Gaussian shaped) and its full width at half maximum (FWHM) is less than approximately 10% of the most probable energy. A larger FWHM will decrease the normalized dose gradient, but will not affect the dose in the build-up region. An asymmetric wedge shaped spectrum with a low-energy extension simultaneously increases the dose in the build-up region and decreases the dose gradient. The relationship between the energy spectral width and the normalized dose gradient is, however, smaller than published analytical expressions indicate. Some well-established energy-range relationships were shown to be accurate for most of the initial beams studied. The energy spectrum at the phantom surface was also derived from a measured depth–dose curve through different methods. The extracted spectrum depends on the beam model and the spectral reconstruction algorithm. Even though the depth–dose curve is fairly independent of initial beam characteristics, a correct description of the low-energy tail of the energy spectrum is important to obtain good agreement between measured and Monte Carlo calculated doses in the build-up region.

Klinische Dosimetrie für schwere geladene Teilchen

Since December 1997, patients are treated with carbon ions at GSI (Gesellschaft für Schwerionenforschung). Dose delivery is performed with the intensity-controlled raster-scanning technique, which allows a highly conformal treatment of the tumor. To meet the special requirements of dosimetry with heavy ion beams, new dosimetric measurement techniques were developed and introduced into clinical application by the DKFZ (Deutsches Krebsforschungszentrum). The techniques comprise calibration of the irradiation monitor, checks of lateral and depth dose profiles, as well as verification of the beam delivery for complex three-dimensional dose distributions. The developed dosimetric methods are now integral part of clinical application and enable safe treatment with carbon ion therapy.

Monte Carlo simulation and dosimetric verification of radiotherapy beam modifiers

Monte Carlo simulation of beam modifiers such as physical wedges and compensating filters has been performed with a rectilinear voxel geometry module. A modified version of the EGS4/DOSXYZ code has been developed for this purpose. The new implementations have been validated against the BEAM Monte Carlo code using its standard component modules (CMs) in several geometrical conditions. No significant disagreements were found within the statistical errors of 0.5% for photons and 2% for electrons. The clinical applicability and flexibility of the new version of the code has been assessed through an extensive verification versus dosimetric data. Both Varian multi-leaf collimator (MLC) wedges and standard wedges have been simulated and compared against experiments for 6 MV photon beams and different field sizes. Good agreement was found between calculated and measured depth doses and lateral dose profiles along both wedged and unwedged directions for different depths and focus-to-surface distances. Furthermore, Monte Carlo-generated output factors for both open and wedged fields agreed with linac commissioning beam data within statistical uncertainties of the calculations (<3% at largest depths). Compensating filters of both low-density and high-density materials have also been successfully simulated. As a demonstration, a wax compensating filter with a complex three-dimensional concave and convex geometry has been modelled through a CT scan import. Calculated depth doses and lateral dose profiles for different field sizes agreed well with experiments. The code was used to investigate the performance of a commercial treatment planning system in designing compensators. Dose distributions in a heterogeneous water phantom emulating the head and neck region were calculated with the convolution–superposition method (pencil beam and collapsed cone implementations) and compared against those from the MC code developed herein. The new technique presented in this work is versatile, DICOM-RT compliant and accurate in the simulation of beam modulators. This paper addresses the need to reduce the sources of error in the modelling of beam modifiers since they remain a viable alternative to the MLC technique in the delivery of IMRT beams.

Measurement of the dose deposition characteristics of x-ray fluoroscopy beams in water.

The purpose of this study was to characterize the x-ray dose distribution of fluoroscopy beams by measuring their percent depth dose curves and lateral dose profiles in a water phantom. Percent depth dose curves were measured near the surface with an Attix parallel plate chamber and deep within the water phantom with a Farmer-type cylindrical chamber. Percent depth dose curves were compared to published data where applicable. Lateral dose profiles were measured at depths of 2, 5, 10, and 15 cm in phantom with a Farmer chamber. Pulsed, 50 mA x-ray beams with peak tube potentials of 60, 80, 100, and 120 kV and half value layers of 1.89, 2.52, 3.20, and 4.09 mm Al, respectively, were investigated.

Monte Carlo modelling of a virtual wedge

Compared with a set of physical photon wedges, a non-physical wedge (virtual or dynamic wedge), realized by a moving collimator jaw, offers an alternative that allows creation of a wedged field with any arbitrary wedge angle instead of the traditional four physical wedges (15°, 30°, 45° and 60°). It is commonly assumed that non-physical wedges do not alter the photon spectrum compared with physical wedges that introduce beam hardening and loss of dose uniformity in the unwedged direction. In this study, we investigated the influence of a virtual wedge on the photon spectra of a 6-10 MV Siemens MD2 accelerator with the Monte Carlo code EGS4/BEAM. Good agreement was obtained between calculated and measured lateral dose profiles at the depth of maximum dose and at 10 cm depth for 20 × 20 cm2 fields for 6 and 10 MV photon beams. By comparing Monte Carlo models of a physical wedge and the virtual wedge that was studied in this work, it is confirmed that the latter has an insignificant effect on the beam quality, whereas the former can introduce significant beam hardening.

Monte Carlo calculation of the dose distributions of two 106Ru eye applicators

Background and purpose: Beta emitting 106Ru applicators are widely used to treat choroidal melanoma. In view of the importance of clinical applications of this radioisotope and the relative lack of knowledge of the dose distributions, three-dimensional dose maps of two concave applicators were calculated by means of Monte Carlo simulation. Materials and methods: Simulations of small CCA and CCB concave applicators manufactured by Bebig were performed using the Monte Carlo code PENELOPE, which allows the description of the structure (geometry and materials) of the applicator in detail. Electrons are emitted from the 106Ru nuclei isotropically, with initial energy randomly sampled from the corresponding Fermi spectra and with initial positions uniformly distributed on the radioactive layer. Primary electrons, as well as the produced δ-rays, are assumed to be absorbed in the medium when they slow down to an energy of 70 keV. Bremsstrahlung photons with energies larger than 7 keV are also simulated. The simulation code has been run on a 166 MHz PENTIUM PC. Results: Three-dimensional dose distributions produced by the CCA and CCB applicators in a water sphere, concentric with the applicator, were evaluated. To minimize the magnitude of statistical uncertainties, advantage has been taken of the cylindrical symmetry of the problem. The relative depth-dose (along the symmetry axis of the applicator) was also evaluated from the applicator surface up to distances larger than 1 cm, with statistical uncertainties of a few percent. Results compare well with data supplied by the manufacturer. Conclusions: We have performed accurate Monte Carlo calculations of three-dimensional dose distributions from CCA and CCB 106Ru applicators. The results, presented in the form of two-dimensional maps, depth-dose distributions along the symmetry axis and lateral dose profiles, provide a detailed description of the dose delivered in treatments of choroidal melanoma.

Monte Carlo dosimetry study of a 6 MV stereotactic radiosurgery unit

Small-field and stereotactic radiosurgery (SRS) dosimetry with radiation detectors, used for clinical practice, have often been questioned due to the lack of lateral electron equilibrium and uncertainty in beam energy. A dosimetry study was performed for a dedicated 6 MV SRS unit, capable of generating circular radiation fields with diameters of 1.25-5 cm at isocentre using the BEAM/EGS4 Monte Carlo code. With this code the accelerator was modelled for radiation fields with a diameter as small as 0.5 cm. The radiation fields and dosimetric characteristics (photon spectra, depth doses, lateral dose profiles and cone factors) in a water phantom were evaluated. The cone factor for a specific cone c at depth d is defined as , where is the reference cone. To verify the Monte Carlo calculations, measurements were performed with detectors commonly used in SRS such as small-volume ion chambers, a diamond detector, TLDs and films. Results show that beam energies vary with cone diameter. For a 6 MV beam, the mean energies in water at the point of maximum dose for a 0.5 cm cone and a 5 cm cone are 2.05 MeV and 1.65 MeV respectively. The values of obtained by the simulations are in good agreement with the results of the measurements for most detectors. When the lateral resolution of the detectors is taken into account, the results agree within a few per cent for most fields and detectors. The calculations showed a variation of with depth in the water. Based on calculated electron spectra in water, the validity of the assumption that measured dose ratios are equal to measured detector readings was verified.