Next Article in Journal
Rice Husk Ash-Based Concrete Composites: A Critical Review of Their Properties and Applications
Next Article in Special Issue
Intensity of Radiative Recombination in the Germanium/Silicon Nanosystem with Germanium Quantum Dots
Previous Article in Journal
Multivariate Analysis Applications to Crystallography
Previous Article in Special Issue
Response of 4H-SiC Detectors to Ionizing Particles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theory of the Thermal Stability of Silicon Vacancies and Interstitials in 4H–SiC

I3N, Department of Physics, University of Aveiro, Campus Santiago, 3810-193 Aveiro, Portugal
Crystals 2021, 11(2), 167; https://doi.org/10.3390/cryst11020167
Submission received: 15 January 2021 / Revised: 2 February 2021 / Accepted: 4 February 2021 / Published: 8 February 2021
(This article belongs to the Special Issue Crystalline Materials for Radiation Detection: A New Perspectives)

Abstract

:
This paper presents a theoretical study of the electronic and dynamic properties of silicon vacancies and self-interstitials in 4H–SiC using hybrid density functional methods. Several pending issues, mostly related to the thermal stability of this defect, are addressed. The silicon site vacancy and the carbon-related antisite-vacancy (CAV) pair are interpreted as a unique and bistable defect. It possesses a metastable negative- U neutral state, which “disproportionates” into V Si + or V Si , depending on the location of the Fermi level. The vacancy introduces a ( / + ) transition, calculated at E c 1.25   eV , which determines a temperature threshold for the annealing of V Si into CAV in n-type material due to a Fermi level crossing effect. Analysis of a configuration coordinate diagram allows us to conclude that V Si anneals out in two stages—at low temperatures ( T 600 °C) via capture of a mobile species (e.g., self-interstitials) and at higher temperatures ( T 1200 °C) via dissociation into V C and C Si defects. The Si interstitial ( Si i ) is also a negative- U defect, with metastable q = + 1 and q = + 3 states. These are the only paramagnetic states of the defect, and maybe that explains why it escaped detection, even in p-type material where the migration barriers are at least 2.7 eV high.

1. Introduction

Silicon carbide (SiC) is currently a mature semiconductor that supports a variety of technologies, most notably in power electronics [1,2]. Among the merits of this material, particularly of its 4H polytype, we find high crystalline quality, a 3.2 eV wide bandgap, a large breakdown field, the possibility to grow both n-type and p-type layers, and exceptional thermal and mechanical stability.
SiC can be grown with a relatively low residual defect concentration ( ~ 10 11 10 13   cm 3 ). Intrinsic defects in SiC have recently received much attention from the community for their negative effects on device performance and for their promising role as building blocks in quantum technologies (e.g., quantum-bit holders, single-photon emitters, or quantum-sensors) [3,4]. Most of these centers involve vacancies; they are usually introduced via electron or ion irradiation, and a precise understanding of their electronic and dynamic properties is of utmost importance, for instance, in order to control their positioning upon dynamic and thermal annealing treatments.
Due to the extreme radiation hardness and the low leakage current of SiC junctions, this material has been proposed for the fabrication of radiation detectors [5,6], including neutrons [7], allowing for operation under harsh conditions, for instance, at high temperatures and under intense radiation fields [8].
During exposure to high-energy radiation or high temperatures, the formation of intrinsic defects can play a decisive role regarding the resulting material specifications. A critical example is the impact of deep carrier traps on the minority carrier lifetime (due to non-radiative recombination at point and extended defects), and therefore on the functioning of bipolar devices.
Defects that are capable of trapping carriers can occur in at least two distinct charge states. These defects are termed acceptors or donors, whether they become negatively or positively charged, respectively. Consequently, trapping defects have transition energy levels, which follow directly from their charge state occupancy against the Fermi level under equilibrium conditions [9]. Traps with transition levels 0.1   eV away from the band edges can be conveniently monitored by deep-level transient spectroscopy (DLTS) or derived techniques [10]. The most conspicuous DLTS signals detected in irradiated high-purity n-type 4H–SiC that are likely to be related to intrinsic defects are Z 1 / 2 , EH 6 / 7 , EH 1 / 3 , EH 4 , EH 5 , M, and S 1 / 2 . They all build up markedly upon high-energy irradiation, although Z 1 / 2 and EH 6 / 7 stand most notably for being present already in as-grown samples with a concentration in the range of 10 12   cm 3 .
Among the above traps, only Z 1 / 2 and EH 6 / 7 gather wide agreement within the community as regard to their origin. Accordingly, the Z 1 / 2 is a two-fold electron emission from double-negative carbon vacancies ( V C ) located on h and k sites ( Z 1 and Z 2 , respectively) of the 4H–SiC lattice [11,12,13,14,15,16,17]. On the other hand, EH 6 / 7 has been assigned to a donor transition of V C . Activation energies for electron emission of Z 1 / 2 and EH 6 / 7 were measured as E a , n = 0.7   eV and 1.5 eV, respectively [11,12,14,16]. V C defects have relatively low formation energy, E f   ~   5   eV , thus explaining their presence in as-grown epi-layers—they are thermally generated and left “frozen” in the layers upon cooling from growth conditions (usually above 1200 °C) [18].
EH 4 and EH 5 traps were recently assigned to a carbon-related antisite-vacancy (CAV) pair [19], an isomer of the silicon vacancy ( V Si ) [20,21,22]. The key arguments for the assignment were that (i) the introduction rates of both EH 4 and EH 5 clearly displayed a linear dependence with the proton fluence, and they were approximately the same as (although slightly and systematically lower than) isolated V C and V Si . This observation suggests single atomic displacements (primary defects) as to the origin of these traps; (ii) the participation of self-interstitials was ruled out based on their low thermal stability as predicted by local-density functional calculations [23]; and (iii) considering that CAV can lie in the lattice along four different orientations, a good match was noted between the calculated CAV ( 0 / + ) levels, around E c 1   eV [22], and the activation energies for electron emission of EH 4 ( E a , n 1.0   eV ) and EH 5 ( E a , n 1.13   eV ). It is noted that EH 5 was connected to one particular alignment of CAV with a calculated ( 0 / + ) level slightly deeper than the other three, all of which accounted for EH 4 .
Although reasonable, the above rationale, in particularly the link of EH 5 to CAV, is still questionable. For instance, if EH 4 and EH 5 are similar defects, essentially V Si isomers, differing in their crystalline orientation, why EH 5 (which allegedly accounts for one of the orientations) is detected when the material is irradiated with electrons whose energy is below the threshold for displacing Si atoms, while EH 4 (representing the other three alignments) is not [24,25,26]? Another puzzle is the fact that EH 4 and EH 5 do not anneal as expected as if they were so similar (identical stoichiometry, structure, and charge state). Whereas the amplitude of EH 4 steadily decreases during 20-min isochronal anneals starting from 300 °C until it vanishes at 1200 °C, EH 5 grows considerably at 400 °C and anneals out at ~ 700 °C [27].
Other traps, namely S 1 / 2 , EH 1 / 3 , and M, have overlapping peak temperatures in the DLTS spectra ( ~ 200   K and ~ 320   K ). Unfortunately, for that reason they have been referred to interchangeably in the literature, thus raising further difficulty in their assignment (see for instance [28]).
S 1 ( E c 0.41   eV ) and S 2 ( E c 0.71   eV ) arise from different charge state transitions of the same defect [29,30]. Their measured capture cross sections at 335   K are 8 × 10 17   cm 2 and 3 × 10 15   cm 2 , thus being consistent with ( q 2 / q 1 ) and ( q 1 / q ) transitions, respectively. The S center anneals out at about 350 °C according to first-order kinetics with an activation energy of E a = 1.8   eV [30]. The frequency of attempt of the process, ~ 1 × 10 11   s 1 , is somewhat low compared to the Debye frequency of SiC, suggesting that it involves a mechanism other than a single atomic jump. S 1 and S 2 have been assigned to V Si based on the calculated energy levels for this defect [14,31], and also from the observed correlation between the magnitude of the S 1 / 2 peaks and the intensity of the V1/V1’ photoluminescence lines (known to arise from transitions of the spin-3/2 state of V Si ) as a function of the proton implantation fluence [32].
Again, this picture is not free of controversy. An obvious question we may ask is how consistent the annealing of S 1 / 2 is with the calculated barriers involving V Si migration and reactions. Hence, during annealing, we either have (i) an unknown defect X that becomes mobile and reacts with V Si , (ii) defect X is immobile and traps a fast-diffusing V Si , or (iii) V Si transforms into a different configuration (e.g., the CAV structure). The fastest among these mechanisms must have a barrier of E a = 1.8   eV . Additionally, in order to comply with the observed first-order kinetics, the concentration of defect X must be at least one order of magnitude higher than that of V Si . For that reason, (i) is hardly consistent with the assignment of X to a fast-diffusing self-interstitial, whereas (ii) and (iii) are not compatible with the large barriers calculated for V Si migration ( E a > 2.6   eV ) and conversion into CAV ( E a > 3.2   eV ) in 4H–SiC [33].
EH 1 ( E c 0.45   eV ) and EH 3 ( E c 0.72   eV ) traps have been detected in 4H–SiC irradiated with low- and high-energy electrons [26,27,34,35,36], protons [30,37], and heavier ions and neutrons [38,39]. In order to clarify whether S 1 / 2 and EH 1 / 3 correspond to the same or distinct defects, Alfieri and Mihaila [36] investigated 116 keV electron irradiated material that had been oxidized in order to suppress the as-grown carbon vacancies. With these low-energy electrons, only carbon-displacements (vacancies, interstitials, and Frenkel pairs) could be introduced. EH 1 and EH 3 annealed out simultaneously at about 350 °C with an activation energy E a = 1.13   eV and 1.17   eV . Their concentration depth profiles corresponded to a one-to-one correlation, thus indicating that they belong to the same defect. The activation energies for annealing are about 0.7   eV lower than those of S 1 and S 2 , thus indicating that EH 1 / 3 and S 1 / 2 are not the same. Based on their observations, Alfieri and Mihaila [36] suggested that EH 1 / 3 is a carbon interstitial related defect.
Due to the injection of interstitials during the oxidation treatment, the concentration of Z 1 / 2 in the pre-irradiated samples of Ref. [36] dropped below the detection limit of DLTS. However, upon irradiation (with a dose of 5 × 10 16   cm 2 ) the introduction of the EH 1 / 3 the peak was about two times more effective than that of Z 1 / 2 . Considering that the latter corresponds to a double emission, the concentration of both EH 1 / 3 and Z 1 / 2 was essentially identical. Although this points to the assignment of EH 1 / 3 to a carbon interstitial, further work is necessary to clarify the introduction efficiency of EH 1 / 3 and Z 1 / 2 , in particular its dose dependency.
The M-center is a bistable defect observed in two configurations (A and B). It becomes evident in the DLTS spectra of electron- and proton-irradiated samples after gentle annealing at ~ 200 °C for 30 mins [40]. Such treatments clean up the crowded spectra of the as-irradiated samples, and it is well possible that the defect is already formed at this stage. Configuration A is formed when the sample temperature is T > 20 °C and reverse bias is applied. This process is thermally activated ( E a = 0.90   eV ) and corresponds to a B A transformation. On the other hand, A B takes place at T > 140 °C with no bias, and the activation energy is now E a = 1.40   eV [41]. Configuration A is responsible for two DLTS peaks, M 1 ( E c 0.42   eV ) and M 3 ( E c 0.83   eV ), while configuration B has one peak referred to as M 2 ( E c 0.63   eV ) [40,41]. The M-center anneals out in the range 310 370 °C according to first-order kinetics with an activation energy of E a = 2.0   eV [41]. The introduction of the M-center is also reported in a paper of this special issue [39], where samples were both ion-implanted and neutron-irradiated, adding further support for its connection to an intrinsic defect.
The bulk of previous theoretical work on migration and reaction mechanisms involving intrinsic defects in SiC was carried out for the simpler 3C allotrope (see for instance [23]). Although elucidative, we now know that these results do not necessarily apply to 4H–SiC [42,43]. In addition, the narrower gap of 3C or the use of too aggressive approximations, such as the local or semi-local treatments of the many-body electronic interactions, implies that they should be revised using more accurate methods. The same applies to results obtained for 4H– and 6H–SiC [33,42,43,44,45]. This issue was recently investigated by Yan et al [46] using hybrid density functional theory (DFT), where it was found a much stronger charge state dependence of migration barriers of intrinsic defects in 4H–SiC than previously thought [23]. For instance, the migration barrier of V C increased from about 2.5 eV in the double-negative charge state to 5.5 eV in the double-positive charge state. Another interesting aspect is the anisotropic effect found by Bathen et al. [47], showing both theoretically and experimentally that the annealing of V C proceeds preferentially via migration of the defect along the basal plane of the 4H–SiC crystal.
The present paper attempts to consolidate some of the most recent findings and to revisit some of the open issues highlighted above, in particular by analyzing the charge-state dependence of transformation, dissociation, and migration of Si-related intrinsic defects in 4H–SiC. The data were obtained from highly accurate hybrid density functional calculations. After a short description of the methodology (Section 2), results are reported for the Si vacancy and self-interstitial (Section 3). In both cases, charge-state-driven bistability and migration mechanisms are discussed (Section 4).

