Jan Seuntjens

Comparison of measured and Monte Carlo calculated dose distributions from the NRC linac

We have benchmarked photon beam simulations with the EGS4 user code BEAM [Rogers et al., Med. Phys. 22, 503-524 (1995)] by comparing calculated and measured relative ionization distributions in water from the 10 and 20 MV photon beams of the NRC linac. Unlike previous calculations, the incident electron energy is known independently to 1%, the entire extra-focal radiation is simulated, and electron contamination is accounted for. The full Monte Carlo simulation of the linac includes the electron exit window, target, flattening filter, monitor chambers, collimators, as well as the PMMA walls of the water phantom. Dose distributions are calculated using a modified version of the EGS4 user code DOSXYZ which additionally allows scoring of average energy and energy fluence in the phantom. Dose is converted to ionization by accounting for the (L/rho)water(air) variation in the phantom, calculated in an identical geometry for the realistic beams using a new EGS4 user code, SPRXYZ. The variation of (L/rho)water(air) with depth is a 1.25% correction at 10 MV and a 2% correction at 20 MV. At both energies, the calculated and the measured values of ionization on the central axis in the buildup region agree within 1% of maximum ionization relative to the ionization at 10 cm depth. The agreement is well within statistics elsewhere. The electron contamination contributes 0.35(+/- 0.02) to 1.37(+/- 0.03)% of the maximum dose in the buildup region at 10 MV and 0.26(+/- 0.03) to 3.14(+/- 0.07)% of the maximum dose at 20 MV. The penumbrae at 3 depths in each beam (in g/cm2), 1.99 (dmax, 10 MV only), 3.29 (dmax, 20 MV only), 9.79 and 19.79, agree with ionization chamber measurements to better than 1 mm. Possible causes for the discrepancy between calculations and measurements are analyzed and discussed in detail.

Fricke dosimetry: the difference between G(Fe3+) for 60Co gamma-rays and high-energy x-rays.

A calibration of the Fricke dosimeter is a measurement of epsilon G(Fe3+). Although G(Fe3+) is expected to be approximately energy independent for all low-LET radiation, existing data are not adequate to rule out the possibility of changes of a few per cent with beam quality. When a high-precision Fricke dosimeter, which has been calibrated for one particular low-LET beam quality, is used to measure the absorbed dose for another low-LET beam quality, the accuracy of the absorbed dose measurement is limited by the uncertainty in the value of G(Fe3+). The ratio of G(Fe3+) for high-energy x-rays (20 and 30 MV) to G(Fe3+) for 60Co gamma-rays, G(Fe3+)MV(Co), was measured to be 1.007(+/-0.003) (confidence level of 68%) using two different types of water calorimeter, a stirred-water calorimeter (20 MV) and a sealed-water calorimeter (20, 30 MV). This value is consistent with our calculations based on the LET dependence of G(primary products) and, as well, with published measurements and theoretical treatments of G(Fe3+).

Monte Carlo study of correction factors for Spencer–Attix cavity theory at photon energies at or above 100 keV

To develop a primary standard for 192Ir sources, the basic science on which this standard is based, i.e., Spencer-Attix cavity theory, must be established. In the present study Monte Carlo techniques are used to investigate the accuracy of this cavity theory for photons in the energy range from 20 to 1300 keV, since it is usually not applied at energies below that of 137Cs. Ma and Nahum [Phys. Med. Biol. 36, 413-428 (1991)] found that in low-energy photon beams the contribution from electrons caused by photons interacting in the cavity is substantial. For the average energy of the 192Ir spectrum they found a departure from Bragg-Gray conditions of up to 3% caused by photon interactions in the cavity. When Monte Carlo is used to calculate the response of a graphite ion chamber to an encapsulated 192Ir source it is found that it differs by less than 0.3% from the value predicted by Spencer-Attix cavity theory. Based on these Monte Carlo calculations, for cavities in graphite it is concluded that the Spencer-Attix cavity theory with delta = 10 keV is applicable within 0.5% for photon energies at 300 keV or above despite the breakdown of the assumption that there is no interaction of photons within the cavity. This means that it is possible to use a graphite ion chamber and Spencer-Attix cavity theory to calibrate an 192Ir source. It is also found that the use of delta related to the mean chord length instead of delta = 10 keV improves the agreement with Spencer-Attix cavity theory at 60Co from 0.2% to within 0.1% of unity. This is at the level of accuracy of which the Monte Carlo code EGSnrc calculates ion chamber responses. In addition, it is shown that the effects of other materials, e.g., insulators and holders, have a substantial effect on the ion chamber response and should be included in the correction factors for a primary standard of air kerma.

