# Modeling Surface Recombination at the p-TypeSi/SiO

Authors: Moustafa Y. Ghannam; Husain A. Kamal
Journal: Advances in Condensed Matter Physics
 (2014)

---

## Abstract

An integral model is proposed for recombination at the silicon/silicon dioxide (Si/SiO2) interface of thermally oxidized p-type silicon via P
b amphoteric centers associated with surface dangling bonds. The proposed model is a surface adaptation of a model developed for bulk recombination in amorphous silicon based on Sah-Shockley statistics which is more appropriate for amphoteric center recombination than classical Shockley-Read-Hall statistics. It is found that the surface recombination via amphoteric centers having capture cross-sections larger for charged centers than for neutral centers is distinguished from Shockley-Read-Hall recombination by exhibiting two peaks rather than one peak when plotted versus surface potential. Expressions are derived for the surface potentials at which the peaks occur. Such a finding provides a firm and plausible interpretation for the double peak surface recombination current measured in gated diodes or gated transistors. Successful fitting is possible between the results of the model and reported experimental curves showing two peaks for surface recombination velocity versus surface potential. On the other hand, if charged and neutral center capture cross-sections are equal or close to equal, surface recombination via amphoteric centers follows the same trend as Shockley-Read-Hall recombination and both models lead to comparable surface recombination velocities.


---

## Body

## 1. Introduction

Recombination at the Si/SiO2 interface has been classically modeled using the two-charge-state Shockley-Read-Hall (SRH) statistics [1] based on Fermi-Dirac occupation probability. Defects at the Si/SiO2 interface, however, have been identified by Electron Paramagnetic and Electron Spin Resonance (EPR, ESR) as dangling bonds, referred to as P
b  centers, precisely located at the Si/SiO2 interface [2–4], and ascertained to be amphoteric [3–6]. Recombination via amphoteric centers is better described by Sah-Shockley recombination statistics [7] which has been used to model bulk recombination due to dangling bond defects in amorphous silicon (a-Si) [8] taking into account the multicharge nature of the centers and the energy correlation between them. Recently, models for recombination via amphoteric dangling bond defects at the a-Si/c-Si interface of heterojunction silicon solar cells using Sah-Shockley statistics have been developed [9, 10], and their application to recombination at the Si/SiO2 interface has been proposed [10].In the present work a model for surface recombination via amphoteric dangling bondP
b centers at the Si/SiO2 interface is developed based on the available model for bulk recombination in a-Si [8]. The surface recombination velocity is calculated using the proposed model and using SRH classical recombination adapted to the surface. In both models, recombination centers are assumed to form a continuum distributed throughout the energy gap and carrier dependence on surface potential at the thermally oxidized Si/SiO2 interface is taken into account. Comparison between the results of both amphoteric center and SRH models is presented. Finally, by fitting previously reported experimental data with the results of the proposed amphoteric center model, an attempt is made to give a convincing explanation for the double peak surface recombination velocity extracted from recombination current in gated diodes and gated transistors.


## 2. Surface Recombination Physical Parameters and Models

Let us assume that minority electrons are injected into a p-type silicon region through a pn junction atx
=
0, travel through the region, and reach its oxide passivated surface at x
=
W where they recombine with holes. The band diagram at the p-type oxide passivated surface sketched in Figure 1 shows a band bending that extends over a shallow surface space charge layer in the c-Si substrate which results in an electric field extending throughout the layer (surface field layer) and in a surface potential ψ
s. The surface potential is positive if the bending is downwards (surface depletion) and negative if the bending is upward (accumulation).Figure 1
Band diagram at a p-type back surface showing band bending, surface potential, and surface field (space charge) layer.The surface carrier concentrations depend on the surface potential such that at low level injection(1a)
p
s
=
p
(
W
)
e
-
q
ψ
s
/
k
T
≈
p
o
e
-
q
ψ
s
/
k
T
=
N
A
e
-
q
ψ
s
/
k
T
,
(1b)
n
s
=
n
(
W
)
e
q
ψ
s
/
k
T
=
[
n
o
+
Δ
n
(
W
)
]
e
q
ψ
s
/
k
T
≈
Δ
n
(
W
)
e
q
ψ
s
/
k
T
=
n
o
e
V
F
/
V
T
e
q
ψ
s
/
k
T
,
(1c)
p
s
×
n
s
=
p
(
W
)
n
(
W
)
≈
p
(
W
)
Δ
n
(
W
)
=
N
A
×
n
o
exp
⁡
(
V
F
V
T
)
=
n
i
2
exp
⁡
(
V
F
V
T
)
,where p
o
=
N
A is the p-type doping concentration and Δ
n
(
W
) is the excess minority electron concentration at the inner edge of the surface space charge region associated with a split V
F in the Fermi Energy. In (1c) bulk recombination in the very thin surface space charge (field) region is neglected which is usually a valid approximation. Hence, the surface recombination rate U
s can be defined at the inner edge of the surface space charge region and is related to the effective surface recombination velocity S
eff through
(2)
S
eff
=
U
s
Δ
n
(
W
)
.


### 2.1. Shockley-Read-Hall (SRH) Surface Recombination Model

