Next Article in Journal
Theoretical Study on Electronic, Optical Properties and Hardness of Technetium Phosphides under High Pressure
Next Article in Special Issue
Experimental and Theoretical Evidence for Surface-Induced Carbon and Nitrogen Fractionation during Diamond Crystallization at High Temperatures and High Pressures
Previous Article in Journal
Structural and Electrical Properties Characterization of Sb1.52Bi0.48Te3.0 Melt-Spun Ribbons
Previous Article in Special Issue
Morphology of Diamond Layers Grown on Different Facets of Single Crystal Diamond Substrates by a Microwave Plasma CVD in CH4-H2-N2 Gas Mixtures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evolution in Time of Radiation Defects Induced by Negative Pions and Muons in Crystals with a Diamond Structure

Moscow Institute of Physics and Technology, Institutsky lane 9, 141700 Dolgoprudny, Russia
Crystals 2017, 7(6), 174; https://doi.org/10.3390/cryst7060174
Submission received: 27 March 2017 / Revised: 29 May 2017 / Accepted: 9 June 2017 / Published: 14 June 2017
(This article belongs to the Special Issue Diamond Crystals)

Abstract

:
Evolution in time of radiation defects induced by negatively-charged pions and muons in crystals with diamond structures is considered. Negative pions and muons are captured by the nucleus and ionize an appropriate host atom, forming a positively-charged radiation defect in a lattice. As a result of an evolution in time, this radiation defect transforms into the acceptor center. An analysis of the full evolution process is considered for the first time. Formation of this acceptor center can be divided into three stages. At the first stage, the radiation defect interacts with a radiation trace and captures electrons. The radiation defect is neutralized completely in Si and Ge for a short time t 10 11 s, but in diamond, the complete neutralization time is very large t 10 6 s. At the second stage, broken chemical bonds of the radiation defect are restored. In Si and Ge, this process takes place for the neutral radiation defect, but in diamond, it goes for a positively-charged state. The characteristic time of this stage is t < 10 8 s for Si and Ge and t < 10 11 s for diamond. After the chemical bonds’ restoration, the positively-charged, but chemically-bound radiation defect in diamond is quickly neutralized because of the electron density redistribution. The neutralization process is characterized by the lattice relaxation time. At the third stage, a neutral chemically-bound radiation defect captures an additional electron to saturate all chemical bonds and forms an ionized acceptor center. The existence of a sufficiently big electric dipolar moment leads to the electron capture. Qualitative estimates for the time of this process were obtained for diamond, silicon and germanium crystals. It was sown that this time is the shortest for diamond (≤ 10 8 s) and the longest for silicon (≤ 10 7 s)

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 E 100 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 τ π ± 2.6 × 10 8 s and τ μ 2.2 × 10 6 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 12 , and a capture of a negative pion gives rise to the boron acceptor:
C 12 + π B 11 + n ,
where n is a decay neutron. The boron nucleus spin is I = 3 / 2 .
In silicon isotope Si 28 (92% in nature, see e.g., [9]) can capture π and than transform into the aluminum acceptor:
Si 28 + π Al 27 + n .
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:
Ge 70 + π Ga 69 + n and Ge 72 + π Ga 71 + n .
Both gallium isotopes possess a spin equal to 3/2. All isotopes B 11 , Al 27 , Ga 69 and Ga 71 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:
μ e + ν μ + ν ˜ e ,
where ν μ and ν ˜ e 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 τ μ 2.2 × 10 6 s. The large lifetime allows investigating with a high precision the processes with a characteristic time t < 10 5 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 1 s -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 Z 1 . After defect neutralization, a replacement impurity is formed, or a muonic atom, similar to an atom isotope with a nuclear charge Z 1 .
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 Z 1 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 ( Z = 6 ), 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 ( Z = 14 ) and Ge ( Z = 32 ), the negative muon is captured by a nucleus forming the pseudo-aluminum μ Al with a nuclear charge equal to Z = 13 and pseudo-gallium μ Ga with a nuclear charge equal to Z = 31 , 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, μ S R 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 ( E 1 keV) is released. A totally ionized positively-charged center and Auger electron ionization environment appear, creating secondary electrons. This process takes short a time, t 10 14 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 t 10 10 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 10 7 s 1 only at a respectively short interval of time 10 10 s. In silicon, the recombination frequency reaches approximately 10 11 s 1 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 t 10 11 s. The probability of the neutralization of a muonic atom μ B is less than 10 3 for the interval t 10 10 s. The recombination frequency in diamond sharply falls for t > 10 10 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, n s and n p are represented with equal probability, where n = 2 , 3 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 Z 1 , 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:
| ψ in = | n s n p 2 .
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 S = 3 / 2 , and a spacial part of its wave function one may represent in the view of equiprobable superposition of three Slater’s determinants:
Ψ s p 2 ( r 1 , r 2 , r 3 ) = 1 3 Ψ + 1 ( r 1 , r 2 , r 3 ) + Ψ 0 ( r 1 , r 2 , r 3 ) + Ψ 1 ( r 1 , r 2 , r 3 ) ,
where:
Ψ M ( r 1 , r 2 , r 3 ) = 1 6 ψ n s ( r 1 ) ψ n p , m ( r 1 ) ψ n p , m ( r 1 ) ψ n s ( r 2 ) ψ n p , m ( r 2 ) ψ n p , m ( r 2 ) ψ n s ( r 3 ) ψ n p , m ( r 3 ) ψ n p , m ( r 3 ) .
Here, M = m + m , and m , m = 0 , ± 1 .
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:
Ψ Cr ( r 1 , r 2 , r 3 ) = 1 6 P ( 1 ) P Ψ μ A ( P ( r 1 , r 2 , r 3 ) ) .
Here, summation is carried out over all permutations P of the valence electrons of the impurity and:
Ψ μ A ( r 1 , r 2 , r 3 ) = 1 2 ( ψ n 1 ( r 1 ) ψ n 2 ( r 2 ) ψ n 3 ( r 3 ) + ψ n 2 ( r 1 ) ψ n 3 ( r 2 ) ψ n 4 ( r 3 ) + ψ n 3 ( r 1 ) ψ n 4 ( r 2 ) ψ n 1 ( r 3 ) + ψ n 4 ( r 1 ) ψ n 1 ( r 2 ) ψ n 2 ( r 3 ) ) .
The unit vectors n a are directed from the impurity to nearest neighbors (along the directions of the chemical bonds). The one-particle functions ψ n a ( r ) are the hybridized states with directed bonds, and they can be written in a form (see e.g., [26,27]):
ψ n a ( r ) = α a ψ n s ( r ) + β n a ψ n p , n a ( r )
and satisfy normalization conditions:
ψ n a * ( r ) ψ n b ( r ) d r = δ a b .
From this condition, we can obtain relations for the coefficients in the superposition (6):
α a 2 + β n a 2 = 1 , a α a 2 = 1 , a β n a 2 = 3 , α a α b + β n a β n b cos ( n a , ^ n b ) = 0 .