2. Theoretical Methods

The electronic structure calculations reported below were carried out within range-separated hybrid density functional theory using the Vienna Ab-initio Simulation Package (VASP) [48,49,50,51]. The many-body exchange-correlation interactions were accounted for using the functional proposed by Heyd, Scuseria, and Ernzerhof (commonly named HSE06 functional) [52]. The projector-augmented wave method was used to treat the core electrons [53], whereas valence states were expanded in plane-waves with a cut-off energy of 420 eV. The latter was solved self-consistently until the total energy between two consecutive steps differed by less than 10 7 eV.
Defects were inserted in 400-atom supercells of 4H–SiC with a hexagonal shape, constructed by replication of 5 × 5 × 2 primitive cells. The calculated basal and axial lattice constants were a = 3.0712   Å and c = 10.051   Å , respectively. Forces acting on atomic nuclei were obtained within the generalized gradient approximation (GGA) to the exchange-correlation potential [54]. They were minimized with respect to the atomic coordinates until the maximum residual force was lower than 0.01 eV/Å. For the calculation of forces and structural relaxation (GGA-level), the Brillouin zone (BZ) was sampled on a Γ -centered 2 × 2 × 2   k -point grid ( Γ - 2 3 ). On the other hand, for total energy calculations (HSE06-level), the sampling was done at the Γ -point (see [47,55] for details).
The formation energy of a Si intrinsic defect in a state described by its net charge and atomistic structure ( q , R ) is obtained as a function of the Fermi energy ( E F ) using the usual expression,
E f ( q , R ; E F ) = E def ( q , R ) E bulk ± μ Si + q ( E v + E F ) ,
where E def is the calculated energy of the defective supercell, E bulk is the energy of a bulk (defect-free) supercell, μ Si is the chemical potential of a Si atom removed from (plus sign) or inserted into (negative sign) the bulk supercell in order to obtain a vacancy or an interstitial, and E v is the valence band maximum energy, here taken as the highest occupied Kohn-Sham state of a bulk calculation. The chemical potential of silicon was calculated for a carbon-rich (silicon-poor) crystal, μ Si = μ Si 0 + Δ H SiC 0 , where μ Si 0 is the energy per silicon atom in pure silicon and Δ H SiC 0 = 0.61   eV the calculated heat of formation of 4H–SiC (to be compared to 0.72   eV from calorimetry measurements [56]). We note that E def includes a charge correction, calculated here according to the recipe of Freysoldt, Neugebauer, and Van de Walle [57].
Migration and transformation mechanisms were investigated using the climbing image nudged elastic band method (CI-NEB) [58]. While already reported elsewhere [47], a summary of the main workflow steps are as follows: the first step consisted of postulating a mechanism, which started and finished at specific minimum-energy configurations (previously obtained). Essentially, it formed a sequence of seven intermediate structures between the end-configurations. On the second step, an exploratory NEB-constrained calculation was performed using the Γ -point for sampling the BZ. This allowed for a relatively fast screening of the configurational space, close to the minimum energy path between the end-structures. In the third step, the relaxed image sequence from the previous step was further refined by turning on the climbing-image mode of the NEB algorithm and improving the forces by increasing the BZ sampling to Γ - 2 3 . Finally, a self-consistent calculation within HSE06-level was carried out using the transition state geometry obtained from the previous step.
With the above specifications, quantities derived from the energy of (local) minimum structures (e.g., formation energies, transition levels) have an estimated numerical accuracy of about 10 meV. On the other hand, transition state calculations (activation energies for migration and structural transformation) have a numerical accuracy of 0.1 eV [47,55].

3. Results

Before moving on to the results, let us introduce a few useful words about the notation employed throughout. A defect state is referred to as D   S   q ( R ) , where D stands for a stoichiometric label ( V Si or Si i for a silicon atom subtracted or added to a perfect crystal, respectively), q = { , 1 ,   0 ,   + 1 , } is a charge state, S = { 0 , 1 / 2 , 1 , } the total spin, and R is a structural label to indicate a specific atomistic geometry. For instance, V   3 / 2 Si 1 ( k ) represents the lack of a pseudo-cubic ( k ) Si atom in 4H–SiC, with a trapped electron (net charge –1) and total spin S = 3 / 2 .
Should a generalization be made, or when the context makes the notation unambiguous, we may reduce the above to D   S   q , D   q , D or simply R . Some defect structures can be described as a combination of elementary defects (e.g., antisites and vacancies). The complex pair formed by first-neighboring carbon-antisite and carbon-vacancy (CAV) can adopt up to four different alignments, depending on the lattice location of its elements. One of its alignments could be represented as C Si ( k ) V C ( h ) . However, because its stoichiometry is identical to that of a silicon vacancy, we simply write V Si ( CAV , k h ) , where first and second lattice indexes refer to the antisite and vacancy sites (not interchangeable).

3.1. The Silicon Vacancy

