Abstract

The structural and optical properties of -quartz cluster with oxygen-deficiency centers (ODCs) defects have been investigated based on the density functional theory (DFT). For cluster models with ODC(I) defect, with the increasing of cluster size and shape, the equilibrium length of Si-Si bond decreases. The excitation peaks of cluster models with ODC(I) defect are from 6.87 eV to 7.39 eV, while the excitation peaks of cluster models with ODC(II) defect are from 5.20 eV to 5.47 eV. We also study the interconversion between ODCs (≡Si-Si≡ bond and divalent Si) induced by UV irradiation. Our study predicted the existence of a metastable structure of ODC(I) for the first time in literature. Our results are in good agreement with the previous results and provide strong theoretical support to the viability of the processes.

1. Introduction

As one of the most important crystalline oxides, -quartz is widely used in microelectronics, piezoelectric devices, optical elements, and geological dating, as well as a support in heterogeneous catalysis and other fields of research and technology. Up to now, a large number of experimental and theoretical studies have been devoted to the characterization of the structure of point defects in -quartz [13]. Among many defect centers in -quartz, oxygen-deficiency centers (ODCs) are believed to play a vital role in the irradiation process and subsequent photostructural changes [4] and thus have attracted much attention in research [510].

The crystal structure and symmetry determine that at least two distinct diamagnetic ODCs can form in SiO2 and they are commonly denoted as ODC(I) and ODC(II), giving rise to the photoabsorption bands at ~7.6 and ~5.0 eV, respectively [4, 11]. ODC(I) is a general consensus of the “relaxed oxygen vacancy,” namely, the ≡Si-Si≡ bond, having a Si-Si bond distance about 2.3 Å [4]. As for ODC(II), there has been a large amount of controversy [4, 1214] on the structural model. At least two alternatives were proposed. One is the “unrelaxed oxygen monovacancy,” namely, the ≡Si⋯Si≡ bond, having a Si⋯Si distance similar to that of regular Si-O-Si bonding (3.1 Å) [11], and the other is the “divalent Si,” having two Si-O bonds and a lone pair of electrons in a Si hybrid orbital [1517].

It has been found that the completely unrelaxed geometry with a Si-Si distance of 3.1 Å is not a stable arrangement and is unlikely for the ground state [17]. Indeed, the observed optical properties of ODC(II) can be reproduced almost satisfactorily by the divalent Si model rather than the unrelaxed oxygen vacancy model [1820]. It seems that ODC(II) is more inclined to the divalent Si model. On the other hand, some researchers have shown that even a puckered unrelaxed oxygen vacancy proposed originally for -quartz has a barrier of only ~0.2–0.3 eV against the relaxation into the stable ≡Si-Si≡ bond [21, 22]. It was suggested that an interconversion may occur between the divalent defect and the relaxed oxygen vacancy upon ultraviolet irradiation, although a detailed mechanism of this interconversion is unknown.

In theoretical studies, it is generally accepted that the main intrinsic defect process in SiO2 is the Frenkel mechanism, that is, formation of a pair of an oxygen monovacancy (Si-Si bond) and an interstitial oxygen atom (O0) from a regular Si-O-Si bond. This mechanism has been verified in amorphous SiO2 [2328], where its disordered structure influences both the defect creation and migration. Furthermore, clear experimental evidence for (1) in the ordered lattice of crystalline SiO2 has been given [29]. However, the remaining main unsolved question about ODC(II) is to explain the observed close relationships between ODC(II) and ODC(I). Both centers accommodate the same amount of oxygen deficiency and it is suggested that at some configurations they can interconvert. Details of this process in -quartz are still not unambiguously established [30]