4. Electron States of a Positively-Charged Radiation Defect in Diamond

A positively-charged radiation defect in diamond has the effective nuclear charge Z = 5 . Its atomic configuration contains only two electrons in an external (unfilled) electron shell, and they are in the “mixed” state [28]:
| ψ in = | 2 s 2 p .
We suppose also that the electronic configuration for the principle quantum number n = 1 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:
Ψ s p ( r 1 , r 2 ) = A 0 Ψ s p ( 0 ) ( r 1 , r 2 ) | 0 , 0 + A 1 Ψ s p ( 1 ) ( r 1 , r 2 ) M S | 1 , M S ,
where:
Ψ s p ( S ) ( r 1 , r 2 ) = 1 2 ψ 2 s ( r 1 ) ψ 2 p ( r 2 ) + ( 1 ) S ψ 2 p ( r 1 ) ψ 2 s ( r 2 ) ,
ψ 2 s ( r ) is the wave function of the 2 s -state. We assume that all p-states with different projections have equal probabilities:
ψ 2 p ( r ) = 1 3 ψ 21 , + 1 ( r ) + ψ 21 , 0 ( r ) + ψ 21 , 1 ( r ) .
S = 0 , 1 are the values of the total electron spin; | S , M S 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:
| A 0 | 2 + 3 | A 1 | 2 = 1 .
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:
Ψ Cr ( S ) ( r 1 , r 2 ) = 1 6 n a , n b Ψ n a n b ( S ) ( r 1 , r 2 ) ,
where summation is carried out over all possible directions of the chemical bond of impurity valence electrons with the nearest neighbors of the lattice,
Ψ n a n b ( S ) ( r 1 , r 2 ) = 1 2 ψ n a ( r 1 ) ψ n b ( r 2 ) + ( 1 ) S ψ n a ( r 2 ) ψ n b ( r 1 ) .
The unit vectors n a are directed from the impurity to the neighbor atoms (along the direction of chemical bonds) like in (5). One-particle wave functions ψ n a ( r ) 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 4 ) 0 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”:
d w = 2 π ψ i n | V ^ rad | ψ C r 2 δ ( E f E i ) d ν f .
The interaction operator is:
V ^ rad = e m c a p ^ a A ( r a ) ,
where A ( r ) is the vector-potential of the free radiation field.
Let us consider now only the term for one electron with a = 1 in the operator (14) to simplify the following calculations. In this case, matrix elements of the perturbation operator could be represented by the expression:
ψ i n | p ^ 1 | ψ Cr s a = ψ n s * ( r 1 ) p ^ 1 ψ n a ( r 1 ) d r 1 ψ n p , m * ( r 2 ) ψ n b ( r 2 ) d r 2 ψ n p , m * ( r 3 ) ψ n c ( r 3 ) d r 3
and:
ψ i n | p ^ 1 | ψ Cr p a = ψ n p , m * ( r 1 ) p ^ 1 ψ n a ( r 1 ) d r 1 ψ n p , m * ( r 2 ) ψ n b ( r 2 ) d r 2 ψ n s * ( r 3 ) ψ n c ( r 3 ) d r 3 ,
where the indexes are a , b , c = 1 , 2 , 3 , 4 .
If we direct the axis z | | n 1 , then the other three p-states in hybridized states (6) turn out as a result of the rotation of the state ψ n 1 , 0 ( r ) in the state with a rotation moment projection equal to zero on the axes n b . In this case, we get the opportunity to calculate easy integrals incoming in Expressions (16) and (17):
ψ n s * ( r ) ψ n a ( r ) d r = α a n s | n s = α a ;
ψ n p , m * ( r ) ψ n a ( r ) d r = β a n p , m | R ^ ( θ a , φ a ) | n 1 , 0 ,
where R ^ ( θ a , φ a ) is the rotation operator, and the matrix elements (19) are determined by the second column of the rotation matrix:
n 1 , ± 1 | R ^ ( θ a , φ a ) | n 1 , 0 = 1 2 sin θ a e i φ a , n 1 , 0 | R ^ ( θ a , φ a ) | n 1 , 0 = cos θ a .
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:
ψ n s * ( r 1 ) p ^ 1 z ψ n a ( r 1 ) d r 1 = β a ψ n s * ( r 1 ) p ^ 1 z R ^ ( θ a , φ a ) ψ n 1 , 0 ( r 1 ) d r 1 = β a cos θ a I s p ( n ) ,
ψ n p , m * ( r 1 ) p ^ 1 z ψ n a ( r 1 ) d r 1 d r 1 = α a ψ n 1 , 0 * ( r 1 ) p ^ 1 z ψ n s ( r 1 ) d r 1 = α a I s p ( n ) * ,
where:
I s p ( n ) = ψ n s * ( r ) p ^ z ψ n 1 , 0 ( r ) d r .
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.
Ψ + 1 | p ^ 1 z | Ψ n 1 , n 2 , n 3 = 1 4 3 β 1 β 2 β 3 sin θ 2 cos θ 3 e i φ 2 cos θ 2 sin θ 3 e i φ 3 I s p ( n ) + 1 4 3 α 1 β 2 α 3 sin θ 2 e i φ 2 α 2 β 3 sin θ 3 e i φ 3 I s p ( n ) * .
The rest of the three matrix elements are determined by the other possible sets of n a for the electrons with coordinates r 1 , r 2 and r 3 :
Ψ + 1 | p ^ 1 z | Ψ n 2 , n 3 , n 4 = 1 4 3 β 2 β 3 β 4 cos θ 2 sin θ 3 cos θ 4 e i φ 3 cos θ 3 sin θ 4 e i φ 4 I s p ( n ) 1 4 3 α 2 α 3 β 4 sin θ 4 e i φ 4 β 3 α 4 sin θ 3 e i φ 3 I s p ( n ) * ;
Ψ + 1 | p ^ 1 z | Ψ n 3 , n 4 , n 1 = 1 4 3 β 4 β 1 β 3 cos θ 3 I s p ( n ) α 1 α 3 I s p ( n ) * sin θ 4 e i φ 4
Ψ + 1 | p ^ 1 z | Ψ n 4 , n 1 , n 2 = 1 4 3 β 2 β 1 β 4 cos θ 4 I s p ( n ) α 1 α 4 I s p ( n ) * sin θ 2 e i φ 2 .
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 α a , β a , θ a and φ a . However, it is necessary to take into account that the system under consideration has a symmetry at less C 3 v . 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:
α a = α = 1 2 , β a = β = 3 2 .
If the vector n 4 lies in the x z plane, the angles θ a , φ a are equal:
θ 1 = 0 , θ 2 = θ 3 = θ 4 = θ , cos θ = 1 3 , φ 2 = φ 3 = 2 π 3 , φ 4 = 0 .
In the case of symmetrical structure, we get the following expressions for the matrix elements:
Ψ + 1 | p ^ 1 z | Ψ n 1 , n 2 , n 3 = i 8 β β 2 sin 2 θ I s p ( n ) α 2 sin θ I s p ( n ) * ,
Ψ + 1 | p ^ 1 z | Ψ n 2 , n 3 , n 4 = 3 8 1 1 3 β sin θ β 2 cos 2 θ I s p ( n ) α 2 I s p ( n ) * ,
Ψ + 1 | p ^ 1 z | Ψ n 3 , n 4 , n 1 = Ψ + 1 | p ^ 1 z | Ψ n 4 , n 1 , n 2 = 1 4 3 β sin θ β 2 cos θ I s p ( n ) α 2 I s p ( n ) * .
Adding up Expressions (30)–(32), we get:
Ψ + 1 | p ^ 1 z | Ψ Cr = 1 i 3 8 3 β 3 sin 2 θ sin 2 θ 2 I s p ( n ) .
The matrix elements for the state with M = 1 are calculated by a similar way:
Ψ 1 | p ^ 1 z | Ψ Cr = 1 + i 3 8 3 β 3 sin 2 θ sin 2 θ 2 I s p ( n ) .
For the state with M = 0 :
Ψ 0 | p ^ 1 z | Ψ Cr = i 2 2 β 3 sin 2 θ cos 2 θ 2 I s p ( n ) .
After substitution of the values o Parameters (28) and (29) and adding up Expressions (33)–(35), we have:
Ψ in | p ^ 1 z | Ψ Cr = β 3 4 3 sin 2 θ sin 2 θ 2 + i 2 sin 2 θ cos 2 θ 2 I s p ( n ) = 1 i 4 72 2 I s p ( n ) .
For the calculation of the integral I s p , we take the appropriate wave functions of the hydrogen-like atom with an effective nuclear charge equal to Z ˜ . In accordance with Slater [30], an effective charge is determined as Z ˜ = Z σ , 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 n = 3 , and appropriate calculations for μ Al give the following results:
I s p ( 3 ) = i 10 Z ˜ 3 6 , in   dimensional   units I s p ( 3 ) = i 10 3 6 Z ˜ a 0 ,
where a 0 is the Bohr radius.
In a germanium lattice, a gallium μ -atom μ Ga with the principle quantum number n = 4 must be formed. The unknown value of the matrix element for μ Ga is equal to:
I s p ( 4 ) = i Z ˜ 5 3 4 2 , correspondingly I s p ( 4 ) = i 5 3 4 2 Z ˜ a 0 .
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 2 S + 1 = 4 :
d w hybr = 2 π 2 ( 2 S + 1 ) e m c 2 2 2 π c 2 ω 3 Ψ in | p ^ 1 z | Ψ Cr 2 cos 2 ϑ δ ( ω i f ω ) k 2 d k d Ω ( 2 π ) 3 .
Here, ω i f 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:
w hybr = 17 × 5 2 / 2 3 3 6 Z ˜ Al 2 α 3 ω i f Si 3 . 8 × 10 8 s 1 for silicon , 17 × 15 / 2 12 Z ˜ Ga 2 α 3 ω i f Ge 4 . 9 × 10 8 s 1 for germanium ,
where α = 1 / 137 is the fine structure constant. We have substituted here the effective charge Z ˜ Al 3.5 and Z ˜ Ga 5.0 in accordance with Slater [30].

