## Revista Brasileira de Ensino de Física

*On-line version* ISSN 1806-9126

### Rev. Bras. Ensino Fís. vol.28 no.2 São Paulo Apr./June 2006

#### http://dx.doi.org/10.1590/S0102-47442006000200010

**ARTIGOS GERAIS**

**Temperature dependence of impurity levels in bulk semiconductors**

**Estudo do efeito térmico sobre os níveis de impurezas em semicondutores**

**M.A. Amato ^{1} **

Instituto de Física, Universidade de Brasília, Brasília, DF, Brasil

**ABSTRACT**

The temperature dependence of impurity levels in semiconductors is a very interesting problem in the domain of deep levels impurities. Its evidence in experimental results is reported by many authors. However, the interpretation on theoretical basis requires sophisticated models. In this paper we propose a very simple model that allows to extract the relevant parameters from the experiment without needing heavy computation.

**Keywords:** bulk semiconductors, impurity levels, thermal properties.

**RESUMO**

O estudo da dependência da posição dos níveis de energia das impurezas profundas em semicondutores com a temperatura é um problema bastante interessante, e sua evidência experimental é discutida por vários autores. Neste artigo propomos um modelo teórico bastante simples, a partir do conhecimento da secção de choque de fotoionização e de alguns parametros associados aos fónons, para o entendimento deste comportamento.

**Palavras-chave:** semicondutores volumétricos, impurezas profundas, propriedades térmicas.

In recent years progress has been made in investigating the physical properties of impurity centres in semiconductors either theoretically or experimentally. For deep level impurities their identification and characterisation remains as a difficult problem and it is still an important and active field in semiconductor research [1]. On the other hand, many of the shallow impurities are by now well understood. The effective mass theory of shallow impurities [2] led to the well-known hydrogenic theory providing both conceptual guidance and, in many cases, quantitative predictions about the binding energies [3]

However, one major difficult to obtain accurate solutions to the problem is the presence of the interaction between the electrons and the lattice. The model generally used for a quantitative treatment of electron-phonon interaction is the so-called configuration co-ordinate model. Its formalism has developed a long time ago and reviews exist on this topic for semiclassical as well as quantum treatment [4, 5, 6]. In its original form this model is worked out for transitions within localised states, and it has been applied successfully to such systems as F-centres in alkali halides and rare-earth impurities in semiconductors. It is also widely accepted for studying transitions related to deep centres, mainly due to a large Franck-Condon shift [7].

It is generally accepted that even were a full theory available, it would be necessary to extract from it an acceptable, simplified model which could be used to define the principal measurable quantities and to provide a useful terminology for describing experimental results. The presence of such a model is of some importance and we attempt in this article to illustrate what experimental features lend themselves most readily to modelling.

Whilst many of these subjects can readily be discussed in a qualitative way in an undergraduate class, it is difficult to go into the theory underlying the thermal process unless the students have a sound background in solid state theory. However, some important quantitative features of the process can be conveyed by elementary methods described in this paper.

The temperature dependence of the energy levels in semiconductors is discussed in Ref. [8] for different types of defects in Si and Ge, and in Refs. [9] and [10] for some deep levels in compound semiconductors. In particular Ref. [11] provides a detailed discussion of thermodinamics functions involved in the thermal processes associated impurities in semiconductors. A rather common technique applied to interpret the experimental results is the convolution method. An exaustive discussion of this method can be find in Ref. [12] and references therein.

In order to proceed for the study of the temperature dependence of the energy levels, we just make an attempt to look and understand these effects in terms of thermal broadening. We also admit that the primary effect of the electron-phonon coupling is to shift the edge of the photoionisation cross section [13]. Thus, we may write

where *E _{I}* is the ionisation energy,

*S*a measure of the electron-phonon coupling (the Huang-Rhys factor), w

*the phonon energy that takes part in the electronic transition, which we admit a transition from a donor-like impurity to the conduction band, and is the the optical ionisation energy in the absence of electron-phonon coupling (*

_{p}*S*= 0) . More realistically, the coupling with the lattice introduces a temperature dependence of photoionization cross section and a shift in the threshold energy. An interesting discussion involving indirect transitions is given in Ref. [14].

To model this effect most simply we will assume that the lattice coupling does not strongly affects the shape of the spectral dependence, but merely shift the threshold energy by an amount proportional to the temperature [15]

