Use of Terbium Doped Phosphate Glasses for High Dose Radiation Dosimetry—Thermoluminescence Characteristics, Dose Response and Optimization of Readout Method

The phosphate glass samples doped with Tb2O3 oxide (general formula: P2O5-Al2O3-Na2O-Tb2O3) were synthesized and studied for usage in high-dose radiation dosimetry (for example, in high-activity nuclear waste disposals). The influence of terbium concentration on thermoluminescent (TL) signals was analyzed. TL properties of glasses were investigated using various experimental techniques such as direct measurements of TL response vs. radiation dose, Tmax–Tstop and VHR (various heating rate) methods, and glow curve deconvolution analysis. The thermoluminescence dosimetry (TLD) technique was used as the main investigation tool to study detectors’ dose responses. It has been proved that increasing the concentration of terbium oxide in glass matrices significantly increases the thermoluminescence yield of examined material. For the highest dose range (up to 35 kGy), the dependence of the integrated thermoluminescent signals vs. dose can be considered as a saturation-type curve. Additional preheating of samples improves linearity of signal vs. dose dependencies and leads to a decrease of the signal loss over time. All obtained data suggest that investigated material can be used in high-dose radiation dosimetry. Additional advantages of the investigated dosimetric system are its potential ability to re-use the same dosimeters multiple times and the fact that reading dosimeters only requires usage of a basic TL reader without any modifications.


Introduction
Fast progress and rising interest in technologies involving interactions of ionizing radiation with matter force the continuous development of dosimetric techniques. The number of dosimetric systems present nowadays on the market is limited, especially in the narrow field of high-dose radiation dosimetry (especially in dose ranges up to 35 kGy). Many of these systems are based on organic compounds and polymers (alanine, poly-methylmethacrylate, etc.) and have numerous limitations in terms of use in harsh environments. One of the potential alternatives that may turn out to be complementary to existing systems is phosphate glasses doped with terbium ions.
Terbium is a lanthanide, belonging to the rare-earth elements (REE) group. REE include lanthanum and the f-block elements, cerium through lutetium. Next to the above, scandium and yttrium are included in this group as well, due to the similar ionic radii with the lighter f-block elements. Moreover, they co-exist together in the same types of ores [1]. Lanthanides exhibit an important feature, namely, they can play the role of light-emitting activators in various matrices. The rare earth phosphors can be divided into two groups: the broad-band-emitting, owing to the 5d→4f transition (Eu 2+, Ce 3+ ), and the narrow-bandemitting, owing to the transition between the 4f levels like Tb 3+ . Terbium ion emits green light due to the four sharp peaks arising from 5 D 4 → 7 F J transitions, where J = 3, 4, or 6, with the strongest emission peak at 550 nm [2,3]. Terbium possesses only one stable isotope Tb-159, characterized by a low nuclear spin value (I = 3/2), and exhibits a relatively simple hyperfine structure of its electronic levels [4].
Terbium ions (Tb 3+ ) are very often used as dopants in order to extend luminescence properties in various materials, e.g., in nitridosilicates, glasses, and ceramics [5][6][7][8]. Terbium ions are also able to play the role of activators in mechano-luminescence light-emission processes, often as inorganic salts such as sulphates, carbonates, nitrates, and chlorides [9,10]. Terbium dopants have been also used in various types of thermoluminescent (TL) materials. For example, magnesium orthosilicate doped with terbium (Mg 2 SiO 4 :Tb) has been investigated in TL dosimetry since 1971, thanks to its high sensitivity and effective atomic number of 11 [11,12]. Unfortunately, TL properties of orthosilicates might be unstable, which are strongly influenced by the preparation method [13]. Another group of lightemitting materials that are based on terbium dopant, improving their TL properties, are ionic crystals, e.g., CsCl:Tb 3+ , terbium-doped strontium pyrophosphates (Sr 2 P 2 O 7 :Tb 3+ ), potassium bromide crystals, and terbium-doped double fluorides [14][15][16][17]. Moreover, the application of terbium oxide as a dopant significantly improved the properties of materials that had been already used in TL dosimetry. A straightforward example of such materials is aluminum oxide doped with terbium and thalium [18]. Lu 2 O 3 :Tb,Ta ceramic storage phosphors are also a good example of how the presence of terbium ions can significantly change the thermoluminescence properties of materials; Tb,Ta ceramics show high overall TL intensity [19]. Despite all the advantages of the crystal-based phosphors, they possess drawbacks as well. For TL dosimeters based on crystalline systems, one can observe an increased light-scattering efficiency, due to their limited transparency. A reasonably good alternative for these phosphors might be found among dosimetric systems based on glasses, for example, borosilicate glasses [20]. One such example is zinc-lithium-borate glass doped with terbium. Introduction of terbium into the material results in an increase of TL yield, enhancement of sensitivity of the phosphor, and significant minimization of fading behavior [21]. BaO-P 2 O 5 glasses are also an interesting example of influence of terbium ions on matrices. The analysis of the TL data of BaO-M 2 O 3 -P 2 O 5 :Tb 2 O 3 glasses suggests that, comparatively, there are high non-radiative losses in In 2 O 3 mixed glasses [22]. Another example of the improvement of dosimetric characteristics by the introduction of terbium into glass hosts is terbium-doped CaO-Al 2 O 3 -B 2 O 3 glasses. As in the previous case, dosimeters were found to exhibit good sensitivity and linearity. Sadly, for these materials one can observe more than a 60% decrease in signal intensity in the first 72 h after irradiation [23].
Phosphate glasses are currently used in various fields of science, technology, and medicine, such as in matrices for high-level radioactive waste (HLW) immobilization, in fast ion conductors, and in bone transplantation [24][25][26]. Phosphate glasses exhibit lower transition temperatures, lower viscosity when in a molten state, higher thermal expansion coefficients, and higher electrical conductivity than borate and silicate glasses [26]. As shown before, phosphate glasses also possess interesting TL properties suitable for highdose ionizing radiation dosimetry [27]. Glasses doped with Tb 3+ are also used in other luminescent materials, for example, lasers. [28]. An obvious advantage of phosphate glass systems is that their synthesis requires significantly lower temperatures and less energy, which makes the process of their preparation economically competitive in comparison with other glasses. In the past, there were some papers dealing with the thermoluminescence of phosphate glasses doped with rare earth elements such as terbium, but none of them investigated in detail thermoluminescence kinetics and the usage of this material in highdose dosimetry [29,30].
The presented work reports the influence of terbium doping on dosimetric properties of phosphate glass by means of thermoluminescence techniques. Here, one investigated the influence of terbium concentration on TL signal characteristics and the dependencies of TL signal intensity vs. dose absorbed in a material. Next to the above, the nature of thermoluminescence signals was investigated using various methods (such as the T max -T stop or VHR methods). Moreover, the paper also proposes the readout method optimal for investigated dosimeters and investigates its influence on TL signals, linearity of signals intensity vs. radiation dose absorbed within the material, and TL signal loss over time.

