Determination of electron beam mean incident energy from d50 (ionization) values

Depth-ionization measurements have been obtained with an air-filled Nordic Association of Clinical Physicists (NACP) design parallel-plate ionization chamber in a water phantom for ten foil-scattered electron beams from two different machines with nominal energies between 6 and 20 MeV and field sizes from 6 X 6 to 25 X 25 cm2. Depths of 50% ionization and practical range have been determined from least-squares fits to both the raw data and values corrected to parallel-beam geometry using measured virtual source distances. Depths of 50% dose have also been obtained from fits to divergence-corrected depth-dose measurements performed under identical conditions using, a p-type silicon diode detector. Utilizing accepted conversion factors between mean incident energy (E0) and depth of 50% dose for parallel incident beams, and taking advantage of the fact that p-type silicon diode detector readings are nearly directly indicative of relative dose, conversion factors between E0 and depth of 50% ionization for divergence-corrected and raw, uncorrected finite source-surface distance depth-ionization data are empirically determined. Those values, obtained using the results of both ETRAN and EGS4 dose calculations as base lines, are compared to values currently recommended for use in clinical dosimetry.     

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.     

Monte Carlo study of si diode response in electron beams

Silicon semiconductor diodes measure almost the same depth-dose distributions in both photon and electron beams as those measured by ion chambers. A recent study in ion chamber dosimetry has suggested that the wall correction factor for a parallel-plate ion chamber in electron beams changes with depth by as much as 6%. To investigate diode detector response with respect to depth, a silicon diode model is constructed and the water/silicon dose ratio at various depths in electron beams is calculated using EGSnrc. The results indicate that, for this particular diode model, the diode response per unit water dose (or water/diode dose ratio) in both 6 and 18 MeV electron beams is flat within 2% versus depth, from near the phantom surface to the depth of R50 (with calculation uncertainty <0.3%). This suggests that there must be some other correction factors for ion chambers that counter-balance the large wall correction factor at depth in electron beams. In addition, the beam quality and field-size dependence of the diode model are also calculated. The results show that the water/diode dose ratio remains constant within 2% over the electron energy range from 6 to 18 MeV. The water/diode dose ratio does not depend on field size as long as the incident electron beam is broad and the electron energy is high. However, for a very small beam size (1 X 1 cm(2)) and low electron energy (6 MeV), the water/diode dose ratio may decrease by more than 2% compared to that of a broad beam.     

Calculation of photon energy deposition kernels and electron dose point kernels in water

Effects of changes in the physics of EGSnrc compared to EGS4/PRESTA on energy deposition kernels for monoenergetic photons and on dose point kernels for beta sources in water are investigated. In the diagnostic energy range, Compton binding corrections were found to increase the primary energy fraction up to 4.5% at 30 keV with a corresponding reduction of the scatter component of the kernels. Rayleigh scattered photons significantly increase the scatter component of the kernels and reduce the primary energy fraction with a maximum 12% reduction also at 30 keV where the Rayleigh cross section in water has its maximum value. Sampling the photo-electron angular distribution produces a redistribution of the energy deposited by primaries around the interaction site causing differences of up to 2.7 times in the backscattered energy fraction at 20 keV. Above the pair production threshold, the dose distribution versus angle of the primary dose component is significantly different from the EGS4 results. This is related to the more accurate angular sampling of the electron-positron pair direction in EGSnrc as opposed to using a fixed angle approximation in default EGS4. Total energy fractions for photon beams obtained with EGSnrc and EGS4 are almost the same within 0.2%. This fact suggests that the estimate of the total dose at a given point inside an infinite homogeneous water phantom irradiated by broad beams of photons will be very similar for kernels calculated with both codes. However, at interfaces or near boundaries results can be very different especially in the diagnostic energy range. EGSnrc calculated kernels for monoenergetic electrons (50 keV, 100 keV, and 1 MeV) and beta spectra (32P and 90Y) are in excellent agreement with reported EGS4 values except at 1 MeV where inclusion of spin effects in EGSnrc produces an increase of the effective range of electrons. Comparison at 1 MeV with an ETRAN calculation of the electron dose point kernel shows excellent agreement.     

