Next Article in Journal
Hydroxyapatite Double Substituted with Zinc and Silicate Ions: Possibility of Mechanochemical Synthesis and In Vitro Properties
Next Article in Special Issue
The d Orbital Multi Pattern Occupancy in a Partially Filled d Shell: The KFeF3 Perovskite as a Test Case
Previous Article in Journal
Effects of Oxide Fragments on Microstructure and Mechanical Properties of AA6061 Aluminum Alloy Tube Fabricated by Thermomechanical Consolidation of Machining Chips
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Difference in Structure and Electronic Properties of Oxygen Vacancies in α-Quartz and α-Cristobalite Phases of SiO2

by
Katherine L. Milton
1,*,
Thomas R. Durrant
2,
Teofilo Cobos Freire
2 and
Alexander L. Shluger
1,2,*
1
London Center for Nanotechnology, University College London, Gower Street, London WC1E 6BT, UK
2
Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK
*
Authors to whom correspondence should be addressed.
Materials 2023, 16(4), 1382; https://doi.org/10.3390/ma16041382
Submission received: 3 January 2023 / Revised: 30 January 2023 / Accepted: 1 February 2023 / Published: 7 February 2023

Abstract

:
α -cristobalite ( α -C) is a polymorph of silica, mainly found in space exploration and geochemistry research. Due to similar densities, α -C is often used as a proxy for amorphous SiO 2 , particularly in computer simulations of SiO 2 surfaces and interfaces. However, little is known about the properties of α -C and its basic oxygen defects. Using density functional theory (DFT) simulations we provide a comprehensive report on the properties of perfect structure and oxygen vacancies in α -C. The calculated properties of α -C are compared with those of the better-characterized α -quartz ( α -Q). Our results demonstrated that the positively charged O vacancy in α -C is most stable in the dimer configuration, in contrast to α -Q, which favors the puckered configuration. A back-projected configuration was also predicted in both polymorphs. We calculated the optical transition energies and isotropic hyperfine constants for O vacancies in both α -Q and α -C, and compared our findings with the results of previous studies and experiments. This work, thus, offers one of the first in-depth investigations of the properties of oxygen vacancies in α -C.

1. Introduction

Oxygen vacancies are among the most basic defects in oxide crystals, and their properties affect and control the performance of oxides in technology. Defects in SiO 2 have long been the focus of many-electron calculations using quantum chemistry [1,2,3] and periodic DFT [4,5,6,7,8], due to their importance in numerous applications, including microelectronics, catalysis, fiber optics, and dosimetry. These calculations demonstrated that oxygen vacancies in SiO 2 distinguish themselves from other oxides in at least two accounts: (1) the number of different configurations they adopt, and (2) the significant distortions they are capable of exerting on the surrounding lattice. Both of these properties are due to the extraordinary flexibility of the SiO 2 network. Oxygen vacancies in silica are best studied experimentally in α -quartz ( α -Q) and amorphous SiO 2 (a-SiO 2 ), where they have been experimentally shown to be stable under many conditions and have at least two (+1 and 0) accessible charge states [9]. We note that, quite loosely, we denote here many types of silica glass, as well as amorphous silica films, as a-SiO 2 . Although we are aware of many subtle differences, these do not affect our discussion.
On the other hand, theoretical calculations predicted the existence of five charge states of the O vacancy in α -Q and a-SiO 2 (with defect charges ranging + 2 0 2 ) [7,10,11,12,13], as well as more than eight different geometric configurations for the +1 charge state vacancy in a-SiO 2 (see, e.g., [3,8]). When transforming between these charge states and configurations, O vacancies in SiO 2 undergo structural changes with atomic displacements often exceeding 0.5 Å [3], in excess of typical displacements in other oxides (0.1–0.3 Å [14]). The abundance of predicted configurations of positively charged vacancies in a-SiO 2 results from different short and medium range order environments surrounding oxygen sites, due to disorder.
Quite extraordinarily, SiO 2 is one of those materials where some of the experimentally suggested defect configurations have never been verified theoretically (e.g., the so called E δ center in a-SiO 2 [15]) and a large number of theoretically predicted defect configurations and charge states have not yet been identified experimentally. This huge gulf between theory and experiment is, in part, caused by high energies and the instability of some defect configurations, and by difficulties in experimental measurements and interpretation. The uncertainty in the verification of theoretical predictions undermines current attempts to employ machine learning to predict defect structures in solids (see e.g., [16]) and, particularly, in amorphous SiO 2 [17].
One of the longest controversies concerns the existence of the so-called “dimer” configuration of the +1 charge state of O vacancy in SiO 2 . This simple configuration was predicted by most theoretical calculations of α -Q and a-SiO 2 but never observed experimentally, despite numerous attempts (see an overview in ref. [15]). This is likely due to the small barrier for its transformation into the more complex dangling bond configuration of O vacancy, also known as the “puckered” configuration (see details in the next section). Since the existence and stability of both configurations is strongly affected by the short and medium range environment of O sites in SiO 2 , one could gain additional insights from investigating these defects in other crystalline SiO 2 polymorphs. However, surprisingly little is known about the properties of O vacancies in other stable silica polymorphs, such as α -cristobalite ( α -C) or tridymite. Samples of these polymorphs are much less abundant and often much more contaminated with impurities than those of α -Q [18,19]. Therefore, theoretical calculations are often the only way to shed light on their properties and provide the information required to identify defects in these polymorphs.
Here, we focused on calculating properties of O vacancies in α -C, which is a high-temperature and low-pressure silica polymorph. It has a less dense structure than α -Q (see Figure 1) with larger rings close to a-SiO 2 structures. There are no theoretical studies of the O vacancies in α -C. By comparing our results for α -C with those for α -Q and a-SiO 2 we hoped to better understand the vacancy structure and guide experimental studies of their identification in α -C. In particular, we were interested in comparing the stable configurations of positively charged O vacancies in both crystals.
Our results demonstrated that the positively charged O vacancy in α -C is most stable in the dimer configuration, in contrast to α -Q, which favors the puckered configuration. A back-projected configuration was also predicted in both polymorphs. We calculated the optical transition energies and isotropic hyperfine constants for O vacancies in both α -Q and α -C, and compared our findings with the results of previous studies using molecular clusters [2] and experiments. This work, thus, offers one of the first in-depth looks into the properties of oxygen vacancies in α -C.

2. Background

2.1. O Vacancy Configurations in SiO 2

Previous calculations demonstrated that a neutral O vacancy in α -Q and a-SiO 2 (in a-SiO 2 this defect is also called the Oxygen Deficiency Center (ODC) I [9,20]) is stabilized by forming a bond between two Si atoms that share two dangling electrons O 3 ≡Si-Si≡O 3 , or by partially, or fully, passivating these bonds with hydrogen atoms O 3 ≡Si-H H-Si≡O 3 [5]. The formation of a bare vacancy in α -Q was accompanied by a very strong displacement of the two Si atoms toward each other, see Figure 2a for the structure of this defect [3,5]. Some authors also suggested the existence of so-called unrelaxed neutral O vacancy in amorphous silica [21]. Hole trapping by neutral O vacancy, or hydrogen de-passivation reactions, can produce positively charged and paramagnetically active defect designated E centers, based on their EPR signals [22,23]. In α -Q, the most stable E center was the E 1 , which was a close relative to the family of E γ centers in a-SiO2, the general structure of which can be seen in Figure 2c. The E center naming convention uses Greek letters for a-SiO2 and Arabic numerals for α -Q. Recent electron paramagnetic resonance (EPR) studies distinguished at least 15 different EPR signals attributed to E centers in different types of α -Q perturbed by different defects and impurities [23,24,25].
Theoretical calculations predicted the existence of two lowest energy relaxed structures of E 1 center in pure α -Q: the so-called dimer, E d i , and puckered, E p u c k , configurations [1,10], see Figure 2 and a detailed description below. More E′ center configurations in a-SiO 2 have been proposed in e.g., references [8,26]. The dimer and puckered configurations differ strongly in terms of the values of isotropic hyperfine constants. The dimer type configurations are characterized by similar hyperfine constants on the two Si ions neighboring the vacancy of around 10 mT. Due to the similarity with observed hyperfine constants, in some papers this configuration has been attributed to the E δ center [3,5]. However, a thorough analysis of the experimental EPR spectra of the E δ center suggested a much more complex configuration for this defect involving four or five Si atoms [15,27]. In the dangling bond puckered configuration, the spin density is almost completely localized on one Si ion, which warrants a very strong hyperfine interaction with the Si of 42 mT, close to the measured experimental value [5,24].
In α -Q, the puckered configuration of the E center was shown to be more stable than the dimer configuration in several quite different calculations [5,28,29], supporting the results of EPR measurements [22]. This is related to the structure of α -Q, which has two non-equivalent positions of Si in the crystal structure and, correspondingly, two slightly different Si-O bonds. The so-called SiS-O and SiL-O bonds, where SiS is the Si atom associated with the shorter bond and SiL—with the slightly longer one. The calculations show that the SiL atom can back-bond to O to form the puckered configuration of the E 1 center, as discussed in more detail below. The most stable configuration of this E center has a forward projection of electrons on the Si S towards the vacant O site. The less stable E d i configuration has not yet been observed experimentally in either α -Q or a-SiO 2 [15].