Glass Samples Synthesis
The phosphate glass samples were synthesized according to the melt-quench technique from raw materials, such as Na 2 CO 3 , NH 4 H 2 PO 4 , Al 2 O 3, and Tb 4 O 7 of highest available purity. The details on the synthesis procedure were discussed in previous papers [27]. Four different glass compositions were prepared for this research. The chemical compositions of each investigated glass system synthesized differed only due to terbium oxide concentration, whereas the main components of the glasses remained practically unchanged. In order to follow the overall influence of the terbium dopant on TL properties of the investigated glasses, for comparison raw (undoped) glass sample was prepared simultaneously (PPG) and investigated in the same manner as terbium-doped glass samples. The detailed compositions of all investigated phosphate glass samples are presented in Table 1. The as-synthesized glass was grounded in a planetary ball mill (Retsch PM100, Haan, Germany) using zirconium oxide (ZrO 2 ) grinding vessel and balls. The obtained powders were classified using laboratory sieves, to remove remaining coarse fraction.

XRD and DSC Analysis of Phosphate Glasses
To confirm the glassy state of the phosphate glass samples XRD (X-Ray Diffraction), analysis was performed using powdered samples. The X-ray diffractometer used was X'Pert Pro MPD (PANanalytical, Almelo, The Netherlands), and measurements were performed within 2theta angle between 5-60 deg. Data acquisition was carried out using Cu Kα radiation and scan rate of 6 deg/min, whereas X-ray tube was operating at 40 kV and 30 mA. Next to XRD analysis, one also performed differential scanning calorimetry analysis of the native and terbium-doped phosphate glasses within 40-550 • C temperature range using Q200 DSC (TA Instruments, New Castle, DE, USA) with ca. 15 mg samples in aluminum crucibles applying a heating rate of 20 deg/min.

Preparation of Glass Pellet Dosimeters
To enhance properties of detectors and to facilitate manual operations, powders of phosphate glass were pressed into pellets. The weighted portions of glass powders (40 mg) were placed in the cylindrical pressing die (hole diameter: 5 mm). Materials were compressed using a hydraulic press, then pellets were sintered in a laboratory furnace (LE6/11/P300, Nabertherm, Lilienthal, Germany) for 1 h at 420 • C. The sintering step was applied to improve the mechanical properties of the detectors and to erase any residual TL signals from the detectors.

Samples Irradiation
The investigated samples were irradiated using γ and β radiation installed in Institute of Applied Radiation Chemistry, Lodz University of Technology. Gamma irradiation was performed using Co-60 radiation source installed in Ob-Servo-D panoramic irradiator Appl. Sci. 2021, 11, 7221 4 of 21 (Izotop, Budapest, Hungary), whereas βirradiation was performed using electron beam from linear electron accelerator (ELU-6E), which was able to produce electrons of energy of 6 MeV. The e-beam irradiation was performed in a single-pulse mode (17 ns or 4 µs duration), depending on the required dose. The range of cumulative doses applied to samples using 17 ns pulses was 0.06-1.0 kGy and 0.33-35 kGy for 4 µs pulses. As a reference, for verification of radiation dose delivered to the samples, e-Scan alanine dosimetry system (Bruker, Karlsruhe, Germany) was used.

Thermoluminescence (TL) Measurements
TLD reader (Mikrolab TL Ra94, Poland) was used for thermoluminescence emission measurements. The measurements were carried out in the analyzer mode, within the temperature range of 40-400 • C using heating rate of 5 deg/s. Due to excessive light emission from samples containing terbium, it was necessary to use a shutter with hole (diameter of 0.5 mm) to limit light flux reaching photomultiplier tube.
To estimate influence of the shutter used on thermoluminescence signal integrals, additional experiment was performed. Commercially available LiF dosimeters (MCP-N) were irradiated with different doses and read out with and without shutter. The slopes of the functions integrals vs. dose were calculated. Ratio of slopes between these two reading conditions determined the ability of the shutter to reduce light. For the shutter used, TL signal intensity was reduced by two orders of magnitude. Due to much lower light-emission intensity in case of undoped (raw) phosphate glass, the shutter was not used.

Results of XRD and DSC Analysis of Synthesized Glasses
All synthesized samples were analyzed using X-ray powder diffraction technique to ensure their glassy state. As one can see from the diffraction patterns presented in Figure 1a, one should consider all phosphate glass samples as fully amorphous materials. The lack of crystalline peaks is observed for all samples, even for PPG_Tb30 sample, containing relatively high contribution (3.0%, w/w) of Tb 2 O 3 oxide. The presented XRD patterns show that for all samples only wide scattering peaks (so called "halo") are observed at ca. 21 • and 31 • 2theta angles. In the low range of 2theta angle (5 • to 16 • ), one can also observe an increased background intensity due to more efficient X-ray scattering in air. Obtained X-ray patterns are similar to those that can be find in literature [31][32][33]. The results of DSC analysis are presented in Figure 1b. Thermal analyses were carried out within the temperature range of 40-550 • C and confirmed the glassy nature of the investigated materials by existing of the glass transitions within the temperature range of 407-423 • C. For pure, undoped phosphate glass, the glass transition temperature was determined at 423 • C, whereas for the terbiumdoped samples the transition temperatures were shifted towards lower values between 408-410 • C. In terms of DSC, for undoped phosphate glass exhibits there was similar T g to that reported in literature [34]. Decreasing of T g temperature for terbium-doped glasses may be caused by a global weakening of the glass network induced by relatively large Tb 3+ ion incorporation, which can be ascribed to the increase in non-bridging oxygens due to the cleavage of P-O-P bonds to form P-O-Tb 2+ [35][36][37][38]. However, this issue requires further research using other research methods, for example, NMR.

Thermoluminescence (TL) Characteristics of Irradiated Glasses
The phosphate glass dosimeters samples were irradiated in different modes in order to investigate the detectors response. For this purpose, both electron beam and gamma photon irradiation were performed on investigated samples.

Electron Accelerator, Pulse Duration 4 μs
At first, the native (undoped) and terbium-doped phosphate glass samples were irradiated using electron accelerator with pulse duration of 4 μs. The range of radiation doses delivered to the samples was 1-30 kGy. The results of these experiments are depicted in Figure 2, where TL glow curves for each sample is presented for different radiation doses. The undoped phosphate glass already exhibits thermoluminescent properties. The irradiated samples of PPG glass show wide, complex TL signals with maximum emission at ca. 125 °C (cf. Figure 2a). For the higher radiation doses, it is also possible to identify higher temperature peak, which in TL curves is located at ca. 180 °C (Figure 2a). For all terbium-doped samples, special shutter was installed to decrease amount of light reaching PMT tube (ca. 100 times). As one can see in Figure 2b, 0.1% addition of terbium oxide increases the yield of thermoluminescence significantly in comparison with native phosphate glass. Higher contribution of terbium results in further increasing of TL intensity, which can be clearly observed in Figure 2c,d. For the glass samples doped with 0.5% of terbium oxide, the maximum of signals appears at ca. 225 °C (Figure 2b,c). Further increasing of terbium concentration in glass, to 3%, is followed by formation of additional signal at higher temperatures. This is clearly presented in Figure 2d, where one can observe low temperature signal at ca. 125 °C with easily recognized signal at 275 °C. Signal at 125 °C is common for all samples irradiated with different doses, whereas peak at 275 °C appears for radiation doses higher than 8 kGy. Moreover, for radiation dose of 31 kGy one can observed partial degradation of the low temperature signal. Additionally PPG_Tb01and PPG_Tb30 shows an extensive change in the peak positions at higher radiation dose than in PPG_Tb05. Such effect may be due to the intensity of glass network modification process caused by introduction of terbium dopant. In case of low content of terbium, one deals only with slight modification of the glass network, whereas when increasing terbium oxide concentration one can expect more significant change of glass matrix structure and weakening of the glass-network-forming bonds. Such processes may lead to formation of additional trap centers, which can be more easily bleached when irradiated with higher doses. For the presented glow curves, one can also observe the slight distortions of the signal shape related to the increasing radiation dose. Nevertheless, the overall character of the presented glow curves remains similar for each glass sample. This