Sensitivity of megavoltage photon beam Monte Carlo simulations to electron beam and other parameters

The BEAM code is used to simulate nine photon beams from three major manufacturers of medical linear accelerators (Varian, Elekta, and Siemens), to derive and evaluate estimates for the parameters of the electron beam incident on the target, and to study the effects of some mechanical parameters like target width, primary collimator opening, flattening filter material and density. The mean energy and the FWHM of the incident electron beam intensity distributions (assumed Gaussian and cylindrically symmetric) are derived by matching calculated percentage depth-dose curves past the depth of maximum dose (within 1% of maximum dose) and off-axis factors (within 2sigma at 1% statistics or less) with measured data from the AAPM RTC TG-46 compilation. The off-axis factors are found to be very sensitive to the mean energy of the electron beam, the FWHM of its intensity distribution, its angle of incidence, the dimensions of the upper opening of the primary collimator, the material of the flattening filter and its density. The off-axis factors are relatively insensitive to the FWHM of the electron beam energy distribution, its divergence and the lateral dimensions of the target. The depth-dose curves are sensitive to the electron beam energy, and to its energy distribution, but they show no sensitivity to the FWHM of the electron beam intensity distribution. The electron beam incident energy can be estimated within 0.2 MeV when matching either the measured off-axis factors or the central-axis depth-dose curves when the calculated uncertainties are about 0.7% at the 1 sigma level. The derived FWHM (+/-0.1 mm) of the electron beam intensity distributions all fall within 1 mm of the manufacturer specifications except in one case where the difference is 1.2 mm.     

Energy spectra, angular spread, fluence profiles and dose distributions of 6 and 18 MV photon beams: results of monte carlo simulations for a varian 2100EX accelerator

The purpose of this study is to provide detailed characteristics of incident photon beams for different field sizes and beam energies. This information is critical to the future development of accurate treatment planning systems. It also enhances our knowledge of radiotherapy photon beams. The EGS4 Monte Carlo code, BEAM, has been used to simulate 6 and 18 MV photon beams from a Varian Clinac-2100EX accelerator. A simulated realistic beam is stored in a phase space data file, which contains details of each particle's complete history including where it has been and where it has interacted. The phase space files are analysed to obtain energy spectra, angular distribution, fluence profile and mean energy profiles at the phantom surface for particles separated according to their charge and history. The accuracy of a simulated beam is validated by the excellent agreement between the Monte Carlo calculated and measured dose distributions. Measured depth-dose curves are obtained from depth-ionization curves by accounting for newly introduced chamber fluence corrections and the stopping-power ratios for realistic beams. The study presents calculated depth-dose components from different particles as well as calculated surface dose and contribution from different particles to surface dose across the field. It is shown that the increase of surface dose with the increase of the field size is mainly due to the increase of incident contaminant charged particles. At 6 MV, the incident charged particles contribute 7% to 21% of maximum dose at the surface when the field size increases from 10 x 10 to 40 x 40 cm2. At 18 MV, their contributions are up to 11% and 29% of maximum dose at the surface for 10 x 10 cm2 and 40 x 40 cm2 fields respectively. However, the fluence of these incident charged particles is less than 1% of incident photon fluence in all cases.     

Beam hardening of 10 MV radiotherapy x-rays: analysis using a convolution/superposition method

Total and primary polyenergetic dose spread arrays (PDSA) have been generated for a high energy 10 MV radiotherapy photon beam using the electron gamma shower (EGS) Monte Carlo code. By considering the attenuation of fluence per energy interval, PDSA have been produced at radiological depths of 0 cm (the surface PDSA) and 40 cm (the beam hardened PDSA). By comparing primary PDSA produced at these different depths, the effect of beam hardening on the PDSA has been quantified. Calculations show that the mean electron range due to the surface primary PDSA is 6.67 mm and the mean electron range of the beam hardened primary PDSA is 8.24 mm. In comparison, a 3 MeV primary monoenergetic dose spread array (MDSA) has a much smaller mean electron range of 4.81 mm. A radiotherapy x-ray beam computation method is introduced which involves a single superposition of the surface generated PDSA or beam hardened PDSA with a polyenergetic TERMA. The mean percentage difference between depth-dose curves obtained using super-position of surface and beam hardened PDSA is only 0.1%. The mean percentage difference from experimental data for these superposition curves is 2.8% down to 40 cm in a homogeneous phantom. The superposition process is shown to be forgiving to spectral differences when calculating the PDSA, but sensitive to the incident photon energy spectrum used to calculate the TERMA.     

