^{1}

^{*}

^{2}

The aim of this paper is to give some simplified expressions related to the peak shape method. The modified equations have been used to calculate the activation energy (
*E*) of commercial thermoluminescent dosimeters (TLD), as well as of ZnO thermoluminescent material produced in laboratory; the values so determined have been compared to the values obtained using the classical expressions of the peak shape method. The modified equations proposed are as a function of peak shape parameters or the peak temperature at the maximum. This expression could be useful to obtain approximated E values in the case of complex glow curves as well, when the peaks are not well resolved but the peak temperature at the maximum may be easily determined.

Luminescence is defined as the emission of light from some solids commonly called phosphors. This emission is the release of energy stored within the solid through some type of prior excitation of the solid electronic system, i.e. by visible, infrared (IR) or ultra violet (UV) light and ionizing radiation [

The luminescence effect can be used in solid state dosimetry for measurement of ionizing radiation dose; the most popular luminescence dosimetric technique is called thermoluminescence (TL): it consists in a transient emission of light from an irradiated solid when heated [

Thermoluminescence is characterized by the following main points:

1) Transient effects, superimposed on an intrinsic background;

2) A memory but not a permanent record of the irradiation dose.

Materials of interest in thermoluminescent dosimetry (TLD) are principally insulators in which the conduction electrons are entirely due to the absorbed radiation energy.