2.2. Properties of α -Cristobalite

Tetragonal α -C has a structure similar to α -Q with rings of the same size in both, but α -C has a much lower density (see Figure 1). It is present in various planetary materials, e.g., lunar and martian rocks [30] and meteorites [31]. On the Earth, it occurs as a metastable phase in many geologic settings, for example as crystals deposited from vapor within the pores of volcanic rocks [32]. Most data on this silica phase have been obtained from geochemical analysis of volcanic rocks, where α -C crystals or paracrystals are found with a mixture of silica phases, such as opal-CT [33,34,35]. Extensive studies on the effect of pressure on α -C have contributed to our understanding of the phase [36,37,38,39,40,41,42]. Several earlier computational studies of α -C focused mainly on producing a complete data set of silica polymorphs and lacked detailed characterization of α -C itself [43,44,45,46].
Quite independently of its geochemical and mineral interest, α -C is popular in theoretical calculations of silica as a mimic of the bulk and surface properties of amorphous silica [47,48,49], as both have very similar densities and structural characteristics (e.g., see discussion in [50]). In particular, Emami et al. [51] showed that the α -C (10 1 ¯ ) surface had a similar silanol group density to a-SiO2. This makes α -C a good model for the a-SiO2 surface structure without the need for computationally demanding melt-quench methods to produce theoretical models of a-SiO2 surface structures. Therefore, interfaces of α -C with water [52,53] and with solids have been extensively calculated [48,54,55].
In α -C, all Si positions are equivalent, so the distinction of SiS and SiL cannot be made. However, as the E′ center is investigated in both, the Si atom involved in the back bonding must be labeled. To distinguish α -C as separate from the α -Q SiS notation, a new labeling was introduced, where Si A is the Si involved in the back bonding (SiL equivalent) and Si B is the other Si at the oxygen vacancy (SiS equivalent). For the sake of comparison, α -Q is also labeled in the same way, see Figure 3.

3. Materials and Methods

In this paper, we compared properties of O vacancies in two SiO 2 polymorphs, α -Q and α -C. Most previous calculations of α -C used molecular cluster models, or local DFT functionals and small periodic cells. In contrast, we used periodic DFT calculations with a hybrid functional and large periodic cells to more accurately predict the structure, electronic, EPR, and optical absorption properties of positively charged, neutral, and negatively charged O vacancies in α -C. These results were compared with those obtained in α -Q and amorphous silica.
Calculations were carried out using DFT implemented in the CP2K software package, using a Gaussian Plane Wave method [56]. A 324 atom 3 × 3 × 3 supercell of α -C and a 260 atom 3 × 3 × 3 supercell of α -Q were used. All calculations were carried out at the Γ -point of the Brillouin zone. Double- ζ Gaussian basis sets [57] were used with the GTH pseudopotential [58]. Plane wave cutoffs were set to 650 Ry (8844 eV) with a relative cutoff of 70 Ry (952 eV). The Broyden–-Fletcher–-Goldfarb–-Shanno algorithm (BFGS) was used to minimize the force on the atoms to 0.0001 a.u (0.005 eV) Å 1 [59].
The local Exchange–Correlation (XC) functional of Perdew–Burke–Ernzerhof (PBE) was used to optimize the cell parameters and atomic positions in all calculations [60]. The optimized structures from the PBE calculations were then used as input for the calculations with the non-local XC functional PBE0–TC–LRC [61]. These PBE0–TC–LRC calculations had an exchange interaction cutoff of 2 Å and a Hartree–Fock exchange of 25%. To reduce the computational cost of non-local functional calculations, the auxiliary density matrix method (ADMM) was used [62]. The same GTH pseudopotentials and plane wave cutoffs were used in both the PBE and PBE0–TC–LRC calculations.
The defect formation energies were calculated using equation [63]:
E f [ X q ] = E t o t [ X q ] E t o t [ b u l k ] + i n i μ i + q ( E V B M + μ e ) + E c o r r .
Here, E f [ X q ] is the defect formation energy, E t o t [ X q ] is the total energy of the defect system, E t o t [ b u l k ] is the total energy of the pristine system, i n i μ i is the chemical potential of the element i where n i atoms of that type are exchanged with a chemical reservoir, q is the defect charge, E V B M is the energy of the valence band maximum, μ e is the Fermi energy or the chemical potential of electrons, and E c o r r is a finite-size electrostatic correction using the anisotropic Makov-Payne correction [64,65]. Nudged Elastic Band (NEB) calculations were performed using the Climbing Image NEB (CI–NEB) method [66].
Optical absorption (OA) energies were calculated using the CP2K implementation of the Time-Dependent Density Functional Theory (TD–DFT), which does not include long range correction [67,68]. The truncation cutoff of the PBE0–TC functional was increased, making the function more like the PBE0 XC functional [69]. This was used with an enhanced TZVP–MOLOPT–GTH basis set, as the DZVP–MOLOPT–GTH basis sets did not accurately calculate the optical absorption peak energies. The local structure of the defect was re-optimized with these new parameters.
Positively charged vacancies are paramagnetic and are characterized using Electron Paramagnetic Resonance (EPR). To calculate the hyperfine EPR constants of defects, all-electron calculations were performed using CP2K and a Gaussian and Augmented Plane Wave basis [70]. For these all-electron calculations, pcJ basis sets were employed [71]. The EPR properties were then calculated within the framework of the density functional perturbation theory [72].

4. Results of Calculations

4.1. Pristine Crystals

First, we compared the properties of α -C and α -Q to understand the underlying structural differences that could contribute to any deviation in the stability of oxygen vacancies. As seen in Table 1, the calculated α -C and α -Q structures in this work were close to the experimental values. The space group of α -C was P4 1 2 1 2 with all lattice angles equaling 90. α -Q has a right-hand crystallization P3 2 21 space group with lattice angles α = β = 90 and γ = 120.
The band gap width of α -Q was studied in more detail than for α -C. The early experimental works by Bart et al. [73,74] showed a range in band gap values in α -Q of 8.6–9.5 eV, depending on the surface structure of α -Q. The bulk experimental band gap values were then determined to be 9.65 eV by Garvie et al. [75,76] A close value was also calculated with high-level theories, such as quasi-particle self-consistent GW [77].
Earlier calculations using the PBE0–TC–LRC XC functional (which was also used in this work) predicted the single-electron α -Q band gap at 8.6 eV [78]. Although this value corresponded well to that by Bart et al. [73,74], it was expected to be an underestimate of the bulk SiO 2 value. This underestimate was also seen in this work with a band gap of 8.5 eV.
The band gap values for α -C are only known from theoretical work and depend on the XC functional used. Higher-level XC functionals increase the accuracy of the band gap. The highest level theory calculation found in the literature used a meta–GGA XC functional to produce a single-electron band gap of 8.54 eV, [44]. Consequently, this value was used for comparison. Using the PBE0–TC–LRC XC functional, we calculated the band gap to be 8.57 eV, similar to the meta-GGA functional. Other lower levels of theory produced band gap values in the range of 5.5–10.3 eV [43,45,79,80,81].
The geometry of the two phases show that the Si–O–Si bond angles of α -C are larger than in α -Q. This indicates that α -C has a more open and flexible structure than α -Q, leading to a lower density when compared to α -Q. This is seen in Figure 1, where the large rings in α -C can be clearly seen. In addition, the Si–O bond lengths of α -C were slightly shorter than in α -Q, which might mean that any distortion introduced into α -C introduces more strain as the bonds are perturbed to a greater extent.
The calculated density of α -C (2.35 g/cm3) and a-SiO2 (2.20 g/cm3) [82] showed a difference of 0.15 g/cm3, and, thus, the band gaps were expected to differ with a lower density indicating a lower band gap value. However, El-Sayed et al. [83] calculated a band gap value of 8.9 eV for a-SiO2 using the hybrid HSE06 XC functional, while experimental band gap values of a-SiO2 ranged from 8.9–9.7 eV [75,84,85,86].
In summary, a comparison of the calculated α -Q band gap to the experimental and high-level theory data showed that the α -Q band gap was likely underestimated in this work. Therefore, it was assumed that the α -C band gap was also underestimated. The geometry of the two phases were close to previously observed and calculated results, with α -C having a less dense structure with larger open rings and shorter Si–O bonds.
Table 1. The structural properties of α -cristobalite and α -quartz compared to the previous experimental and theoretical literature. The geometric structure is also shown, with the band gap separated at the bottom of the table. All bond lengths and lattice vectors are in Å, all angles are in degrees (°), the density is in g/cm3, and the band gap is in eV.
Table 1. The structural properties of α -cristobalite and α -quartz compared to the previous experimental and theoretical literature. The geometric structure is also shown, with the band gap separated at the bottom of the table. All bond lengths and lattice vectors are in Å, all angles are in degrees (°), the density is in g/cm3, and the band gap is in eV.
Parameter α -Cristobalite α -Quartz
This WorkLiterature DataThis WorkLiterature Data
Lattice Vectors
a = b5.054.97 [87]4.934.91 [88]
c7.086.93 [87]5.435.40 [88]
Bond length
Si A -O1.604–1.6091.600–1.607 [87]1.6081.604 [88]
Si B -O1.604–1.6091.600–1.607 [87]1.6121.613 [88]
Bond Angles
O-Si-O108–111.8108.1–111.3 [89]109.2–110.5109.0–109.5 [88]
Si-O-Si150–153146.6 [87]144.8–145.1143.7 [88]
Density2.352.18–2.37 [44,45]2.412.47–2.70 [44,45,90]
Band gap8.578.54 [44]8.59.65 [75]
References [44,45,90] are the results of theoretical calculations. References [75,87,88,89] are experimental data.