Relative electron beam measurements: scaling depths in clear polystyrene to equivalent depths in water

For purposes of mean incident energy determination and the accrual of consistent treatment planning data, measurements of relative ionization or dose made in clear polystyrene must be scaled in depth to produce depth-ionization or depth-dose curves equivalent to what would have been measured in water. Recommendations from various protocols for clear polystyrene to water depth scaling factors differ by as much as 5%. Here, central axis measurements of relative ionization as a function of depth have been made with parallel-plate chambers both in a popular clear polystyrene phantom of density 1.045 g/cm3 and in water. Perturbation and displacement corrections were thus minimized. Both depth-ionization and depth-dose curves were formed for electron beams with nominal incident energies between 6 and 20 MeV and field sizes from 6 cm X 6 cm to 15 cm X 15 cm. Comparisons of the depths for 50% relative reading and practical range between the two phantoms yield average empirical scaling factors of 0.990 and 1.002, respectively.     

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.     

Monte Carlo simulation of small electron fields collimated by the integrated photon MLC

In this study, a Monte Carlo (MC)-based beam model for an ELEKTA linear accelerator was established. The beam model is based on the EGSnrc Monte Carlo code, whereby electron beams with nominal energies of 10, 12 and 15?MeV were considered. For collimation of the electron beam, only the integrated photon multi-leaf-collimators (MLCs) were used. No additional secondary or tertiary add-ons like applicators, cutouts or dedicated electron MLCs were included. The source parameters of the initial electron beam were derived semi-automatically from measurements of depth-dose curves and lateral profiles in a water phantom. A routine to determine the initial electron energy spectra was developed which fits a Gaussian spectrum to the most prominent features of depth-dose curves. The comparisons of calculated and measured depth-dose curves demonstrated agreement within 1%/1?mm. The source divergence angle of initial electrons was fitted to lateral dose profiles beyond the range of electrons, where the imparted dose is mainly due to bremsstrahlung produced in the scattering foils. For accurate modelling of narrow beam segments, the influence of air density on dose calculation was studied. The air density for simulations was adjusted to local values (433?m above sea level) and compared with the standard air supplied by the ICRU data set. The results indicate that the air density is an influential parameter for dose calculations. Furthermore, the default value of the BEAMnrc parameter 'skin depth' for the boundary crossing algorithm was found to be inadequate for the modelling of small electron fields. A higher value for this parameter eliminated discrepancies in too broad dose profiles and an increased dose along the central axis. The beam model was validated with measurements, whereby an agreement mostly within 3%/3?mm was found.     

Electron fluence correction factors for various materials in clinical electron beams

Relative to solid water, electron fluence correction factors at the depth of dose maximum in bone, lung, aluminum, and copper for nominal electron beam energies of 9 MeV and 15 MeV of the Clinac 18 accelerator have been determined experimentally and by Monte Carlo calculation. Thermoluminescent dosimeters were used to measure depth doses in these materials. The measured relative dose at dmax in the various materials versus that of solid water, when irradiated with the same number of monitor units, has been used to calculate the ratio of electron fluence for the various materials to that of solid water. The beams of the Clinac 18 were fully characterized using the EGS4/BEAM system. EGSnrc with the relativistic spin option turned on was used to optimize the primary electron energy at the exit window, and to calculate depth doses in the five phantom materials using the optimized phase-space data. Normalizing all depth doses to the dose maximum in solid water stopping power ratio corrected, measured depth doses and calculated depth doses differ by less than +/- 1% at the depth of dose maximum and by less than 4% elsewhere. Monte Carlo calculated ratios of doses in each material to dose in LiF were used to convert the TLD measurements at the dose maximum into dose at the center of the TLD in the phantom material. Fluence perturbation correction factors for a LiF TLD at the depth of dose maximum deduced from these calculations amount to less than 1% for 0.15 mm thick TLDs in low Z materials and are between 1% and 3% for TLDs in Al and Cu phantoms. Electron fluence ratios of the studied materials relative to solid water vary between 0.83+/-0.01 and 1.55+/-0.02 for materials varying in density from 0.27 g/cm3 (lung) to 8.96 g/cm3 (Cu). The difference in electron fluence ratios derived from measurements and calculations ranges from -1.6% to +0.2% at 9 MeV and from -1.9% to +0.2% at 15 MeV and is not significant at the 1sigma level. Excluding the data for Cu, electron fluence correction factors for open electron beams are approximately proportional to the electron density of the phantom material and only weakly dependent on electron beam energy.     

