Blog

First-principles studies of defect behaviour in bismuth germanate | Scientific Reports

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  12, Article number: 15728 (2022 ) Cite this article Bi4ge3o12 Scintillator

First-principles studies of defect behaviour in bismuth germanate | Scientific Reports

Intrinsic defects are known to greatly affect the structural and electronic properties of scintillators thereby impacting performance when these materials are in operation. In order to overcome this effect, an understanding of the defect process is required for the design of more stable materials. Here we employed density functional theory calculations and the PBE0 hybrid functional to study the structural, electronic,defect process and optical properties of \(\hbox {Bi}_4\hbox {Ge}_3\hbox {O}_{{12}}\) (BGO), a well know material used as scintillator. We examined possible intrinsic defects and calculated their formation energy and their impact on the properties that affect the scintillation process. Furthermore, we investigated the effect and role of rare earth element (REE = Nd, Pr, Ce and Tm) doping on the properties of the BGO system. While the PBE functional underestimated the band gap, the PBE0 was found to adequately describe the electronic properties of the system. Out of all the defects types considered, it was found that \(\hbox {Bi}_{{Ge}}\) antisite is the most favourable defect. Analysis of the effect of this defect on the electronic properties of BGO revealed an opening of ingap states within the valence band. This observation suggests that the \(\hbox {Bi}^{3+}\) could be a charge trapping defect in BGO. We found that the calculated dopant substitution formation energy increases with increase in the size of the dopant and it turns out that the formation of O vacancy is easier in doped systems irrespective of the size of the dopant. We analyzed the optical spectra and noted variations in different regions of the photon energy spectra.

Shiyuan Gao, Hsiao-Yi Chen & Marco Bernardi

Wahib Aggoune, Alberto Eljarrat, … Claudia Draxl

Liang Zhao, Zhifang Shi, … Jinsong Huang

Scintillators are materials that convert high energy rays such as X-rays and \(\gamma\) rays to light. This characteristic is desirable in so many fields. Over the years there has been increased interest in them especially in fast time measurement in nuclear physics1, for precision calorimetry in high-energy physics2 and for positron emission tomography in medical physics3. The interest has resulted in intense efforts channelled towards discovery, research and development of inorganic scintillator materials1,4. To be considered efficient, a scintillator is required to be stable upon exposure to radiation, posses high light yield, fast response, and high efficiency in absorbing radiation. The scintillation properties are closely linked to the structure of the material used. Moreover, the scintillation efficiency is controlled by the presence of defect and crystallographic properties, isotropic propagation of light in scintillation crystals notwithstanding2. Similarly, the scintillation yield, transport and luminescence yield are all dependent on the crystal structure of the scintillator3. Specifically, the energy transfer in the scintillator is a structure sensitive phenomena governed by carrier capture in deep and shallow traps, as well as other radiation-dependent defects5. Defects serve as trap for electrons and holes, interrupting energy transfer in the process. Indeed, previous investigations confirmed the existence of traps in scintillators, although a complete understanding of the energetics of these defects in most materials is still scarce. In view of the strong performance-structure relationship, it is important to understand the defect chemistry of the material, especially those that can be induced when the material is in operation6,7. This is required to improve existing and in design of resilient materials.

One of the most studied material used as scintillator is the Bismuth germanate \(\hbox {Bi}_4\hbox {Ge}_3\hbox {O}_{{12}}\) (BGO)8. First synthesized in 1957, BGO is a cubic crystal with space group I43d possessing eulytine-type structure with \(\hbox {BiO}_6\) octahedron and \(\hbox {GeO}_4\) tetrahedron, an arrangement that can accommodate defects1,9. It has been extensively studied as a scintillator for various applications and found to emit light in the visible region upon exposure to ionizing radiation10. It is also known to possesses high density (7.13 g/\(\hbox {cm}^3\) ) and high light yield (9000 photons/MeV), in addition to impeccable characteristics such as a short radiation length and absence of hygroscopicity11. These characteristics have resulted in the material finding application in tomography scanners and high-precision calorimeters used in the detection of electromagnetic radiation12. This advantages notwithstanding, its wide spread application in fields such as high energy physics is limited due to low response time (about 300 ns). Another issue is the presence of germanium in the system, which raises question about its cost11. Most of the issues encountered in this material is related to its crystal structure13. The arrangement of atoms in BGO is such that a number of charge-trapping sites exist. An understanding of the defect behaviour and formation in the material is required to enhance its performance and minimize deficiencies. Previously, thermo-luminescence experiments have been employed to characterize intrinsic defects in BGO. It was reported that the relative intensities of the glow peaks observed above room temperature depends on radiation dose and the presence of impurities. Certain defect types are suggested as trapping sites with further analysis revealing a range of trapping levels in pristine and doped BGO14. Atomistic simulations employing empirical pair-potential was used to calculate the formation energy of basic defects in BGO. Obtained results supports experiment observation of charge trapping defects in BGO15.

Theoretical methods been have applied successfully to describe various material characteristics related to defect and to calculate defect energetics of materials13,16. Specifically, density functional theory has proved to useful in determining band gaps and defect properties of scintillators17. Studies of this nature consider deviation from stoichiometry resulting from formation of intrinsic defects, whose presence determines the stability of the material in operation. The incorporation of dopant ions into perovskite and similar structures in a wide range of concentrations has been reported to improve properties and applicability of materials18,19. Specifically, the introduction of rare-earth elements (REE) has received immense attention due to its ability to modify electronic properties and luminescence in scintillator materials20. Moreover, the doping of BGO attracted attention due to the ability of its photons to interact with the material effectively and combine to form new photons with doubled energy and frequency21. Indeed, REE are interesting dopants for enhancing the properties of BGO. Among REE ions, \(\hbox {Pr}^{3+}\) , \(\hbox {Nd}^{3+}\) , \(\hbox {Tm}^{3+}\) and \(\hbox {Ce}^{3+}\) have received attention as a result of offering remarkable activator ion for luminescence12,19. Jazmati and coworkers22 investigated BGO: Ce samples implanted at linear no-threshold model at 77 K with He ions for manufacturing waveguides. They observed a phase change in the BGO, modifying its cubic structure to an anisotropic guide layer generated from the ’stress’ of the He beam deployment and, at the same time modifying the optical activity. Besides, Nd doped BGO demonstrates the properties that allow its use in the construction of solid-state lasers23. The advantages reported for these resultant materials notwithstanding, their practical applicability has been hindered by lack of detailed information about their microstructure and the position of the dopant atom in the system. The choice of the REE dopant employed for our investigation is guided by experimental findings. Different REE dopants have been reported to improve scintillation performance12,22. For instance, it has been show that the radiation resistance of BGO crystal was improved by Eu doping leading to faster induced absorption recovery24. In the same vein, Ce doping has been found to lead to occurrence of thermo-luminescence (TSL) peaks around room temperature (RT)25. Similarly, Nd, Tm and Ce are attractive dopant as they have been found to posses emission lines due to 4f–4f transmission from visible to near-infrared wavelength, hence are known as luminescence centers26.

(a) Crystal structure of optimized BGO unit cell and (b) \(2 \times 2\times\) 2 supercell used for defect simulations.

Previously, density functional theory (DFT)13,27,28,29 has been employed to study the electronic and optical properties of BGO. Also, attempts have been made to study the doping of BGO with rare earth elements in order to improve its applicability using DFT12. These investigations employed different flavours of the DFT within Generalized gradient approximation (GGA). It was reported that the major part of incoming energy is absorbed by oxygen (O) P electrons and transferred to bismuth (Bi) p states. All of these investigations reported band gaps that are below the experimental reported values. In order to accurately characterize the defect behaviour and optical properties of a system, it is important to obtain electronic properties that are close to experiment30. This work employs the hybrid PBE0 functional to study BGO. The hybrid functional, specifically the PBE0 has been used successfully to study systems where GGA has been found to underestimate band gaps31. In this work, we investigate the defect process in BGO and predict the stability of the different defect type, the effect of the prominent defect on factors affecting scintillation using density functional theory. Although this work focuses on BGO, inferences drawn from this investigation will aid understanding other materials used for similar applications. Moreover, it is expected that understanding of these kind of defects will help to optimize the efficiency of scintillators.

The BGO system crystallizes in the cubic symmetry, comprising of a regular arrangement of \(\hbox {GeO}_4\) tetrahedra sharing vertices with distorted \(\hbox {BiO}_6\) octahedra, see Fig. 1a. The primitive unit cell contains two formula unit of 38 atoms13. Rodriguez and coworkers28 reported the crystal structure of BGO with experimental lattice parameter a = b = c = 10.513 Å. We begin our investigation by calculating the lattice parameter of pristine BGO and obtained a lattice parameter of 10.6 Å which is in good agreement with values reported from previous theoretical work12 and experimentally28,32. The BGO structure is such that two different Bi–O bonds exist due to the distortion of the the \(\hbox {BiO}_6\) octahedron. The bond lengths for the pristine system are presented in Table 1. Observation reveals an agreement with values reported by experiment32. After validating our model for the calculation of the structural properties of defect-free BGO, we proceeded to calculate the electronic structure of the pristine system. The density of states in Fig. 2 offers insight into the chemical bonding in pristine BGO. Figure 2a presents the GGA calculated density of states, where we obtained a band gap of 3.4 eV comprising of a O-p states dominated valence band maximum (VBM) and a conduction band that comprises of hybridization of Bi-p and O-p states. Note that the calculated band gap is smaller than the experimental band gap of 4.58–4.7 eV33,34. Previous DFT investigations using GGA12,28 reported a band gap of 3.5 eV. Similar calculations using full-potential linearized augmented plane wave, taking into account spin-orbit interaction yielded similar band gap13. All previous work were unable to reproduce the experimental band gap. Correct prediction of the electronic properties especially fundamental properties of the gap, its width and states is important in determining the defect energetics30, optical properties and characterization of scintillation process13. In view of this, there is a need to consider another functional in a bid to improve the result.The DFT-PBE0 method with 12.5% exchange was employed and the obtained band gap of 4.6 eV is in better agreement with experiment than the one obtained with DFT-PBE. Figure 2b presents PBE0 calculated DOS for pristine BGO. The same features obtained for DFT-PBE are reproduced here with the only difference resulting from increase in the size of the band gap.The band structures (Fig. 2c) calculated along calculated along \(\Gamma\) \(\rightarrow\) H \(\rightarrow\) N \(\rightarrow\) \(\Gamma\) \(\rightarrow\) P \(\rightarrow\) H|P high symmetry path of the Brillouin zone further supports the information obtained from the density of state plot shown in Fig. 2b. It is worth noting that our calculated band gap using the PBE0 significantly improves the band gap previously obtained using PBE-GGA functionals (see Table 1).

The electronic band gap of a material is a valuable feature that provides a deep understanding of its electronic, defect and optical properties. From our investigation, it was found that, the PBE0 approximation improves greatly the value of the band gap energy. Infact, PBE0 approximation improves, significantly, the calculated gap value better than the conventional GGA approximation. Llalic and coworkers13 have previously carried out first-principles calculations, including spin orbit coupling on BGO, while certain features were found to be improved in comparison to standard DFT, the band gap was underestimated. Using PBE0, key features of the electronic states are reproduced while a much improved band gap is also obtained. An accurate prediction of the band gap is important for the meaningful characterization of defect energetics35,36. Moreover, several DFT calculations36,37,38 on Bi-containing oxides have been successfully carried out, including investigations involving defect without the employing spin-orbit coupling.

Total and partial density of states BGO calculated using (a) DFT-GGA (b) DFT-PBE0 and (c) band structure of pristine BGO calculated with DFT-PBE0.

Total and partial density of states of (a) pristine BGO and BGO (b) antisite defect (c) O vacancy defect (d) Bi vacancy defect (e) Ge vacancy defect.

Several important optical and luminescent characteristics are influenced by energy levels formed within band width of scintillators31. These levels are due to activator ions, impurities and point defects. Defects in materials are important as they control the physical, thermodynamic and electronic properties31,39. Intrinsic defects are disorders that can be thermally generated and not from doping or reaction with the environment31. They have been studied previously in similar systems and for similar applications31,40. For this investigation, we considered only isolated single defects for simplicity and excluded defect clustering. Similar approach was adopted in studies of defect for similar applications40,41.

The intrinsic defects in BGO include Bi, O, and Ge vacancy/interstitial. We modelled intrinsic defects in the BGO system by adding/removing an atom (of Bi, Ge, or O) to/from a \(2 \times 2\times\) 2 supercell (Fig. 1b). The formation energy of the different intrinsic defects considered are calculated thus:

\(E_{defective}\) and \(E_{Pristine}\) are the total energies of defective and defect-free supercells, respectively. \(\mu\) is the chemical potential of the atoms considered. The sign in equation is \(+\) for vacancy and − for interstitial. A combination of one vacancy and one interstitial defect makes a Frenkel pair. To calculate, the Frenkel energy, we first investigate formation energy of the different possible O vacancy sites and find the most favourable site, considering O vacancies in the \(\hbox {BiO}_6\) and \(\hbox {GeO}_4\) coordination. Generally, we find that the O vacancy formation energies are slightly lower for Bi coordinated O atoms than Ge coordinated ones. We obtained O vacancy formation energy of 2.32 and 2.30 eV for Ge and Bi coordinated O atoms respectively. The same approach used in determination of stable site for O vacancy formation is adopted in finding the most favorable interstitial site. The O interstitial formation energy is calculated to vary between 1.47 and 1.61 eV for the possible scenarios considered. From the most energetically favorable configurations for the O vacancy and interstitial, we obtained the O Frenkel energy for BGO to be between 0.83 eV (2.30–1.47 eV) and 0.71 eV ( 2.32–1.61 eV) respectively for the Ge and Bi sublattices. From the calculation of the O defect formation energies, we find that the \(\hbox {GeO}_4\) sublattice is stable, hence most defects are associated with the \(\hbox {BiO}_6\) sublattice. In Table 2, we list the calculated formation energies of O, Ge and Bi vacancies and interstitials.