6. Formation of Neutral Center ( μ B C 4 ) 0 in Diamond

Formulae (16) and (17) must be modified for diamond in accordance with Equations (11) and (13) as follows:
ψ i n | p ^ 1 | ψ Cr s a = ψ 2 s * ( r 1 ) p ^ 1 ψ n a ( r 1 ) d r 1 ψ 2 p * ( r 2 ) ψ n b ( r 2 ) d r 2
and:
ψ i n | p ^ 1 | ψ Cr p a = ψ 2 p * ( r 1 ) p ^ 1 ψ n a ( r 1 ) d r 1 ψ 2 s * ( r 2 ) ψ n b ( r 2 ) d r 2 ,
To derive a common expression now, we consider that the Z-axis does not coincide with any of the bond directions n a . Therefore, as for (18) and (19), we have:
ψ 2 s * ( r ) ψ n b ( r ) d r = α b 2 s | 2 s = α b ;
ψ 2 p * ( r ) ψ n b ( r ) d r = 1 3 β b ( 2 cos φ b i 2 sin θ b sin φ b + cos θ b ) .
Without the reduction of the generality of the calculations, we consider a matrix element only for the z- projection of a momentum operator.
ψ 2 s * ( r 1 ) p ^ 1 z ψ n a ( r 1 ) d r 1 = β a 3 m 2 , 1 , 0 | R ^ ( n a ) | 2 , 1 , m ψ 2 s * ( r 1 ) p ^ 1 z ψ 21 , 0 ( r 1 ) d r 1 = β a 3 cos θ a I s p ( 2 ) ,
ψ 2 p * ( r 1 ) p ^ 1 z ψ n a ( r 1 ) d r 1 d r 1 = α a 3 ψ 21 , 0 * ( r 1 ) p ^ 1 z ψ 2 s ( r 1 ) d r 1 = α a 3 I s p ( 2 ) * ,
where I s p ( 2 ) is determined by Equation (23) for n = 2 .
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 S = 1 , 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 C 2 v of the ( μ B C 4 ) + cluster, we direct axes as shown in Figure 2b and introduce the following designations: θ 1 = θ 2 θ 1 , θ 3 = θ 4 θ 3 , φ 1 = π / 2 , φ 2 = 3 π / 2 , φ 3 = 0 , φ 4 = π . We have, respectively, β 1 = β 2 β 1 , α 1 = α 2 α 1 , β 3 = β 4 β 3 , α 3 = α 4 α 3 . After summation over all bonds, one obtains the expression:
ψ i n | p ^ 1 z | ψ C r = A 0 3 2 1 3 I s p * β 1 2 cos 2 θ 1 + 4 β 1 β 3 cos θ 1 cos θ 3 + β 3 2 cos 2 θ 3 + α 1 2 + 4 α 1 α 3 + α 3 2 I s p .
We can determine superposition coefficients in (6). Taking into account the relations (8) as α a 2 = 1 β a 2 and, respectively, α 1 = α 2 , β 1 = β 2 , α 3 = α 4 , β 3 = β 4 . Therefore, we have:
β 1 , 3 = 1 2 sin θ 1 , 3 , α 1 , 3 = 1 2 1 cot 2 θ 1 , 3 .
Here, the condition 2 α 1 2 + 2 α 3 2 = 1 and, respectively, cos ( n 1 , ^ n 3 ) = cos θ 1 cos θ 3 were used.
Finally, the matrix element (47) can be written in the form:
ψ i n | p ^ 1 z | ψ C r = A 0 3 2 1 3 I s p * cot 2 θ 1 + 4 cot θ 1 cot θ 3 + cot 2 θ 3 + 2 cot 2 θ 1 cot 2 θ 3 + 4 cos 2 θ 1 cos 2 θ 3 sin θ 1 sin θ 3 I s p .
For the appropriate wave functions of the hydrogen-like atom with an effective nuclear charge equal to Z ˜ , the integral I s p ( 2 ) is equal to:
I s p ( 2 ) = i Z ˜ , or   in   dimension   units I s p ( 2 ) = i Z ˜ m e 2 = i Z ˜ a 0 .
An effective charge for a boron atom is equal to Z ˜ B 2.6 [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:
w hybr = 16 3 | A 0 | 2 f 2 ( θ 1 , θ 3 ) Z ˜ 2 α 3 ω s p .
Here, ω s p 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 4 ) + .
In accordance with the matrix element (49) and the result of the matrix element (50) calculation, a configuration factor f ( θ 1 , θ 3 ) is equal to:
f ( θ 1 , θ 3 ) = 1 + 3 3 cot 2 θ 1 + cot 2 θ 3 + 4 3 cot θ 1 cot θ 3 4 cos 2 θ 1 cos 2 θ 3 sin θ 1 sin θ 3 2 .
A transition frequency ω s p and angles for a configuration factor f ( θ 1 , θ 3 ) 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 5 H* 12 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 5 H* 12 ( T d ) fragment. As a result of a structural relaxation, the defect transits into the final hybridized state | ψ C r described by the lowest in energy triplet structure of the [BC 4 H* 12 ] + ( T , C 2 v ) cluster (Figure 2b). The energy of such a transition is calculated from the difference between the total energies of levels | ψ i n and | ψ C r to be equal to 1.17 eV.
The angles in the cluster are equal to θ 1 65 , θ 3 128 . Calculations of a spin density at the center of the cluster give the value of the superposition parameter A 0 = 7 / 4 . The effective charge for 2 s and 2 p states of the boron atom is Z ˜ = 2.6 . 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:
w hybr = τ + 1 1.7 × 10 11 s 1 .
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).

