Reduction of uncertainties in radiotherapy assessed by Monte Carlo simulation: spectral analysis applied to absorbed dose correction

Abstracts

OBJETIVO: Determinar, por simulação Monte Carlo, os espectros de feixes de cobaltoterapia em profundidade na água e fatores de correção para doses absorvidas em dosímetros termoluminescentes de fluoreto de lítio. MATERIAIS E MÉTODOS: As simulações dos espectros secundários da fonte clínica de cobalto-60 foram realizadas com o código Monte Carlo PENELOPE, em diversas profundidades na água. Medidas experimentais de dose profunda foram obtidas com dosímetros termoluminescentes e câmara de ionização em condições de referência em radioterapia. Os fatores de correção para os dosímetros termoluminescentes foram obtidos através da razão entre as absorções relativas ao espectro de baixa energia e ao espectro total. RESULTADOS: A análise espectral em profundidade revelou a existência de espectros secundários de baixa energia responsáveis por uma parcela significativa da deposição de dose. Foram observadas discrepâncias de 3,2% nas doses medidas experimentalmente com a câmara de ionização e com os dosímetros termoluminescentes. O uso dos fatores de correção nessas medidas permitiu diminuir a discrepância entre as doses absorvidas para, no máximo, 0,3%. CONCLUSÃO: Os espectros simulados permitem o cálculo de fatores de correção para as leituras de dosímetros termoluminescentes utilizados em medidas de dose profunda, contribuindo para a redução das incertezas associadas ao controle de qualidade de feixes clínicos em radioterapia.

Simulação Monte Carlo; Radioterapia; Controle de qualidade; Espectrometria; Dosimetria termoluminescente; TLD

ORIGINAL ARTICLES

Reduction of uncertainties in radiotherapy assessed by Monte Carlo simulation: spectral analysis applied to absorbed dose correction^*

Tatiana Marques^I; Mirko Alva-Sánchez^I; Patrícia Nicolucci^II

^IPhD Fellow degree, Researchers at Departamento de Física e Matemática da Universidade de São Paulo (DFM-USP), Ribeirão Preto, SP, Brazil

^IIPhD, Professor at Departamento de Física e Matemática da Universidade de São Paulo (DFM-USP), Ribeirão Preto, SP, Brazil

Mailing address

ABSTRACT

OBJECTIVE: To calculate spectra of cobalt-60 beam at water depth and correction factors for absorbed dose measurements obtained with lithium fluoride thermoluminescent dosimeters using Monte Carlo simulation.

MATERIALS AND METHODS: The simulations of secondary spectra of clinical cobalt-60 sources were performed with the PENELOPE Monte Carlo code at different water depths. Experimental measurements of deep doses were obtained with thermoluminescent dosimeters and ionization chamber under reference conditions for radiotherapy. Correction factors for the thermoluminescent dosimeters detectors were obtained through the ratio between the relative energy absorption for the low energy spectrum and the total spectrum.

RESULTS: Deep spectral analysis has demonstrated the presence of secondary low-energy spectra responsible for a significant portion of the dose deposition. Discrepancies of 3.2% were observed among the doses measured with ionization chamber and thermoluminescent dosimeters. The adoption of correction factors has allowed a reduction in the discrepancy among absorbed doses to a maximum of 0.3%.

CONCLUSION: Simulated spectra allow the calculation of correction factors for reading of thermoluminescent dosimeters utilized in the measurement of deep doses, contributing for the reduction of uncertainties associated with quality control of clinical beams in radiotherapy.

Keywords: Monte Carlo simulation; Radiotherapy; Spectrometry; Quality control; Thermoluminescent dosimetry; TLD.

INTRODUCTION

Quality assurance in radiotherapy must associate high accuracy in the planning and prescribed doses with reproducibility of the planned technique, detailed documentation and a careful dosimetry of both the planned situation and treatment^(1–3). Characteristics such as absorbed dose homogeneity and accuracy are absolutely relevant for the success of clinical treatments in cobalt therapy units, particularly in whole body irradiations, which are normally combined with intensive chemotherapy and bone marrow transplantation. Small errors in the determination of depth doses may compromise vital organs, for example, the lungs^(4,5).