Commissioning of a medical accelerator photon beam Monte Carlo simulation using wide-field profiles

A method for commissioning an EGSnrc Monte Carlo simulation of medical linac photon beams through wide-field lateral profiles at moderate depth in a water phantom is presented. Although depth-dose profiles are commonly used for nominal energy determination, our study shows that they are quite insensitive to energy changes below 0.3 MeV (0.6 MeV) for a 6 MV (15 MV) photon beam. Also, the depth-dose profile dependence on beam radius adds an additional uncertainty in their use for tuning nominal energy. Simulated 40 cm x 40 cm lateral profiles at 5 cm depth in a water phantom show greater sensitivity to both nominal energy and radius. Beam parameters could be determined by comparing only these curves with measured data.     

Influence of initial electron beam parameters on Monte Carlo calculated absorbed dose distributions for radiotherapy photon beams

Our aim in the present study was to investigate the effects of initial electron beam characteristics on Monte Carlo calculated absorbed dose distribution for a linac 6 MV photon beam. Moreover, the range of values of these parameters was derived, so that the resulted differences between measured and calculated doses were less than 1%. Mean energy, radial intensity distribution and energy spread of the initial electron beam, were studied. The method is based on absorbed dose comparisons of measured and calculated depth-dose and dose-profile curves. All comparisons were performed at 10.0 cm depth, in the umbral region for dose-profile and for depths past maximum for depth-dose curves. Depth-dose and dose-profile curves were considerably affected by the mean energy of electron beam, with dose profiles to be more sensitive on that parameter. The depth-dose curves were unaffected by the radial intensity of electron beam. In contrast, dose-profile curves were affected by the radial intensity of initial electron beam for a large field size. No influence was observed in dose-profile or depth-dose curves with respect to energy spread variations of electron beam. Conclusively, simulating the radiation source of a photon beam, two of the examined parameters (mean energy and radial intensity) of the electron beam should be tuned accurately, so that the resulting absorbed doses are within acceptable precision. The suggested method of evaluating these crucial but often poorly specified parameters may be of value in the Monte Carlo simulation of linear accelerator photon beams.     

Depth dose enhancement of electron beams subject to external uniform longitudinal magnetic fields: a Monte Carlo study

We studied the dose distributions from electron beams subjected to a longitudinal magnetic field applied to them before they reach the phantom. We found that specific combinations of the length and intensity of the magnetic field produced enhancement of the peaks of the central-axis depth-dose distributions. The EGS4 Monte Carlo system was used in this study. In the simulations, a uniform axial magnetic field parallel to the electron beam direction was applied to the air gap between the collimation and the phantom. We extensively studied the simplified case of an 18 MeV electron beam point source. Dose deposition was calculated for various magnetic field strengths, distances through which the magnetic field was applied, collimation sizes, and source to collimation distances. The magnetic field strengths varied from 0 to 3 T, the source-to-collimation distances studied were 50 and 95 cm, the collimation sizes studied were 10 x 10 and 20 x 20 cm2, and the distance through which the field was applied ranged from 10 to 20 cm. Specific combinations of these variables resulted in as much as a 70% enhancement of the peak dose relative to the surface dose. Finally, to determine how the geometry of a real accelerator affects the resulting dose distribution, we performed a full simulation of an Elekta SL20 linear accelerator and compared the results with the ideal case.     