In order to reveal preferable structure of -quartz cluster with ODC(I) defect and the intrinsic mechanism of structural conversion from ODC(I) to ODC(II) in -quartz clusters, we apply the first-principles calculation to different cluster models. The size and shape effect of cluster models is included. Our study predicted the existence of a metastable structure of ODC(I) for the first time in literature. Our results also suggested the mechanism of structural conversion from ODC(I) to ODC(II) in -quartz clusters. We applied the first-principle calculation based on density functional theory (DFT) and time dependent density functional theory (TDDFT) implemented in ORCA program to the electronic structure properties and optical properties. Our paper is organized as follows. In Section 2, the models and methods are described. The results and discussions are presented in Section 3. Finally, a summary of the main conclusions of this work is given in Section 4.

2. Computational Methods

We have studied clusters with various size. The formulas are Si2H2, Si2O6H6, Si5O15H12, Si8O24H18, SiO2H2, and Si3O8H6, as shown in Figure 1.

In order to avoid the spurious states caused by the defects located the boundary of the cluster, we use hydrogen atoms to passivate all clusters. The distance between the capping hydrogen atoms and the oxygen atoms is 0.80 Å [31]. In the process of optimization, the interior atoms were allowed to move freely, while the terminal hydrogen atoms and surface oxygen atoms were frozen. The main geometrical parameters are given for each cluster in Table 1.

The ORCA quantum chemistry program was used in all calculations. DFT and the hybrid density functional B3LYP method [32, 33] were employed to perform the geometry optimization. Following geometry optimization, TDDFT (nroots = 16) was used for the sake of calculating the electronic excited states. Triple--quality basis sets with one set of polarization functions (TZVP) were employed for the Si, O, and H atoms, respectively [34]. In the self-consistent field calculations, tight scf convergence and a dense integration grid were selected. The maximum element of the direct inversion of iterative subspaces (DIIS) error vector is set as . The structural optimization is allowed to be relaxed until the maximum force on each atom becomes less than 0.01 eV/Å and maximum energy change between two steps is smaller than  Eh. The density change is no more than  Eh [35, 36].

3. Results and Discussion

3.1. Preferable Structure of α-Quartz Cluster with ODC(I) Defect

ODC(I) defect pair is created by shifting an interstitial oxygen atom (O0) from a regular Si-O-Si bond and constitutes a pair of an oxygen monovacancy (Si-Si bond) (see (1)). In (1), bonds with three oxygen atoms were represented by ≡. For the purpose of finding out the cluster models which have the stable structure, prediction of the apposite Si-Si band length becomes a key object. The Si-Si band length increased from 2.0 Å to 3.5 Å, which is the main structure change, as shown in Figures 1(a)1(d). Figure 2 shows the relationship between the total energy and Si-Si bond length in singlet state of four clusters with ODC(I) defects. The total energy of the -quartz clusters will have the lowest value along with the increasing of the Si-Si bond length. The cluster size can represent a bulk structure. In order to describe the singlet state adequately, relaxation processes were performed with four cluster models.

To improve the accuracy of the calculations, we consider the parent molecule disilane, Si2H6 cluster model (see Figure 1(a)). The cluster model contains a direct Si-Si bond. As a result of the structural scan and geometry optimization, we have found that, in Figure 2(a), the equilibrium bond length in singlet state of Si2H6 cluster model is 2.34 Å, similar to the ordinary Si-Si bond length of 2.36 Å [31]. However, Figures 2(b)2(d) show that the singlet state equilibrium bond lengths of Si2O6H6, Si5O15H12, and Si8O24H18 cluster models are 2.67 Å, 2.53 Å, and 2.43 Å. With the increasing of cluster size and shape, the equilibrium Si-Si bond length decreases.

3.2. Electronic Structure Properties and Optical Properties of ODC(I)

We compared the electronic density of states (DOS) of the four clusters with ODC(I) defects (see Figure 1) so that the geometric effect on the electronic structure can be seen clearly. Figure 3 shows the total electronic DOS for Si2H6, Si2O6H6, Si5O15H12, and Si8O24H18 clusters. In a nutshell, the electronic DOS of all the clusters shows notable structural sensitivity. The DOS for Si2H6 cluster reveals more uniform distributions than the other three cases, because of the small size. When the valence band states participate in optical transitions of point defects, the structural dependence of electronic state will result in different OA spectra.