4.2. Oxygen Vacancies

Three charge states of oxygen vacancies were investigated in both α -Q and α -C in the following different configurations: neutral Si-Si dimer ( V O 0 ), negative Si-Si dimer ( V O ), positive Si-Si dimer ( E d i ), puckered ( E p u c k ) and back-projected ( E b p ) configurations. The structures and highest occupied molecular orbitals (HOMOs) of these defects are shown in Figure 4 and Figure 5. The short-range geometric parameters of the four configurations are summarized in Table 2 in conjunction with Figure 4 and Figure 5. Below we discuss their geometric, electronic and optical properties in more detail.

4.2.1. Geometric Structure and Stability

Neutral Vacancies

Neutral O vacancies in SiO 2 are characterized by formation of a Si–Si bond sharing two electrons, which is accompanied by very significant displacement of two Si atoms in the direction of the vacant O site (see Figure 4a and Table 2). Hence, large changes are observed in the Si A -Si B distance with the introduction of oxygen vacancies in both α -C and α -Q. Compared to pristine structures, the formation of V O 0 was accompanied by a decrease in the distance between the two Si atoms by 0.81 Å and 0.69 Å for α -C and α -Q, respectively. The character of electron localization between the two Si atoms, forming the bond is seen in Figure 4a.

Positively Charged Vacancies

Some of the ways of creating positively charged vacancies in SiO 2 are via hole trapping or electronic excitation and ionization of neutral vacancies. Both are complex processes that can lead to significant reorganization of the surrounding lattice due to strong electron–phonon coupling. In static calculations presented here, positive charging of a periodic cell led to hole trapping on the Si–Si bond and the Si A -Si B distance increased, compared to the V O 0 defect, by 0.46 Å and 0.59 Å for α -C and α -Q, respectively (see Table 2). If no further reorganization occurs, this configuration of the positively charged defect is called a dimer, E d i . The Si–Si bond was significantly weakened, due to the loss of an electron and the consequent reduction in bonding order between the two silicon atoms (see Figure 4b), and, due to asymmetry of the lattice structure, more defect configurations are possible.
One such an asymmetric configuration in α -Q has been discovered in references [1,10]. where Si A in Figure 3b relaxed through the plane of its three back bonded oxygen ions and made a bond with the so-called back oxygen O B , as shown in Figure 5a. This Si atom carried a hole and Si B carried a strongly localized unpaired electron. This configuration was also predicted to be stable in some sites of a-SiO 2 (see discussion in ref. [26]) and is commonly called the puckered configuration. Therefore, we labeled it as E p u c k .
The E p u c k configuration had the largest Si A -Si B distance, due to the puckering of Si A through the plane of O atoms, as seen in Figure 5a, and this distance was close to the literature values [3,5,91]. This configuration was by 0.23 eV lower in energy than E d i . The CI–NEB calculated adiabatic barrier between the two configurations was 0.26 eV. Similar results were also found in [4,29]. In particular, the results presented in reference [29] demonstrated the dependence of the barrier height on the size of the crystal region allowed to relax. The barrier for the reverse process of conversion from E p u c k configuration back to the dimer was 0.38 eV, which was in good agreement with previous calculations [5].
The structure of α -C is more symmetric than that of α -Q with Si A and Si B being equivalent and O B further away from them in the corresponding rings (see Figure 3a). Therefore, the E p u c k configuration proved to be unstable in α -C. However, we found another stable configuration where the Si B , carrying the unpaired electron, was projected through the plane of its three back bonded oxygen ions, as shown in Figure 5b. A similar, so-called “back-projected” configuration, E b p , was found in DFT calculations of amorphous SiO 2 [8,92]. We note that in this configuration Si A carrying the hole relaxed into the plane of its nearest oxygen ions. The energy of the E b p configuration was 0.83 eV higher than that of the dimer. The NEB calculation to determine a barrier between E d i and E b p did not fully converge, but the barrier exceeded 1.0 eV. These results suggested that E d i was the dominant configuration of the positively charged O vacancy in α -C. Below, we calculated its optical and EPR characteristics, which could be used for identifying it experimentally.
The E b p configuration was also stable in α -Q (see Figure 5c), but in this case Si A formed a bond with the back oxygen O B . Therefore, the energy difference between the more stable E p u c k and E b p configurations was smaller than in α -C at 0.44 eV. Thus, there were at least three different configurations of positively charged vacancy in α -Q. The calculated EPR and optical properties of the E p u c k configuration were in agreement with the experimental data for the E 1 center in α -Q [2,10] (see also our results below) and this configuration was often associated with this defect. Below we compared the calculated optical and EPR parameters of this and other configurations in α -Q and α -C.

Negatively Charged Vacancies

The flexibility of the SiO 2 structure suggests that the lattice relaxation could promote trapping of an extra electron by a neutral vacancy, as suggested in earlier calculations [10] and, then, further investigated in [11,12]. The results presented here for α -Q and α -C were broadly in line with the previous calculations. The extra electron occupied the anti-bonding state, shown in Figure 4c. The Si A -Si B distance remained close to that in the neutral vacancy, but the bonds with the nearest O ions were elongated, due to repulsion with the extra electron, as shown in Table 2.

4.2.2. Charge Transition Levels

The charge transition levels (CTL) of the O vacancy in α -C and α -Q are compared in Figure 6. In α -C, the E d i center was the most stable defect at Fermi levels close to the valence band (VB). The transition from E d i to V O 0 (+1 to 0) occurred at 1.64 eV above the VBM, while that from the less stable E b p occurred at 0.81 eV. The transition from V O 0 to V O (0 to −1) occurred at 8.34 eV, close to the conduction band minimum (CBM) energy.
On the contrary, the α -Q CTL diagram shows that the E p u c k configuration was the most stable defect near the VB. The +1/0 transition from E p u c k to V O 0 occurred at 2.23 eV above the VBM, whereas the transition from E d i occurred at 1.97 eV. E b p was less stable than E d i , with a transition at 1.79 eV. The negatively charged V O was more stable in α -Q than in α -C, with a 0/−1 CT at 7.73 eV. The fact that the puckered configuration of the positively charged vacancy was more stable than the dimer configuration agreed with the results of previous calculations in α -Q, e.g., [1,5,10]. However, in α -C positive oxygen vacancies had not been previously investigated, so the contrast in the stability of E d i highlights the differences between α -C and α -Q structures.

4.2.3. Optical Absorption

To understand how these oxygen vacancy defects affect the observable properties of the two polymorphs, the optical transitions of α -C and α -Q were investigated using TD–DFT. Due to the CP2K implementation of TD–DFT, the results presented in this section were calculated using the PBE0–TC XC functional, in contrast with the other sections where the long-range GGA correction (LRC) was used. PBE0–TC uses a constant fraction of exact exchange in the short range until a specified cut-off distance, at which it sets exact exchange to zero. To check whether the TD–DFT spectra were sensitive to the cut-off used in truncating the exact exchange, apart from using the 2 Å cut-off, we also used the maximum cut-off for each silica phase, according to the cell size, 6.4 Å for α -Q and 7.5 Å for α -C, respectively. This ensured that the lack of LRC did not introduce artifacts and also meant that the XC functional was almost identical to the classic PBE0 XC functional.
The results shown in Table 3 demonstrate the effect of the cut-off radius on transition energies. We first considered the results for the E p u c k configuration in α -Q. In Table 3 we only show transitions with oscillator strengths exceeding 0.1. The transitions at 6.18 eV for the 6.4 Å cut-off (and at 6.43 eV for the 2.0 Å cut-off) had a strong contribution from the localized orbitals on the Si B and atoms in the surrounding rings. The unpaired electron mixed with the Si p d and Op states of the surrounding rings, and this formed the bottom of the conduction band. The transition energies were close to the previously calculated values of 5.8–6.1 eV [2,93].
The optical absorption in the 5.8–6.0 eV region in SiO 2 was usually associated with the E center and was attributed to its E p u c k configuration. Experimental optical absorption spectra in α -Q are quite scarce. Guzzi et al. [94] attributed the 6.08 eV peak measured at 77K to the E 1 center in neutron-irradiated α -Q. Our results at the 6.4 Å cut-off supported this attribution. The nature of electronic excitations responsible for optical absorption of the E p u c k center was discussed in reference [95] where it was proposed that it was caused by electron transfer from Si B to Si A carrying the hole (see Figure 5a). This model was confirmed by DFT calculations in cluster models [2,93]. However, we observed significant contribution of the charge transfer transition only at the small cut-off. The difference in the nature of excitation with the results of [2,93] could be attributed to the much larger size of the periodic cell and more delocalized nature of the CBM states in our calculations, compared to the relatively small clusters used by [2,93]. We observed that optical transition energies were quite robust with respect to the exchange cut-off value and method of calculation.
Similar trends, with respect to the cut-off radius, can be seen also for the V O 0 center in α -C, where the transition energies with the 6.4 Å radius were reduced by about 0.25 eV with respect to those at 2.0 Å. The main transitions at 7.62 and 7.72 eV were in very good agreement with the absorption peak at 7.6 eV [9,94] usually attributed to the ODC(I) defect, where the proposed structure was a relaxed diamagnetic oxygen vacancy with two electrons localized between two silica [96,97,98] in quartz and SiO 2 glass samples. The main transition was between the bonding and anti-bonding σ states localized at the Si–Si bond. The strongly localized nature of the bond might explain why the predicted transitions for α -C are so similar to those in α -Q and a-SiO 2 , as the long-range structure of the silica might not affect this transition. Two further peaks were seen at 8.2 and 8.48 eV, which could be related to those observed experimentally [94] and were due to localized high-energy π defect states in the conduction band.
The calculated transition energies for the E d i configuration in α -C at 7.56 eV and 7.72 eV (see Table 3) overlapped strongly with those for V O 0 but had lower oscillator strengths. Therefore other data, such as luminescence and EPR spectra, are required for their identification. We noted the similarity with F and F + centers in MgO and CaO in this respect, where the absorption spectra of these defects also strongly overlapped (see e.g., ref. [99]). Below we present the results of our calculations for the EPR parameters of E d i and E p u c k configurations.