Surface recombination has been traditionally treated using Shockley-Read-Hall (SRH) statistics which uses Fermi-Dirac occupation probability function (occupancy),f
(
E
t
), for recombination centers at energy E
t. At equilibrium f
o
(
E
t
) is given by(3a)
f
o
(
E
t
)
=
1
1
+
e
(
E
t
-
E
F
)
/
k
T
which states that the probability of occupation of a recombination center at energy E
t is smaller than 1/2 if E
t is larger than the Fermi Energy E
F, greater than 1/2 if E
t
<
E
F, and exactly equal 1/2 when E
t
=
E
F.The SRH surface recombination rate is expressed by an expression similar to the SRH bulk recombination rate with replacing the bulk electron and hole concentrationsn and p by the surface potential dependent surface carrier concentrations n
s and p
s. Unlike previous reports the interface recombination centers in the present work do not occur at a single discrete energy level with a density N
i
t (/cm2) but form a continuum in the energy bandgap with a state density D
i
t
(
E
t
) (/cm2/eV). Consequently, the SRH surface recombination rate U
s is to be determined from the integration over all energies in the bandgap from the valence band edge E
v up to the conduction band edge E
c
(3b)
U
s
,
SRH
=
∫
E
v
E
c
v
th
D
i
t
(
E
t
)
(
p
s
n
s
-
n
i
2
)
d
E
t
(
n
s
+
n
i
e
(
E
t
-
E
i
)
/
k
T
)
/
σ
p
+
(
p
s
+
n
i
e
(
E
i
-
E
t
)
/
k
T
)
/
σ
n
.
Using (2), the effective surface recombination velocity can then be obtained from
(3c)
S
SRH
=
U
s
,
SRH
Δ
n
(
W
)
≈
∫
E
v
E
c
v
th
D
i
t
(
E
t
)
p
o
d
E
t
(
n
s
+
n
i
e
(
E
t
-
E
i
)
/
k
T
)
/
σ
p
+
(
p
s
+
n
i
e
(
E
i
-
E
t
)
/
k
T
)
/
σ
n
,where v
th is the thermal velocity (≈107 cm/s), n
i is the intrinsic carrier concentration, E
i is the intrinsic Fermi Energy (midgap energy level), k is Boltzmann constant, T is the absolute temperature, and σ
n and σ
p represent the electron and hole capture cross-sections of the centers.


### 2.2. Surface Recombination Amphoteric Dangling Bond Model

It has been argued that the two-state charge Fermi-Dirac occupation function has some limitations in a-Si abundantly hosting amphoteric dangling bonds [8, 11] and that multicharge Sah-Shockley correlated electron statistics [7] is more appropriate for this case. This may also apply to recombination at the Si/SiO2 where defects are mainly due to dangling bond amphoteric P
b centers. Hence we propose to use the amphoteric center recombination model developed for bulk recombination in a-Si after proper adaptation to the surface as done for SRH surface recombination.Since the amphoteric recombination model is not very well known to many readers, the main formulation aspects detailed in [8] are repeated in the appendix for convenience. Besides the flaw distribution density D
i
t
(
E
t
), the model parameters include electron capture cross-sections of positive and neutral centers σ
n
+ and σ
n
0 and hole capture cross-sections of negative and neutral centers σ
p
- and σ
p
0. The total surface recombination rate via a continuum of amphoteric centers in the energy gap is obtained by integrating the recombination rate at a single level U
s
,
amph
(
E
t
) over all the energies from E
v to E
c:(4a)
U
s
,
amph
=
∫
U
s
,
amph
(
E
t
)
·
d
E
t
,
and the effective surface recombination velocity S
amph follows from (2):
(4b)
S
amph
=
U
s
,
amph
Δ
n
(
W
)
.


## 2.1. Shockley-Read-Hall (SRH) Surface Recombination Model

Surface recombination has been traditionally treated using Shockley-Read-Hall (SRH) statistics which uses Fermi-Dirac occupation probability function (occupancy),f
(
E
t
), for recombination centers at energy E
t. At equilibrium f
o
(
E
t
) is given by(3a)
f
o
(
E
t
)
=
1
1
+
e
(
E
t
-
E
F
)
/
k
T
which states that the probability of occupation of a recombination center at energy E
t is smaller than 1/2 if E
t is larger than the Fermi Energy E
F, greater than 1/2 if E
t
<
E
F, and exactly equal 1/2 when E
t
=
E
F.The SRH surface recombination rate is expressed by an expression similar to the SRH bulk recombination rate with replacing the bulk electron and hole concentrationsn and p by the surface potential dependent surface carrier concentrations n
s and p
s. Unlike previous reports the interface recombination centers in the present work do not occur at a single discrete energy level with a density N
i
t (/cm2) but form a continuum in the energy bandgap with a state density D
i
t
(
E
t
) (/cm2/eV). Consequently, the SRH surface recombination rate U
s is to be determined from the integration over all energies in the bandgap from the valence band edge E
v up to the conduction band edge E
c
(3b)
U
s
,
SRH
=
∫
E
v
E
c
v
th
D
i
t
(
E
t
)
(
p
s
n
s
-
n
i
2
)
d
E
t
(
n
s
+
n
i
e
(
E
t
-
E
i
)
/
k
T
)
/
σ
p
+
(
p
s
+
n
i
e
(
E
i
-
E
t
)
/
k
T
)
/
σ
n
.
Using (2), the effective surface recombination velocity can then be obtained from
(3c)
S
SRH
=
U
s
,
SRH
Δ
n
(
W
)
≈
∫
E
v
E
c
v
th
D
i
t
(
E
t
)
p
o
d
E
t
(
n
s
+
n
i
e
(
E
t
-
E
i
)
/
k
T
)
/
σ
p
+
(
p
s
+
n
i
e
(
E
i
-
E
t
)
/
k
T
)
/
σ
n
,where v
th is the thermal velocity (≈107 cm/s), n
i is the intrinsic carrier concentration, E
i is the intrinsic Fermi Energy (midgap energy level), k is Boltzmann constant, T is the absolute temperature, and σ
n and σ
p represent the electron and hole capture cross-sections of the centers.


## 2.2. Surface Recombination Amphoteric Dangling Bond Model

