Band-gap
widening
in
heavily
Sn-daped
In203
I.
Hamberg andC.
G.
GranqvistPhysics Department, Chalmers University
of
Technology, S-412 96Goteborg, SwedenK.
-F.
Berggren,B.E.
Sernelius, andL.
EngstromTheoretical Physics Group, Department
of
Physics and Measurement Technology,Linkoping University, S-581 83Linkoping, Sweden
(Received 17 April 1984)
Films ofpure and Sn-doped semiconducting Inz03 were prepared by reactive e-beam evaporation. The spectral absorption coefficient was evaluated by spectrophotometry in the (2
—
6)-eV range. The extracted band gap increases with electron density (n,) approximately as n, for n,(10
' cm This result isinterpreted within an effective-mass model for n-doped semiconductors well above the Mott critical density. Because ofthe high degree ofdoping, the impurities are ionized and the asso-ciated electrons occupy the bottom ofthe conduction band in the form ofan electron gas. Themodel accounts for a Burstein-Moss shift aswell as electron-electron and electron-impurity
scatter-ing treated in the random-phase approximation. Experiments and theory were reconciled by assum-ing a parabolic valence band with an effective mass
-0.
6m. Earlier work on doped oxide semicon-ductors are assessed inthe light ofthe present results.I.
INTRODUCTION AND SUMMARYIn this paper we investigate the optical properties
of
evaporated filmsof
doped semiconducting In203 in the (2—
6)-eV range,i.
e., around the fundamental band gap.The study serves two main purposes: to elucidate basic properties
of
aheavily n-doped semiconductor, and to im-prove our understandingof
a technologically important material which is widely used when transmittanceof
visi-ble or solar radiation needs to be combined with good electrical conduction orlow thermal emittance.We are interested in semiconductors which are n doped sothat the Mott critical density is exceeded and for which electrons occupy the host conduction band in the form
of
an electron gas. The energy gap is shifted as a resultof
the doping. The magnitudeof
the shift is determined by two competing mechanisms. There is aband-gap nar~om-ing which is a consequence'of
many-body effects on the conduction and valence bands. This shrinkage is coun-teracted by the Burstein-Moss effect3 which gives a band-gap widening as a resultof
the blockingof
the lowest states in the conduction band.For
doped In203, the net effect istoincrease the gap.We report below in Sec.
II
on the production and analysisof
evaporated In203 films with up to 9 mol%%uoSn02. The spectral absorption coefficient is evaluated for films with electron densities up to
—
10'cm;
the data in-dicate a band-gap widening by as much as-0.
8 eV.To
understand this widening from basic principles, we outline in Sec.III
a theory which includes the Burstein-Moss shift and self-energies due to electron-electron and electron-impurity scattering. The calculations are per-formed within the frameworkof
the random-phase ap-proximation (RPA), along the lines given in an earlier pa-per.A proper comparison
of
theory and experiments must account for the shapeof
the spectral absorption coeffi-cient around the energy gap. The analysis procedure is presented in Sec.IV.
We are able to reconcile theory and experiments by having effective electron massesof
-0.
3m for the conduction band and-0.
6m for the valence band (where m is the free-electron mass). Our techniques for evaluating and interpreting the band gaps go beyond what is normally done for oxide semiconduct-ors. Hence itisof
interest to consider the earlier work in the lightof
the present results. This is done for In203, CdO, and Sn02 in Sec. V. As far as we know, all earlier work on these materials have neglected the self-energies.It
is demonstrated that this may lead to aqualitatively in-correct interpretationof
the shifted band gaps.The present analysis provides aconsistent model forthe optical properties
of
doped In203 around its fundamental band gap.It
embraces a dopingof
a host semiconductor to achieve a high densityof
electrons described within theRPA
—
and accounts for the ensuing ionized impuri-ties.It
is gratifying that the same model can be used 'for the optical performance in the infrared, where the properties are governed by a degenerate electron gas with a plasma energy
—
1 eVfor the highest doping levels. We are thus establishing a model capableof
explaining the key propertiesof
transparent conductors and transparent heat mirrors: a high transmittance between a properly shifted and broadened energy gap and the plasma wave-length, and a high reflectance and concomitant high con-ductivity beyond the plasma wavelength. Coatings with these properties, applied to substratesof
glass and plastic, are extensively used for alarge numberof
applications; someof
the most important are front-surface electrodes on solar cells and display devices, and low-emittance coat-ings forwindows.II.
FILM FABRICATION AND OPTICAL DATA The thin films studied in this work were produced by reactive e-beam evaporationof
pure In203 andof
In203 with up to 9mo1%
Sn02 onto substratesof
CaFz in a sys-tem with accurate process controls. The deposition ma-terials were hot-pressed pellets with99.
99%
purity sup-plied by Kyodo International, Japan. The evaporation rate was kept at a constant value in the0.
2—
0.
3nm/s in-terval by using the output from a vibrating quartz micro-balance to feedback control the electrical powerof
the e-beam source. A constant oxygen pressure in the (5—
8)X10
-Torr range was maintained by continuous gas inlet through a precision valve. The substrate, posi-tioned 3Scm above the source, was kept at-300'C
dur-ing the evaporation. These parameters are known toyield particularly good optical performance. Film thicknesses were monitored on the vibrating quartz micro-balance during the evaporation and were subsequently determined to an accuracy of,typically,
+2
nm by acom-100
80—
V c C$ V O ~I
60-'D C lg O v40—
c cj ~~ E m 20— CI-bination
of
optical interference methods. The analyzed films were between 75 and 1200nm thick.Optical properties
of
the films were determined by spectrophotometry. Normal transmittance (T) and near-normal reflectance(R)
were recorded at room tempera-ture as a functionof
photon energy (fico) on a Beckman ACTA MVII double-beam instrument interfaced to a computer. Figure 1 shows typical data; (a) refers to a 75-nm-thick film made by evaporationof
pure In203, and (b) refers to a 110-nm-thick film made by evaporationof
In203+9
mo1%
Sn02. A large decrease inT
as the pho-ton energy is increased signifies the band gap. Oscilla-tions inR
indicate optical interference.Spectral transmittance and reflectance data were used
to coinpute the complex refractive index,
n+ik,
of
the films by useof
Fresnel's equations' and known propertiesof
the CaF2 substrates."
Spectrophotometric data were selected in such a way that the ensuing optical constants were accurate to within prescribed limits, and so that spurious solutions were avoided.It
was also verified that the Kramers-Kronig relation was obeyed. These aspects, which fall outside the scopeof
this paper, are discussed in some detail elsewhere. ' The optical constants n and kwere derived at
—
100 equally spaced energies in the (2—
6)-eV interval. Figure 2 shows smooth curves drawn through these data points for the two films reported inFig. 1. For
filmsof
pure In203 (solid curves), we find that kgoes up rapidly atduo-3.
8eV and that n is rough-ly 2 with a small peak at Acu-4 eV.For
Sn-doped In203 (dashed curves), the increase in k is displaced towards higher energy and now occurs at%co-4.
3 eV; the corre-sponding curveof
n again shows a small peak atfico-4
eV.The absorption coefficient
a
will be analyzed in detail below. This quantity is defined bya=4irk/A,
, 0 2.0 100 I 3.0 I 4.0 Energy (eV) I5.0 6.0 where A. is the wavelength. We computed
a
at—
100equally spaced energies for each sample, and joined the
in-80—
I
V C tg VI
60—
I
'a c lOI
40—
C lg O~ E th 20— f-110nm 3.0 c2.5—
I
'a C ~~ o 2.0—
V o 1.5—
0 lg ~ 1.0—
I
K 0.5—
—
0.8~
I
C ~~—
0.6~I
V lg LeI
V—
0.4 o tO CL i lg —0.2-
Ql 0 2.0 3.0 4.0 5.0 6.0 Energy (eV)FIG.
1. Spectral normal transmittance and near-normal re-flectance for films of(a)pure In&03 and (b)Sn-doped In2O3 on substrates of CaF2. Estimated experimental errors and the reproducibility among measurements on the same sample are considerably lessthan one percent unit.0
2.0 3.0 4.0 5.0 6.0
Energy (eV)
FIG.
2. Optical constants, n and k,versus photon energy for films ofpure In203 (solid curves) and Sn-doped Inq03 (dashed4, $ I I I 1 i i i i ) i I Sample Material n~(10 crn 3) (b)
E3-Q CI
.22— O 0 O C 0 ~~ CL 0 tD E~(k)=*'E
o rnc E'„(k)= 2FAy Ego jl~)I
kFE,
(k)= E',(k)+tZ,(k) E„(k)= E'„(k)+RE„(k) k 0—— I 3 4 5 Energy (eV)FIG.
3. Absorption coefficient versus photon energy for films of pure In203 and Sn-doped In203. Solid curves were drawn between individual data points whose scatter was lessthan the width ofthe lines. The curve denoted Dwas obtained for athick film which became virtually opaque for Ago
)
4.7 eV.Inset table shows the electron densities in the various films.
dividual points by smooth curves. Figure 3 shows such data for four different films. Curves
3
andC
pertain to the filmsof
pure In203 and Sn-doped In203, respectively, which were reported earlier in Figs. 1 and2. It
is ap-parent that the band gaps lie at different energies, and that the onsetof
strong absorption is rather gradual.For
low energies it was found that 1na versus ficodisplayed an approximately linear relation. 'The data in
Fig.
3 can be systematized by considering the electron density,n„
in the various films. Its magni-tude was obtained from determinationsof
the plasma en-ergy (i.e.,the energy for which the real partof
the dielec-tric function equals zero), as described in earlier papers. 'For
filmsof
pure Inz03, we have n,=(0.
4+0.
05)&&10
cm as a result
of
doubly charged oxygen vacancies. 'For
Sn-doped In203, we have electron densities between(1.7+0.
2) and(8.
0+0.
5)&&10 cm primarily as aconsequence
of
the Sn atoms which enter substitutionally as Sn + on In + sites, and hence act as singly charged donors. ' Actual valuesof
n, are given in the insetof
Fig.
3.
It
is concluded that increasing electron densities lead to aprogressive enhancementof
the energy gap.III.
THEORY OF SHIFTED BAND GAPS As a starting point for discussing the relation between optical band gap and electron density we consider the band structureof
undoped In203. The band structure is unknown in most respects. The only available informa-tion concerns the direct and indirect band gaps' and the region around the bottomof
the conduction band, which is thought to be parabolic with an effective mass(m,
')
of'
-0.
3m, where m is the free-electron mass. We can only make an assumption regarding the shapeof
the valence band; we take it to be parabolic and characterized by an effective mass (m,*)
of
unknown magnitude. Figure 4(a) illustrates this band structure. With the topof
theFIG.
4. (a) shows the assumed band structure of undopedIn203 in the vicinity ofthe top ofthe valence band and the bot-tom ofthe conduction band. (b)describes the effects ofSn dop-ing: The valence band isshifted upward by many-body effects
while the conduction band is shifted downward. Shaded areas denote occupied states. Band gaps, Fermi wave number, and dispersion relations are indicated.
valence band as reference energy, the dispersions for the unperturbed valence and conduction bands are
E„(k)=
—
A' k /2m,'
(2)E,
(k)=Ego+fi
k /2m,*,
(2')respectively. Ego is the band gap
of
the undoped semicon-ductor, k is the wave number, and superscript0
denotes unperturbed bands. Single crystalline platesof
In203 have' E~o3.
75 eV.Po——
lycrystalline thin films may have a somewhat different Eso, its actual magnitude is depen-dent on the detailed preparation conditions. A detailed discussionof
these shifts is not possible, but we note that the band gap can be altered by local strain induced by im-purities, point defects, and poor crystallinity.In the doped material we have to consider three dif-ferent effects: First, the shapes
of
the valence and con-duction bands may not be accounted for by precisely the same effective masses as in the undoped material. Indeed, it has been found''
that m, isweakly dependent on the electron concentration and goes up to &0.
4m at n,&3&10
cm.
A corresponding variation for m,* cannot be ruled out, but nothing is known. Second, above the Mott critical density' the partial fillingof
the con-duction band leads to a blockingof
the lowest states and hence a wideningof
the optically observed band gap. This is the well-known Burstein-Moss (BM) shift. Third, again above the Mott critical density the valence and con-duction bands are shifted in energy as a resultof
electron-electron and electron-impurity scattering. In In&03, these tend to partially compensate theBM
shift. Figure 4(b) shows schematically the rolesof
the second and third effects.We first neglect the role
of
electron-electron and electron-impurity scattering. The energy gap for direct transitions in the doped material is then given in termsof
the unperturbed bands asEg E—
—
,
(kF)—
E„(kF
),
(3)where kF
—
—
(3m. n,)'
is the Fermi wave number. Alter-natively, we may write0 BM
Eg
—
—
Egp+
AEgwhere the BMshift isgiven by g2
gEBM
(3+
)2/32muc
with the reduced effective mass
1 1 1
~+
Pl~ Plv mc
(3')
E,
(k,co)=E,
(k)+iriX,
(k,co),
(6')where
iiiX„and
iiiX, are self-energies due to electron-electron and electron-electron-impurity scattering. We now obtain, insteadof
Eqs. (3) and(3'),
the shifted optical gapEs
E,
(kF,co)E,
(—
—
kF,co),
—
or, alternatively,(7) Equation (4) predicts an energy-gap shift proportional to ~2/3
With the purpose
of
giving a more correct theoretical model for the shifted band gaps, we now include electron interactions and impurity scattering. The free electrons in the doped material cause a downward shiftof
the conduc-tion band as a resultof
their mutual exchange and Coulomb interactions. This shift is further accentuated by the attractive impurity scattering. The valence band is infiuenced in the opposite way. The effectof
the various interactions can be described simply by replacing the bare-band dispersions in Eqs. (2) and(2')
by the corre-sponding quasiparticle dispersionsE„(k,
co)=E„(k)+iriX„(k,
co)and
where Rg runs over a random distribution
of
Sn + on In + sites. We estimate V;„by
the Heine-Abarenkovpseudopotentials' appropriate to the two ionic species. Figure 5 shows the Fourier transform
of
the unscreened difference Vs„(r
)—
V,„(
r
) divided by the bare Coulomb potential—
4me /q for aunit point charge. This ratio is close to unity for k&2kF, which is the pertinent range in the computations, thus proving that the scattering centers behave as screened point charges. We stress that this simplifying feature would break down for electron con-centrations much higher than thoseof
the present films as well as for doping elements whose pseudopotentials are drastically different from thatof
the substituted ion. Ex-plicit results forfiX"
,
and A'X",
were obtained by using m, fromRef.
17(together with areasonable extrapolation to-ward high electron densities). The computationsof
A'X„"and A'X,
"
require in principle that the full frequencydependence
of
the background is included.For
simplicity we suppress this dependence. In all calculations we set the dielectric constant for the In203 host equal to thestat-ic
value. 'The computations now proceed as discussed at length in
Ref. 2.
Thus the screening properties were included by using the random-phase approximation. As commonly done, irido in the expressions for the self-energies was putequal to fPkF /Zm, ~,i The in.teractions do not only shift
the positions
of
the valence and conduction bands but also distort them to some degree, and hence the k dependenceof
the self-energies should be retained in principle.As an application
of
the above formulas, we now report on computed band-gap shifts versus n, The eff.ective valence-band mass remains as an unknown parameter. Figure 6contains results for three different valuesof
m„.
It
is found that an increasing electron density leads to an enhanced gap shift, and that at constant electron density the gap shift is largest for the smallest m„*.The shifted band gap is further elucidated in
Fig. 7. It
shows b,Ez form„'=0.
6m as well as its three contribu-tions:Eg
—
—
Eg0+
AEg,
with(7')
b,
Es
bEg+AX,
(k—
—
F,
co) fiX, (kF,co).
—
Electron-electron (ee) and electron-impurity
(ei)
scatter-ing are taken as additive processes within our perturbation treatment,i.
e.,iriX,
(k,
co)=iiiX„''(k,co)+A'X",(k,co)V;m~(r)
=+[Vs„(r
—
Rg)—
V&„(r—
Rg)],
(10) AX,(k,co)=AX,
"(k,
co)+A'X,"(k,
co).
We refer to
Ref.
2for details.Wefirst consider fiX'„'and fiX,
".
The Sn+
ions—
which replace substitutionally someof
the In + ions—
act effec-tively as singly charged scattering centers. Their scatter-ing potential can be written0.
5
I ) I I
1.
0
1.
5
a(au)
FIG.
5.F(q)
as a function of q in atomic units (a.u.). F(q)isthe ratio between the Fourier transform ofthe bare scattering potential Vs„(r)
—
Vi„(r) and the Fourier transform of the Coulomb potential for a point charge.F(q)
was evaluated for n,=
10 'cm,
but is not sensitive tothe electron concentration. Arrow points at k~for n,=7)&10 cmIV. COMPARISON OF THEORY AND EXPERIMENTS ~ 10 E V O Ol C CO Q 'U C 0 4 Q LLI and gEBM g A'X"(kF)=I)i&,
"(kF)
—
IIi&",(kF),
2—
I I I I I I I I l I I I I I I I I I 00.
20.
40.
6 0.8Energy gap shift hEg (eV}
FIG.
6. AE~ versus n, computed for three magnitudes ofthe effective valence-band mass.
, X
I(I
II'If&
I'
(N
—
CO,f
)+T
where
i.
accounts for the broadeningof
the initial andfi-nal states,
P
is the probability that the state is occupied, and cof;=(Ef
E;
)/fi.—
In the limiti
~
acEq. (11)
goes over to the usual golden-rule expression.We now identify the initial states with the filled valence band and the final states with the partially filled conduc-tion band.
It
is straightforward to prove thatR
ccI
dx(x+fico
—
W)'i
(1
P,
),
—
Xp
x
+I
(12)where wehave introduced the notation
The experimental data on 0,versus %co show a gradual onset
of
strong absorption, and hence it is not obvious how to locate a unique optical band gap which can be compared with Eg as derived in the preceding section. In order to understand the experimental broadening, we con-sider the quantum-mechanical transition rateR
for transi-tions between an initial(i)
and a final(f)
state.Accord-ing to time-dependent perturbation theory, we have
fiX"(kF):
—
fiX,"(kF)
—
I)i'",(kF).
x:—
Ak
+
O'—
Ac@,2&i (13)
The
BM
shift is dominating except for the lowest electron concentrations. AEgBM versus n,2/3 does not give a straight line, which is aconsequenceof
the empirical rela-tion' between m, and n, [and (Ref. 21)between thestat-ic
dielectric constant and n,]
The s.elf-energies are seen to vary approximately as n, ; this dependence is rather fortuitous and is the net resultof
several competing ef-fects.xo
=
4Eg+
8'
—
Ace,I
=&/w.
An analogous expression for
R
has been given by Finken-rath who used classical arguments. The minimum dis-tance between the valence and conduction band in the doped material—
denoted WinFig.
4(b)—
is given by the approximate relationW=EsP+fiX,
(kF)—
I)iX„(kF).
10— CV I E O ~8—
I I I I I I I I I I I I IIn writing R as in
Eq.
(12) we have also ignored the k dependenceof
I .
Thermal excitations above the Fermi energy are represented by aFermi function according to
A' k
P,
=
exp,
—
p
+1.
2m,'
C6—
'Q C 0 Le V O LLlwhere k&T is Boltzmann's constant times the tempera-ture, and
p
is the chemical potential. Atkg
T «A'
kF/2m,':
—
eF we have2' I I I I I I I I I I I I I I I
0 0.5 1.0
Energy (eV)
FIG.
7. AE~ versus n, computed for m„*=0.6m, andanalogous plots for the contributions to EE~from the Burstein-Moss shift (EEL ),electron-impurity scattering [fiX"(kF)],aud electron-electron scattering [A'X"(kF)
].
2
kmT
P
~QF3 2CF
Inthe limit
I
~0
andT~O,
Eq. (12)goes to(fico W)'~ for
fico) W+b,
E—
0
for~&
g
+kg
(19) (19')
(20) or, tolinear order,
R
oc1—
(W+bE
fici)).
—
~l
g (21)According to
Eq.
(20) the inflection point is given byiricoo—
=
W+b,
Eg.
A numerical studyof
the full expres-sion inEq.
(12)shows that this estimate is indeed good as long asI
is not too large. A first estimateof
the absorp-tion edge is thus found simply by locating the maximumof
da/d(irico). Drawing the tangent in the linear regionof
the measured absorption, the intercept with the Ace axisgives
W+
hE~—
(irl/2),
and hence afirst approximate value forI
.
The numerical studyof
Eq. (12)shows that this construction somewhat overestimatesI'.
The final values forW+EEg
andI
shouldof
course be deter-mined more accurately by fittingEq.
(12) to measured quantities, as will be described shortly. Our approximate procedure yields, however, values forW+b,
Ez™
andI
which are most useful in connection with such a fine tuningof
the parameters.Figure 8 shows comparison between theory and experi-ments for two
of
the samples, whose spectral absorption coefficients were given earlier inFig. 3.
Open and filled circles correspond to evaluations based on spectropho-tometric measurements. Solid curves denote theoretical data, which have been fitted by selecting appropriate mag-nitudesof
8'+AEg
andI
.
In doing this fitting we have primarily regarded the middle partsof
the curves.It
is seen fromFig.
8 that theory and experiments can be brought in good agreement around the steepest partsof
the curves, whereas the"tails"
toward high and low ener-gy cannot be reproduced to an equal precision. Oneof
the possible reasons for this is that we have ignored the k dependenceof
I
and AX. The dashed curves pertain to al-ternative valuesof
I
.
It
is apparent that the computedThis relation has formal similarities to its well-known counterpart for the undoped semiconductor (forwhich we substitute
W~Ego
and EEg~0).
The effectof
I
is thus to smear the sharp absorption edge inherent in Eqs. (19)and(19').
Putting
R
o:a
(Ref. 26) gives a complete scheme for computing the spectral absorption coefficient in termsof
two parameters. Oneof
these is taken asW+bEg
.
It
gives the energy around which the transition from low to high absorption is centered. This parameter is convenient since it can be directly compared with Eg as computed in the preceding section. The second parameter isI,
which gives the widthof
the transition.For
qualitativegurposes, and for making rough esti-matesof
W+&Ez
andI,
it may be useful to approxi-mateEq.
(12)in the following way.For
small valuesof
I
the Lorentzian is a sharply peaked function. Hence, we may consider the square-root expression as slowly varying and take it outside the integral sign. In the regionof
the band-gap we therefore have, with xo inserted in the square root arid puttingT
=0,
I
BM
R
cc(AE
gM))g2 1—
—
2tan )W+AEg
r
I i I I I I I I I t I I i l f I I I i ( I I l I I 3— 0 I Theory CI
2— V ~~ V O 0 V C O f CL 0 V) r (eV 0.25 0.15 0—
0.05 '" I s i i i I I I I i I i & I I I i i i i I i 3 3.5 4 4.5 5 5.5 Energy (eV)FIG.
8. Absorption coefficient versus photon energy for films ofSn-doped In203. Circles refer to experimental data (cf. Fig. 3). Solid and dashed curves were computed by using8'+EEg
=4.
04and 4.46eVfor the samples, together withI
(ineV)as shown in the figure.
TABLE
I.
Parameters used for fitting theoretical curves to experimental spectral absorption coefficients. The band gap corresponds to8'+DE~,
and the width in the curve isgoverned by
I .
Sample
W+
aE,
'
(eV)3.88 4.04 4.46
-4.
55r
(eV) 0.07 0.13 0.37-0.
35 curves are strongly dependent onI
and, conversely, that reasonably uniqueI"
s can be extracted from the experi-ments. When including a finite temperature we have usedT=300 K.
This parameter has a marginal influence on the curve shapes, since the degeneracy temperatureof
the electron gas is much higher than 300K.
Table
I
containsW+EEg
andI'
as obtained from "bestfits"
between theory and experiments. Improved ac-curacy was achieved by comparing theoretical and experi-mental results on da/d(irico) versus fico. Data for sampleD
are less certain than for the others, since we have access toa
only for fico&4.
7 eV.It
is inferred thatW+b,
Eg andI
both become larger as n, is increased(cf.
the inset table in Fig.3).
In
Fig.
9 we compare the evaluationsof
8'+EEg
with computationsof Eg.
From Fig. 6 it is evident that the slope in the theoretical curves depends on m„'.It
is then clear that m,* should be chosen in the range(-0.
6—
0.
7)m.For
m,*=0.
7m we obtain best agreement by setting Ego—
—
3.
82eV.For
m,*=0.
6m the correspond-ing energy is3.
75 eV. Thus itappears that we have to as-sume aband gap for the undoped material which is either the same or somewhat higher than the value3.
75 eV per-tainingto'
single crystalline In203. An uncertaintyof
fQ ) i [ ~ ~ i [ i f i f I ) I f I I 8 Ol E Q O 6 C 4 O C 0 lm V O Ill 2 The mmmm Q+v 0 I i I i I i I i I .i I i I i I i I i 3.8 4.0 4.2 4.4 4.6
Energy gap (eV)
FIG.
9. Electron density to the power 3 versus energy gap 2for In203 with different degrees ofdoping. Circles represent ex-perimental data (cf.Table I). Vertical bars indicate the uncer-tainty in extracting n, from observed plasma frequency; they ac-count for conceivable errors in effective conduction-band mass and dielectric constant of undoped In203. The uncertainty in
the determination of
8'+AEg
is less easy to estimate withconfidence. Solid and dashed curves indicate the results of com-putations using the shown values of
m„and
Ego. Curves are confined to n,)
3&10' cm since the theory only applies to electron densities well above the Mott critical density (cf.Ref.18).
this order
of
magnitude is by no means surprising. garding the value for m,, there are unfortunately no band-structure calculations with which tocompare.The parameter
I
can be understood in the following way. Becauseof
scattering, initial and final states are col-lision broadened.For
ionized impurity scattering one has for the conduction-band electronsn,
f
i V(q) i5(E,
(k)
—
E,
(k+q)),
(22) where V(q) is the dielectrically screened scattering poten-tial for asingle impurity.
For
the valence states one findsv„'(k) =v;
'(k)m„*/m,*.
The electron-impurity scattering therefore gives rise tothe broadeningI
=
—
=
—
fi[r, (k)+~„(k)]
) i=
—
~,
(k)(1+m„/m,
).
A'
c v
2 c
(23) Putting m,
=0.
6m and m, as given by Ohhata et al.,'we find for
k=kF
thatI
=0.
11 eV for n,=1.
7~10
cm andI'=0.
17eV for n,=
.
6X210 cm.
Accord-ing toTableI
the experimental estimates are0.
13 and0.
37 eV, respectively. Impurity scattering thus gives the correct orderof
magnitude forI
.
We recall that the k dependenceof
I
was neglected inEq.
(12).For
this reason the agreement above may be considered as rather satisfactory.The electron lifetime for the conduction-band electrons
inay also be estimated from transport measurements. The expression for this transport time
r,
„(k)
is the same asfor
r,
(k) inEq.
(22) except for a factor(1
—
cosO), whereL9 is the angle between the vectors k and
k+q.
Sinceionized impurity scattering isbelieved to be the most im-portant damping mechanism, we have calculated the in-fluence
of
this factor at different impurity concentrations.For
low values (—
10'
cm )the scattering is highly an-isotropic implying thatr,
t',(kz)
&«,
'(kF).
Increasing the impurity concentration reduces the anisotropy so that at high doping (—
10 cm )we haver,
,
',(kF)&~'(kz),
i.
e., nearly isotropic (s-wave) scattering. Hence we may user,
„
to estimate the broadeningI
by meansof Eq.
(23). From measurementsof
the frequency-dependent resistivity, and hencer,
„,
we obtain the followingI'
values for the different samples listed in TableI:
A,0.
05;8,
0.
08; C,0.
13; D,0.
16eV. CasesB
andC
agree quite well with our previous theoretical estimates, and as before, there is, within about a factorof
two, an overall agree-ment withI
as estimated from the measured absorption coefficient. With increasing concentration there is a marked tendencyof
an extra contribution toI
.As we have demonstrated above, ionized impurity scattering is the most important contribution to
I,
but other effects may also come into play. AtT=O
K,
electron-electron interactions give no lifetime broadening at the Fermi level, but since the folding inEq.
(12)also in-volves states away from the Fermi energy,I
may be af-fected. Effectsof
electron-phonon interactions are simi-lar. There may also be smearing effects due to local de-formationsof
the lattice becauseof
the randoinly intro-duced Sn ions. The electric fieldof
these ionized donors and the difference in atomic sizeof
Sn relative to In affect the absolute valueof
the band gap for each microcellof
the crystal. Hence the band gap should be statistically distributed around the value8'.
In our theoretical expres-sions we have taken into account the overall band-gap narrowing due to the electron-impurity and electron-electron interactions, but neglecting these possible fluctua-tionsof
the energy levels. Clearly this latter effect should increase with the doping concentration. We shall not pur-sue such effects or electron-electron interactions any fur-ther.For
our purposes it suffices to consider ionized im-purity scattering only.V. CRITICAL ASSESSMENT OF EARLIER WORK
Band-gap widening associated with doping has been ob-served in earlier work on doped In203 (Refs. 8, 15, 16,and 27
—
35), ZnO (Refs. 36 and 37), CdO (Ref. 23), Sn02 (Refs. 30, and 38—
43),and Cd2Sn04 (Refs. 44—
46). The evaluationof
shifted band gaps, as well as the interpreta-tionof
these data, are at variance with the results given above, and hence an assessmentof
the earlier work is ap-propriate.For
doped In2Q3 there have been numerous' ' 'at-tempts toevaluate band gaps by plotting
a
versus Aceand associating the band gap with the energy obtained through a linear extrapolation toward zero.It
appears that this procedure has its root in a mistaken interpretationof Eq.
(19),which would yield W and notW+b,
E&™.
Besidesi I I I I I I I I I I I I I I 10— E 0 co 8 0 O 8 C4
~6—
tO e Experiment: C 0I
0a
4 This wo P e~IP
~ Vainshte o o Qhhata etal~ 2— SSIn203 (weihsr aLeyi
p. i
l
i I~ I i I i I i I i I i I i I i I i I i I i I3.4 3.6 3.8 4.0 4.2 4.4 4.6
Energy gap (eV)
FIG.
11. Electron density to the power 3 versus energy gapfor In203 with different degrees ofdoping. Approximate energy
gaps were obtained from extrapolations in plots ofcz versus fico.
Symbols denote results from the present work on pure and Sn-doped films prepared by reactive e-beam evaporation, as well as for sputtered films studied by Vainshtein and Fistul (Ref.27), Ohhata et al. (Ref.17), and Smith and Lyu (Ref. 33). The elec-tron density was determined from measurements ofthe Hall ef-fect in Refs. 27and 33,and from measurements ofthe plasma
energy in Ref. 17. Error bars on n, are the same as those in
Fig. 9. Arrow points at the energy gap for single-crystalline
In203 asobserved by Weiher and Ley (Ref. 15). Solid curve was
computed from our theoretical model with Ego
—
—
3.75 eV andpl~
=
1.0m. OI5—
E Q O O 4 Ol cI
~~3 Q ~~I
0 Q C2— 0 ~~ 0 rk in &Fistul'Smith & Lyu
I 4.5 I 5.0 I I I 4.0 Energy {eV)
FIG.
10. Absorption coefficient squared versus photon ener-gyfor the samples A—
Dreported earlier in Fig.3. Dashed linesindicate linear extrapolations. The zero energy intercepts lie at
3.81,3.91, 4.23,and 4.30eV.
I
3.0
band gaps by using the theory outlined in Sec.
III.
Best agreement was obtained with m,=1.
0mandEg0-3.
75 eV. This valueof
m,'
is significantly different from that found by the detailed analysis in Sec. IV,which points to the importanceof
using a proper evaluationof
the band gap. Correspondence with the dataof
Ohhata et al. ' could be obtained with their band gap Ego——3.
67 eV. Correspondence with the dataof
Vainshtein and Fistul' andof
Smith and Lyu requiredEso-3.
55eV.The experimental observation that the band gap scales with n, in doped In203 has led to the erroneous conclusion
—
iterated in numerous papers'
'
'—
thatthe
BM
shift alone would determine the band-gap widen-ing.If,
for the sakeof
argument, we neglect the self-energy effects, we can obtain agreement between theory and experiments only by assuming a negative magnitudeof
m, somewhere between—
—
0.
6 and—
—
1.
0m. Thiswould imply that the valence band is curved in the same direction as the conduction band. Inclusion
of
the self-energies leads to apositive valueof
m„*, as discussed inSec.
IV.
This proves the importanceof
considering electron-electron and electron-impurity scattering, and that qualitative differences can occurif
they are neglected.As far as we know, self-energy effects have been neglected in almost all earlier analyses
of
band-gap shifts in doped oxide semiconductors.For
three materi-als, ZnO, CdO, and Sn02, there are band-structure calcu-lations which allow—
at least in principle—
an assess-this, we have shown thatEq.
(19)is inapplicable becauseof
lifetime broadening effects. Nevertheless—
evenif
a strict theoretical backing is absent—
it turns out that the simple extrapolation procedure yi.elds values which, at least for our present samples, are not grossly different from those obtained from the more elaborate analysis given in the preceding section. This point is illustrated in Fig. 10,where we replot the absorption coefficients given earlier inFig.
3asa
versus fico. The shifted band gaps, derived from the indicated extrapolations, lie below those given in TableI
by amagnitude ranging from0.
07eVfor the lowest electron density to0.
25 eV for the highest elec-tron density. %'e may understand this result by consider-ing the approximate relation inEq.
(21).If
squared, andif
only terms ug to linear order are retained, the root equalsW+AEs
irIl.
4,i.
e.,t—
hea
procedure leadstoband gaps which are too small by an amount
—
I
.
Thevalues for
I
in TableI
are obviously consistent with such arguments. Other techniques to obtain shifted band gapsare based on extrapolations towards zero in plots
of
(Ref. 30)a
versus fico, and determinationsof
the energy corre-sponding to (Ref. 28)10%
or (Refs. 29 and 31)50%
transmittance through the samples or (Ref. 32)a=10
cmIt
is interesting to compare our band gaps with earlier evaluations. In order to have a consistent setof
data, we focus on works in which approximate band gaps have been obtained by evaluatinga
versus Ace. Figure 11 shows a collectionof
such data plotted against n,It.
appears that our results (solid triangles) are in rather good agreement with thoseof
Ohhata etal.
' (open circles), whereas the resultsof
Vainshtein and Fistul' (solid squares) andof
Smith and Lyu (solid circles) are dis-placed toward lower energy. Still larger displacements have been reported byRay et al. We find, approximate-ly, that the shiftof
the energy gap scales withn,
.
This kindof
relation has been reported also by others. ' 'ment
of
the earlier analyses. Sn02 doped byF
and Sbhas been investigated recently by Shanti et al. ' Interpreting the band-gap variation asa BM shift, they concluded thatmU' would be
——
0.
9m. A similar analysisof
CdO, per-formed much earlier by Finkenrath, also gave anegative value of m„*. Neitherof
these predictions are corroborat-ed by recent band-structure calculations for (Ref. 48) Sn02 and CdO. Inclusionof
self-energy effects wouldtend to improve the correspondence between the theoreti-cal analyses
of
shifted band gaps and band-structure data, but no detailed results are available at present.ACKNOWLEDGMENTS
This work was financially supported by grants from the Swedish Natural Science Research Council and the ¹
tional Swedish Board forTechnical Development.
R.
A.Abram, G.J.
Rees, andB. L.
H.Wilson, Adv. Phys. 27, 799 (1978).K.-F.
Berggren andB. E.
Sernelius, Phys. Rev.B
24, 1971 (1981).3E.Burstein, Phys. Rev. 93,632(1954);
T.
S.Moss, Proc.Phys. Soc. London Ser.B67,775(1954).4I. Hamberg and C.G.Granqvist, Thin Solid Films 105, L83 (1983).
5I. Hamberg and C. G. Granqvist, Proc. Soc.Photo-Opt. In-strum. Engr. 428,2(1983);Appl. Phys. Lett. 44,721{1984}.
6C.G.Granqvist, Appl. Opt. 20, 2606 (1981);Proc. Soc. Photo-Opt. Instrum. Engr. 401,330 (1983).
H. Kostlin, in Festkoerperprobleme, edited by P. Grosse (Vieweg, Braunschweig, 1982), Vol. 22,p. 229.
SK.L.Chopra, S.Major, and
D. K.
Pandya, Thin Solid Films 102,1(1983}.I.
Hamberg, A. Hjortsberg, and C. G.Granqvist, Appl. Phys. Lett. 40, 362 (1982); Proc. Soc. Photo-Opt. Instrum. Engr. 342, 31(1982};I.
Hamberg and C. G.Granqvist, Appl. Opt. 22, 609 (1983).ioM. Born and
E.
Wolf, Principlesof
Optics, 6th ed. (Pergamon, Oxford, 1980}.~~Results obtained from the Harshaw Chemical Company. 2A. Hjortsberg, Appl. Opt. 20, 1254(1981);
T.
S.Eriksson andA.Hjortsberg, Proc. Soc.Photo-Opt. Instrum. Engr. 428, 135 (1983).
This phenomenon is known as "Urbach tails"; cf. M. V. Kurik, Phys. Status Solidi 8, 9(1971}.
G.Frank and H.Kostlin, Appl. Phys. A27, 197 (1982).
R.
L.Weiher andR.
P.Ley,J.
Appl. Phys. 37, 299 (1966). Z. M.Jarzebski, Phys. Status Solidi A71,13(1982).Y.Ohhata,
F.
Shinoki, and S.Yoshida, Thin Solid Films 59, 255(1979).N.
F.
Mott, Metal-Insulator Transitions (Taylor and Francis, London, 1974). The critical density n, forthe metal-insulator transition in a doped semiconductor can be estimated from Mott's criterion n,'a~-0.
25. The effective Bohr radius isatt
—
—
eoirt /m,e,
where eo is the static dielectric constant ofthe host material [P. P. Edwards and M.
J.
Sienko, Phys. Rev. B 17,2575 (1978)]. For n type In203 w-e thus obtain n,=4&(10' cm.
One should recall, however, that our e-beam evaporated In203 films have n,=4&&10' cm as are-sult ofdoubly charged oxygen vacancies.
A. O.
E.
Animalu and V.Heine, Philos. Mag. 12,1249(1965); M.L.
Cohen and V. Heine, in Solid State Physics, edited byH. Ehrenreich,
F.
Seitz, andD.
Turnbull (Academic, NewYork, 1970),Vol.24.
~oRecent calculations ofthe frequency-dependent resistivity also
show that the impurities may betreated as point charges toa
high degree ofaccuracy [1.. Engstrom and
I.
Hamberg(un-published)].
The static dielectric constant is adequate forelectron-impurity scattering; cf.
B.
E.
Sernelius andK.-F.
Berggren, Philos. Mag. 43, 115 (1981). Similar to the case of m,*,
the static dielectric constant isrelated to n,.
We usethe dependencein-herent inthe results by Ohhata et al. (Ref.17).
22See, for example, G. H. Wannier, Elements
of
Solid StateTheory {Cambridge University Press, Cambridge, England, 1959), p. 212.
2 H.Finkenrath, Z.Phys. 159, 112 (1960).
Equation (16}would be exact in arigid-band model (i.e.,
as-suming only vertical displacements ofthe energy bands).
2sSee, for example,
F.
Wooten, Optical Propertiesof
Solids (Academic, New York, 1972),p. 115.In practice it makes little difference whether weequate
8
witha
(-
k/A, )or@~co(-nk/1,
), since n isnearly constant. V. M.Vainshtein and V.I.
Fistul', Fiz.Tekh. Poluprovodn. 1,135(1967)[Sov. Phys.
—
Semicond. 1,104(1976)].2 H.Kostlin,
R.
Jost,and W.Lems, Phys. Status Solidi A29,87(1975).
~W.G.Haines and
R.
H.Bube,J.
Appl. Phys. 49, 304 (1978).J.
-C.Manifacier,L.
Szepessy,J. F.
Bresse, M.Perotin, andR.
Stuck, Mater. Res.Bull. 14, 163(1979).'A.
J.
Steckl and G. Mohammed,J.
Appl. Phys. 51, 3890 (1980).M.Mizuhashi, Thin Solid Films 70, 91 (1980).
3
F.
T.
J.
Smith and S.L.
Lyu,J.
Electrochem. Soc. 128, 2388(1981).
3~J.Szczyrbowski, A. Dietrich, and H. Hoffman, Phys. Status Solidi A69, 217 (1982};78,243(1983).
35S. Ray, R. Banerjee, N. Basu, A.
K.
Batabyal, and A.K.
Barua,J.
Appl. Phys. 54, 3497 (1983}.A. P. Roth,
J.
B.
Webbs, and D.F.
Williams, Solid State Commun. 39,1269(1981);Phys. Rev.B
25, 7836 (1982). O.Carporaletti, Sol. Energy Mater. 7,65(1982). 33T. Arai,J.
Phys. Soc. Jpn. 15,916(1960).39H.Koch, Phys. Status Solidi, 7,263(1964).
~S.
P.Lyashenko and V.K.
Miloslavskii, Opt. Spektrosk. 19, 108(1965) [Opt.Spectrosc. (USSR) 19,55(1965)].E.
Shanti, A. Banerjee, V.Dutta, andK.
L.Chopra,J.
Appl. Phys. 51,6243 (1980);53, 1615 (1982);E.
Shanti, A.Banerjee,and
K.
L.Chopra, Thin Solid Films 88, 93 (1982).K.
B.
Sundaram and G.K.
Bhagavat,J.
Phys. D 14, 921 (1981).K.
Suzuki and M.Mizuhashi, Thin Solid Films 97, 119 (1982).~A.
J.
Nozik, Phys. Rev.8
6, 453 (1972).~5N. Miyata,
K.
Miyake,K.
Koga, andT.
Fukushima,J.
Elec-trochem. Soc.127,918 (1980).~
E.
Leja,K.
Budzynska,T.
Pisarkiewicz, andT.
Stapinski,4~Work on ZnO by Roth et al. (Ref. 36) has considered
self-energy effects,
4~J.Robertson,
J.
Phys. C 12, 4767 (1979).4~J.C.Boettger and A.
B.
Kunz, Phys. Rev.B27, 1359 (1983).5OD. M.Kolb and H.-J.Schulz, Curr. Top. Mater. Sci. 7,227 (1981).