A new method for the measurement of bremsstrahlung spectra

An improved physics-based approach for unfolding megavoltage bremsstrahlung spectra using transmission analysis

PURPOSE

To develop a physics-based approach to improve the accuracy and robustness of the ill-conditioned problem of unfolding megavoltage bremsstrahlung spectra from transmission data.

METHODS

Spectra are specified using a rigorously-benchmarked functional form. Since ion chambers are the typical detector used in transmission measurements, the energy response of a Farmer chamber is calculated using the EGSnrc Monte Carlo code, and the effect of approximating the energy response on the accuracy of the unfolded spectra is studied. A proposal is introduced to enhance spectral sensitivity by combining transmission data measured with multiple detectors of different energy response and by combining data from multiple attenuating materials. Monte Carlo methods are developed to correct for nonideal exponential attenuation (e.g., scatter effects and secondary attenuation). The performance of the proposed methods is evaluated for a diverse set of validated clinical spectra (3.5-25 MV) using analytical transmission data with simulated experimental noise.

RESULTS

The approximations commonly used in previous studies for the ion-chamber energy response lead to significant errors in the unfolded spectra. Of the configurations studied, the one with best spectral sensitivity is to measure four full transmission curves using separate low-Z and high-Z attenuators in conjunction with two detectors of different energy response (the authors propose a Farmer-type ion chamber, once with a low-Z, and once with a high-Z buildup cap material), then to feed the data simultaneously to the unfolding algorithm. Deviations from ideal exponential attenuation are as much as 1.5% for the smallest transmission signals, and the proposed methods properly correct for those deviations. The transmission data with enhanced spectral sensitivity, combined with the accurate and flexible spectral functional form, lead to robust unfolding without requiring a priori knowledge of the spectrum. Compared with the commonly-used methods, the accuracy is improved for the unfolded spectra and for the unfolded mean incident electron kinetic energy by at least factors of three and four, respectively. With simulated experimental noise and a lowest transmission of 1%, the unfolded energy fluence spectra agree with the original spectra with a normalized root-mean-square deviation, %Δ(ψ), of 2.3%. The unfolded mean incident electron kinetic energies agree, on average, with the original values within 1.4%. A lowest transmission of only 10% still allows unfolding with %Δ(ψ) of 3.3%.

CONCLUSIONS

In the presence of realistic experimental noise, the proposed approach significantly improves the accuracy and robustness of the spectral unfolding problem for all therapy and MV imaging beams of clinical interest.

[Improvement of the accuracy of the Laplace transform method for the determination of radiotherapy spectra of clinical linear accelerators]

The present study focused on the reconstruction of the bremsstrahlung spectrum of a clinical linear accelerator from the measured transmission curve, with the aim of improving the accuracy of this method. The essence of the method was the analytic inverse Laplace transform of a parameter function fitted to the measured transmission curve. We tested known fitting functions, however they resulted in considerable fitting inaccuracy, leading to inaccuracies of the bremsstrahlung spectrum. In order to minimise the fitting errors, we employed a linear combination of n equations with 2n-1 parameters. The fitting errors are now considerably smaller. The measurement of the transmission function requires that the energy-dependent detector response is taken into account. We analysed the underlying physical context and developed a function that corrects for the energy-dependent detector response. The factors of this function were experimentally determined or calculated from tabulated values.

Monte Carlo calculation of nine megavoltage photon beam spectra using the BEAM code

A recent paper analyzed the sensitivity to various simulation parameters of the Monte Carlo simulations of nine beams from three major manufacturers of commercial medical linear accelerators, ranging in energy from 4-25 MV. In this work the nine models are used: to calculate photon energy spectra and average energy distributions and compare them to those published by Mohan et al. [Med. Phys. 12, 592-597 (1985)]; to separate the spectra into primary and scatter components from the primary collimator, the flattening filter and the adjustable collimators; and to calculate the contaminant-electron fluence spectra and the electron contribution to the depth-dose curves. Notwithstanding the better precision of the calculated spectra, they are similar to those calculated by Mohan et al. The three photon spectra at 6 MV from the machines of three different manufacturers show differences in their shapes as well as in the efficiency of bremsstrahlung production in the corresponding target and filter combinations. The contribution of direct photons to the photon energy fluence in a 10 x 10 field varies between 92% and 97%, where the primary collimator contributes between 0.6% and 3.4% and the flattening filter contributes between 0.6% and 4.5% to the head-scatter energy fluence. The fluence of the contaminant electrons at 100 cm varies between 5 x 10(-9) and 2.4 x 10(-7) cm(-2) per incident electron on target, and the corresponding spectrum for each beam is relatively invariant inside a 10 x 10 cm2 field. On the surface the dose from electron contamination varies between 5.7% and 11% of maximum dose and, at the depth of maximum dose, between 0.16% and 2.5% of maximum dose. The photon component of the percentage depth-dose at 10 cm depth is compared with the general formula provided by AAPM's task group 51 and confirms the claimed accuracy of 2%.

Computer algebra for x-ray spectral reconstruction between 6 and 25 MV

A previously presented algorithm for the reconstruction of bremsstrahlung spectra from transmission data has been implemented into MATHEMATICA. Spectra vectorial algebra has been used to solve the matrix system A * F = T. The new implementation has been tested by reconstructing photon spectra from transmission data acquired in narrow beam conditions, for nominal energies of 6, 15, and 25 MV. The results were in excellent agreement with the original calculations. Our implementation has the advantage to be based on a well-tested mathematical kernel. Furthermore it offers a comfortable user interface.