The determination of uncertainties associated with radiotherapy is absolutely relevant for the assurance of efficiency of clinical treatments. A considerable fraction of such uncertainty resides in intrinsic characteristics of the methods and instrumental apparatuses utilized in the determination of dosimetric parameters, which are directly applied in the calculation for prescription of absorbed dose. In order to minimize the uncertainty associated with dosimetric parameters, it is possible to utilize the Monte Carlo simulation of ionizing radiation interaction with matter in order to identify, quantify and correct inaccuracies from the dosimeters utilized in radiotherapy quality control.

Thermoluminescent dosimeters (TLDs) are among the dosimeters most frequently utilized in radiotherapy, with characteristics of effective atomic number (Z_effective = 8.2) close to that of water, high spatial detector resolution (0.9 mm) and low cost associated with their use^(6–8). TLDs with lithium fluoride (LiF) present energetic dependence for photons with energies < 662 keV⁽⁹⁾, which implies the necessity to correct readings performed with such dosimeters as they are exposed to low energy radiation. In measurements of depth absorbed dose⁽¹⁰⁾, for example, those applied in quality control of whole body irradiation⁽⁸⁾, the spectrum of therapy beams at depths of clinical interest is composed by high energy fractions, to which LiF presents a linear response, and also by low energy fractions, in which there is energetic dependence⁽⁹⁾.

The simulation of clinical treatments based on the Monte Carlo method has been widely utilized in radiotherapy mainly due to the accuracy associated with its use⁽¹¹⁾. The Monte Carlo PENELOPE simulation code is a well established and efficient tool for spectrometry and clinical beam dosimetry, allowing the faithful geometrical representation of several types of radiotherapy apparatuses⁽¹²⁾. The PENELOPE code uses material cross section libraries based on international standards⁽¹³⁾, thus assuring a faithful representation of characteristics interaction between materials of dosimetric interest.

The present study presents the energy spectra of a cobalt-60 (Co-60) clinical beam simulated in water depth, using the Monte Carlo PENELOPE simulation code. By means of quantitative analysis of dose contribution of primary and secondary components of such beam, it is possible to calculate correction factors for LiF-100 TLDs readings at several water depths. The contribution of the use of correction factors in absorbed doses in TLDs can be attained by comparison between the depth dose rate (DDR) curves built with experimental data collected with an ionization chamber and with TLDs, whose responses were corrected by these factors.

MATERIALS AND METHODS

Monte Carlo simulation

The computer simulations of ionizing radiations interaction with matter were performed by using the Monte Carlo PENELOPE simulation package, which comprises the PENELOPE particles transport code and the geometry package PENGEOM. In these simulations, the source was positioned at 80 cm from the phantom, with an aperture of 3.57°, which provided a radiation field of 10 × 10 cm² on the surface of the phantom. The primary photon spectrum is composed by two peaks of equal occurrence probability, with 1.17 MeV and 1.33 MeV⁽¹²⁾ energies, respectively. The phantom was a homogeneous water cube with a 30 cm edge.

The photon spectra at depths were obtained by means of virtual impact detectors: virtually drawn pellets, simulating the TLDs experimental geometry at depth inside the phantom, with 150 energy channels uniformly distributed between 10 KeV and 1.33 MeV. The impact detectors have the function of counting the photons that reach the depth in which the detector is positioned, separating into channels and storing their respective energies⁽¹³⁾. So, the response generated at the end of each simulation, comprises the probability density function associated with the photon occurrence, separated by energy, at each depth. Based on the data in the response files, the total and low energy spectra were built at several depths of the phantom. Also, a simulated DDR curve in the central axis of the phantom was obtained, for determining the PENELOPE code accuracy in this type of simulation, by means of comparison with experimental data of depth dose.

Experimental data

The DDRs were experimentally obtained with a Farmer type 0.6 cm³ cylindrical ionization chamber with LiF-100 TLDs, both positioned at depth, within a cubic, homogeneous water phantom, with an edge of 50 cm, analogous to that in the computer simulation. The field projected on the surface of the phantom was 10 × 10 cm², and was at a source-surface distance (SSD) of 80 cm.