It has been argued that the two-state charge Fermi-Dirac occupation function has some limitations in a-Si abundantly hosting amphoteric dangling bonds [8, 11] and that multicharge Sah-Shockley correlated electron statistics [7] is more appropriate for this case. This may also apply to recombination at the Si/SiO2 where defects are mainly due to dangling bond amphoteric P
b centers. Hence we propose to use the amphoteric center recombination model developed for bulk recombination in a-Si after proper adaptation to the surface as done for SRH surface recombination.Since the amphoteric recombination model is not very well known to many readers, the main formulation aspects detailed in [8] are repeated in the appendix for convenience. Besides the flaw distribution density D
i
t
(
E
t
), the model parameters include electron capture cross-sections of positive and neutral centers σ
n
+ and σ
n
0 and hole capture cross-sections of negative and neutral centers σ
p
- and σ
p
0. The total surface recombination rate via a continuum of amphoteric centers in the energy gap is obtained by integrating the recombination rate at a single level U
s
,
amph
(
E
t
) over all the energies from E
v to E
c:(4a)
U
s
,
amph
=
∫
U
s
,
amph
(
E
t
)
·
d
E
t
,
and the effective surface recombination velocity S
amph follows from (2):
(4b)
S
amph
=
U
s
,
amph
Δ
n
(
W
)
.


## 3. Numerical Results

A MATLAB code is developed and used for calculating the effective surface recombination velocityS
SRH and S
amph, according to (3c) and (4b), at the interface of an oxide passivated p-type Si surface having a doping concentration of 1.5
×
10
16/cm3. The recombination centers, although initially having a U shape distribution [5], are assumed to be uniformly distributed in the energy gap as it happens to be after annealing in hydrogen or forming gas [12–14]. After this technological step the midgap center density drops from the 1011/cm2/eV range to the 1010/cm2/eV and the correlation energy between the centers is in the range 0 to 0.5 eV. The capture cross-sections of amphoteric centers are reported to be in the range 10−16 cm2 [15, 16] and it is widely accepted that σ
n can be up to 100 times larger than σ
p [17]. Values for D
i
t, E
corr, σ
n, and σ
p close to these typical values are used in the calculations carried out here.As detailed in [8] and repeated in the appendix, the correlation energy E
corr between amphoteric centers influences the occupation probability which affects the recombination rate and velocity. Therefore two extreme E
corr values will be considered: (
1
)  E
corr
=
0 and E
corr
=
0.4 eV. In the first case and according to (A.7), the equilibrium occupation probability of centers considered to be positively charged is high for centers having E
t
>
E
F, equals 25% for centers with E
t
=
E
F, and vanishes for centers with E
t more than 3
k
T below E
F. The probability of occupation for centers considered to be negatively charged is high for centers having E
t
<
E
F, equals 25% for centers with E
t
=
E
F, and vanishes for centers with E
t more than 3
k
T above E
F. Finally, the probability of occupation for centers considered to be neutral is symmetric around E
F, equals 50% for E
t
=
E
F, and vanishes when E
t is more than 6
k
T below or above E
F. For the extreme case E
corr
=
0.4 eV the center of symmetry of the occupation probability is moved from E
t
=
E
F to E
F
-
E
corr
/
2, that is, to 0.2 eV below E
F. In that case, the probability of occupation of positive centers is 1/3 at E
t
=
E
F, is higher than 1/3 for centers having E
t
>
E
F, and vanishes for E
t more than 6
k
T below E
F. The probability of occupation of negative centers is equal to 1
/
3 for centers with E
t
=
E
F
-
0.4
eV, is higher for centers with lower energies, and vanishes for centers with E
t more than 4
k
T above E
F. Finally, the probability of occupation of neutral centers is 2/3 for E
t
=
E
F and E
t
=
E
F
-
0.4
eV and almost one between E
t
=
E
F
-
0.3
eV and E
F
-
0.1
eV.


### 3.1. Effective Surface Recombination Velocity versus Excess Carrier Concentration

The effective surface recombination velocityS
eff is calculated using SRH and amphoteric center recombination models and plotted as a function of the excess electron concentration Δ
n
(
W
) at the edge of the surface field region. Since nonnegligible positive fixed charges are present in the passivating oxide layer, the p-type surface would be depleted which creates a surface field region, as shown in Figure 1, and results in band bending at the surface and a positive surface potential. Therefore calculations of S
eff are carried out at flat band condition (surface potential ψ
s
=
0) and for a positive surface potential ψ
s
=
0.1 and ψ
s
=
0.2 V. The results in Figure 2(a) represent the case where the correlation energy E
corr is equal to zero, D
i
t is uniform and equal to 1010/cm2/eV throughout the energy gap, and all capture cross-sections are equal to 10−16cm2. It is always assumed in the present analysis that the capture cross-sections of the SRH centers and of neutral amphoteric centers are equal such that σ
n
=
σ
n
0 and σ
p
=
σ
p
0. The impact of asymmetry in other capture cross-sections (σ
n
0
/
σ
p
0
≠
1, σ
n
+
/
σ
n
0
≠
1, σ
p
-
/
σ
p
0
≠
1) and of positive correlation energy will be discussed in later sections.(a) SRH and amphoteric center surface recombination velocityS
SRH and S
amph as a function of excess electron concentration Δ
n
(
W
) assuming all capture cross-sections are equal, E
corr
=
0, for ψ
s
=
0, 0.1 V and 0.2 V. (b) Amphoteric center surface recombination velocity S
amph as a function of excess electron concentration Δ
n
(
W
) assuming all capture cross-sections are equal for E
corr
=
0 and 0.4 eV, for ψ
s
=
0, 0.1 V and 0.2 V.
(a)
(b)By inspecting Figure2(a) it appears that S
SRH and S
amph follow the same trend versus excess electron concentration at the surface, Δ
n
(
W
). The value of S
SRH is larger than S
amph but stays in the same order of magnitude. Both values get closer as the surface potential is increased and S
amph maybecome larger at higher surface potential or higher correlation energy as will be seen in later sections. For several orders of magnitude in the lower Δ
n range S
eff exhibits its highest value that stays constant, then starts to decay gradually in the higher range, and continues to decay to reach very small values close to unity at high values of Δ
n. Such a behavior can be predicted from (3c) since, for a constant value of the surface potential, for several orders of magnitude in the low Δ
n
(
W
) range the n
s term present in the denominator is negligible compared to the p
s term which results in a constant and high value of S
SRH. As Δ
n
(
W
) increases n
s may become significantly boosted such that it starts to affect the value of the denominator of (3c) which results in a gradual decrease in the value of S
SRH as depicted in Figure 2(a). At positive surface potential 0.1 and 0.2 V, which are typical values at the p-type surface due to the presence of positive fixed oxide charges in the SiO2 passivating layer, the value of n
s is exponentially boosted and that of p
s exponentially reduced as predicted by (1a) and (1b) which shifts the onset of S
eff degradation to smaller values of Δ
n
. Moreover, due to the reduced value of p
s the value of the denominator becomes smaller which boosts the maximum constant value of S
eff at low Δ
n.


