Roberto Capote

Ionization chamber dosimetry of small photon fields: a Monte Carlo study on stopping-power ratios for radiosurgery and IMRT beams

Absolute dosimetry with ionization chambers of the narrow photon fields used in stereotactic techniques and IMRT beamlets is constrained by lack of electron equilibrium in the radiation field. It is questionable that stopping-power ratio in dosimetry protocols, obtained for broad photon beams and quasi-electron equilibrium conditions, can be used in the dosimetry of narrow fields while keeping the uncertainty at the same level as for the broad beams used in accelerator calibrations. Monte Carlo simulations have been performed for two 6 MV clinical accelerators (Elekta SL-18 and Siemens Mevatron Primus), equipped with radiosurgery applicators and MLC. Narrow circular and Z-shaped on-axis and off-axis fields, as well as broad IMRT configured beams, have been simulated together with reference 10 x 10 cm2 beams. Phase-space data have been used to generate 3D dose distributions which have been compared satisfactorily with experimental profiles (ion chamber, diodes and film). Photon and electron spectra at various depths in water have been calculated, followed by Spencer-Attix (delta = 10 keV) stopping-power ratio calculations which have been compared to those used in the IAEA TRS-398 code of practice. For water/air and PMMA/air stopping-power ratios, agreements within 0.1% have been obtained for the 10 x 10 cm2 fields. For radiosurgery applicators and narrow MLC beams, the calculated s(w,air) values agree with the reference within +/-0.3%, well within the estimated standard uncertainty of the reference stopping-power ratios (0.5%). Ionization chamber dosimetry of narrow beams at the photon qualities used in this work (6 MV) can therefore be based on stopping-power ratios data in dosimetry protocols. For a modulated 6 MV broad beam used in clinical IMRT, s(w,air) agrees within 0.1% with the value for 10 x 10 cm2, confirming that at low energies IMRT absolute dosimetry can also be based on data for open reference fields. At higher energies (24 MV) the difference in s(w,air) was up to 1.1%, indicating that the use of protocol data for narrow beams in such cases is less accurate than at low energies, and detailed calculations of the dosimetry parameters involved should be performed if similar accuracy to that of 6 MV is sought.

An EGSnrc Monte Carlo study of the microionization chamber for reference dosimetry of narrow irregular IMRT beamlets

Intensity modulated radiation therapy (IMRT) has evolved toward the use of many small radiation fields, or beamlets, to increase the resolution of the intensity map. The size of smaller beamlets can be typically about 1-5 cm2. Therefore small ionization chambers (IC) with sensitive volumes < or = 0.1 cm3 are generally used for dose verification of IMRT treatment. The dosimetry of these narrow photon beams pertains to the so-called nonreference conditions for beam calibration. The use of ion chambers for such narrow beams remains questionable due to the lack of electron equilibrium in most of the field. The present contribution aims to estimate, by the Monte Carlo (MC) method, the total correction needed to convert the IBA-Wellhöfer NAC007 micro IC measured charge in such radiation field to the absolute dose to water. Detailed geometrical simulation of the microionization chamber was performed. The ion chamber was always positioned at a 10 cm depth in water, parallel to the beam axis. The delivered doses to air and water cavity were calculated using the CAVRZ EGSnrc user code. The 6 MV phase-spaces for Primus Clinac (Siemens) used as an input to the CAVRZnrc code were derived by BEAM/EGS4 modeling of the treatment head of the machine along with the multileaf collimator [Sánchez-Doblado et al., Phys. Med. Biol. 48, 2081-2099 (2003)] and contrasted with experimental measurements. Dose calculations were carried out for two irradiation geometries, namely, the reference 10x10 cm2 field and an irregular (approximately 2x2 cm2) IMRT beamlet. The dose measured by the ion chamber is estimated by MC simulation as a dose averaged over the air cavity inside the ion-chamber (Dair). The absorbed dose to water is derived as the dose deposited inside the same volume, in the same geometrical position, filled and surrounded by water (Dwater) in the absence of the ionization chamber. Therefore, the Dwater/Dair dose ratio is a MC direct estimation of the total correction factor needed to convert the absorbed dose in air to absorbed dose to water. The dose ratio was calculated for several chamber positions, starting from the penumbra region around the beamlet along the two diagonals crossing the radiation field. For this quantity from 0 up to a 3% difference is observed between the dose ratio values obtained within the small irregular IMRT beamlet in comparison with the dose ratio derived for the reference 10x10 cm2 field. Greater differences from the reference value up to 9% were obtained in the penumbra region of the small IMRT beamlet.

