Open access peer-reviewed chapter - ONLINE FIRST

Thermoluminescence Dosimetry Technique for Radiation Detection Applications

Written By

Hossam Donya

Submitted: December 29th, 2021 Reviewed: January 17th, 2022 Published: March 25th, 2022

DOI: 10.5772/intechopen.102728

Dosimetry Edited by Thomas J. FitzGerald

From the Edited Volume

Dosimetry [Working Title]

Dr. Thomas J. FitzGerald

Chapter metrics overview

28 Chapter Downloads

View Full Metrics


Due to the risk of radiation exposure, radiation dosimetry is performed regularly to ensure the occupational safety of personnel and radiation workers. Therefore, various dosimeters are widely used to detect neutrons, gamma, X-ray, and proton irradiation fields. As an example, in medical applications, routine personal dosimetry is used to monitor and limit workers’ long-term occupational exposure. Radiation workers who undertake X-ray diagnostic, radiotherapy operations, in clinical and industrial application. Although, the overheads of running an in-house TLD (Thermoluminescent dosimetry) service for monitoring doses to eyes, pacemakers and so on seems rather high for the benefits conferred, however, it is still widely used for reporting doses accurately in various medical centers over the world. TLD also is widely used for measuring entrance doses on a handful of patients to validate a new LINAC/TPS combination. As well as in the industrial field as if petroleum, companies or nuclear reactor, RSO (radiation safety officer) used TLD badges to report delivered doses. In this chapter, we focus on the TLD technique for measuring doses of various ionizing radiation detection. Different methods for evaluations of TL Kinetics are covered. Modern TLD applications in the clinical field are also investigated. Some recommendations on advance dosimetry failure of TLD are concluded.


  • TLD
  • TL kinetics
  • radiotherapy
  • hybrid-functional density theory
  • modern clinical applications

1. Introduction

TLD method is considered an important technique as it can store radiation in trap centers for long period. Glow peaks of thermoluminescence dosimeters are later measured and discussed based on some models related to the physical changes in the band structure of dosimeter because of ionizing radiation exposure. A wide range of substances exhibits thermoluminescence (TL) phenomena after being exposed to nuclear radiation such as activated LiF and CaSO4. Thermoluminescent dosimeter (TLD) emits light when heated up after being irradiated. Due to this special property, TLD could be used as a radiation dosimeter. TLD has many advantages and sensitive to different types of radiation. A dosimeter of higher TL response to thermal neutrons is most commonly used in mixed radiation fields (neutron and gamma ray). The sensitivity of TLD to neutrons depends on the detector compound type, environment and neutron energy. For neutron dosimetry purposes, the neutron and gamma ray dose contribution must be separated by using two different detector types of TLD. The first one should be sensitive to gamma and the other should be sensitive to neutrons plus gamma (as LiF-700 and LiF-600) [1, 2].

The response of fast neutrons depends on the cross-section for the interaction in TLD material and the relative TL efficiency, which depends on the linear energy transfer (LET) of the reaction products in the first place. The response to intermediate-energy depends mainly on the cross-section of the reaction, which may take place with the composite material of the TLD.

1.1 TLD applications in neutron and gamma ray dosimetry

Generally, there are three types of TLD used for neutron dosimetry as follow:

1.1.1 Albedo neutron dosimeter

A considerable fraction of intermediate and fast neutrons can be slowed down to epithermal neutron energy and backscattered in the human body, interacting with the sensitive TL material. An albedo neutron dosimeter is a type of neutron monitor and is typically used in the neutron energy range of 0.2 eV to around 0.5 MeV. The slow neutrons interact with TL material, usually through 6Li (n,α) 3H reaction, and the resulting induced charged particles to stimulate the TL material. Recently, some of albedo TLD dosimeters depend on 10B (n, α) 7Li reactions. Because neutron TL sensitive material responds to gamma radiation, and neutrons are accompanied by this gamma radiation, another TLD is usually utilized in conjunction with TLD with a gamma ray.

The neutron albedo dosimeter measures (a) direct fast neutrons, (b) direct thermal neutrons, and (c) albedo neutrons reflected from the body. This type of dosimeter uses Lexan polycarbonate and/or CR-39 foils, as well as two 10B (n, γ) 7Li converters in a cadmium cover, to efficiently measure the three neutron dosage components independently [3, 4, 5]. Fast neutron dose is assessed in CR-39 by counting proton recoil tracks, while thermal neutron dose is determined by counting α particles created during the process. Because the albedo dosimeter has a sensitivity range of 0.3–30 mSv, it is advised that it be used as a backup dosimeter to assist in the assessment of high dose values in the event of accidents or patients receiving neutron therapy.

In another application, the 10B (n, α) 7Li reactions with the backscattered albedo neutrons employed with Electret’s ionization chamber proposed by Seifert et al. [6, 7]. In this chamber, induced 7Li from the ionization of the gas in the chamber worn on the body’s surface in the above reaction instance. Under saturation conditions, produced charge carriers with the corresponding polarity travel to the surface of the electret. As a result, the change in the electrets voltage is a direct measure of albedo neutron fluence and an indirect estimate of primary neutron fluence. In general, the advantages of albedo TLD dosimeter are: they are relatively inexpensive and can be reused, easily fabricated, lightweight to wear, Readout is simple and can be automated, Insensitive to humidity.

While their disadvantages are: Some of TLD exhibit fading, TLD is sensitive to gamma-ray, they must be worm properly or serious errors can be resulted, the measured values of TLD does not give permanent record as the track detectors, their sensitivity is highly dependent on the angle and energy of incidence radiation.

1.1.2 Hydrogenous radiator TLDs

In this type of dosimeters, the fast neutrons knock out protons from hydrogenous material mixed with the phosphor, and the protons dispel their energy in the dosimeter. In this method, the hydrogenous substances are called proton radiators [8]. This technique has demonstrated that TL materials mixed with hydrogenous material can detect fast neutrons, but the sensitivity needs to be improved by one order of magnitude before using in personnel neutron dosimetry.