We investigated several configurations for V Si in 4H–SiC, of which the most stable are represented in Figure 1. The left-most structure represents a pristine piece of crystal for the sake of reference. White and dark grey spheres represent Si and C atoms, respectively. The first configuration is the site vacancy structure, V Si ( k ) , which is simply the lack of a Si atom, in this case, removed from the pseudo-cubic site (represented by a translucent white sphere). The second structure is V Si ( CAV ) , a carbon antisite-vacancy first neighboring pair (in this case, with an off-axis k h alignment). This structure “communicates” with V Si ( k ) via jumping of the C Si atom into the V C site. Another stable structure investigated, and also shown in Figure 1, is V Si ( CAV 2 ) . It is also a carbon antisite-vacancy pair, but V C is on the second neighboring C-shell with respect to C Si (both lying on pseudo-cubic sub-lattice sites).
The relative stability of different V Si configurations are reported in Table 1. For each net charge of the defect, q , zero energy indicates the most stable state considering both atomistic and spin configurations. In agreement with previous results [22], the simple site configurations ( { k , h } ) are more stable in the negative charge states, while the CAV structure is more favorable when positively charged. Both configurations have four unsaturated radicals, and their electronic structure share an identical shape, i.e., [ a 1 ] a 1 e for V Si 0 within the C 3 v point group (the singlet state within brackets is resonant with the valence band, while the other states are located in the bandgap). However, while site vacancies have four nearly equivalent and highly electronegative carbon radicals with states in the lower part of the gap, the CAV defect has a doublet localized on three nearly equivalent Si atoms with states high in the gap. Hence, the deep C-radicals of { k , h } site structures can trap more electrons than CAV, whereas the latter can trap more holes. It is noted that some defect structures mentioned above, notably V   1 Si 0 ( k ) ( { k , h } ) , V   1 / 2 Si 3 ( k ) ( { k , h } ) , V   1 Si 0 ( CAV , { k k , h h } ) and V   1 / 2 Si 1 ( CAV , { k k , h h } ) do not retain the full C 3 v symmetry and display small monoclinic Jahn–Teller distortions of the order of tens of meV. In this case, the doublet splits into close a ′ and a ″ singlets.
The site vacancy in 4H–SiC was found to be stable in charge states from 0 to 3 . The strong localization of the C-radicals favors high spin states [59], namely a 1 e , a 1 e , a 1 e and a 1 e , for q = 0 , 1 , 2 and 3 , respectively, where we dropped the resonant [ a 1 ] level for the sake of simplicity. Conversely, a 1 0 e 0 , a 1 e 0 , a 1 e , and a 1 e were found for charge states q = + 2 , + 1 , 0 , and 1 of CAV, respectively, where the spin S = 1 of the neutral state originates from a very small overlap of the a 1 level (strongly localized on the unique C atom) and the e level (localized on the Si atoms and nodal on the C atom).
The CAV2 configuration, comprising a more remote C Si V C defect pair (see right-hand side of Figure 1) can be thought of a carbon vacancy perturbed by the tensile strain field of a nearby inert C Si defect. We actually found that like V C [14,16], CAV2 is stable in charge states from + 2 to 2 with a negative- U ordering of the acceptor levels. However, CAV2 is always metastable with respect to both CAV and { k , h } configurations by about 1   eV and it is not expected to be found under equilibrium conditions.
The energy of a neutral uncorrelated pair of C Si and V C defects were also calculated as 1.1–1.2 eV higher than that of V Si 0 ( CAV , h h ) , depending on the sub-lattice occupation of the defect components. This figure is representative of the binding energy of CAV. We finally note that V Si structures involving Si antisites, e.g., C Si Si C V Si obtained by displacing a Si atom into V C in CAV (see Figure 1) were not stable and relaxed back to CAV.

3.1.1. Bistability of the Silicon Vacancy

We now look at the relative stability of the several vacancy states in equilibrium at a specific temperature T . Figure 2a depicts the formation energy for the relevant states as obtained from Equation (1). Formation energies for site and CAV structures are shown in red and blue, respectively. Solid and dashed lines represent the energies of defects located on hexagonal and cubic sites, respectively. There is a total of four possible alignments for CAV. However, we only considered the axial defects, V Si ( CAV , h h ) and V Si ( CAV , k k ) . The lower horizontal axis in Figure 2a shows the Fermi energy with respect to the valence band top, E F E v . The upper horizontal axis marks the location of the Fermi level in the material (assumed to be n-type), for the temperatures indicated. For that, a donor concentration N D = 10 14   cm 3 was assumed. Accordingly, for an extrinsic semiconductor [60],
E F ( T ) = E i ( T ) + k B T   log ( n n i ) ,
where E i is the intrinsic Fermi level,
E i ( T ) = E v + E c 2 + 1 2 k B T   log [ g v g c ( m h * m e * ) 3 / 2 ] ,
n is the free electron density,
n = N D 2 + ( N D 2 ) 2 + n i 2 ,
and n i is the intrinsic carrier density,
n i = N c N v   exp ( E g 2 k B T ) .
In the above, E c , v are the energies of the conduction band bottom and valence band top denoted by c and v subscripts, respectively. These correspond to a density of states N c , v = 2 g c , v   ( 2 π   m e , h *   k B T / h 2 ) 3 / 2 , where g c , v are band edge degeneracy factors ( g c = 3 and g v = 1 in 4H–SiC), k B is the Boltzmann constant and h the Plank constant. In the density of states and Equation (3), m e , h * are the free electron and hole effective mass.
The results of Figure 2a are in line with previous reports using hybrid density functional theory [22]. While the silicon vacancy adopts the CAV structure in p-type material, it displays the “simpler” site configuration in n-type material. The defect is amphoteric and strongly compensating. It can capture up to three electrons in n-type 4H–SiC, whereas in p-type material it captures up to two holes. The site configuration has three acceptor transitions at E ( / 0 ) = E c 2.13 ± 0.02   eV , E ( = / ) = E c 0.74 ± 0.13   eV and E ( 3 / 2 ) = E c 0.33 ± 0.01   eV . These are average values obtained from the calculated levels of cubic and hexagonal configurations (c.f. Table 2). On the other hand, the CAV configuration has donor levels at E ( + / + + ) = E c 1.98 ± 0.03   eV and E ( 0 / + ) = E c 1.03 ± 0.01   eV , and one acceptor transition at E ( / 0 ) = E c 0.44 ± 0.03   eV .
The formation energy of the V Si + 2 ( CAV ) state is rather low in p-type material (about 3–4 eV). Of course, this is still high enough to prevent the thermal generation of these defects at room temperature. However, the formation of CAV defects should be preferred in detriment of the site configuration, for instance in irradiated or ion-implanted p-type SiC. This was actually predicted experimentally by Konstantinov [61], who postulated the formation of antisite-type defects in the context of the formation of shallow and deep boron acceptors in Si and C sites, respectively.
It is noteworthy from Figure 2a that the neutral vacancy is not stable (under thermodynamic conditions), irrespective of is configuration and location of the Fermi level. The reason is that the following disproportionation reaction is exothermic.
2 V Si 0 ( CAV ) V Si + ( CAV ) + V Si ( CAV ) V Si + ( CAV ) + V Si
Although the energy of the intermediate state with respect to reactants is positive, U = 0.63   eV (which is the Hubbard correlation energy of the neutral CAV state, in this case for the k k alignment), this figure is lower in magnitude than the energy release of the final step due to the V Si ( CAV ) V Si reconfiguration. In the end the effective U = 0.46   eV . This feature is what distinguishes a negative- U defect [9], thus implying the existence of a ( / + ) transition, here calculated at E c 1.25 ± 0.02   eV .
From the upper horizontal axis of Figure 2a we arrive at the conclusion that in n-type material and above ~ 600 °C the CAV structure becomes the most favorable state (positively charged). This suggests that an annealing/quenching process could be used in order to create a non-equilibrium population of CAV defects in n-type 4H–SiC. We note that although being exothermic, the kinetics of Reaction (6) strongly depends on the emission/capture rates of carriers, as well as on transformation barriers along the process. This will be discussed further below. In p-type 4H–SiC the site structure of V Si is invariably unstable against CAV.
The left-hand side of Figure 3 depicts a configuration coordinate diagram of the silicon vacancy in n-type 4H–SiC. It shows the potential along the transformation coordinate V Si ( k ) V Si ( CAV , k k ) only, although that will suffice for our needs. The barriers involving transformations between V Si ( h ) and V Si ( CAV , h h ) systematically exceeded those of Figure 3 by ~ 0.5   eV . Barriers for all charge states are indicated on the trail of the arrows along the potential energy surfaces. Vertical arrows indicate the energy difference between the respective energy minima. Note that V Si 1 and V Si 1 ( CAV ) have a different spin in their ground states ( S = 3 / 2 and 1 / 2 , respectively). This means that a spin-flipping event will take place somewhere along the conversion V 1 / 2 Si 1 ( CAV ) V 3 / 2 Si 1 . Our calculations assume that single atomic jumps occur on a much shorter time scale than that involving the change of the spin state [62]. Hence, transformations are expected to conserve the spin of the initial state.
The V Si V Si ( CAV ) transformation essentially involves a jump of a carbon atom from the edge of the Si vacancy into the empty site. This may occur either in the neutral or negative charge states. Only the former leads to an exothermic reaction. Other charge states, namely V Si 3 and V Si 2 must emit electrons into the conduction band or capture holes from the valence band before making the jump. The barriers for jumping into the CAV structure are rather high, namely 4.3 eV and 3.6 eV for charge states q = 1 and q = 0 , respectively.
Another distinct feature of a negative- U defect is the double emission (or double capture) of carriers during transients upon changing the Fermi level. For instance, and according to Figure 3, if starting from the negative charge state in the k configuration in a reverse-biased n-type diode, the ionization into the most stable neutral state, which adopts the CAV structure, involves a net change in the energy of Δ E ( / 0 ) = E c E ( / 0 ) = 1.51   eV . If the temperature of the sample is high enough for this step to proceed within an observable timescale, an additional ionization into charge state q = + 1 will follow at an even faster rate because it involves an electron emission with Δ E ( 0 / + ) = 1.04   eV only. Of course, the experimental verification of this effect depends on the intermediate transformation barriers, and as we will see further ahead, on the competition with the migration and annealing out of the defect itself.
The thermodynamic drive that controls the V Si V Si ( CAV ) reaction direction is the Fermi energy. From Figure 2, we find that in n-type 4H–SiC and up to T ~ 600 °C, the vacancy is negatively charged ( q = 3 ,   2 , or 1 ) in the site configuration. Above this temperature, the positive CAV state is more favorable. Supposing that the Fermi level is located at mid-gap and the vacancy is in the V   3 / 2 Si 1 (non-equilibrium) state, Figure 3 tells us that the CAV structure is only reachable if a jump over a barrier of about 4.3   eV is performed. Alternatively, it could first emit an electron into the conduction band (with an activation energy of about 2.2   eV ), followed by a jump in the neutral charge state over a 3.6   eV barrier. Neither of these alternatives is likely to occur below T ~ 1000 °C in a timescale of seconds (compare with the migration of the carbon vacancy, which is observed above 1200 °C and has a measured barrier of 3.6   eV [47]).
On the other hand, if a V Si + ( CAV ) defect is found in n-type 4H–SiC without an external bias, it will be “eager” for electrons. First of all, it will quickly capture two electrons and become negatively charged, V 1 / 2 Si 1 ( CAV ) . Further energy lowering is only possible by overcoming a 3.3 eV barrier via V 1 / 2 Si 1 ( CAV ) V 3 / 2 Si 1 . After that, additional electrons can be captured depending on the Fermi level.
Finally, it is noted that previous semi-local calculations gave a qualitatively similar picture for the silicon vacancy, including the acceptor character and the preference for the CAV configuration in p-type material [59]. However, and despite its merits, the semi-local approach used to describe the exchange-correlation functional lacks accuracy, making it difficult to draw quantitative conclusions. For instance, the calculated electronic transition levels reported in [59] suffer from the same insufficiency that underestimates the bandgap width by about 40%. By employing the experimental bandgap in the construction of the formation energy diagram, the calculated acceptor levels came out too deep (with respect to E c ) and actually, a ( 4 / 3 ) transition was artificially predicted. For the same reason, the calculated V Si V Si ( CAV ) conversion barriers were about 30–50% smaller than those reported here.