The phantom was irradiated in a Siemens Gammatron II S-80 cobalt therapy system, at the Radiotherapy Unit of Hospital das Clínicas da Faculdade de Medicina de Ribeirão Preto da Universidade de São Paulo. The TLDs were treated and read at the Center of Instrumentation, Dosimetry and Radioprotection of Universidade de São Paulo (Cidra-USP).

Determination of correction factors by energetic dependence

For each depth, the summatories of spectrum energies multiplied by the ratio between the respective LiF and water mass absorption coefficients (µ_en/ρ)⁽¹⁴⁾ were calculated in two parts: one with energies in the 10 keV to 670 keV interval and another comprising all the energies in the spectrum, i.e., from 10 keV to 1.33 MeV. The calculation of the correction factors was made as presented in equation 1.

In the equation 1 (µ_en(E_i)/ρ) corresponds to the ratio between the LiF and water mass absorption coefficients (µ_en/ρ) for an energy E_i; P(E_i) is the probability density of the occurrence of particles with energy E_i at the x depth, relative to the total spectrum; P_sec(E_i) is the probability density of occurrence of particles with energy E_i, renormalized for the secondary spectrum at the x depth; FC_x is the correction factor for the TLD reading calibrated in dose at the x depth of the phantom.

RESULTS

Spectra were obtained at 24 different depths in the phantom, at 2 mm intervals in the first centimeter to observe the behavior at the build-up region, and at every 1 cm in depth after the build-up region, scanning the whole phantom extent. The total spectrum simulated at a depth of 0.5 cm shown on Figure 1, presents a low energy spectral component with a peak at 220 keV (2.5%) with a total expression rate of 5.5%.

The secondary spectra simulated for the depths of 5.5 cm, 10.5 cm 20.5 cm and 25.5 cm along the phantom central axis are presented on Figure 2. The relative intensities of the simulated spectra reach a maximum at the 5.5 cm depth, where the low energy spectra peak represents 7% of the total spectrum. After this depth, the contribution of the energy peak, centered at 220 keV, goes to 5% at 10.5 cm, decreasing to 1.4% and 1.2% at 20.5 cm and 25.5 cm respectively. However, the contribution in dose deposition increases, as the relative expression of energies > 220 keV grows.

In order to obtain the FC_x, the depths at which there were DDR data with both dosimeters (TLD and ionization chamber), were considered. Figure 3A graphically presents the DDR values measured by the TLDs, with and without response correction, together with the DDR measurements performed with the ionization chamber and the DDR curve calculated by means of integration of simulated energy spectra with the PENELOPE code. The percentage differences between the absorbed dose readings with the TLDs, with and without the application of the FC_x, and measurements with the ionization chamber are presented on Figure 3B.

The chart presented on Figure 3A demonstrates the PENELOPE code as appropriate and accurate for simulating spectra in water depth, co-validating the spectral analysis presented in the present study. Based on the analysis of Figure 3, it is possible to observe how much the readings with TLDs become less accurate at greater depths, because of the presence of low energy radiation, which confirms the need to apply the FC_x. On Figure 3B, it is clear that the FC_x approach the measurements performed with the TLDs to those performed with the ionization chamber, reducing discrepancies from up to 3.3% to a maximum of 0.38%. The numerical data of the FC_x can be observed on Table 1.

DISCUSSION

The clinical practice of radiotherapy with high quality indexes is primarily based on rigorous, periodic and reproducible reference dosimetry^(14,15) as corroborated by the utilization of TLDs as a complementary alternative to reference dosimetry, allowing more frequent quality control tests with relative simplicity, without excessive cost and interference in the clinical routine⁽¹⁶⁾.

The reduction of uncertainties and errors associated with dosimetry in radiotherapy is the central theme of several scientific and clinical studies, all of them facing the challenge the fact that precision and accuracy are not easily found in a single dosimeter⁽¹⁷⁾. The use of Monte Carlo simulation, combined with experimental measurements with TLDs for directly determining dosimetric parameters is a precise and viable alternative to face the need of reducing uncertainty and errors in radiotherapy. With the correction factors presented in the present study, it is possible to utilize LiF TLDs for DDR measurements in water, with an error interval relative to the ionization chamber coherent with the indexes recommended by the reference literature on quality control of radiotherapy beams.