1.1.3 LET-dependent deep trap TLD glow peaks

The fast neutron interacts directly with the TL material as calcium fluoride (CaF2:Tm) which is commercially called TLD-300. This type has a glow curve with two glow peaks and the peak temperature Tm centered 150 and 250°C, respectively. The higher temperature peak (250°C) has a greater response to the fast neutrons. TLD-300 dosimeter CaF2: Tm (0.35 Mol. %) showed a lower detection limit of about 0.3 mSv from 241Am-Be source.


2. Characteristic of TLD phosphors

2.1 The glow curve

The term “Glow curve” refers to the graph of TL as a function either of temperature or of time as shown in Figure 1.

Figure 1.

TLD glow curve and time–temperature profile (TTP).

Glow curves have the following features:-

  • The glow curve of a certain phosphor probably best characterizes that phosphor. For example, the appearance of glow peaks only at low temperatures implies that the phosphor loses its stored TL with time, and therefore would be unsuitable for long-term measurements.

  • A peak at very high temperatures indicates a phosphor that will produce infrared radiation at the temperature necessary to release the TL. This is a problem for the instrumentation.

  • A glow curve without well-defined peaks makes the selection of the appropriate end for integration difficult. Ideally then, a glow curve should reveal only a single thermoluminescent peak, which occurs at a temperature high enough to ensure room temperature stability but not so high as to present instrumental problems. However, the temperature at which a TL peak appears is quite affected by the heating rate.

The following factors may affect the shape of the glow curve:

2.2 TL sensitivity

The sensitivity of TLD should be evaluated for competitor’s materials to determine the dose linearity, sub linearity or supralinearity behavior of phosphors. The sensitivity and thermal stability of glass samples were found to be dependent on both the starting materials and the method of preparation in Refs. [9, 10, 11, 12, 13].

2.3 Dose rate dependence

TL dosimeters, in general, have demonstrated no dose rate effects over a wide range.

2.4 Stability

TL is the release, in the form of visible light, of energy absorbed from previous irradiation. The rate at which this energy is released is dependent upon the phosphor temperature and increases sharply at a higher temperature. Even though the concept of “glow peak temperature” that temperature at which the maximum of the glow peak occurs, is both useful and easily demonstrated. It should be remembered that a finite rate of loss of stored energy exists even for much lower temperatures. If a phosphor shows an insignificant loss of TL at room temperature, it is said to have good stability.


3. TL kinetics

The physical process leading to the emission of TL from a sample is related in most cases to the traffic of charge carriers, usually electrons and holes, between different imperfection states in the solid sample. Studying the kinetics of the TL process means the investigation of electron–hole transitions between energy states in cases of both the irradiation of the TL sample and the readout processes. Although, in most experimental situations the TL curve consists of several overlapping peaks, it is appropriate to start the discussion by dealing with a single peak to understand the basic process. For most purposes, it is not necessary to assume that the glow curve consists of only one peak. The analysis of a single peak may just be valid if a series of peaks occur, provided that the peak of interest is sufficiently separated from others, either because it appears separately or because we have an efficient method to isolate it from the rest of the curve.

3.1 First kinetics order

If n is the number of trapped electrons in the sample, which is maintained at constant temperature T, n decreases with time t as:


From Eq. (1) we get:


The rate of photon emission, and hence the rate of release of electrons from traps to their rate of arrival at luminescence centers, determine the strength of the TL glow peak [14].


Where Cis a luminescence efficiency constant.

When the dosimeter is heated with rate β=dT/dt. Then we may write dn/dtas:


Substitute in Eq. (2), we get:


Where no denotes the number of electrons in the trap at a given time and temperature, to and To, respectively.


Substitute in Eq. (3)


It should be noted that I(T) in Eq. (6) dependents on two physical parameters, the activation energy E, and frequency factor S, and the heating rate β. The activation energy is the minimum energy required to release the electrons from their traps. Differentiation of Eq. (6) with respect to the temperature gives:


Where, fT=SβTexpEKTdT

At T=TmIT=o


where Tmis peak position or the temperature at maximum intensity.

Equation (7) describes the condition of the occurrence of the maximum intensity and the determination of the corresponding temperature, which we call, Tm. The reduction in the second exponential function is faster than the growth in the first exponential function above this temperature, and the product function decreases until the traps are fully depopulated. This accounts for the end of the peak. A theoretical (calculated) glow peak plotted using Eq. (6) is shown in Figure 2. The main feature of the first-order peak is that the asymmetric, is such that at temperatures over Tm, the reduction is faster than the rise at low temperatures.

Figure 2.

Theoretical glow peak plotted using the first kinetics order equation.

The initial concentration noappears in the first kinetics order acts only as a constant multiplying the temperature-dependent factors. In this particular case of the first kinetics order, changing the initial concentration no has no effect on the curve’s form because adjusting the intensity at each temperature has the same proportional effect. Figure 3 shows several glow peaks with different no. One of the aspects of this fact is that Tmis independent of the initial concentration no.

Figure 3.

Glow curves plotted using the first-order kinetics equation for differentno.

This appears well in the condition of the Tmdescribed by Eq. (7), where no does not appear in the equation. This property of the independent of Tmon nois specified to the first-order case, and will not occur for most of the other kinetics possibilities [13]. Eq. (7) can be written in the following form:


We see that changing the heating rate βmust change Tmin a such way that equality still holds. The term Tm2expEKTis monotonically increasing with Tm, therefore increasing the rate βwill immediately cause Tmto increase. Since Tm2expEKTis a very rapidly increasing function of Tm, only a small change of Tmmay accompany a large variation in the heating rate β, this variation is usually rather easily observable.

3.1.1 Second kinetics order

One assumption made up by Randall and Wilkins [14] which led to the first kinetics order was that once a charge carrier is thermally elevated into the band, it is bound to recombine rather quickly with an opposite sign carrier trapped in a recombination center. Gralick and Gibson [15] considered another case in which the free carriers may re-trap with equal retrapping recombination probabilities with the further assumption that the concentration of electrons in traps and holes in recombination centers are equal during the entire process. Denoting the total number of traps of the given type (free electrons or holes) by N, they found the kinetics equation:


