Séverine Rossomme

Fluence correction factors and stopping power ratios for clinical ion beams

Abstract Background. In radiation therapy, the principal dosimetric quantity of interest is the absorbed dose to water. Therefore, a dose conversion to dose to water is required for dose deposited by ion beams in other media. This is in particular necessary for dose measurements in plastic phantoms for increased positioning accuracy, graphite calorimetry being developed as a primary standard for dose to water dosimetry, but also for the comparison of dose distributions from Monte Carlo simulations with those of pencil beam algorithms. Material and methods. In the conversion of absorbed dose to phantom material to absorbed dose to water the water-to-material stopping power ratios (STPR) and the fluence correction factors (FCF) for the full charged particle spectra are needed. We determined STPR as well as FCF for water to graphite, bone (compact), and PMMA as a function of water equivalent depth, zw, with the Monte Carlo code SHIELD-HIT10A. Simulations considering all secondary ions were performed for primary protons as well as carbon, nitrogen and oxygen ions with a total range of 3 cm, 14.5 cm and 27 cm as well as for two spread-out Bragg-peaks (SOBP). STPR as a function of depth are also compared to a recently proposed analytical formula. Results. The STPR are of the order of 1.022, 1.070, and 1.112 for PMMA, bone, and graphite, respectively. STPR vary only little with depth except close to the total range of the ion and they can be accurately approximated with an analytical formula. The amplitude of the FCF depends on the non-elastic nuclear interactions and it is unity if these interactions are turned off in the simulation. Fluence corrections are of the order of a percent becoming more pronounced for larger depths resulting in dose difference of the order of 5% around 25 cm. The same order of magnitude is observed for SOBP. Conclusions. We conclude that for ions with small total range (zw-eq ≤3 cm) dosimetry without applying FCF could in principle be performed in phantoms of materials other than water without a significant loss of accuracy. However, in clinical high-energy ion beams with penetration depths zw-eq ≥3 cm, where accurate positioning in water is not an issue, absorbed dose measurements should be directly performed in water or accurate values of FCF need to be established.

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.

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.

Under-response of a PTW-60019 microDiamond detector in the Bragg peak of a 62 MeV/n carbon ion beam.

To investigate the linear energy transfer (LET) dependence of the response of a PTW-60019 Freiburg microDiamond detector, its response was compared to the response of a plane-parallel Markus chamber in a 62 MeV/n mono-energetic carbon ion beam. Results obtained with two different experimental setups are in agreement. As recommended by IAEA TRS-398, the response of the Markus chamber was corrected for temperature, pressure, polarity effects and ion recombination. No correction was applied to the response of the microDiamond detector. The ratio of the response of the Markus chamber to the response of the microDiamond is close to unity in the plateau region. In the Bragg peak region, a significant increase of the ratio is observed, which increases to 1.2 in the distal edge region. Results indicate a correlation between the under-response of the microDiamond detector and high LET values. The combined relative standard uncertainty of the results is estimated to be 2.38% in the plateau region and 12% in the distal edge region. These values are dominated by the uncertainty of alignment in the non-uniform beam and the uncertainty of range determination.

Development and application of a water calorimeter for the absolute dosimetry of short-range particle beams.

In this work, we describe a new design of water calorimeter built to measure absorbed dose in non-standard radiation fields with reference depths in the range of 6-20 mm, and its initial testing in clinical electron and proton beams. A functioning calorimeter prototype with a total water equivalent thickness of less than 30 mm was constructed in-house and used to obtain measurements in clinical accelerator-based 6 MeV and 8 MeV electron beams and cyclotron-based 60 MeV monoenergetic and modulated proton beams. Corrections for the conductive heat transfer due to dose gradients and non-water materials was also accounted for using a commercial finite element method software package. Absorbed dose to water was measured with an associated type A standard uncertainty of approximately 0.4% and 0.2% for the electron and proton beam experiments, respectively. In terms of thermal stability, drifts were on the order of a couple of hundred µK min-1, with a short-term variation of 5-10 µK. Heat transfer correction factors ranged between 1.021 and 1.049. The overall combined standard uncertainty on the absorbed dose to water was estimated to be 0.6% for the 6 MeV and 8 MeV electron beams, as well as for the 60 MeV monoenergetic protons, and 0.7% for the modulated 60 MeV proton beam. This study establishes the feasibility of developing an absorbed dose transfer standard for short-range clinical electrons and protons and forms the basis for a transportable dose standard for direct calibration of ionization chambers in the users beam. The largest contributions to the combined standard uncertainty were the positioning (⩽0.5%) and the correction due to conductive heat transfer (⩽0.4%). This is the first time that water calorimetry has been used in such a low energy proton beam.