4.2.4. EPR Parameters

The pseudopotentials used so far in this paper were unsuitable for the calculation of EPR parameters, and instead computationally expensive all-electron (AE) calculations had to be performed with an expanded basis set suitable for EPR calculations. Benchmarking calculations demonstrated that, as the quality of the AE basis set was systematically increased, the AE optimized geometry tended to the pseudopotential geometries already obtained. Therefore, the previously calculated defect geometries were maintained, and, then, an AE calculation of the wave function was performed for these geometries. The pcJ basis sets were employed [71], which systematically increased in quality from pcJ-0 to pcJ-2.
The results of these calculations are summarized in Table 4. The isotropic hyperfine coupling constants were calculated for the considered defects, which were dominated by Fermi contact interactions and neglected the subtler anisotropic splittings. As can be seen in Table 4, there were significant differences in the calculated hyperfine couplings between the two configurations E d i (in α -C and α -Q) and E p u c k (in α -Q), due to the differing degree of electronic localization in the two cases. In the E p u c k configuration, the unpaired electron associated with the defect state was localized on a single Si, and, hence, a strong hyperfine coupling was observed. In contrast, in the E d i configuration in both phases, the defect state was delocalized across two Si atoms, and, hence. a much weaker coupling was observed. Although there were important structural differences between α -C and α -Q, our results showed that they did not significantly affect the resultant EPR parameters. The very minor difference in the calculated values for the two Si atoms in the E d i defect occurred because the introduction of the oxygen vacancy lowered the symmetry and the two Si sites were no longer symmetrically equivalent.
We note that the hyperfine coupling constant obtained here for the E p u c k configuration in α -Q is in good agreement with 40.15–42.5 mT attributed to the E 1 center [24,100]. These values for the E d i configuration are similar to the 10 mT seen in the literature for E δ center [101], as well as in calculations of dimer vacancy configurations [11]. However, the analysis of experimental EPR spectra of the E δ center in silica glass concluded [15,27] that this center had a more complex structure than the E d i configuration described here and in previous calculations. Since E d i was the only stable configuration of positively charged vacancy in α -C, EPR measurements for the positively charged O vacancy in α -C offered a unique opportunity to identify this configuration and to verify the theoretical predictions.

5. Discussion and Conclusions

Predicting properties of defects in solids is a growing activity involving advanced DFT calculations and, increasingly, Machine Learning. Still, the properties of even basic intrinsic defects observed in binary oxides, such as SiO 2 , are not fully understood. Creating and proving models of defects in solids is still a major challenge. In this work we used the O vacancy in α -quartz to calibrate our computational procedure and test previous predictions using less accurate methods. We then compared the structures and properties of O vacancies in two SiO 2 crystalline phases, α -quartz and α -cristobalite, to investigate how structural difference affects these properties and to find ways to identify the dimer configuration of the positively charged O vacancy. The optical absorption spectra of O vacancies in α -C and α -Q were calculated using TD–DFT and isotropic hyperfine constants were calculated for positively charged configurations to aid in experimental identification.
The results for neutral and negatively charged vacancies are similar in both phases and agree with previous calculations. However, there is a significant difference in the number and relative stability of configurations of positively charged vacancies. We found three configurations in α -Q, with E p u c k being the most stable, followed by the E d i and E b p configurations. There were only two stable configurations in α -C, with E d i being much more stable than E b p . This difference was attributed to the difference in crystal structures and short-range environment of Si atoms in both phases. The back-projected, E b p , configuration of the positively charged vacancy was predicted in previous calculations of O vacancies in amorphous SiO 2 , but had not been considered in crystalline phases. To the best of our knowledge, it has not been attributed, so far, to any experimental signal.
The optical and EPR characteristics calculated for the E p u c k configuration in α -Q were close to the experimental data for the best-studied E 1 center. Periodic cells used in our calculations were much larger than 72-atom cells typically used in the past, as well as molecular clusters. They better accounted for the long-range defect-induced lattice distortions predicted by embedded cluster models [3]. However, the optical transition energies were overestimated with respect to the maxima of absorption spectra observed in α -Q (6.08 eV [94]) and a-SiO 2 (5.8 eV [9]). TD–DFT showed contributions of both charge-transfer excitations and those into the states more localized on one Si and the surrounding atoms into this absorption peak. We noted that, although our results, in terms of the structure, relative energies, optical and EPR properties were in line with many previous studies of the E p u c k configuration in α -Q, the accuracy of these calculations, as well as experimental optical measurements, should still be improved to make a fully confident attribution.
On the other hand, both the E d i and E b p configurations have never been experimentally identified in either α -Q or a-SiO 2 . By comparison of hyperfine constants, earlier theoretical simulations associated the E d i configuration with the E δ center observed in a-SiO 2 [5]. However, several experimental EPR studies ruled out this attribution [15,27]. The predicted stability of the dimer configuration in α -C might provide an opportunity to finally observe and investigate it using the optical transition energies and hyperfine splittings calculated here.
We note that the fact that E d i and E b p configurations had not yet been identified experimentally in a-SiO 2 , in spite of numerous theoretical predictions, may indicate that the O vacancy did not serve as the main precursor for formation of E center. Experiments using field-dependent recombination of holes trapped in thermally grown a-SiO 2 films with injected electrons [102] revealed that the paramagnetic state of the E center did not always correlate with the entity bearing the positive charge. It has been suggested that the positive charge is protonic in origin. Consequently, the O 3 ≡Si–H entity in a-SiO 2 is suggested to be a possible E center precursor [102], whereupon hole trapping hydrogen dissociates in the form of a proton leaving behind a neutral paramagnetic E center, a hypothesis later corroborated by theoretical calculations [103].
To summarize, these results may help defect characterization in silica phases for geochemistry and space exploration. They demonstrate how subtle differences in the local structure of silica polymorphs strongly affect the structure and stability of basic intrinsic defects. Although the surface of α -C is a good mimic of that for a-SiO 2 , positively charged O vacancies in the bulk show significant differences in the structure and relative stability of configurations. This may affect modeling of charge transfer processes at interfaces simulated using α -C.

Author Contributions

Conceptualization, A.L.S. and K.L.M.; methodology, K.L.M. and T.R.D.; validation, K.L.M. and T.R.D.; formal analysis, K.L.M. and A.L.S.; investigation, K.L.M., T.R.D. and T.C.F.; data curation, K.L.M.; writing—original draft preparation, K.L.M. and A.L.S.; writing—review and editing, T.R.D. and A.L.S.; visualization, K.L.M.; supervision, A.L.S.; project administration, A.L.S.; funding acquisition, A.L.S. All authors have read and agreed to the published version of the manuscript.

Funding