where (S/N) is a constant having units of m3s−1, which we may denote by S′. Then we have


where S′is called “pre-exponential factor” which does not have the same meaning of “frequency factor” as was in the first kinetics order.

For linear heating rate β, we have:


Substitute in Eq. (10),


where Eq. (12) represents the intensity of a glow peak according to the second kinetics order model. At high temperature, the second decreasing function dominates so that the product function is decreasing. Somewhere between two regions the glow curve, therefore, reaches its maximum. Figure 4 Displays a hypothetical glow peak plotted using Eq. (12).

Figure 4.

Theoretical glow peak plotted using the second-order kinetics equation.

The condition of the maximum is found by setting the derivative of Eq. (12) to zero (dI/dT = zero) [16], then we may find:


Multiply by 1+S´βnoTmexpEKTmdT3and rearrange, one gets


Then Eq. (13) represents the condition of the peak maximum according to the second kinetics order. As can see noappears in the equation and therefore we expect that Tmwill depend on no. It can be shown numerically or analytically, that increasing nocauses Tmto decrease. An exception to this rule of the shift of a second-order peak with nocan be found by Wrzesinska [17], who writes Eq. (10) with S=Sno. The resulting peak has all the regular features of a second-order peak (e.g., symmetry properties) except one can write Sinstead of noS’ and thus Eq. (10) turns out to be independent of no. The ensuring Tmis, therefore independent of no. It is not clear, however, what physical circumstances result in S′ being equal to S/no[17]. Other aspects of the dependence of the glow curve on the initial concentration noare paramount importance when we are interested in a TL as a dosimetric tool. In many cases, one associated the initial concentration with the imparted dose and then the dependence of different parts of the glow peak on nois important. In the first kinetics order, since the intensity at each point is multiplied by the same factor while changing no, the total area varies with the same amount so that the total area is proportional to no. Its occurrence in second-order peak can be illustrated by integrating Eq. (9) with respect to time from zero to infinity;


Both in the first order and second order, as well as other cases, n is zero and therefore the integral, which represents the area under the glow peak is equal (in appropriate units) to no.

Now we can consider the dependence of different portions of the second-order peak on no. First, we shall study the dependence of I on no for a given temperature T. In the initial rise range, Eq. (12) reduces to:


This shows immediately that for a given temperature in this range the dependence of Ion nois superlinear, namely I α no. It is to be emphasized that it is true only in the initial rise region; as already shown the total area is proportional to noand different dependencies are expected on other portions of the curve. Using the maximum condition equation and approximation to TexpE/KTdT, it can be shown that the two terms in the brackets in Eq. (12), namely unity and noSβToexpE/KTdTare more or less equal at T = Tm. At higher temperature, the latter term increases substantially and the unity can be neglected so that we obtain:


The main point in Eq. (18) is that the term includes nocancel. This means that at a higher temperature range the TL intensity is independent of nofor any given temperature [17].

Figure 5 shows plotted glow peaks using Eq. (12) for different no. In the low-temperature range, the TL intensity appears to depend on no. As noincreases, Tmdecreases which makes the peaks appear to be shifted to the low-temperature side. As the temperature increases the effect of noon the peak starts to decrease which makes the peaks approach each other’s on the high-temperature side.

Figure 5.

Plotted glow peaks using the second-order kinetics equation for differentno.

3.1.2 A single TL peak analysis

As seen in Figure 6, the concentration of the trapping state is denoted by N(m−3), with n(t) (m−3) being filled by electrons at time t(s). These electrons can be thermally elevated into the conduction band by crossing an energy barrier of E (eV) at a rate proportional to exp.(−E/kT), resulting in a concentration of free electrons nc(t) (m−3). Following that, these can be retrapped in a similar trap with a re-trapping probability Anor recombined with a trapped hole in a center with a recombination center probability Am, generating a photon with the recombination center energy h. A set of three simultaneous differential equations governs this operation. The following factors influence the recombination process:


where n, m, and ncare the trapped electron, hole in the center, and free-electron concentrations, respectively, and (dm/dt) is the recombination rate. This means that the amount of light emitted is proportional to the pace at which m decreases. The rate of recombination is proportional to both the instantaneous concentration of free electrons ncand the concentration of hole centers m, the proportional constant Am (m3 s−1). The product of cross-section recombination σ(m2) and thermal velocity is commonly used to calculate this value (m.s−1). The second equation is concerned with the movement of electrons that have been thermally liberated from the trapped condition. The rate of release of these electrons –dn/dtis proportional to the trapped electron concentration n (m−3) and the Boltzman constant exp.(−E/KT), with S serving as the proportional constant (s−1).

Figure 6.

A general treatment of the charge carriers’ transitions in the TL sample.

However, the actual rate of change of n is also related to the retraping term. The rate of retraping is proportional to the concentration of free electrons nc, and the unoccupied trapping states N-n, the proportional factor being the recombination probability An(m3s−1). Thus, the second combined equation is given by:


The third equation is that of charge neutrality. In its simplest form, it should read m = n + nc. Taking the first derivative with respect to time, the charge neutrality condition can be written as:


This equation has been given by Adirovitch [18] for phosphorescence and by Halperin and Braner [19].

Now let us discuss the kinetics of the process in more general terms and see how the simplified cases of first, second, and more general cases emerge from Eqs. (17)(20). Two simplifying assumptions were first made by Adirovitch [18] and later by many other investigators [19, 20, 21, 22, 23]. These are related to the relation between the concentration of the electrons in the conduction band and in traps and to the rate of change of these concentrations, namely:


Although, it seems to be the same connection between these two conditions, basically they are two separate relations and the occurrence of one does not necessarily imply the other. With these assumptions, Helperin and Braner [19] found the expression:


Since this equation contains two unknown functions, n(t) and m(t), it cannot be solved without further assumption. As mentioned, Randall and Wilkins [14] wrote their first-order equation assuming strong recombination. This can be expressed in more specific terms. If we assume with relation to Eq. (22) that:


The condition of Eq. (22) is the relation between functions rather than parameters. It is, therefore, possible that at the low-temperature range of a glow peak, the strong inequality holds, and at higher temperatures where m and n decreases, the inequality “weakens” may be inverted. This may result in a shift from first-order behavior to non-first-order behavior within the same peak [24].

Then, one can say that:


Then Eq. (22) will take the following form:


Then we see that Eq. (23) takes the same form of Eq. (3). For linear heating rate function, the general solution of Eq. (23) is given by Eq. (24):


Then from Eq. (22) with Randall and Wilkins [14] assumptions, we reached the first kinetics order equation.

The abovementioned second kinetics order, resulting from different assumptions associated with Eq. (22). In one set of assumptions, one can take n(t) = m(t) which is not very different from the parametric equality no = moonce the assumption ncnis made.

In addition, we have to assert the retraping dominates [15]


We also suppose that the trap is far from being saturated, i.e., the retrapping duration.


Then, from Eq. (26) one can write:


Using the condition of Eq. (26) in Eq. (27) one gets:


and since we have assumed that n(t) = m(t),


Then, apply to Eq. (22), we get:


Alternatively, one can assume, in addition to the concentration equality, that An = Am[18] which yields:


Then Eq. (22) takes the same form of Eq. (8) which is found by Gralick and Gibson [15]. Where Eq. (30) sums up both these possibilities by employing the parameter S′ (m3s−1), the pre-exponential factor that replaces AmSAnNin one case and S/N in the other. The solution of Eq. (30) is given by Eq. (32)


It should be emphasized that two cases discussed so far, namely first and second kinetics order, are only special cases in a sense, extreme cases and the general case described by equations Eq. (17) through Eq. (19) may be neither first nor second order even if the simplifying conditions of Eq. (21) are assumed to be general. The resulting Eq. (20) consists of many intermediate cases that do not have a distinct kinetics order. Although, some researchers still attempt to determine for every TL peak a first or second kinetics order [25].

Several attempts [16, 26] have been made to add a third parameter to the two basic ones, the activation energy Eand the pre-exponential constant S′ (or S), all the attempts extend the “order parameter” implied when talking about first or second-order peak. The order parameters considered so far as a discrete magnitude assuming the value of 1 and 2 can be extended to be a continuous parameter. It is to be noted, however, that the addition of a third parameter is in principle one step in the right direction since the general treatment should include eight parameters (E, S, Am, An, N, no, mo, nc). The best-known way of including the third parameter is that of general kinetics order, b, according to which one can assume that the glow peak is governed by [25].


The kinetics order, b, is normally considered to be between 1 and 2, but it can occasionally exceed this range [13]. The rationale behind writing Eq. (33) is as follows: it is readily seen that a first-order peak is asymmetric, where a second order peak is nearly symmetric. Following Halperin and Braner [19] and Chen [16] we can define the symmetry factor μgas:


where δ=T2Tm,ω=T2T1as it is shown in Figure 7, and T1 and T2 are the low and high temperatures on half- maximum intensity, respectively. It has been shown [16] that for the first kinetics order, μg0.42and the second kinetics order, μg0.52.

Figure 7.

Parameters used in the calculation of the symmetry factor.

Of course, intermediate symmetries represented by different values of μgare found and the simplest way to present them by taking 1<b<2in Eq. (30). Chen [16] has shown that μgchanges from 0.42 to 0.52 as b increasing from 1 to 2. The solution of Eq. (33) for linear heating rate β, is given by:


where S=Snob1.Eq. (35) represents glow peak intensity according to the general kinetics order.

A few words of caution are in order with respect to this treatment. First, although Eq. (34) has been shown to quite accurately described measured TL peaks [27, 28], it is to be noted that in most cases it is only an empirical presentation and is not based on the three differential equations [Eqs. (17) up to (19)], seem to be more physically significant. However, the general order case is still important because it can handle intermediate circumstances and smooth the first and second-order cases as b1 and b2, respectively.

3.1.3 General-order kinetics

May and Partridge supposed the empirical equation that has been suggested to explain the thermoluminescence glow peak if the first or second-order kinetics do not describe the glow peak. The equation is namely the general- order kinetics and written by:


Hence s=sNn0is called the pre-exponential factor, bthe order of kinetics and the rang supposed between 1 and 2 but sometimes this rang has able to be greater than those. The pre-exponential factor sis constant for given the dose, however, it differs with changing the absorbed dose withn0.

3.1.4 Trap parameters evaluation techniques Empirical methods

We can deduce that the higher the peak temperature Tm, the higher the activation energy Urbach [29], and Urbach [30] found empirically for KCl crystals:


This can also be written as E = 23KTm and it differs according to the types of the sample. Halperin [19] deduced E = 38 KTm for NACL samples, and Miller and Bube [31] arrived at E = 39 KTm for LiF.

The maximum intensity of the peak, according to Randall and Wilkins [12, 13], occurs around the temperature where the electron escape probability is 1 s-1. As a result of Eq. (1), we have:

P=SexpEKTm=1yieldsE=KTmlnSE37 Initial rise method

According to Eqs. (6), (12) and (3), we can say that at the start of the glow peak (initial rise region) the TL intensity is proportional to expE/kT, irrespective of whether the first kinetics order is obeyed or not [32]. This temperature relationship persists until the quantity of trapped electrons is drastically reduced. Hence, by plotting Log (I) versus 1/T, the value of Ecan be obtained from the slope of the straight line obtained. As a result, using the equation: it is possible to calculate Ewithout knowing the frequency factor S:


From Eq. (6), we see that when Tis slightly greater than Tm, the argument of the second exponential is very small and therefore the value of the exponential function is close to unity and varies very slowly with temperature. The temperature dependence of I(t) is therefore dominated by the first exponential function, however the second exponential function decreases with increasing temperature and at higher temperatures it decreases very rapidly [13].

Therefore, the range of the initial rise must be chosen in which the second exponential function has minimum influence on the TL intensity temperature dependence. Therefore, it is necessary to restrict the temperature range such that the TL intensity does not exceed one-tenth of the maximum intensity [32].