#### 3.1.1. Impact of the Correlation Energy between Amphoteric Centers

A major property of recombination via amphoteric centers is that the recombination centers are energy correlated. The correlation energy affects the occupation probability and hence would certainly affect the recombination rate and recombination velocity. In the previous section the correlation energy was assumed to be zero which would logically mean a reduced recombination activity. As depicted in Figure2(a), the surface recombination velocity via amphoteric centers indeed increases with increasing correlation energy. For a correlation energy E
corr
=
0.4 eV close to the maximum value, the surface recombination velocity S
amph displayed in Figure 2(b) is more than doubled compared to its values plotted in Figure 2(a) for E
corr
=
0 and exceeds the value of S
SRH at ψ
s
=
0.1 V and 0.2 V. On the other hand, E
corr does not seem to affect the trend nor the ranges characterizing the dependence of S
amph on Δ
n
(
W
).


#### 3.1.2. Asymmetric Electron and Hole Capture Cross-Sections

It has been reported that the capture cross-section of electrons may be up to 100 times larger than the capture cross-sections of holes [17]. Therefore S
eff is calculated assuming that the electron capture cross-sections of SRH centers and of the neutral amphoteric centers are equal (σ
n
=
σ
n
0) and are 100 times larger than their hole capture cross-sections such that σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
100. Equality is also maintained between capture cross-sections of charged and neutral center (R
n
=
σ
n
+
/
σ
n
0
=
1 and R
p
=
σ
p
-
/
σ
p
0
=
1). The impact of the asymmetry seems to affect only the value of S
eff since both S
SRH and S
amph are boosted by approximately a factor of 100 compared to the results in Figure 2(a) which reflects the ratio between σ
n and σ
p but not the general trend it has versus Δ
n, as depicted in Figure 3.Figure 3
Effective surface recombination velocity as a function of excess electron concentration for both SRH and amphoteric center recombination forE
corr
=
0: (a) at flat band (
ψ
s
=
0
), (b) ψ
s
=
0.1 V, (c) ψ
s
=
0.2 V. The ratio σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
100.


#### 3.1.3. Asymmetric Carrier Capture Cross-Sections of Neutral and Charged Centers

Another fundamental difference between SRH statistics and Sah-Shockley statistics is the possibility in the latter to have different capture cross sections for neutral and charged centers. This would lead to a discrepancy betweenS
SRH and S
amph which was not found when capture cross sections of neutral and charged states were assumed to be equal. Indeed, a significant discrepancy is found between S
SRH and S
amph when R
n and R
p are not equal to one, as depicted in Figure 4(a) obtained for ψ
s
=
0.2 V and E
corr
=
0.4 eV. These results are obtained assuming equal capture cross-sections of electrons and holes of SRH and neutral centers σ
n
0
=
σ
p
0
=
σ
n
=
σ
p. When R
n > 1 (σ
n
+
>
σ
n
0) the constant value of S
amph in the low Δ
n range is significantly enhanced and the range itself is strongly reduced which is the result of a boost in the product σ
n
+
n
s strongly affecting the electron capture rate. On the other hand, S
amph is boosted in the high voltage range when R
p
>
1 (σ
p
-
>
σ
p
0) and the range is extended to higher Δ
n as a result of a boost in the product σ
p
-
p
s strongly affecting the hole capture rate. In the presence of both asymmetries R
n
>
1 and R
p
>
1 (σ
n
+
>
σ
n
0 and σ
p
-
>
σ
p
0) S
amph exhibits a double step behavior resulting from the superposition of the low voltage and high voltage trends displayed in Figure 4(a).(a) Effective surface recombination velocity as a function of excess electron concentration for both SRH and amphoteric center recombination forψ
s
=
0.2 V, E
corr
=
0.4 eV with σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
1 for several combinations of R
n
=
σ
n
+
/
σ
n
0 and R
p
=
σ
n
-
/
σ
p
0. (b) Effective surface recombination velocity as a function of excess electron concentration for both SRH and amphoteric center recombination for ψ
s
=
0.2 V, E
corr
=
0.4 eV with σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
1, R
n
=
R
p
=
10, and R
n
=
R
p
=
100.
(a)
(b)


### 3.2. Effective Surface Recombination Velocity versus Surface Potential