Another type of defect considered is the antisite defect. Antisite defects are a common defect type in many oxide based compounds and have been reported experimentally and theoretically35,42. To generate the antisite defect, we placed Ge in a site that is originally occupied by Bi and vice versa. Similar approach has been adopted in study of antisite defects in oxides35,43. The antisite defect formation energy \(\hbox {E}_{antisite_formation}\) is defined as the energy difference between energy of configuration with an antisite defect (\(\hbox {E}_{{antisite}}\) ) and the energy of the pristine system43.

We obtained the energy required to form the \(\hbox {Bi}_{{Ge}}\) and \(\hbox {Ge}_{{Bi}}\) antisites as − 0.48 eV and 4.30 eV respectively. A comparison of the Frenkel and antisite defect energetics shows that the \(\hbox {Bi}_Ge\) antisite defect are much easier to form than the O Frenkel. This observation suggests that the \(\hbox {Bi}^{3+}\) is likely to be charge trap site and the main optically active constituent of the BGO.

The electronic structure of scintillators is a crucial factor in the luminescence properties of these components. In order to achieve the best possible efficiency, the band gap needs to be narrowed13. To investigate the origin of the electronic structure modification due to intrinsic defects, we analyse the density of states as shown in Fig. 3 and compare with some of the prominent defects. As mentioned earlier, our calculated band gap for pristine BGO is 4.6 eV (shown again in Fig. 3a). Our analysis of the electronic structure was carried out on the most favourable defect types. Here we consider the different vacancies and antisite defect. The presence of vacancy introduces defects states within the band. Figure 3b–e shows the densities of state for the BGO with antisite defect, O vacancy, Bi vacancy and Ge vacancy respectively. The plots show that the position of the defect state induced depends on the type of defect. Generally, the overall shape of the DOS and composition of the bands are unaltered for the vacancy defects. However, for the case of system with antisite defect, there is a downward shift in the conduction band into the band gap when compared with the pristine case. Similar downward shift was observed for the valence band culminating in a split into discrete bands. The states introduced are mainly as a result of Bi p states. This suggests that charge carrier thermalization can be suppressed thereby leading to intra band luminescence of carriers.

(a) Change in volume of BGO upon doping with REE (b) Formation energy of doped REE systems.

The trivalent Bi cation has a suitable ionic size to accommodate the REE ion and has been found to influence properties of scintillators44. This makes them suitable dopants for the BGO system22. To substitute the REE atom at the Bi site, we considered all the non-equivalent Bi sites and proceeded with our calculation with the the most stable site for doping. For our investigation, we considered REE = Ce, Pr, Nd, and Tm as dopant elements;all considered in previous experimental studies22,24. The REE atoms are considered to replace Bi atom culminating in a REE/Bi ratio of 1/64 in the supercell. For the considered dopants, the size of the ionic radii decreases from Ce (1.01 Å) to Tm (0.88 Å) resulting in a decrease in REE-O bond length. Overall, the doping of BGO with REE affects the structural properties of BGO minimally. Figure 4a shows the percentage change in volume of relaxed structure of the BGO:REE. We find that the volume is not hugely changed upon doping with REE, this is because the ionic radii of the Bi and the REE ions are similar and the concentration of REE introduced is small. Our simulation revealed a maximum change of 1.62% (Tm). The minimal change in volume is found for BGO:Ce. This is expected as the ionic radii of Ce (1.01Å) is similar to that of Bi (1.03Å)45. Overall, the volume of the system decreases as the ionic radii of REE decreases (REE doping leads to reduction in volume). Our investigation proceeded with the determination of the stability of the dopant systems by evaluating the formation energy of the doping with REE. We calculated the formation energies for the different dopant systems using16

where \(\hbox {E}_{{formation}}\) , \(\hbox {E}_{{REE}}\) , \(\hbox {E}_{{pristine}}\) are the formation energy of substitution of REE, energy of system with REE dopant and energy of pristine system respectively. \(\mu _{Bi}\) and\(\mu _{REE}\) are the chemical potential of the Bi, REE and oxygen atoms respectively. The value calculated quantifies the ease of substitution of the dopant element in the BGO system. Figure 4b shows the calculated dopant substitution formation energy as a function of the different dopant elements. We observed that the value of \(\hbox {E}_{{formation}}\) gradually reduces as the size of the REE substituted. However, the Nd doped system does not follow this trend. To shed light into the reason for this deviation, the REE-O bond length is examined. We find that although the REE–O length increases as the ionic radii45 of REE increases, Nd–O bond length is similar to Tm–O bond length. Jazmati et al.22 studied the role of rare earth on properties of BGO and found unique properties in Nd doped BGO. It was found to modify its cubic structure of BGO to an anisotropic guide layer thereby exhibiting non linear optical behaviour.

O vacancy formation energy of pristine and REE doped BGO.

The existence of oxygen vacancy in BGO is inherent and its formation in pristine BGO has been discussed earlier. Moreover, the local distortion caused by the REE dopant (as a result of difference in size) may affect the stability of surrounding O atoms16. The effect of location distortion induced as result of REE doping has been reported in oxides16,46. Hence it is important to quantify the stability of the O vacancy in the system. Upon substitution of the REE dopant, the O vacancy formation energy is calculated using Eq. (5)

where \(\hbox {E}_{V^1_o}\) , \(\hbox {E}_{{REE}}\) and \(\hbox {E}_{(REE+Vac)}\) are the energy of formation of O vacancy in a system doped with REE, energy of system with REE dopant at the Bi site only and energy of system with REE dopant with O vacancy, respectively and \(\mu _{O}\) is the chemical potential of O. The O vacancy formation energy as a function of dopant is presented in Fig. 5. We considered O sites close to the REE dopant (see Fig. 6). Our investigation shows that all of the considered dopant systems are less stable in the presence of O vacancy than the pristine BGO system with the Ce-doped system being the most stable of all the considered REE dopants. It is worth noting that while the O vacancy formation energies obtained for Tm, Nd and Ce follow a trend that mimicks the ionic radii of the dopant (see Fig. 4a), Pr does not follow this trend. O vacancy formation energy is formed more easily in Pr doped system compared to the other dopants. Overall, the REE with largest ionic radii posses the highest O vacancy formation energy, see Fig. 5. The very small difference in O vacancy formation energy in the dopant system suggests similar behavior in the system, and shows that even with low concentration of dopant, the stability of the BGO system can be altered in the presence of defect like O vacancy.

Arrangement of atom around REE atom (a) before O vacancy (b) after O vacancy. The yellow ball depicts Pr atom and marked O is the atom removed.

In this section, we discuss the optical properties of BGO within the photon energy 0–30 eV, calculated using the PBE0 functional. Here, we calculate the absorption coefficient, reflectivity, extinction coefficient and refractive index. The ability of a material to store electrical charges is often quantified by its dielectric constant (\(\epsilon\) ). It is a significant optical property of a material. It comprises of real and imaginary part as described above. The imaginary part is directly proportional to absorption spectrum. The absorption spectra is shown Fig 7a. Analysis of the absorption spectra characterized the highest intensity to be within 4–13 eV with the peak at 8 eV. Above 13 eV, we notice a medium intensity with a smaller peak at 15.5 eV. Above 20 eV there is negligible intensity. Where there is zero absorption intensity, it implies that at the corresponding energy there is absence of dispersion resulting in maximum absorption. The nature of reflectance of incident radiation on the BGO is described by its reflectivity. The reflective spectra is shown in Fig. 7b. Similar to the absorption spectra, it is characterized by three regimes with the peak attaining maximum at about 12 eV. The obtained spectrum reproduced features observed experimentally47. Not only is the shape of the spectrum similar, the peaks are situated at comparable energies. The extinction properties are shown in Fig. 7c, it elucidates the absorption losses at particular ranges of incident electromagnetic spectrum. From our calculated spectrum, it can be seen that the extinction coefficient increases gradually and is maximum at 6.4 eV from which its value declines rapidly further. We obtained a refractive index of 2.6 at 3.5 eV, see Fig. 7d. This value is in agreement refractive index obtained in experiment47.

The calculated of (a) absorption coefficient (b) reflectivity, (c) extinction coefficient and (d) refractive index of pristine BGO. The lower panel presents (d) absorption coefficient (e) reflectivity, (f) extinction coefficient and (g) refractive index of Ce-doped BGO.

We present the absorption, reflective, extinction and refractive spectra of the most stable dopant system (Ce-doped BGO) in Fig. 7e–h. A comparison of the spectra with those obtained for pristine BGO shows some similarities. For example, the peaks for the four spectra are of the same order of magnitude and obtained at almost the same energy. The similarities notwithstanding, we observed contrast in certain features. For example, while the absportion spectra has no intensity above 20 eV for pristine BGO, we observed intensity upto 25 eV for Ce-doped BGO. Similar observation is made for reflectivity. The fact that upon doping the BGO system with very low concentration of Pr we observe contributions to the optical adsorption suggests that the low light yield and slow response which has plagued the applicability of the BGO crystals as scintillators for positron electron tomograhy can be improved by defect engineering as well as improved crystal growth.

We have investigated the energetics of the formation and electronic properties of isolated defects in BGO using first-principles method. The investigation is essential to identify routes to improve scintillator performance. The calculation employed the hybrid PBE0 functional with 12.5% Hatree Fock exchange component to reproduce experimental band gap since PBE functional has been found to underestimate this very important property. \(\hbox {B}_{{Ge}}\) antisite defect was predicted to be the most favorable defect of all considered defect types and introduces states within the band gap, introduces states within the valence band while shifting the bands relative to the states obtained for the pristine BGO, confirming the presence of trap sites in BGO. Furthermore, our investigation examined different rare earth dopants and their formation energies. We found O vacancy formation to be easily formed in REE doped system than in the pristine system. We also investigated the optical properties using approach that best reproduce the band gap of pristine BGO and noted variations in different regions of the photon energy spectra. Our investigation shed light on the role of defects and how using cation doping can be a route to control the stability of BGO.

All calculations were carried out using the ab initio density functional theory using the projector augmented wave method of the Vienna Ab initio Simulation Package48. A plane wave basis set with a cut-off energy of 520 eV and a \(4 \times 4\times\) 4 Monkhorst-Pack generated sets of k-points were used for optimizing \(2 \times 2\times\) 2 supercell of BGO. We tested the convergence of the cutoff energy and K-points and found that increasing both produced negligible difference in the calculated results. (see Figs. S1 and S2 in Supplementary information. Both the lattice parameter and atomic positions relaxed with the energies and the forces of each ion were converged within \(1.0 \times 10^{-5}\)  eV/atom and 0.01 eV/Å, respectively. In order to reproduce the electronic properties as reported by previous experiment, hybrid functional, PBE0 (a scheme that combines a fraction of Hartree-Fock exchange with PBE exchange) calculations was carried out after Perdew-Burke-Ernzerhof (PBE)49,50 was unable to correctly predict the band gap of BGO. Defects are simulated by removing an O atom from a \(2 \times 2\times\) 2 supercell constructed from the cubic structure. The size of cell used for the simulation of O vacancy is large enough to prevent defect-defect interactions.A minimal deviation in formation energy (about 1 meV) was found when a larger \(3 \times 3\times\) 3 supercell was used.

Furthermore, we calculated the optical properties within the random phase approximation (RPA) by evaluating the complex dielectric tensor \(\epsilon\) . The dielectric tensor comprises of real and imaginary part, with the imaginary part being directly proportional to the optical spectrum of the system. Also, the dielectric function can be expressed in terms of refractive index n, extinction coefficient k as follows13:

Knowledge of the energy dependence of this component will help to quantify the response of the material to incident energy and the subsequent decay of that light propagating through the absorbing medium51.

The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.

Weber, M. J. Inorganic scintillators: Today and tomorrow. J. Lumin. 100, 35–45 (2002).

Lecoq, P., Annenkov, A., Gektin, A., Korzhik, M. & Pedrini, C. Influence of the crystal structure defects on scintillation properties. Inorg. Scintill. Detect. Syst. Phys. Princ. Cryst. Eng. 123–173 (2006).

Van Eijk, C. Inorganic scintillators for the next generation of neutron beam facilities. In Proceedings of the Fifth International Conference on Inorganic Scintillators and Their Applications, 22–32 (Moscow State University, 2000).

McDonald, K. A., McDonald, M. R., Bailey, M. N. & Schweitzer, G. K. Parametric study on the production of the gagg: Ce and lso: Ce multicomponent oxide scintillator materials through use of a planetary ball mill. Dalton Trans. 47, 13190–13203 (2018).

Article  CAS  PubMed  Google Scholar 

Robbins, D. On predicting the maximum efficiency of phosphor systems excited by ionizing radiation. J. Electrochem. Soc. 127, 2694 (1980).

Article  ADS  CAS  Google Scholar 

Blasse, G. Luminescent materials: Is there still news?. J. Alloy. Compd. 225, 529–533 (1995).

Stanek, C., McClellan, K., Levy, M. & Grimes, R. Defect behavior in rare earth re al o 3 scintillators. J. Appl. Phys. 99, 113518 (2006).

Article  ADS  CAS  Google Scholar 

Casey, M. & Nutt, R. A multicrystal two dimensional bgo detector system for positron emission tomography. IEEE Trans. Nucl. Sci. 33, 460–463 (1986).

Tekin, H. et al. Newly developed bgo glasses: synthesis, optical and nuclear radiation shielding properties. Ceram. Int. 46, 11861–11873 (2020).

Weber, M. J. & Monchamp, R. R. Luminescence of bi4 ge3 o12: Spectral and decay properties. J. Appl. Phys. 44, 5495–5499 (1973).

Article  ADS  CAS  Google Scholar 

Kobayashi, M., Ishii, M., Harada, K. & Yamaga, I. Bismuth silicate bi4si3o12, a faster scintillator than bismuth germanate bi4ge3o12. Nucl. Instrum. Methods Phys. Res. Sect. A 372, 45–50 (1996).

Article  ADS  CAS  Google Scholar 

Chen, F. et al. Structure and luminescence properties of a nd 3+ doped bi 4 ge 3 o 12 scintillation crystal: New insights from a comprehensive study. J. Mater. Chem. C 5, 3079–3087 (2017).

Lalic, M. & Souza, S. The first-principles study of electronic and optical properties of bgo and bso scintillators. Opt. Mater. 30, 1189–1192 (2008).

Article  ADS  CAS  Google Scholar 

Raymond , S. , Luff , B. , Townsend , P. , Feng , X. & Hu , G. Thermoluminescence spectra of doped bi4ge3o12 .Radiation.Meas.23, 195–202 (1994).

Jackson, R. & Valerio, M. Computer modelling of radiation generated defects in bgo scintillators. Nucl. Instrum. Methods Phys. Res. Sect. B 218, 145–147 (2004).

Article  ADS  CAS  Google Scholar 

Zhu, J., Gu, M., Jia, L. & Song, G. Structural properties of lu2sio5 doped with rare-earth elements. Mater. Lett. 256, 126410 (2019).

Liu, B. et al. First-principles study of oxygen vacancies in lu2sio5. J. Phys. Condens. Matter 19, 436215 (2007).

Article  ADS  CAS  Google Scholar 

Zulueta, Y., Lim, T. & Dawson, J. Defect clustering in rare-earth-doped batio3 and srtio3 and its influence on dopant incorporation. J. Phys. Chem. C 121, 23642–23648 (2017).

Gökçe, M. Development of eu3+ doped bismuth germanate glasses for red laser applications. J. Non-Cryst. Solids 505, 272–278 (2019).

Article  ADS  CAS  Google Scholar 

He, L. et al. Insight into the characteristics of 4f-related electronic transitions for rare-earth-doped klus2 luminescent materials through first-principles calculation. J. Phys. Chem. C 124, 932–938 (2019).

Cates, J. W. & Levin, C. S. Electronics method to advance the coincidence time resolution with bismuth germanate. Phys. Med. Biol. 64, 175016 (2019).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Jazmati, A. & Townsend, P. Photoluminescence from re doped bgo waveguides. Nucl. Instrum. Methods Phys. Res. Sect. B 166, 597–601 (2000).

Feng, X.-Q. et al. Growth, laser and magneto-optic properties of nd-doped bi4ge3o12 crystals. Mater. Sci. Eng. B 23, 83–87 (1994).

Arslanlar , YT , Kotan , Z. , Kibar , R. , Canımoğlu , A. & Can , N. Rare earth photoluminescence in bismuth-germanate crystals .Spectrosc.Lett.Rev. 46, 590–596 (2013).

Article  ADS  CAS  Google Scholar 

Mares, J. et al. Time development of scintillating response in ce-or pr-doped crystals. Phys. Status Solidi c 4, 996–999 (2007).

Article  ADS  CAS  Google Scholar 

Frindell, K. L. et al. Visible and near-ir luminescence via energy transfer in rare earth doped mesoporous titania thin films with nanocrystalline walls. J. Solid State Chem. 172, 81–88 (2003).

Article  ADS  CAS  Google Scholar 

Sousa, O., Lima, A. & Lalic, M. New insights into the electronic and optical properties of the bi4m3o12 (m= si or ge) scintillators. Opt. Mater. 73, 642–646 (2017).

Article  ADS  CAS  Google Scholar 

Rodriguez, J. R. et al. Bismuth germanate (bi 4 ge 3 o 12), a promising high-capacity lithium-ion battery anode. Chem. Commun. 54, 11483–11486 (2018).

Isik, M., Surucu, G., Gencer, A. & Gasanly, N. M. First principles study of bi12geo20: Electronic, optical and thermodynamic characterizations. Mater. Today Commun. 27, 102299 (2021).

Akande, S. O., Chroneos, A., Vasilopoulou, M., Kennou, S. & Schwingenschlögl, U. Vacancy formation in moo 3: Hybrid density functional theory and photoemission experiments. J. Mater. Chem. C 4, 9526–9531 (2016).

Ning, L., Cheng, W., Zhou, C., Duan, C. & Zhang, Y. Energetic, optical, and electronic properties of intrinsic electron-trapping defects in yalo3: A hybrid dft study. J. Phys. Chem. C 118, 19940–19947 (2014).

Milenov, T., Rafailov, P., Petrova, R., Kargin, Y. F. & Gospodinov, M. X-ray diffraction study of a bi4ge3o12 crystal. Mater. Sci. Eng. B 138, 35–40 (2007).

Santana, G. C., De Mello, A. C. S., Valerio, M. E. & Macedo, Z. S. Scintillating properties of pure and doped bgo ceramics. J. Mater. Sci. 42, 2231–2235 (2007).

Article  ADS  CAS  Google Scholar 

Polosan, S., Galca, A. & Secu, M. Band-gap correlations in bi4ge3o12 amorphous and glass-ceramic materials. Solid State Sci. 13, 49–53 (2011).

Article  ADS  CAS  Google Scholar 

Akande, S. O., Chroneos, A. & Schwingenschlögl, U. O vacancy formation in (pr/gd) baco 2 o 5.5 and the role of antisite defects. Phys. Chem. Chem. Phys. 19, 11455–11459 (2017).

Li, M., Li, J., Guo, C. & Zhang, L. Doping bismuth oxyhalides with indium: A dft calculations on tuning electronic and optical properties. Chem. Phys. Lett. 705, 31–37 (2018).

Article  ADS  CAS  Google Scholar 

Kweon, K. E. & Hwang, G. S. Hybrid density functional study of the structural, bonding, and electronic properties of bismuth vanadate. Phys. Rev. B 86, 165209 (2012).

Article  ADS  CAS  Google Scholar 

Iori, F. et al. Bismuth iron garnet: Ab initio study of electronic properties. Phys. Rev. B 100, 245150 (2019).

Article  ADS  CAS  Google Scholar 

Kuganathan , N. & Chroneos , A. Defects and dopants in cafesi2o6: Classical and dft simulations.Energies 13, 1285 (2020).

Liu, X.-Y., Pilania, G., Talapatra, A. A., Stanek, C. R. & Uberuaga, B. P. Band-edge engineering to eliminate radiation-induced defect states in perovskite scintillators. ACS Appl. Mater. Interfaces 12, 46296–46305 (2020).

Article  CAS  PubMed  Google Scholar 

Hu, C. et al. Antisite defects in nonstoichiometric lu3al5o12: Ce ceramic scintillators. Phys. Status Solidi (b) 252, 1993–1999 (2015).

Article  ADS  CAS  Google Scholar 

Kempaiah Devaraju, M., Duc Truong, Q., Hyodo, H., Sasaki, Y. & Honma, I. Synthesis, characterization and observation of antisite defects in linipo4 nanomaterials. Sci. Rep. 5, 1–8 (2015).

Dey, U., Chatterjee, S. & Taraphder, A. Antisite-disorder engineering in la-based oxide heterostructures via oxygen vacancy control. Phys. Chem. Chem. Phys. 20, 17871–17880 (2018).

Article  CAS  PubMed  Google Scholar 

Zhydachevskyy, Y. et al. Band gap engineering and trap depths of intrinsic point defects in ralo3 (r= y, la, gd, yb, lu) perovskites. J. Phys. Chem. C 125, 26698–26710 (2021).

Shannon, R. D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides. Acta Crystallogr. Sect. A Cryst. Phys. Diffr. Theor. Gen. Crystallogr. 32, 751–767 (1976).

El Sayed, M. E., Naji, S., Murshed, M. N. & Samir, A. Cation substitution for tunable electronic, optical and scintillation properties of pb1-xcaxwo4 materials: A density functional theory study. Results Phys. 30, 104826 (2021).

Antonangeli, F., Zema, N., Piacentini, M. & Grassano, U. Reflectivity of bismuth germanate.Phys.Rev. B 37, 9036 (1988).

Article  ADS  CAS  Google Scholar 

Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169 (1996).

Article  ADS  CAS  Google Scholar 

Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).

Article  ADS  CAS  PubMed  Google Scholar 

Ernzerhof, M. & Scuseria, GE Assessment of the perdew-burke-ernzerhof exchange-correlation functional.J. Chem.Phys.110, 5029–5036 (1999).

Article  ADS  CAS  Google Scholar 

Dar, S. A., Sharma, R., Srivastava, V. & Sakalle, U. K. Investigation on the electronic structure, optical, elastic, mechanical, thermodynamic and thermoelectric properties of wide band gap semiconductor double perovskite ba 2 intao 6. RSC Adv. 9, 9522–9532 (2019).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

The advanced computing facility of Texas A & M University at Qatar is used for all calculations. This work is supported by an internal grant from TAMU Qatar and a joint agreement from Qatar Computing Research Institute. Open Access funding provided by the Qatar National Library.

Texas A & M University at Qatar, PO Box 23874, Doha, Qatar

Salawu Omotayo Akande - I'm Not Afraid (Official Music Video)

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

S.O. performed all the calculations. O.B. supervised overall and provided funding. All authors took part in discussion, analysis of result and in finalizing the manuscript.

Correspondence to Salawu Omotayo Akande.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Akande, S.O., Bouhali, O. First-principles studies of defect behaviour in bismuth germanate. Sci Rep 12, 15728 (2022). https://doi.org/10.1038/s41598-022-18586-x

DOI: https://doi.org/10.1038/s41598-022-18586-x

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

First-principles studies of defect behaviour in bismuth germanate | Scientific Reports

Gagg Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.