Total skin electron therapy treatment verification: Monte Carlo simulation and beam characteristics of large non-standard electron fields

Total skin electron therapy (TSET) is a complex technique which requires non-standard measurements and dosimetric procedures. This paper investigates an essential first step towards TSET Monte Carlo (MC) verification. The non-standard 6 MeV 40 x 40 cm2 electron beam at a source to surface distance (SSD) of 100 cm as well as its horizontal projection behind a polymethylmethacrylate (PMMA) screen to SSD = 380 cm were evaluated. The EGS4 OMEGA-BEAM code package running on a Linux home made 47 PCs cluster was used for the MC simulations. Percentage depth-dose curves and profiles were calculated and measured experimentally for the 40 x 40 cm2 field at both SSD = 100 cm and patient surface SSD = 380 cm. The output factor (OF) between the reference 40 x 40 cm2 open field and its horizontal projection as TSET beam at SSD = 380 cm was also measured for comparison with MC results. The accuracy of the simulated beam was validated by the good agreement to within 2% between measured relative dose distributions, including the beam characteristic parameters (R50, R80, R100, Rp, E0) and the MC calculated results. The energy spectrum, fluence and angular distribution at different stages of the beam (at SSD = 100 cm, at SSD = 364.2 cm, behind the PMMA beam spoiler screen and at treatment surface SSD = 380 cm) were derived from MC simulations. Results showed a final decrease in mean energy of almost 56% from the exit window to the treatment surface. A broader angular distribution (FWHM of the angular distribution increased from 13 degrees at SSD = 100 cm to more than 30 degrees at the treatment surface) was fully attributable to the PMMA beam spoiler screen. OF calculations and measurements agreed to less than 1%. The effect of changing the electron energy cut-off from 0.7 MeV to 0.521 MeV and air density fluctuations in the bunker which could affect the MC results were shown to have a negligible impact on the beam fluence distributions. Results proved the applicability of using MC as a treatment verification tool for complex radiotherapy techniques.

Analytical expressions for the dispersive contributions to the nucleon-nucleus optical potential

Mahaux and co-workers @1‐5# have shown how the study of the nuclear mean field may benefit from the use of dispersion relations. These are mathematical expressions that link certain contributions to the real and imaginary components of the optical model potential ~OMP!. The constraint imposed by these dispersion relations helps in reducing ambiguities in the construction of phenomenological potentials from fits to the experimental data. We refer specifically to the so-called dispersive contribution DV, which adds dynamical content to the otherwise static ~and real! Hartree-Fock potential term VHF . Under favorable conditions of analyticity in the complex E plane, the real part DV can be constructed from the knowledge of the imaginary part W of the mean field on the real axis through the dispersion relation

Dispersion relations in the nuclear optical model

We have developed a code to calculate integrals arising from the dispersion relation between the real and the imaginary parts of the nuclear optical model potential (OMP). Both, analytical solution for the most common dispersive OMP, and the general numerical solution, are included. In the numerical integration, fast convergence is achieved by means of the Gauss–Legendre integration method, which offers accuracy, easiness of implementation and generality for dispersive optical model calculations. The numerical method is validated versus analytical solution. The use of this package in the OMP parameter search codes allows for an efficient and accurate dispersive analysis. uf6d9 2003 Elsevier Science B.V. All rights reserved.

A standard dosimetry procedure for 192Ir sources used for endovascular brachytherapy