Between temperatures T1 and T2 (both < Tm) corresponding to values equal to a1Imand a2Imrespectively as in Figure 8, where:

Figure 8.

Extracted parameters from “Christodoulides expression” are to correct the value of the activation energy evaluated by the initial rise method.


On the temperature scale, a series of points were taken at equal intervals and plotted as ln(I) versus (1/T). The value Ec can then be calculated from the slope of the straight line as the energy determined by the initial rise technique; this value is smaller than the real activation energy Eby the amount that grows as a1 and a2 increase. Christodoulides [33] devised the following expression for the corrected energy Ein terms of the measured values Ec, a1, and a2:


The range of applicability of this equation is restricted by:

10EKTm100E41 Peak shape method

Grossweiner [34] established the first peak shape approach for first-order peaks, writing:


Where: Tmis the temperature at the maximum intensity, T1 is the temperature at the half of the maximum intensity in low-temperature side, τ=TmT1as in Figure 9. Grossweiner used the coefficient 1.51, which was later [20] amended to 1.41. Lushchik [35] developed a method for evaluating the activation energy by utilizing the high-temperature half width δ=T2T1for first peaks he suggested:

Figure 9.

Peak shape method used to calculate the activation energy.


and for second-order peak:


Chen [16] improved these equations by adding a factor of 0.976 in front of the former and replacing the factor 2 by 1.71 in the latter.

Halperin and Braner [19] have derived their equations for both first [Eq. (45)] and second kinetics orders [Eq. (46)]:


Chen [16] managed to establish expressions for general kinetics order, which is dependent on the geometry factor of the glow peak which is defined by Eq. (35):


Where φstands for τ,δ,ωand the values of Cφand bφfor the three methods are:


whereμgis geometrical shape factor that equalδω. Various heating rates method

As mentioned above about Eq. (7), Tm changes with the heating rateβ, writing Eq. (7) twice for heating rate β1 and β2 with maximum temperatures Tm1 and Tm2 we get [36]:


The activation energy that will be evaluated from Eq. (54) will be of course in accord with the first kinetics order only. However, Chen and Winer [37], Chen and Kirsh [38] showed that it can be used as a very good approximation for nonfirst-order cases as well.

The maximum condition, Eq. (16), can also take the following form:


According to this equation, Hoogenstraaten [39] suggested using several heating rates, a plot of ln(β/T2m) vs. (1/Tm) should yield a straight line of slope E/K, so that the activation energy is evaluated. Extrapolation to 1/Tm → 0 gives the value of ln(E/SK) from which the frequency factor is immediately found. It was shown that a plot of ln(Im) versus 1/Tmfor various heating rates usually yields a straight line too and the activation energy can be extracted similarly. It is to be noted from the theoretical point of view that β should be varied in as board a range as possible. However, this may cause various experimental difficulties. At very low heating rates, the maximum intensity will be low and in fact, the peak smeared, thus not allowing effective extraction of the experimental parameters. At high heating rates, a delay between the sample temperature and that of the measuring device impairs the temperature measurement. Moreover, temperature gradients within the sample usually occur at high heating rates which result in a smearing effect of a different kind. In practice, one should therefore compromise on a relatively narrow range of heating rates [10].

3.1.5 Three points method

A new technique was developed by Rasheedy [25], to evaluate the trap parameters from the measured glow curve according to the general kinetics order.

The behavior of a phosphor’s TL intensity is determined by the following equation, [40], for generic kinetics order.


Where Iis the intensity of the TL, n(cm−3), is the electron concentration trapped at time t(s), N(cm−3) is the traps concentration and K(eV/oK) is the Boltzman constant. Eq. (55) is more general than the two equations describing the first and second kinetics orders.

Eq. (55) is a modification of Eq. (33) in which the pre-exponential factor is defined as: S=SNb1instead of =Snob1.The solution of Eq. (55) is given by Rasheedy [40]:


Where the pre-exponential factor S″ = S(no/N)b−1 which is constant for a given dose but it varies with changes in the absorbed dose, i.e., with n0.

This method is based on the proportional of the concentration of populated traps during the running of the TL to the area under the glow peak.

Ixis the TL intensity at temperature Txat any portion of the glow peak as shown in Figure 10, then Eq. (55) becomes:

Figure 10.

Three points method used by Rasheedy [25] to investigate the equations used to calculate the trap parameters.


Where Axis the area under the glow peak between the temperatures Txand Tf(the final temperature of glow peak). Similarly, we have:


Where Iyand Izare the TL intensities at temperatures Tyand Tz, respectively.

From Eq. (57) and Eq. (58), we shall get


And from Eq. (57) and Eq. (59), we shall get:


The order of kinetics, b can be obtained using Eqs. (60) and (61) which leads to:


Then, the order of kinetics b can be obtained from Eq. (62). Once the order of kinetics b is determined, the activation energy E(eV) can be determined by using Eq. (60) or Eq. (61).

Since, at T=TmyieldsIT=0

From Eq. (56) and using Eq. (59) leads to the following expression [41]:


A simple analytical method has been developed to obtain the relative value of noin the case of general kinetics order [41]:


where Tm, and Imcan be obtained from the shape of the glow peak.

Thus, by calculating the kinetics order b, the activation energy E, and the initial trapped electrons number for many points that cover sufficient range on the glow peak, and taking the average value for each parameter, one can determine the trap parameters according to the general kinetics order.

3.1.6 Glow curve analysis (peak shape methods)

A review of the expression used in an intercomparison of glow curve analysis computer programs to evaluate TLD-100 glow curve is given in Ref. [42] where I(T) is written in the following form:


where: A = area (counts); b = kinetics order; E = activation energy; I = intensity (counts per s, counts per K); S = frequency factor (s−1).

On the other hand, Eq. (66) is based on first order kinetics which was used by Puchalska, [43], to develop glow-curve analysis software, in the following form:-


Where the parameter α is defined in Eq. (67) as