3.1.2. Migration of the Silicon Vacancy

As reported by Bockstedte et al. [63], migration of a silicon vacancy in SiC must occur via exchange with Si second neighbors. The right-hand side of Figure 3 depicts the calculated migration barriers of the vacancy in 4H–SiC within the cubic sub-lattice. All charge states are included. Solid and dashed lines represent potential energy surfaces for migration along the basal and axial directions, respectively. We note that each axial mechanism involves a total of four consecutive jumps. Values in Figure 3 correspond to the energy of the highest saddle point along the way with respect to the initial and final states. Two features are readily noticeable—(i) the migration barrier increases as the defect emits electrons and becomes less negative and (ii) basal migration barriers are invariably lower than their axial counterparts (for the same charge state). This suggests that, like the carbon vacancy, diffusion of silicon vacancies could be highly anisotropic.
Before jumping to the conclusions, we should bear in mind several effects that compete with the jumps of the vacancy, namely, electron emission and capture, in addition to conversion into CAV and subsequent dissociation into a remote C Si V C pair via migration of V C . We know that V C migrates only at T 1200   ° C , at which point the Fermi level is close to the middle of the gap, and for which the ground state of the Si vacancy is actually V Si + ( CAV ) . From here, it could capture a free-electron (releasing 1 eV) and convert it into the metastable V Si 0 species by overcoming a 4.2 eV high barrier and finally perform a migration jump (over a 4 eV barrier). Considering V Si + ( CAV ) as the zero-energy reference, the overall energy barrier for migration could be as low as 3.7 eV if the energy that results from the electron capture was not dissipated via phonon emission, but rather converted into kinetic energy and assist the Si nuclei to perform the migration jump. However, if the energy drop from the capture is dissipated, the migration barrier increases to 4.7 eV. It is also noted that at least two kinetic effects are expected to hinder V Si migration. Firstly, the capture of a carrier depends on their availability and cross section, and secondly, the V Si 0 V Si 0 ( CAV ) return barrier is considerably lower than the migration barrier.
Alternatively, and again assuming that the Fermi energy is close to mid-gap (high temperatures), V Si + ( CAV ) can convert into V Si + ( CAV 2 ) (show in Figure 1), and from there V C could migrate in the positive or neutral charge state. Considering that the energy of uncorrelated neutral C Si and V C defects are about 1 eV higher than V Si ( CAV ) , and adding this to the 3.6 eV migration barrier of V C [47], we arrive at a dissociation barrier for the CAV configuration of about 4.6 eV.

3.2. Silicon Interstitials

Now we turn to the silicon self-interstitials ( Si i ). Several structures were investigated, among which we have Si atoms at empty sites (cage structures), Si–Si and Si–C split interstitials. The cage structures are denoted by the sublattice layer where the interstitial Si atom is located. Four of these sites are represented graphically in Figure 4 (left). Two types of Si–Si split structures were found, namely those with the Si–Si dimer direction possessing a component along the hexagonal axis of the crystal (referred to as axial or A-structures), and those whose Si–Si dimer is parallel to the basal plane of the crystal (referred to as basal or B-structures). Examples of these are depicted in the middle and right-hand side of Figure 4, respectively. Si–C split structures were found to be unstable and relaxed to cage configurations. These findings were already reported in previous studies [46,63,64].
Split-interstitial structures A and B are also divided into two groups, namely those occupying k and h sublattice Si sites. Figure 4 shows Si i ( A , h ) and Si i ( B , k ) configurations for the sake of exemplification.
Table 3 reports the relative stability of the lower energy Si i defects in 4H–SiC for several charge states. As it will be shown below, acceptor states are not stable, and therefore only neutral and positive charge states are reported. The energy of the most stable configuration in each charge state is set to zero. Energies for k and h configurations are not included in the table. These structures are either unstable (in q = 0 and + 1 charge states) or their relative energies are way too high ( E > 2.5   eV for q = + 2 ,   + 3 , and + 4 states) with respect to other competing structures. All cage structures retain the maximum site symmetry ( C 3 v ), except the Si i 0 ( h ) state, which undergoes a distortion upon symmetry-unconstrained relaxation, possibly due to a pseudo-Jahn–Teller effect. For this reason, it is referred as h * , and its relative energy in Table 3 is also labeled with a star.
We can readily draw a few conclusions from Table 3. Overall, cage structures are more stable in positive charge states, while split-interstitial forms are more favorable in the neutral state. The Si i defect almost exclusively inhabits the h layer of 4H–SiC. The exception is the Si i + 4 state, which stabilizes in the k configuration in p-type material. As the defect becomes more positively charged (moving from left to right in Table 3), the Si i ( A , h ) state becomes less stable in favor of Si i ( h ) , and the latter becomes less stable in favor of Si i ( k ) . We can rationalize these findings in the following way: the fewer electrons bound to the Si–Si split defect, the less stable the bonds with ligands become, and with further ionization, the internal crystal field across the large h -cage tends to stabilize interstitial Si cations towards the more electronegative carbon anion in the 4H–SiC lattice (upper left region of Figure 4).
Blank cells in Table 3 mean that the respective states are not stable. Hence, cage structures can emit up to four electrons into the conduction band bottom (or capture four holes from the valence band top). On the other hand, split interstitial structures can only emit up to two electrons. The Si i 2 + ( A , h ) state is not stable and relaxes into the Si i 2 + ( h ) cage form.

3.2.1. Bistability of the Silicon Interstitial

We can get a better perspective of the above results by looking at the right-hand side of Figure 2. Here we depict a diagram with formation energies of the cage (blue lines), split-axial (red lines), and split-basal (green lines) configurations as a function of the Fermi level (bottom horizontal axis). Solid lines and dashed lines represent the formation energies of defects located at h and k sublattice sites. The black line labeled h * stands for the neutral distorted Si i 0 ( h * ) as described above. The energy of the non-distorted trigonal Si i 0 ( h ) configuration is 0.9   eV higher than that of Si i 0 ( h * ) . The top horizontal axis represents the Fermi level location at the temperatures indicated in a sample with a free-electron density of n = 1 × 10 14   cm 3 .
In general, the formation energy of the Si vacancies is considerably lower than that of Si self-interstitials. This is particularly noticeable in n-type material and especially relevant under intrinsic conditions (high temperatures) as it has an important impact on self-diffusion. With the Fermi level at mid-gap, Si vacancies and interstitials are most stable in the positive V Si + ( CAV ) and double-positive Si i + 2 ( h ) states, respectively. Hence, mutual repulsion will hider their annihilation should any of them become mobile. However, in n-type material in the range of T 400 600 °C, the vacancy adopts the site structure V Si and could be an efficient trap for mobile Si i + 2 defects.
In n-type 4H–SiC the Si interstitial adopts the split Si i 0 ( A , h ) state (solid red horizontal line in Figure 2). Among competing configurations, at about 0.3 eV higher in the energy scale we find the split-basal structures and Si i 0 ( h * ) . On the other hand, in p-type 4H–SiC, Si i 4 + ( h ) is the most favorable state. Hence, while Si vacancies compensate n-type material, the deep donors of the Si interstitial are expected to strongly compensate p-type 4H–SiC as well.
The carbon self-interstitial has two negative- U metastable charge states, respectively; q = + 1 and q = + 3 . These are only expected to form under non-equilibrium conditions, e.g., under illumination or thermal excitation. The effect is clearly illustrated in Figure 5, where a configuration coordinate diagram is presented for the Si i defect in n-type 4H–SiC. Only states that are relevant for transformations and migration of Si i are included in the diagram. Split-basal configurations are metastable in the q = 0 , + 1 , and + 2 charge states.
In equilibrated n-type material, the defect is found in the Si i 0 ( A , h ) state. If the defect is provided with at least 1.45   eV (thermally or optically) so that it becomes ionized as Si i 0 ( A , h ) Si i + ( A , h ) + e , a subsequent transformation into a slightly lower energy state, Si i + ( h ) , followed by further emission Si i + ( h ) Si i 2 + ( h ) + e will take place immediately. This is because the transformation barrier is only 0.2 eV high and the second ionization costs 0.5 eV, which must less than the first ionization. Likewise, hole emission from Si i 4 + ( k ) Si i 3 + ( k ) + h + involves the absorption of about 1.6   eV (which is the position of the E ( + 3 / + 4 ) E v transition involving the cage interstitial in the k site). However, the ionized Si i 3 + ( k ) species is not the ground state anymore, and an eventual transformation into the lower energy Si i 3 + ( h ) state leaves the defect in a state with a hole ionization energy of only 0.75 eV (which is the position of the E ( + 2 / + 3 ) E v transition involving the cage interstitial in the h site). It is noted that in this case the structural transformation involves overcoming a large barrier of 2.1 eV, and as we will see, that is higher than the activation energy for migration in this particular charge state.

3.2.2. Migration of the Silicon Interstitial

The reporting of migration mechanisms and barriers of silicon interstitials is now carried out for n-type and intrinsic conditions. Results for p-type conditions are presented at the end of this section. Similar to the vacancy, we assume that the interstitials do not exchange carriers with the crystalline states in the course of elementary jumps. This may only occur when the defect “lands” at intermediate (meta-)stable states along the mechanisms, whose lifetime is much longer than the charge carrier emission/capture periods. Figure 2 (right) suggests that up to about 400   ° C in n-type 4H–SiC, the Si i defect is expected to adopt the Si i 0 ( A , h ) state. The lowest energy path for migration of this species was found to be
Si i 0 ( A , h ) 1.54   eV Si i 0 ( h * ) 1.21   eV Si i 0 ( B , k ) 1.19   eV Si i 0 ( A , h ) ,
which accounts for both axial and basal migration. This is because the first and final configurations can differ by a basal lattice vector. Along the hexagonal direction, they differ by half lattice vector only. However, the remaining half path is identical to the first half by symmetry if we consider a rotation of the Si–Si dimer in Si i 0 ( A , h ) by 2 π / 3 (passing by Si i 0 ( B , h ) ), which has a barrier of 1.3   eV . The transition-state energies indicated over the forward arrows in Reaction (7) are all relative to the Si i 0 ( A , h ) ground state. Overall, the migration barrier is E a , m = 1.54   eV (limited by the first step), and this is summarized on the left-hand side of the configuration coordinate diagram of Figure 5.
In intrinsic material (with the Fermi level close to mid-gap) the Si i + 2 ( h ) state is the most favorable. In this case, we found that “cage-jumping” leads to the lowest transition state energies. In this case, the mechanism is
Si i 2 + ( h ) 2.63   eV Si i 2 + ( k ) 2.69   eV Si i 2 + ( h ) ,
which accounts for axial and basal migration as well. The end-configurations differ by one basal lattice vector and half axial vector. The saddle point is actually a structure close to the Si i 2 + ( k ) state. Again, energies over reaction barriers are relative to the initial (final) ground state. A simplified potential energy surface for the mechanism is depicted in the upper left part of Figure 5.
As shown in Figure 2, in p-type material the ground state is Si i + 4 ( k ) . For these conditions, the migration mechanism is still of cage-jumping type and it is now
Si i 4 + ( k ) 3.71   eV Si i 4 + ( h ) 3.56   eV Si i 4 + ( k ) .
Overall the migration barrier is E a , m = 3.71   eV high. The first step involves a purely axial jump of the Si i 4 + ion and it is the limiting one. The second step, in addition to an axial component, involves the displacement of Si i 4 + by a basal lattice vector. Hence, the mechanism accounts for both basal and axial migration.

4. Discussion

The calculated E ( 3 / 2 ) and E ( = / ) levels of V Si in the site configuration at about E c 0.33   eV and E c 0.74   eV , are rather close to the depth of S 1 and S 2 traps measured at E c 0.41   eV and E c 0.71   eV , respectively. The large barriers that were found for the migration of V Si in different charge states (>3 eV) and those for the conversion into the CAV form (>3.5 eV) are not compatible with the observed first-order annealing kinetics with activation energy E a = 1.8   eV for S 1 / 2 . Such low activation energy suggests an annihilation mechanism due to the capture of mobile self-interstitials by the vacancies. The figure is not far from the calculated migration barrier of neutral Si i defects ( 1.54   eV ). However, in order to account for the first-order kinetics, the concentration of interstitials has to be much larger than that of vacancies. One possibility is that carbon interstitials (which are easier to be displaced from the lattice in irradiated material) are the main culprits for the annihilation of V Si during anneals at above 300 °C.
In n-type material, at higher temperatures ( T 600 800 °C), any available V Si defect that escaped annihilation at lower temperatures, is now more stable in the positive CAV form. This suggests that an annealing stage driven by a Fermi level effect should occur in this temperature range.
According to Figure 3, migration of V Si could occur via the following mechanism,
V Si + ( CAV ) + e ( 1 ) V Si 0 ( CAV ) ( 2 ) V Si 0 ( 3 ) V Si 0 ,
which is strongly hampered by the need of capturing an electron (or emitting a hole) in step 1, plus the fact that the barrier for the back-reaction of step 2 is lower than the forward barrier (3.58 eV against 4.22 Ev; c.f. Figure 3). The overall barrier of Reaction 10 is estimated as high as E a , m = 4.7   eV . Hence, considering the binding energy of C Si and V C of 1 eV with respect to CAV, plus the migration barrier of V C (3.6 eV), it is likely that the high-temperature stage for the annealing of the Si vacancy occurs via dissociation of CAV with an activation barrier of about E a , d = 4.6   eV . Annealing via capture of a mobile species (e.g., self-interstitials from dissolving interstitial aggregates) is also possible.
The CAV defect has calculated donor levels at E ( 0 / + ) = E c 1.03   eV and E ( + / + + ) = E c 1.98   eV . It also has a calculated acceptor transition at E ( / 0 ) = E c 0.44   eV . As recently suggested by Karsthof et al. [19], the first donor level is rather close to the EH4 and EH5 electron traps. Other potential DLTS signals that could be assigned to CAV ( + / + + ) and CAV ( 0 / + ) transitions are HK4 and EM1 [65]. These are hole and electron traps found in p-type 4H–SiC with levels at E v + 1.56   eV (HK4) and E v + 2.26   eV (EM1), respectively. Considering the calculated band gap of 3.3 eV, we place the calculated CAV donor levels at E ( + / + + ) = E v + 1.32   eV and E ( 0 / + ) = E v + 2.27   eV . HK4 and EM1 can actually be found in as-grown p-type 4H–SiC with a concentration of the order of 10 12   cm 3 [65], thus being consistent with the low formation energy of CAV in the p-type material. Considering their relative locations in the gap, it is possible that EM1 (minority carrier trap in p-type material) and EH4 (majority carrier trap in n-type material) could originate from the same defect transition. In any case, we are still left with an additional CAV level of acceptor type, predicted at E c 0.44   eV (which should be easy to detect DLTS), but for which there is no experimental match with compatible intensity and thermal stability to the other peaks.
Interestingly, in n-type 4H–SiC, silicon vacancies have several acceptor states, whereas interstitials have several donor states. Their multiple opposite charges are expected to drive a strong Coulomb attraction, and therefore to promote annihilation. On the other hand, in p-type and intrinsic material, both vacancies (with the CAV structure) and interstitials are positively charged, and annihilation should be hindered.
During epitaxial growth and higher temperatures ( T 1400 °C), the material is intrinsic. Under these conditions, the formation energy of the Si vacancy is about 6.2 eV and that of the Si interstitial is about 8.7 eV. Assuming that (i) V Si ( CAV ) dissociates before being able to migrate in the form of a site vacancy and (ii) the migration barrier of Si i + 2 is E a , m = 2.7   eV , we arrive at an estimate for the activation energy for self-diffusion of silicon, which is solely due to the motion of silicon interstitials, of E a , sd Si = 11.4   eV . Within this picture, V Si defects do not contribute to self-diffusion. This exceeds the measured analogous quantity for carbon self-diffusion by about 4 eV [66,67], and it is consistent with the observation that silicon diffusivity is more than two orders of magnitude lower than that of carbon [68,69,70].
Regarding the migration of neutral Si interstitials in n-type 4H–SiC, the proposed mechanism is in line with that found for the migration of Si i in 3C–SiC by Bockstedte et al. [63], where a ⟨ 110 ⟩-kick-out sequence was found to be the most favorable path. It is noted that a strictly identical sequence is not possible in 4H–SiC due to the lower symmetry of the host crystal. That is why the Si i 0 ( h * ) state (which is not far from a split configuration), appears along the way. In Ref. [63], the kick-out migration mechanism of the neutral split-interstitial in 3C–SiC was reported with an activation barrier E a , m = 1.4   eV . Our results indicate that in 4H–SiC the barrier is about 0.1–0.2 eV higher.
For p-type 3C–SiC, the authors of Ref. [63] found a migration barrier for Si i 4 + of E a , m = 3.5   eV . Again, this figure is about 0.2 eV lower than the analogous barrier found in this work for 4H–SiC. In a more recent study, Yan et al. [46] found migration barriers for the Si i defect lower than 0.9   eV for the neutral state and about 1.8   eV for the q = + 4 charge state. These results, which employed methodologies identical to those in this paper, imply a considerably faster diffusivity and lower temperature for annealing of Si i . Unfortunately, the details provided in Ref. [46] regarding the mechanisms involved are rather limited, and therefore it is rather difficult to provide a comment other than suggesting that this issue needs to be revisited.
The formation of interstitial clusters in SiC is an important issue. Hornos et al. [71] found that silicon interstitials can aggregate and form stable and electrically active complexes. They suggested that at high temperatures, the Si clusters can dissolve and emit interstitials with significant consequences to the formation of other defects. Our calculated migration barriers indicate that the kinetics of these processes should be much faster in n-type than in p-type material.

5. Conclusions

This paper presents a theoretical study of the transformation and migration mechanisms of Si vacancies and self-interstitials in 4H–SiC using first-principles methods. Particular attention was given to the investigation of charge state effects on the potential barriers.
Considering that the vacancy has two distinct configurations (site and CAV structures) in n-type and p-type material, it is shown that the neutral state is metastable, irrespectively of the position of the Fermi level. That is, V Si has a negative- U in the neutral state, giving rise to a ( / + ) transition estimated at E c 1.25 ± 0.02   eV . This level is likely to define the threshold for the annealing of V Si into the CAV structure in n-type material due to a displacement of the Fermi level deep into the gap. V Si is also an amphoteric center that strongly compensates both p- and n-type material.
The calculated electronic transitions and thermal stability of CAV are consistent with the DLTS emission energy and annealing of EH 4 . Accordingly, EH4 is stable up to 1200   ° C , which is interpreted as the migration of V C away from C Si . The dissociation barrier of CAV is estimated as E a , d = 4.6   eV . Si vacancies are unlikely to anneal out via migration of the site structure. That route was shown to be hindered by the need for the emission or capture of a carrier. Additionally, at high temperatures (intrinsic conditions), Si vacancies and interstitials are most stable in the positive V Si + ( CAV ) and double-positive Si i + 2 ( h ) states, respectively. Hence, mutual repulsion will hamper their annihilation should any of them become mobile. However, in n-type 4H–SiC in the range of T 400 600 °C, the vacancy adopts the site structure V Si and could be an efficient trap for mobile Si i + 2 defects as well as carbon interstitials.
The Si I defect shows a strong charge-state dependent migration barrier. In n-type material the defect is neutral, and the mechanism consists of a sequence of kick-out jumps involving Si–Si split structures with an overall barrier of 1.5 eV. On the other hand, in intrinsic and p-type material the stable charge states are q = + 2 and q = + 4 , respectively. In that case, the migration mechanism involves a jump of Si i between neighboring cages with barriers of 2.7 eV and 3.7 eV, respectively.
We find that the Si i is a negative- U defect with metastable q = + 1 and q = + 3 states. These are the only paramagnetic states of the defect and that could explain why it escaped detection by electron paramagnetic resonance, even in p-type material where the migration barrier is at least 2.7 eV high. Of course, other techniques are in principle sensitive to the presence of Si i in 4H–SiC, including those related to the absorption or luminescence involving electronic transition within gap states.

Funding

The author thanks the support of the i3N projects, Refs. UIDB/50025/2020 and UIDP/50025/2020, financed by the Fundação para a Ciência e a Tecnologia in Portugal. This work was supported by the Science for Peace and Security NATO Program through project SPS G5674.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kimoto, T. Material Science and Device Physics in SiC Technology for High-Voltage Power Devices. Jpn. J. Appl. Phys. 2015, 54, 040103. [Google Scholar] [CrossRef]
  2. She, X.; Huang, A.Q.; Lucia, O.; Ozpineci, B. Review of Silicon Carbide Power Devices and Their Applications. IEEE Trans. Ind. Electron. 2017, 64, 8193–8205. [Google Scholar] [CrossRef]
  3. Castelletto, S.; Rosa, L.; Johnson, B.C. Silicon Carbide for Novel Quantum Technology Devices. In Advanced Silicon Carbide Devices and Processing; Saddow, S.E., Via, F.L., Eds.; InTech: London, UK, 2015; ISBN 978-953-51-2168-8. [Google Scholar] [CrossRef] [Green Version]
  4. Lohrmann, A.; Johnson, B.C.; McCallum, J.C.; Castelletto, S. A Review on Single Photon Sources in Silicon Carbide. Rep. Prog. Phys. 2017, 80, 034502. [Google Scholar] [CrossRef] [PubMed]
  5. Ruddy, F.H.; Seidel, J.G.; Chen, H.; Dulloo, A.R.; Ryu, S.-H. High-Resolution Alpha-Particle Spectrometry Using Silicon Carbide Semiconductor Detectors. In Proceedings of the IEEE Nuclear Science Symposium Conference Record, Fajardo, Puerto Rico, 23–29 October 2005; Volume 3, pp. 1231–1235. [Google Scholar]
  6. Coutinho, J.; Torres, V.J.B.; Capan, I.; Brodar, T.; Ereš, Z.; Bernat, R.; Radulović, V.; Ambrožič, K.; Snoj, L.; Pastuović, Ž.; et al. Silicon Carbide Diodes for Neutron Detection. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2021, 986, 164793. [Google Scholar] [CrossRef]
  7. Dulloo, A.R.; Ruddy, F.H.; Seidel, J.G.; Adams, J.M.; Nico, J.S.; Gilliam, D.M. The Thermal Neutron Response of Miniature Silicon Carbide Semiconductor Detectors. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2003, 498, 415–423. [Google Scholar] [CrossRef]
  8. Szalkai, D.; Ferone, R.; Issa, F.; Klix, A.; Lazar, M.; Lyoussi, A.; Ottaviani, L.; Tutto, P.; Vervisch, V. Fast Neutron Detection With 4H-SiC Based Diode Detector up to 500 °C Ambient Temperature. IEEE Trans. Nucl. Sci. 2016, 63, 1491–1498. [Google Scholar] [CrossRef]
  9. Coutinho, J.; Markevich, V.P.; Peaker, A.R. Characterisation of Negative-U Defects in Semiconductors. J. Phys. Condens. Matter 2020, 32, 323001. [Google Scholar] [CrossRef] [PubMed]
  10. Peaker, A.R.; Markevich, V.P.; Coutinho, J. Tutorial: Junction Spectroscopy Techniques and Deep-Level Defects in Semiconductors. J. Appl. Phys. 2018, 123, 161559. [Google Scholar] [CrossRef] [Green Version]
  11. Dalibor, T.; Pensl, G.; Matsunami, H.; Kimoto, T.; Choyke, W.J.; Schöner, A.; Nordell, N. Deep Defect Centers in Silicon Carbide Monitored with Deep Level Transient Spectroscopy. Phys. Status Solidi (A) 1997, 162, 199–225. [Google Scholar] [CrossRef]
  12. Hemmingsson, C.; Son, N.T.; Kordina, O.; Bergman, J.P.; Janzén, E.; Lindström, J.L.; Savage, S.; Nordell, N. Deep Level Defects in Electron-Irradiated 4H SiC Epitaxial Layers. J. Appl. Phys. 1997, 81, 6155–6159. [Google Scholar] [CrossRef]
  13. Hemmingsson, C.G.; Son, N.T.; Ellison, A.; Zhang, J.; Janzén, E. Negative-U Centers in 4H Silicon Carbide. Phys. Rev. B 1998, 58, R10119–R10122. [Google Scholar] [CrossRef]
  14. Hornos, T.; Gali, A.; Svensson, B.G. Large-Scale Electronic Structure Calculations of Vacancies in 4H-SiC Using the Heyd-Scuseria-Ernzerhof Screened Hybrid Density Functional. MSF 2011, 679–680, 261–264. [Google Scholar] [CrossRef] [Green Version]
  15. Son, N.T.; Trinh, X.T.; Løvlie, L.S.; Svensson, B.G.; Kawahara, K.; Suda, J.; Kimoto, T.; Umeda, T.; Isoya, J.; Makino, T.; et al. Negative-U System of Carbon Vacancy in 4H-SiC. Phys. Rev. Lett. 2012, 109, 187603. [Google Scholar] [CrossRef] [Green Version]
  16. Coutinho, J.; Torres, V.J.B.; Demmouche, K.; Öberg, S. Theory of the Carbon Vacancy in 4H-SiC: Crystal Field and Pseudo-Jahn-Teller Effects. Phys. Rev. B 2017, 96, 174105. [Google Scholar] [CrossRef] [Green Version]
  17. Capan, I.; Brodar, T.; Pastuović, Ž.; Siegele, R.; Ohshima, T.; Sato, S.; Makino, T.; Snoj, L.; Radulović, V.; Coutinho, J.; et al. Double Negatively Charged Carbon Vacancy at the H- and k-Sites in 4H-SiC: Combined Laplace-DLTS and DFT Study. J. Appl. Phys. 2018, 123, 161597. [Google Scholar] [CrossRef] [Green Version]
  18. Ayedh, H.M.; Nipoti, R.; Hallén, A.; Svensson, B.G. Elimination of Carbon Vacancies in 4H-SiC Employing Thermodynamic Equilibrium Conditions at Moderate Temperatures. Appl. Phys. Lett. 2015, 107, 252102. [Google Scholar] [CrossRef]
  19. Karsthof, R.; Bathen, M.E.; Galeckas, A.; Vines, L. Conversion Pathways of Primary Defects by Annealing in Proton-Irradiated n-Type 4H-SiC. Phys. Rev. B 2020, 102, 184111. [Google Scholar] [CrossRef]
  20. Umeda, T.; Son, N.T.; Isoya, J.; Janzén, E.; Ohshima, T.; Morishita, N.; Itoh, H.; Gali, A.; Bockstedte, M. Identification of the Carbon Antisite-Vacancy Pair in 4H-SiC. Phys. Rev. Lett. 2006, 96, 145501. [Google Scholar] [CrossRef] [Green Version]
  21. Steeds, J.W. Photoluminescence Study of the Carbon Antisite-Vacancy Pair in 4H- and 6H-SiC. Phys. Rev. B 2009, 80, 245202. [Google Scholar] [CrossRef]
  22. Szász, K.; Ivády, V.; Abrikosov, I.A.; Janzén, E.; Bockstedte, M.; Gali, A. Spin and Photophysics of Carbon-Antisite Vacancy Defect in 4H Silicon Carbide: A Potential Quantum Bit. Phys. Rev. B 2015, 91, 121201. [Google Scholar] [CrossRef] [Green Version]
  23. Bockstedte, M.; Mattausch, A.; Pankratov, O. Ab Initio Study of the Annealing of Vacancies and Interstitials in Cubic SiC: Vacancy-Interstitial Recombination and Aggregation of Carbon Interstitials. Phys. Rev. B 2004, 69, 235202. [Google Scholar] [CrossRef] [Green Version]
  24. Sasaki, S.; Kawahara, K.; Feng, G.; Alfieri, G.; Kimoto, T. Major Deep Levels with the Same Microstructures Observed in N-Type 4H–SiC and 6H–SiC. J. Appl. Phys. 2011, 109, 013705. [Google Scholar] [CrossRef] [Green Version]
  25. Kawahara, K.; Thang Trinh, X.; Tien Son, N.; Janzén, E.; Suda, J.; Kimoto, T. Investigation on Origin of Z1/2 Center in SiC by Deep Level Transient Spectroscopy and Electron Paramagnetic Resonance. Appl. Phys. Lett. 2013, 102, 112106. [Google Scholar] [CrossRef] [Green Version]
  26. Storasta, L.; Bergman, J.P.; Janzén, E.; Henry, A.; Lu, J. Deep Levels Created by Low Energy Electron Irradiation in 4H-SiC. J. Appl. Phys. 2004, 96, 4909–4915. [Google Scholar] [CrossRef]
  27. Alfieri, G.; Monakhov, E.V.; Svensson, B.G.; Linnarsson, M.K. Annealing Behavior between Room Temperature and 2000 °C of Deep Level Defects in Electron-Irradiated n-Type 4H Silicon Carbide. J. Appl. Phys. 2005, 98, 043518. [Google Scholar] [CrossRef]
  28. Beyer, F.C.; Hemmingsson, C.; Pedersen, H.; Henry, A.; Janzén, E.; Isoya, J.; Morishita, N.; Ohshima, T. Annealing Behavior of the EB-Centers and M-Center in Low-Energy Electron Irradiated n-Type 4H-SiC. J. Appl. Phys. 2011, 109, 103703. [Google Scholar] [CrossRef] [Green Version]
  29. Åberg, D.; Storasta, L.; Hallén, A.; Svensson, B.G. Implantation Temperature Dependent Deep Level Defects in 4H-SiC. MSF 2001, 353–356, 443–446. [Google Scholar] [CrossRef]
  30. David, M.L.; Alfieri, G.; Monakhov, E.M.; Hallén, A.; Blanchard, C.; Svensson, B.G.; Barbot, J.F. Electrically Active Defects in Irradiated 4H-SiC. J. Appl. Phys. 2004, 95, 4728–4733. [Google Scholar] [CrossRef]
  31. Gali, A. Defects in SiC: Theory. MSF 2011, 679–680, 225–232. [Google Scholar] [CrossRef]
  32. Bathen, M.E.; Galeckas, A.; Müting, J.; Ayedh, H.M.; Grossner, U.; Coutinho, J.; Frodason, Y.K.; Vines, L. Electrical Charge State Identification and Control for the Silicon Vacancy in 4H-SiC. NPJ Quantum Inf. 2019, 5, 111. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, X.; Zhao, M.; Bu, H.; Zhang, H.; He, X.; Wang, A. Formation and Annealing Behaviors of Qubit Centers in 4H-SiC from First Principles. J. Appl. Phys. 2013, 114, 194305. [Google Scholar] [CrossRef]
  34. Castaldini, A.; Cavallini, A.; Rigutti, L.; Pizzini, S.; Le Donne, A.; Binetti, S. Diffusion Length and Junction Spectroscopy Analysis of Low-Temperature Annealing of Electron Irradiation-Induced Deep Levels in 4H-SiC. J. Appl. Phys. 2006, 99, 033701. [Google Scholar] [CrossRef]
  35. Hazdra, P.; Vobecký, J. Radiation Defects Created in N-Type 4H-SiC by Electron Irradiation in the Energy Range of 1–10 MeV. Phys. Status Solidi A 2019, 216, 1900312. [Google Scholar] [CrossRef]
  36. Alfieri, G.; Mihaila, A. Isothermal Annealing Study of the EH1 and EH3 Levels in N-Type 4H-SiC. J. Phys. Condens. Matter 2020, 32, 465703. [Google Scholar] [CrossRef]
  37. Pastuović, Ž.; Siegele, R.; Capan, I.; Brodar, T.; Sato, S.; Ohshima, T. Deep Level Defects in 4H-SiC Introduced by Ion Implantation: The Role of Single Ion Regime. J. Phys. Condens. Matter 2017, 29, 475701. [Google Scholar] [CrossRef] [Green Version]
  38. Capan, I.; Brodar, T.; Yamazaki, Y.; Oki, Y.; Ohshima, T.; Chiba, Y.; Hijikata, Y.; Snoj, L.; Radulović, V. Influence of Neutron Radiation on Majority and Minority Carrier Traps in N-Type 4H-SiC. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms 2020, 478, 224–228. [Google Scholar] [CrossRef]
  39. Brodar, T.; Bakrač, L.; Capan, I.; Ohshima, T.; Snoj, L.; Radulović, V.; Pastuović, Ž. Depth Profile Analysis of Deep Level Defects in 4H-SiC Introduced by Radiation. Crystals 2020, 10, 845. [Google Scholar] [CrossRef]
  40. Martin, D.M.; Kortegaard Nielsen, H.; Lévêque, P.; Hallén, A.; Alfieri, G.; Svensson, B.G. Bistable Defect in Mega-Electron-Volt Proton Implanted 4H Silicon Carbide. Appl. Phys. Lett. 2004, 84, 1704–1706. [Google Scholar] [CrossRef]
  41. Nielsen, H.K.; Hallén, A.; Svensson, B.G. Capacitance Transient Study of the Metastable M Center in n-Type 4 H-Si C. Phys. Rev. B 2005, 72, 085208. [Google Scholar] [CrossRef]
  42. Defo, R.K.; Zhang, X.; Bracher, D.; Kim, G.; Hu, E.; Kaxiras, E. Energetics and Kinetics of Vacancy Defects in 4 H-SiC. Phys. Rev. B 2018, 98, 104103. [Google Scholar] [CrossRef] [Green Version]
  43. Rauls, E.; Lingner, T.; Hajnal, Z.; Greulich-Weber, S.; Frauenheim, T.; Spaeth, J.-M. Metastability of the Neutral Silicon Vacancy in 4H-SiC. Phys. Status Solidi (B) 2000, 217, R1–R3. [Google Scholar] [CrossRef]
  44. Torpo, L.; Marlo, M.; Staab, T.E.M.; Nieminen, R.M. Comprehensive Ab Initio Study of Properties of Monovacancies and Antisites in 4H-SiC. J. Phys. Condens. Matter 2001, 13, 6203–6231. [Google Scholar] [CrossRef]
  45. Wiktor, J.; Jomard, G.; Bertolus, M. Electronic Structure Calculations of Positron Lifetimes in SiC: Self-Consistent Schemes and Relaxation Effect. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms 2014, 327, 63–67. [Google Scholar] [CrossRef]
  46. Yan, X.; Li, P.; Kang, L.; Wei, S.-H.; Huang, B. First-Principles Study of Electronic and Diffusion Properties of Intrinsic Defects in 4H-SiC. J. Appl. Phys. 2020, 127, 085702. [Google Scholar] [CrossRef]
  47. Bathen, M.E.; Coutinho, J.; Ayedh, H.M.; Hassan, J.U.; Farkas, I.; Öberg, S.; Frodason, Y.K.; Svensson, B.G.; Vines, L. Anisotropic and Plane-Selective Migration of the Carbon Vacancy in SiC: Theory and Experiment. Phys. Rev. B 2019, 100, 014103. [Google Scholar] [CrossRef] [Green Version]
  48. Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 47, 558–561. [Google Scholar] [CrossRef] [PubMed]
  49. Kresse, G.; Hafner, J. Ab Initio Molecular-Dynamics Simulation of the Liquid-Metal–Amorphous-Semiconductor Transition in Germanium. Phys. Rev. B 1994, 49, 14251–14269. [Google Scholar] [CrossRef] [PubMed]
  50. Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef]
  51. Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  52. Heyd, J.; Scuseria, G.E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207–8215, Erratum in J. Chem. Phys. 2006, 124, 219906, doi:10.1063/1.2204597. [Google Scholar] [CrossRef] [Green Version]
  53. Blöchl, P.E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Gouveia, J.D.; Coutinho, J. Can We Rely on Hybrid-DFT Energies of Solid-State Problems with Local-DFT Geometries? Electron. Struct. 2019, 1, 015008. [Google Scholar] [CrossRef] [Green Version]
  56. Greenberg, E.; Natke, C.A.; Hubbard, W.N. The Enthalpy of Formation of Silicon Carbide by Fluorine Bomb Calorimetry. J. Chem. Thermodyn. 1970, 2, 193–201. [Google Scholar] [CrossRef]
  57. Freysoldt, C.; Neugebauer, J.; Van de Walle, C.G. Fully Ab Initio Finite-Size Corrections for Charged-Defect Supercell Calculations. Phys. Rev. Lett. 2009, 102, 016402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Henkelman, G.; Uberuaga, B.P.; Jónsson, H. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901–9904. [Google Scholar] [CrossRef] [Green Version]
  59. Mattausch, A. Ab Initio-Theory of Point Defects and Defect Complexes in SiC. Doctoral Thesis, Friedrich-Alexander-Universität (FAU), Erlangen-Nürnberg, Germany, 2005. [Google Scholar]
  60. Sze, S.M.; Ng, K.K. Physics of Semiconductor Devices; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2006; ISBN 978-0-470-06832-8. [Google Scholar]
  61. Konstantinov, A.O. Nature of the Point Defects Generated during the Diffusion of Acceptor Impurities in Silicon Carbide. Soviet Phys. Semicond. 1992, 26, 151. [Google Scholar]
  62. Harvey, J.N. Understanding the Kinetics of Spin-Forbidden Chemical Reactions. Phys. Chem. Chem. Phys. 2007, 9, 331–343. [Google Scholar] [CrossRef]
  63. Bockstedte, M.; Mattausch, A.; Pankratov, O. Ab Initio Study of the Migration of Intrinsic Defects in 3C-SiC. Phys. Rev. B 2003, 68, 205201. [Google Scholar] [CrossRef] [Green Version]
  64. Liao, T.; Roma, G.; Wang, J. First-Principles Study of Neutral Silicon Interstitials in 3C- and 4H-SiC. Philos. Mag. 2009, 89, 2271–2284. [Google Scholar] [CrossRef] [Green Version]
  65. Alfieri, G.; Kimoto, T. Minority Carrier Transient Spectroscopy of As-Grown, Electron Irradiated and Thermally Oxidized p-Type 4H-SiC. MSF 2014, 778–780, 269–272. [Google Scholar] [CrossRef]
  66. Linnarsson, M.K.; Janson, M.S.; Zhang, J.; Janzén, E.; Svensson, B.G. Self-Diffusion of 12C and 13C in Intrinsic 4H-SiC. J. Appl. Phys. 2004, 95, 8469–8471. [Google Scholar] [CrossRef]
  67. Rüschenschmidt, K.; Bracht, H.; Stolwijk, N.A.; Laube, M.; Pensl, G.; Brandes, G.R. Self-Diffusion in Isotopically Enriched Silicon Carbide and Its Correlation with Dopant Diffusion. J. Appl. Phys. 2004, 96, 1458–1463. [Google Scholar] [CrossRef]
  68. Hon, M.H.; Davis, R.F.; Newbury, D.E. Self-Diffusion of Si-30 in Polycrystalline Beta-SiC. J. Mater. Sci. 1980, 15, 2073–2080. [Google Scholar] [CrossRef]
  69. Hong, J.D.; Davis, R.F. Self-Diffusion of Carbon-14 in High-Purity and N-Doped Alpha-SiC Single Crystals. J. Am. Ceram. Soc. 1980, 63, 546–552. [Google Scholar] [CrossRef]
  70. Bathen, M.E.; Linnarsson, M.; Ghezellou, M.; Ul Hassan, J.; Vines, L. Influence of Carbon Cap on Self-Diffusion in Silicon Carbide. Crystals 2020, 10, 752. [Google Scholar] [CrossRef]
  71. Hornos, T.; Son, N.T.; Janzén, E.; Gali, A. Theoretical Study of Small Silicon Clusters in 4H-SiC. Phys. Rev. B 2007, 76, 165209. [Google Scholar] [CrossRef]
Figure 1. Ball-and-stick models of a bulk region of 4H–SiC and three structures of a silicon vacancy, namely the k -site vacancy (k), the first neighboring carbon antisite-vacancy (CAV) pair, and a second neighboring carbon antisite-vacancy (CAV2) pair. Silicon and carbon atoms are shown in white and gray, respectively. Missing atoms are semi-transparent.
Figure 1. Ball-and-stick models of a bulk region of 4H–SiC and three structures of a silicon vacancy, namely the k -site vacancy (k), the first neighboring carbon antisite-vacancy (CAV) pair, and a second neighboring carbon antisite-vacancy (CAV2) pair. Silicon and carbon atoms are shown in white and gray, respectively. Missing atoms are semi-transparent.
Crystals 11 00167 g001
Figure 2. Formation energy as a function of the Fermi level (with respect to the valence band top) for (a) the silicon vacancy and (b) the silicon self-interstitial in 4H–SiC. For the vacancy, site (red) and CAV (blue) structures are represented. For the interstitial, cage (blue and black), split axial (red), and split basal (green) are represented. Solid and dashed lines represent the energies of defects located on hexagonal and cubic sites. The upper axis of each diagram indicates the position of the Fermi level for a specific temperature, considering a free-carrier (n-type) density of n = 1 × 10 14   cm 3 .
Figure 2. Formation energy as a function of the Fermi level (with respect to the valence band top) for (a) the silicon vacancy and (b) the silicon self-interstitial in 4H–SiC. For the vacancy, site (red) and CAV (blue) structures are represented. For the interstitial, cage (blue and black), split axial (red), and split basal (green) are represented. Solid and dashed lines represent the energies of defects located on hexagonal and cubic sites. The upper axis of each diagram indicates the position of the Fermi level for a specific temperature, considering a free-carrier (n-type) density of n = 1 × 10 14   cm 3 .
Crystals 11 00167 g002
Figure 3. Configuration coordinate diagram silicon vacancy in n-type 4H–SiC on the pseudo-cubic site. The transformation between site- and CAV-configurations is shown along the c -axis (left-hand side). Migration barriers for different charge states are depicted on the right-hand side. Here, solid and dashed lines refer to vacancy motion along the basal and axial directions, respectively. All values are in eV.
Figure 3. Configuration coordinate diagram silicon vacancy in n-type 4H–SiC on the pseudo-cubic site. The transformation between site- and CAV-configurations is shown along the c -axis (left-hand side). Migration barriers for different charge states are depicted on the right-hand side. Here, solid and dashed lines refer to vacancy motion along the basal and axial directions, respectively. All values are in eV.
Crystals 11 00167 g003
Figure 4. Ball-and-stick models of silicon self-interstitials in 4H–SiC. On the left we indicate the location of the four stable cage forms ( k , k , h , and h ), whereas the middle- and right-structures depict two examples of split forms, specifically with (pseudo-)axial and basal aligned Si–Si dimers on h and k sites, respectively. Silicon and carbon atoms are shown in white and gray, respectively.
Figure 4. Ball-and-stick models of silicon self-interstitials in 4H–SiC. On the left we indicate the location of the four stable cage forms ( k , k , h , and h ), whereas the middle- and right-structures depict two examples of split forms, specifically with (pseudo-)axial and basal aligned Si–Si dimers on h and k sites, respectively. Silicon and carbon atoms are shown in white and gray, respectively.
Crystals 11 00167 g004
Figure 5. Configuration coordinate diagram silicon self-interstitial in n-type 4H–SiC. Only the most stable configurations are shown. The left-hand side of the diagram shows the migration mechanism of positive charge states. The right-hand side of the diagram shows the migration mechanism of the neutral charge state. In the q = + 1 potential curve, solid and dashed lines refer to vacancy motion along the basal and axial directions, respectively. Anisotropy was not found for migration barriers in other charge states (see text). All values are in eV.
Figure 5. Configuration coordinate diagram silicon self-interstitial in n-type 4H–SiC. Only the most stable configurations are shown. The left-hand side of the diagram shows the migration mechanism of positive charge states. The right-hand side of the diagram shows the migration mechanism of the neutral charge state. In the q = + 1 potential curve, solid and dashed lines refer to vacancy motion along the basal and axial directions, respectively. Anisotropy was not found for migration barriers in other charge states (see text). All values are in eV.
Crystals 11 00167 g005
Table 1. Energies (in eV) of different silicon-vacancy structures in 4H–SiC, within the same charge state, with respect to the corresponding ground state (zero energy). Site ( k , h ) and axial CAV configurations ( k k , h h ) were considered. The total spin of each state is given between parentheses.
Table 1. Energies (in eV) of different silicon-vacancy structures in 4H–SiC, within the same charge state, with respect to the corresponding ground state (zero energy). Site ( k , h ) and axial CAV configurations ( k k , h h ) were considered. The total spin of each state is given between parentheses.
Structure q = 3 q = 2 q = 1 q = 0 q = + 1 q = + 2
k 0.00 (1/2)0.00 (1)0.02 (3/2)0.73 (1)
h 0.23 (1/2)0.24 (1)0.00 (3/2)0.67 (1)
CAV, k k 1.12 (1/2)0.09 (1)0.12 (1/2)0.05 (0)
CAV, h h 0.96 (1/2)0.00 (1)0.00 (1/2)0.00 (0)
CAV2, k k 2.50 (0)2.08 (1/2)1.16 (0)1.63 (1/2)1.18 (0)
CAV2, k h 2.09 (0)1.72 (1/2)0.91 (0)1.50 (1/2)1.37 (0)
Table 2. Depth of the transition levels with respect to the conduction band minimum, Δ E ( q 1 / q ) = E c E ( q 1 / q ) , of silicon-vacancy structures in 4H–SiC. Site ( k , h ) and axial CAV configurations ( k k , h h ) were considered. All values are in eV.
Table 2. Depth of the transition levels with respect to the conduction band minimum, Δ E ( q 1 / q ) = E c E ( q 1 / q ) , of silicon-vacancy structures in 4H–SiC. Site ( k , h ) and axial CAV configurations ( k k , h h ) were considered. All values are in eV.
Structure Δ E ( 3 / 2 ) Δ E ( = / ) Δ E ( / 0 ) Δ E ( 0 / + ) Δ E ( + / + + )
k 0.330.882.15
h 0.340.612.11
CAV, k k 0.411.041.95
CAV, h h 0.481.012.01
Table 3. Energies (in eV) of different silicon self-interstitials in 4H–SiC, within the same charge state, with respect to the corresponding ground state (zero energy). Cage-like ( { k , k , h } ), split-axial (A, { k , h } ) and split-basal (B, { k , h } ) structures were considered. Empty cells indicate that the respective states were unstable. All states have a low-spin configuration ( S = 0 and S = 1 / 2 for even and odd charge states, respectively).
Table 3. Energies (in eV) of different silicon self-interstitials in 4H–SiC, within the same charge state, with respect to the corresponding ground state (zero energy). Cage-like ( { k , k , h } ), split-axial (A, { k , h } ) and split-basal (B, { k , h } ) structures were considered. Empty cells indicate that the respective states were unstable. All states have a low-spin configuration ( S = 0 and S = 1 / 2 for even and odd charge states, respectively).
Structure q = 0 q = + 1 q = + 2 q = + 3 q = + 4
h 0.32 *0.000.000.001.49
k 2.541.341.510.210.00
A, h 0.000.07
A, k 1.150.731.07
B, h 0.370.612.43
B, k 0.240.532.26
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Coutinho, J. Theory of the Thermal Stability of Silicon Vacancies and Interstitials in 4H–SiC. Crystals 2021, 11, 167. https://doi.org/10.3390/cryst11020167

AMA Style

Coutinho J. Theory of the Thermal Stability of Silicon Vacancies and Interstitials in 4H–SiC. Crystals. 2021; 11(2):167. https://doi.org/10.3390/cryst11020167

Chicago/Turabian Style

Coutinho, José. 2021. "Theory of the Thermal Stability of Silicon Vacancies and Interstitials in 4H–SiC" Crystals 11, no. 2: 167. https://doi.org/10.3390/cryst11020167

APA Style

Coutinho, J. (2021). Theory of the Thermal Stability of Silicon Vacancies and Interstitials in 4H–SiC. Crystals, 11(2), 167. https://doi.org/10.3390/cryst11020167

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