Dose properties of a laser accelerated electron beam and prospects for clinical application

Laser wakefield acceleration (LWFA) technology has evolved to where it should be evaluated for its potential as a future competitor to existing technology that produces electron and x-ray beams. The purpose of the present work is to investigate the dosimetric properties of an electron beam that should be achievable using existing LWFA technology, and to document the necessary improvements to make radiotherapy application for LWFA viable. This paper first qualitatively reviews the fundamental principles of LWFA and describes a potential design for a 30 cm accelerator chamber containing a gas target. Electron beam energy spectra, upon which our dose calculations are based, were obtained from a uniform energy distribution and from two-dimensional particle-in-cell (2D PIC) simulations. The 2D PIC simulation parameters are consistent with those reported by a previous LWFA experiment. According to the 2D PIC simulations, only approximately 0.3% of the LWFA electrons are emitted with an energy greater than 1 MeV. We studied only the high-energy electrons to determine their potential for clinical electron beams of central energy from 9 to 21 MeV. Each electron beam was broadened and flattened by designing a dual scattering foil system to produce a uniform beam (103%>off-axis ratio>95%) over a 25 x 25 cm2 field. An energy window (deltaE) ranging from 0.5 to 6.5 MeV was selected to study central-axis depth dose, beam flatness, and dose rate. Dose was calculated in water at a 100 cm source-to-surface distance using the EGS/BEAM Monte Carlo algorithm. Calculations showed that the beam flatness was fairly insensitive to deltaE. However, since the falloff of the depth-dose curve (R10-R90) and the dose rate both increase with deltaE, a tradeoff between minimizing (R10-R90) and maximizing dose rate is implied. If deltaE is constrained so that R10-R90 is within 0.5 cm of its value for a monoenergetic beam, the maximum practical dose rate based on 2D PIC is approximately 0.1 Gy min(-1) for a 9 MeV beam and 0.03 Gy min(-1) for a 15 MeV beam. It was concluded that current LWFA technology should allow a table-top terawatt (T3) laser to produce therapeutic electron beams that have acceptable flatness, penetration, and falloff of depth dose; however, the dose rate is still 1%-3% of that which would be acceptable, especially for higher-energy electron beams. Further progress in laser technology, e.g., increasing the pulse repetition rate or number of high energy electrons generated per pulse, is necessary to give dose rates acceptable for electron beams. Future measurements confirming dosimetric calculations are required to substantiate our results. In addition to achieving adequate dose rate, significant engineering developments are needed for this technology to compete with current electron acceleration technology. Also, the functional benefits of LWFA electron beams require further study and evaluation.     

A diamond detector in the dosimetry of high-energy electron and photon beams.

A diamond detector type 60003 (PTW Freiburg) was examined for the purpose of dosimetry with 4-20 MeV electron beams and 4-25 MV photon beams. Results were compared with those obtained by using a Markus chamber for electron beams and an ionization chamber for photon beams. Dose distributions were measured in a water phantom with the detector connected to a Unidos electrometer (PTW Freiburg). After a pre-irradiation of about 5 Gy the diamond detector shows a stability in response which is better than that of an ionization chamber. The current of the diamond detector was measured under variation of photon beam dose rate between 0.1 and 7 Gy min(-1). Different FSDs were chosen. Furthermore the pulse repetition frequency and the depth of the detector were changed. The electron beam dose rate was varied between 0.23 and 4.6 Gy min(-1) by changing the pulse-repetition frequency. The response shows no energy dependence within the covered photon-beam energy range. Between 4 MeV and 18 MeV electron beam energy it shows only a small energy dependence of about 2%, as expected from theory. For smaller electron energies the response increases significantly and an influence of the contact material used for the diamond detector can be surmised. A slight sublinearity of the current and dose rate was found. Detector current and dose rate are related by the expression i alpha Ddelta, where i is the detector current, D is the dose rate and delta is a correction factor of approximately 0.963. Depth-dose curves of photon beams, measured with the diamond detector, show a slight overestimation compared with measurements with the ionization chamber. This overestimation is compensated for by the above correction term. The superior spatial resolution of the diamond detector leads to minor deviations between depth-dose curves of electron beams measured with a Markus chamber and a diamond detector.     