The wave function of highest occupied molecular orbital (HOMO) and the electron localization function (ELF) are used to analyze the electronic structure of cluster with ODC(I) defect. We calculate the ELF and the wave function of HOMO on four clusters (see Figures 1(a)1(d)) to clarify the charge distribution around ODC(I) defect. From the electron density, we can approximately derive ELF, introduced by Becke and Edgecombe [37]. The ELF pictures are shown in Figure 4. For the singlet state of the clusters, the ELF was, respectively, computed using the DFT/B3LYP, at the stable geometry of the singlet state. The nature of bonding between Si-Si atoms and charge transfer process is shown in the ELF. The transition can be depicted as an excitation from the top valence band to the ODC(I) defect level in the gap. The excited electron of the ODC(I) defect is located mainly on the central silicon atom. Figures 4(e)4(h) show the wave function of HOMO. The excited electron transfer in the ODC(I) defect can be found in here too.

Table 2 and Figure 5 show the optical properties calculated from four clusters with ODC(I) defects by TDDFT-B3LYP, neglecting other paramagnetic defects. The singlet-to-singlet () excitation energy of Si2H6 cluster was calculated to be 7.39 eV, which agrees well with the previously reported excitation energies (~7.6 eV) [36, 38]. On the other hand, the calculated excitation energy of Si2O6H6, Si5O15H12, and Si8O24H18 clusters was 6.87~7.04 eV. These excitation energy values are only slightly smaller than the experimental transition energy (~7.6 eV), because the basis sets (def2-TZVP) may not be flexible for highly energy. Analyzing these data, the nature of the 7.6 eV absorption band cannot be ascribed to a single point defect in -quartz with ODC(I) defects. It can be explained as a manifestation of the localized states of the disordered structure of silica modified by an oxygen deficit.

3.3. Electronic Structure Properties and Optical Properties of ODC(II)

The ODC(II) defect has been the subject of considerable interest since the discovery of the photorefractive effect in SiO2, and consequently, their microscopic structure and optical properties have been extensively studied using the cluster approach. Figure 6 shows the DOS of the -quartz cluster models with ODC(II) defects computed by DFT-B3LYP. We calculate the total electronic DOS for SiO2H2 and Si3O8H6 clusters (see Figures 1(e) and 1(f)), so that we can explore the geometric effect on the electronic structure. On the whole, the electronic DOS of the cluster shows remarkable structural sensitivity.  We note that the details of the electronic structure are essentially dependent on the -quartz cluster size and shape. Such structural dependence of electronic state will result in different excitation energy when these valence band states participate in optical transition of point defects.

To further illustrate the electronic structure of cluster with ODC(II) defect, we also present the electron localization function (ELF) (see Figures 7(a) and 7(b)) and the wave function (see Figures 7(c) and 7(d)) of HOMO on the two clusters. In order to research the charge distribution around ODC(II) defect, we plotted the ELF. For the singlet state of the clusters, the ELF was, respectively, computed using the DFT/B3LYP. The nature of divalent Si and charge transition process can be obtained in the ELF picture. The transition can be described as the lowest singlet excited state corresponding to a single electron promotion from the Si lone pair orbital to a vacant orbital on the same atom. The excited electron of the ODC(II)defect is located mainly on the central siliconatom. The wave function of HOMO is shown in Figures 7(c) and 7(d). The excited electron transfer in the ODC(II) defect can also be found.