LET dependence of the response of a PTW-60019 microDiamond detector in a 62 MeV proton beam

This study was initiated following conclusions from earlier experimental work, performed in a low-energy carbon ion beam, indicating a significant LET dependence of the response of a PTW-60019 microDiamond detector. The purpose of this paper is to present a comparison between the response of the same PTW-60019 microDiamond detector and an IBA Roos-type ionization chamber as a function of depth in a 62MeV proton beam. Even though proton beams are considered as low linear energy transfer (LET) beams, the LET value increases slightly in the Bragg peak region. Contrary to the observations made in the carbon ion beam, in the 62MeV proton beam good agreement is found between both detectors in both the plateau and the distal edge region. No significant LET dependent response of the PTW-60019 microDiamond detector is observed consistent with other findings for proton beams in the literature, despite this particular detector exhibiting a substantial LET dependence in a carbon ion beam.

SU-E-T-408: Determination of KQ,Q0-Factors From Water and Graphite Calorimetry in a 60 MeV Proton Beam

PURPOSE To reduce the uncertainty of the beam quality correction factor kQ,Q0, for scattered proton beams (SPB). This factor is used in dosimetry protocols, to determine absorbed dose-to-water with ionization chambers. For the Roos plane parallel chambers (RPPICs), the IAEA TRS-398 protocol estimates kQ,Q0-factor to be 1.004(for a beam quality Rres=2 g.cm2 ), with an uncertainty of 2.1%. METHODS A graphite calorimeter (GCal), a water calorimeter (WCal) and RPPICs were exposed, in a single experiment, to a 60 MeV non-modulated SPB. RPPICs were calibrated in terms of absorbed dose-to-water in a 20 MeV electron beam. The calibration coefficient is traceable to NPLs absorbed dose standards. Chamber measurements were corrected for environmental conditions, recombination and polarity. The WCal corrections include heat loss, heat defect and vessel perturbation. The GCal corrections include heat loss and absorbed dose conversion. Except for heat loss correction and its uncertainty in the WCal system, all major corrections were included in the analysis. Other minor corrections, such as beam profile non-uniformity, are still to be evaluated. Experimental kQ,Q0-factors were derived by comparing the results obtained with both calorimeters and ionometry. RESULTS The absorbed dose-to-water from both calorimeters was found to be within 1.3% with an uncertainty of 1.2%. kQ,Q0-factor for a RPPIC was found to be 0.998 and 1.011, with a standard uncertainty of 1.4% and 0.9% when the dose is based on the GCal and the WCal, respectively. CONCLUSION Results suggest the possibility to determine kQ,Q0-values for PPICs in SPB with a lower uncertainty than specified in the TRS-398 thereby helping to reduce uncertainty on absorbed dose-to-water. The agreement between calorimeters confirms the possibility to use GCal or WCal as primary standard in SPB. Because of the dose conversion, the use of GCal may lead to slightly higher uncertainty, but is, at present, considerably easier to operate.

Ion recombination correction factor in scanned light-ion beams for absolute dose measurement using plane-parallel ionisation chambers