Thermoluminescence is observed under condition of steadily increasing temperature. In the usual thermoluminescence experiments, the TL material is irradiated at room temperature (RT) and later heated through a temperature range where the luminescence is bright, until a temperature level at which all the charges have been thermally excited out of their metastable levels and the luminescence completely disappears [

If the light intensity is plotted as a function of temperature (or time), the resulting graph is called glow-curve (

In the field of TL dosimetry, the analysis of multipeak glow curves includes the evaluation of the kinetics parameters E (activation energy) and s (frequency factor) related to each TL peak. Several methods to determine these kinetics parameters, mainly the activation energy, have been proposed, based on one-trap models (related to a single peak), and hence the applicability of these methods to multi-peak glow curves is limited for the complexity of the glow curves shape [

parameters of a glow curve [

The aim of this paper is to propose some simplified expressions related to the peak shape method. The modified equations have been used to calculate the activation energy, when the peaks are not well resolved but the peak temperature at the maximum may be easily determined.

This paper presents simulated glow peaks obtained using the computation of the first order kinetics equation, i.e. the Randall-Wilkins equation [

I ( T ) = n 0 s exp ( − E k T ) exp [ − s β ∫ T 0 T exp ( − E k T ′ ) d T ′ ] (1)

where, n_{0} is the initial concentration of filled traps; s is a pre-exponential frequency factor; E is the activation energy or trap depth; k is Boltzmann’s constant; T is the absolute temperature and β is the heating rate.

For the second order equation, i.e. the Garlick-Gibson equation [

I ( T ) = n 0 s exp ( − E k T ) [ 1 + s β ∫ T 0 T exp ( − E k T ′ ) d T ′ ] 2 (2)

and, finally, the general order equation, i.e. the May-Partridge equation [

I ( T ) = s n 0 exp ( − E k T ) [ 1 + s ( b − 1 ) β ∫ T 0 T exp ( − E k T ′ ) d T ′ ] − b b − 1 (3)

where, b the kinetic order, a parameter with values typically between 1 and 2.

The shape of the glow peak is different according to the kinetic process, as is observed in

An isolated TL glow peak, obtained using a linear heating rate, without the kinetic order it follows, can be characterized by some geometrical parameters T_{M}, T_{1}, T_{2}, τ, δ and θ, as it can be seen in _{M}, T_{1}, T_{2}: are respectively the peak temperature at the maximum and the temperatures on either side of the temperature at the maximum, corresponding to half intensity. The geometrical parameters from a peak are evaluated as follows:

- τ = T M − T 1 : is the half-width at the low temperature side of the peak,

- δ = T 2 − T M : is the half-width towards the fall-off of the glow peak, (4)

- ω = T 2 − T 1 : is the total half-width,

- μ g = δ / ω : is the so called symmetrical geometrical factor.

Considering the so given characteristic temperatures of a glow peak, glow peaks of both first and second order kinetics have been generated by computer using Equations (1) and (2) and the following relations among the peak temperatures have been obtained:

- In case of first and second order peaks: T 1 = 0.95029 T M

- For a first order peak: T 2 = 1.03132 T M (5)

- For a second order peak: T 2 = 1.05648 T M

It has also to be noted that:

- According to the asymmetrical property on a glow curve shape, for a first-order peak τ is almost 40% bigger than δ, due the re-trapping process involved.

- The geometrical factor µ_{g} is equal to 0.42 for first-order kinetics, and 0.52 in the case of second-order.

- µ_{g} is practically independent of E, in the range from 0.1 to 1.6 eV.

- µ_{g} is strongly dependent on the kinetics order, b, in the range 0.7 £ b £ 2.5,

- Another factor, namely γ = δ/τ, is ranging from 0.7 to 0.8 for a first-order peak, and from 1.05 to 1.20 for a second-order.

An important and widely used method to investigate the trapping levels in crystals is based on the geometrical characterization of a TL glow peak, the well-known peak shape methods. In fact, for calculating the activation energy of the trapping level corresponding to a peak in the glow curve, one needs to measure three temperature values on the peak itself: T_{M}, T_{1}andT_{2}. The equations for determining the activation energy according to the peak shape methods are the followings:

Grossweiner was the first to use the shape of the glow curve to calculate the trap depth and frequency factor [_{1}. The expressions to evaluate E, as revised by Chen [

1^{st} order ( E G ) τ = 1.41 k T 1 T M τ (6)

2^{nd} order ( E G ) τ = 1.68 k T 1 T M τ (7)

Lushchik also proposed a method based on the shape for first- and second-order kinetics introducing the parameter δ [

1^{st} order ( E L ) δ = 0.976 k T M 2 δ (8)

2^{nd} order ( E L ) δ = 1.706 k T M 2 δ (9)

where the coefficients have been computed by Chen [

Halperin & Braner gave a different approach by using both T_{1} and T_{2}. They consider the luminescence emission as mainly due to different kinds of recombination processes [

1^{st} order ( E H B ) τ = 1.72 k T M 2 τ ( 1 − 1.58 Δ M ) (10)

2^{nd} order ( E H B ) τ = 2 k T M 2 τ ( 1 − 3 Δ M ) (11)

where Δ M = 2 k T M E .

Both the previous Equations (9) and (10) require an iterative process owing the presence of the quantity D_{M}. To avoid this difficulty, Chen modified them as following [

1^{st} order ( E H B ) τ = 1.52 k T M 2 τ − 3.16 k T M (12)

2^{nd} order ( E H B ) τ = 1.813 k T M 2 τ − 4 k T M (13)

Chen also gave two expressions, using numerical approximations and without knowledge of the kinetics order which is directly found from the peak shape. Based on the Lushchik, Halperin and Braner previous studies, he proposed the following equations based on parameter ω [

1^{st} order ( E H B ) ω = 2 k T M ( 1.25 T M ω − 1 ) (14)

2^{nd}order ( E H B ) ω = 2 k T M ( 1.76 T M ω − 1 ) (15)

Balarin also gave expressions based on the quantity ω, [

1^{st} order ( E B ) ω = T M 2 ω ⋅ 4998 (16)

2^{nd} order ( E B ) ω = T M 2 ω ⋅ 3542 (17)

Considering the three temperature values characterizing a TL peak, i.e. T_{M}, T_{1} and T_{2}, and the asymmetrical and symmetrical characteristics of the glow peaks following first or second order kinetics respectively, synthetic peaks have been generated, giving the temperature at the maximum, T_{M}, ranging from 300 K to 750 K.

In this way, the original equations used in the peak shape method, from Equation (6) to Equation (17), have been modified using a little algebra and the following new equations have been obtained as a function of the peak temperatures already defined trough Equations (4) and (5).

The modified equations were evaluated as a function of the peak temperature parameters and of the peak temperature at the maximum, as follow:

Also, the modified equations were evaluated as a function of the peak temperature at the maximum, as follow:

The following _{2}: Tm (TLD-300) [

With the aim to apply the modified equations on a glow curve having peaks not well resolved, the values of E of a ZnO thermoluminescent material have been calculated, and then compared to the values obtained using the original peak shape expression, reported on a previous work [

In

Phosphor | Reference | E (eV) 1st peak | E (eV) 2nd peak | E (eV) 3rd peak | E (eV) 5th peak |
---|---|---|---|---|---|

CaF2:Tm (TLD-300) | [ | 0.69 ÷ 0.74 | 1.07 ÷ 1.31 | 1.68 ÷ 1.90 | |

This work | 0.91 ÷ 1.04 | 0.91 ÷ 1.18 | 1.35 ÷ 1.45 | ||

BeO | [ | 0.99 ÷ 1.55 | 1.16 ÷ 1.35 | ||

This work | 1.17 ÷ 1.38 | 1.49 ÷ 1.72 |

Method | Kinetic Order | Activation energy, E (eV) 1st Peak | Activation energy, E (eV) 2nd Peak | Activation energy, E (eV) 3rd Peak | |
---|---|---|---|---|---|

Grossweiner | 1st Order | 0.54 | 0.74 | 0.99 | |

2nd Order | 0.65 | 0.85 | 1.18 | ||

Lushchik | 1st Order | -- | -- | 1.07 | |

2nd Order | -- | -- | 1.88 | ||

Halperin- Braner | 1st Order | 0.55 | 0.77 | 1.03 | |

2nd Order | 0.62 | 0.88 | 1.17 | ||

Balarin | 1st Order | -- | -- | 0.96 | |

2nd Order | -- | -- | 1.04 | ||

experimental glow curve (

As it can be observed from the Tables, the E values calculated using the classical methods, i.e. Equations (6) to (17) are not much different from the ones obtained using the modified equations here proposed.

For the expressions as a function of the peak temperature parameters, the E values for low temperature peaks are higher than the values obtained from the

Method | Kinetic Order | Activation energy, E (eV) 1st Peak | Standard Deviation | Activation energy, E (eV) 2nd Peak | Standard Deviation | Activation energy, E (eV) 3rd Peak | Standard Deviation | |
---|---|---|---|---|---|---|---|---|

Grossweiner | 1st Order | 1.48 | 0.66 | 1.84 | 0.77 | 1.75 | 0.53 | |

2nd Order | 1.76 | 0.78 | 2.2 | 0.95 | 2.09 | 0.64 | ||

Lushchik | 1st Order | -- | -- | -- | -- | 1.68 | 0.43 | |

2nd Order | -- | -- | -- | -- | 1.31 | 0.40 | ||

Halperin-Braner | 1st Order | 2.34 | 1.26 | 2.92 | 1.52 | 2.78 | 1.23 | |

2nd Order | 1.47 | 0.60 | 1.84 | 0.67 | 1.74 | 0.40 | ||

Balarin | 1st Order | -- | -- | -- | -- | 1.73 | 0.54 | |

2nd Order | -- | -- | -- | -- | 1.67 | 0.44 | ||

Method | Kinetic Order | Activation energy, E (eV) 1st Peak | Standard Deviation | Activation energy, E (eV) 2nd Peak | Standard Deviation | Activation energy, E (eV) 3rd Peak | Standard Deviation |
---|---|---|---|---|---|---|---|

Grossweiner | 1st Order | 0.92 | 0.26 | 1.19 | 0.31 | 1.34 | 0.24 |

2nd Order | 1.09 | 0.31 | 1.42 | 0.40 | 1.60 | 0.29 | |

Lushchik | 1st Order | 1.06 | -- | 1.38 | -- | 1.55 | 0.33 |

2nd Order | 1.03 | -- | 1.34 | -- | 1.50 | 0.26 | |

Halperin- Braner | 1st Order | 0.93 | 0.26 | 1.21 | 0.31 | 1.36 | 0.23 |

2nd Order | 1.10 | 0.33 | 1.43 | 0.38 | 1.61 | 0.31 | |

Balarin | 1st Order | 1.02 | -- | 1.32 | -- | 1.48 | 0.36 |

2nd Order | 1.02 | -- | 1.33 | -- | 1.49 | 0.31 |

original equations. As the T_{M} value increases the standard deviation decreases and, furthermore, the E values are more closed to the values obtained before.

In any case, the expressions proposed could be useful to obtain approximated E values in the case of complex glow curves, when the peaks are not well resolved but the peak temperature at the maximum may be easily determined; furthermore, the approximated E values so determined may be inserted in a deconvolution programmer which allows a more precise determination of the kinetics parameters.

Borbón-Nuñez, H.A. and Furetta, C. (2017) Activation Energy of Modified Peak Shape Equations. World Journal of Nuclear Science and Technology, 7, 274-283. https://doi.org/10.4236/wjnst.2017.74021