where b* _{T}* is the temperature coefficient of the apparent shift of energy with temperature caused by thermal broadening, and

*E*(0) is the threshold energy

_{I} Essentially, what we have done is to associate thermal broadening with the apparent shift of threshold energy with temperature. Its origin is in the temperature dependence of emitted phonons, *p*, which maximizes the transition rate *W _{p}* [13]

where |*M*| ^{2} is the electronic matrix element, and r(*E*_{k}) is the electronic density of states, *J*(*E*_{p}) is the square of the matrix element involving the overlap of displaced oscillator wavefunctions [16] for the case of emitting an energy *E _{p}* in phonons. It also accomodates the case of absorption, in this case,

*E*being negative.

_{p}*J*(

*E*) maximizes at all temperatures when

_{p}*p = S*. We illustrate this behaviour in Fig. 1, for

*S*= 3.0 and an effective temperature

*T*

^{*}= 1.0 in units of w

*.*

_{p}/k_{B}

However, due to the dependence of *W _{p}* and r(

*E*

_{k}) on the final states with energy

*E*

_{k}the maximum of

*W*is not at

_{p}*p = S*but at a value of

*p*less than

*S*, since for phonon emission process the transition rate increases with

*E*

_{k}, and

*E*

_{k}=

*h*n-

*E*-

_{I}*pw*, with hn being the energy of the incident radiation [17]. Moreover, this value of

_{p}*p*decreases with increasing temperature as the dependence of

*J*(

*E*) on

_{p}*p*become less pronounced, and it is this temperature-dependence which causes the apparent threshold to vary.

As an illustration, this behaviour of the photoionisation cross section [18] for a donor type impurity with photon energy at different temperatures is shown in Fig. 2, for a trap state with ionisation energy = 20*w _{p}* [19].

Although Eq. (2) is an empirical equation, we have to provide a way to obtain b* _{T}*. This temperature coefficient is derived graphically from plots similar to those shown in Fig. 4 for various values of

*S*,

*w*and . The procedure for obtaining this coefficient is depicted in the following paragraphs.

_{p}

A plot of [ s(*h*n)·*h*n] ^{2/3} *vs. h*n for three different values of temperature in units of *w _{p}/k_{B}* is depicted in Fig. 3. It shows an apparent change in the threshold of the ionisation energy, and as we get near this threshold the curve bends upwards due tho the presence of phonons in the transition. A linear fitting of this curve provides the threshold ionisation energy due to photon absorption, or the optical ionisation energy as a function of temperature

*E*(

_{I}*T*). So, for different temperatures one gets different values of

*E*(

_{I}*T*). However, the photoionisation cross section is sensitive to different values of S, and so should be this threshold energy. This behaviour is shown in Fig. 4.

A linear fit of these curves provides b_{T} and *E _{I}*( 0)

Yet the coefficient b_{T} in the above equation depends on , *S* and *w _{p}*. For electrons making a transition between |sñ -like states only the zone edge LA and LO phonon energies are involved [20], so

*w*is in principle known. However, neither

_{p}*S*or is obtained separately from this observation of effective threshold. To circumvent this problem one has to obtain a relation between (b

_{T}) and . Figure 5 provides a way to obtain such relation. From the graph we obtain:

with b_{0} __~__ 3.2*k _{B}*, where

*k*is the Boltzmann constant and b

_{B}_{T}is proportional to

*S*. As expected, inserting this into Eq. (5) one recovers Eq. (3).

Measured values of cross sections, or more usually, of a quantity proportional to s(*h*n) such as photosensitivity - the change in photoconductivity per photon, which are used to obtain the effective threshold in the way described in this paper, span typically an order of magnitude near the shoulder of the spectral-dependence curve. Below this region lies a spectral tail which extends through and below the effective threshold, spanning two or more orders of magnitude. The form of this tail is sensitive to temperature and to lattice-coupling strength, and provided it is not influenced by lower energy optical transitions, its form and temperature dependence can be analysed to provide a measure of S. One may wonders about the effect of the band-gap energy as it also depends on temperature on the determination of *E*(*T*). Rather than going into the calculations in order to obtain the relevant thermodynamics parameter, such as the Gibbs free energy and thus the entropy, we suggest that it may additive, so we should include into Eq. (2) a term similar to one encountered in Varshini equation [21].