Verification of absorbed doses determined with thimble and parallel‐plate ionization chambers in clinical electron beams using ferrous sulphate dosimetry

Absorbed dose values determined with the commonly applied NACP and PTW/Markus parallel-plate chambers and the cylindrical NE2571 Farmer chamber were compared to values obtained with ferrous sulphate dosimetry in a number of electron beams. For the ionometry with the parallel-plate chambers the dose-to-air chamber factor ND (or Ngas) was derived from a 60Co beam calibration free in air with an additional buildup layer of 0.54 g cm-2 graphite as proposed by the protocol for electron dosimetry published by the Netherlands Commission on Radiation Dosimetry. For the product kattkm in this calibration geometry values of 0.980 +/- 0.003 [1 standard deviation (s.d.)] and 0.993 +/- 0.004 (1 s.d.) were obtained for the parallel-plate NACP and PTW/Markus chambers, respectively. The behavior of the fluence perturbation correction factor pf versus the mean electron energy at depth was deduced for the flat PTW/Markus and cylindrical NE2571 chamber by comparison with the NACP chamber, for which pf was assumed unity. Our results show a small but significant energy dependence of pf for the PTW/Markus chamber. The absorbed dose values, determined ionometrically with the different chambers considered in the study using the experimentally determined kattkm and pf values, are systematically 0.5% higher than those obtained with ferrous sulphate dosimetry adopting 352 x 10(-6) m-2 kg-1 Gy-1 for epsilon mG. The performed comparative study confirms also that for the NACP chamber pf is unity independent of the electron energy down to a mean energy at depth of 2 MeV.

Conversion factor f for X-ray beam qualities, specified by peak tube potential and HVL value

Conversion factors of exposure to absorbed dose f in compact bone, muscle, water and adipose tissue are deduced from the measured photon fluence spectra of x-ray beams with potentials from 50 to 250 kV. The results show that accurate determination of f requires beam quality specification not only by the HVL value but also by a second parameter, e.g. the peak tube potential. The differences between the conversion factors f obtained in this work and those deduced from the monoenergetic equivalent energy are evaluated.

Water calorimetry and ionization chamber dosimetry in an 85-MeV clinical proton beam

In recent years, the increased use of proton beams for clinical purposes has enhanced the demand for accurate absolute dosimetry for protons. As calorimetry is the most direct way to establish the absorbed dose and because water has recently been accepted as standard material for this type of beam, the importance of water calorimetry is obvious. In this work we report water calorimeter operation in an 85-MeV proton beam and a comparison of the absorbed dose to water measured by ionometry with the dose resulting from water calorimetric measurements. To ensure a proper understanding of the heat defect for defined impurities in water for this type of radiation, a relative response study was first done in comparison with theoretical calculations of the heat defect. The results showed that pure hypoxic water and hydrogen-saturated water yielded the same response with practically zero heat defect, in agreement with the model calculations. The absorbed dose inferred from these measurements was then compared with the dose derived from ionometry by applying the European Charged Heavy Particle Dosimetry (ECHED) protocol. Restricting the comparison to chambers recommended in the protocol, the calorimeter dose was found to be 2.6% +/- 0.9% lower than the average ionometry dose. In order to estimate the significance of chamber-dependent effects in this deviation, measurements were performed using a set of ten ionization chambers of five different types. The maximum internal deviation in the ionometry results amounted to 1.1%. We detected no systematic chamber volume dependence, but observed a small but systematic effect of the chamber wall thickness. The observed deviation between calorimetry and ionometry can be attributed to a combination of the value of (Wair/e)p for protons, adopted in the ECHED protocol, the mass stopping power ratios of water to air for protons, and possibly small ionization chamber wall effects.

Calibration of low activity 192Ir brachytherapy sources in terms of reference air kerma rate with large volume spherical ionization chambers

An air kerma rate calibration method for low dose rate 192Ir brachytherapy sources was elaborated using three different large volume spherical ionization chambers (PTW LS-10, NE 2551, Exradin A6). To this end these chambers were calibrated for X-ray qualities in the energy range 35 keV to 305 keV and for 137Cs and 60Co gamma rays. The results of these measurements are used to derive mean weighted calibration factors for the photon radiation of 192Ir brachytherapy sources. Furthermore, the effect of the finite chamber size on the effective measuring point and the correction for the contribution of scattered radiation where studied. Measurements on the same sets of sources show that the ionization chambers and the methods used yield results for the reference air kerma rates of 192Ir sources which agree within 0.42%.

Correction factors for a cylindrical ionization chamber used in medium-energy X-ray beams

Correction factors are derived for a cylindrical NE2571 ionization chamber, for absorbed dose determinations in medium-energy X-rays. These new correction factors are proposed as weighted mean values of the factors derived from two methods. The first method is based on a comparative study of dosimetry of medium-energy X-rays with the cylindrical ionization chamber and a water calorimeter. The calorimetric results show that when the ionisation chamber, calibrated free in air, is used to measure the absorbed dose in a water phantom, solely applying mass energy absorption coefficient ratios for the conversion, the ratio of the absorbed dose calorimetry to ionometry is 1.007+or-0.015 at 250 kV (HVL: 2.5 mm Cu), 1.026+or-0.018 at 150 kV (HVL: 0.82 mm Cu) and 1.040+or-0.020 at 100 kV (HVL: 4.54 mm AI). The second method is based on an experimental investigation of the NE2571 ionization chamber in connection with Monte Carlo calculations. Using this method three components in the overall correction factor are investigated: the effect of the displaced volume, the combined angular-energy dependence of the chamber response within the water phantom and the stem correction.

Dependence of overall correction factor of a cylindrical ionization chamber on field size and depth in medium-energy x-ray beams

In this paper we examine the depth and field size dependence of the overall correction factor kch for in-phantom dose determinations in orthovoltage x-ray beams. The overall correction factor is considered to be composed of three contributions, i.e., (1) a contribution from the angular dependence of the chamber response free-in-air, derived based on the measured directional response of the NE2571 for different energies combined with Monte Carlo calculations; (2) a displacement effect and (3) a stem effect, both calculated using the Monte Carlo method for different field sizes and depths. The results show a variation of, at most, 2.2% at the lowest photon energies (29.8-keV average photon energy) when going from 2 cm to 5 cm for a small circular 20-cm2 field. In the medium-energy range (> or = 100 kV), variations are limited to, at most, 1.5% for 120 kV-150 kV when comparing the most extreme variations in field size and depth (i.e., 2-cm depth; 20-cm2 area compared to 5 cm depth; 200-cm2 area). Depth variations most importantly affect the overall correction factor by hardening of the photon fluence spectrum, whereas field diameter variations affect the factor by increase or decrease of contributions of photon scattering. The work shows that taking into account the uncertainties adopted in the recent review of data and methods recommended in the IAEA code of practice, the dependence of the overall correction factor on depth and field size is insignificant for the radiation qualities between 100 kV (HVL 0.17-mm Cu, average energy: 52 keV) and 280 kV (HVL 3.41-mm Cu, average energy: 144 keV).

Kinetics of [methyl-11C]thymidine in patients with squamous cell carcinoma of the head and neck

Thymidine labelled with a position emitter and imaging by positron emission tomography has been suggested to measure tumour proliferation in vivo non-invasively. The present study presents a compartmental model which describes the behavior of [methyl-11C]thymidine in tumour tissue and allows quantification of the data using a proliferation parameter (PP). In nine patients with squamous cell carcinoma of the head and neck, 11C-activity over the tumour and the adjacent blood vessels was measured. The median PP without correction for blood pool activity was found to be 8.61 10(-4) min(-1). With correction this parameter changed from -14% to +40%. We found no correlation between either sex or tumour differentiation and PP. The described model can be used to quantitate [methyl-11C]thymidine uptake in human tumours. Indivdual correction for tumour blood volume is mandatory.