Thermoluminescence (TL) Characteristics of Irradiated Glasses
The phosphate glass dosimeters samples were irradiated in different modes in order to investigate the detectors response. For this purpose, both electron beam and gamma photon irradiation were performed on investigated samples.

Electron Accelerator, Pulse Duration 4 µs
At first, the native (undoped) and terbium-doped phosphate glass samples were irradiated using electron accelerator with pulse duration of 4 µs. The range of radiation doses delivered to the samples was 1-30 kGy. The results of these experiments are depicted in Figure 2, where TL glow curves for each sample is presented for different radiation doses. The undoped phosphate glass already exhibits thermoluminescent properties. The irradiated samples of PPG glass show wide, complex TL signals with maximum emission at ca. 125 • C (cf. Figure 2a). For the higher radiation doses, it is also possible to identify higher temperature peak, which in TL curves is located at ca. 180 • C (Figure 2a). For all terbiumdoped samples, special shutter was installed to decrease amount of light reaching PMT tube (ca. 100 times). As one can see in Figure 2b, 0.1% addition of terbium oxide increases the yield of thermoluminescence significantly in comparison with native phosphate glass. Higher contribution of terbium results in further increasing of TL intensity, which can be clearly observed in Figure 2c,d. For the glass samples doped with 0.5% of terbium oxide, the maximum of signals appears at ca. 225 • C (Figure 2b,c). Further increasing of terbium concentration in glass, to 3%, is followed by formation of additional signal at higher temperatures. This is clearly presented in Figure 2d, where one can observe low temperature signal at ca. 125 • C with easily recognized signal at 275 • C. Signal at 125 • C is common for all samples irradiated with different doses, whereas peak at 275 • C appears for radiation doses higher than 8 kGy. Moreover, for radiation dose of 31 kGy one can observed partial degradation of the low temperature signal. Additionally PPG_Tb01and PPG_Tb30 shows an extensive change in the peak positions at higher radiation dose than in PPG_Tb05. Such effect may be due to the intensity of glass network modification process caused by introduction of terbium dopant. In case of low content of terbium, one deals only with slight modification of the glass network, whereas when increasing terbium oxide concentration one can expect more significant change of glass matrix structure and weakening of the glass-network-forming bonds. Such processes may lead to formation of additional trap centers, which can be more easily bleached when irradiated with higher doses. For the presented glow curves, one can also observe the slight distortions of the signal shape related to the increasing radiation dose. Nevertheless, the overall character of the presented glow curves remains similar for each glass sample. This effect may be related to the possible thermal and geometrical conditions during readout and/or imperfect grain size distribution of the material used for detectors preparation, which was already observed previously [39]. As one can see in Figure 2, introduction of terbium ions into the glass matrix significantly modifies thermoluminescence signals, both in terms of their shapes and intensity. Thermoluminescence signals at 125 and 180 • C may have their origins from pure glass lattice. For doped glasses, one can observe additional signals at 225 and 275 • C that are clearly caused by presence of Tb. The more detailed explanation for these effects is presumably related to the facts-that aluminophosphate glasses contain PO 4 tetrahedra as well as AlO 4 units, which are able to induce point defects owing to the fact of mismatching in bonding as well as formation links with neighboring phosphorous elements. Due to the presence of point defects, the simultaneous occurrence of localized energy levels within the band gap is possible [40].
Appl. Sci. 2021, 11, x FOR PEER REVIEW 6 of 21 effect may be related to the possible thermal and geometrical conditions during readout and/or imperfect grain size distribution of the material used for detectors preparation, which was already observed previously [39]. As one can see in Figure 2, introduction of terbium ions into the glass matrix significantly modifies thermoluminescence signals, both in terms of their shapes and intensity. Thermoluminescence signals at 125 and 180 °C may have their origins from pure glass lattice. For doped glasses, one can observe additional signals at 225 and 275 °C that are clearly caused by presence of Tb. The more detailed explanation for these effects is presumably related to the facts-that aluminophosphate glasses contain PO4 tetrahedra as well as AlO4 units, which are able to induce point defects owing to the fact of mismatching in bonding as well as formation links with neighboring phosphorous elements. Due to the presence of point defects, the simultaneous occurrence of localized energy levels within the band gap is possible [40]. Rare earth metals, such as terbium, induce higher sensitivity in TL materials [21]. From literature it is known that matrix distortions due to presence of large dopant ions such as Tb 3+ may form complex defect sites, which can be stabilized by Na + and intrinsic oxygen vacancies [41]. The possible mechanism of TL emission from terbium-doped phosphors includes few steps. It is considered that trivalent terbium ions Tb 3+ may be effectively reduced by free electrons formed in the matrix after β-and γ irradiation. Further, the formed Tb 2+ ion is efficient hole scavenger, forming excited Tb 3+ ions in this process, and TL emission is considered to follow the reaction expressed by Equation (1).
The holes are supplied by traps via thermal excitation [42]. An elegant explanation for these mechanisms based on the simultaneous participation of rare earth metals ions and interstitial oxygen atoms is presented in the Abdel-Kader paper [43]. The complex nature of energy transitions between excited states and ground level of Tb3+ ions involves different transitions between D and F levels (especially 5 D4→ 7 F6, 5 D4→ 7 F5, 5 D4→ 7 F4 and Rare earth metals, such as terbium, induce higher sensitivity in TL materials [21]. From literature it is known that matrix distortions due to presence of large dopant ions such as Tb 3+ may form complex defect sites, which can be stabilized by Na + and intrinsic oxygen vacancies [41]. The possible mechanism of TL emission from terbium-doped phosphors includes few steps. It is considered that trivalent terbium ions Tb 3+ may be effectively reduced by free electrons formed in the matrix after βand γ irradiation. Further, the formed Tb 2+ ion is efficient hole scavenger, forming excited Tb 3+ ions in this process, and TL emission is considered to follow the reaction expressed by Equation (1).
The holes are supplied by traps via thermal excitation [42]. An elegant explanation for these mechanisms based on the simultaneous participation of rare earth metals ions and interstitial oxygen atoms is presented in the Abdel-Kader paper [43]. The complex nature of energy transitions between excited states and ground level of Tb3+ ions involves different transitions between D and F levels (especially 5 D 4 → 7 F 6 , 5 D 4 → 7 F 5 , 5 D 4 → 7 F 4 and 5 D 4 → 7 F 3 ), which are then responsible for enhanced photon emission from terbium-doped materials [44].
As one can see in Figure 3, dependencies of signal integral vs. dose cannot be considered as linear and are clearly saturation-type curves, which resemble functions of exponential type. The integrated luminescence intensity (N int ) may be approximated by a simple summation, according to Equation (2) [39]: where Np is the photon count intensity at point temperature T. For the approximation of the experimental data presented in Figure 2, the function expressed by Equation (3) can be used: where N inf is the pre-exponential factor dependent on the material properties and related to the saturation intensity, k is the thermoluminescence yield per unit dose factor or material sensitivity coefficient, in kGy-1, and D is the radiation dose in kGy. It is obvious that increasing the Tb 3+ ions concentration has the most pronounced influence on the lightemission efficiency, what is also confirmed by the data provided in the further parts of this work. The sensitizing effect of rare earth elements was already discussed in previous papers, where other glass samples containing such types of dopants were investigated [27,[45][46][47]. It seems that terbium is one of the most efficient sensitizing agent among rare earth metals. However, for the terbium ions the strongest enhancement of the light-emission efficiency is observed for the low dopant concentration range to upmost few percent and decrease above this value. Above relatively low concentrations of ca. 0.3%, one observed the quenching effects due to simultaneous effects of multipolar transitions or these related to the so called inner effect filter, which were observed already for terbium concentration starting from 6 (mole) % [46,47]. The data obtained and presented in this study correlate well with those provided in literature, which is demonstrated by the results provided in Figure 3, such as the significant increase of TL intensity when shifting from PPG_Tb01 sample to PPG_Tb05 and lower efficiency for further transition to PPG_Tb30 glass.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 7 of 21 5 D4→ 7 F3), which are then responsible for enhanced photon emission from terbium-doped materials [44]. As one can see in Figure 3, dependencies of signal integral vs. dose cannot be considered as linear and are clearly saturation-type curves, which resemble functions of exponential type. The integrated luminescence intensity (Nint) may be approximated by a simple summation, according to Equation (2) [39]: where Np is the photon count intensity at point temperature T. For the approximation of the experimental data presented in Figure 2, the function expressed by Equation (3) can be used: where Ninf is the pre-exponential factor dependent on the material properties and related to the saturation intensity, k is the thermoluminescence yield per unit dose factor or material sensitivity coefficient, in kGy-1, and D is the radiation dose in kGy. It is obvious that increasing the Tb 3+ ions concentration has the most pronounced influence on the lightemission efficiency, what is also confirmed by the data provided in the further parts of this work. The sensitizing effect of rare earth elements was already discussed in previous papers, where other glass samples containing such types of dopants were investigated [27,[45][46][47]. It seems that terbium is one of the most efficient sensitizing agent among rare earth metals. However, for the terbium ions the strongest enhancement of the light-emission efficiency is observed for the low dopant concentration range to upmost few percent and decrease above this value. Above relatively low concentrations of ca. 0.3%, one observed the quenching effects due to simultaneous effects of multipolar transitions or these related to the so called inner effect filter, which were observed already for terbium concentration starting from 6 (mole) % [46,47]. The data obtained and presented in this study correlate well with those provided in literature, which is demonstrated by the results provided in Figure 3, such as the significant increase of TL intensity when shifting from PPG_Tb01 sample to PPG_Tb05 and lower efficiency for further transition to PPG_Tb30 glass. Figure 3. Dependencies of TL signal integral vs. dose absorbed for glass dosimeter irradiated with beta radiation (4 μs); for PPG sample, TL glow curves were recorded without shutter. Figure 3. Dependencies of TL signal integral vs. dose absorbed for glass dosimeter irradiated with beta radiation (4 µs); for PPG sample, TL glow curves were recorded without shutter. Considering samples containing terbium, one can observe that with increasing amount of Tb 2 O 3 in the glass, the yield of luminescence increases. As can be found from data included in Table 2, despite the shutter, samples of terbium-doped glass possess similar or even higher coefficients than undoped glass that was read without shutter, with observable tendency to increase with the terbium concentration (except of k coefficient value for the sample PPG_Tb01).  For verification of the behavior of the investigated glasses for the lower dose ranges, samples were irradiated with radiation doses up to 1000 Gy. During this experiment, samples were irradiated by electron accelerator with pulse duration of 17 ns.
As shown in Figure 4a, the raw phosphate glass possesses maximum of TL signal at 125 • C. For radiation dose of ca. 2 kGy, one can easily find additional maximum, which appears at ca. 210 • C; however, after closer analysis this signal is also observable for the lower doses. For the glass samples containing terbium (Figure 4b-d), one can observe wide TL signals from 100 to 250 • C, whose shapes suggest more complex structure. Probably next to the main signals at 125 • C (signal from glass matrix) and at 200 • C, other TL peaks also contribute to the overall emission from terbium-doped glass. This will be further discussed in Section 3.3.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 8 of Considering samples containing terbium, one can observe that with increasin amount of Tb2O3 in the glass, the yield of luminescence increases. As can be found fro data included in Table 2, despite the shutter, samples of terbium-doped glass possess sim ilar or even higher coefficients than undoped glass that was read without shutter, wi observable tendency to increase with the terbium concentration (except of k coefficie value for the sample PPG_Tb01).  For verification of the behavior of the investigated glasses for the lower dose range samples were irradiated with radiation doses up to 1000 Gy. During this experiment, sam ples were irradiated by electron accelerator with pulse duration of 17 ns.
As shown in Figure 4a, the raw phosphate glass possesses maximum of TL signal 125 °C. For radiation dose of ca. 2 kGy, one can easily find additional maximum, whic appears at ca. 210 °C; however, after closer analysis this signal is also observable for th lower doses. For the glass samples containing terbium (Figure 4b-d), one can observ wide TL signals from 100 to 250 °C, whose shapes suggest more complex structure. Pro ably next to the main signals at 125 °C (signal from glass matrix) and at 200 °C, other T peaks also contribute to the overall emission from terbium-doped glass. This will be fu ther discussed in Section 3.3.  As one can see in Figure 5, all samples exhibit linear dependencies of TL signal intensity vs. dose in the dose range from 50 to 2000 Gy. Analyzing data included in Table 3, one can see that yield of TL emission is directly linked with amount of dopant in the glass, analogical as in case of higher dose ranges, for which the results were presented in Section 3.2.1. With increasing concentration of rare earth oxide, the total TL emission increases significantly.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 21 Figure 5. Dependencies of TL signal integral vs. dose absorbed for glass dosimeter irradiated with beta radiation (17ns); for PPG sample, TL glow curves were recorded without shutter. Gamma irradiation was performed using Co-60 source (Ob-Servo-D panoramic irradiator). The samples were placed on sources housing surrounding Co-60 radiation sources.
As one can see in Figure 6a, raw phosphate glass possesses maximum of TL signal at 140 °C, that is, at temperature of ca. 15 °C higher in comparison with electron beam irradiated samples. Additional peaks appear at 210 °C and 280 °C. For samples containing terbium (Figure 6b-d), one can observe a wide TL complex signal from 100 to 300 °C. For samples containing 3% of terbium oxide, additional maximum at 300 °C is observed for doses above 5 kGy. Comparing the data provided in Figures 4 and 6, one can easily find that native PPG glasses exhibit weak or moderate light-emission ability, whereas it is not observed for the terbium-doped samples. Taking this observation under analysis, it is possible to explain such behavior by two effects. On one hand, the PPG glasses were measured without light reducing shutter, due to significantly lower light-emission efficiency as compared to Tb 3+ doped glasses, so even weaker signals may be easily observed, especially for higher radiation doses. On the other hand, as a reasonable explanation for the lack of light emission in case of terbium-doped glasses, it is possible that incorporation of the Tb 3+ ions into glass matrix modifies the trap structure in such a way that deep traps existing in native PPG glass are converted into traps, which are emptied almost completely at lower temperatures.  Gamma irradiation was performed using Co-60 source (Ob-Servo-D panoramic irradiator). The samples were placed on sources housing surrounding Co-60 radiation sources.
As one can see in Figure 6a, raw phosphate glass possesses maximum of TL signal at 140 • C, that is, at temperature of ca. 15 • C higher in comparison with electron beam irradiated samples. Additional peaks appear at 210 • C and 280 • C. For samples containing terbium (Figure 6b-d), one can observe a wide TL complex signal from 100 to 300 • C. For samples containing 3% of terbium oxide, additional maximum at 300 • C is observed for doses above 5 kGy. Comparing the data provided in Figures 4 and 6, one can easily find that native PPG glasses exhibit weak or moderate light-emission ability, whereas it is not observed for the terbium-doped samples. Taking this observation under analysis, it is possible to explain such behavior by two effects. On one hand, the PPG glasses were measured without light reducing shutter, due to significantly lower light-emission efficiency as compared to Tb 3+ doped glasses, so even weaker signals may be easily observed, especially for higher radiation doses. On the other hand, as a reasonable explanation for the lack of light emission in case of terbium-doped glasses, it is possible that incorporation of the Tb 3+ ions into glass matrix modifies the trap structure in such a way that deep traps existing in native PPG glass are converted into traps, which are emptied almost completely at lower temperatures.  As can be seen in Figure 7, all samples doped with terbium oxide exhibit close to linear dependencies of TL signal intensity vs. dose in the dose range from 1 to 8 kGy. One can see in Table 4 that also for gamma irradiation, yield of TL emission is linked with amount of dopant in glass. For all samples, goodness of function fitting is near 1, except for samples of PPG (undoped) glass-this is probably caused by incorrect position of these kinds of detectors during irradiation. As can be seen in Figure 7, all samples doped with terbium oxide exhibit close to linear dependencies of TL signal intensity vs. dose in the dose range from 1 to 8 kGy. One can see in Table 4 that also for gamma irradiation, yield of TL emission is linked with amount of dopant in glass. For all samples, goodness of function fitting is near 1, except for samples of PPG (undoped) glass-this is probably caused by incorrect position of these kinds of detectors during irradiation. As can be seen in Figure 7, all samples doped with terbium oxide exhibit close to linear dependencies of TL signal intensity vs. dose in the dose range from 1 to 8 kGy. One can see in Table 4 that also for gamma irradiation, yield of TL emission is linked with amount of dopant in glass. For all samples, goodness of function fitting is near 1, except for samples of PPG (undoped) glass-this is probably caused by incorrect position of these kinds of detectors during irradiation.

Analysis and Structure of TL Emission from Phosphate Glasses
As one can find, from the previous data, changing concentration of terbium oxide in phosphate glasses leading to increase of their thermoluminescence yield and also influencing the shapes of TL curves. Because of this, it was necessary to investigate TL mechanism in synthesized glasses. To follow and analyze the impact of terbium dopant and its concentration on TL signal and traps structure, one can apply verified approach based on the mixed experimental methods and computerized deconvolution of the TL curves. The substantial question was the selection of the proper traps distribution model, which should be selected for glow curves deconvolution procedure. Thus, at first one performed T max -T stop measurements for both pure PPG and terbium-doped glass samples irradiated with 10 kGy radiation dose, within the temperature ranges of 50-360 • C and changing temperature in 2-6 • C steps. The results obtained from these experiments are presented in Figure 8. From the presented data, one can easily determine the number and nature of traps present in the investigated glasses. For pure PPG glass sample (Figure 8a), the dependence of position of the first maximum present in a glow curve (T max ) obtained during reading the sample after irradiation and preheating to T stop temperature vs. T stop revealed the existence of five trapping centers. The first peak represented by a flat line in region of 50-90 • C was characterized as a single first order kinetics (FOK) signal, next three signals were identified as defects with continuous distribution of energy (straight increasing lines in regions 90-140, 140-190, and 190-310 • C), whereas the last peak visible as a flat region followed by an increasing tail gave presumption for existence of a single general order kinetics (GOK) peak in the region of 280-360 • C [48]. Similar analysis of the terbium-doped glass delivered evidence for existence of one FOK peak in the region of

Analysis and Structure of TL Emission from Phosphate Glasses
As one can find, from the previous data, changing concentration of terbium oxid phosphate glasses leading to increase of their thermoluminescence yield and also i encing the shapes of TL curves. Because of this, it was necessary to investigate TL me nism in synthesized glasses. To follow and analyze the impact of terbium dopant an concentration on TL signal and traps structure, one can apply verified approach base the mixed experimental methods and computerized deconvolution of the TL curves substantial question was the selection of the proper traps distribution model, w should be selected for glow curves deconvolution procedure. Thus, at first one perfor Tmax-Tstop measurements for both pure PPG and terbium-doped glass samples irrad with 10 kGy radiation dose, within the temperature ranges of 50-360 °C and chan temperature in 2-6 °C steps. The results obtained from these experiments are prese in Figure 8. From the presented data, one can easily determine the number and natu traps present in the investigated glasses. For pure PPG glass sample (Figure 8a), th pendence of position of the first maximum present in a glow curve (Tmax) obtained du reading the sample after irradiation and preheating to Tstop temperature vs. Tstop reve the existence of five trapping centers. The first peak represented by a flat line in regio 50-90 °C was characterized as a single first order kinetics (FOK) signal, next three sig were identified as defects with continuous distribution of energy (straight increasing in regions 90-140, 140-190, and 190-310 °C), whereas the last peak visible as a flat re followed by an increasing tail gave presumption for existence of a single general o kinetics (GOK) peak in the region of 280-360 °C [48]. Similar analysis of the terb doped glass delivered evidence for existence of one FOK peak in the region of 50-9 and five peaks with a continuous trap energy distribution in regions of 90-110, 110-130-160, 160-190, and 190-320 °C. Direct determination of the exact positions of signals contributing to a complex g curve with strongly overlapping peaks is always a difficult task, so next to the prese above Tmax-Tstop measurements, one also performed variable heating rate method to itate the resolving of the glow curves structure. The method was applied both for phosphate glass (PPG) and for terbium-doped samples within the heating rates betw Direct determination of the exact positions of signals contributing to a complex glow curve with strongly overlapping peaks is always a difficult task, so next to the presented above T max -T stop measurements, one also performed variable heating rate method to facilitate the resolving of the glow curves structure. The method was applied both for pure phosphate glass (PPG) and for terbium-doped samples within the heating rates between 0.5-10 • C on samples irradiated with 10 kGy radiation dose. Results of these measurements are presented in Figure 9. In case of terbium-doped samples, shutter-reducing light intensity had to be applied. From the VHR measurements, the shifting of TL emission maxima is evident. For PPG sample (cf. Figure 9a), the most prominent peak maximum moves from 98 to 145 • C, upon changing heating rate from 0.5 to 10 • C/s. Unfortunately, due to the strong overlapping of peaks, this is the only signal whose position may be unambiguously determined. Nevertheless, the overall shapes of the TL emission curves remain similar, regardless of heating rate applied.

PEER REVIEW 12 of 21
intensity had to be applied. From the VHR measurements, the shifting of TL emission maxima is evident. For PPG sample (cf. Figure 9a), the most prominent peak maximum moves from 98 to 145 °C, upon changing heating rate from 0.5 to 10 °C/s. Unfortunately, due to the strong overlapping of peaks, this is the only signal whose position may be unambiguously determined. Nevertheless, the overall shapes of the TL emission curves remain similar, regardless of heating rate applied. For terbium-doped glass (0.5% Tb2O3, PPG_Tb05, cf. Figure 9b) and the lowest heating rate of 0.5 °C/s, one can easily determine the low temperature peak position at 101 °C/s and roughly the position of the other, higher temperature peak at ca. 175 °C. Moreover, for terbium-doped glass the much stronger dependence of the overall TL emission curves vs. heating rate can be easily found. Increasing the heating rate, the intensity of the first peak decreases progressively; however, it can be still easily found even in case of the highest heating rate of 10 °C/s on the lower temperature shoulder of the most prominent TL emission peak, whose maximum is located at ca. 200 °C/s. Very similar results based on the VHR method were observed also for the 0.1 and 3.0% Tb2O3-doped phosphate glasses (PPG_Tb01 and PPG_Tb30).
Due to the significant complexity of the presented TL emissions, it is very difficult to completely analyze the structure of the glow curves and contribution of the particular TL peaks taking into account only results obtained from Tmax-Tstop and VHR measurements. Nevertheless, they can be a starting point for the deconvolution analysis using specialized software with implemented appropriate TL models. Based on the data presented above, one can propose models suitable for resolving observed TL glow curves. For the PPG glass, the model involves mixed contribution of FOK peak, three continuous energy distribution peaks, and GOK signal, according to Equation (4): where ITL and IPeak (1)(2)(3)(4)(5) are total and fractional TL intensities due to the particular peaks For terbium-doped glass (0.5% Tb 2 O 3 , PPG_Tb05, cf. Figure 9b) and the lowest heating rate of 0.5 • C/s, one can easily determine the low temperature peak position at 101 • C/s and roughly the position of the other, higher temperature peak at ca. 175 • C. Moreover, for terbium-doped glass the much stronger dependence of the overall TL emission curves vs. heating rate can be easily found. Increasing the heating rate, the intensity of the first peak decreases progressively; however, it can be still easily found even in case of the highest heating rate of 10 • C/s on the lower temperature shoulder of the most prominent TL emission peak, whose maximum is located at ca. 200 • C/s. Very similar results based on the VHR method were observed also for the 0.1 and 3.0% Tb 2 O 3 -doped phosphate glasses (PPG_Tb01 and PPG_Tb30).
Due to the significant complexity of the presented TL emissions, it is very difficult to completely analyze the structure of the glow curves and contribution of the particular TL peaks taking into account only results obtained from T max -T stop and VHR measurements. Nevertheless, they can be a starting point for the deconvolution analysis using specialized software with implemented appropriate TL models. Based on the data presented above, one can propose models suitable for resolving observed TL glow curves. For the PPG glass, the model involves mixed contribution of FOK peak, three continuous energy distribution peaks, and GOK signal, according to Equation (4): where I TL and I Peak (1)(2)(3)(4)(5) are total and fractional TL intensities due to the particular peaks contributions. Such a combined approach to TL glow curve analysis has already been applied and reported in literature [49]. The glow peaks following first (FOK) and general order kinetics (GOK) were fitted using Equations (5) and (6) [48,50]: where I and I m are TL intensities for a given temperature (T, in K) or at peak maximum (Tm) and temperature linearly increases with time (t) from T 0 to T, T(t) = T 0 + βt, β ( • C/s) is the linear heating rate; E (eV) is the trap activation energy; k (eV·K −1 ) is the Boltzmann constant; and b is the kinetic order parameter. In case of continuous energy distribution of traps present in material, usually exponential or gaussian functions are used, which can be described by Equations (7) and (8) [49][50][51]: where n and no are number of trapping centres of a given energy (E) and total density of trapped charges (cm −3 ), E 0 is either the shallowest energy in the distribution (exponential model) or the most populated trapping level (gaussian distribution), and σ is the parameter related to the width of the distribution. Based on these distributions, two models are generally in use, describing dependencies of TL emission intensity vs. energy and temperature according to Equations (9) and (10): (10) where s (s −1 ) is the frequency coefficient related to the type of defect and other symbols are as described above. The models described according to Equations (5)-(10) are implemented in a dedicated glow curve deconvolution software and are described in several papers [49,50]. The software was used for analysis of glow curves obtained during measurement of PPG sample, and graphical representation of the obtained results are presented in Figure 10. The deconvolution procedure was performed in order to recognize the structure of the TL glow curve, structure, and types of defects present in the material and to estimate the energy levels of the trapping centers related to the particular TL emissions from the material. Verification of the proper localization of particular peaks was based on the consecutive analysis of the experimental data and results obtained from deconvolution analysis using the GCA (Glow Curve Analysis) software [50]. The selection between exponential and gaussian distributions for continuous type defects was made taking as a criterion the possibly lowest value of FOM (Figure of Merit) parameter, which indicates the difference between experimental data and computer-simulated glow curves based on the parameters obtained during deconvolution procedure. The commonly accepted values of FOM are usually below 5%. The results obtained during deconvolution analysis of glow curves obtained from VHR measurements of PPG glass are presented in Table 5. √2πσ where s (s −1 ) is the frequency coefficient related to the type of defect and other symbols are as described above. The models described according to Equations (5)-(10) are implemented in a dedicated glow curve deconvolution software and are described in several papers [49,50]. The software was used for analysis of glow curves obtained during measurement of PPG sample, and graphical representation of the obtained results are presented in Figure 10. The deconvolution procedure was performed in order to recognize the structure of the TL glow curve, structure, and types of defects present in the material and to estimate Figure 10. Results of deconvolution procedure of TL signal from PPG samples irradiated with 10 kGy radiation dose (electron accelerator, pulse duration 4 µs, heating rate β = 0.5 • C/s). Table 5. Fitting results and traps parameters obtained for PPG glass samples irradiated with radiation dose of ca. 10 kGy and heating rates (β) of 0.5, 1.0, and 2.0 • C/s. In Table 5, one presented only data obtained during deconvolution of TL emission curves for three lowest heating rates of 0.5, 1.0, and 2.0 • C/s. The general conclusion from these results is that all sets of data correlate quite well each other for particular peaks. For Peak 2, one can find relatively high value of the frequency factor; however, this could be due to the exponential expression used to obtain this value and based on energy of the trapping center, so even small variation in E value may lead to large variation and uncertainty in calculated value of s. It was impossible to correlate all the data obtained from deconvolution procedure with the results from VHR measurements, especially for the higher heating rates above 2 • C/s. It was probably due to a limited thermal conductivity of the samples (ca. 40 mg in weight, 5 mm in diameter, and ca. 1 mm thickness) and in consequence delayed response of the TL emission vs. temperature measured during samples reading. Nevertheless, the careful analysis of the available data allowed for determination of the positions of the peaks in TL emission curves for which temperature shifts were monitored in relation to the applied heating rate. Analysis was performed according to the known VHR based correlation of the TL peak position and heating rate, which can be described by Equation (11) [49]:

Peak Type
where all symbols are as described previously. Equation (11) can be easily linearized, and the obtained dependencies of ln T 2 m /β = f(1/T) for PPG glass samples are presented in Figure 11. As mentioned above, due to limited thermal conductivity of the samples, it was very difficult to correlate results obtained from the higher heating rates, for which calculated values of departed significantly from linearity. Such prepared data allowed for experimental verification of the values obtained during deconvolution procedure. For all presented dependencies, one can observe purely linear correlations, whose slopes allowed for experimental determination of the energies of trapping centers, presented in Table 6 together (for comparison) with those obtained from deconvolution procedure. In most cases, the presented values are fully convergent, with only higher discrepancy observed for Peak 3, in which case the value obtained from VHR method was ca. 15% lower than that from deconvolution analysis.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 15 of 21 very difficult to correlate results obtained from the higher heating rates, for which calculated values of departed significantly from linearity. Such prepared data allowed for experimental verification of the values obtained during deconvolution procedure. For all presented dependencies, one can observe purely linear correlations, whose slopes allowed for experimental determination of the energies of trapping centers, presented in Table 6 together (for comparison) with those obtained from deconvolution procedure. In most cases, the presented values are fully convergent, with only higher discrepancy observed for Peak 3, in which case the value obtained from VHR method was ca. 15% lower than that from deconvolution analysis.  Taking into account the data from Tmax-Tstop and VHR experiments, it seems that structure of the TL signals in terbium-doped phosphate glasses should exhibit similar behavior as observed in case of pure PPG glass. Because of a vast number of experiments required and calculation works to be performed, one decided that due to the complexity of the thermoluminescence signals from terbium-doped phosphate glasses, this topic deserves wider discussion and analysis in another separate publication, for which the results  Taking into account the data from T max -T stop and VHR experiments, it seems that structure of the TL signals in terbium-doped phosphate glasses should exhibit similar behavior as observed in case of pure PPG glass. Because of a vast number of experiments required and calculation works to be performed, one decided that due to the complexity of the thermoluminescence signals from terbium-doped phosphate glasses, this topic deserves wider discussion and analysis in another separate publication, for which the results presented in this work provide a great introduction. Because of the presence of large number of signals at lower temperatures, one can expect that the signal may be degraded during time due to the random recombinations and emission processes. After verification, this supposition was positively confirmed in experiment, whose results are shown further in Section 3.4.2.
Thus, we decided to check the effect of preheating detectors prior to readout and its influence on the obtained TL curves. The preheating and TL signal stability experiments were performed on the PPG_Tb30 samples. Glass detectors were preheated directly before readout. Two temperatures were used for preheating, namely, 120 • C and 200 • C. The results of these experiments are presented in Figure 12. From the data presented in Figure 8, it is evident that preheating erases part of the TL signal, which cannot be obtained again. In case of the not preheated sample PPG_Tb30, at first sight one can observe three main peaks, at ca. 120 • C, 180 • C, and 260 • C. During preheating at 120 • C, two or even three low-temperature peaks, as presented in Figure 8b, are cleaned by emptying traps presented in the material. It is clear that removal of the lower temperature signals leads to situation, in which, after annealing, they do not contribute to the overall TL emission curve. Similarly, preheating to higher temperature of 200 • C consequently must lead to emptying of the five traps with lower energies existing in the material, as shown in Figure 8. Because of the presence of large number of signals at lower temperatures, one can expect that the signal may be degraded during time due to the random recombinations and emission processes. After verification, this supposition was positively confirmed in experiment, whose results are shown further in Section 3.4.2.
Thus, we decided to check the effect of preheating detectors prior to readout and its influence on the obtained TL curves. The preheating and TL signal stability experiments were performed on the PPG_Tb30 samples. Glass detectors were preheated directly before readout. Two temperatures were used for preheating, namely, 120 °C and 200 °C. The results of these experiments are presented in Figure 12. From the data presented in Figure  8, it is evident that preheating erases part of the TL signal, which cannot be obtained again. In case of the not preheated sample PPG_Tb30, at first sight one can observe three main peaks, at ca. 120 °C, 180 °C, and 260 °C. During preheating at 120 °C, two or even three low-temperature peaks, as presented in Figure 8b, are cleaned by emptying traps presented in the material. It is clear that removal of the lower temperature signals leads to situation, in which, after annealing, they do not contribute to the overall TL emission curve. Similarly, preheating to higher temperature of 200 °C consequently must lead to emptying of the five traps with lower energies existing in the material, as shown in Figure  8. The overall influence of the preheating procedure on the radiation dose and integrated TL intensity response is presented in Figure 13. Obviously, the TL signal integrals are more intensive for the non-preheated samples. This is a logical consequence of partial emptying of traps through preheating. Interestingly, preheating also has an influence on The overall influence of the preheating procedure on the radiation dose and integrated TL intensity response is presented in Figure 13. Obviously, the TL signal integrals are more intensive for the non-preheated samples. This is a logical consequence of partial emptying of traps through preheating. Interestingly, preheating also has an influence on the dependencies of TL signal integrals vs. radiation dose absorbed. For the non-preheated samples, the dose response is of exponential character (saturation-type curve, triangles in Figure 13), for which signal saturation is observed above 15 kGy. For the sample preheated at 120 • C, the dose response is of linear character for lower dose range (up to 15 kGy), and for higher doses exhibits it exponential character (circles in Figure 13). Finally, after preheating at 200 • C, the presented in Figure 13 data representing the dependence of TL signal integral vs. dose absorbed remain linear even for the highest dose up to ca. 30 kGy. For comparison purposes, in Figure 13 both exponential fittings as well as linear approximations are shown for a clear presentation of the observed differences and effects of samples preheating at different conditions. After this analysis it is clear that preheating efficiently erases part of the lower energy traps, that, as one can observe, affects significantly the linearity of the dependence of integrated TL signal vs. radiation dose. With increasing of the absorbed dose, the low temperature defects saturate more effectively in comparison with those corresponding to higher temperatures. In such case, high-energy defects play a major role in TL signal generation. heated samples, the dose response is of exponential character (saturation-type curve, triangles in Figure 13), for which signal saturation is observed above 15 kGy. For the sample preheated at 120 °C, the dose response is of linear character for lower dose range (up to 15 kGy), and for higher doses exhibits it exponential character (circles in Figure 13). Finally, after preheating at 200 °C, the presented in Figure 13 data representing the dependence of TL signal integral vs. dose absorbed remain linear even for the highest dose up to ca. 30 kGy. For comparison purposes, in Figure 13 both exponential fittings as well as linear approximations are shown for a clear presentation of the observed differences and effects of samples preheating at different conditions. After this analysis it is clear that preheating efficiently erases part of the lower energy traps, that, as one can observe, affects significantly the linearity of the dependence of integrated TL signal vs. radiation dose.
With increasing of the absorbed dose, the low temperature defects saturate more effectively in comparison with those corresponding to higher temperatures. In such case, highenergy defects play a major role in TL signal generation.

Effect of Preheating on Signal Loss in Time
Another important aspect in radiation dosimetry is the stability of signal generated upon irradiation within the detector used. To investigate the signal loss in time, samples of PPG_Tb30 were irradiated using LINAC with dose of 3.65 kGy. After irradiation, samples were stored in dark and preheated right before readout. To compare the obtained results in a quantitative way, Equation (12) was used: where SL is the TL signal loss in time, in %, and Nint and Nint(to) are the integrated TL signal intensities at a given time and just after irradiation. As one can find from data presented in Figure 14, for the non-preheated samples, Figure 13. Dependencies of TL signal integral vs. dose absorbed for preheated PPG_Tb30 glass dosimeter irradiated with beta radiation (4 µs pulse duration).

Effect of Preheating on Signal Loss in Time
Another important aspect in radiation dosimetry is the stability of signal generated upon irradiation within the detector used. To investigate the signal loss in time, samples of PPG_Tb30 were irradiated using LINAC with dose of 3.65 kGy. After irradiation, samples were stored in dark and preheated right before readout. To compare the obtained results in a quantitative way, Equation (12) was used: where SL is the TL signal loss in time, in %, and N int and N int(to) are the integrated TL signal intensities at a given time and just after irradiation. As one can find from data presented in Figure 14, for the non-preheated samples, nearly 40% of the dosimetric signal is lost during first 24 h after irradiation. Finally, at the end of the monitored period of time, the signal loss after ca. 160 h (7 days) is roughly 63%. For the samples preheated at 120 • C, one observes significant improvement, and the observed signal loss can be estimated for ca. 20-30% during first 160 h. What is interesting, for the samples preheated at highest temperature (200 • C), the observed signal loss is the slightest and during 340 h (14 days) was found to be of ca. 19%. Moreover, when analyzing the shape of signal loss vs. time dependency, one can find that it exhibits a purely linear character. This observation confirms that for the signal loss only recombination of the lower energy traps are responsible, after the emptying of which the main contributors to the TL signal are relatively stable traps activated at higher temperatures. Both, the dose-dependence of the signal and the decrease of the TL signal intensity over time indicate that the most optimal method of reading samples is their preheating at 200 • C, which provides signal linearity for high doses and relative stability of the dosimetric signal over time. Additionally, in case of the readout method based on preheating at 200 • C, one can easily apply the time-related correction of the dose measurement.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 18 of 21 end of the monitored period of time, the signal loss after ca. 160 h (7 days) is roughly 63%. For the samples preheated at 120 °C, one observes significant improvement, and the observed signal loss can be estimated for ca. 20-30% during first 160 h. What is interesting, for the samples preheated at highest temperature (200 °C), the observed signal loss is the slightest and during 340 h (14 days) was found to be of ca. 19%. Moreover, when analyzing the shape of signal loss vs. time dependency, one can find that it exhibits a purely linear character. This observation confirms that for the signal loss only recombination of the lower energy traps are responsible, after the emptying of which the main contributors to the TL signal are relatively stable traps activated at higher temperatures. Both, the dosedependence of the signal and the decrease of the TL signal intensity over time indicate that the most optimal method of reading samples is their preheating at 200 °C, which provides signal linearity for high doses and relative stability of the dosimetric signal over time. Additionally, in case of the readout method based on preheating at 200 °C, one can easily apply the time-related correction of the dose measurement.

Discussion
The alumino-phosphate glasses doped with terbium with different amounts of REO (Rare Earth Oxide) were synthesized by melt-quench method. Glass powders were sieved and pelletized in the form of discs. Samples were irradiated using electron accelerator and irradiation chamber (γ and β-) with radiation dose ranges up to 35 kGy. For samples containing terbium oxide, special shutter had to be used due too significantly higher lightemission efficiency. Addition of terbium oxide, even in amounts of only 0.1%, improves TL intensity significantly. The observed effect is a concentration-dependent phenomenon, as introduction of higher amounts of terbium results in further increasing of thermolumi-

Discussion
The alumino-phosphate glasses doped with terbium with different amounts of REO (Rare Earth Oxide) were synthesized by melt-quench method. Glass powders were sieved and pelletized in the form of discs. Samples were irradiated using electron accelerator and irradiation chamber (γ and β-) with radiation dose ranges up to 35 kGy. For samples containing terbium oxide, special shutter had to be used due too significantly higher lightemission efficiency. Addition of terbium oxide, even in amounts of only 0.1%, improves TL intensity significantly. The observed effect is a concentration-dependent phenomenon, as introduction of higher amounts of terbium results in further increasing of thermoluminescence yield. During the performed research, one established that TL signals from irradiated native and terbium-doped phosphate glasses consist of a combination of different defect types such as first order kinetics (FOK), general order kinetics (GOK), and signals of continuous energy distribution of trapping centers. Verification of the proper localization of particular peaks was based on the consecutive analysis of the experimental data, and the results were obtained from deconvolution analysis using GCA (Glow Curve Analysis) software. Due to the extensive research related to the exact determination of signals deconvolution for glasses studied in this work, one may conclude that this is a very interesting topic that deserves further research and analysis. Additionally, it was found that simple preheating of dosimeters before readout influences dosimetric signals and signal loss and allows for optimization of the readout method of the investigated dosimeters, making the TL signal vs. radiation dose response more predictable. The performed research proved that preheating of the investigated samples at 200 • C improves the linearity of integral TL intensity vs. dose absorbed for doses higher than 15 kGy. Moreover, preheating allows one to minimize relative signal loss in samples over time. An additional advantage of the investigated dosimetric system is its ability to read dosimeters using a basic TL reader without any modifications. All obtained data suggest that the investigated material has high potential for use in high-dose ionizing radiation dosimetry.