1. Introduction
Radiation defects in diamond and silicon are examined actively because these semiconductors are widely used as detectors and some other devices in high energy physics. The main problem of these investigations is connected with their radiation hardness (see e.g., some recent works [
1,
2,
3,
4,
5,
6]). Radiation defects induced by protons, neutrons and heavier particles with kinetic energies
MeV are studied in most works. After slowing down to kinetic energies
E less than ionization energies of host atoms in crystals, these impinging particles stop in the lattice, damaging it. If the recoil energy is higher than the lattice binding energy, a host atom will be displaced from its site. Numerical modeling of these processes is carried out, e.g., in [
3,
5]. In [
7], many types of these radiation defects in diamond are well described. The second problem is implantation of ions in crystals for preparing necessary impurity atoms. Ion implantation is a commonly-used method for modifying properties of materials in the field of microelectronics. The application for diamond is represented, e.g., in [
8].
For a number of reasons, radiation defects induced by light negatively-charged particles used to be out of interest in high energy physics experiments. First, these particles are secondary particles, as a rule, and, second, they cannot inflict many lattice damages compared to protons. Nevertheless, the interaction of these particles with crystals can be very important for many different applications of electronic devices.
We will consider in this article radiation defects induced by light negatively-charged particles like pions (
-mesons) and muons in crystals with a diamond structure. These particles do not destroy the crystal structure, like heavy particles, but can create specific defects in the lattice. Indeed, they are captured by a nucleus creating an impurity atom and, thus, can change the electronic properties of a crystal. Negative pions and muons are respectively long-lived particles: the lifetimes of charged pion and muon are
s and
s, respectively. Pions are born as usual when high energy protons are stopped in a target, and muons are born after the decay of pions:
where
is a muon antineutrino for the negative muon and a neutrino for the positive one. This picture can be observed in cosmic rays. Negatively-charged pions and muons are stopped in a media very effectively because of the capture by nuclei.
The capture mechanism differs for pions and muons, but the result manifests in the same way in electronic properties, because they form finally the same acceptor impurity. Consider a capture of negative pions by stable nuclei of the main semiconductors: C, Si and Ge. In diamond, we have only one stable isotope C
, and a capture of a negative pion gives rise to the boron acceptor:
where
n is a decay neutron. The boron nucleus spin is
.
In silicon isotope Si
(92% in nature, see e.g., [
9]) can capture
and than transform into the aluminum acceptor:
Processes in germanium are more complicated, because it has only two stable isotopes with an atomic number equal to 70 and 72 with 21.2% and 22% in nature, respectively [
9], which can capture a negative pion and decay to an appropriate gallium isotope. Therefore, we have:
Both gallium isotopes possess a spin equal to 3/2. All isotopes B, Al, Ga and Ga are stable, and the capture of the negative pions irreversibly changes the concentration of the main acceptor impurities in semiconductors.
A capture process of negatively-charged muon in crystals strongly differs from the negative pion capture process. Consider this difference in more detail. Positively- and negatively-charged muons (
and
) are widely used for research of condensed matter in many different areas, for the simulation of the behavior of hydrogen-like light element impurities and chemical processes with atomic hydrogen (see e.g., [
10]). The application of muons for materials’ investigation has become possible due to a well-developed
SR-technique based on the possibilities of supervision for a muon magnetic moment in the sample. Negatively- and positively-charged muons (
) are unstable leptons with spin 1/2.
The negatively-charged muon (
) decays according to the scheme:
where
and
are muonic neutrino and electronic antineutrino, respectively. The escape probability of a decay electron depends on the angle between the electron momentum direction and the average muon spin
s, due to what appears to be the possibility to study local fields of a target. A muon has a relatively high decay time of
s. The large lifetime allows investigating with a high precision the processes with a characteristic time
s, which provides the opportunity for a
SR-technique for the material property studies, well comparable with the possibilities of the widely-applied methods of NMR and ESR.
The behavior of and in a medium is radically different. From the chemical point of view, the positively-charged muon is a light element impurity modeling a light hydrogen isotope. The negatively-charged muon cascades into the ground -state forming a muonic atom (-atom). The mass of a muon equals 207-times the mass of an electron, and therefore, its binding energy with an atomic nucleus is 207-times larger than that of the electron. After a muon capture, much energy is released, leading to a high ionization of a target atom due to the emission of Auger electrons. Further, the target Auger electrons are captured by the positively-charged radiation-induced defect. Due to a high muon mass value, the negative muon screens a nuclear charge Z, which is effectively becoming . After defect neutralization, a replacement impurity is formed, or a muonic atom, similar to an atom isotope with a nuclear charge .
This fact was well known since the initial stage of the muon research (see e.g., [
11]) and gave rise to the foundation of the muon method of materials research (
SR). Systematic study of impurities’ formation with a nuclear charge equal to
in condensed matter was carried out at the early stages of the
SR-research [
12,
13,
14]. The muonic atom formed inside a semiconductor lattice models an acceptor center. For example, in diamond (
), the negative muon, as a result of capture by a nucleus, forms a pseudo-boron, or muonic boron, which can be designated as
B. In Si (
) and Ge (
), the negative muon is captured by a nucleus forming the pseudo-aluminum
Al with a nuclear charge equal to
and pseudo-gallium
Ga with a nuclear charge equal to
, respectively. These chemical elements are the main acceptor impurities in silicon and germanium semiconductors. Therefore, a radiation defect induced by a negative muon is unstable and disappears after its decay. However, nevertheless, this kind of defect is very interesting because it provides the possibility to study the evolution in time of the processes considered above.
The study of acceptor center properties using
was suggested in [
15]. The possibility to extract valuable information about the hyperfine structure and interactions with a lattice of acceptor centers in different semiconductors with the help of negative muons was shown in the works [
16,
17,
18,
19]. Recently,
research of synthetic diamond crystals were carried out in [
20,
21,
22,
23]. We will examine evolution in time of radiation defects induced by negative muons in the following, taking in mind that they are the same for negative pions, as well.
The total process of this kind of radiation defect formation can be separated into two principally different stages. At the first stage, a center with a large positive charge appears after the stopping of a negative muon or pion. This center interacts with electrons of a trace, created by the charged particle, when it is decelerated in a crystal. As a result of this interaction, the positively-charged center partially compensates its charge or becomes neutralized if it is possible. At the second stage, the center with a compensated charge restores broken chemical bonds with a lattice and then forms an acceptor center.
Now, we will outline briefly the main results of the first stage and show the difference between diamond and other diamond structure crystals following [
24]. After that, we will consider the second stage in more detail.
2. Interaction with the Trace
When a negative pion is captured by a nucleus or a negative muon is captured to the
K-shell of the muonic atom, substantial energy (
keV) is released. A totally ionized positively-charged center and Auger electron ionization environment appear, creating secondary electrons. This process takes short a time,
s. After the ionization, free electrons lose their energy; for a while, it will be of the order of the forbidden gap energy. This is a diffusion process, when ionized impurities and host atoms are neutralized, and it takes respectively a long time
s. All of these findings were obtained as a result of a numerical modeling of a neutralization process of muonic atoms in a kinetic approximation for diamond and silicon crystals with different concentrations of impurities [
24].
A different situation is observed in diamond and silicon already at this stage. The capture of a negative muon on a silicon nucleus creates the number of free charge carriers approximately at two orders more with respect to diamond. This result is connected with a difference in the number of Auger electrons, in ionization energies of impurities and a forbidden energy band for these two crystals.
Numerical calculations have shown that the recombination frequency of electrons with positively-charged ions in diamond reaches approximately s only at a respectively short interval of time s. In silicon, the recombination frequency reaches approximately s for the same interval of time. As a result, all ionized impurities in silicon including a muonic atom Al are neutralized for a very short time s. The probability of the neutralization of a muonic atom B is less than for the interval s. The recombination frequency in diamond sharply falls for s, and the neutralization time in this process becomes more than both a muon lifetime and a characteristic time for chemical bonds’ restoration.
Thus numerical modeling has shown that a positively-charged radiation defect, created by a negative muon in silicon and germanium, must be quickly neutralized before chemical bonds with the lattice can be restored. In diamond, we observe other behavior. Namely, the radiation defect must restore chemical bonds with the lattice to be positively charged. Therefore, we need to consider different initial states of the radiation defect at the second stage for diamond and other crystals with the diamond structure.
The second stage of the radiation defect formation consists of a few steps that lead to an acceptor center formation. A neutral defect with restored chemical bonds is not an acceptor center yet, because there exist unsaturated chemical bonds. Therefore, we need to consider at least three steps of an acceptor center formation:
- (1)
restoration of broken chemical bonds and neutralization of the radiation defect;
- (2)
capture of a missing electron and saturation of chemical bonds of a neutral radiation defect in the lattice; formation of an ionized acceptor center;
- (3)
formation of an acceptor center in the ground state.
The first two steps are discussed in this article.
3. Electron States of a Neutral Radiation Defect in Si and Ge
The muonic impurity atom is in an exited state just after formation because its chemical bonds with host atoms are broken. In accordance with the standard idea of quantum chemistry, only electrons with the same principal quantum number can create a chemical bond if they were on an unfilled energy level of the atom. In this case, they form hybridized states. For lattices with diamond structure electron states, and are represented with equal probability, where and 4 relate to C, Si and Ge, respectively. Hybridized states are formed in atomic time, but chemical bonds’ formation is determined by exchange interactions that are weaker than Coulomb interactions, which form the appropriate atomic configuration.
When a chemical bond is formed, a significant energy (of the order of some eV) can be emitted. In gasses and liquids, this excess energy can be transferred to the third body. This kind of energy transfer in crystal must be connected with a phonon emission. One-phonon emission with the energy of ≥1 eV in covalent crystal is impossible. Therefore, a transfer of this energy value could be realized in the case of a multi-phonon process. This kind of process have a very small probability. Thus, a radiation transition with a photon emission seems to us more preferable with respect of the other processes.
In this section, we consider this process for a neutralized radiation defect in Si and Ge when three electrons are in hybridized states [
25]. Consider an impurity atom with the nuclear charge
, which is formed as a result of a neutralization process in an atomic time and has an atomic configuration where electrons at the external shell are in the “mixed”, but not in the ground, state:
We assume that the electron configuration with the principal quantum number less then
n is completely occupied, and the state of such electrons is described by the unperturbed wave function of the free atom. The initial state of the radiation defect in a silicon lattice is sketched in
Figure 1a.
The mixed state (
1) forms a chemical bond, if it possesses by maximal spin
, and a spacial part of its wave function one may represent in the view of equiprobable superposition of three Slater’s determinants:
where:
Here, , and .
A spacial part of the wave function of the defect in a final state can be represented as a superposition of three hybridized states forming the chemical bond with host atoms of the lattice:
Here, summation is carried out over all permutations
of the valence electrons of the impurity and:
The unit vectors
are directed from the impurity to nearest neighbors (along the directions of the chemical bonds). The one-particle functions
are the hybridized states with directed bonds, and they can be written in a form (see e.g., [
26,
27]):
and satisfy normalization conditions:
From this condition, we can obtain relations for the coefficients in the superposition (
6):
4. Electron States of a Positively-Charged Radiation Defect in Diamond
A positively-charged radiation defect in diamond has the effective nuclear charge
. Its atomic configuration contains only two electrons in an external (unfilled) electron shell, and they are in the “mixed” state [
28]:
We suppose also that the electronic configuration for the principle quantum number
is completely filled, and electronic states of the external atomic shell are described by unperturbed wave functions of a free atom. The initial state of the radiation defect in a diamond lattice is shown schematically in
Figure 2a.
The mixed state (
9) forms a chemical bond with the nearest host atoms of the lattice. We express a wave function of the initial state (
9) in distinguish from the state (
1) in the form of superposition with all possible spin states:
where:
is the wave function of the
-state. We assume that all
p-states with different projections have equal probabilities:
are the values of the total electron spin;
is the appropriate spin-state vector. Spin states with different projections are considered as having equal probabilities; so coefficients in the superposition (
10) satisfy the following condition:
The space part of the defect wave function in the final state must correspond to the determined value of the total electron spin
S, and this can be represented in the form of the superposition of hybridized states providing a chemical bond with the lattice host atoms:
where summation is carried out over all possible directions of the chemical bond of impurity valence electrons with the nearest neighbors of the lattice,
The unit vectors
are directed from the impurity to the neighbor atoms (along the direction of chemical bonds) like in (
5). One-particle wave functions
of hybridized states with directed bonds are determined by Equations (
6)–(
8).
The final state of the charged impurity is described by a function similar to the superposition (
10) where the wave functions with a determined spin must be replaced by Expression (
13).
5. Formation of the Neutral Center (A A in Si and Ge
A lifetime of exited states (
1) and (
10) is determined by a rate of a radiation transition in a bond state and can be calculated by using Fermi’s “golden rule”:
The interaction operator is:
where
is the vector-potential of the free radiation field.
Let us consider now only the term for one electron with
in the operator (
14) to simplify the following calculations. In this case, matrix elements of the perturbation operator could be represented by the expression:
and:
where the indexes are
.
If we direct the axis
, then the other three
p-states in hybridized states (
6) turn out as a result of the rotation of the state
in the state with a rotation moment projection equal to zero on the axes
. In this case, we get the opportunity to calculate easy integrals incoming in Expressions (
16) and (
17):
where
is the rotation operator, and the matrix elements (19) are determined by the second column of the rotation matrix:
We put in Equation (
18) the obvious expressions of the matrix elements of the rotation operator (see e.g., [
29]).
Without reduction of the generality of the calculations, we consider a matrix element only for the
z- projection of a momentum operator. Therefore, we have:
where:
For clarity, we show some intermediate calculations. The total matrix element in Expression (
14) consists of 72 different items corresponding to different matrix elements between states of the superpositions (
2) and (
4). However, it is enough to calculate only four of them. We give them below.
The rest of the three matrix elements are determined by the other possible sets of
for the electrons with coordinates
and
:
It is easy to see that permutations of electrons in the superposition (
5) do not change expressions for the matrix elements (
24)–(27). Therefore, the number of permutations in the state (
5) with similar expressions reduces the total number of items by six-times.
We get a very cumbersome expression for the arbitrary values of the parameters
and
. However, it is necessary to take into account that the system under consideration has a symmetry at less
. In this case, the result could be essentially simplified. We examine the simplest case at first, when the system has a tetrahedral symmetry and all parameters in the hybridized states (
6) are equal to each other:
If the vector
lies in the
plane, the angles
are equal:
In the case of symmetrical structure, we get the following expressions for the matrix elements:
Adding up Expressions (
30)–(32), we get:
The matrix elements for the state with
are calculated by a similar way:
For the state with
:
After substitution of the values o Parameters (
28) and (
29) and adding up Expressions (
33)–(
35), we have:
For the calculation of the integral
, we take the appropriate wave functions of the hydrogen-like atom with an effective nuclear charge equal to
. In accordance with Slater [
30], an effective charge is determined as
, where
Z is the real nuclear charge and
is a screening constant.
In a silicon crystal, an aluminum
-atom
Al is formed. It has the principle quantum number
, and appropriate calculations for
Al give the following results:
where
is the Bohr radius.
In a germanium lattice, a gallium
-atom
Ga with the principle quantum number
must be formed. The unknown value of the matrix element for
Ga is equal to:
In the calculation of a transition probability per unit time, we will take into account that at least three electrons participate in the matrix element of the operator (19), and the number of spin states in Determinant (
3) is
:
Here, appropriates the transition frequency of a neutral radiation defect from the energy level of the corresponding free atom state on the energy level corresponding to a hybridized state in a lattice.
After integration over the wavevector of photons and averaging over all angles, we get:
where
is the fine structure constant. We have substituted here the effective charge
and
in accordance with Slater [
30].
6. Formation of Neutral Center (B C in Diamond
Formulae (
16) and (
17) must be modified for diamond in accordance with Equations (
11) and (
13) as follows:
and:
To derive a common expression now, we consider that the
Z-axis does not coincide with any of the bond directions
. Therefore, as for (
18) and (19), we have:
Without the reduction of the generality of the calculations, we consider a matrix element only for the
z- projection of a momentum operator.
where
is determined by Equation (
23) for
.
The interaction operator (
15) conserves the total spin, and we need to calculate only matrix elements for the superpositions (
10) and (
12) between states with equal total spins. For states with a total electron spin
, the matrix elements of the interaction operator in Expression (
14) are equal to zero because of the symmetry of the two-particle states (
11) and (
13). Therefore, contrary to calculations performed for neutral defects in Si and Ge, it is necessary to calculate matrix elements only for singlet states. Accordingly, a radiation transition in dipole approximation for triplet-states (with a maximum spin value of the (
B)
defect) is forbidden. Note that singlet spin states constitute only 1/4 part of all spin states under the assumption of an equal probability of a population of all spin states. This fact leads to the reduction of a total probability transition by nine-times with respect to the probability transition between triplet states.
Taking into account a local symmetry
of the
B C
cluster, we direct axes as shown in
Figure 2b and introduce the following designations:
. We have, respectively,
. After summation over all bonds, one obtains the expression:
We can determine superposition coefficients in (
6). Taking into account the relations (
8) as
and, respectively,
. Therefore, we have:
Here, the condition and, respectively, were used.
Finally, the matrix element (
47) can be written in the form:
For the appropriate wave functions of the hydrogen-like atom with an effective nuclear charge equal to
, the integral
is equal to:
An effective charge for a boron atom is equal to
[
30].
In the calculation of the probability of a transition per unit time, we shall take into account that two electrons participate in the matrix element of the operator (
15). Carrying out integration over a wavevector of photons and averaging over all angles, we obtain:
Here, corresponds to a transition frequency of the charged radiation defect from the energy level of the free ion (B) to the energy level of the hybridized (bound) state in the lattice B C.
In accordance with the matrix element (
49) and the result of the matrix element (
50) calculation, a configuration factor
is equal to:
A transition frequency
and angles for a configuration factor
were calculated numerically in [
28] by the quantum-chemical methods. The crystalline chemical environment of clusters in
Figure 2 has been taken into account by the procedure [
31] based on a passivation of unnatural valences on a border of the cluster by hydrogen atoms. The variation of the geometrical position of the H* atoms, if it is possible, ensures the stoichiometry of the charge distribution on the carbon atoms of the model C
H*
fragment (
Figure 2a, where the B atom is substituted by the central C atom).
The initial state of the radiation (
B)
defect in a diamond lattice is modeled by the tetrahedral structure of
Figure 2a where the central carbon atom is substituted by the B atom, and the non- equilibrium length of four B–C bonds coincides with the equilibrium length of 1.523 Å for C–C bonds found by us earlier after geometry optimization of the central structural C
H*
fragment. As a result of a structural relaxation, the defect transits into the final hybridized state
described by the lowest in energy triplet structure of the [BC
H*
cluster (
Figure 2b). The energy of such a transition is calculated from the difference between the total energies of levels
and
to be equal to 1.17 eV.
The angles in the cluster are equal to
. Calculations of a spin density at the center of the cluster give the value of the superposition parameter
. The effective charge for
and
states of the boron atom is
. Substituting the calculated values into Formulas (
51) and (
52), we obtain a numerical estimate of the radiation transition rate of the impurity center
B
into a hybridized state:
The obtained value confirms the validity of the assumption on the kinetics of a charged radiation defect
B
thermalization reported in [
24]. The hybridization rate (
53) is two orders of magnitude less than the rate of a non-hybridized charged center formation.
The hybridized charged center (
B)
quickly, during characteristic lattice times, transfers into a neutral state. Therefore, the neutralization time of a charged defect formed by a negative muon in a diamond lattice is determined by the value of (
53), and at least two order higher than that for silicon and germanium (
40).
8. Discussion
We have considered the total process of an acceptor center formation in crystals with the diamond structure, which appears as a radiation defect induced by negative pions or muons. It was shown that the evolution of this kind of radiation defects can be divided into two physically different stages. At the first stage, the negatively-charged particle is stopped in the crystal and captured by a nucleus with the charge Z (in the case of ) or at the K-shell of a muonic atom (in the case of ). Both negative pions and muons create a host nucleus with an effective charge . This strongly-charged center interacts with trace electrons and captures them. This stage of radiation defect neutralization exists because of Coulomb interaction. Numerical calculations show that this stage of neutralization strongly differs in diamond and other crystals. Namely, the radiation defect is completely neutralized in Si and Ge for a relatively short time s. In diamond, this radiation defect can be completely neutralized for a long time s.
The second stage of evolution of this radiation defect is connected with restoration of chemical bonds with the lattice as the first step and formation of an appropriate acceptor center as the final step. Chemical bonds in Si and Ge are restored for neutral radiation defects and in C chemical bonds can be formed for a single-fold charged center. Formation of a chemically-bound radiation defect is accompanied with sufficiently large energy emission. Therefore, the process of chemical bonds’ formation can be described with the help of radiation transition. Our estimates gave rather a long time for this step: s for Si and s for Ge. Charged radiation defect in C forms chemical bonds very quickly: s. We can see that this time is many orders shorter than . Therefore, the radiation defect in diamond is neutralized in the chemically-bound state. This time cannot be strictly estimated, but we can suppose that it is determined by characteristic electronic times in the lattice and must be of the order of s. We can conclude that the first step of the formation of a chemically-bound neutral radiation defect is approximately two orders shorter in diamond with respect to silicon and germanium.
The second step finishes the formation of an acceptor center in the ionized state. This step is similar in all of the above-mentioned crystals, but very complicated for obtaining good quantitative results. We can see that the chemically-bound neutral radiation defect is a cluster with unsaturated chemical bonds. This unsaturated chemical bond can be saturated if an electron of a valence band of the crystal is captured by the cluster. What kind of interaction gives an opportunity to capture a missing electron? The neutral cluster possesses a relatively large electric dipole moment, which interacts with lattice electrons. This interaction can be qualitatively described in the dipolar approximation. Our analytical and numerical calculations show that:
The total time for the acceptor center formation in the ionized state as a result of a radiation defect induced by negative pions and muons is the sum of times for all steps:
and is determined by the longest item. The final step is slower with respect to the first two steps, and we can conclude that the formation time of the ionized acceptor center is the shortest for diamond (≤
s) and the longest for silicon (≤
s). These values are comparable with characteristic times in semiconductor devices.
The ionized acceptor center is neutralized through the mechanism of Coulomb capture of the hole from the valence band. This process is well studied in many articles (see e.g., [
39,
40,
41]), and we will not concern ourselves with this problem here.