7. Formation of an Ionized Acceptor Center

In this section, we will consider the process of an electron capture on the neutral radiation defect with totally restored chemical bonds and the formation of an acceptor center in the ionized state.
According to our cluster calculations, the neutral [ ( μ B)C 4 ] 0 defect has C 3 v symmetry and creates the substantial electric dipole moment, which is equal to 1.08D, in a diamond lattice. Here, D is Debye, the unit of an electric dipole moment in the atomic system of units (D = 10 18 CGSE). The dipole moment is directed along the symmetry axis. Therefore, we can suppose that any neutral cluster of a type [ ( μ A) A 4 ] 0 in crystals with a diamond structure possesses an electric dipole moment of the order of 1D.
This electric dipole moment gives rise to an interaction necessary to capture a lattice electron and the form of an ionized acceptor center. The neutral center μ A has unsaturated chemical bonds because a crystal lattice turns out to be deformed. This deformation is a reason to change a phonon spectrum and the local phonon mode appearance. Chemical bounds are saturated after the missing electron capture. The new cluster is an ionized acceptor center and possesses a local crystal symmetry. When the ionized acceptor center is formed, an appropriate phonon of the local mode is radiated, and crystal deformations disappear. Therefore, the problem is very similar to the problem of the thermalization of molecular ions in molecular crystals and cryocrystals of noble atoms (see e.g., [32,33,34]). An exact solution of the problem taking into account a crystal symmetry is scarcely possible. However, qualitative estimations can be obtained in some approximations [35,36].

7.1. Effective Hamiltonian and Interaction Operator