Due to the exponential dependence of the surface carrier concentrationsn
s and p
s on the surface potential ψ
s, the surface recombination activity is also a strong function of ψ
s. Such a dependence is very well established and has been investigated thoroughly in previous work [18–20] based on SRH recombination statistics. To our knowledge this work is the first rigorous analysis that treats surface recombination via amphoteric centers at the Si/SiO2 interface, and hence the dependence of S
amph on the surface potential needs to be investigated and compared to the classical SRH dependence. Typical results for the effective surface recombination velocity S
SRH and S
amph versus surface potential are plotted in Figure 5. A single peak in the SRH recombination rate and SRH surface recombination velocity dependence on ψ
s is predicted by (3b) and (3c), respectively. This peak occurs when n
s
σ
n
=
p
s
σ
p, which, when using (1b) and (1a), would occur at a surface potential given by (5a)
ψ
s
,
peak
_
SRH
=
V
T
2
[
ln
⁡
N
A
Δ
n
(
W
)
+
ln
⁡
σ
p
σ
n
]
.
For T
=
300 K, σ
n
=
σ
p, and the given typical values N
A
=
1.5
×
10
16/cm3 and Δ
n
(
W
)
=
10
10/cm3 and according to (5a), S
SRH would peak at ψ
s
,
peak
=
0.184 V which is confirmed in Figure 5. Further increase in surface potential results in a decrease of S
SRH due to the dominance of n
s in the denominator of (3c).Figure 5
Effective surface recombination velocity versus surface potential for an injection levelΔ
n
(
W
)
=
10
10/cm3, E
corr
=
0.4 eV, assuming σ
n
=
σ
n
0 and σ
p
=
σ
p
0 and different combinations of R
n and R
p.On the other hand, our analysis proves thatS
amph versus surface potential is distinguished from S
SRH by the appearance of two peaks, one to the left and one to the right of the SRH peak, as depicted in Figure 5. The formation of these peaks is strongly dependent on the ratios R
n and R
p since, like SRH, they correlate to the relationship between the terms involving the products n
s
σ
n and p
s
σ
p. In this case, however, due to the presence of neutral, positively charged, and negatively charged centers, amphoteric center recombination involves four products instead of two, namely, p
s
σ
p
0, n
s
σ
n
+, p
s
σ
p
-, and n
s
σ
n
0. Assuming σ
n
0
=
σ
p
0 the competition between p
s
σ
p
0 and n
s
σ
n
+ at large values of R
n
=
σ
n
+
/
σ
n
0 leads to a first peak when p
s
σ
p
0
=
n
s
σ
n
+ occurring at ψ
s
=
ψ
s
,
peak
_
R
n given by
(5b)
ψ
s
,
peak
_
R
n
=
V
T
2
[
ln
⁡
N
A
Δ
n
(
W
)
-
ln
⁡
R
n
]
,
which, for example, for R
n
=
50 results in the low surface potential S
amph peak occurring at ψ
s
,
peak
_
R
n
=
0.133 V which is confirmed in Figure 5. On the other hand, the competition between n
s
σ
n
0 and p
s
σ
p
- at large values of R
p
=
σ
p
-
/
σ
p
0 leads to a first peak when n
s
σ
n
0
=
p
s
σ
p
- occurring at ψ
s
=
ψ
s
,
peak
_
R
p given by
(5c)
ψ
s
,
peak
_
R
p
=
V
T
2
[
ln
⁡
N
A
Δ
n
(
W
)
+
ln
⁡
R
p
]
,which, again for R
p
=
50, results in the high surface potential S
amph peak occurring at ψ
s
,
peak
_
R
p
=
0.234 V which is also confirmed in Figure 5. Note that, when R
n or R
p is equal to one, the magnitude of the corresponding S
amph peak is small and tends to be close to the SRH value, as explained in Section 3.1, and occurs exactly at the same surface potential Ψ
s
_
peak
_
SRH as predicted by (5b) and (5c) when R
=
1. Note that the two peak behavior is pronounced only if both R
n and R
p are large and comparable. On the other hand, if R
n or R
p are asymmetric, the peak associated with the larger is maintained while the peak associated with the smaller disappears or shows as a hump or a shoulder at the peak position determined by (5b) and (5c). It is clear that according to (5b) and (5c), the whole plot in Figure 5 may be shifted to the right or to the left by changing the excess electron concentration Δ
n, and the peak positions may be stretched away from each other or brought closer to each other by modifying the ratios R
n and R
p. The two-peak behavior and the peak-hump behavior have been previously reported for surface recombination velocity deduced from measurements of the recombination current at the Si/SiO2 interface in gated diodes and gated transistors (e.g., [18, 21, 22]). The presence of mobile ions in the oxide layer [21], nonuniformities at the Si/SiO2 interface [18], interface states having different capture cross-sections [22] are previous interpretations for the two-peak behavior which in our opinion still need to be justified.Using the proposed amphoteric center recombination model a successful fitting is demonstrated in Figure6 for an experimental double peak surface recombination—versus n
s
/
p
s plot previously reported in the literature for a p-type surface doping concentration N
A
=
5.5
×
10
16/cm3 and excess electron concentration Δ
n
=
3
×
10
16/cm3 [22]. It is possible to obtain such fitting using several combinations with reasonable values for the parameters D
i
t, σ
n
0, σ
p
0, R
n, and R
p and therefore there is no need to insist on a specific combination. It is impossible, however, to fit these data using SRH recombination unless mathematical oriented assumptions with no physical justifications are proposed such as splitting the centers into groups with different capture cross-sections. On the other hand, the present analysis firmly attributes the two-peak surface recombination velocity to surface recombination via amphoteric centers with larger capture cross-sections of charged centers than of neutral centers.Figure 6
Fitting of experimentally measured surface recombination velocityS
gen
=
S
eff
×
(
Δ
n
/
p
)
0.5 versus surface electron/hole concentration ratio n
s
/
p
s.Circles: experimental results [22],line: results of proposed model.


## 3.1. Effective Surface Recombination Velocity versus Excess Carrier Concentration