Table 3 shows the optical properties calculated from two clusters with ODC(II) defects by TDDFT-B3LYP, neglecting other paramagnetic defects. transition energies of the SiO2 with ODC(II) defects were calculated by Trukhin [10], Lü et al. [8], and Adelstein et al. [7] using different cluster models and different levels of theory. These studies yield almost the same calculated results concerning . In our work, the transition energies of Si2H6 and Si3O8H6 clusters are, respectively, calculated to be 5.47 eV and 5.20 eV, which agree well with the previously reported excitation energies. Analyzing these data, the nature of the 5.0 eV absorption band can be ascribed to a single electron promotion from the Si lone pair orbital to a vacant orbital on the same atom.

3.4. Photoinduced Interconversions of ODC Defects in α-Quartz

A theoretical study of the interconversion models of amorphous silicon dioxide has been performed by Uchino et al. [39] who proposed a mechanism for the ODC(I) ODC(II) transformation based on photoionization of ODC(I). However, a detailed mechanism of this interconversion in -quartz is unknown. To get better knowledge about the detailed photoinduced interconversion mechanism of ODC defects in -quartz, we here employ quantum-chemical program (ORCA) using Si5O15H12 cluster model.

Starting from the ground state geometry of Si5O15H12 cluster model with ODC(I) (see Figure 1(c)), we have excited the model into the triplet state which is 4.11 eV higher in single point energy [31]. The cluster model relaxes under the modified potential energy surface, the Si-Si disconnects and expands up to 3.02 Å (see Figure 8(a)), which is slightly longer than the equilibrium bond length (2.53 Å) calculated for Figure 1(c), thus inducing a large distortion surrounding the defect. Following the reaction path toward the ODC(II) (see Figure 8(c)), we find the transition state (see Figure 8(b)) and the computed relaxation energy is 0.46 eV higher than the state of Si5O15H12 cluster model with ODC(I). The transition state cannot be stable.

By further decreasing the reaction coordinate, Si5O15H12 cluster transforms spontaneously into the ODC(II) configuration and the computed relaxation energy is 0.30 eV lower than the transition state of Si5O15H12 cluster model. This implies that the transition state configuration corresponds to a metastable minimum and is possible to relax into a lower energy configuration shown in Figure 8(c) during thermal processes. Among the many possible local configurations of the ODC(I) defect in -quartz, only a small fraction is expected to transform into the ODC(II) upon excitation. Although our used cluster model is small, other local configurations different from ours may be more favorable for the interconversion.

4. Conclusions

The DFT method was employed to study the structures of six different cluster models for -quartz with ODCs defects. We predict the apposite Si-Si band length in order to search for the stable structure of cluster models with ODC(I) defects. For Si2O6H6, Si5O15H12, and Si8O24H18 cluster models with ODC(I) defect, with the increasing of cluster size and shape, the equilibrium Si-Si bond length decreases. The OA spectroscopic properties and excitation energies are calculated based on TDDFT. For clusters with ODC(I), the nature of the 7.6 eV absorption band can be explained as a manifestation of the localized states of the disordered structure of silica modified by an oxygen deficit. However, for clusters with ODC(II), the nature of the 5.0 eV absorption band can be ascribed to a single electron promotion from the Si lone pair orbital to a vacant orbital on the same atom. We also study the detailed photoinduced interconversion mechanism of ODCs defects in -quartz. We predicted the existence of a metastable structure of ODC(I) for the first time. Our results provide strong theoretical support to the viability of the processes.

Conflicts of Interest

There are no conflicts of interest related to this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (nos. 61377097, 61675032, and 61671085), the National Basic Research Program of China (973 Program) under Grant no. 2014CB643900, the Open Program of State Key Laboratory of Functional Materials for Informatics, National Natural Science Foundation for Theoretical Physics Special Fund “Cooperation Program” (no. 11547039), Shaanxi Institute of Scientific Research Plan Projects (no. SLGKYQD2-05), and the Fund of State Key Laboratory of Information Photonics and Optical Communications (BUPT) (no. IPOC2015ZT05). The authors acknowledge the computational support from the Beijing Computational Science Research Center (CSRC).