The spectral analysis based on total and secondary photon spectra at water depth, simulated with the PENELOPE code^(18,19), allows the calculation of correction factors that can be applied in a relatively simple manner. Photons spectra in water depth allow to observe the presence of low energies at the maximum dose depth, distributed between 30 and 350 keV, region in which the LiF response has a non-linear behavior alternating between a relatively descending response at energies from 30 to 100 keV and mildly increasing at energies in the 100 to 300 keV interval, approximately.

In all simulated spectra up to 15.5 cm within the phantom, the low energy peak is kept at 220 keV and the expression of portions with energies from 30 to 300 keV is predominant, which does not happen at depths below the middle of the phantom, at which the probability of occurrence of low energy particles are equally distributed in a larger interval, from 30 to 670 keV.

As the depth increases, the total spectra present a higher number of lower energy particles, which dislocates the mean energy of the simulated beam reaching that depth to a lower value than the mean energy of the incident Co-60 beam, that would be 1.25 MeV, commonly considered for the calculation of absorbed energy at water depth.

It is possible to observe that the area of secondary spectra with energies between 10 keV and 670 keV remains approximately the same after the build-up depth, therefore dislocating the effective energy of the secondary spectrum which increases the dose in each detector. Such an increase is not duly observed in experimental measurements performed with LiF TLDs, as a function of their response energetic dependence in this energy range, which can be observed in the DDR curves obtained with these dosimeters (see Figure 3A).

In spite of the decrease in energy with depth be easily visualized from the observation of simulated spectra, a more precise study of absorbed dose at a given depth must consider the (µ_en/ρ) corresponding to each discrete set of particles with a given energy, that are, by their turn, differently absorbed by water and by the LiF, generating different deposited doses in each one of these mediums. It is exactly by observing the low energy portions and their significance at each depth that it is possible to relate the differences in absorption among the medium of interest, the water, and the dosimeter material, LiF, applying the FC_x determined in the present study.

When compared, the PDDR curves obtained with an ionization chamber and with TLDs presented a percentage difference that increases with depth, from 0.2% to 3.2%, and remaining approximately constant, at around 2.5% for depths > 10 cm. Also, when observing that the relative probability of the low energy particles increase with depth in the phantom (analysis on Figure 2), and considering that the total number of simulated particles in each interaction remains the same, it is easy to observe that the FC_x must increase with depth, as shown on Table 1.

CONCLUSION

The spectral analysis of a clinical Co-60 beam based on the Monte Carlo simulation allows the calculation of satisfactory correction factors for the absorbed dose readings in LiF TLDs. The application of the correction factors FC_x in TLD readings contributes for the increase in accuracy in the determination of dosimetric parameters, intimately associated with the quality control of clinical beams in radiotherapy.

The results of the present study may be extended for several clinical protocols, with the objective of making the use of TLDs more accurate and common.

Acknowledgments

Authors thanks to Setor de Radioterapia do Hospital das Clínicas da Faculdade de Medicina de Ribeirão Preto (HCFMRP-USP) and to Centro de Instrumentação, Dosimetria e Radioproteção da Universidade de São Paulo (CIDRA-USP) for essencial technical support.

REFERENCES

• 1. Kawa-Iwanicka A, Lobodziec W, Iwanicki T, et al. Dose uniformity in the total body irradiation technique using 15 MV photon beam. Physica Medica. 2004;20(Suppl 1):1446.
• ²
  International Atomic Energy Agency. Absorbed dose determination in external beam radiotherapy. An International Code of Practice for Dosimetry Based on Standards of Absorbed Dose to Water. Technical Reports Series No. 398. Vienna: International Atomic Energy Agency; 2000.
• ³
  International Commission on Radiation Units and Measurements. Measurement of absorbed dose in a phantom irradiated by a single beam of X or gamma rays. ICRU Report 23. Bethesda: International Commission on Radiation Units and Measurements; 1973.