Based on international reference dosimetry protocols for light-ion beams, a correction factor (k s) has to be applied to the response of a plane-parallel ionisation chamber, to account for recombination of negative and positive charges in its air cavity before these charges can be collected on the electrodes. In this work, k s for IBA PPC40 Roos-type chambers is investigated in four scanned light-ion beams (proton, helium, carbon and oxygen). To take into account the high dose-rates used with scanned beams and LET-values, experimental results are compared to a model combining two theories. One theory, developed by Jaffé, describes the variation of k s with the ionization density within the ion track (initial recombination) and the other theory, developed by Boag, describes the variation of k s with the dose rate (volume recombination). Excellent agreement is found between experimental and theoretical k s-values. All results confirm that k s cannot be neglected. The solution to minimise k s is to use the ionisation chamber at high voltage. However, one must be aware that charge multiplication may complicate the interpretation of the measurement. For the chamber tested, it was found that a voltage of 300 V can be used without further complication. As the initial recombination has a logarithmic variation as a function of 1/V, the two-voltage method is not applicable to these scanned beams.

Fluence correction factor for graphite calorimetry in a clinical high-energy carbon-ion beam

The aim of this work is to develop and adapt a formalism to determine absorbed dose to water from graphite calorimetry measurements in carbon-ion beams. Fluence correction factors, [Formula: see text], needed when using a graphite calorimeter to derive dose to water, were determined in a clinical high-energy carbon-ion beam. Measurements were performed in a 290 MeV/n carbon-ion beam with a field size of 11  ×  11 cm2, without modulation. In order to sample the beam, a plane-parallel Roos ionization chamber was chosen for its small collecting volume in comparison with the field size. Experimental information on fluence corrections was obtained from depth-dose measurements in water. This procedure was repeated with graphite plates in front of the water phantom. Fluence corrections were also obtained with Monte Carlo simulations through the implementation of three methods based on (i) the fluence distributions differential in energy, (ii) a ratio of calculated doses in water and graphite at equivalent depths and (iii) simulations of the experimental setup. The [Formula: see text] term increased in depth from 1.00 at the entrance toward 1.02 at a depth near the Bragg peak, and the average difference between experimental and numerical simulations was about 0.13%. Compared to proton beams, there was no reduction of the [Formula: see text] due to alpha particles because the secondary particle spectrum is dominated by projectile fragmentation. By developing a practical dose conversion technique, this work contributes to improving the determination of absolute dose to water from graphite calorimetry in carbon-ion beams.

Response of synthetic diamond detectors in proton, carbon, and oxygen ion beams

Purpose: In this work, the LET‐dependence of the response of synthetic diamond detectors is investigated in different particle beams. Method: Measurements were performed in three nonmodulated particle beams (proton, carbon, and oxygen). The response of five synthetic diamond detectors was compared to the response of a Markus or an Advanced Markus ionization chamber. The synthetic diamond detectors were used with their axis parallel to the beam axis and without any bias voltage. A high bias voltage was applied to the ionization chambers, to minimize ion recombination, for which no correction is applied (+300 V and +400 V were applied to the Markus and Advanced Markus ionization chambers respectively). Results: The ratio between the normalized response of the synthetic diamond detectors and the normalized response of the ionization chamber shows an under‐response of the synthetic diamond detectors in carbon and oxygen ion beams. No under‐response of the synthetic diamond detectors is observed in protons. For each beam, combining results obtained for the five synthetic diamond detectors and considering the uncertainties, a linear fit of the ratio between the normalized response of the synthetic diamond detectors and the normalized response of the ionization chamber is determined. The response of the synthetic diamond detectors can be described as a function of LET as (−6.22E‐4 ± 3.17E‐3) • LET + (0.99 ± 0.01) in proton beam, (−2.51E‐4 ± 1.18E‐4) • LET + (1.01 ± 0.01) in carbon ion beam and (−2.77E‐4 ± 0.56E‐4) • LET + (1.03 ± 0.01) in oxygen ion beam. Combining results obtained in carbon and oxygen ion beams, a LET dependence of about 0.026% (±0.013%) per keV/μm is estimated. Conclusions: Due to the high LET value, a LET dependence of the response of the synthetic diamond detector was observed in the case of carbon and oxygen beams. The effect was found to be negligible in proton beams, due to the low LET value. The under‐response of the synthetic diamond detector may result from the recombination of electron/hole in the thin synthetic diamond layer, due to the high LET‐values. More investigations are required to confirm this assumption.