K.L.M. was supported by the the Engineering and Physical Sciences Research council (EPSRC) grant EP/T517793/1. T.D. is grateful to EPSRC (grant number EP/R034540/1) for financial support. Via our membership of the UK’s HEC Materials Chemistry Consortium, which is funded by EPSRC (EP/R029431), this work used the ARCHER2 UK National Supercomputing Service (http://www.archer2.ac.uk, accessed on 2 Feburary 2023).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data required to reproduce these findings are available from the corresponding author upon reasonable request.

Acknowledgments

We would like to thank Jack Strand for his assistance with Python scripts, as well as Matthew Watkins and Jack Strand for their support in CP2K applications. We would also like to acknowledge Anna Regoutz and her group for valuable discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DFTDensity Functional Theory
α -C α -cristobalite
α -Q α -quartz
a-SiO 2 Amorphous silicon dioxide
CTCharge transition
ODCOxygen deficient center
EPRElectron paramagnetic resonance
Opal-CTOpal cristobalite tridymite
LDALocal density approximation
GTHGoedecker-Teter-Hutter
BFGSBroyden–Fletcher–Goldfarb–Shanno
XCExchange-correlation
PBEPerdew-Burke-Ernzerhof
ADMMAuxiliary density matrix method
TD-DFTTime-dependent density functional theory
TC-LRCTruncated coulomb long-range correction
NEBNudged elastic band
CI-NEBClimbing image nudged elastic band
GGAGeneralized gradient approximation
HSEHeyd-Scuseria-Ernzerhof
VBValence band
CBConduction band
CBMConduction band minimum
CTLCharge transition level
TSTransition state
OAOptical absorption
AEall-electron

References

  1. Feigl, F.J.; Fowler, W.B.; Yip, K.L. Oxygen Vacancy Model for the E′-center in SiO2. Solid State Commun. 1974, 14, 225–229. [Google Scholar] [CrossRef]
  2. Pacchioni, G.; Ierano, G.; Márquez, A.M. Optical absorption and nonradiative decay mechanism of E′ center in silica. Phys. Rev. Lett. 1998, 81, 377–380. [Google Scholar] [CrossRef]
  3. Sushko, P.V.; Mukhopadhyay, S.; Mysovsky, A.S.; Sulimov, V.B.; Taga, A.; Shluger, A. Structure and properties of defects in amorphous silica: New insights from embedded cluster calculations. J. Phys. Condens. Matter 2005, 17, S2115. [Google Scholar] [CrossRef]
  4. Boero, M.; Pasquarello, A.; Sarnthein, J.; Car, R. Structure and Hyperfine Parameters of E1′ Centers in α-Quartz and in Vitreous SiO2. Phys. Rev. Lett. 1997, 78, 887–890. [Google Scholar] [CrossRef]
  5. Blöchl, P.E. First-principles calculations of defects in oxygen-deficient silica exposed to hydrogen. Phys. Rev. B 2000, 62, 6158–6179. [Google Scholar] [CrossRef]
  6. Lu, Z.Y.; Nicklaw, C.J.; Fleetwood, D.; Schrimpf, R.D.; Pantelides, S. Structure, Properties, and Dynamics of Oxygen Vacancies in Amorphous SiO2. Phys. Rev. Lett. 2002, 89, 285505. [Google Scholar] [CrossRef]
  7. Chadi, D.J. Negative-U property of the oxygen vacancy defect in SiO2 and its implication for the E′ center in α-quartz. Appl. Phys. Lett. 2003, 83, 437–439. [Google Scholar] [CrossRef]
  8. Giacomazzi, L.; Martin-Samos, L.; Boukenter, A.; Ouerdane, Y.; Girard, S.; Richard, N. EPR parameters of E′ centers in v-SiO2 from first-principles calculations. Phys. Rev. B 2014, 90, 014108. [Google Scholar] [CrossRef]
  9. Skuja, L. Optically active oxygen-deficiency-related centers in amorphous silicon dioxide. J.-Non-Cryst. Solid. 1998, 239, 16–48. [Google Scholar] [CrossRef]
  10. Rudra, J.K.; Fowler, W.B. Oxygen vacancy and the E′1 center in crystalline SiO2. Phys. Rev. B 1987, 35, 8223–8230. [Google Scholar] [CrossRef]
  11. Sushko, P.V.; Mukhopadhyay, S.; Stoneham, A.M.; Shluger, A.L. Oxygen vacancies in amorphous silica: Structure and distribution of properties. Microelectron. Eng. 2005, 80, 292–295. [Google Scholar] [CrossRef]
  12. Kimmel, A.; Sushko, P.; Shluger, A.; Bersuker, G. Positive and Negative Oxygen Vacancies in Amorphous Silica. ECS Trans. 2009, 19, 3–17. [Google Scholar] [CrossRef]
  13. Munde, M.S.; Gao, D.Z.; Shluger, A.L. Diffusion and aggregation of oxygen vacancies in amorphous silica. J. Phys. Condens. Matter 2017, 29, 245701. [Google Scholar] [CrossRef] [PubMed]
  14. Ganduglia-Pirovano, M.V.; Hofmann, A.; Sauer, J. Oxygen vacancies in transition metal and rare earth oxides: Current state of understanding and remaining challenges. Surf. Sci. Rep. 2007, 62, 219–270. [Google Scholar] [CrossRef]
  15. Jivanescu, M.; Stesmans, A.; Afanas’ev, V.V. Multifrequency ESR analysis of the Eδ′ defect in a-SiO2. Phys. Rev. B 2011, 83, 094118. [Google Scholar] [CrossRef]
  16. Arrigoni, M.; Madsen, G.K.H. Evolutionary computing and machine learning for discovering of low-energy defect configurations. Npj Comput. Mater. 2021, 7, 71. [Google Scholar] [CrossRef]
  17. Milardovich, D.; Jech, M.; Waldhoer, D.; El-Sayed, A.M.B.; Grasser, T. Machine Learning Prediction of Defect Structures in Amorphous Silicon Dioxide. In Proceedings of the IEEE 51st European Solid-State Device Research Conference (ESSDER2021), Grenoble, France, 13–22 September 2021; pp. 239–242. [Google Scholar] [CrossRef]
  18. Qiao, Z.; Liu, Q.; Zhang, S.; Wu, Y. The mineralogical characteristics between opaline silica in bentonite and α-cristobalite. Solid State Sci. 2019, 96, 105948. [Google Scholar] [CrossRef]
  19. Önal, M.; Kahraman, S.; Sarikaya, Y. Differentiation of α-cristobalite from opals in bentonites from Turkey. Appl. Clay Sci. 2007, 35, 25–30. [Google Scholar] [CrossRef]
  20. Griscom, D.L. Trapped-electron centers in pure and doped glassy silica: A review and synthesis. J.-Non-Cryst. Solid. 2011, 357, 1945–1962. [Google Scholar] [CrossRef]
  21. Imai, H.; Arai, K.; Hosono, H.; Abe, Y. Two types of oxygen-deficient centers in synthetic silica glass. Phys. Rev. B 1988, 38, 12773–12775. [Google Scholar] [CrossRef]
  22. Weeks, R.A. Paramagnetic resonance of lattice defects in irradiated quartz. J. Appl. Phys. 1956, 27, 1376–1381. [Google Scholar] [CrossRef]
  23. Weeks, R.A.; Magruder, R.H.; Stesmans, A. Review of some experiments in the 50 year saga of the E′ center and suggestions for future research. J.-Non-Cryst. Solid. 2008, 354, 208–216. [Google Scholar] [CrossRef]
  24. Mashkovtsev, R.I.; Pan, Y. Nature of paramagnetic defects in alpha-quartz: Progresses in the first decade of the 21st century. In New Developments in Quartz Research: Varieties, Crystal Chemistry and Uses in Technology; Nova Science Publishers: New York, NY, USA, 2013; pp. 66–104. [Google Scholar]
  25. Mashkovtsev, R.I.; Nepomnyashchikh, A.I.; Zhaboedov, A.P.; Paklin, A.S. EPR study of the E′ defects in optical glasses and cristobalite. Europhys. Lett. 2021, 133, 14003. [Google Scholar] [CrossRef]
  26. Mukhopadhyay, S.; Sushko, P.V.; Stoneham, A.M.; Shluger, A.L. Modeling of the structure and properties of oxygen vacancies in amorphous silica. Phys. Rev. B 2004, 70, 195203. [Google Scholar] [CrossRef]
  27. Buscarino, G.; Agnello, S.; Gelardi, F.; Parlato, A. Electron paramagnetic resonance investigation on the hyperfine structure of the Eδ′ center in amorphous silicon dioxide. J.-Non-Cryst. Solid. 2007, 353, 518–521. [Google Scholar] [CrossRef]
  28. Fowler, W.B.; Edwards, A.H. Theory of defects and defect processes in silicon dioxide. J.-Non-Cryst. Solid. 1997, 222, 33–41. [Google Scholar] [CrossRef]
  29. Mysovsky, A.S.; Sushko, P.V.; Mukhopadhyay, S.; Edwards, A.H.; Shluger, A.L. Calibration of embedded-cluster method for defect studies in amorphous silica. Phys. Rev. B 2004, 69, 085202. [Google Scholar] [CrossRef]
  30. Kayama, M.; Nagaoka, H.; Niihara, T. Lunar and martian silica. Minerals 2018, 8, 8070267. [Google Scholar] [CrossRef]
  31. Brearley, A.J.; Jones, R.H. Chondritic meteorites. In Planetary Materials; Reviews in Mineralogy; De Gruyter: Berlin, Germany, 1998; Volume 36, pp. 3.1–3.398. [Google Scholar]
  32. Schipper, C.I.; Rickard, W.D.; Reddy, S.M.; Saxey, D.W.; Castro, J.M.; Fougerouse, D.; Quadir, Z.; Conway, C.; Prior, D.J.; Lilly, K. Volcanic SiO2-cristobalite: A natural product of chemical vapor deposition. Am. Mineral. 2020, 105, 510–524. [Google Scholar] [CrossRef]
  33. Smith, D.K. Opal, cristobalite, and tridymite: Noncrystallinity versus crystallinity, nomenclature of the silica minerals and bibliography. Power Diffr. 1998, 13, 2–19. [Google Scholar] [CrossRef]
  34. Correcher, V.; Garcia-Guinea, J.; Bustillo, M.A.; Garcia, R. Study of the thermoluminescence emission of a natural α-cristobalite. Radiat. Eff. Defects Solid. 2009, 164, 59–67. [Google Scholar] [CrossRef]
  35. Dağlar, S.; Kahya, N.D.; Ustunisik, G.; Onal, M. Thermal Crystallization Kinetics of an Opal-like Biogenic Silica. Silicon 2022, 14, 7211–7217. [Google Scholar] [CrossRef]
  36. Dera, P.; Lazarz, J.D.; Prakapenka, V.B.; Barkley, M.; Downs, R.T. New insights into the high-pressure polymorphism of SiO2 cristobalite. Phys. Chem. Miner. 2011, 38, 517–529. [Google Scholar] [CrossRef]
  37. Salvadó, M.A.; Pertierra, P.; Morales-García, A.; Menéndez, J.M.; Recio, J.M. Understanding chemical changes across the α-cristobalite to stishovite transition path in silica. J. Phys. Chem. C 2013, 117, 8950–8958. [Google Scholar] [CrossRef]
  38. Parise, J.B.; Yeganeh-Haeri, A.; Weidner, D.J.; Jorgensen, J.D.; Saltzberg, M.A. Pressure-induced phase transition and pressure dependence of crystal structure in low (α) and Ca/Al-doped cristobalite. J. Appl. Phys. 1994, 75, 1361–1367. [Google Scholar] [CrossRef]
  39. Hatch, D.; Ghose, S. The α-β phase transition in cristobalite, SiO2. Phys. Chem. Miner. 1991, 17, 554–562. [Google Scholar] [CrossRef]
  40. Garg, N.; Sharma, S.M. Classical molecular dynamical simulations of high pressure behavior of alpha cristobalite (SiO2). J. Phys. Condens. Matter 2007, 19, 456201. [Google Scholar] [CrossRef]
  41. Cernok, A.; Marquardt, K.; Caracas, R.; Bykova, E.; Habler, G.; Liermann, H.P.; Hanfland, M.; Mezouar, M.; Bobocioiu, E.; Dubrovinsky, L. Compressional pathways of α-cristobalite, structure of cristobalite X-I, and towards the understanding of seifertite formation. Nat. Commun. 2017, 8, 15647. [Google Scholar] [CrossRef]
  42. Shelton, H.; Bi, T.; Zurek, E.; Smith, J.; Dera, P. The Ideal Crystal Structure of Cristobalite X-I: A Bridge in SiO2 Densification. J. Phys. Chem. C 2018, 122, 17437–17446. [Google Scholar] [CrossRef]
  43. Güler, E.; Uğur, G.; Uğur, Ş.; Güler, M. A theoretical study for the band gap energies of the most common silica polymorphs. Chin. J. Phys. 2020, 65, 472–480. [Google Scholar] [CrossRef]
  44. Warmbier, R.; Mohammed, F.; Quandt, A. Optical and other material properties of SiO2 from ab initio studies. Opt. Eng. 2014, 53, 071808. [Google Scholar] [CrossRef]
  45. Sevik, C.; Bulutay, C. Theoretical study of the insulating oxides and nitrides: SiO2, GeO2, Al2O3, Si3N4, and Ge3N4. J. Mater. Sci. 2007, 42, 6555–6565. [Google Scholar] [CrossRef]
  46. Xu, Y.N.; Ching, W.Y. Electronic and optical properties of all polymorphic forms of silicon dioxide. Phys. Rev. B 1991, 44, 11048–11059. [Google Scholar] [CrossRef] [PubMed]
  47. Akiyama, T.; Kawamoto, K.; Kageshima, H.; Uematsu, M.; Nakamura, K.; Ito, T. A first-principles study of O2 incorporation and its diffusion in compressively strained high-density silicon oxides. Thin Solid Films 2006, 508, 311–314. [Google Scholar] [CrossRef]
  48. Nagao, K.; Neaton, B.; Ashcroft, W. First-principles study of adhesion at Cu/SiO2 interfaces. Phys. Rev. B 2003, 68, 125403. [Google Scholar] [CrossRef]
  49. Rimola, A.; Costa, D.; Sodupe, M.; Lambert, J.F.; Ugliengo, P. Silica surface features and their role in the adsorption of biomolecules: Computational modeling and experiments. Chem. Rev. 2013, 113, 4216–4313. [Google Scholar] [CrossRef]
  50. Huang, L.; Kieffer, J. Amorphous-amorphous transitions in silica glass. I. Reversible transitions and thermomechanical anomalies. Phys. Rev. B 2004, 69, 224203. [Google Scholar] [CrossRef]
  51. Emami, F.S.; Puddu, V.; Berry, R.J.; Varshney, V.; Patwardhan, S.V.; Perry, C.C.; Heinz, H. Erratum: Force Field and a Surface Model Database for Silica to Simulate Interfacial Properties in Atomic Resolution. Chem. Mater. 2016, 26, 2647–2658. [Google Scholar] [CrossRef]
  52. Yang, J.; Meng, S.; Xu, L.; Wang, E.G. Water adsorption on hydroxylated silica surfaces studied using the density functional theory. Phys. Rev. B 2005, 71, 035413. [Google Scholar] [CrossRef] [Green Version]
  53. Valencia, E. Calculation of the adsorption potential energy of water vapour on α-cristobalite. J. Chem. Soc. Faraday Trans. 1994, 90, 2555–2559. [Google Scholar] [CrossRef]
  54. Shan, T.R.; Devine, B.D.; Phillpot, S.R.; Sinnott, S.B. Molecular dynamics study of the adhesion of Cu/SiO2 interfaces using a variable-charge interatomic potential. Phys. Rev. B 2011, 83, 115327. [Google Scholar] [CrossRef]
  55. Tang, C.; Cai, L. The effect of surface heterogeneity between α-quartz and α-cristobalite on adsorption behaviors toward Cu2+ solution. Colloids Surfaces A Physicochem. Eng. Asp. 2021, 609, 125651. [Google Scholar] [CrossRef]
  56. Vandevondele, J.; Krack, M.; Mohamed, F.; Parrinello, M.; Chassaing, T.; Hutter, J. Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach. Comput. Phys. Commun. 2005, 167, 103–128. [Google Scholar] [CrossRef]
  57. VandeVondele, J.; Hutter, J. Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases. J. Chem. Phys. 2007, 127, 114105. [Google Scholar] [CrossRef] [PubMed]
  58. Goedecker, S.; Teter, M. Separable dual-space Gaussian pseudopotentials. Phys. Rev. B 1996, 54, 1703–1710. [Google Scholar] [CrossRef] [PubMed]
  59. Pfrommer, B.G.; Côté, M.; Louie, S.G.; Cohen, M.L. Relaxation of crystals with the quasi-Newton method. J. Comput. Phys. 1997, 131, 233–240. [Google Scholar] [CrossRef]
  60. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef]
  61. Guidon, M.; Hutter, J.; VandeVondele, J. Robust periodic Hartree-Fock exchange for large-scale simulations using Gaussian basis sets. J. Chem. Theory Comput. 2009, 5, 3010–3021. [Google Scholar] [CrossRef]
  62. Guidon, M.; Hutter, J.; Vandevondele, J. Auxiliary density matrix methods for Hartree-Fock exchange calculations. J. Chem. Theory Comput. 2010, 6, 2348–2364. [Google Scholar] [CrossRef] [Green Version]
  63. Zhang, S.B.; Northrup, J.E. Chemical Potential Dependence of Defect Formation Energies in GaAs: Application to Ga Self-Diffusion. Phys. Rev. Lett. 1991, 67, 2339–2342. [Google Scholar] [CrossRef]
  64. Makov, G.; Payne, M.C. Periodic boundary conditions in ab initio calculations. Phys. Rev. B 1995, 51, 4014–4022. [Google Scholar] [CrossRef] [PubMed]
  65. Murphy, S.T.; Hine, N.D. Anisotropic charge screening and supercell size convergence of defect formation energies. Phys. Rev. B 2013, 87, 094111. [Google Scholar] [CrossRef]
  66. Henkelman, G.; Uberuaga, B.P.; Jónsson, H. Climbing image nudged elastic band method for finding saddle points and minimum energy paths. J. Chem. Phys. 2000, 113, 9901–9904. [Google Scholar] [CrossRef]
  67. Iannuzzi, M.; Chassaing, T.; Wallman, T.; Hutter, J.; Hutter, J. Ground and Excited State Density Functional Calculations with the Gaussian and Augmented-Plane-Wave Method. Chimia 2005, 59, 499–503. [Google Scholar] [CrossRef]
  68. Strand, J.; Chulkov, S.K.; Watkins, M.B.; Shluger, A.L. First principles calculations of optical properties for oxygen vacancies in binary metal oxides. J. Chem. Phys. 2019, 150, 044702. [Google Scholar] [CrossRef]
  69. Adamo, C.; Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys. 1999, 110, 6158–6170. [Google Scholar] [CrossRef]
  70. Lippert, G.; Hutter, J.; Parrinello, M. The Gaussian and augmented-plane-wave density functional method for ab initio molecular dynamics simulations. Theor. Chem. Acc. 1999, 103, 124–140. [Google Scholar] [CrossRef]
  71. Jensen, F. The optimum contraction of basis sets for calculating spin–spin coupling constants. Theor. Chem. Acc. 2010, 126, 371–382. [Google Scholar] [CrossRef]
  72. Weber, V.; Iannuzzi, M.; Giani, S.; Hutter, J.; Declerck, R.; Waroquier, M. Magnetic linear response properties calculations with the Gaussian and augmented-plane-wave method. J. Chem. Phys. 2009, 131, 014106. [Google Scholar] [CrossRef] [Green Version]
  73. Bart, F.; Gautier, J.; Henriot, M. (011¯0) α-quartz Surface: A LEED, XANES and ELS Study. Surf. Sci. 1992, 274, 317–328. [Google Scholar] [CrossRef]
  74. Bart, F.; Gautier, M.; Jollet, F.; Duraud, J. Electronic structure of the (0001) and (101¯0) quartz surfaces and of their defects as observed by reflection electron energy loss spectroscopy (REELS). Surf. Sci. 1994, 306, 342–358. [Google Scholar] [CrossRef]
  75. Garvie, L.A.J.; Rezb, P.; Alvarezb, J.R.; Buseck’, P.R. Interband Transitions of Crystalline and Amorphous SiO2: An Electron Energy-Loss Spectroscopy (EELS) Study of the Low-Loss Region. Solid State Commun. 1998, 106, 303–307. [Google Scholar] [CrossRef]
  76. Garvie, L.A.J.; Rez, P.; Alvarez, J.R.; Buseck, P.R.; Craven, A.J.; Brydson, R. Bonding in α-quartz (SiO2): A view of the unoccupied states. Am. Mineral. 2000, 85, 732–738. [Google Scholar] [CrossRef]
  77. Kresse, G.; Marsman, M.; Hintzsche, L.E.; Flage-Larsen, E. Optical and electronic properties of Si3N4 and α-SiO2. Phys. Rev. B 2012, 85, 045205. [Google Scholar] [CrossRef]
  78. El-Sayed, A.M.; Watkins, M.B.; Afanas’ev, V.V.; Shluger, A.L. Nature of intrinsic and extrinsic electron trapping in SiO2. Phys. Rev. B 2014, 89, 125201. [Google Scholar] [CrossRef]
  79. Maj, S. Energy gap and density in SiO2 polymorphs. Phys. Chem. Miner. 1988, 15, 271–273. [Google Scholar] [CrossRef]
  80. Li, Y.P.; Ching, W.Y. Band structures of all polycrystalline forms of silicon dioxide. Phys. Rev. B 1985, 31, 2172–2179. [Google Scholar] [CrossRef]
  81. Gnani, E.; Reggiani, S.; Colle, R.; Rudan, M. Band-Structure Calculations Of SiO2 by Means of Hartree-Fock and density-functional techniques. IEEE Trans. Electron Devices 2000, 47, 1795–1803. [Google Scholar] [CrossRef]
  82. Mozzi, R.L.; Warren, B.E. Structure of Vitreous Silica. J. Appl. Crystallogr. 1969, 2, 164–172. [Google Scholar] [CrossRef]
  83. El-Sayed, A.M.; Watkins, M.B.; Shluger, A.L.; Afanas’ev, V.V. Identification of intrinsic electron trapping sites in bulk amorphous silica from ab initio calculations. Microelectron. Eng. 2013, 109, 68–71. [Google Scholar] [CrossRef]
  84. DiStefano, T.; Eastman, D. The Band Edge of Amorphous SiO2 by Photoinjection and Photoconductivity Measurements. Solid State Commun. 1971, 9, 2259–2261. [Google Scholar] [CrossRef]
  85. Nithianandam, V.J.; Schnatterly, S.E. Soft-x-ray emission and inelastic electron-scattering study of the electronic excitations in amorphous and crystalline silicon dioxide. Phys. Rev. B 1998, 38, 5547–5553. [Google Scholar] [CrossRef]
  86. Weinberg, Z.A.; Rubloff, G.W.; Bassous, E. Transmission, photoconductivity, and the experimental band gap of thermally grown SiO2 films. Phys. Rev. B 1979, 19, 3107–3117. [Google Scholar] [CrossRef]
  87. Lee, S.; Xu, H. Using powder XRD and pair distribution function to determine anisotropic atomic displacement parameters of orthorhombic tridymite and tetragonal cristobalite. Acta Crystallogr. Sect. B Struct. Sci. Cryst. Eng. Mater. 2019, 75, 160–167. [Google Scholar] [CrossRef]
  88. Baur, W.H. In search of the crystal structure of low quartz. Z. Krist. 2009, 224, 580–592. [Google Scholar] [CrossRef]
  89. Pluth, J.J.; Smith, J.V.; Faber, J. Crystal structure of low cristobalite at 10, 293, and 473 K: Variation of framework geometry with temperature. J. Appl. Phys. 1985, 57, 1045–1049. [Google Scholar] [CrossRef]
  90. Anderson, O.L.; Schreiber, E. The Relation between Refractive Index and Density of Minerals Related to the Earth’s Mantle. J. Geophys. Res. 1965, 70, 1463–1471. [Google Scholar] [CrossRef]
  91. Capron, N.; Carniato, S.; Lagraa, A.; Boureau, G.; Pasturel, A. Local density approximation and generalized gradient approximation calculations for oxygen and silicon vacancies in silica. J. Chem. Phys. 2000, 112, 9543–9548. [Google Scholar] [CrossRef]
  92. Mukhopadhyay, S.; Sushko, P.V.; Edwards, A.H.; Shluger, A.L. Calculation of relative concentrations of E′ centres in amorphous silica. J.-Non-Cryst. Solid. 2004, 345-346, 703–709. [Google Scholar] [CrossRef]
  93. Raghavachari, K.; Ricci, D.; Pacchioni, G. Optical properties of point defects in SiO2 from time-dependent density functional theory. J. Chem. Phys. 2002, 116, 825–831. [Google Scholar] [CrossRef]
  94. Guzzi, M.; Piot, F.; Spinolo, G.; Vedda, A.; Azzoni, C.B.; Paleari, A. Neutron irradiation effects in quartz: Optical absorption and electron paramagnetic resonance. J. Phys. Condens. Matter 1992, 4, 8635–8648. [Google Scholar] [CrossRef]
  95. Griscom, D.; Fowler, W.B. Electron-transfer model for E′-center optical absorption in SiO2. In The Physics of MOS Insulators, Proceedings of the International Topical Conference on the Physics of MOS Insulators; Pergamon Press: Oxford, UK, 1980; pp. 97–101. [Google Scholar]
  96. Trukhin, A.N. Radiation processes in oxygen-deficient silica glasses: Is ODC(I) a precursor of E′-center? J.-Non-Cryst. Solid. 2006, 352, 3002–3008. [Google Scholar] [CrossRef]
  97. Griscom, D.L. Optical properties and structure of defects in silica glass. J. Ceram. Soc. Jpn. 1991, 99, 923–942. [Google Scholar] [CrossRef]
  98. Hosono, H.; Abe, Y.; Imagawa, H.; Imai, H.; Arai, K. Experimental evidence for the Si-Si bond model of the 7.6-eV band in SiO2 glass. Phys. Rev. B 1991, 44, 12043–12045. [Google Scholar] [CrossRef] [PubMed]
  99. Henderson, B.; Chen, Y.; Sibley, W.A. Temperature Dependence of Luminescence of F’ and F Centers in CaO. Phys. Rev. B 1972, 6, 4060–4068. [Google Scholar] [CrossRef]
  100. Perlson, B.D.; Weil, J.A. Electron paramagnetic resonance studies of the E′ centers in alpha-quartz. Can. J. Phys. 2008, 86, 871–881. [Google Scholar] [CrossRef]
  101. Zhang, L.; Leisure, R.G. The Eδ′ and triplet-state centers in x-irradiated high-purity amorphous SiO2. J. Appl. Phys. 1996, 80, 3744–3749. [Google Scholar] [CrossRef]
  102. Afanas’ev, V.; Stesmans, A. Charge state of paramagnetic E′ centre in thermal SiO2 layers on silicon. J. Phys. Condens. Matter 2000, 12, 2285–2290. [Google Scholar] [CrossRef]
  103. Ling, S.; El-Sayed, A.M.; Lopez-Gejo, F.; Watkins, M.B.; Afanas’ev, V.V.; Shluger, A.L. A computational study of Si–H bonds as precursors for neutral E0 centres in amorphous silica and at the Si/SiO2 interface. Microelectron. Eng. 2013, 109, 310–313. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Pristine structures of SiO 2 in different phases, viewed along the a axis. Yellow atoms correspond to Si, red atoms correspond to oxygen atoms. (a) α -cristobalite structure; (b) α -quartz structure. Both phases are on the same scale as α -quartz has a smaller unit cell than α -cristobalite.