• 4. Vrtar M. A dosimetric method of total body irradiation. Cell Mol Biol Lett. 2002;7:33740.
• 5. Zabatis Ch, Koligliatis T, Xenofos S, et al. Dosimetry in translation total body irradiation technique: a computer treatment planning approach and an experimental study concerning lung sparing. J Buon. 2008;13:25362.
• 6. Bloemen-van Gurp EJ, Mijnheer BJ, Verschueren TA, et al. Total body irradiation, toward optimal individual delivery: dose evaluation with metal oxide field effect transistors, thermoluminescence detectors, and a treatment planning system. Int J Radiat Oncol Biol Phys. 2007;69:1297304.
• 7. Kron T. Thermoluminescence dosimetry and its applications in medicinePart 1: Physics, materials and equipament. Australas Phys Eng Sci Med. 1994;17:17599.
• 8. Giordani AJ, Segreto HRC, Segreto RA, et al. Verificação das doses de radiação absorvidas durante a técnica de irradiação de corpo inteiro nos transplantes de medula óssea, por meio de dosímetros termoluminescentes. Radiol Bras. 2004;37:3439.
• 9. Reyes FE. Procedimento de calibração de dosímetros termoluminescentes em dose absorvida na água para fontes de irídio-¹⁹²I de alta taxa de dose [dissertação de mestrado]. Rio de Janeiro: Instituto de Radioproteção e Dosimetria; 2004.
• 10. Almeida CE, Affonseca M, Calcina CSG, et al. Rastreabilidade das referências metrológicas em dose absorvida na água do Programa de Qualidade em Dosimetria. Radiol Bras. 2005;38:2058.
• 11. Garnica-Garza HM. Monte Carlo-derived TLD cross-calibration factors for treatment verification and measurement of skin dose in accelerated partial breast irradiation. Phys Med Biol. 2009;54: 162131.
• 12. Panettieri V, Sempau J, Andreo P. Chamber-quality factors in ⁶⁰Co for three plane-parallel chambers for the dosimetry of electrons, protons and heavier charged particles: PENELOPE Monte Carlo simulations. Phys Med Biol. 2008;53: 591726.
• 13. International Commission on Radiation Units and Measurements. Tissue substitutes in radiation dosimetry and measurement. ICRU Report 44. Bethesda: International Commission on Radiation Units and Measurements; 1989.
• 14. Morávek Z, Rickhey M, Hartmann M, et al. Uncertainty reduction in intensity modulated proton therapy by inverse Monte Carlo treatment planning. Phys Med Biol. 2009;54:480319.
• 15. Guckenberger M, Krieger T, Richter A, et al. Potential of image-guidance, gating and real-time tracking to improve accuracy in pulmonary stereotactic body radiotherapy. Radiother Oncol. 2009;91:28895.
• 16. McGarry CK, Cosgrove VP, Fleming VA, et al. An analysis of geometric uncertainty calculations for prostate radiotherapy in clinical practice. Br J Radiol. 2009;82:1407.
• 17. Budrukkar A, Dutta D, Sharma D, et al. Comparison of geometric uncertainties using electronic portal imaging device in focal three-dimensional conformal radiation therapy using different head supports. J Cancer Res Ther. 2008;4:706.
• 18. Salvat F, Fernández-Varea JM, Acosta E, et al. PENELOPE a code system for Monte Carlo simulation of electron and photon transport. Paris: Nuclear Energy Agency; 2008.
• 19. Salvat F, Fernández-Varea JM. Overview of physical interaction models for photon and electron transport used in Monte Carlo codes. Metrologia. 2009;46:S11238.

*

Trabalho realizado no Hospital das Clínicas da Faculdade de Medicina de Ribeirão Preto da Universidade de São Paulo (HCFMRP-USP), Ribeirão Preto, SP, Brasil. Suporte financeiro: Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) e Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Capes).

Publication Dates

• Publication in this collection
  07 May 2010
• Date of issue
  Apr 2010

History

• Received
  26 Mar 2009
• Accepted
  15 Oct 2009