The experimental dosimetry of a high dose rate (HDR) 192Ir source used for the brachytherapy of peripheral vessels is reported. The direct determination of the reference air kerma rate Kr agrees, within the experimental uncertainty, with the results obtained by a well ionization chamber calibrated at the NIST and the manufacturers certification. A highly sensitive (HS) radiochromic film (RCF), that presents only one active layer, was used for the source dosimetry in a water phantom. The adopted experimental set-up, with the source in its catheter positioned on the RCF plane, seems to have given better accuracy of the RCF optical density measurements. The agreement between the measurement of the dose rate constant DKr (10 mm, pi/2) and the literature data confirmed the coherence of the HS RCF calibration obtained by the kerma in air measurements. The RCF measurements supplied dosimetric information about the dose to water per reference air kerma rate D(r, theta)/Kr along the source transverse bisector axis, the radial dose function g(r) and the anisotropy function F(r, theta). The value D(2 mm, pi/2)/Kr = 22.4 +/- 1.2 cGy h(-1)/(microGy h(-1)) is supplied with a dose uncertainty that is essentially due to the indeterminacy of the source position in the catheter. The data of the radial and anisotropy functions have been compared with Monte Carlo determinations reported in the literature.

A general numerical solution of dispersion relations for the nuclear optical model

A general numerical solution of the dispersion integral relation between the real and the imaginary parts of the nuclear optical potential is presented. Fast convergence is achieved by means of the Gauss-Legendre integration method, which offers accuracy, ease of implementation and generality for dispersive optical model calculations. The use of this numerical integration method in the optical model parameter search codes allows for a fast and accurate dispersive analysis.

Experimental dosimetry of a 32P catheter-based endovascular brachytherapy source.

The experimental dosimetry in a water phantom of a 32P linear source, 20 mm in length, used for the brachytherapy of coronary vessels is reported. The source content activity, A, was determined by means of a calibrated well ion-chamber and the value was compared with the contained activity reported in the manufacturers certification. In this field of brachytherapy dosimetry, radiochromic film supplies a high enough spatial resolution. A highly sensitive radiochromic film, that presents only one active layer, was used in this work for the source dosimetry in a water phantom. The radiochromic film was characterized by electron beams produced by a clinical linac. A Monte Carlo calculation of beta spectra in water at different distances along the source transverse bisector axis allowed to take into account the low dependence of film response from the electron beam energy. The adopted experimental set-up, with the source in its catheter positioned on the film plane inside the water phantom, supplies accurate dosimetric information. The measured dose rate to water per unit of source activity at reference distance, D(r0, theta0)/A, in units of cGy s(-1) GBq(-1), was in agreement with the value reported in the manufacturers certification within the experimental uncertainty. The radial dose function, g(r), is in good agreement with the literature data. The anisotropy function F(r, theta) is also reported. The analysis of the dose profile obtained at 2 mm from the source longitudinal axis shows that the uniformity is within 10% along 75% of the 20 mm treatment length. The adopted experimental set-up seems to be adequate for the quality control procedure of the dose homogeneity distribution in the water medium.

The wall correction factor for a spherical ionization chamber used in brachytherapy source calibration

The effect of wall chamber attenuation and scattering is one of the most important corrections that must be determined when the linear interpolation method between two calibration factors of an ionization chamber is used. For spherical ionization chambers the corresponding correction factors A(w) have to be determined by a non-linear trend of the response as a function of the wall thickness. The Monte Carlo and experimental data here reported show that the A(w) factors obtained for an Exradin A4 chamber, used in the brachytherapy source calibration, in terms of reference air kerma rate, are up to 1.2% greater than the values obtained by the linear extrapolation method for the studied beam qualities. Using the Aw factors derived from Monte Carlo calculations, the accuracy of the calibration factor N(K,Ir) for the Exradin A4, obtained by the interpolation between two calibration factors, improves about 0.6%. The discrepancy between the new calculated factor and that obtained using the complete calibration curve of the ion-chamber and the 192Ir spectrum is only 0.1%.

Dispersive spherical optical model of neutron scattering from 27 Al up to 250 MeV

A spherical optical model potential (OMP) containing a dispersive term is used to fit the available experimental database of angular distribution and total cross section data for n + Al27 covering the energy range 0.1- 250 MeV using relativistic kinematics and a relativistic extension of the Schroedinger equation. A dispersive OMP with parameters that show a smooth energy dependence and energy independent geometry are determined from fits to the entire data set. A very good overall agreement between experimental data and predictions is achieved up to 150 MeV. Inclusion of nonlocality effects in the absorptive volume potential allows to achieve an excellent agreement up to 250 MeV.