Detailed calculations of an ionized acceptor center formation in diamond, silicon and germanium were carried out in [35,36], and here, we will present the main results. The electric dipole moment creates a scalar potential, and the interaction energy with lattice electrons is U = e φ = e ( d r ) / ( ε r 3 ) , where ε is a dielectric penetration. Therefore, taking into account the displacement u ( r r + u ) and neglecting the changing of the denominator of the potential, we can write an electron-phonon interaction operator:
V ^ e ph e ( d u ^ ) ε r 3 = e d ε r 3 cos θ u ^ ,
The operator of radial displacements can be determined in the approach of an isotropic elastic media. For this reason, we need to study vibrations of a sphere at the center of which is placed an electric dipolar moment d . This dipolar moment creates an electric induction D , and a displacement vector u satisfies the following equation:
u ¨ i c 2 Δ u i ϰ ε 2 ρ D i D k u k = 0 ,
where c is the longitudinal sound velocity, ϰ is the dielectric susceptibility and ε = 1 + 4 π ϰ , ρ is the density of a media.
To solve Equation (55), we will take into account only radial vibrations of a deformed crystal lattice with respect to the center of the sphere. The center of this sphere coincides with our radiation defect. Therefore, we can introduce as usual for problems with central symmetry a radial displacement χ = u r . In this case, we have the more simple equation for χ :
χ ¨ c 2 + 4 ϰ d 2 ε 2 ρ 1 r 6 χ = 0 ,
where χ is the second derivative on the radial variable r. Making the one-dimensional Fourier transformation for the function χ ( r , t ) :
χ ω , k = 2 d t e i ω t 0 d r sin ( k r ) χ ( r , t ) ,
we obtain a dispersion relation in a long wavelength approximation:
( ω 2 k 2 c 2 Ω 2 ) χ k = 0 , Ω = d ε R 1 4 ϰ 6 π ρ ,
where R 1 is the radius of the first coordination sphere. More detailed calculations are represented in Appendix A.
The numerical estimates for the frequency Ω in diamond, silicon and germanium crystals are presented in Table 1 (the parameters oft the crystals were taken from [37,38]). Since a dipole moment d is considered as an unknown parameter, the numerical estimates are presented in units of Debye.
Now, we can construct an effective Hamiltonian describing the radial vibrations of a lattice and determine the operator of radial displacements u ^ in Operator (54). Consider radial vibrations in a sphere with the radius R D with boundary conditions u ( R D ) = 0 . Keeping a dependence on time t, we need to use discrete Fourier amplitudes χ n instead of Equation (57):
χ n = 2 0 R D d r sin ( k n r ) χ ( r , t ) ,
where k n = π n / R D .
Therefore, χ n and ρ χ ˙ n are generalized coordinates and the momentum of the system under consideration, respectively. The effective Hamiltonian for a system of independent oscillators can be written as:
H ^ p h = n 1 2 ω n P ^ n 2 + Q ^ n 2 .
Here,
ω n = c 2 k n 2 + Ω 2 ,
and the dimensionless generalized momentum and coordinate are defined as usual:
P ^ n = ρ χ ˙ ^ n / p 0 n , Q ^ n = χ ^ n / q 0 n ,
while the units of the generalized momentum and coordinate are equal to:
q 0 n = ω n ρ R D , p 0 n = ω n ρ R D .
Introducing, as usual, the annihilation and creation operators:
b ^ n = 1 2 Q ^ n + i P ^ n , b ^ n = 1 2 Q ^ n i P ^ n , [ b ^ n , b ^ n ] = δ n n ,
we obtain the operator of radial displacements of a lattice, which must be substituted in Equation (54), in the form:
u ^ = χ ^ r = 1 r n 2 ρ ω ( k n ) R D sin ( k n r ) b n + b n .
It is easy to see that the interaction operator (54) has a strong singularity at r 0 , which leads to the divergence of matrix elements. However, Operator (54) is obtained in a dipolar approximation; its expression is valid for respectively large values of the radius vector and is not applicable at r 0 . In this case, as usual (see e.g., [39]), the potential is taken as constant at r < R 0 , where R 0 is a certain characteristic distance, i.e.,
V ^ e ph e d / ( ε r 3 ) cos θ u ^ , for r > R 0 , e d r / ( ε R 0 4 ) cos θ u ^ , for r < R 0 .
In our case, we can take R 0 r 0 , where r 0 is the length of chemical bonds in the lattice.

7.2. Electron Capture Rate