The effective surface recombination velocityS
eff is calculated using SRH and amphoteric center recombination models and plotted as a function of the excess electron concentration Δ
n
(
W
) at the edge of the surface field region. Since nonnegligible positive fixed charges are present in the passivating oxide layer, the p-type surface would be depleted which creates a surface field region, as shown in Figure 1, and results in band bending at the surface and a positive surface potential. Therefore calculations of S
eff are carried out at flat band condition (surface potential ψ
s
=
0) and for a positive surface potential ψ
s
=
0.1 and ψ
s
=
0.2 V. The results in Figure 2(a) represent the case where the correlation energy E
corr is equal to zero, D
i
t is uniform and equal to 1010/cm2/eV throughout the energy gap, and all capture cross-sections are equal to 10−16cm2. It is always assumed in the present analysis that the capture cross-sections of the SRH centers and of neutral amphoteric centers are equal such that σ
n
=
σ
n
0 and σ
p
=
σ
p
0. The impact of asymmetry in other capture cross-sections (σ
n
0
/
σ
p
0
≠
1, σ
n
+
/
σ
n
0
≠
1, σ
p
-
/
σ
p
0
≠
1) and of positive correlation energy will be discussed in later sections.(a) SRH and amphoteric center surface recombination velocityS
SRH and S
amph as a function of excess electron concentration Δ
n
(
W
) assuming all capture cross-sections are equal, E
corr
=
0, for ψ
s
=
0, 0.1 V and 0.2 V. (b) Amphoteric center surface recombination velocity S
amph as a function of excess electron concentration Δ
n
(
W
) assuming all capture cross-sections are equal for E
corr
=
0 and 0.4 eV, for ψ
s
=
0, 0.1 V and 0.2 V.
(a)
(b)By inspecting Figure2(a) it appears that S
SRH and S
amph follow the same trend versus excess electron concentration at the surface, Δ
n
(
W
). The value of S
SRH is larger than S
amph but stays in the same order of magnitude. Both values get closer as the surface potential is increased and S
amph maybecome larger at higher surface potential or higher correlation energy as will be seen in later sections. For several orders of magnitude in the lower Δ
n range S
eff exhibits its highest value that stays constant, then starts to decay gradually in the higher range, and continues to decay to reach very small values close to unity at high values of Δ
n. Such a behavior can be predicted from (3c) since, for a constant value of the surface potential, for several orders of magnitude in the low Δ
n
(
W
) range the n
s term present in the denominator is negligible compared to the p
s term which results in a constant and high value of S
SRH. As Δ
n
(
W
) increases n
s may become significantly boosted such that it starts to affect the value of the denominator of (3c) which results in a gradual decrease in the value of S
SRH as depicted in Figure 2(a). At positive surface potential 0.1 and 0.2 V, which are typical values at the p-type surface due to the presence of positive fixed oxide charges in the SiO2 passivating layer, the value of n
s is exponentially boosted and that of p
s exponentially reduced as predicted by (1a) and (1b) which shifts the onset of S
eff degradation to smaller values of Δ
n
. Moreover, due to the reduced value of p
s the value of the denominator becomes smaller which boosts the maximum constant value of S
eff at low Δ
n.


### 3.1.1. Impact of the Correlation Energy between Amphoteric Centers

A major property of recombination via amphoteric centers is that the recombination centers are energy correlated. The correlation energy affects the occupation probability and hence would certainly affect the recombination rate and recombination velocity. In the previous section the correlation energy was assumed to be zero which would logically mean a reduced recombination activity. As depicted in Figure2(a), the surface recombination velocity via amphoteric centers indeed increases with increasing correlation energy. For a correlation energy E
corr
=
0.4 eV close to the maximum value, the surface recombination velocity S
amph displayed in Figure 2(b) is more than doubled compared to its values plotted in Figure 2(a) for E
corr
=
0 and exceeds the value of S
SRH at ψ
s
=
0.1 V and 0.2 V. On the other hand, E
corr does not seem to affect the trend nor the ranges characterizing the dependence of S
amph on Δ
n
(
W
).


### 3.1.2. Asymmetric Electron and Hole Capture Cross-Sections

It has been reported that the capture cross-section of electrons may be up to 100 times larger than the capture cross-sections of holes [17]. Therefore S
eff is calculated assuming that the electron capture cross-sections of SRH centers and of the neutral amphoteric centers are equal (σ
n
=
σ
n
0) and are 100 times larger than their hole capture cross-sections such that σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
100. Equality is also maintained between capture cross-sections of charged and neutral center (R
n
=
σ
n
+
/
σ
n
0
=
1 and R
p
=
σ
p
-
/
σ
p
0
=
1). The impact of the asymmetry seems to affect only the value of S
eff since both S
SRH and S
amph are boosted by approximately a factor of 100 compared to the results in Figure 2(a) which reflects the ratio between σ
n and σ
p but not the general trend it has versus Δ
n, as depicted in Figure 3.Figure 3
Effective surface recombination velocity as a function of excess electron concentration for both SRH and amphoteric center recombination forE
corr
=
0: (a) at flat band (
ψ
s
=
0
), (b) ψ
s
=
0.1 V, (c) ψ
s
=
0.2 V. The ratio σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
100.


### 3.1.3. Asymmetric Carrier Capture Cross-Sections of Neutral and Charged Centers

Another fundamental difference between SRH statistics and Sah-Shockley statistics is the possibility in the latter to have different capture cross sections for neutral and charged centers. This would lead to a discrepancy betweenS
SRH and S
amph which was not found when capture cross sections of neutral and charged states were assumed to be equal. Indeed, a significant discrepancy is found between S
SRH and S
amph when R
n and R
p are not equal to one, as depicted in Figure 4(a) obtained for ψ
s
=
0.2 V and E
corr
=
0.4 eV. These results are obtained assuming equal capture cross-sections of electrons and holes of SRH and neutral centers σ
n
0
=
σ
p
0
=
σ
n
=
σ
p. When R
n > 1 (σ
n
+
>
σ
n
0) the constant value of S
amph in the low Δ
n range is significantly enhanced and the range itself is strongly reduced which is the result of a boost in the product σ
n
+
n
s strongly affecting the electron capture rate. On the other hand, S
amph is boosted in the high voltage range when R
p
>
1 (σ
p
-
>
σ
p
0) and the range is extended to higher Δ
n as a result of a boost in the product σ
p
-
p
s strongly affecting the hole capture rate. In the presence of both asymmetries R
n
>
1 and R
p
>
1 (σ
n
+
>
σ
n
0 and σ
p
-
>
σ
p
0) S
amph exhibits a double step behavior resulting from the superposition of the low voltage and high voltage trends displayed in Figure 4(a).(a) Effective surface recombination velocity as a function of excess electron concentration for both SRH and amphoteric center recombination forψ
s
=
0.2 V, E
corr
=
0.4 eV with σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
1 for several combinations of R
n
=
σ
n
+
/
σ
n
0 and R
p
=
σ
n
-
/
σ
p
0. (b) Effective surface recombination velocity as a function of excess electron concentration for both SRH and amphoteric center recombination for ψ
s
=
0.2 V, E
corr
=
0.4 eV with σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
1, R
n
=
R
p
=
10, and R
n
=
R
p
=
100.
(a)
(b)