Are neutrons responsible for the dose discrepancies between Monte Carlo calculations and measurements in the build-up region for a high-energy photon beam?

This study presents measured neutron dose using a neutron dosimeter in a water phantom and investigates a hypothesis that neutrons in a high-energy photon beam may be responsible for the reported significant dose discrepancies between Monte Carlo calculations and measurements at the build-up region in large fields. Borated polyethylene slabs were inserted between the accelerator head and the phantom in order to remove neutrons generated in the accelerator head. The thickness of the slab ranged from 2.5 cm to 10 cm. A lead slab of 3 mm thickness was also used in the study. The superheated drop neutron dosimeter was used to measure the depth-dose curve of neutrons in a high-energy photon beam and to verify the effectiveness of the slab to remove these neutrons. Total dose measurements were performed in water using a WELLHOFER WP700 beam scanner with an IC-10 ionization chamber. The Monte Carlo code BEAM was used to simulate an 18 MV photon beam from a Varian Clinac-2100EX accelerator. Both EGS4/DOSXYZ and EGSnrc/DOSRZnrc were used in the dose calculations. Measured neutron dose equivalents as a function of depth per unit total dose in water were presented for 10 x 10 and 40 x 40 cm2 fields. The measured results have shown that a 5-10 cm thick borated polyethylene slab can reduce the neutron dose by a factor of 2 when inserted between the accelerator head and the detector. In all cases the measured neutron dose equivalent was less than 0.5% of the photon dose. In order to study if the ion chamber was highly sensitive to the neutron dose, we have investigated the disagreement between the Monte Carlo calculated and measured central-axis depth-dose curves in the build-up region when different shielding materials were used. The result indicated that the IC-10 chamber was not highly sensitive to the neutron dose. Therefore, neutrons present in a high-energy photon beam were unlikely to be responsible for the reported discrepancies in the build-up region for large fields.     

The antiproton depth-dose curve in water

We have measured the depth-dose curve of 126 MeV antiprotons in a water phantom using ionization chambers. Since the antiproton beam provided by CERN has a pulsed structure and possibly carries a high-LET component from the antiproton annihilation, it is necessary to correct the acquired charge for ion recombination effects. The results are compared with Monte Carlo calculations and were found to be in good agreement. Based on this agreement we calculate the antiproton depth-dose curve for antiprotons and compare it with that for protons and find a doubling of the physical dose in the peak region for antiprotons.     

Monte Carlo calculations of correction factors for plane-parallel ionization chambers in clinical electron dosimetry

Recent standard dosimetry protocols recommend that plane-parallel ionization chambers be used in the measurements of depth-dose distributions or the calibration of low-energy electron beams with beam quality R50 <4 g/cm2. In electron dosimetry protocols with the plane-parallel chambers, the wall correction factor, Pwall, in water is assumed to be unity and the replacement correction factor, Prepl, is taken to be unity for well-guarded plane-parallel chambers, at all measurement depths. This study calculated Pwall and Prepl for NACP-02, Markus, and Roos plane-parallel chambers in clinical electron dosimetry using the EGSnrc Monte Carlo code system. The Pwall values for the plane-parallel chambers increased rapidly as a function of depth in water, especially at lower energy. The value around R50 for NACP-02 was about 10% greater than unity at 4 MeV. The effect was smaller for higher electron energies. Similarly, Prepl values with depth increased drastically at the region with the steep dose gradient for lower energy. For Markus Prepl departed more than 10% from unity close to R50 due to the narrow guard ring width. Prepl for NACP-02 and Roos was close to unity in the plateau region of depth-dose curves that includes a reference depth, dref. It was also found that the ratio of the dose to water and the dose to the sensitive volume in the air cavity for the plane-parallel chambers, Dw/[Dair]pp, at d(ref) differs significantly from that assumed by electron dosimetry protocols.