The probability of electron capture per unit time can be calculated with the well-known Fermi golden rule:
d w i f d w capt = 2 π f | V ^ e ph | i 2 δ ( E f E i ) d ν f d k ( 2 π ) 3 .
We will consider a case of low temperature, and the initial state | i = | k | 0 ph corresponds to the free electron in the valence band with the wave-vector k and the absence of exited radial phonons. The final state is determined by the captured electron to the | n s or | n p hybridized state of the cluster and by the excitation of the radial phonon with the wavevector k ph :
| f = | n s | n p | k ph .
Here, n = 2 , 3 and 4 for C, Si and Ge, respectively. Equation (66) takes into account that all electrons having wavevectors in the interval from k to k + d k can be captured by the cluster.
The volume element of the final states in the case of one-dimensional motion is equal to:
d ν f = R D d k ph 2 π d o ,
where d o is the element of the solid angle into which a phonon is emitted. The energies of the initial and final states are:
E i = E 0 + 2 k 2 2 m * , E f = E + ω ( k ph ) , E 0 E = ,
m * is the effective mass of the electron. To obtain numerical results, the ionization energy of the cluster (acceptor) is taken for .
Integration of (66) over k gives:
d w capt = m * ( 2 π ) 2 2 3 ρ e d ε 2 e f F ( r ) e i k ph r k 0 2 k 0 ω ( k ph ) d k ph d Ω ,
where a vector of the final state | f e takes into account electron states only, and according to Equation (65),
F ( r ) = r 4 , for r > R 0 , R 0 4 , for r < R 0 .
Here:
k 0 2 = 2 m * ( ω ph ) .
For the further analysis, it is convenient to introduce the dimensionless parameters and variables:
x = Z ˜ n a 0 r , k ˜ ph = n a 0 Z ˜ k ph , k ˜ 0 = n a 0 Z ˜ k 0 ,
where n = 2 , 3 and 4 for C, Si and Ge, respectively. The estimated characteristic parameters for considered crystals are summarized in Table 2.
All considered crystals are anisotropic, and the sound velocity c in the dispersion relation (58) depends on the direction of a phonon propagation. To take this into account, we will use, as usual, the average value of the longitudinal velocity of sound c , neglecting the effects of anisotropy.
We define a characteristic frequency:
ω 0 = c Z ˜ n a 0
and corresponding dimensionless frequencies ω ˜ = ω ph / ω 0 , ˜ = / ω 0 and Ω ˜ = Ω / ω 0 . Taking into account the dispersion relation of the radial phonons (58), we can write the electron capture rate in the form:
w capt = G n d Ω 4 π Ω ˜ ˜ k ˜ 0 A ( k ˜ ph , k ˜ 0 ) 2 ( ˜ ω ˜ ) ( ω ˜ 2 Ω ˜ 2 ) d ω ˜ ,
where:
G n = 1 6 π ρ m * 2 m * c 3 / 2 e d ε 2 Z ˜ n a 0 13 / 2 .
The dimensionless matrix element A is determined by Function (68):
A ( k ˜ ph , k ˜ 0 ) = e f F ( x ) e i k ˜ ph x k ˜ 0 .
Equation (72) shows that the electron can be captured both in s- and p-states of the hybridized state of the cluster. As was shown in [36], the capture rate to the p-state w p is approximately two orders less than the capture rate to s-state w s . Therefore, in the following, we can consider the electron capture to the s-state. In this case, the matrix element (74) is given by the expression:
A n s ( k ˜ ph , k ˜ 0 ) = 1 2 π k ˜ 0 0 F ( x ) R n 0 ( x ) e i k ˜ ph x sin ( k ˜ 0 x ) x d x = k ˜ 0 1 ( I 1 ( n ) ( ω ˜ ) i I 2 ( n ) ( ω ˜ ) ) ,
where I 1 ( ω ˜ ) and I 2 ( ω ˜ ) are the real and imaginary parts of the matrix element, respectively. Finally, the capture rate (72) has the form:
w capt w s = G n Ω ˜ ˜ ( I 1 ( n ) ( ω ˜ ) ) 2 + ( I 2 ( n ) ( ω ˜ ) ) 2 ( ˜ ω ˜ ) ( ω ˜ 2 Ω ˜ 2 ) d ω ˜ ,
Obtained Formula (76) determines the formation rate of the ionized acceptor center through the capture of the electron of the medium on the neutral radiation defect induced by the negative muon or pion in crystals. Analytical expressions depend on well-known characteristic parameters of the medium except for two parameters of the cluster, namely, the electric dipolar moment d and the parameter R 0 in Equation (65). The dependence on d of the integral in Equation (75) is weak and can assume that w capt d 2 . Unfortunately, the dependence of the results on the parameter R 0 is more critical, because the matrix elements exponentially depend on the lower integration limit. Nevertheless, the parameter R 0 cannot be smaller than the length of the chemical bound in the lattice. We can determine the upper limit also R 0 R 1 , the radius of the first coordination sphere.
Numerical calculation were performed for several different values of the uncertain parameters d and R 0 . The results are summarized in Table 3.
According to the obtained results in the considered range of the uncertain parameters R 0 and d, the spread of the estimates of the capture rates is about two orders of magnitude and seems at first unsatisfactory. However, it is not surprising because both the dipole approximation and the dispersion relation for the radial phonons (58) are quite rough for r R 1 . Nevertheless, the considered interaction mechanism is fairly justified and describes the process of the electron capture on the neutral defect of the lattice with the formation of the ionized state of the acceptor center.

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 Z 1 . 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 τ n 10 11 s. In diamond, this radiation defect can be completely neutralized for a long time τ n > 10 6 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: τ hybr 2.6 × 10 9 s for Si and τ hybr 2.0 × 10 9 s for Ge. Charged radiation defect in C forms chemical bonds very quickly: τ hybr 0.6 × 10 11 s. We can see that this time is many orders shorter than τ n . 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 10 10 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:
0 . 5 × 10 10 s τ capt 2 × 10 8 s for diamond , 1 . 4 × 10 8 s τ capt 2 × 10 7 s for silicon , 0 . 5 × 10 9 s τ capt 1 . 7 × 10 8 s for germanium .
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:
τ ac = τ n + τ hybr + τ capt ,
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 (≤ 2 × 10 8 s) and the longest for silicon (≤ 2 × 10 7 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.

9. Conclusions

The obtained results will be useful both for μ SR experiments and research of different radiation defects in semiconductors. The considered approach can be applied to crystals with the sphalerite-type structure (A I I I B V semiconductors, e.g., GaAs, InSb, CdS), which are widely used in electronic devices. Unfortunately, these crystals are more complicated for analysis because of the large variety of possible impurity centers. In addition, the model of chemical bonds must be modified for some calculations.

Acknowledgments

I acknowledge my colleagues Leonid E. Fedichkin and Leonid P. Sukhanov for fruitful discussions. This work is supported by the Ministry of Education and Science of the Russian Federation 16.7162.2017/8.9.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The function χ ( r , t ) can be represented in the form of the “one-dimensional” Fourier-transform:
χ ( r , t ) = 2 d ω d k ( 2 π ) 2 χ ω , k e i ω t sin ( k r ) ,
χ ω , k = 2 d t e i ω t 0 d r sin ( k r ) χ ( r , t ) .
An integral equation can be obtained for the Fourier amplitude given by Equation (A1) taking into account that the expression for the force in Equation (54) is valid for r > R 1 , where R 1 is the radius of the first coordination sphere, i.e.,
( ω 2 k 2 c 2 ) χ k = 4 ϰ d 2 ε 2 ρ R 1 d r e i k r 1 r 6 χ ( r , t ) .
After passage of the Fourier amplitude on the right-hand side of Equation (A3), this equation is reduced to the form:
( ω 2 k 2 c 2 ) χ k = 4 ϰ d 2 π ε 2 ρ × 0 d k k 2 χ k R 1 d r r 6 cos ( k k ) r cos ( k + k ) r .
Another serious simplification should be done in Equation (A4). The maximum contribution to the second integral comes from wavevectors k k . For this reason, the term cos ( k + k ) r can be neglected, and the Fourier amplitude χ k can be taken at the point k. In particular, this means that that the main contribution to the integral comes from the values k r 1 . Correspondingly, in the accepted approximation, the integral on the right-hand side can be modified to the form:
0 d k k 2 χ k R 1 d r r 6 cos ( k k ) r χ k R 1 d r r 6 0 r 1 k 2 d k = χ k 24 R 1 8 .
After that, the simple dispersion relation (58) is obtained.

References

  1. Honniger, F.; Fretwurst, E.; Lindstrom, G. DLTS measurements of radiation induced defects in epitaxial an MCz silicon detectors. Nucl. Instrum. Methods Phys. Res. Sect. A 2007, 583, 104–108. [Google Scholar] [CrossRef]
  2. Grilj, V.; Skukan, N.; Jakšić, M. Irradiation if thin diamond detectors and radiation hardness tests using MeV protons. Nucl. Instrum. Methods Phys. Res. Sect. B 2013, 306, 191–194. [Google Scholar] [CrossRef]
  3. Guthoff, M.; de Boer, W.; Muller, S. Simulation of beam induced lattice defects of diamond detectors using FLUKA. Nucl. Instrum. Methods Phys. Res. Sect. A 2014, 735, 223–228. [Google Scholar] [CrossRef]
  4. Sato, Y.; Shimaoka, T.; Kaneko, J.H. Radiation hardness of a single crystal CVD diamond detector for MeV energy protons. Nucl. Instrum. Methods Phys. Res. Sect. A 2015, 784, 147–150. [Google Scholar] [CrossRef]
  5. Grilj, V.; Skukan, N.; Jakšić, M. The evaluation of radiation damage parameter for CVD diamond. Nucl. Instrum. Methods Phys. Res. Sect. B 2016, 372, 161–164. [Google Scholar] [CrossRef]
  6. Makgato, T.H.; Sideras-Haddad, E.; Ramos, M.A. Magnetic properties of point defects in proton irradiated diamond. J. Magn. Magn. Mater. 2016, 413, 76–80. [Google Scholar] [CrossRef]
  7. Vins, V.G. New radiation induced defects in HPHT synthetic diamonds. Diam. Relat. Mater. 2005, 14, 364–368. [Google Scholar] [CrossRef]
  8. Dresselhaus, M.S.; Kalish, R. Ion implantation in Diamond, Graphite and Related Materials; Springer: Berlin, Germany, 1992. [Google Scholar]
  9. Vonsovsky, S.V. Magnetizm of Microparticles; Nauka: Moscow, Russia, 1973; pp. 1–197. (In Russian) [Google Scholar]
  10. Smilga, V.P.; Belousov, Y.M. The Muon Method in Science; Nova Science: New York, NY, USA, 1994; pp. 1–420. [Google Scholar]
  11. Waisenberg, A.O. Mu-Meson; Nauka: Moscow, Russia, 1964; pp. 1–400. (In Russian) [Google Scholar]
  12. Dzhuraev, A.A.; Evseev, V.S. Mechanism of mu-meson depolarization in molecular condensed media. Sov. Phys. JETP 1972, 35, 615–819. [Google Scholar]
  13. Dzhuraev, A.A.; Evseev, V.S.; Myasishcheva, G.G. Depolarization of negative muons in solids. Sov. Phys. JETP 1972, 35, 748–752. [Google Scholar]
  14. Dzhuraev, A.A.; Evseev, V.S. Depolarization of negative muons in condensed molecular media. Sov. Phys. JETP 1974, 39, 207–211. [Google Scholar]
  15. Gorelkin, V.N.; Smilga, V.P. On the theory of precession of the polarization vector of μ-mesons in mu-nucleonic atoms. Sov. Phys. JETP 1974, 66, 586–591. [Google Scholar]
  16. Gorelkin, V.N.; Grebinnik, V.G.; Gritsaj, K.I. μSR-study of the behavior of impurity atoms in silicon. Phys. At. Nucl. 1993, 56, 1316–1319. [Google Scholar]
  17. Mamedov, T.N.; Duginov, V.N.; Grebinnik, V.G. Investigation of the behavior of the impurity atoms in Si by μSR-method. Hyperfine Interact. 1994, 86, 717–722. [Google Scholar] [CrossRef]
  18. Mamedov, T.N.; Chaplygin, I.L.; Duginov, V.N. Anomalous frequency shift of negative muon spin precession in n-type silicon. Hyperfine Interact. 1997, 105, 345–349. [Google Scholar] [CrossRef]
  19. Mamedov, T.N.; Gritsaj, K.I.; Stoykov, A.V. μSR investigations in silicon. Phys. B 2000, 289–290, 574–577. [Google Scholar] [CrossRef]
  20. Mamedov, T.N. Study of Boron Acseptor Center in Synthetic Diamond by μSR-Technique; JINR Preprint P14-2007-12; JINR: Dubna, Russia, 2007. (In Russian) [Google Scholar]
  21. Mamedov, T.N.; Baturin, A.S.; Blank, V.D. Non-equilibrium charge carrier dynamics in synthetic diamond studied by μSR-method. Diam. Relat. Mater. 2008, 17, 1221–1224. [Google Scholar] [CrossRef]
  22. Mamedov, T.N.; Baturin, A.S.; Gritsaj, K.I. Muonic atom as an acceptor centre in diamond. J. Phys. Conf. Ser. 2014, 551, 012046. [Google Scholar] [CrossRef]
  23. Baturin, A.S.; Gorelkin, V.N.; Mamedov, T.N.; Soloviev, V.R. Absolute negative mobility of charge carriers in diamond and interpretation of μSR experiments. Phys. B 2006, 374–375, 340–346. [Google Scholar] [CrossRef]
  24. Antipov, S.A.; Belousov, Y.M.; Soloviev, V.R. Dynamics of charge carriers at the place of the formation of a muonic atom in diamond and silicon. JETP 2012, 115, 866–875. [Google Scholar] [CrossRef]
  25. Belousov, Y.M. Formation of radiation defect induced by negative muons in crystals with diamond structure. J. Phys. Conf. Ser. 2012, 343, 012013. [Google Scholar] [CrossRef]
  26. Pauling, L. Natureof the Chemical Bond, 3rd ed.; Cornell University Press: Ithaca, NY, USA, 1960. [Google Scholar]
  27. Coulson, C.A. Valence; Oxford University Press: Oxford, UK, 1961. [Google Scholar]
  28. Belousov, Y.M.; Sukhanov, L.P. Formation of chemically bound positively charged radiation defect induced by negative muon in diamond crystals. Diam. Relat. Mater. 2015, 58, 10–15. [Google Scholar] [CrossRef]
  29. Landau, L.D.; Lifshitz, E.M. Quantum Mechanics. Nonrelativistic Theory, 3rd ed.; Pergamon: Oxford, UK; New York, NY, USA, 1977. [Google Scholar]
  30. Slater, J.C. Atomic Shielding Constants. Phys. Rev. 1930, 36, 57–64. [Google Scholar] [CrossRef]
  31. Groshev, G.E.; Sukhanov, L.P. Quantum-chemical cluster models of the Si(111)/SiO2 interface. Theor. Exp. Chem. 1990, 26, 268–275. (In Russian) [Google Scholar] [CrossRef]
  32. Belousov, Y.M.; Smilga, V.P. Theory of muon spin relaxation in condensed phases of hydrogen isotopes. JETP 1994, 79, 811–818. [Google Scholar]
  33. Belousov, Y.M. Theoretical description of the muon spin depolarization in the crystalline phase of 3He. JETP 2007, 104, 215–229. [Google Scholar] [CrossRef]
  34. Belousov, Y.M.; Smilga, V.P. Theory of muon spin relaxation in condensed phases of hydrogen isotopes. Hyperfine Interact. 1997, 106, 19–25. [Google Scholar] [CrossRef]
  35. Belousov, Y.M. Process of an acceptor center μB formation in diamond crystals. J. Phys. Conf. Ser. 2014, 551, 012047. [Google Scholar] [CrossRef]
  36. Belousov, Y.M. Process of negative-muon-induced formation of an ionized acceptor center (μA) in crystals with diamond structure. JETP 2016, 123, 1160–1169. [Google Scholar] [CrossRef]
  37. Madelung, O. Semiconductors—Basic Data., 3rd ed.; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  38. Adachi, S. Properties of Group-IV, III-V and II-VI Semiconductors. Wiley Series in Materials for Electronic & Optoelectronic Applications; John Wiley & Sons: Chichester, UK, 2005. [Google Scholar]
  39. Abakumov, V.N.; Perel, V.I.; Yasievich, I.N. Nonradiative Recombination in Semiconductors; SPb.Institute of Nuclear Physics RAS: Amsterdam, The Netherlands, 1997. [Google Scholar]
  40. Lax, M. Cascade capture of electron in solids. Phys. Rev. 1960, 119, 1502–1523. [Google Scholar] [CrossRef]
  41. Abakumov, V.N.; Yasievich, I.N. Cross section for recombination of an electron with a positively charged center in a semiconductor. Sov. Phys. JETP 1976, 44, 345–349. [Google Scholar]
Figure 1. Radiation defect μ Al state in a silicon lattice: (a) all bonds with host atoms are broken in the initial state, and the defect state is determined by the “mixed” function 3 s 3 p 2 ; (b) three electrons of the impurity form chemical bonds with host atoms in the final state, and an unsaturated (broken) bond is equiprobable for four nearest neighbors of the cluster (AlSi 4 ) 0 .
Figure 1. Radiation defect μ Al state in a silicon lattice: (a) all bonds with host atoms are broken in the initial state, and the defect state is determined by the “mixed” function 3 s 3 p 2 ; (b) three electrons of the impurity form chemical bonds with host atoms in the final state, and an unsaturated (broken) bond is equiprobable for four nearest neighbors of the cluster (AlSi 4 ) 0 .
Crystals 07 00174 g001
Figure 2. Structure of the radiation ( μ B) + defect in a diamond lattice: (a) the ( μ BC 4 ) + cluster has a T d symmetry in the initial triplet state; (b) it has a lower C 2 v symmetry with angles θ 1 = 65 and θ 3 = 128 in the final triplet state. The internuclear distances are given in Å. The spin densities on atoms are marked by bold type.
Figure 2. Structure of the radiation ( μ B) + defect in a diamond lattice: (a) the ( μ BC 4 ) + cluster has a T d symmetry in the initial triplet state; (b) it has a lower C 2 v symmetry with angles θ 1 = 65 and θ 3 = 128 in the final triplet state. The internuclear distances are given in Å. The spin densities on atoms are marked by bold type.
Crystals 07 00174 g002
Table 1. Physical parameters of the C, Si and Ge crystals and the estimate of the frequency Ω .
Table 1. Physical parameters of the C, Si and Ge crystals and the estimate of the frequency Ω .
CrystalDensity ρ , g cm 3 R 1 , 10 8 cmDielectric Constant Ω , 10 9 s 1 Ionization Energy of the Acceptor, eV
Diamond3.513.575.758.1 d μ B, 0.37
Silicon2.335.4311.971.4 d μ Al, 0.069
Germanium5.3235.6616.00.66 d μ Ga, 0.011
Table 2. Characteristic parameters estimated for C, Si and Ge.
Table 2. Characteristic parameters estimated for C, Si and Ge.
Crystal r 0 , 10 8 cm Z ˜ r ˜ 0 ω 0 , 10 14 c 1 k 0 , max , 10 7 cm 1 k ˜ 0 , max r ˜ 0 k ph , max , 10 8  cm 1 k ˜ ph , max r ˜ 0
C1.542.63.783.323.10.484.215.9
Si2.343.55.151.051.340.311.182.78
Ge2.445.05.751.240.540.130.320.77
Table 3. Estimate of the formation rate of ionized acceptor centers ( μ A ) in the C, Si and Ge crystals.
Table 3. Estimate of the formation rate of ionized acceptor centers ( μ A ) in the C, Si and Ge crystals.
d x 0 , C w capt , C · 10 8 s 1 x 0 , Si w capt , Si · 10 8 s 1 x 0 , Ge w capt , Ge · 10 8 s 1
1.03.7822.05.120.715.7614.3
0.53.785.95.120.195.763.8
0.23.780.985.120.0345.760.65
1.05.52.27.50.0957.52.4
0.55.50.587.50.0247.50.64
0.25.50.097.52.6 · 10 3 7.50.11

Share and Cite

MDPI and ACS Style

Belousov, Y.M. Evolution in Time of Radiation Defects Induced by Negative Pions and Muons in Crystals with a Diamond Structure. Crystals 2017, 7, 174. https://doi.org/10.3390/cryst7060174

AMA Style

Belousov YM. Evolution in Time of Radiation Defects Induced by Negative Pions and Muons in Crystals with a Diamond Structure. Crystals. 2017; 7(6):174. https://doi.org/10.3390/cryst7060174

Chicago/Turabian Style

Belousov, Yury M. 2017. "Evolution in Time of Radiation Defects Induced by Negative Pions and Muons in Crystals with a Diamond Structure" Crystals 7, no. 6: 174. https://doi.org/10.3390/cryst7060174

APA Style

Belousov, Y. M. (2017). Evolution in Time of Radiation Defects Induced by Negative Pions and Muons in Crystals with a Diamond Structure. Crystals, 7(6), 174. https://doi.org/10.3390/cryst7060174

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