## 3.1.1. Impact of the Correlation Energy between Amphoteric Centers

A major property of recombination via amphoteric centers is that the recombination centers are energy correlated. The correlation energy affects the occupation probability and hence would certainly affect the recombination rate and recombination velocity. In the previous section the correlation energy was assumed to be zero which would logically mean a reduced recombination activity. As depicted in Figure2(a), the surface recombination velocity via amphoteric centers indeed increases with increasing correlation energy. For a correlation energy E
corr
=
0.4 eV close to the maximum value, the surface recombination velocity S
amph displayed in Figure 2(b) is more than doubled compared to its values plotted in Figure 2(a) for E
corr
=
0 and exceeds the value of S
SRH at ψ
s
=
0.1 V and 0.2 V. On the other hand, E
corr does not seem to affect the trend nor the ranges characterizing the dependence of S
amph on Δ
n
(
W
).


## 3.1.2. Asymmetric Electron and Hole Capture Cross-Sections

It has been reported that the capture cross-section of electrons may be up to 100 times larger than the capture cross-sections of holes [17]. Therefore S
eff is calculated assuming that the electron capture cross-sections of SRH centers and of the neutral amphoteric centers are equal (σ
n
=
σ
n
0) and are 100 times larger than their hole capture cross-sections such that σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
100. Equality is also maintained between capture cross-sections of charged and neutral center (R
n
=
σ
n
+
/
σ
n
0
=
1 and R
p
=
σ
p
-
/
σ
p
0
=
1). The impact of the asymmetry seems to affect only the value of S
eff since both S
SRH and S
amph are boosted by approximately a factor of 100 compared to the results in Figure 2(a) which reflects the ratio between σ
n and σ
p but not the general trend it has versus Δ
n, as depicted in Figure 3.Figure 3
Effective surface recombination velocity as a function of excess electron concentration for both SRH and amphoteric center recombination forE
corr
=
0: (a) at flat band (
ψ
s
=
0
), (b) ψ
s
=
0.1 V, (c) ψ
s
=
0.2 V. The ratio σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
100.


## 3.1.3. Asymmetric Carrier Capture Cross-Sections of Neutral and Charged Centers

Another fundamental difference between SRH statistics and Sah-Shockley statistics is the possibility in the latter to have different capture cross sections for neutral and charged centers. This would lead to a discrepancy betweenS
SRH and S
amph which was not found when capture cross sections of neutral and charged states were assumed to be equal. Indeed, a significant discrepancy is found between S
SRH and S
amph when R
n and R
p are not equal to one, as depicted in Figure 4(a) obtained for ψ
s
=
0.2 V and E
corr
=
0.4 eV. These results are obtained assuming equal capture cross-sections of electrons and holes of SRH and neutral centers σ
n
0
=
σ
p
0
=
σ
n
=
σ
p. When R
n > 1 (σ
n
+
>
σ
n
0) the constant value of S
amph in the low Δ
n range is significantly enhanced and the range itself is strongly reduced which is the result of a boost in the product σ
n
+
n
s strongly affecting the electron capture rate. On the other hand, S
amph is boosted in the high voltage range when R
p
>
1 (σ
p
-
>
σ
p
0) and the range is extended to higher Δ
n as a result of a boost in the product σ
p
-
p
s strongly affecting the hole capture rate. In the presence of both asymmetries R
n
>
1 and R
p
>
1 (σ
n
+
>
σ
n
0 and σ
p
-
>
σ
p
0) S
amph exhibits a double step behavior resulting from the superposition of the low voltage and high voltage trends displayed in Figure 4(a).(a) Effective surface recombination velocity as a function of excess electron concentration for both SRH and amphoteric center recombination forψ
s
=
0.2 V, E
corr
=
0.4 eV with σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
1 for several combinations of R
n
=
σ
n
+
/
σ
n
0 and R
p
=
σ
n
-
/
σ
p
0. (b) Effective surface recombination velocity as a function of excess electron concentration for both SRH and amphoteric center recombination for ψ
s
=
0.2 V, E
corr
=
0.4 eV with σ
n
/
σ
p
=
σ
n
0
/
σ
p
0
=
1, R
n
=
R
p
=
10, and R
n
=
R
p
=
100.
(a)
(b)


## 3.2. Effective Surface Recombination Velocity versus Surface Potential