where the constants a0, a1 … and b0, b1 … are listed in the followings: -


Equation (68) will be used throughout our results which give better fitting to the resultant deconvoluted peaks. Different software was developed by Ratovonjanahary et al. [32], which uses the first kinetics order with an approximation of the second kinetics order. In this software the following equation was used:




Such a technique was also developed to analyze the glow curve using Eq. (53) by Rasheedy [41], which used the value of the trap parameters obtained by the three points method.


4. Modern clinical applications of TLD

TLD is widely used in various clinical fields for different purposes. The key reasons are undoubtedly their widespread availability, well-studied dosimetric characteristics, and applicability across a broad dose range. Imaging and Radiation Oncology Core-Houston IROC-H conducts remote dosimetry audits on MV photon and electron beams. IROC-H usually used integration between TLD-100 and other dosimetry system like nanoDot or diode systems for achieving the dose commissioning and calibrating dosimetry systems in an acrylic mini-phantom [44]. The failure rate was recorded in dose curves after modeling of the TPS (RayStation-Elekta Inc.) using phantom tests, which was not observed by patient-specific IMRT QA. Such failure was related to little changes in the MLC leaf-tip offset rather than leaf-tip width. Koger and his team [45] in IROC-H prosed four labeled TLD distributed in an anthropomorphic head-and-neck phantom for correcting such failure, (see Figure 11). It was utilized a 3D diode array were used in addition to assess the detectability of modeling mistakes [45].

Figure 11.

TLDs were labeled in head-and-neck phantom at IROC-H [45].

Another crucial issue is to increase the staff’s awareness about radiation safety and enhance radiation protection against unnecessary radiation doses. For such purpose, TLD-100 was recently used to validate occupational doses both inside and outside the nuclear medicine department, radiation protection purposes as well as the dose rate distribution around the positron emission tomography or computed tomography (PET/CT) [46].

Some recent studies were envisaged to see how the department compared to reports from other centers across the world in terms of the annual number of procedures and exposure limits, and to see if there was an opportunity for further radiation protection enhancements. As an example, personal TLD was calibrated to estimate the personal equivalent dose Hp (10) and Hp (0.07) at PET/CT. It was used for assessing the employee’s exposure [47]. On the other hand, TLD rings personal dosimeters were worn by surgeons in their fingers through sentinel node biopsy procedure to measure personal doses Hp(10) and Hp(0.07), as well as ambiental dose for operating theater and during injection [48]. This will assure that personal equivalent doses are within the acceptable annual determined limits [49].

Other important recent TLD application in diagnostics is using an anthropomorphic phantom that modeled the reference person to get a conversion coefficient connecting dose area product (DAP) to effective patient dosage. They concluded that the effective dosage at the clinical dark-field radiography system, which generates both attenuation and dark-field pictures, is within the range of chest radiography standard dose values [50].

TLDs showed to be an excellent choice for skin dosimetry. Omojola et al. [51] used TLD in measurements of 3D skin dosimetry and verify their results using TPS planning verification at specific spots in the phantom. A full perspective of the dose distribution was achieved; however, they revealed that regions outside the PTV require special attention [52, 53, 54].

In addition, in the field of proton therapy, a novel tissue-equivalent TLD-sheet of manganese doped lithium triborate showed a valuable and effective dosimetry technique. It may also be a great in vivo skin dosimetry instrument for proton treatment due to its flexible and reusable properties. Despite the presence of significant energy dependences in the Bragg peak region, the response properties studied in this work, including as reproducibility, fading effects, dosage linearity and dose homogeneity are acceptable [55].

Monte Carlo (MC) simulation is considered a good tool to understand well the TLD [56]. Some algorithm methods based on MC as if pencil beam could be involved in accurate dose in MV radiotherapy calculations. It could be useful to calculate the spectrum inside the detector based on four categories primary photon and electrons and secondary photon and electrons [57].

On the other hand, Low-energy (100 keV) photons (x-rays and gamma) have been widely employed in biological research and medical applications for more than a century, including mammography, fluoroscopy, general radiography, computed tomography, and brachytherapy treatment, among others. The majority of electrons created by low photon energy beams have energies below 10 keV, according to research. The physical processes through which these low-energy electrons interact with matter, on the other hand, are still unknown. Furthermore, it is commonly thought that all energy put within a dosimeter-sensitive volume is converted into a response. However, this assumption could be inaccurate because some of the deposited energy could be utilized to build flaws or damages at the molecular and atomic levels [58].

The hybrid-functional density theory (H-DFT) has shown to be a promising tool for localizing secondary electrons within a dosimeter volume and calculating the energy spent on creating defects or colors centers, among other things, when it comes to the relationship between the energy deposited and the response of a dosimeter. Following that, the quantity of energy that can be truly turned into a dosimeter response following exposure to ionizing radiation would be more accurately determined.


5. Conclusion

The aim of this chapter is concerned with TLD materials, measurements and recent various applications in clinical and industrial fields. TL kinetics are also covered in details due to their importance in knowing traps parameters and band structure-related phenomena that are responsible for TL phenomena. Modern clinical applications of TLD are also covered like quality assurance purposes for proton, x-ray and gamma radiotherapy based on phantom tests. In addition, we shed spot on using TLD for recent accurate methods for skin dose evaluation under IMRT/VMAT radiotherapy. Special attention should be oriented to hybrid-functional density theory Monte Carlo simulation to model TL dosimeters. Recent studies proved a promising tool for localizing secondary electrons within a dosimeter volume and calculating the energy spent on creating defects or colors centers, among other things, when it comes to the relationship between the energy deposited and the response of a dosimeter. Such methods could give knowledge about misunderstanding behaviors of some TLD and could eliminate its disadvantages like missing TL signal or fading; angle and energy of incidence ionizing radiation. In general, the properties of TLD like its inexpensive cost and reusability; easily fabricated, lightweight to wear, readout is simple and can be automated, insensitive to humidity make it advantageous in different clinical and radiation safety applications.


  1. 1. Furuta Y, Tanaka S. Response of 6LiF and 7LiF thermoluminescence dosimeters to fast neutrons. Nuclear Instruments and Methods. 1972;104:36-374
  2. 2. Moreira Ribeiro R, Souza-Santos D. Monte Carlo characterization of an individual albedo neutron monitor. Brazilian Journal of Radiation Sciences. 2021;9(2C):1-8. DOI: 10.15392/bjrs.v9i2C.1657
  3. 3. El-Faramawy N, Chopra V, Rawash S, El-Hafez AA, Dhoble SJ. Response of TLD-600/TLD-700 and CR-39 to neutrons for medical dosimetry. Luminescence. 2021;36(5):1257-1264
  4. 4. Gibson AB, Piesch E. Technical Reports Series No. 252. Neutron Monitoring for Radiological Protection. Vienna: International Atomic Energy Agency. 1985
  5. 5. ICRU. Determination of Dose Equivalents Resulting from External Radiation Sources. In: Report 39, International Commission on Radiation Protection and Measurements. Bethesda, MD. 1985
  6. 6. Seifert H, Dörschel B, Pawelke J, Hahn T. Comparison of calculated and measured neutron sensitivities of an electret albedo dosemeter. Radiation Protection Dosimetry. 1991;37(1):13-18
  7. 7. Seifert H, Dorschel B, Pawelke J. Neutron dosimetry using an electret albedo dosimeter. In: 1991 Proceedings 7th International Symposium on Electrets (ISE 7). London: IEEE; 1991. pp. 753–758
  8. 8. Fellinger J, Hahn T, Henniger J, Hübner K, Schmidt P. Fast neutron sensitivity of TL detectors using proton radiator techniques. Isotopenpraxis Isotopes in Environmental and Health Studies. 1991;27(7):342-346. DOI: 10.1080/10256019108622562
  9. 9. El-Adawy A, Khaled NE, El-Sersy AR, Hussein A, Donya H. TL dosimetric properties of Li2O–B2O3 glasses for gamma dosimetry. Applied Radiation and Isotopes. 2010;68(6):1132-1136
  10. 10. Polymeris GS, Çoskun S, Tsoutsoumanos E, Konstantinidis P, Aşlar E, Şahiner E, et al. Dose response features of quenched and reconstructed, TL and deconvolved OSL signals in BeO. Results in Physics. 2021;25:104222
  11. 11. Begum M, Rahman AM, Abdul-Rashid HA, Yusoff Z, Nawi SN, Khandaker MU, et al. Photonic crystal fibre as a potential medium for radiotherapy dosimetry. Applied Radiation and Isotopes. 2021;174:109771
  12. 12. Saray AA, Kaviani P, Shahbazi-Gahrouei D. Dosimetric characteristics of lithium triborate (LiB3O5) nanophosphor for medical applications. Radiation Measurements. 2021;140:106502
  13. 13. Horowitz YS, Yossian D. Computerised glow curve deconvolution: Application to thermoluminescence dosimetry. Radiation Protection Dosimetry. 1995;60(1):1-114
  14. 14. Randall JT, Wilkins MHF. Phosphorescence and electron traps II. The interpretation of long-period phosphorescence. Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences. 1945;184(999):390-407
  15. 15. Garlick GFJ, Wilkins MHF. Short period phosphorescence and electron traps. Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences. 1945;184(999):408-433
  16. 16. Chen R. On the calculation of activation energies and frequency factors from glow curves. Journal of Applied Physics. 1969;40(2):570-585
  17. 17. Wrzesińska A. Their production and thermoluminescence curves. Acta Physica Polonica. 1956;15:151
  18. 18. Adirovitch EI. La formule de Becquerel et la loi élémentaire du déclin de la luminescence des phosphores cristallins. Journal de Physique et le Radium. 1956;17(8-9):705-707
  19. 19. Halperin A, Braner AA, Ben-Zvi A, Kristianpoller N. Thermal activation energies in NaCl and KCl crystals. Physical Review. 1960;117(2):416
  20. 20. Dussel GA, Bube RH. Theory of thermally stimulated conductivity in a previously photoexcited crystal. Physical Review. 1967;155(3):764
  21. 21. Saunders IJ. The thermally stimulated luminescence and conductivity of insulators. Journal of Physics C: Solid State Physics. 1969;2(12):2181
  22. 22. De Muer D. Development of a universal method for calculating the thermoluminescence parameters. Physica. 1970;48(1):1-12
  23. 23. Bräunlich P, Kelly P, Fillard JP. Thermally stimulated luminescence and conductivity. In: Thermally Stimulated Relaxation in Solids. Berlin, Heidelberg: Springer; 1979. pp. 35-92
  24. 24. Moharil SV. On the general‐order kinetics in thermoluminescence. Physica Status Solidi A: Applications and Materials Science. 1982;73(2):509-514
  25. 25. Rasheedy MS. A new evaluation technique for analyzing the thermoluminescence glow curve and calculating the trap parameters. Thermochimica Acta. 2005;429(2):143-147
  26. 26. May CE, Partridge JA. Thermoluminescent kinetics of alpha‐irradiated alkali halides. The Journal of Chemical Physics. 1964;40(5):1401-1409
  27. 27. Prokić M. Analysis of the thermoluminescence glow curves of natural barite. Journal of Physics and Chemistry of Solids. 1977;38(6):617-622
  28. 28. De Blasi C, Gallassini S, Manfredotti C, Micocci G, Ruggiero L, Tepore A. Trapping levels in PbI2. Solid State Communications. 1978;25(3):149-153
  29. 29. Urbach R. Zur lumineszenz der alkalihalogenide. Sitzungsberichte Akad. der Wiss. Wien. 1930;139:363-372
  30. 30. Urbach F. Storage and Release of Light by Phosphors. Vol. 115. New York: John Wiley and Sons; 1948
  31. 31. Miller LD, Bube RH. Luminescence, trapping, and F centers in lithium fluoride crystals. Journal of Applied Physics. 1970;41(9):3687-3697
  32. 32. Ratovonjanahary AJF, Raboanary R, Andriambololona R. Quartz Glow-Peaks Lifetime Analysis: TL Glow-Curve Deconvolution Functions for First Order of Kinetic Compared to Initial Rise Method. In: HEPMAD 04 Conference, Madagascar, 27 Sep-01 Oct. 2004
  33. 33. Christodoulides C. Errors involved in the determination of activation energies in TL and TSDC by the initial rise method. Journal of Physics D: Applied Physics. 1985;18(8):1665
  34. 34. Grossweiner LI. A note on the analysis of first‐order glow curves. Journal of Applied Physics. 1953;24(10):1306-1307
  35. 35. Lushchik CB. The investigation of trapping centers in crystals by the method of thermal bleaching. Soviet Physics Jetp-USSR. 1956;3(3):390-399
  36. 36. Booth AH. Calculation of electron trap depths from thermoluminescence maxima. Canadian Journal of Chemistry. 1954;32(2):214-215
  37. 37. Chen R, Winer SAA. Effects of various heating rates on glow curves. Journal of Applied Physics. 1970;41(13):5227-5232
  38. 38. Chen R, Kirsh Y. Analysis of Thermally Stimulated Process. Oxford: Pergamon Press; 1981
  39. 39. Hoogenstraaten W. Electron traps in zinc sulphide phosphors. Philips Research Reports. 1958;13:515-693
  40. 40. Rasheedy MS. On the general-order kinetics of the thermoluminescence glow peak. Journal of Physics: Condensed Matter. 1993;5(5):633
  41. 41. Rasheedy MSRMS. A complete system for obtaining the trap parameters of thermoluminescence glow peak. Japanese Journal of Applied Physics. 1996;35(2R):634
  42. 42. Bos AJJ, Piters TM, Gómez-Ros JM, Delgado A. An intercomparison of glow curve analysis computer programs: I. Synthetic glow curves. Radiation Protection Dosimetry. 1993;47(1-4):473-477
  43. 43. Puchalska M, Bilski P. GlowFit—A new tool for thermoluminescence glow-curve deconvolution. Radiation Measurements. 2006;41(6):659-664
  44. 44. Alvarez P, Kry SF, Stingo F, Followill D. TLD and OSLD dosimetry systems for remote audits of radiotherapy external beam calibration. Radiation Measurements. 2017;106:412-415
  45. 45. Koger B, Price R, Wang D, Toomeh D, Geneser S, Ford E. Impact of the MLC leaf-tip model in a commercial TPS: Dose calculation limitations and IROC-H phantom failures. Journal of Applied Clinical Medical Physics. 2020;21(2):82-88
  46. 46. Nilsson I, Himmelman J, Khan J, Dalmo J. The potential to use Tld measurements to validate the occupational radiation protection at the Department of Nuclear Medicine. Radiation Protection Dosimetry. 2021;195(3-4):355-362
  47. 47. Pavičar B, Davidović J, Petrović B, Vuleta G, Trivić S, Šajinović V, et al. Nuclear medicine staff exposure to ionising radiation in 18F-FDG PET/CT practice: A preliminary retrospective study. Arhiv za Higijenu Rada i Toksikologiju. 2021;72(3):216-223
  48. 48. Petrovic B, Vicko F, Radovanovic D, Samac J, Tot A, Radovanovic Z, et al. Occupational radiation dose of personnel involved in sentinel node biopsy procedure. Journal of Medical Physics. 2021;91:117-120
  49. 49. Ali W, Sulieman A, Tamam N, Boshara N, Aldhebaib A, Alkhorayef M, et al. Estimation of patients organ doses and staff exposure during bone scan examination. Radiation Physics and Chemistry. 2021;188:109693
  50. 50. Frank M, Urban T, Willer K, Noichl W, De Marco F, Schick R, et al. Dosimetry on first clinical dark-field chest radiography. Medical Physics. 2021;48(10):6152-6159
  51. 51. Omojola AD, Akpochafor MO, Adeneye SO, Akala IO, Agboje AA. Chest X-rays of newborns in a medical facility: Variation between the entrance skin dose measurements using the indirect and direct methods for clinical dose audit. Japanese Journal of Radiology. 2021;40(2):1-7
  52. 52. Moradi F, Khandaker MU, Mahdiraji GA, Ung NM, Bradley DA. Dose mapping inside a gamma irradiator measured with doped silica fibre dosimetry and Monte Carlo simulation. Radiation Physics and Chemistry. 2017;140:107-111
  53. 53. Moradi F, Ung NM, Mahdiraji GA, Khandaker MU, Entezam A, See MH, et al. Angular dependence of optical fibre thermoluminescent dosimeters irradiated using kilo- and megavoltage X-rays. Radiation Physics and Chemistry. 2017;135:4-10
  54. 54. Moradi F, Mahdiraji GA, Dermosesian E, Khandaker MU, Ung NM, Mahamd Adikan FR, et al. Influence of dose history on thermoluminescence response of Ge-doped silica optical fibre dosimeters. Radiation Physics and Chemistry. 2017;134:62-70
  55. 55. Kato T, Sagara T, Komori S, Kato R, Takeuchi A, Narita Y. Dosimetric properties of a newly developed thermoluminescent sheet-type dosimeter for clinical proton beams. Journal of Applied Clinical Medical Physics. 2021;22(4):158-165
  56. 56. Donya H, Seniwal B, Darwesh R, Fonseca TC. Prospective Monte Carlo simulation for choosing high efficient detectors for small-field dosimetry. In: Theory, Application, and Implementation of Monte Carlo Method in Science and Technology. London: IntechOpen; 2019
  57. 57. Donya H. Pencil-beam fluence evaluation based on Monte Carlo simulations algorithm of high energetic treatment photons. Journal of Medical Signals and Sensors. 2018;8(2):81
  58. 58. Massillon-JL G. Future directions on low-energy radiation dosimetry. Scientific Reports. 2021;11:10569. DOI: 10.1038/s41598-021-90152-3

Written By

Hossam Donya

Submitted: December 29th, 2021 Reviewed: January 17th, 2022 Published: March 25th, 2022