In this paper one tries to provide a simple method which could be easily applied to experimental results in order to extract relevant parameters and replace heavy computation, along the same lines we have proposed in Ref. [22]. In an experiment one does not know the charge of the centre, and one does not need to know the charge to obtain something useful, so its effect is not include in this simple model. Also, the size of the centre is not included. As a step forward one proposes an empirical formula for temperature dependence of the energy threshold associated with optical transitions involving impurities in semiconductors. In summary, the model described in this paper is undoubtedly over-simplistic, and should not be regarded as anything more than a conceptual aid. On the other side, it does provide an easily understood introduction to the subject, and may enable some useful quantitative results to be obtained without a detailed quantum mechanical calculations.

**Referências**

[1] S. Pantelides, in *Deep Centers in Semiconductors: A State of the Art Approach*, S. Pantelides (ed.), Perspectives in the Past, Present, and Future of Deep Centers (Gordon and Breach, New York, 1997), p. 1-85. [ Links ]

[2] J.M. Luttinger and W. Kohn. Phys. Rev. **97**, 969 (1955). [ Links ]

[3] W. Kohn, in *Solid State Physics*, F. Seitz and D. Turnbull (eds.), *Shallow Impurity States in Silicon and Germanium, 5* (Academic Press, New York, 1957), p. 257-320. [ Links ]

[4] P. Vogl, H.P. Hjalmarson and J.D. Dow, J. Phys. Chem. Solids **44**, 365 (1983). [ Links ]

[5] K. Huang and A. Rhys, Proc. Royal Soc. A **204**, 406 (1950). [ Links ]

[6] M. Lax, J. Chem. Phys. **20**, 1752 (1952). [ Links ]

[7] J.J. Markham, Rev. Mod. Phys. **31**, 956 (1959). [ Links ]

[8] J.A. Van Vechten and C.D Thurmond, Phys. Rev. B **14**, 3539 (1976). [ Links ]

[9] R. Pässler, J. Appl. Phys. **96**, 715 (2004); id. **97**(11) 3533 (2005). [ Links ]

[10] L. Samuelson and B. Monemar, Phys. Rev. B **18**, 830 (1978). [ Links ]

[11] C.D. Thurmond, J. Eletrochem. Soc. **122**, 1133 (1975). [ Links ]

[12] R. Pässler, H. Petterson, H.G. Grimmeis and K.R. Schulz, Phys. Rev. B **55**, 4312 (1997). [ Links ]

[13] B.K. Ridley and M.A. Amato, J. Phys. C: Solid State Phys. **14**, 1255 (1981). [ Links ]

[14] H. Dai, M.A. Gundersen, C.W. Myles and P.G. Snyder, Phys. Rev. B **37**, 1205 (1988). [ Links ]

[15] Some authors refer to *E*(*T)* as zero-phonon binding energy.

[16] See footnote 16

[17] See footnote 17

[18] M.A. Amato and B.K. Ridley, J. Phys. C: Solid State Phys. **13**, 2027 (1980). [ Links ]

[19] The theoretical model for obtaining s(*h*n) is given in Ref [18]. It includes the charge of the centre *Z* by inserting in the conduction band wavefunction the Coulomb factor. However, in this paper one just consider the case *Z* = 0.

[20] M.A. Amato and B.K. Ridley, Phys. Letter A **78**, 170 (1980). [ Links ]

[21] Y.P. Varshini, Physica **34**, 149 (1967). [ Links ] The Varshini relation is quoted in Jacques I. Pankove, in *Optical Processes in Semiconductors* (Dover, New York, 1975), p. 27. [ Links ]

[22] M.A. Amato, Rev. Bras. Ens. Fís. **24**, 1 (2002). [ Links ]

Recebido em 24/10/2005; Aceito em 24/2/2006

1 E-mail: maamato@unb.br.

[16]* J*(*E _{p}*) is the square of the matrix element involving the overlap of displaced oscillator wavefunctions for the case of emitting an energy

*E*in phonons. A quantum mechanical calculation gives

_{p}*J*(

*E*) = exp[-2

_{p}*S*(

*n*+

^{})]exp

*I*{2

_{p}*S*[

*n*(

*n*+1)]

^{1/2}}, where

*n*is the Bose-Einstein phonon population, and

*I*(

_{p}*z*) is the modified Bessel function.

[17] For phonon absorption processes one has

*E*

_{k}=

*h*n-

*E*+

_{I}*p*

*w*.

_{p}