Due to the exponential dependence of the surface carrier concentrationsn
s and p
s on the surface potential ψ
s, the surface recombination activity is also a strong function of ψ
s. Such a dependence is very well established and has been investigated thoroughly in previous work [18–20] based on SRH recombination statistics. To our knowledge this work is the first rigorous analysis that treats surface recombination via amphoteric centers at the Si/SiO2 interface, and hence the dependence of S
amph on the surface potential needs to be investigated and compared to the classical SRH dependence. Typical results for the effective surface recombination velocity S
SRH and S
amph versus surface potential are plotted in Figure 5. A single peak in the SRH recombination rate and SRH surface recombination velocity dependence on ψ
s is predicted by (3b) and (3c), respectively. This peak occurs when n
s
σ
n
=
p
s
σ
p, which, when using (1b) and (1a), would occur at a surface potential given by (5a)
ψ
s
,
peak
_
SRH
=
V
T
2
[
ln
⁡
N
A
Δ
n
(
W
)
+
ln
⁡
σ
p
σ
n
]
.
For T
=
300 K, σ
n
=
σ
p, and the given typical values N
A
=
1.5
×
10
16/cm3 and Δ
n
(
W
)
=
10
10/cm3 and according to (5a), S
SRH would peak at ψ
s
,
peak
=
0.184 V which is confirmed in Figure 5. Further increase in surface potential results in a decrease of S
SRH due to the dominance of n
s in the denominator of (3c).Figure 5
Effective surface recombination velocity versus surface potential for an injection levelΔ
n
(
W
)
=
10
10/cm3, E
corr
=
0.4 eV, assuming σ
n
=
σ
n
0 and σ
p
=
σ
p
0 and different combinations of R
n and R
p.On the other hand, our analysis proves thatS
amph versus surface potential is distinguished from S
SRH by the appearance of two peaks, one to the left and one to the right of the SRH peak, as depicted in Figure 5. The formation of these peaks is strongly dependent on the ratios R
n and R
p since, like SRH, they correlate to the relationship between the terms involving the products n
s
σ
n and p
s
σ
p. In this case, however, due to the presence of neutral, positively charged, and negatively charged centers, amphoteric center recombination involves four products instead of two, namely, p
s
σ
p
0, n
s
σ
n
+, p
s
σ
p
-, and n
s
σ
n
0. Assuming σ
n
0
=
σ
p
0 the competition between p
s
σ
p
0 and n
s
σ
n
+ at large values of R
n
=
σ
n
+
/
σ
n
0 leads to a first peak when p
s
σ
p
0
=
n
s
σ
n
+ occurring at ψ
s
=
ψ
s
,
peak
_
R
n given by
(5b)
ψ
s
,
peak
_
R
n
=
V
T
2
[
ln
⁡
N
A
Δ
n
(
W
)
-
ln
⁡
R
n
]
,
which, for example, for R
n
=
50 results in the low surface potential S
amph peak occurring at ψ
s
,
peak
_
R
n
=
0.133 V which is confirmed in Figure 5. On the other hand, the competition between n
s
σ
n
0 and p
s
σ
p
- at large values of R
p
=
σ
p
-
/
σ
p
0 leads to a first peak when n
s
σ
n
0
=
p
s
σ
p
- occurring at ψ
s
=
ψ
s
,
peak
_
R
p given by
(5c)
ψ
s
,
peak
_
R
p
=
V
T
2
[
ln
⁡
N
A
Δ
n
(
W
)
+
ln
⁡
R
p
]
,which, again for R
p
=
50, results in the high surface potential S
amph peak occurring at ψ
s
,
peak
_
R
p
=
0.234 V which is also confirmed in Figure 5. Note that, when R
n or R
p is equal to one, the magnitude of the corresponding S
amph peak is small and tends to be close to the SRH value, as explained in Section 3.1, and occurs exactly at the same surface potential Ψ
s
_
peak
_
SRH as predicted by (5b) and (5c) when R
=
1. Note that the two peak behavior is pronounced only if both R
n and R
p are large and comparable. On the other hand, if R
n or R
p are asymmetric, the peak associated with the larger is maintained while the peak associated with the smaller disappears or shows as a hump or a shoulder at the peak position determined by (5b) and (5c). It is clear that according to (5b) and (5c), the whole plot in Figure 5 may be shifted to the right or to the left by changing the excess electron concentration Δ
n, and the peak positions may be stretched away from each other or brought closer to each other by modifying the ratios R
n and R
p. The two-peak behavior and the peak-hump behavior have been previously reported for surface recombination velocity deduced from measurements of the recombination current at the Si/SiO2 interface in gated diodes and gated transistors (e.g., [18, 21, 22]). The presence of mobile ions in the oxide layer [21], nonuniformities at the Si/SiO2 interface [18], interface states having different capture cross-sections [22] are previous interpretations for the two-peak behavior which in our opinion still need to be justified.Using the proposed amphoteric center recombination model a successful fitting is demonstrated in Figure6 for an experimental double peak surface recombination—versus n
s
/
p
s plot previously reported in the literature for a p-type surface doping concentration N
A
=
5.5
×
10
16/cm3 and excess electron concentration Δ
n
=
3
×
10
16/cm3 [22]. It is possible to obtain such fitting using several combinations with reasonable values for the parameters D
i
t, σ
n
0, σ
p
0, R
n, and R
p and therefore there is no need to insist on a specific combination. It is impossible, however, to fit these data using SRH recombination unless mathematical oriented assumptions with no physical justifications are proposed such as splitting the centers into groups with different capture cross-sections. On the other hand, the present analysis firmly attributes the two-peak surface recombination velocity to surface recombination via amphoteric centers with larger capture cross-sections of charged centers than of neutral centers.Figure 6
Fitting of experimentally measured surface recombination velocityS
gen
=
S
eff
×
(
Δ
n
/
p
)
0.5 versus surface electron/hole concentration ratio n
s
/
p
s.Circles: experimental results [22],line: results of proposed model.


## 4. Conclusions

A model is proposed for surface recombination via amphoteric defects at the Si/SiO2 interface of thermally oxidized p-type silicon surface. The model is an adaptation to the surface of the model developed for the bulk recombination in amorphous silicon based on Sah-Shockley multicharge statistics for energy correlated amphoteric dangling bond defects. The results indicate that surface recombination via amphoteric centers behaves in a very similar way to SRH recombination when capture cross-sections of charged and neutral centers are equal or close to equal. On the other hand, if electron or hole capture cross-sections are larger for charged centers than for neutral centers the surface recombination velocity is enhanced at small or at high excess electron concentration. If both asymmetries are present, the superposition of both trends leads to a double step behavior of S
amph versus excess electron concentration. In addition, in the presence of such asymmetries the surface recombination velocity exhibits two peaks versus surface potential, which relate to a larger electron capture cross-section of positively charged centers and to a larger hole capture cross-section for negatively charged centers than for neutral centers. Expressions are derived for the surface potentials at which the two peaks occur. Recombination via amphoteric centers with asymmetric capture cross-sections for charged and neutral centers is a firm and plausible interpretation for the double peak surface recombination velocity extracted from measurements of surface recombination current in gated diode and gated transistor versus gate voltage.



---

*Source: 857907-2014-03-12.xml*