1. Introduction
The
Z1 and
Z2 defect centers were detected for the first time by Hemmingsson et al. in
n-type 4H-SiC epitaxial layers, grown through chemical vapor deposition (CVD), with the net donor concentration in the range of 2 × 10
14–5 × 10
15 cm
−3 [
1]. It is worth stressing that due to the small difference between the thermal electron emission rates, the time constants separation of the two exponential components associated with the electron emission from the
Z1 and
Z2 defects occurring in the capacitance relaxation waveforms measured in the range of 260–340 K was not possible through the correlation procedure applied in conventional deep level transient spectroscopy (DLTS). As a result of the strong overlapping of two DLTS signals, only one peak in the spectrum was observed, and to mark that the two kinds of defects contribute to the thermal electron emission, this peak was named
Z1/2. However, it was possible to resolve the two exponential components induced by the thermal emission of electrons from the
Z1 and
Z2 centers by direct fitting the capacitance relaxation waveforms with the sum of exponential functions [
1]. In this way, the thermal activation energies for electron emission from the
Z1 and
Z2 centers were found to be 0.72 and 0.76 eV, respectively [
1].
Additional measurements performed by Hemmingsson et al. indicated that when the sample was illuminated with 2.64 eV photons before each short (50 ns) voltage pulse filling the defect centers in the space charge region, the two new peaks, associated with the two shallower donor levels
Z10/+ and
Z20/+, appeared in the DLTS spectrum, while merely the trace signal corresponding to the
Z1/2 peak was recorded [
1]. The activation energies for the thermal electron emission leading to the
Z10/+ and
Z20/+ charge state changes were found to be 0.52 and 0.45 eV, respectively [
1]. On the grounds of the DLTS results, a model attributing negative-
U properties to the
Z1 and
Z2 defects was proposed [
1]. It is worth adding that a defect has the negative-
U properties if it can capture two electrons from the conduction band and the second is more strongly bound than the first [
1,
2,
3]. A possible explanation for this phenomenon is that the energy gained through electron pairing combined with a lattice relaxation and/or a change in the atomic configuration in the defect’s vicinity might overcome the Coulombic repulsion at the site. In this way, the effective net attraction needed to give a defect negative-
U properties would arise [
1,
2,
3]. According to the model, each of the
Z1 and
Z2 defects can capture two electrons, and the thermal emission of these electrons to the conduction band is reflected in the DLTS spectrum through the
Z1/2−/+ peak [
1,
2]. This model also assumes that outside the area of the space charge layer, the
Z1 and
Z2 defects are singly positively ionized donors, and after capturing one electron, they become neutral [
1]. As a result of capturing two electrons, they behave as negatively ionized acceptors. The illumination with 2.64 eV photons results in removing one electron, enabling, in this way, in the DLTS experiment, the transitions
Z10/+ and
Z20/+ to be observed [
1,
2]. It has been suggested that the
Z1 and
Z2 centers originate from the same point defect located at the hexagonal (
h) and quasi-cubic (
k) lattice sites, without specifying the defect type and its charge state [
1].
On the grounds of the other model, based on the Laplace-transform DLTS (LDLTS) results combined with the calculations made using the density functional theory (DFT), the
Z1 and
Z2 centers are believed to be acceptors attributed to the carbon vacancy (
VC) located at the
h and
k sites of the 4
H-SiC lattice, respectively [
4]. Before filling with electrons in the LDLTS experiment, the
VC(
h) and
VC(
k) acceptors are neutral, and by capturing one or two electrons, they become singly negatively ionized,
VC(
h)
− and
VC(
k)
−, or doubly negatively ionized,
VC(
h)
2− and
VC(
k)
2−, respectively. The energy levels of the
Z1−/0 and
Z2−/0 transitions involving the one-electron thermal emission, identified with the
VC(
h)
−/0 and
VC(
k)
−/0 charge state changes, are located at 0.48 and 0.41 eV below the conduction band minimum (CBM), respectively [
4]. The deeper energy levels of the
Z12−/0 and
Z22−/0 transitions, involving the two-electron thermal emission and identified with the
VC(
h)
2−/0 and
VC(
k)
2−/0 charge state changes, are located at 0.59 and 0.67 eV, respectively [
4].
For over 25 years, the
Z1 and
Z2 defects have been observed in
n-type epitaxial 4
H-SiC through conventional DLTS as a single electron trap, labeled
Z1/2, with the activation energies 0.63–0.68 eV and concentrations in the range of 10
12–10
13 cm
−3 [
5,
6,
7,
8,
9]. As a result of extensive investigations, several experimental facts have been established. Firstly, by using low-energy electron (100–200 keV) irradiation, it has been found that the
Z1/2 trap can be created through primary displacements of carbon atoms [
8,
9]. Secondly, it has been shown that in the 4
H-SiC epilayers grown under C-rich conditions, the
Z1/2 trap concentration is clearly lower [
5]. Thirdly, the
Z1/2 center has been found to be a thermally stable defect, as no changes in its concentration during annealing up to 1600 °C have been observed [
7]. Furthermore, during the CVD growth at temperatures ranging from 1550 to 1750 °C and the C/Si ratio of 1.5, the
Z1/2 center concentration has been found to increase from ~5 × 10
11 to ~1 × 10
13 cm
−3, respectively, indicating that the interstitials’ involvement in its formation is unlikely [
7,
9]. Fourthly, the possibility of removing the
Z1/2 center by implanting C
+ ions followed by annealing has been found [
6]. It should be stressed that the most important results aiming at finding the atomic configuration of the
Z1/2 center were obtained by Son et al. [
10], who used electron paramagnetic resonance (EPR) to determine the energy levels of
VC and revealed the defect’s negative-
U properties. By performing combined studies using EPR and DLTS, it was found that the most common defects occurring in epitaxial 4
H-SiC, namely, the
Z1/2 and
EH7 centers, are related to the double acceptor level of
VC, corresponding to the (2−/0) charge state change, and to the single donor level of
VC, corresponding to the (0/+) transition, respectively [
10]. However, all of these experimental results do not allow for attributing unambiguously the
Z1/2 center to
VC or presenting an alternative atomic configuration for the center [
5,
6,
7,
8,
9,
10].
To obtain deeper insight into the
Z1 and
Z2 defects properties, we investigated the thermal electron emission from these defects in high-purity semi-insulating (HPSI) 4
H-SiC wafers after filling them with the excess charge carriers generated by 3.31 eV photons. These unique studies were conducted using the Laplace-transform photoinduced transient spectroscopy (LPITS), which enabled the four deep levels for both defects to be resolved. In view of the LDLTS results, the
Z1 and
Z2 electron traps represent the defects located at the
h- and
k-lattice sites, which means the activation energies for electron emission are slightly different and can be distinguished only through high-resolution spectroscopic techniques [
4]. Apart from the activation energies and capture cross-sections, the concentrations of the
Z1 and
Z2 defects contributing to the one-electron thermal emission as well as the two-electron thermal emission are determined. The defects’ properties and concentrations are compared for the pristine and annealed HPSI 4
H-SiC samples. For comparison, the electronic properties and concentrations of the
Z1 and
Z2 defects are also studied in a nitrogen-doped epitaxial layer of 4
H-SiC with a net donor concentration of ~1.2 × 10
15 cm
−3. The results indicating the possible identification of the
Z1 and
Z2 defects with divacancies
VCVSi involving the neighboring carbon and silicon vacancies located at
h- and
k-sites of the 4
H-SiC lattice are discussed.
2. Materials and Methods
The
Z1 and
Z2 defects were characterized in both conducting 4
H-SiC epitaxial layers and semi-insulating wafers originating from a bulk crystal of HPSI 4
H-SiC. The epitaxial layers were grown in an Epigress VP508 horizontal hot-wall chemical vapor deposition (CVD) reactor (Epigress, Lund, Sweden) using silane and propane as precursors. The layers with a thickness of ~10 µm were deposited in a flow of hydrogen, used as a carrier gas, at a temperature of 1600 °C. The epitaxial growth was performed on standard epi-ready n-type 4
H-SiC wafers with a net donor concentration of ~2 × 10
18 cm
−3, oriented 4° off of the (0001) plane toward the [
11,
12,
13,
14,
15,
16,
17,
18,
19,
20] direction. The substrates of 15 × 15 mm
2 in size were cut out from the wafers with a 4-inch diameter. The samples for the LDLTS measurements were in the form of Schottky diodes with a 1 mm diameter, made by evaporating a 20 nm layer of Ti and a 150 nm layer of Ni on the epitaxial layers’ surface for the Schottky contacts and by evaporating 4 nm of Ti and a 150 nm layer of Ni on the substrate back side for the Ohmic contacts. After the electrodes deposition, the samples were annealed at 500 °C. To determine the net donor concentration (
n =
ND −
NA) at room temperature (RT), the samples were initially characterized through capacitance–voltage (
C-
U) measurements.
For investigations of the
Z1 and
Z2 defects’ properties and concentrations, the capacitance relaxation waveforms induced by the thermal electron emission were measured by means of the state-of-the-art experimental system built using a Boonton 7200 capacitance meter, a home-made compensator of the steady-state capacitance, a Janis nitrogen cryostat, and a LakeShore 340 temperature controller. The measurements were performed at a temperature range of 270–370 K with a temperature increment of 3 K. They were carried out at a constant reverse bias of −5 V and a filling pulse amplitude of 4 V. The filling pulse width was 100 µs, and the filling pulse repetition period was 100 ms. To extract the
Z1 and
Z2 defects’ electronic properties from the temperature-induced changes in the capacitance relaxation waveforms, the two-dimensional (2D) analysis were used to transform the waveforms into the 2D spectra by means of both the correlation procedure and the numerical calculations based on the inverse Laplace transformation algorithm (ILT) implemented in the CONTIN code [
11]. In this way, the temperature dependences of the emission rate (
eT) for the detected defects were found and visualized in the 3D space as the ridgelines of the folds in the surfaces created above the plane defined by the absolute temperature (
T) and emission rate (
eT) axes. The temperature and emission rate values were taken from the ridgelines to draw the Arrhenius plots illustrating the dependences of ln(
T2/
eT) versus 1/
kBT, where
kB denotes the Boltzmann constant. The plots were in the form of straight lines, and the values of the electron emission activation energy
Ea and pre-exponential factor
A in the Arrhenius equation were calculated from the slope and intercept of each line, respectively, by means of linear regression. The electron capture cross-sections for the
Z1 and
Z2 defects were determined as
σn =
A/
γn, where
γn is the electron effective mass-dependent material constant, which, for 4
H-SiC, was calculated to be equal to 2.51 × 10
21 cm
−2 s
−1 K
−2. For the convenience of the results presentation, the folds associated with the thermal electron emission were projected onto the plane defined by the (
T,
eT) axes, and the ridgelines illustrating the temperature dependences of the emission rate for the
Z1 and
Z2 defects were depicted. The defects’ concentrations were derived from the amplitudes of the two exponential components of the capacitance relaxation waveforms using the method typical of the DLTS technique [
4,
5]. However, assuming that the two-electron emission occurs from the
Z1 and
Z2 defects [
1,
4], the concentration values obtained in this way were divided in half.
The LPITS studies of the Z1 and Z2 defects’ properties and concentrations were carried out using chips of 1 × 1 cm2 in size cut out from a HPSI 4H-SiC wafer with a thickness of 500 µm and a diameter of 100 mm supplied by Wolfspeed. The wafer was of “epi-ready” type with both sides perpendicular to the c-axis (<0001> direction). The surfaces of both sides were polished; however, the front surface of the wafer lying in the (0001) plane (Si-face) was especially prepared through mechano-chemical polishing (MCP) so that it could be directly used for epitaxy. The LPITS studies were performed using the pristine chips as well as the chips subjected to annealing at 1400 °C for 3 h in argon ambience at a pressure of 100 mbar. The samples were cooled from 1400 to 700 °C for 0.5 h and for 1 h from 700 °C to room temperature. It is worth adding that during the annealing, the morphology evolution of the (0001) plane surface took place. According to the results of atomic force microscopy (AFM) observation, the regularly stepped terraces were formed. The width, length, and height of adjacent terraces were (1–3) µm, (8–10) µm, and ~1.3 nm, respectively. The surface roughness values (Ra), determined for the regions of 20 × 20 µm2 on the pristine and annealed samples surfaces, were 2.09 ± 0.2 nm and 2.3 ± 0.3 nm, respectively.
For the LPITS measurements, the arrays of co-planar contact pairs with an inter-distance of 0.7 mm were evaporated by using an electron-beam. The contacts, in the shape of the squares of 2.5 × 2.5 mm2 in size, were made using a 20 nm layer of Cr and a 200 nm layer of Au. After the deposition of contacts, the chips were diced into samples with dimensions of 4 × 7 mm2, which were used for the photocurrent transient measurements as well as for the measurements of the dark current and mobility–lifetime product temperature dependences. The photocurrent transients were measured in the range of 300–400 K with a temperature increment of 3 K. The pulses of 3.31 eV photons generating the excess electron–hole pairs were emitted using a semiconductor laser produced by Power Technology. The transients were amplified using a Keithley 428 fast current amplifier, operating as a conductance–voltage converter, and then digitized with a 12-bit amplitude resolution and a 1 µs time resolution. In order to improve the signal to noise ratio, the digitally recorded photocurrent transient resulted from averaging of 500 waveforms. The photon flux, the duration of the excitation pulses, and their repetition period were ~1.9 × 1018 cm−2s−1, 50 ms, and 500 ms, respectively. The voltage applied between the two co-planar contacts separated by 0.7 mm was 20 V.
The properties of
Z1 and
Z2 defects were extracted from the temperature-induced changes in the photocurrent relaxation waveforms observed after switching off the UV pulse. The relaxation waveforms were analyzed under the assumption that the decay time of each waveform is much longer than the charge carrier’s lifetime [
12,
13]. According to the assumed model, each relaxation waveform is described by the sum of the exponential components, whose the temperature-dependent time constants are equal to the reciprocals of the thermal emission rate of charge carriers for the defects contributing to the thermal emission process [
12,
13]. Therefore, the temperature dependences of the electron emission rate for the
Z1 and
Z2 defects, and, derived from these dependences, the defects’ electronic properties, were experimentally found by means of the 2D analysis of the photocurrent relaxation waveforms, performed in a similar way as in the case of the capacitance relaxation waveforms. The concentrations of the
Z1 and
Z2 defects were determined from the amplitudes of the exponential components of the photocurrent relaxation waveforms, assuming that the amplitudes are proportional to the concentrations’ electrons trapped by the defects when the optical excitation pulse was switched off. The detailed description of the method used for deriving the concentrations of defect centers from LPITS measurements was published earlier [
12].
The measurements of dark current as a function of temperature were carried out to compare the activation energy of dark conductivity (
EADC) values for the pristine and annealed HPSI 4
H-SiC samples. Because the resistivity of the samples at room temperature was of the order of 10
11 Ωcm, the measurements were carried out in the broad temperature range of 300–700 K at an applied voltage of 20 V. The
EADC values were determined from the slope of linear parts of the dark current characteristics plotted in coordinates of log(
T 3/2/
I) against 1000/
T, where
T and
I denote the absolute temperature and electrical current, respectively. The temperature dependences of the
µτ product for both kinds of samples were also determined in the range of 300–700 K using the transient photocurrent method (TPM), based on the assumption that the
µτ value at a given temperature is proportional to the amplitude of the photocurrent pulse when the optical pulse generating the excess charge carriers is terminated [
14]. The
µτ measurements were carried out at a voltage of 20 V and a possibly low flux of 3.31 eV photons, being ~9.7 × 10
14 s
−1cm
−2, to minimize the effect of lifetime shortening by a higher concentration of the optically generated charge carriers [
12].
4. Discussion
The positions of the
Z1 and
Z2 defects’ energy levels found in this work in comparison with those determined by Hemmingsson et al. [
1] and Capan et al. [
4] are illustrated in
Figure 6. It is worth stressing that the LPITS results for the T1, T2, T1A, and T2A traps represent the first demonstration of the
Z1 and
Z2 defects’ charge state transitions in HPSI 4
H-SiC. This new finding is of great importance in terms of understanding their atomic configurations and formation mechanisms. The data shown in
Figure 6 significantly enlarge the current status of the knowledge on the defects’ properties and encourage theoreticians to further calculate the defects’ transition levels using the currently available powerful computational tools [
21,
22,
23]. The results clearly reveal not only that the energy levels of the
Z1 and
Z2 defects lie in negative-
U ordering, but also the fact that that the levels are closely spaced. The small differences between the ionization energies of the same kind of point defects, below 100 meV, are typical of 4
H-SiC crystals due to the complex lattice structure in which each substitutional defect (vacancy, antisite, impurity atom) can take up two inequivalent sites, namely hexagonal and quasi-cubic, being in an environment with the slightly different electric field distribution [
21,
22,
23]. For example, the ionization energy for nitrogen atoms, commonly dopants used in 4
H-SiC as shallow donors, are 60 and 110 meV when they are located in the
h and
k sites, respectively [
22].
According to the model presented by Capan et al. [
4], the
Z1 and
Z2 defects in 4
H-SiC can adopt three charge states, 2−, 1−, and 0, similarly to the carbon vacancies. This means that before trapping electrons, the defects are neutral. Such a model was earlier investigated for the oxygen negative-
U center in GaAs [
3], and it was found that the defect’s atomic configurations for the two negatively charged states are very similar; however, those of the neutral charge state clearly differ from the former. From the recent calculations, which employed hybrid density functionals, it has been found that the configuration coordinate diagrams for the negative and doubly negative charge states of
VC, occupying the
h- and
k-sites in the 4
H-SiC lattice are also very similar but substantially differ from that for the neutral
VC. In other words, the capture of one or two electrons vitally changes the arrangement of the four Si atoms in the defect’s vicinity. The interpretation of the charge state changes of
Z1 and
Z2 defects given by Hemmingsson et al. [
1] is unlike that proposed in Ref. [
4]. Although it has also been assumed that each of the defects can be in the three charge states, these states are 1−, 0, and 1+. This means that before the electron capture, the defects are positively charged and behave like typical donors, becoming neutral after the capture of one electron. On the other hand, after capturing two electrons, they behave like acceptors, being negatively charged. Such a model allows us to understand why in the normal DLTS experiment only the two-electron thermal emission is observed, and special modifications in the traps’ filling procedure are needed to avoid the accumulation of the centers in the negative charge state. In other words, generating the capacitance relaxation waveforms corresponding to the
Z10/+ and
Z20/+ transitions is difficult, even in the LDLTS experiment. This is because the Coulombic interaction between the positively charged defects and free electrons enhances the capture of two electrons from the conduction band. In the experiment described in Ref. [
1], to reveal the one-electron emission, the sample was illuminated with 2.64 eV photons before each filling pulse of a 50 ns width. In this way, the first electron was optically released from the negatively charged
Z1 and
Z2 defects, and the second was thermally emitted from the neutral defects. In the experiment described in Ref. [
4], when using LDLTS, the one-electron emission signals were detected by applying filling pulses of a 100 ns width after quickly cooling the sample from room temperature down to 220 K at a reverse bias of −5V.
The negative-
U properties, similar to those proposed for the
Z1 and
Z2 defects in Ref. [
1], were established experimentally for interstitial boron in silicon [
2]. According to the EPR results and the model proposed by Watkins and Troxell, interstitial boron in this material exists in three charge states, 1−, 0, and 1+, giving rise to two energy levels lying in the upper half of the silicon bandgap [
2]. These are the donor level at
EC – 0.13 eV and the acceptor level at
EC− 0.37 eV. It should be added that the former was detected only through the modified DLTS experiment in which an optical trap-filling pulse was used to prevent the B
i atom from getting into the negative charge state [
2]. In our LDLTS experiment aimed at the studies of the
Z1 and
Z2 defects’ properties and concentrations in epitaxial 4
H-SiC, only the electrical trap-filling pulse was used, and the defects were predominantly transferred to the negative charge state during each pulse. Therefore, the measured capacitance relaxation waveforms were induced by the thermal emission of two electrons from the defects’ acceptor levels. In the LPITS experiment, the
Z10/+ and
Z20/+ transitions, involving the thermal one-electron emission from the donor levels, and the
Z1−/+ and
Z2−/+ transitions, involving the thermal two-electron emission from the acceptor levels, are perfectly separated. It is worth noting that all of the transitions with the various emission rates were detected in the same temperature range, 300–400. According to the Laplace spectral fringes shown in
Figure 4, the photocurrent relaxation waveform induced by the thermal electron emission at 340 K was composed of the four exponential components related to the
Z20/+,
Z10/+,
Z1−/+ and
Z2−/+ transitions with the time constants of around 0.11, 0.28, 1.4, and 5.6 ms, respectively. The time constants of the exponential components induced by these transitions and resolved in the waveform recorded at 360 K were 0.044, 0.11, 0.35, and 1.8 ms. The fact that the charge state transitions of the particular defects could be well-determined indicates that the defects’ donor and acceptor levels had been almost fully filled electrons before the thermal emission began. It should be added that the excess electrons trapped by the defects were excited from the valence band by illuminating the samples with 3.31 eV photons when the optical excitation pulse was turned on. During the illumination, two processes took place. Firstly, two electrons were captured by the positively charged
Z1 and
Z2 defects, placing the defects in the negative charge states, and then the second electron was by the optical emission transferred back to the conduction band, leaving the defects in the neutral charge states [
1,
2]. When the optical trap-filling pulse was turned off, the thermal emission of single electrons began, and the electrons released from these states were recaptured by the positively charged defects, which again appeared. This process is possible because the electron capture rate is much higher than the rate of either one-electron or two-electron thermal emission, which is strongly temperature-dependent [
1,
2].
The question of whether the thermal one-electron emission follows the electron capture by the donor or acceptor levels of the
Z1 and
Z2 defects is still unsolved [
1,
4,
10]. However, the results of this work shown in
Figure 6 are more in line with the model of Hemmingsson et al. [
1] than with that of Capan et al. [
4]. The values of the T2A and T1A traps’ activation energy are very close to that of the former and clearly higher than that of the latter. The diminished values of the activation energy reported in Ref. [
4] may result from the effect of the electric field on the emission rate, known as the Poole–Frenkel effect [
2,
24]. During the measurements described in Ref. [
4], the thermal emission took place in the depletion region of the Schottky diodes, where the electric field was around 5 × 10
4 V/cm. On the one hand, the photocurrent transients generated in the HPSI 4
H-SiC samples were measured at the electric field of ~3 × 10
2 V/cm. According to the Poole–Frenkel effect, the electric field enhances the emission rate in the case of charge carrier emission from a Coulombic attractive center due to lowering the potential barrier by an amount proportional to the square root of the electric field. It should be noted that the Arrhenius plots shown in Ref. [
4] for the one-electron emission from the
Z1 and
Z2 defects are shifted towards lower temperatures by ~50 K with respect to those plotted in
Figure 5, and their slopes are accordingly smaller. These facts unambiguously indicate that the lower activation energy values obtained by Capan et al. [
4] compared to those obtained in this work for the T1A and T2A traps are due to the Poole–Frenkel effect [
2,
24]. Thus, the assumption made in Ref. [
4] that the one-electron emission occurs from the acceptors levels of
Z1 and
Z2 defects is in contradiction with the LPITS results and requires further experimental verification. It is worth stressing that the activation energy values for the T1 and T2 traps related to the two-electron emission from the defects’ acceptor levels are in good agreement with those determined in Ref. [
4]. This result confirms the theoretical predictions that in the case of thermal electron emission from acceptor levels, the Poole–Frenkel effect is inactive [
2].
The experimental findings of this work show that the [
Z2]/[
Z1] values determined for the epitaxial 4
H-SiC as well as for the pristine and annealed semi-insulating materials are 4.38, 0.57, and 1.23, respectively. Simultaneously, the [
Z1−]/[
Z10] values after filling the defects with electrons for the pristine and annealed HPSI 4
H-SiC are 1.11 and 1.29, respectively, and the [
Z2−]/[
Z20] values are, accordingly, 1.25 and 2.7. The view that
Z1 and
Z2 defects are associated with specific point defects residing at the
h and
k sites has been, so far, widely assumed [
1,
4,
10]. However, the question arises as to whether these defects can be identified with an individual point defect, like a carbon vacancy being in a certain charge state, e.g., (
VC−(
h),
VC−(
k)), occupying the different lattice sites [
4,
9,
10], or a complex defect, such as the divacancy, which may be neutral, negatively charged, or positively charged and consists of the nearest-neighbor carbon and silicon vacancies located at the mixed lattice sites, e.g., (
VC(
h)
VSi(
k))
− [
15,
25,
26]. It should be added that in the 4
H-SiC lattice, the
VCVSi complex can take on four possible configurations, namely
hh,
kk,
hk, and
kh, and in each of them it can be in various charge states [
23]. Moreover, the
VCVSi axis can be parallel or inclined at the angle of ~109.5° to the
c-axis and, as a result, the defect can exhibit either the high, axial (C
3v) or the low, orthorhombic (C
1h), symmetry [
23]. The defects with the former and latter symmetries are referred to as the axial and basal plane divacancies, respectively [
23]. The primitive cell of 4
H-SiC with the divacancies formed by the
VC and
VSi occupying the same and mixed lattice sites is schematically illustrated in
Figure 7.
The discrepancy between the [
Z2]/[
Z1] values determined for the three kinds of 4
H-SiC crystals subjected to various technological treatments can be discussed by taking into account how far from the thermodynamical equilibrium these materials were during these processes [
17]. Under the equilibrium conditions, the
h and
k sites are likely to be equally occupied by the isolated
VC and
VSi. This postulate is theoretically supported by the close values of the formation energy (
Ef) for neutral
VC(
h) and
VC(
k). For stoichiometric 4
H-SiC, these values are 4.21 and 4.07 eV, respectively, which means that the concentration ratio of [
VC(
k)]/[
VC(
h)] is 1.15 [
15]. The
Ef values for the neutral
VSi(
h) and
VSi(
k) in stoichiometric 4
H-SiC are 8.26 and 8.37 eV, respectively, which gives the [
VSi(
k)]/[
VSi(
h)] ratio of 0.9 [
15]. Furthermore, the calculated energy barriers for
VC migration in 4
H-SiC along axial (
h-
h,
k-
k) and basal (
h-
k-
h-
k) paths are approximately the same [
16,
27]. It has been stated that even ion implantation, which is an extremely non-equilibrium process, generates approximately equal concentrations of
VC(
h) and
VC(
k), irrespective of the defect distribution between the sites prior to implantation [
16]. Thus, the
VC(
k)/
VC(
h) ratio considerably higher than 1 for the epitaxial material seems to be unlikely [
16,
23,
27]. In other words, assuming that the
Z1 and
Z2 defects originate from the two point defects involving different atomic configurations is fully justified.
It is worth noting that in the pristine and annealed semi-insulating materials, the concentrations of the
Z1 defect in the two charge states (0) and (−), arising after the capture of one or two electrons, respectively, are approximately equal, which indicates that all of the
Z1 defects present in these materials are taking part in the charge state transitions. However, in the case of the
Z2 defect, the [
Z2−]/[
Z20] ratio is close to 1 only for the pristine material, being equal to 2.7 for the annealed material. This result fully corresponds to the data shown in
Figure 5, illustrating the temperature dependences of the excess charge carriers’ lifetime. According to the Laplace spectral fringes (
Figure 3), the
Z2+/− transitions occur at temperatures close to 380 K, at which the carriers’ lifetime in the annealed material is ~7 × 10
−9 s, while the
Z2+/0 transitions are observed at temperatures near 330 K, when the lifetime is by a factor of ~2.5 shorter. This phenomenon results in the accordingly lower excess electron concentration to be captured by the
Z2+ defects [
20].
The most important result of this work is the demonstration that as a result of the heat treatment at 1400 °C, the
Z1 and
Z2 defects’ concentrations in the semi-insulating material go up from 2.1 × 10
13 to 2.2 × 10
14 cm
−3 and from 1.2 × 10
13 to 2.7 × 10
14 cm
−3. This finding is in line with the postulate that the
Z1 and
Z2 defects can be attributed to the
VCVSi complexes consisting of the nearest-neighbor
VC and
VSi residing at
k and
h as well as
h and
k sites, respectively. This postulate is based on the results of first-principles calculations for divacancy defects in 4
H-SiC, revealing their formation energies and stability, their ionization levels, and symmetry point groups corresponding to neutral and charged states [
15,
23]. According to these results, the divacancies possess a remarkably high binding energy of ~4 eV, and there is a strong dependence of the formation energies on the sites in which the nearest-neighbor silicon and carbon monovacancies are located [
15]. It is worth stressing that these results, derived from the calculations based on the density functional theory (DFT), prove that a negative-
U behavior leading to the 1+/1− charge state transitions occurs only for the divacancies involving the
VC and
VSi located in the nearest neighborhood at the
hk and
kh lattice sites, respectively [
15]. Furthermore, it has been shown that the formation energy of the
VC(
h)-
VSi(
k) complex is ~0.3 eV lower than that of the
VC(
k)
VSi(
h) one [
15]. By assuming that the
Z2 defect is attributed to the
VC(
h)
VSi(
k), it is possible to understand why this defect is predominant in the
n-type epitaxial 4
H-SiC used for the studies described in this work as well as in Ref. [
4]. Moreover, the impact of the C/Si ratio during the epitaxial growth on the
Z2 defect’s concentration can also be explained. For a long time, this concentration has been hidden in that of the artificial
Z1/2 defect [
5,
7]. At the low C/Si values (0.4–1), the [
VC] dominates, and the
Z2 defect concentration is limited by the [
VSi]. On the other hand, at the most frequently used C/Si values in the middle range of (1.5–3), the [
VSi] increases, and the [
VC] contributing to the
Z2 defect formation may be additionally increased by the irradiation with low-energy particles [
7,
8,
9]. At a C/Si value of six, the [
Z2] was found to drop by a factor of four compared to that at the C/Si = 3 [
5]. In this case, the [
VC] significantly decreases, and
VSi is likely to be transformed into the carbon-related antisite–vacancy (CAV) pair, being an isomer of the
VSi in the form of the C
SiVC complex [
21].
Attributing the
Z1 and
Z2 defects to the
VC(
k)
VSi(
h) and
VC(
h)
VSi(
k) divacancies, respectively, allows us to understand two facts established in this work. The first says that as a result of the heat treatment, the defects’ concentrations in the HPSI 4
H-SiC increase by an order of magnitude. During annealing, both
VC and
VSi become mobile, and, to be bound, their migration paths should exceed the distance between them. Assuming that the activation energy for
VC migration is 3.6 eV and the pre-factor
D0 = 0.54 cm
2/s [
16,
27], the
VC diffusion coefficient at 1400 °C (
D) equals 7.8 × 10
−12 cm
2/s, and the diffusion length
LD = 2(
Dt)
1/2 for
t = 3 h is ~6 µm. According to the reported data [
19], the [
VC] and [
VSi] concentrations in HPSI 4
H-SiC crystals are ~1 × 10
15 cm
−3, which means that the average distance between the vacancies is ~0.08 µm. Thus, the interactions leading to the divacancy formation are very likely [
24,
26]. The second phenomenon is that the
Z1 and
Z2 defects’ concentrations after the heat treatment become approximately equal. This means that due to the mutual interactions between the isolated vacancies, the [
VC] and [
VSi] tend to equalize, and the
h and
k lattice sites become equally occupied by the vacancies [
23].
It is worth stressing that there are interesting results obtained through photoluminescence (PL), EPR, and positron annihilation spectroscopy (PAS) showing the unique properties of divacancies in HPSI 4
H-SiC [
25,
26,
28,
29,
30,
31]. The neutral (
VCVSi)
0 divacancy optical properties have been widely investigated through PL measurements, and the photon energies corresponding to the zero-phonon lines (ZPLs), labeled PL1, PL2, PL3, and PL4, have been precisely determined [
28,
29]. The energies indicated by the distinct PL peaks’ positions for these lines are 1.095, 1.096, 1.119, and 1.150 eV, respectively [
29]. By correlating these energies with theoretical values derived from hybrid density functional theory, equal to 1.056, 1.044, 1.081, and 1.103 eV, respectively, the PL1 and PL2 lines have been attributed to the divacancy in the axial
hh and
kk configurations, respectively, while the PL3 and PL4 lines have been assigned to the divacancy in the basal
kh and
hk configurations, respectively [
28,
29]. Furthermore, the PL1–PL4 lines have also been attributed to the P6/P7 centers (S = 1) detected in the EPR spectra measured for as-grown and electron-irradiated HPSI 4
H-SiC [
25,
26]. The spectra were measured at 77 K using the magnetic field parallel to the
c axis and the sample illumination with photons, whose energy was in the 2.0–2.8 eV range. Based on the first principles calculations results, which enabled the parameters of hyperfine (HF) interactions to be found, as well as the results of the angular dependence and intensity of the HF lines’ measurements, the P6/P7 centers were identified with the (
VCVSi)
0 divacancy in the four configurations, including the
hh and
kk ones, with the
C3v symmetry, as well as the
kh and
hk ones, with the
C1h symmetry [
25,
26]. The divacancies in the former two configurations were observed in the EPR spectra as the P6b and P6′b centers, respectively [
25,
26]. The defects in the latter two configurations are known as the P7′b and P7b centers, respectively [
25,
26]. Thus, the PL1, PL2, PL3, and PL4 zero-phonon lines attributed to the (
VCVSi)
0 in the
hh,
kk,
kh, and
hk configurations are also assigned to the P6b, P6′b, P7′b, and P7b centers observed in the EPR spectra, respectively [
25,
26,
29]. The positron annihilation studies of vacancy-type defects in SI 4
H-SiC have shown that the material contains vacancy clusters, and the positron trapping to the clusters is enhanced by annealing for 1 h at 1600 °C in H
2 ambient [
30]. Interesting results obtained through the positron lifetime measurements are presented in Ref. [
31]. It has been shown that after annealing the SI 4
H-SiC samples for 1 h at 1600 °C in H
2 ambient, the material resistivity drops from 2 × 10
9 to 1 × 10
8 Ωcm, but the concentration of
VSi-related defects, being in the pristine material ~6 × 10
16 cm
−3, increases by nearly 10%, while the concentration of Si vacancy clusters decreases from 8 × 10
15 to 6 × 10
15 cm
−3 and the number of vacancies in the cluster rises from 5 to 16 [
31]. This fact indicates that silicon vacancies are mobile during the annealing and they diffuse either outwards of the cluster or inside the cluster. The mobile Si vacancies can also form more stable complexes, such as
VCCSi or divacancies, which act as compensating centers after the thermal treatment [
31].