Figure 1. Pristine structures of SiO 2 in different phases, viewed along the a axis. Yellow atoms correspond to Si, red atoms correspond to oxygen atoms. (a) α -cristobalite structure; (b) α -quartz structure. Both phases are on the same scale as α -quartz has a smaller unit cell than α -cristobalite.
Materials 16 01382 g001
Figure 2. Schematic representation of oxygen vacancy defects in SiO 2 [8]. Yellow atoms are silicon and red atoms are oxygen, the blue oval represents electron localization. (a) Neutral dimer configuration, such as ODC(I). (b) Positive dimer configuration, such as E d i with elongated Si-Si bonds. (c) Puckered E 1 center, showing the back bond between Si L and the back oxygen in the surrounding ring, O B . The labeling is consistent with α -quartz and is described in the text.
Figure 2. Schematic representation of oxygen vacancy defects in SiO 2 [8]. Yellow atoms are silicon and red atoms are oxygen, the blue oval represents electron localization. (a) Neutral dimer configuration, such as ODC(I). (b) Positive dimer configuration, such as E d i with elongated Si-Si bonds. (c) Puckered E 1 center, showing the back bond between Si L and the back oxygen in the surrounding ring, O B . The labeling is consistent with α -quartz and is described in the text.
Materials 16 01382 g002
Figure 3. Local pristine structures of (a) α -cristobalite with two connected tetrahedrons and (b) α -quartz with two connected tetrahedrons with a ring containing back bonding O B . Yellow atoms are silicon and red atoms are oxygen. In α -quartz, Si A and Si B are not equivalent, but in α -cristobalite they are equivalent. The Si A and Si B notations are, therefore, used to distinguish the Si atoms in defects which are later introduced to the pristine structure. More detail is provided in the text.
Figure 3. Local pristine structures of (a) α -cristobalite with two connected tetrahedrons and (b) α -quartz with two connected tetrahedrons with a ring containing back bonding O B . Yellow atoms are silicon and red atoms are oxygen. In α -quartz, Si A and Si B are not equivalent, but in α -cristobalite they are equivalent. The Si A and Si B notations are, therefore, used to distinguish the Si atoms in defects which are later introduced to the pristine structure. More detail is provided in the text.
Materials 16 01382 g003
Figure 4. Square modulus of the wavefunction of the highest occupied state in the studied dimer oxygen defects in α -C shown in purple, the isosurface shown is | 0.1 | . Yellow atoms are silicon and red atoms are oxygen. (a) V O 0 neutral dimer configuration. (b) E d i , positive dimer. (c) V O , negative dimer. (ac) are general HOMOs applicable to either α -cristobalite or α -quartz.
Figure 4. Square modulus of the wavefunction of the highest occupied state in the studied dimer oxygen defects in α -C shown in purple, the isosurface shown is | 0.1 | . Yellow atoms are silicon and red atoms are oxygen. (a) V O 0 neutral dimer configuration. (b) E d i , positive dimer. (c) V O , negative dimer. (ac) are general HOMOs applicable to either α -cristobalite or α -quartz.
Materials 16 01382 g004
Figure 5. Square modulus of the wavefunction of the highest occupied state in the studied oxygen E defects shown in purple, the isosurface shown is | 0.1 | . Yellow atoms are silicon and red atoms are oxygen. (a) E p u c k , α -quartz puckered, forward projected; (b) E b p , α -cristobalite, back projected; (c) E b p , α -quartz, puckered, back projected.
Figure 5. Square modulus of the wavefunction of the highest occupied state in the studied oxygen E defects shown in purple, the isosurface shown is | 0.1 | . Yellow atoms are silicon and red atoms are oxygen. (a) E p u c k , α -quartz puckered, forward projected; (b) E b p , α -cristobalite, back projected; (c) E b p , α -quartz, puckered, back projected.
Materials 16 01382 g005
Figure 6. Charge transition level diagram of oxygen vacancy defects in (a) α -cristobalite and (b) α -quartz. 0 eV is the valence band maximum and the graph cuts off at the conduction band minimum.
Figure 6. Charge transition level diagram of oxygen vacancy defects in (a) α -cristobalite and (b) α -quartz. 0 eV is the valence band maximum and the graph cuts off at the conduction band minimum.
Materials 16 01382 g006
Table 2. Bond length (Å) between the atoms surrounding the oxygen vacancy: Si A , Si B , O B and O surrounding the Si A / B atoms in α -C and α -Q. For some Si A / B -O bond lengths a range is reported as the bond lengths slightly differ due to asymmetries in the relaxation around the defect. The atoms and labeling conventions are described in Figure 3, Figure 4 and Figure 5. α -cristobalite and α -quartz are compared to understand how they are perturbed between each defect.
Table 2. Bond length (Å) between the atoms surrounding the oxygen vacancy: Si A , Si B , O B and O surrounding the Si A / B atoms in α -C and α -Q. For some Si A / B -O bond lengths a range is reported as the bond lengths slightly differ due to asymmetries in the relaxation around the defect. The atoms and labeling conventions are described in Figure 3, Figure 4 and Figure 5. α -cristobalite and α -quartz are compared to understand how they are perturbed between each defect.
Bond α -Cristobalite α -Quartz
Pristine V O 0 E di E bp V O Pristine V O 0 E di E puck E bp V O
Si A -Si B 3.112.382.884.592.393.092.413.004.455.382.47
Si A -O B 5.315.545.224.835.633.664.103.711.831.844.07
Si A -O1.611.61–1.631.57–1.581.62–1.631.67–1.701.611.631.581.621.61–1.641.69–1.73
Si B -O1.611.62–1.631.57–1.581.54–1.551.67–1.711.611.621.571.581.57–1.591.66–1.68
Table 3. Optical absorption of oxygen vacancies calculated using different exact exchange cutoff radii. 2 Å is the cutoff used in geometry optimization. The values 6.4 Å and 7.5 Å were the highest cutoff points for the phases of α -quartz and α -cristobalite, respectively. Transition types were determined by the state that had the largest contribution to the excitation. In the spin-polarized calculations, i.e., positively charged defects, the α spin channel contributed to a large extent to the transitions. Transition types are notated by the symmetry of the states, the Si. → Si. indicates a promotion into the singly occupied dangling bond on the same Si atom.
Table 3. Optical absorption of oxygen vacancies calculated using different exact exchange cutoff radii. 2 Å is the cutoff used in geometry optimization. The values 6.4 Å and 7.5 Å were the highest cutoff points for the phases of α -quartz and α -cristobalite, respectively. Transition types were determined by the state that had the largest contribution to the excitation. In the spin-polarized calculations, i.e., positively charged defects, the α spin channel contributed to a large extent to the transitions. Transition types are notated by the symmetry of the states, the Si. → Si. indicates a promotion into the singly occupied dangling bond on the same Si atom.
Cutoff Radius (Å)PhaseDefectPeak Energy (eV)Oscillator StrengthTransition Type
7.5 α -C V O 0 7.620.22 σ π
7.720.20
8.130.16
8.40.13
E d i 7.560.12 σ π
7.720.11 σ π *
6.4 α -Q E p u c k 6.180.12Si. → deloc. ring
2 α -C V O 0 7.860.27 σ π
7.940.19
8.200.13
8.480.12
E d i 6.270.15 σ σ *
7.860.12 σ π
α -Q E p u c k 6.380.14Si B → Si A
6.430.12Si. → deloc. ring
Table 4. Calculated isotropic hyperfine coupling constants employing the pseudopotential structures and the all-electron wavefunction with the PBE0 functional, as the quality of the pcJ basis sets systematically improved.
Table 4. Calculated isotropic hyperfine coupling constants employing the pseudopotential structures and the all-electron wavefunction with the PBE0 functional, as the quality of the pcJ basis sets systematically improved.
AE Basis α -Quartz (mT) α -Cristobalite (mT)
E di E puck E di
SilicaSi A Si B Si A Si A Si B
pcJ-012.2910.1142.2511.1610.78
pcJ-111.148.7038.9610.4410.01
pcJ-210.268.0637.469.559.17
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Milton, K.L.; Durrant, T.R.; Cobos Freire, T.; Shluger, A.L. Difference in Structure and Electronic Properties of Oxygen Vacancies in α-Quartz and α-Cristobalite Phases of SiO2. Materials 2023, 16, 1382. https://doi.org/10.3390/ma16041382

AMA Style

Milton KL, Durrant TR, Cobos Freire T, Shluger AL. Difference in Structure and Electronic Properties of Oxygen Vacancies in α-Quartz and α-Cristobalite Phases of SiO2. Materials. 2023; 16(4):1382. https://doi.org/10.3390/ma16041382

Chicago/Turabian Style

Milton, Katherine L., Thomas R. Durrant, Teofilo Cobos Freire, and Alexander L. Shluger. 2023. "Difference in Structure and Electronic Properties of Oxygen Vacancies in α-Quartz and α-Cristobalite Phases of SiO2" Materials 16, no. 4: 1382. https://doi.org/10.3390/ma16041382

APA Style

Milton, K. L., Durrant, T. R., Cobos Freire, T., & Shluger, A. L. (2023). Difference in Structure and Electronic Properties of Oxygen Vacancies in α-Quartz and α-Cristobalite Phases of SiO2. Materials, 16(4), 1382. https://doi.org/10.3390/ma16041382

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop