D Shipley

Perturbation correction factors for the NACP-02 plane-parallel ionization chamber in water in high-energy electron beams.

Recent dosimetry protocols for clinical high-energy electron beams recommend measurements of absorbed dose-to-water with a plane-parallel or cylindrical ionization chamber. For well-guarded plane-parallel ionization chambers, the ionization chamber perturbation factor in water, p(Q), has a recommended value of unity in all protocols. This assumption was investigated in detail in this study for one of the recommended ionization chambers in the protocols: the Scanditronix NACP-02 plane-parallel ionization chamber. Monte Carlo (MC) simulations of the NACP-02 ionization chamber with the EGSnrc code were validated against backscatter experiments. MC simulations were then used to calculate p(wall), p(cav) and p(Q) perturbation factors and water-to-air Spencer-Attix stopping powers in 4-19 MeV electron beams of a calibration laboratory (NPL), and in 6-22 MeV clinical electron beams from a Varian CL2300 accelerator. Differences between calculated and the currently recommended (Burns et al 1996 Med. Phys. 23 383-8) stopping powers, water-to-air, were found to be limited to 0.9% at depths between the reference depth z(ref) and the depth where the dose has decreased to 50% of the maximum dose, R50. p(wall) was found to exceed unity by 2.3% in the 4 MeV NPL calibration beam at z(ref). For higher energy electron beams p(wall) decreased to a value of about 1%. Combined with a p(cav) about 1% below unity for all energies at z(ref), this was found to cause p(Q) to exceed unity significantly for all energies. In clinical electron beams all three perturbation factors were found to increase with depth. Our findings indicate that the perturbation factors have to be taken into account in calibration procedures and for clinical depth dose measurements with the NACP-02 ionization chamber.

Validation of a Monte Carlo model of a NACP-02 plane-parallel ionization chamber model using electron backscatter experiments

The accuracy of Monte Carlo (MC) simulation results relies on validating the MC models used in the calculations. In this work, a MC model for the NACP-02 plane-parallel ionization chamber was built and validated against megavoltage electron backscatter experiments using materials of water, graphite, aluminium and copper. Electron energies ranged between 6-18 MeV and the chambers air cavity was at the depth of maximum dose, z(max). A chamber model based on manufacturers specifications resulted in systematic discrepancies of several percents between measured and simulated backscatter factors. Tuning of the MC chamber model against backscatter factors to improve agreement increased the chambers front window mass thickness by 35% over the reported value of 104 mg cm(-2) in the IAEAs TRS-398 absorbed dose protocol. The large increase in chamber window mass thickness was verified by measurements on a disassembled NACP-02 chamber. The new backscatter factor results based on the tuned MC NACP-02 chamber model matched the experimental results within 1-2 standard deviations. We conclude therefore that for MC simulations near z(max), tuning of the NACP-02 chamber model against experimental backscatter measurements is an acceptable method for validating the chamber model.

Evaluation of factors to convert absorbed dose calibrations from graphite to water for the NPL high-energy photon calibration service

The National Physical Laboratory (NPL) provides a high-energy photon calibration service using 4-19 MV x-rays and 60Co gamma-radiation for secondary standard dosemeters in terms of absorbed dose to water. The primary standard used for this service is a graphite calorimeter and so absorbed dose calibrations must be converted from graphite to water. The conversion factors currently in use were determined prior to the launch of this service in 1988. Since then, it has been found that the differences in inherent filtration between the NPL LINAC and typical clinical machines are large enough to affect absorbed dose calibrations and, since 1992, calibrations have been performed in heavily filtered qualities. The conversion factors for heavily filtered qualities were determined by interpolation and extrapolation of lightly filtered results as a function of tissue phantom ratio 20,10 (TPR20,10). This paper aims to evaluate these factors for all mega-voltage photon energies provided by the NPL LINAC for both lightly and heavily filtered qualities and for 60Co y-radiation in two ways. The first method involves the use of the photon fluence-scaling theorem. This states that if two blocks of different material are irradiated by the same photon beam, and if all dimensions are scaled in the inverse ratio of the electron densities of the two media, then, assuming that all photon interactions occur by Compton scatter the photon attenuation and scatter factors at corresponding scaled points of measurement in the phantom will be identical. The second method involves making in-phantom measurements of chamber response at a constant target-chamber distance. Monte Carlo techniques are then used to determine the corresponding dose to the medium in order to determine the chamber calibration factor directly. Values of the ratio of absorbed dose calibration factors in water and in graphite determined in these two ways agree with each other to within 0.2% (1sigma uncertainty). The best fit to both sets of results agrees with values determined in previous work to within 0.3% (1sigma uncertainty). It is found that the conversion factor is not sensitive to beam filtration.

Fluence correction factors for graphite calorimetry in a low-energy clinical proton beam: I. Analytical and Monte Carlo simulations

The conversion of absorbed dose-to-graphite in a graphite phantom to absorbed dose-to-water in a water phantom is performed by water to graphite stopping power ratios. If, however, the charged particle fluence is not equal at equivalent depths in graphite and water, a fluence correction factor, kfl, is required as well. This is particularly relevant to the derivation of absorbed dose-to-water, the quantity of interest in radiotherapy, from a measurement of absorbed dose-to-graphite obtained with a graphite calorimeter. In this work, fluence correction factors for the conversion from dose-to-graphite in a graphite phantom to dose-to-water in a water phantom for 60 MeV mono-energetic protons were calculated using an analytical model and five different Monte Carlo codes (Geant4, FLUKA, MCNPX, SHIELD-HIT and McPTRAN.MEDIA). In general the fluence correction factors are found to be close to unity and the analytical and Monte Carlo codes give consistent values when considering the differences in secondary particle transport. When considering only protons the fluence correction factors are unity at the surface and increase with depth by 0.5% to 1.5% depending on the code. When the fluence of all charged particles is considered, the fluence correction factor is about 0.5% lower than unity at shallow depths predominantly due to the contributions from alpha particles and increases to values above unity near the Bragg peak. Fluence correction factors directly derived from the fluence distributions differential in energy at equivalent depths in water and graphite can be described by kfl = 0.9964 + 0.0024·zw-eq with a relative standard uncertainty of 0.2%. Fluence correction factors derived from a ratio of calculated doses at equivalent depths in water and graphite can be described by kfl = 0.9947 + 0.0024·zw-eq with a relative standard uncertainty of 0.3%. These results are of direct relevance to graphite calorimetry in low-energy protons but given that the fluence correction factor is almost solely influenced by non-elastic nuclear interactions the results are also relevant for plastic phantoms that consist of carbon, oxygen and hydrogen atoms as well as for soft tissues.

Analysis of dose perturbation factors of a NACP-02 ionization chamber in clinical electron beams

For well-guarded plane-parallel ionization chambers, international dosimetry protocols recommend a value of unity for electron perturbation factors in water. However, recent data published by various groups have challenged this. Specifically for the NACP-02 chamber, non-unity electron perturbation factors have already been published by Verhaegen et al (2006 Phys. Med. Biol. 51 1221-35) and Buckley and Rogers (2006 Med. Phys. 33 1788-96). Recently it was found that the mass thickness of the front chamber window can be 35% greater than is listed in the IAEAs TRS-398 absorbed dose protocol (Chin et al 2008 Phys. Med. Biol. 53 N119-26). This study therefore recalculated NACP-02 electron perturbation correction factors for energies 4-18 MeV at depths z(ref) and R(50) to determine the effect of the chamber model change. Results showed that perturbation factors at z(ref) are fairly stable for similar chamber models but become highly sensitive to small changes at deeper depths. The results also showed some dependence on using 1 keV versus 10 keV for the transport cut-off. Additional investigations revealed that the wall perturbation factor, p(wall), is strongly influenced by the chamber back wall at z(ref) and at larger depths small changes in the positioning of the effective point of measurement cause large fluctuations in the final value. Finally, the cavity perturbation factor, p(cav), was found to be primarily influenced by electron backscatter.

Water equivalence of some plastic-water phantom materials for clinical proton beam dosimetry.

Plastic-water phantom materials are not exactly water equivalent since they have a different elemental composition and different interaction cross sections for protons than water. Several studies of the water equivalence of plastic-water phantom materials have been reported for photon and electron beams, but none for clinical proton beams. In proton beams, the difference between non-elastic nuclear interactions in plastic-water phantom materials compared to those in water should be considered. In this work, the water equivalence of Plastic Water® (PW)1, Plastic Water® Diagnostic Therapy (PWDT)1 and solid water (WT1)2 phantoms was studied for clinical proton energies of 60 MeV and 200 MeV. This was done by evaluating the fluence correction factor at equivalent depths; first with respect to water and then with respect to graphite by experiment and Monte Carlo (MC) simulations using FLUKA. MC simulations showed that the fluence correction with respect to water was less than 0.5% up to the entire penetration depth of the protons at 60 MeV and less than 1% at 200 MeV up to 20 cm depth for PWDT, PW and WT1. With respect to graphite the fluence correction was about 0.5% for 60 MeV and about 4% for 200 MeV. The experimental results for modulated and un-modulated 60 MeV proton beams showed good agreement with the MC simulated fluence correction factors with respect to graphite deviating less than 1% from unity for the three plastic-water phantoms.

Conversion from dose-to-graphite to dose-to-water in an 80 MeV/A carbon ion beam.

Based on experiments and numerical simulations, a study is carried out pertaining to the conversion of dose-to-graphite to dose-to-water in a carbon ion beam. This conversion is needed to establish graphite calorimeters as primary standards of absorbed dose in these beams. It is governed by the water-to-graphite mass collision stopping power ratio and fluence correction factors, which depend on the particle fluence distributions in each of the two media. The paper focuses on the experimental and numerical determination of this fluence correction factor for an 80 MeV/A carbon ion beam. Measurements have been performed in the nuclear physics laboratory INFN-LNS in Catania (Sicily, Italy). The numerical simulations have been made with a Geant4 Monte Carlo code through the GATE simulation platform. The experimental data are in good agreement with the simulated results for the fluence correction factors and are found to be close to unity. The experimental values increase with depth reaching 1.010 before the Bragg peak region. They have been determined with an uncertainty of 0.25%. Different numerical results are obtained depending on the level of approximation made in calculating the fluence correction factors. When considering carbon ions only, the difference between measured and calculated values is maximal just before the Bragg peak, but its value is less than 1.005. The numerical value is close to unity at the surface and increases to 1.005 near the Bragg peak. When the fluence of all charged particles is considered, the fluence correction factors are lower than unity at the surface and increase with depth up to 1.025 before the Bragg peak. Besides carbon ions, secondary particles created due to nuclear interactions have to be included in the analysis: boron ions ((10)B and (11)B), beryllium ions ((7)Be), alpha particles and protons. At the conclusion of this work, we have the conversion of dose-to-graphite to dose-to-water to apply to the response of a graphite calorimeter in an 80 MeV/A carbon ion beam. This conversion consists of the product of two contributions: the water-to-graphite electronic mass collision stopping power ratio, which is equal to 1.115, and the fluence correction factor which varies linearly with depth, as k(fl, all) = 0.9995 + 0.0048(zw-eq). The latter has been determined on the basis of experiments and numerical simulations.

Transfer of the UK absorbed dose primary standard for photon beams from the research linac to the clinical linac at NPL

An Elekta Synergy clinical linac facility is now in routine use at the National Physical Laboratory (NPL). For the purpose of therapy-level dosimetry, this has replaced the NPL research linac, which is over 40 years old, and in which the NPL absorbed dose primary standard for high-energy photons was established. This standard has been disseminated to clinical beams by interpolation of the calibration factor as a function of tissue phantom ratio TPR20/10. In this work the absorbed dose standard has been commissioned in all the beams produced by the Elekta Synergy linac. Reference standard ionization chambers have been calibrated in terms of absorbed dose to graphite and this calibration has been converted to one in terms of absorbed dose to water. The results have been combined with the calibration in 60Co γ-rays to obtain measured values for the quality-dependent correction, kQ, for these reference standard chambers used in the Elekta beams. The resulting data are consistent with the interpolated kQ to within 0.4%, which is less than the combined standard uncertainty of kQ, 0.56%.

Evaluation of the water-equivalence of plastic materials in low- and high-energy clinical proton beams

The aim of this work was to evaluate the water-equivalence of new trial plastics designed specifically for light-ion beam dosimetry as well as commercially available plastics in clinical proton beams. The water-equivalence of materials was tested by computing a plastic-to-water conversion factor, [Formula: see text]. Trial materials were characterized experimentally in 60 MeV and 226 MeV un-modulated proton beams and the results were compared with Monte Carlo simulations using the FLUKA code. For the high-energy beam, a comparison between the trial plastics and various commercial plastics was also performed using FLUKA and Geant4 Monte Carlo codes. Experimental information was obtained from laterally integrated depth-dose ionization chamber measurements in water, with and without plastic slabs with variable thicknesses in front of the water phantom. Fluence correction factors, [Formula: see text], between water and various materials were also derived using the Monte Carlo method. For the 60 MeV proton beam, [Formula: see text] and [Formula: see text] factors were within 1% from unity for all trial plastics. For the 226 MeV proton beam, experimental [Formula: see text] values deviated from unity by a maximum of about 1% for the three trial plastics and experimental results showed no advantage regarding which of the plastics was the most equivalent to water. Different magnitudes of corrections were found between Geant4 and FLUKA for the various materials due mainly to the use of different nonelastic nuclear data. Nevertheless, for the 226 MeV proton beam, [Formula: see text] correction factors were within 2% from unity for all the materials. Considering the results from the two Monte Carlo codes, PMMA and trial plastic #3 had the smallest [Formula: see text] values, where maximum deviations from unity were 1%, however, PMMA range differed by 16% from that of water. Overall, [Formula: see text] factors were deviating more from unity than [Formula: see text] factors and could amount to a few percent for some materials.

Development of a primary standard for absorbed dose from unsealed radionuclide solutions

Currently, the determination of the internal absorbed dose to tissue from an administered radionuclide solution relies on Monte Carlo (MC) calculations based on published nuclear decay data, such as emission probabilities and energies. In order to validate these methods with measurements, it is necessary to achieve the required traceability of the internal absorbed dose measurements of a radionuclide solution to a primary standard of absorbed dose. The purpose of this work was to develop a suitable primary standard. A comparison between measurements and calculations of absorbed dose allows the validation of the internal radiation dose assessment methods.The absorbed dose from an yttrium-90 chloride (90YCl) solution was measured with an extrapolation chamber. A phantom was developed at the National Physical Laboratory (NPL), the UKs National Measurement Institute, to position the extrapolation chamber as closely as possible to the surface of the solution. The performance of the extrapolation chamber was characterised and a full uncertainty budget for the absorbed dose determination was obtained. Absorbed dose to air in the collecting volume of the chamber was converted to absorbed dose at the centre of the radionuclide solution by applying a MC calculated correction factor. This allowed a direct comparison of the analytically calculated and experimentally determined absorbed dose of an 90YCl solution.The relative standard uncertainty in the measurement of absorbed dose at the centre of an 90YCl solution with the extrapolation chamber was found to be 1.6% (k = 1). The calculated 90Y absorbed doses from published medical internal radiation dose (MIRD) and radiation dose assessment resource (RADAR) data agreed with measurements to within 1.5% and 1.4%, respectively. This study has shown that it is feasible to use an extrapolation chamber for performing primary standard absorbed dose measurements of an unsealed radionuclide solution. Internal radiation dose assessment methods based on MIRD and RADAR data for 90Y have been validated with experimental absorbed dose determination and they agree within the stated expanded uncertainty (k = 2).