McGuire, G.E. (Eds.) - Semiconductor Materials and Process Technology Handbook-William Andrew Publishing - Noyes (1988)
McGuire, G.E. (Eds.) - Semiconductor Materials and Process Technology Handbook-William Andrew Publishing - Noyes (1988)
and
Process Technology
Handbook
for
Very Large Scale Integration (VLSI)
and
Ultra Large Scale Integration (ULSI)
Edited by
Gary E. McGuire
Microelectronics Center of North Carolina
Research Triangle Park, North Carolina
I 1
L-l "P
NOYES PUBLICATIONS
Westwood, New Jersey, U.S.A.
Copyright 01988 by Noyes Publications
No part of this book may be reproduced in any form
without permission in writing from the Publisher.
Library of Congress Catalog Card Number: 87-31529
ISBN: 08155-l 150-7
Printed in the United States
1098765
Bibliography: p.
Includes index.
1. Integrated circuits--Very large scale
integration--Design and construction--Handbooks,
manuals, etc. I. McGuire, G.E.
TK7874S4178 1986 621.395 87-31529
ISBN 08155-1150-7
MATERIALS SCIENCE AND PROCESS TECHNOLOGY SERIES
Editors
Related Titles
INTRODUCTION. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .I
Gary E. McGuire
References. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .7
vii
viii Contents
INDEX..............................................669
Introduction
Gary E. McGuire
Tektronix, Inc.
Beaverton, Oregon
1
2 Semiconductor Materials
lo- 0
4K DRAM
1 K SRAM
7‘ .O
5- 6K DRA: &? 4K SRAM
4. 16K DRAM
3. 00 16K SRAM
Figure 1: Projected decrease in feature size for dynamic and static random access
memories as a function of the year of introduction.
2 s 4 S 6 7 8
CHIP LENGTH (MM)
Figure 2: Plot of transistor channel length versus overall chip size showing the
trend toward larger chips even though the feature sizes have decreased.
Introduction 3
27
Eo lOOO-
b
rn
0) 600-
d
3 600 -
%
5 400 -
0 1 2 3 4 5 6
3-
2-
l-
0 I I I I
0 1 2 3 4 5
I I I I
1960 1970 1860
VDS = 0.1 v
DRAIN
CURRENT
(IDS)
bATE
GATE VOLTAGE
Figure 6: Plot of the drain current versus gate voltage illustrating the lack of a
sharp cut-off voltage for VLSI devices. Reference: VLSI Technology and Design,
IEEE, 0. Folberth and W. Grobman (1984).
JL (nAlcm2) @
P/A (X 10 cm - 1)
REFERENCES
William C. O’Mara
Aeolus Laboratory
Palo Alto, California
1. INTRODUCTION
This is the silicon age. Just as previous historical periods were named
by the characteristic material, we may assume that silicon will be seen as
paramount in importance to the era to which some refer as the second
industrial revolution. Certainly the transformation in the way people both
workandrelaxis beingchangedsignificantlybyelectronicdevicessuchas
computers, control systems, and audioandvideo products. These electronic
devices are all based on silicon, especially as used for integrated circuits.
Although a large amount of technical literature exists on silicon devices,
comparatively little has been written on the material itself. This is especially
true of material of an introductory nature. Thischapterattemptsasurveyof
the way silicon is made, and includes information on material properties,
especially as modified bythe presence of small amounts of oxygen. Silicon
turns out to be a fascinating substance, and readers of thisvieware invited
to turn to the references for further information.
The method for single-crystal silicon growth was invented by Teal and
Beuhler in 1951.’ It was an extension of earlier work by Teal and Little in
which single crystals of germanium were prepared for transistor manufac-
turing.Germanium,likesilicon,occupiesaGroupIVposition in theperiodic
table, and has semiconductor properties that facilitated the initial device
manufacturing. However, the superior properties of silicon were soon
realized and the method was extended to this element. Pure material was
melted in a quartz crucible in an inert ambient, and a seed crystal was
lowered to begin controlled freezing of the melt. Dopants were added as
needed to control the electrical properties of the material. Reference 1
describes this early work in detail.
Production of silicon ingots today follows this same method, although
extensive improvements have been made in equipment, starting material
and process control. This section describes current practice, with emphasis
on material perfection andcontrolled impurityincorporation.Some people
refer to this method as Czochralski silicon growth, after an earlier experi-
mental method. However, this method was not designed for, nor did it
produce, single-crystal material.2 The ability of Teal and Little to make
single crystals repeatably was crucial to the growth of the solid-state-
device industry.
Other processes have been developed forsilicon-crystal growth. The
most important of these is the floating-zone method, in which thecrystal is
solidified from a small molten zone resting on the crystal itself. A polycrys-
tallinefeed rod is lowered from above into an RF induction coil which melts
its lower end. The rate of lowering is matched to the rate of withdrawal of
the crystal from below in order to maintain a constant melt volume. This
method differs from that of Teal and Little in that no crucible is employed;
the melt is suspended and maintained by surface tension. Because the
melt is not in contact with quartz, no oxygen is incorporated into thesilicon.
Some devices, such as high-voltage, high-power transistors, rectifiersand
thyristors require oxygen-free starting material. The majority of devices,
however, including virtually all integrated circuits, benefit from the presence
of dissolved oxygen and therefore require silicon grown from a quartz
crucible.
Several other processes have been investigated for the manufacture
of silicon devices.3 These include sheet-growth methods using dendrites
or fast-growing silicon-crystal forms, growth from dies or free-form crystal-
lization. Casting has also been employed to prepare ingots of large-grain
polycrystalline material. One product of these novel growth methods is
infrared window blanks for various applications. The main thrust of work in
this area, however, has been to prepare low-cost starting material for
photovoltaic applications. Currently, casting of polycrystafline ingots is
the most common way to prepare photovoltaic substrates.
Figure
Figure 1:
1: Photograph
Photograph of
ofaamodern
modern silicon-crystal
silicon-crystal puller
puller (courtesy Kayex-Hamco).
(courtesy Kayex-Hamco).
Silicon Materials Technology 11
HEATER
HEATER
CONNECTOR
ELECTRODE
introducing dislocations into the lower end or tang of the crystal. These
dislocations, if introduced, could propagate upward and destroy crystal
perfection in much of the ingot.
An important improvement in crystal growth was made in the late
1950s by Dash,4 which allowed the production of ingots free of dislocations,
termed zero-D growth. The process, represented schematically in Figure
3, involves special growth conditions during the initial seeding process.
The seed is a single crystal of silicon, usually oriented along a <l OO> or
<l 1 1> direction. Although it is a single crystal, in general it will contain
dislocations or extended disruptions of the lattice. As material isadded to
the seed by freezing, the dislocationswill propagate. By reducing the seed
diameterto mm or half the initial diameter, and making useof the fact that
dislocationsvirtuallyalwaysmakeat leastasmallanglewith respecttothe
vertical axis, the seed can be grown to the point at which all dislocations
havereacheditssurface.Onceadislocationisatthesurface,itis“pinned”,
and substantial energy is required to initiate a new one. Subsequent
growth of the crystal is routinely maintained in the dislocation-free condi-
tion, and all silicon substrates are supplied in this state. The method of
Dash was crucial to the production of ingots of three-inch diameter and
larger, avoiding the tendency of large dislocated crystals to become
polycrystalline. Figure4 shows the stages of crystal growth in a production
puller, while Figure 5 shows a completed ingot ready for further processing.
I
I
L-SEED
I
I I
C, = C,k(l - g)k-’
Electrical Dopants
Boron
Phosphorus
Arsenic
Antimony
Ubiquitous Impurities
Oxygen 1.25
Carbon 0.07
Nitrogen 7 x 10-4
Metals
Iron B x 1O-6
Nickel 2.7 x 1O-6
Copper 4 x 10-4
Gold 6 x 10-6
Aluminum 2 x 10-3
[concentro
06
001
0 01 02 03 04 05 06 07 08 (
Froctlon solldlfled,l
P-TYPE (100)
100mm DIAMETER
FRACTION GROWN
N TYPE (111)
3 INCH DIAMETER
I I I I
0 5 10 15 Xl
=-1
LI
52
;=
4.5
(, ‘..~‘-..~,_:...,‘-...:
A,, 111’
. ..~ . . .. ” . ..-.. ,, _ ,,,_,._ I . “.,_, _.,:______:..,_‘-2
‘-
>2 6
Zi 5
Ei; A,, r2w, ... . . .-, ....> _.
;;‘ 4 .-.,a.. _,“,.._
__.,_./_._
#., ...,... ..A.” _.,..-‘~‘*
_.,.‘_‘ . . .. .;‘: ., .. .: : .;_.. ..
z
o 3
J CLOCHRALSKI GROWN N FWOSPHORUS HIII 3” WAFER
SEED 3 6 9 12 15 18 TAIL
Figure 10: Axial gradient of oxygen in a silicon crystal. Oxygen measured accord-
ing to ASTM F 121430.
Semiconductor Materials
a _ j . - . SEED
0 10 20 30 40
.2 .4 .6 .0 1.0
FRACTION GROWN
Figure 12: Carbon concentration versus fraction of melt solidified for three con-
ditions of crystal growth: circles represent the maximum level of carbon due to
remelt in the original charge, triangles represent crystals grown with standard
polysilicon at atmospheric pressure and squares present data for material grow.n
with standard polysilicon but a reduced pressure of 30 torr.
22 Semiconductor Materials
-I-
PRIMARY f? PRIMARY
-l FLAT
I
FLAT
SECONDARY
FLAT
11 111 n-TYPE
(3 111 P-TYPE
FLAT
( 1 100 ~-TYPE
{
100) p - TYPE
i Blade
i Axial Blade
l4Movement
3
:.
blade
/
3.3 Polishing
Front-surface polishing of the wafer leaves it with a mirror finish
needed for device fabrication. Wafers are mounted on carriers for this
operation, held either by a thin wax layeror byafriction bond to thecarrier.
These carriers are pressed onto a rotating polishing pad of polymeric
material while a polishing slurry isapplied. Generally the processconsists
of two steps on two different polishing machines. In the first step, perhaps
0.001”oronemilofsiliconisremovedfromthesurfaceinaprocesstermed
stock removal. The carrier is then moved to a machine with a smoother pad
surface for final, mirror-finish polishing.
The polishing slurryconsistsof a solution of colloidal silica maintained
atapHof 11 forstockremovalandapHof9forfinishpolish.Thefunctionof
the silica in the polishing process is not understood. It is probable that the
product of the polishing process is silica itself. That is, the surface of the
wafer is oxidized and hydrated to form a silica that can be wiped away by
the pad. It is curious that one must add the reaction product to initiate the
polishing process. The silica acts as some sort of catalyst, and the primary
polishing agent is probably water!15
3.4 Cleaning
Thefinalstepinwaferpreparationisacarefulc1eaningofthesurface.A
sequence of acids and bases is used to remove any contaminants, including
wax residues, if any, and metallic contaminants in the polishing medium.
The cleaning process is assisted by adding hydrogen peroxide to oxidize
these materials, and by heating the liquid baths. The basis for chemical
cleaning has been established by Kern.16
The final bath in chemical cleaning process is usually chosen so that
the wafersurface is hydrophilic, orwater-loving. This means that thewater
wets the surface as the wafer is withdrawn from the bath, and can be
removed uniformly. Acombination of ammonium hydroxide and hydrogen
peroxide is an example of such a bath. On the other hand, a wafer
withdrawn from a hydrofluoric acid bath will be hydrophobic, and waterwill
bead on its surface. This beading can lead to spots on the surface afterthe
water dries. The difference between a hydrophilic and hydrophobic surface
consists of the difference in the native oxide thickness. Silicon will instantly
grow a thin oxide upon emerging from a chemical bath. In the hydrophobic
case, the oxide thickness is 1OA or so, while for the hydrophilic case the
oxide thickness can range from 155OA, depending on bath characteristics.
This can lead to variations in thickness of a subsequent thermal oxide
layer. This variation is important for thin gate-oxide-layer growth in MOS
device fabrication.
Throughout the fabrication process, wafers are inspected for a variety
of parameters including physical dimensions, electrical resistivity, flatness
and surface perfection. Afterfinal inspection, most wafers are packaged in
Next Page
26 Semiconductor Materials
Class 100 clean rooms in such a fashion that the cleanliness will be
maintained until use by the customer.
4. MATERIAL PROPERTIES
Bruce E. Deal
Research Center
Fairchild Semiconductor Corporation
Palo A/to, California
. COMPONENTS IN DEVICES
. CORROSION PROTECTION
. DEVICE ISOLATION
. DOPANT DIFFUSION SOURCE
. GETTER IMPURITIES
. INCREASE BREAKDOWN VOLTAGE
This chapter deals primarily with the thermal oxidation of silicon. The
kinetics of silicon thermal oxidation is first reviewed, with emphasis on the
general oxidation relationship and the thin oxide regime. Important oxide
properties are then summarized. Following is the most important part of
the chapter, a description of process variable-oxidation reaction inter-
dependencies. Some of the variables included are effects on silicon
surface properties (dopant redistribution, charges), surface property
effects on the oxidation process (orientation, doping, cleaning), and
ambient effects (type, chlorine addition, pressure). After a discussion of
oxidation mechanisms, the chapter concludes with a section on other
oxidation processes and an indication of thermal oxidation applications
and trends with respect to future semiconductor devices.
Si 4. O2 -c Si02 I 11
1.0 10.0
OXIDATION TIME (hr)
Figurel: Oxide thickness vs oxidation time for silicon oxidation in dry oxygen
at various temperatures (after Hess and Deal”).
2.0
1.5
1.0
-z 0.5
,i
4n
2 0.2
5
u
z 0.1
8
z 0.05
0
0.02
0.01
0.1 0.2 0.5 1.0 2 5 10 20 30
OXIDATION TIME hr)
02
(a)
(a) Transfer of the oxidant from the gas phase to the oxide
outer surface:
Fl = h (C’ - Co) c 31
= kCi I 53
*3
I 6bl
= Bt 1103
X02
X o = B/A(t+r) 1111
B/A is the linear rate constant.
From Equations 7 and 8 it can be noted that when the oxidation
process is controlled primarily by the parabolic rate constant (at high
temperatures or thick oxides), the kinetics are affected by changes in the
diffusion process or oxidant solubility in the oxide. The latter is proportional
to ambient pressure. On the other hand, at low temperatures or for thin
oxides, where the linear rate constant predominates, the oxidation is also
sensitive to oxidant solubility in the oxide (and ambient pressure) but
depends on those factors affecting the interface rate constants h and k.
These effects will be discussed in more detail later.
Special mention should be made of the correction factor r which is
related to initial oxide thickness x, in Equation 9. It has been noted that for
oxidation in dry O,, an initial thickness region does not appear to be
satisfied by the general relationship Equation 6. Rather, the plot of oxide
thickness versus oxidation time tends to extrapolate through the thickness
52 Semiconductor Materials
axis at about x= 15OA. Thus, in the absence of a model for oxidation in this
region, the practice has been to assign a value of r corresponding to xi=
15OA. Further discussion on the mechanism of thermal oxidationforthese
very thin oxides is presented in the next section.
The thermal oxidation of silicon can be represented by Equation Gfora
wide range of temperatures, oxide thicknesses, orientations, and oxidation
ambients, provided the dependence of the rate constants 6 and B/A as a
function of these variables is known. Values of the rate constants have
been determined by rearranging the general relationship Equation6 intoa
linear expression, plotting x0 vs (t+r)lx, and extracting B as the slope and
-A as the intercept from the resulting plots. Arrhenius expressions of the
form
c s +E’kT 1123
have been used in plotting log B and log B/A vs l/T. Such plots are
presentedin Figures4and5,andvaluesoftheconstantsfortheArrhenius
expression are tabulated in Table 2. These data can be used to determine
any thickness-time relationship for a given set of oxidation conditions.
Similar data are incorporated in the SUPREM program13 and related
computer process modeling programs.
1.0 1 I I 1 1
‘;:
f
P
2
” O.Ol-
O.OOl
0.6 0.7 0.6 09 1.0
1000/T (OK)
1o.c
z
‘E l.O-
s
s
I- Ol-
i
w O.Ol-
5
%
y O.ool-
2
O.cool I I I 1
Figure 5: Dependence of the linear rate constant B/A on temperature for the
thermal oxidation of silicon in pyrogenic steam l-640 Torr) and dry oxygen
(after Deal” 1. Reprinted by permission of the publisher, The Electrochemical
Society, Inc.
PARABOLIC B = C, ,qWkT
(111) SILICON
DRY 02 Cl ~7.72~10~ pm’/hr
Cp 26.23 x lo6 pm/hr
El = 1.23 eV
E2 = 2.0 eV
ax [I31
Ois
dt & + Cl-p
0
In this expression, the first term on the right-hand side is the contribution
from the original linear-parabolic model. The second term incorporating L,
is possibly related to effects of residue left on the silicon surface from the
cleaning treatment (X, 5 15A). The contribution of L, has not yet been
explained. Subsequently, Han and Helms18 have proposed that a mechanism
based on parallel diffusion reactions provide even a better fit to oxidation
data over the entire thickness range.
It is important for future applications of devices having sub-micrometer
feature sizes and film thicknesses in the nanometer range (especially
MOS gate and capacitoroxides), that reaction mechanisms be understood
and characterized for oxides in the very initial stages of formation. The
investigations described are a good step in that direction.
Physical Properties
Formula Si02
f4alecularweight 60.08
Electrical Properties
Energy gap -8 cv
optical Properties
Chemical Properties
GENERAL RELATIONSHIP
TRANSITION
REGION
y INTERFACE
TRAPPED CHARGE, Q,t
\ t
\
\
0 N$
MOBILE IONIC
t
CHARGE Om ++++ +
---- +
i
OXIDE TRAPPED + .FIXED OXIDE
+ CHARQE. Of
CHARGE. Oet
Figure 7: Names and location of charges associated with the thermally oxidized
silicon structure.
58 Semiconductor Materials
Cl51
c/co
I I
0
VG
Figure 8: The determination of fixed oxide charge density Qf and mobile ionic
charge density Q, in thermal silicon dioxide using the MOS capacitance-voltage
technique.
OXIDE SILICON
A. m< 1 0. mc 1
DIFFUSION IN DIFFUSION IN
OXIDE SLOW OXIDE FAST
IBORON) (BORON IN Hz)
C. m >l D. m >l
DIFFUSION IN DIFFUSION IN
OXIDE SLOW OXIDE FAST
~PHOSPflORuS. ARSENIC, (GALLIUM)
ANTIMONY 1
Aspredictedfromthegeneralrelationship,theparabolicrateconstant
is relatively independent of silicon orientation. However, there is some
indicationthatforlowertemperaturesorthinneroxidesthesiliconsubstrate
may cause some oxide structural effect, which in turn might result in an
orientation dependence of the parabolic rate constant.
4.2.2 Dopant Concentration. It was observed some time ago that
areas of semiconductor devices having surface dopant concentrations
greater than 101gcm-3, i.e. emitter regions in bipolar transistors, exhibit
higher oxidation rates than ajoining lightly doped silicon.37*38 Experiments
indicated that the effect for n-type dopants is more pronounced at lower
temperatures or thinner oxides, while for boron doping, an oxidation
increase is noted to some extent over the entire temperature range. These
results implied that high concentrations of phosphorus affect the oxidation
process primarily through the surface reaction rate constant WA, while for
boron both B and B/A contribute to the increase in oxidation rate.
More recently Ho and co-workers3g investigated in more detail the
effect of phosphorus concentration on the thermal oxidation process.
They have attributed this effect to vacancy generation resulting from high
phosphorus concentrations in silicon. These vacancies provide a driving
force for increased interface reaction ratesand relate to the more recently
proposed atomic model of oxidation to be discussed in a later section of
this chapter. Boron tends to segregate into the oxide and will therefore
tend to affect the parabolic rate constant which is more important at higher
temperatures. As boron concentration increases, however, that at the
silicon surface will also increase and therefore raise the interface oxidation
reaction rate. Typical oxide thickness- time data are shown in Figure 10 for
phosphorus.3g Note that the increase of oxidation rate is more pronounced
at lower temperature (800°C) and for thinner oxides, which reflects the
greater contribution of the linear rate constant B/A under these conditions.
’ ’ “I
Figure 10: Oxide thickness vs oxidation time for oxidation time for silicon oxi-
dation in dry oxygen at 800” and 1100°C using (111) silicon substrates doped
with phosphorus up to solid solubility (after Ho et a139).
64 Semiconductor Materials
1400 I I I I 1 I I I I
1200-
3 lOOO-
%
$j 800- NH40H : H202: Hz0
z
: 600-
x
fz 400- DRY 02
1OOO’C
(100) si
200- n-TYPE, 2-8&I cm
0 I I I I I I I I I
0 20 40 60 80 100 120 140 160 180 200
OXIDATION TIME (mid
Figure 11: Effect of pre-oxidation cleaning process on oxide growth rate (after
Schwettmann et ala).
Thermal Oxidation 65
occurs. it follows that Cl, must be the primary chlorine species incorporated
in the oxide and is driven to the Si-SiO, interface by a field in the oxide
during oxidation. If water is added to the ambient, the above chemical
reaction is driven to the left and less chlorine is incorporated in the oxide.
This is supported by the fact that little chlorine is observed in steam
produced oxides. Figure 12 includes an Auger profile of an 11 OOCchiorine-
66 Semiconductor Materials
50
I 1 I I I I
OO 500 loo0 1500 2000
x.J (11
B. CHLORINE CONCENTRATION VS OXIDE THICKNESS
Figure 12: Auger sputter profile (A) and chlorine concentration vs oxide thick-
ness (6) for thermal oxide prepared in 5% HCl/02 ambient at 1100°C using (100)
silicon (after Rouse et a143).
Figure 13: Oxide thickness vs oxidation time for silicon oxidation in pyrogenic
steam (-640 Torr) at 900°C and various pressures (after Razouk et als3). Re-
printed by permission of the publisher, The Electrochemical Society, Inc.
5. OXIDATION MECHANISM
Si+O2--SiO2
02 -DIFFUSION
I
902 I SI
1
Figure 14: Proposed mechanisms occurring at the Si-Si02 interface during sili-
__
con thermal oxidation (after Plummer’“). These figures were originally presented
at the Spring 1981 Meeting of The Electrochemical Society, Inc., held in Minne-
apolis Minnesota.
dt
;if [181
from these various types of analysis have been correlated with those
obtained using improved electrical techniques such as quasistatic C-V,
DLTS, and conductance-voltage measurements. All of these plus actual
device measurements have provided considerable insight into the exact
nature of the Si-SiO, interface. It is not possible here to reference even a
small number of the papers concerning the evaluation of the Si-SiO,
interface; however, some of the more comprehensive reviews on the
subject of surface and interface analysis may be consulted.63-65
The current understanding of the nature of the Si-SiO, interface in
thermallyoxidizedsilicon may be summarizedasfollows. First, it isgenerally
agreed that the transition region between silicon and the bulk oxide is no
more than 1Odor even one or two monolayers. In this region, the composi-
tion changes rapidly from Si to SiO,; and the oxide isapparentlycrystalline
in nature immediately adjacent to the silicon. As a result, the physical,
electrical, and chemical properties of the oxide in this transition region are
markedly different from those of amorphous SiO, and affect the net
properties of oxides up to 200Aor more. There is also a good possibility
that the Si-0 bond angles in the transition region and beyond (up to 50&
are strained, which can also affect oxide properties.
Depending on the oxidation conditions and the silicon orientation, a
limited number of silicon atoms at the silicon surface (as few as one in 1 05)
might not be bonded to oxygen and thus could act as trapping sites (Cl,&
Similarly, some of the silicon atoms on the oxide side of the interface might
be disconnected from adjacent oxygen ions (or certain oxygen atoms
might be missing) and these silicon species could also act as charge or
trapping sites (Cl,). Although these specific defects or trapping sites have
yet to be positively observed, the evidence for their presence is fairly
conclusive as a result of recent investigations.18~1g~66~6gA proposed cross
section structure of the Si-SiO, interface region is presented in Fig. 15,
which includes the possible origin of the four types of oxide charges. This
‘-\
,I-.\ I ’ \
:,
:‘O\ \I
0 ( Na" 0
I I Qf, /’
THERMAL I\ +Nou
S102 -Si - 0 -4, SitO-Si-O-Si-
I I “\I_/’ I I
--LO 0 0 0
I- ”
TRANSITION
I
- ’ -
I
Si - 0 -1
,/I-‘\
Si+ +- 0 -
I
Si-
REQION ‘\’ Qf// I
-L ,
SILICON
-9i - ffi - si - Fji - qi - si - Cfi -
concept of the Si-SiO, interface structure is the most likely to date. As more
sophisticatedanalysisequipment isdeveloped, it is reasonable toassume
that a more accurate description of the Si-SiO, interface will emerge.
Under any consideration, this interface will play a most important role in
future semiconductor devices.
2 34ooc
p”:’ ;
,” 160
120
60
40 1
20 I I I I I I I
0 1 2 3 4 5 6 7 0
TIME (hrs)
Figure 16: Oxide thickness vs oxidation time for silicon oxidation in dry oxy-
gen plasma (30 mTorr, 1 kW, 0.5 MHz) (after Ray and Reisman%). Reprinted by
permission of the publisher, The Electrochemical Society, Inc.
Thermal Oxidation 73
SiO2
Poly-Si
POLYCRYSTALLINE SILICON
502
Si+O2 --SiO2
Si
SiO2
Si3N4
SILICON NITRIDE
B.
SiO2 Si3N4 + 3 02 --r3 SiO2 + 2 Np
Si
i i
SiO;,
To52
Si
reaction. Analysis of the data indicates that the rate determining step is
primarily diffusion of the oxidizing species through the oxide. Values of B,
the parabolic rate constant are almost identical to those obtained for
conventional silicon oxidation, while B/Avalues are much higher.77-7g If no
silicon is present beneath the silicide, the resulting oxides are mixtures of
refractory metal and silicon oxides and are generally not stable or repro-
ducible.
Typical thickness-time data for the thermal oxidation of tantalum
silicide deposited over polycrystalline silicon are shown in Figure 1 8.77
Single crystal silicon oxidation dataare included in thefigurefor comparison.
1.0 I I I I I I I I
0 (100) Si
0.5 A TaSi2 on Poly-Si
7. FUTURE TRENDS
conductors, e.g. GaAs, GaAIP, etc., will be used for various applications,
including integrated circuits. It is reasonably certain, however, that both
silicon and thermal silicon dioxide will continue to play major roles in
semiconductor technology.
REFERENCES
26. P. Balk and J.M. Eldridge, Proc. IEEE 57: 1558-63 (1969).
27. For a discussion of trapping in silicon oxides, see: special issues on device
radiation effects, /EEE Trans. Nucl. Sci., Dec. issues, Vols. NS 21-28,1974-81;
C.T. Sah, /EEE Trans. Nucl. Sci. NS-23: 1563-68 (1976); W.R. Dawes, Jr., G.F.
Derbenwick and B.L. Gregory, IEEE J. Solid-State Circuits SC-1 1: 459-65
(1976).
28. RR. Razouk and B.E. Deal, DARPA Final Technical Report No. NR 322-080,
Contract No. NO001 4-79-C-0297, April 1982.
29. D.R. Young, J. Appl. Phys. 52: 4090-4 (1981).
30. D.J. Bartelink, in: integrated Circuit Process Mode/s (J.D. Meindl and KC.
Saraswat, eds.), Chapt. 15, Englewood Cliffs, NJ: Prentice-Hall, Inc., to be
published.
31. A Goetzberger, E. Klausmann and M.J. Schulz, CRC Crit. Reviews So/id State
Science 6: l-43 (1976).
32. A.S. Grove, B.E. Deal, E.H. Snow and CT. Sah, Solid-State Electronics 8: 145-63
(1965).
33. KH. Zaininger and F.P. Heiman, Solid-State Tech. 13 (5): 49-55 (1970); 13 (6):
46-55 (1970).
34. R.R. Razouk and B.E. Deal, J. Elecfrochem. Sot. 129: 806-810 (1982).
35. A.S. Grove, 0. Leistiko, Jr. and C.T. Sah, J. Appl. Phys. 35: 2695-701 (1964).
36. J.R. Ligenza, J. Phys. Chem. 65: 201 l-l 4 (1961).
37. W.A. Pliskin, ISM J. Rsch. Dev. 10: 198-206 (1966).
38. B.E. Deal, and M. Sklar, J. Electrochem. Sot. 1 12: 430-35 (1965).
39. C.P. Ho and J.D. Plummer, J. Necfrochem. Sot. 125: 665-71 (1978); 126: 1516-
30 (1979).
40. F.N. Schwettmann, K.L. Chaing and W.A. Brown, Paper No. 276 in The Spring
Meeting of The Electrochemical Society, Seattle, WA, May 21-26, 1978.
41. R.J. Kriegler, Y.G. Cheng and D.R. Colton, J. Electrochem. Sot. 119: 388-96
(1972).
42. B.R. Singh and P. Balk, J. Electrochem. Sot. 125: 453-61 (1978).
43. J.W. Rouse, C.R. Helms, B.E. Deal and R.R. Razouk, J. Electrochem. Sot.,
13 1: 887-894 (1984).
44. S.I. Raider, R.A. Gdula and J.R. Petrak, Appl. Phys. Left. 27: 150-52 (1975).
45. T. Ito,S. Hijiya,T. Nozaki, H.Arakawa,M.ShinodaandY. Fukukawa, J. Electrochem.
Sot. 125: 448-52 (1978).
46. T. Ito, I. Kato,T. Nozaki, T. Nakamuraand H. Ishikawa,App/. Phys. Lett.38:370-72
(1981).
47. T. Ito, T. Nozaki and H. Ishikawa, J. Electrochem. Sot. 127: 2053-57 (1980).
48. S.S. Wong,C.G. Sodini,T.W. Ekstedt, H.R.Grinolds, K.H. Jackson.S.H. Kwanand
W.G. Oldham, J. Electrochem. Sot. 130: 1 139-44 (1983).
49. P.T. Panousis and M. Schneider, Paper No. 53 presented at The Spring Meeting
of The Electrochemical Society, Chicago, IL, May 13-18, 1973.
50. R. Champagne and M. Toole, So/id State Tech. 20(12): 61-63 (1977).
51. N. Tsubouchi, H. Miyoshi, A. Nishimoto and H. Abe, Jap. J. Appl. Phys. 16: 855-56
(1977).
52. R.J. Zeto, N.O. Korolkoff and S. Marshall, SolidState Tech. 22(7): 62-69 (1979).
53. R.R. Razouk, L.N. Lie and B.E. Deal, J. Electrochem. Sot. 128: 2214-20 (1981).
54. L.N. Lie, R.R. Razouk and B.E. Deal, J. Electrochem. Sot. 129: 2828-34 (1982).
55. J.D. Plummer, in Semiconductor Silicon 7987 (H.R. Huff, R.V. Kriegler and Y.
Takeishi, eds.) pp 445-54, The Electrochemical Society, Pennington, NJ
(1981).
56. R.B. Fair, J. Electrochem. Sot. 128: 1360-68 (1981).
57. A. Lit-t, D.A. Antoniadis and R.W. Dutton, Appl. Phys. Left. 35: 799-801 (1979);
J. Electrochem. Sot. 128: 1 131-37 (1981).
Thermal Oxidation 79
Kenneth E. Bean
Texas Instruments Incorporated
Dallas, Texas
INTRODUCTION
. SYNTHESIS
. THIN FILMS
EPITAXY
POLY
. OXIOES
. NITRIDES
. CARBIDES
. SILICIDES
80
Chemical Vapor Deposition 81
. SILICON SYNTHESIS
Sic4 + Hz
SiHCI3 + Hz
CVDlFLUlD BED
SiH4 + Hz I
the process by which the technology is formed, and the uses of this
technology in silicon manufacturing. In the synthesis of ultra-high purity
elemental silicon for today’s semiconductor manufacturing, we may begin
the process with the hydrogen reduction of an ultra-high purity silicon
halide, such as (SiCI,) or (SiHClJ1-4 by CVD of the elemental silicon, on a
high purity silicon rod such as that shown in (a) of Figure 1. This high
temperature reduction takes place in an all quartz system under very
precisely controlled high purity gas flow conditions. When this CVD reaction
has reached completion we will have obtained a polycrystalline rod similar
to the section shown in (b) of Figure 1. This high purity polycrystalline
silicon material is then broken into small pieces, placed in a high purity
quartz liner or crucible which is then heated to the melting point of silicon,
1420”C.Afterthethermalstabilityofthemoltensiliconpoolisestablished,
a carefully oriented seed, cut from single crystal silicon of the desired
crystal orientation, is dipped into this molten silicon and then slowly
rotated and withdrawn to grow the single crystal of the desired diameter.5
This melt is carefully doped to provide the desired conductivity type and
resistivityfortheslicesorsubstrates. Figure lcshowsthetop,orseedend,
of a single crystal of silicon and Figure 1 d shows a sawed slice from such a
crystal. The standard diameter of the silicon slice used by most silicon
device or integrated circuit manufacturers today is 125 mm plus or minus
25 mm while 200 mm is being developed. After the crystal is grown and
sliced by the use of diamond saws the slices are ground, lapped and
chemically/mechanically (chem/mech) polished to remove all surface
damage introduced by the sawing, lapping and polishing operation.
EPITAXIAL DEPOSITION
Table 4: Epitaxy
.The sil icon epitaxial layer, or layers, is usually the only active
semiconductor material in the device or circuit.
Epitaxial Film
Substrate
Figure 2: Epitaxy.
J OXIDE PATTERN
- POLISHING DAMAGE LAYER
- SINGLE CRYSTAL SUBSTRATE
(111)3.5” OFF, (100) 0’ OFF.
EPITAXIAL FILM
OUF
SUBSTRATE
In the early days of silicon epitaxy most of the work or efforts were
aimed towards the deposition of thin films on (1 11)” silicon substrates by
the hydrogen reduction of silicon tetrachloride (SICI,), or the silicon tetra-
bromide (SiBr,). In other attempts silicon tetraiodide (Sil,) was also used.
Table 5 lists the silicon bearing halides and silicon hydrides in the order of
use historically and also in the order of descending energy or temperature
required for the reduction. In production today most people use silicon
dichlorosilane (SiH,CI,). This material readily decomposes at about 1050°C.
Silicon hydride(SiH,) decomposes at an even lowertemperature. However
there are problems in the epitaxial deposition of thick films ,using (SiH,)
dueto thermal decomposition in thevapor phase. When thisoccursabove
the epitaxial substrate particles form in the gas phase which fall on the
substrate resulting in the formation of spurious nucleation sites. It should
also be noted that trichlorosilane(SiHClJ, tribromosilane(SiHBr,), Silane
86 Semiconductor Materials
l Silicon Source-Halides-Hydride
l SiCL, -Early
SiBr4
Sit4
l SiHC13
SiHBr3
l SiH4
l SiHzClz
l Reduction Source-H2
l Vapor Etching-HCI
l Sic14
l SF6
. Hz0
’ Hz
Shortlyafterthebeginningofattemptstodosiliconepitaxyintheearly
1960s it was learned that the substrate was of prime importance and must
be extremely clean and free of defects at the beginning of epitaxial
nucleation.Allattemptstoclean thesubstratepriortoplacing thesubstratein
the epitaxial reactor met with high density defect levels in the epitaxial
material. Due to this problem the in situ HCI chemical vapor etching (CVE)
process was developed in 1963. In this process the substrates were
initially cleaned and then placed in the epitaxial reactor. The final cleanup,
removal of contaminants and crystallographic defects, is then done at
approximately 1150-l 250°C by high temperature in-situ etching, using
HCI as the etchant.‘*-l3 Other etches which have been experimentally
used are listed in Table 5. However, the development of in-situ HCI vapor
etching was a key development in the history of silicon epitaxial technology
and allowed the production of large volume, low defect epitaxial material.
With the development of today’s better cleaning and polishing processes,
Chemical Vapor Deposition 87
DEPOSITION
12
11
2
1
0 '______---_--_-_t-----\--- \
-1 ” ” ” ” ’ ” ’ “I’
4 2 4 6 9 10 12 14 16 18 20 22 24 26 29 30 32
ETCHING MOLEXSiHALlDE
Figure 4: Deposition and etch rate vs mol % silicon halide. The deposition rate
is also affected by reactor design. Curves A and B are from a vertical reactor.
Curve C is from a multiple slice vertical reactor in which each slice rotates on its
own susceptor. Curve Cl is data from Henry Therur of Bell Labs using a single
slice vertical reactor.
MOS PROCESS BIPOLAR PROCESS
REGROWN
;
L
DRAIN
fiREFLOW
J
1) SELF ALIGNEO.SOURCE-
OXIOE
1
I
!izz= N*
J
OXIDE
/ el BASE DIFFUSIDN
3
ibS
g) PHOS OXIDE
0 EMITTER DIFFUSION
cl OXIDATION
c) N.EPI GRDWTH
I
0) CONTACT OR
L-
d) ISOLATION DIFFUSION
a
a) ACTIVE AREA DEFINITION
ilW I
hl METALLIZATIDN
EMISSIVITY UNCORRECTED
0.048 pm
1 t I I
1000 1050 1100 1150 1200 1250
TEMPERATURE. “C
POLISHING
etchrateasafunctionoftemperatureforanyratiodesired.Forexample2%
HCI at approximately 1220°C is an ideal etching condition in order to
remove the impurities and damaged region of the substrate surface prior
to epitaxial deposition. We are now ready to perform theepitaxial deposition.
Silicon bearing halides, or the hydride, may be used in orderto perform our
epitaxial deposition. These materials are shown in Figure 8 with silicon
tetrachloride being the silicon source requiring the highest temperature
for deposition, trichlorosilane (SiHCIJ requiring a medium deposition
temperature and dichlorosilane (SiH,CI,) or the silicon hydride (SiHJ
requiring the lowest temperature for deposition. In the deposition reaction,
the deposition rate increases as a function of temperature. The cross-
hatched bar zone shown in Figure8 delineates the change-over point from
a kinetically controlled reaction to a diffusion controlled reaction. As in the
case of HCI vapor etching, one should operate in the flat portion of the
curve forwell controlled epitaxial deposition. The importance of hydrogen
and chlorine, in the reaction cannot be overlooked. Silicon tetrachloride
decomposes in the presence of hydrogen and thermal energy from the
heated susceptor. This material (SiClJ does not thermally decompose as
does the other halides(SiHCI,) or(SiH,CI,) orthe hydride(SiH,), without the
presence of hydrogen as a reducing agent. The halides and the hydride all
have a hydrogen atom attached to the molecule, while(SiCIJ does not. The
presence of chlorine in the silicon halide is of importance in preferential
deposition.14-l5 The by-product chlorine and/or HCI gives the silicon atom
a vehicle in which it may be transported from one nucleation site to
another, which may be a more desirable crystallographic site. The by-
product of silicon hydride, does not contain chlorine making this material
somewhat lessdesirablefor preferential deposition unless chlorine is also
added to the reaction.14’16 The effect of silicon atomic mobility ortransport
will be discussed in a later section.
"C
1300 1100 900 900 700 600
GETTERING
Intheeverincreasingpackingdensityfortoday’stechnologiesofVLSl
and ULSI circuitry, silicon epitaxy plays a very important role.‘O MOS
circuits built in epitaxial material are in general superior in performance to
those built in bulk silicon material. In the past two years considerable
interest has been placed on improving the quality of epitaxial material for
MOS devices. One of the recent developments which has improved the
quality of epitaxial material has been the denuding of defects from the
substrate surface area prior to epitaxial deposition. This denuding is
accomplished byathermal processwhich in mostcases,consistsof ahigh
temperature cycle (approximately 1 1OOC), followed by a lower or inter-
mediate temperature cycle (approximately 650°C) then followed by a
high temperature (approximately 1OOPC) thermal cycle. This denuding
process is highly dependent upon the oxygen and/orcarbon concentration
in the original bulk silicon substrate. The initial high temperature cycle
provides sufficient energy to dissolve the oxygen precipitates that are
present and to cause out diffusion of the oxygen from the surface areas of
the substrate. The intermediate thermal cycle provides energy for renuclea-
tion and growth to stability of oxygen precipitates in the center, non out-
diffused, region of theslice.Thefollowing high temperaturecycle provides
energyforthese precipitates to grow and getter oxygen, heavy metals, and
other defects, during device processing, to these precipitate sites. See
Figures 9 and 10.
Figure 9 shows a cross-sectional view of a denuded substrate after
epitaxial deposition. The high bulk defect density is evident as is the
denuded zone just below the epitaxial silicon interface. In this figure there
are defects at the interface between the epitaxial film and the original
substrate. This indicates that no HCI vapor etching was carried out in this
process. It also shows that the original substrate surface acts as a trap to
hold these defects, which may be bulkstackingfaults, oxygen precipitates,
heavy metals, or carbon atoms. Figure 9 is a photograph of a cleaved
silicon slice. No polishing or potting was required prior to the Wright
Jenkinsetchinordertobringoutthedefectsinthedenudedareaaswellas
the epitaxial film. It is also interesting to note that the defects at the
substrate/epitaxy interface do not cause defects in the epitaxial film at a 1:l
ratio. It appears that only approximately one out of 10 defects at the
interface actually cause stacking faults or defects in the epitaxial film.
Another method of slice/wafer processing that improves device per-
formance is that of backside gettering prior to epitaxial deposition. There
are several methods of backside gettering including sandblasting or
abrading of the back surface of the substrate to produce defects or traps,
and oxide or nitride films on the backside surface to produce strain fields.
However the most successful appears to be that of depositing a thin CVD
film, approximately 1 micron thick, of polycrystalline silicon across the
back surface of the epitaxial substrate prior to final polishing of the front
side. The high density of defects produced by the grain boundaries and
dislocations, due to the polycrystalline film, provide avery effectivegetter-
ing mechanism. During device processing heavy metals may be gettered
92 Semiconductor Materials
u”
,m
E
.-
Chemical Vapor Deposition 93
Figure 10: (a) 90° cleaved cross section, epitaxy on denuded substrate with back-
side-gettering poly Si layer. 198X P.C. 5 minute W.J. etch. (b) (c) Misfit disloca-
tion extrinsic gettering.
94 Semiconductor Materials
all the way from the front surface to the back surface where they are
trapped at these defect (gettering) sites.
Figure 10a shows a cross-section of a cleaved slice which has a
gettering region at the back surface, a denuded zone just above the
backside gettering media, the high density bulk defect area of the slice,
and then at the top surface a denuded zone just below the epitaxial-
substrate interface. Again, a row of defects is noticable at the epi substrate
interface indicating that this substrate had insufficient or no HCI vapor
etching, in situ, prior to the epitaxial deposition. This also shows that one
could build in intrinsic gettering at desired positions immediately below
theactivedeviceregionofthesemiconductorcircuit.Thistypeofgettering
immediately below the active surface area is very effective18and can also
be designed in discretionarily to provide gettering only at the desired
circuit areas. Figure 1Oa is a cleaved cross-sectional view which required
no polishing or potting prior to the Wright-Jenkins etch to delineate the
defect region, thedenudedzone,andtheepitaxiallayeraswellasthe back
side gettering polycrystalline silicon film.
Asmentionedabove,in-situHCIvaporetchingwillremovethedamage
sites and/or surface traps prior to epitaxy. However, one may wish to leave
or form a new damage layer, for low temperature processing intrinsic
gettering,innearproximitytotheactivedeviceregionofthestructure.Ifso,
a film can be deposited between the substrate and the epitaxial film which
is intentionally doped with, for example germanium, to produce a built-in
misfit dislocation strain field. 10.17,18Figure 10B shows a cleaved cross-
sectional view of a single layer misfit dislocation, extrinsic gettering strain
field, and a single layer epitaxial film. Figure 1Oc shows experimental
multiple layers of strain field/epitaxial silicon films with increasing Ge
doping in the strain fields as they were deposited. This increase in Ge
doping causes a noticeable increase in the density of misfit dislocations
within the strain field layers. Wright-Jenkins etch was used to reveal these
damage sites.
SELECTIVE DEPOSITION
the oxide pinhole. Note also in this figure that the surrounding area is
completely free of polycrystalline silicon deposition or nucleation on the
oxide. This indicates that the silicon atoms, above the oxide area, moved to
a preferred site in the open silicon area. We can take advantage of this
atomic mobility for preferential deposition of silicon at desired sites on an
otherwise oxide or nitride covered masksubstrate. Such a mask is shown
in Figure 12. Shown in this figure are five micron diameter circles, on 25
micron centers, opened up through an oxide or a nitride mask on a silicon
substrate. The goal in this experiment is to preferentially deposit epitaxial
silicon in the five micron diameter circles but have no poly silicon nucleate
on the oxide. An epitaxial diode will be formed at the interface of the
original p substrate with the n epi deposit. After the epitaxial growth has
proceeded up through the mask, in this case oxide, lateral spreading
occurs over the oxide to form a large area “epi top” for electron beam
charging. This process provides a very small area p-n junction diode with
very low parasitic capacitance but with a large “epi top” beam collection
area for the production of Vidicon type detectors. In this process the
preferential deposition must be very complete in that single crystal silicon
is nucleated in the open areas. No spurious nucleation of polycrystalline
silicon can be tolerated on the oxide which would bridge across two diode
“epi tops” thus causing a defect in the array. The diode density in this array
is one million diodes per square inch.
Figure 13 is a top view of such an array after preferential deposition.
The (100) structure is clearly evident in the epitaxial “epi tops”. “Epi tops”
deposited on (111) substrates show an equilateral triangle structure
whereas the (100) substrategives the perfect square‘lepi top”orientation.
Figures 14a and 14b show two SEM photographs with 14a being a low
angle SEM of the cross-sectioned substrate/epi structure and 14b being a
near 90 degree cross-section after the oxide had been etched away. In
14a, the original diode area can be seen as well as the mask oxide which
has been brokenawaywiththecleavageoftheslice. Inthecross-sectionat
14b the original diode structure or size can be seen at the substrate
interface. The lateral spreading, in all directions, over the oxide is approxi-
mately equal to the diameter of the original diode, thus a 3X increase in
diameter and >9X increase in area.
Other examples of preferential deposition making use of the atomic
mobility of the silicon atom are shown in Figures 15 through 7 8. Figure 15
shows the preferential epitaxial deposition of silicon in the vertical lines
across the bottom portion of theslice. In thisgrating there are three micron
wide lines of oxide with two micron wide areas of open silicon between
them. Single crystal silicon is nucleated in the open silicon areas and
grows up through the oxide and then spreads laterally as shown in Figure
18. Also shown in Figure 15 in the top portion of the photograph is an area
with continuous oxide mask with polycrystalline silicon nucleated only on
the top half of this oxide area. The lower half of the oxide area is completely
clean and free of spurious nucleation of polycrystalline silicon. The silicon
atoms have enough mobility, to move to preferred sites in the open silicon
area or to deposit out (at super saturation) as polycrystalline silicon over
the oxide. In this experiment, the deposition temperature was 1 150°C and
Next Page
Chemical
Vapor Deposition
97
Figure 12: Tivicon mask-5µm open circles on 25 µm centers.
Chemical Etching and Slice
Cleanup of Silicon
Kenneth E. Bean
Texas instruments Incorporated
Dallas, Texas
INTRODUCTION
Slice cleaning and wet chemical etching have been key semiconductor
processing technologiessince the beginning ofsemiconductorfabrication in
the late 1940’s and early 1950’s. The demand for cleanliness, control of
purity: and freedom from defects becomes more stringent with each
advance in device and circuit complexity. This chapter discusses wet
chemical etching of silicon from the standpoint of planar etching, orienta-
tion dependent etching (ODE), concentration dependent etching, and
defect delineation etching. It also discusses the etch composition, the
masking materials used for preferential etching, mask alignment, and
applications for the above etching technologies. We will also discuss
silicon slice cleanup procedures and effects thereof. Table 1 summarizes
the subjects to be discussed, in order of discussion.
Planar etch is a solution that etches silicon in all crystallographic
directions at the same rate. A common formulation is made up by mixing
hydrofluoricacid(HF)=8% byvolume,nitricacid(HNO,)=75%,andacetic
acid (C,H,O,) = 17%. At 25°C this solution etches silicon slices(wafers) at
approximately 5 pm per minute, (see Table 2). Orientation dependent
etches (ODE) have been developed which etch much faster in one crystallo-
graphic direction than in another. For example, a solution of potassium
hydroxideandwater(KOH+H,O) inequalparts(50%-50%weight)at80”C
etchessiliconinthe<l lO>direction=700timesfasterthaninthe<l 1 l>
126
Chemical Etching and Slice Cleanup 127
. (lll)TRACE-FLAT SECCO
. ETCH SOLUTIONS . (100) WRIGHT-JENKINS
. ALIGNMENT
1 SCHIMMEL
. APPLICATIONS
SIRTL
. (111) DASH
LEO’S
. STANDARD CLEANUPS
. CHOLINE CLEANUPS
(1101 ODE ETCHES [1101 600 X I1111 KOH - ii20 ~0.8 &MlN AT 80°C
5-I
OIRECTION SO 50 VOL. IN 11101 SILICON
direction. SeeTable 2 and Figure 28. If we add normal propanol to the KOH
and H,O etch, we can also etch silicon in the <lOO> direction approxi-
mately 100 times faster than in the <l 1 1 > direction, at 80°C. See Table 2
and Figure 10.
If we mix the same mineral acids used in the planar etch solution in the
ratios of one part HF, three parts HNO,, and 10 parts C,H,O,, we have an
etch commonly known as Dash etch, which etches p+ silicon or n+ silicon,
>7X101g carrier concentration, much faster than p- or n- silicon. In
contrast, the KOH-propanol-water etch “stops”(slows down by- 20X) at a
p+ interface. The ethylenediamine (EDA) etch, made up of EDA, pyrocate-
chol and water (see Table 2) is also both orientation dependent and
128 Semiconductor Materials
@ (110)
@ (221)
0 (111)
@ (334)
0 (112)
@ (114)
@ (100)
the (100) surface in the pattern shown. These (111) planes intersect the
(100) surface plane at an angel of 54.74 degrees. If we lookat the cleavage
plane or the edge of the cleaved border of the slice we see a crystallo-
graphic plane that is inclined to the surface at this 54.74 degree angle. In
the case of the (111) slice there are three cleavage planes 120 degrees
apart. When they extend all the way across the slice, they break into pie-
shaped segments with 60 degree angles. The (111) cleavage planes
intersect the (1 11) surface of the slice at 70.53 degrees. If one again
cleaves the sections that have already been cleaved, they will continue to
cleave,inthecaseof(lOO)intosquaresorrectangles,andinthecaseofthe
(1 11) into 60 degree triangular shapes. If we cleave a(1 IO) silicon slice it
will cleave 90 degrees to the (110) surface. If we continue to cleave these
sections we will see that they form rhombic shapes as shown in Figure 6.
However, ifwelookattheedgeorthecleavedsurfacewewillseethat inall
casestheyare90degreestothe(l lO)surface.Thesearethe(lll)planes.
In today’s silicon semiconductor processing there is great interest in
MOS-type structures. For this type of structure the (100) slice orientation
is usually used due to the low surface state density at the (100) silicon
surface-silicon dioxide interface.
Figure 7 is a stereographic projection of the standard (001) or (100)
face centered cubic crystal structure. In this projection we are looking
directly at the (100) surface as we would in a (100) silicon slice. Note that
the four (1 11) planes which intersect the( 100) surface are slightly greater
than half way out to the periphery of the projection (slice) in <l lO>
directions. In otherwords, there isa(l 10) planeperpendiculartothe(l00)
plane and tangent to the periphery at this point. The (111) planes are
actually coming into the (100) surface at angles of 54.74 degrees. The
(1 11) planes are also at 90 degree angles to each other. Those planes
designated by the Miller indices at the periphery of the projection are
known as directions. Starting at the bottom, orthe periphery closest tothe
observer, is a (100) plane, indicating that this is a <l OO> direction. Those
Miller indices indicating planes between this (100) plane and the (100)
plane at the center of the projection are lying in this <lOO> direction.
Moving to the right of the bottom center we have a <310> direction.
Moving up along the periphery we find a <210> direction, a <320>
direction, and then the <l 1 O> direction, in which the (11 1) predominant
plane lies. Further examination of this (100) projection shows a four-fold
symmetry in this (100) plane. All four quadrants are exactly alike, and the
(1 1 1) planes are 90 degrees to each other, as are all other families in this
projection. Note also that in this (100) projection are a <31 O> direction to
the right of the <lOO> direction at the bottom of the projection and a
<310> direction to the left of the <l OO> direction at the bottom of the
projection. These same two (310) planes are also to the right and to the left
of the <lOO> direction at the top of the projection. Likewise, they are
above and below the <l OO> direction at the right and at the left of this
projection. This shows there is double four-fold symmetry of <310>
direction planes in the (100) projection. In this double four-fold symmetry
of <310> directions the predominant planes are the (31 1) and the (331).
These planes etch and deposit rapidly in processing. The effect of these
134 Semiconductor Materials
Chemical Etching and Slice Cleanup 135
136 Semiconductor Materials
Ii101 FLAT
(110)FLAT
the traces of the two (1 1 1) sidewall planes meet at the bottom of the etch
moat the orientation dependent etch stops etching. Toverify this statement,
some (100) silicon slices were etched using a mask designed to etch 50
micronsdeepusingthe(l00)ODE.Thismaskalsohascornercompensation
designed to give square or right-angle corners at a depth of 50 microns.
Figure 17a shows a top-focused view of one of these slices, which was
etched for 38 minutes at 80°C. The etch depth is 38 microns. In the left-
hand picture the focus is at the top of the (100) slice. In the right-side of
Figure 17a the focus is at the bottom of the etch moat. Note that the bottom
is flat at the etch front and is a (100) plane, since we have only etched 38
microns deep. Careful measurement of the etch width at the silicon
surface/oxide mask interface shows that the etch moat is 97.5 microns
wide with no measurable undercutting. Also, note that in Figure 17a, at the
lefttopfocusviewthereisaslighttipofsiliconstickingoutatthecornersin
allfourquadrantsoftheetchmoat.Rememberthecornercompensationis
designed for 50 microns etch depth and we have etched only38 microns
deep.Thesliceswerethenplaced backintheODEsolutionandetchedfor
an additional 30 minutes, for a total etch time of 68 minutes. We have now
over-etched by 18 minutes for the corner compensation design. Careful
measurement of the top of the etch moat again show that the etch moat
width at the oxide is 97.5 microns, indicating that there is no undercutting
oretchingafterthetraceofthe(111) planesreachestheedgeoftheoxide
mask. The oxide corner compensation of the mask is clearly visible in both
the top focus and the bottom focus views of Figure 17b. The bottom focus
shows a completely v’d-out etch front with no (100) remaining. To further
prove this statement, the slices were again placed in the ODE solution for
an additional 30 minutes, making a total of 48 minutes over-etching forthe
Chemical Etching and Slice Cleanup 139
Figure 10: Top view and cross-sectional view of ODE etched (100) silicon.
140 Semiconductor Materials
Chemical Etching and Slice Cleanup 141
Figure 14: Corner faceting due to (100) ODE etching using different etches.
Next Page
W.C. Dautremont-Smith
Richard A. Gottscho
R. J. Schutz
AT & T Bell Labora tories
Murray Hill, New Jersey
1. INTRODUCTION
191
192 Semiconductor Materials
2. FUNDAMENTALASPECTS
l I Td = lo-” V cm-z
Plasma Processing 195
plasma. The electrons rapidly drift toward the positive electrode leaving
behind a net positive restoring force, which prohibits further electron
depletion. Asteadystate is achieved when the plasma potential is sufficiently
positive that electron and ion loss rates become equal. The time it takes
electron plasma frequency,17vlg
we = (ne2/tome)lA, (1)
6, = ne*/m,v, (34
IP
RP = -
%AP)
196 Semiconductor Materials
ct?
k
P
Plasma Processing 197
IO' \
\
\\ \\
\ \ \\
106 \ \ \ \
1 1P
\ XP’WCP N WCo
- *P
\ ‘\J
IO"
3 104
N
IO3
IO2
JP meve Pp
RP’(T -( -I--
e P nee2 *P
IO' I I I I
10-3 10-2 10-i 1 10 IO2
w/2-rr (MHz)
Figure 3: Plasma impedances vs. frequency. Resistances, Rs and Rp, are deter-
mined by ion and electron conductivities, Ui = 5.5 X lo1 mhos cm-’ and ue =
lo* mhos cm-’ , sheath and plasma thicknesses, Is = 0.25 cm and Ip = 1 .O cm,
respectively, and the electrode area, AS = Ap = 45.6 cm2. The capacitive impe-
dances, XS and Xp, are determined by the thicknesses, areas, and frequency.
From Figure 3, we see for frequencies below 100 MHz, the plasma body
impedance is predominantly resistive. However, the sheath impedance
changes around the ion plasma frequency (oi = 1 MHz), from being pri-
marily resistive at lower frequencies to being primarily reactive at higher
frequencies. In other words, the sheath capacitor becomes a current
shunt.Abovew, the ionscan no longer respond to the instantaneousvalue
of thefield(see below) and sodisplacement instead of conduction current
dominates.
2.1.2.7. Experimental Verification. Recently, experimental diagnostic
techniques have been developed which allow the concentrations of free
radicals and ions as well as electric field amplitudes to be measured in situ
and non-intrusively.21~22~30-47F or a recent review, see Reference 37. These
techniques allow us to see the extent to which theequivalent circuit model
is appropriate.
2.7.2.7.7 Voltage distribution. How does the electric field vary across
the electrode gap? Acc.ording to the equivalent circuit model, we expect
the field to be largest in the sheaths over the frequency range of interest to
plasma processing (seeTable 1 and Figure 3). In situ electric field measure-
ments are consistent with the model. The local field is plotted as a function
of position for rf discharges through BCI, in Figure 5. For the range of
frequencies studied, 50 kHz to 14 MHz, the field is always greatest in the
sheaths.*’
2.7.2.7.2 Diode behavior. How does the local electric field vary with
time in the electrode sheaths? This has been measured by spectrally
resolving laser-inducedfluorescencefrom parity mixed rotational levels of
the BCI radical formed in rf discharges through BC13,21~46~46 The different
parity levels are mixed by the local electric field; the extent of mixing is
dictated by the field strength as well as the excited state dipole moment
and zero-field energy level splitting. 46.48 Parity mixing is detected by
recording the intensities of transitions which would be “forbidden” in the
absence of an electric field and whose line intensity is a direct measure of
the electric field amplitude. The technique is illustrated in Figure4, where
the field has been sampled at two different times during the rf cycle by
firing the laser synchronously with the applied rf.** Note that the change
signals for the “forbidden” and “allowed” lines in Figure 4 are equal and
opposite insign becausethe”forbidden”component has borrowed intensity
from the“allowed” component asa result of the field-induced mixing. Both
measurementsare made one mm from the powered electrodesheath. The
upper trace is obtained at a time when the powered electrode is the
momentary anode (applied voltage a maximum); the“forbidden” line in the
center is weak compared to the “allowed” lines on either side. Thus, the
field is small during this part of the cycle. However, when the laser is fired
Plasma Processing 199
Q(5) R(4)
ANODE
CATHODE
1 I I
2721 2720 2719
WAVELENGTH 6,
during the cathodic part of the cycle (lower trace), the “forbidden” line is
comparable in amplitude to the “allowed” lines. During this part of the
cycle, the electric field is strong. Thus, the applied field is rectified in the
sheaths. As a result of the difference in electron and ion mobilities the
plasma potential is”tied” to the anode potential. This behavior is accounted
for in the equivalent circuit (Figure 2) by placing diodes in parallel with the
sheath resistors and capacitors.
2.7.2.7.3 Frequency response. As the frequency is varied an interesting
transition is seen to occur in Figure 5. The local field decreases by roughly
a factor of two above 5 MHz. The peak voltage which must be applied
across the plates in order to maintain constant power also decreases by
this amount (Figure 6a). From the equivalent circuit model we see that this
transition corresponds to a change in the sheath impedance from being
predominantly resistive below wit0 predominantlycapacitive above wr The
total impedance decreases above 5 MHz as the sheath capacitive impe-
dance becomes less than the ion resistive impedance (Fisure _ 3). At
constant power the ratio of voltages at high and low frequency is! gicren by,
2
2 4
POSITION (mm)
(a)
-
0
E I I IIrlrli I I1111111 I I Ilrllll Illll~
3 . PIE
& 6 0 LIF (b)
2
g 5
S
FREQUENCY(HZ)
Figure 6: (a) Peak voltage for discharge through BC13 as a function of frequency
at a pressure of 0.3 Torr and a power density of 0.13W cm3. (b) BCI radical densi-
ties as a function of frequency. PIE refers to plasma-induced emission, i.e. ex-
cited state radicals, while LIF refers to laser-induced fluorescence, i.e. ground
state radicals (from Reference 21).
202 Semiconductor Materials
(5)
which is in very good agreement with the values in Figures 5 and 6a. This
transition is also apparent when one examines the voltage and current
waveforms in Figure 7.21 Not only does the current increase and the
voltage decrease above 5 MHz but the phase shift between the two
increases toward 90” as the capacitive component becomes an important
sheath current shunt.
2.7.2.7.4 Power dissipation. From the above discussions of voltage
distribution and frequency effects we can see how and where power is
dissipated in rf dischargesasafunction offrequency.At lowfrequency, the
sheaths are primarily resistive and the sheath resistance is much greater
than the plasma resistance. Thus, we expect power dissipation to occur
primarily in the sheaths. Ions accelerated by the sheath field can dissipate
their energy in basically two ways. Collisions with neutrals can result in
ionization, excitation, chemical reactions and/or heating; collisions with
electrodes can result in any combination of surface damage, sputtering,
secondary electron emission and/or heating with either implantation or
reflection of the incident ion.
The response of ions to the sheath field and the consequences of this
response at low frequency can be seen very clearly in Figures 8 and 9.*’ In
Figure 8, the density of Cl,+ measured in the sheath by laser-induced
fluorescence is plotted as a function of time. During the positive part of the
cycle when the local field is very small (Figure 4), the ion concentration
builds as a result of both diffusion of ions from the p.lasma into the sheath
and ionization by electron impact. During the negative part of thecycle, the
ion concentration decreases precipitously as a result of the large cathodic
fields (Figure 4) which sweep the ions out of the sheath toward the
electrode. The extraction of high energy ions at low frequency causes
ionization, excitation, and secondary emission of electrons. This can be
seen in Figure 9 where the uv emission intensity from BCI radicals is
recorded as a function of position across the electrode gap at different
timesduringtherfcycle.*‘Whentheionsareextractedduringthecathodic
cycle the emission is brightest because the ions and secondary electrons
collide with neutrals and produce a cascade of ionization, excitation
(Figure 9), and dissociation as they are accelerated across the sheath.
This periodic build-upand high-energy extraction of ions in thesheath
has no dc or high frequency analog. In dc anode sheaths, ions build up to
some steady-state level but are never extracted with high energy. In dc
cathode sheaths, the ion concentration can never build up to a large value
owing to the large, extracting sheath field. While the time-averaged flux
may be similar in the dc and low frequency rf discharges, the pulsed ion
bombardment in the latter may make a difference in heterogeneous
reaction rates. In high frequency(i.e. above or) discharges, the ions respond
only to the average field, which is less than at lower frequencies owing to
the resistive to capacitive transition discussed above. The net result at
203
0.2 C-
0.1
-0.4
-0.2
0.4
0.2
9 -0.2
AZ
3 -0.4
2 0.4
5
(3 0.2
-0.2
-0.4
0.4
0.2
- 0.2
-0.4
Figure 7: Current and voltage waveforms for discharges through ECIS at several
different frequencies and a pressure of 0.3 Torr and a power density of 0.13 W
cmm3(from Reference 21). One unit of n corresponds to a 1/2u, where v is the rf
frequency in Hz.
204 Semiconductor Materials
I ” i I”’ I’
0 1 2 3 4
TIME ( UNITS OF 7r 1
Figure 8: Clz+ ion density vs. time in the sheath of a 55 kHz discharge through
Cl? at 0.3 Torr and 0.6 W cm3. One unit of 71corresponds to 18.2 psec.
high frequency is that the ions experience a smaller extraction force and
again build up to some steady state concentration.
The response of ions to the instantaneous field and the change in the
amplitude of this field with frequency at constant power affect the energy
with which ions impact electrode or device surfaces. This is evident in the
ion energy distributions measured by Bruce4g as a function of frequency
(Figure 10). At low frequency, the ions are accelerated to the full sheath
potential, which is approximately the full applied potential (see Figure 5)
on every half cycle as they traverse the sheath. To the extent that there is
ionization and energy loss in the sheath, the ion energy distribution will be
skewed toward lower energies. At 40 Pa(0.30 Torr) of Cl, and 100 kHz, the
Cl+ and Cl,+ ion energies on average are significantly less than the full
sheathpotentialbutthemaximumionenergyisapproximatelyequaltothe
full sheath potential (Figure 10). 4g,50At 13.7 MHz, the ion energy distri-
butions are much narrower and the maximum ion energies are much less
than the peak sheath potential.
Because of the transition from resistive to capacitive sheaths above wi,
power dissipation must shift from the sheaths to the plasma. Since the
current is conducted primarily by electrons in the plasma, the dissipation
mechanisms must involve electron-neutral collisions: ionization, disso-
ciation, and excitation. The shift in power dissipation is evident in Figure
11, where the time-averaged concentrations of excited and ground state
BCI radicals are plotted as a function of position across the gap. Three
different frequencies are displayed; the shifts in emission intensity and
radical density are indicative of the shifts in where power is dissipated.
When these profiles are spatially integrated we can learn not only where
but also how power dissipation changes with frequency. As frequency
increases, the total, spatially-integrated densities of both excited and
Plasma Processing 205
0 4 0 12 16 20
AXIAL POSITION (MM)
104
r
n CP, 0.3 Torr
0.6 W/cm2
SS ELECTRODES
cm SPACING
I I I I
--- PIE
4
3 -
I
.
C
4 2.5 M HZ
4 0.25 MHZ
b.
3
._
0 4 8 IZ I6
o$p~~; , s’
3 0.01 0.1 I 10 100
EXCITATION
FREUUENCY
(MHz)
3- CC!,:Ar/l:l
0.4Torr
75 w
2-
SS ELECTRODES Al
I cmSPACING
l-
I
0.01 0.1 I IO 100
EXCITATION
FREUUENCY
(MHz)
Figure 12: (a) Sietch rate vs.frequency in a discharge through CC14. (b) Al etch
rate vs. frequency in a discharge through Ccl4 (from Reference 52, reprinted by
permission of the publisher, The Electrochemical Society Inc.).
VI A2
‘-= -9 (64
V2 4
Vl
-iE
v2
PERCENT 02
Figure 13: Si etch rate vs. O2 concentration in a CF4/02 discharge (after Refer-
ence 63).
bution can be seen in Figure 14 where the emission intensities from F and
Arare plotted as a function of the 0, feedstock concentration in a CF,/O,/Ar
discharge. As the oxygen concentration increases, the Ar emission intensity
decreasesindicatingthatthenumberofelectronswithenergiesabovethe
7504A line threshold must also decrease. In other words, the electron
energy distribution appears to cool. Note that the F atom emission lines at
7037Aand 6856A, which have similar excitation thresholds, increase with
0, concentration. The overall increase in F atom production by free radical
reactions (see above and Figure 14) more than compensates for the
decrease in excitation efficiency.
Two explanations for the cooling of the electron energy distribution
with the addition of 0, are plausible: (1) the enhanced production of F
atomsand the introduction of oxygen, both species being electronegative
compared to CF, compounds, results in electron cooling by electron
capture processes; or, (2) the lower ionization potential of 0 and 0, relative
to Fand CF,66 results in electron cooling by inelastic ionization of neutral 0
and 0,. Of these two, the second explanation is more viable because
electronegative gases do not attach high energy electrons very effectively
(see Reference 67 and references therein). On the other hand, the intro-
duction of a lower ionization potential gas necessarily reduces the average
electron energy since a lower energy collision channel has been opened.
212 Semiconductor Materials
4 I I
I I I
(a) ’ CF4 + O2
o Ar (7504%) INTENSITY
3 - EXCITATION EFFICIENCY
FOR F ( 7037% 1
.
2
0 F ATOM DENSITY N
0 I I I I I I
0 20 40 60
OXYGEN PERCENTAGE
Theimportanceoffeedstockionizationpotentialsalsoshowsupinthe
dynamics of ion production and loss. 22 For example, small additions of 0
to an Ar discharge cause dramatic changes in the ion composition. 68
Similarly, when Cl, is added to a low frequency discharge through N,, the
ion composition changes drastically. The concentration of Cl, in the
feedstock need be only 10% forthe N,+concentration to drop by an order
of magnitude. The concentration of Cl*+ is virtually the same in this mixed
plasma as it is in a pure Cl, discharge. 22This effect is illustrated in Figure
15, where the time-dependent concentrations of N2+ and Cl,+ are plotted
for various feedstock compositions. The extent of modulation in the ion
concentration waveforms is a measure of theground-state ion lifetime. For
the N,/CI, mixture, the N,+ lifetime is determined by the rate of charge
exchange with Cl222 (Figure 15d). Since the ionization potentials for N,
(15.58 ev) and Cl, (11.48 ev) are so different, the exchange is essentially
irreversible. The charge exchange reaction must proceed rapidly for the
steady-state ion concentration to be so greatly affected by small changes
in feedstock composition. In this particular case, the reaction rate is
probably enhanced by exchange through the excited state of Cl,+,
followed by,
time, which is simply proportional to P/e. This can be seen very clearly in
themassspectrometricdataofTruesdaleetal.6gwhichisreproduced here
as Figure 16.Theystudied the plasmadecomposition of C,F,asafunction
of both pressure and flow-rate and examined the stable products down-
stream from the dischargewith aquadrupole massspectrometer.Thefinal
concentrations of C,F,, CF,, and C,F,vary both with flow-rate at constant
pressure and with pressure at constant flow-rate. When the ratio of pressure
to flow-rate is held constant, however, the final product concentrations
also remain constant (Figure 16).
(cl
-. F = ZOCCIMIN
P/F = 0.04
c2F6
60-
A A ,
a-”
w - __ A A
d 8 m
z
t I
’ ‘CF 4
20
2 .4 6 8 I.0 2 .4 .6 .8 I.(
The best way to vary residence time, i.e. with minimal effect on other
plasma parameters, is to vary the flow-rate at constant pressure since
variation of pressure can have pronounced effects on the discharge
physicalproperties(see below). What happenstothegas-phasechemistry
when the residence time is varied? If we considerthe rate of decomposition
216 Semiconductor Materials
and radical production to be fixed for a given power density, then it is easy
to see that the less time the feedstock gas spends in the plasma volume,
the less extensive the degree of dissociation. Thus, at high flow-rates, or
short residence times, the radical density will decrease with increasing
flow-rate. At low flow-rates, or long residence times, the extent of dissociation
may reach a limiting value and a dynamic equilibrium will be established
such that the radical concentrations become flow-independent.
Let us consider the effects of flow-rate on heterogeneous chemistry
when the heterogeneous rates are not rate-limiting: i.e. when the surface
reaction probability is large. If the radical products of feedstock decom-
position are the primary surface reactants, then the low flow-rate radical
concentration may be negliblesince the heterogeneous reaction actsas a
radicalsinkandtheoverallreactionratemaybesmall becauseofthesmall
reactant flux to the surface. At high flow rates, the feedstock may be
insufficiently decomposed to provide radicals for the surface reaction and
again the heterogeneous rate may be small.Theoverall dependence of the
surface reaction rate on residence time, or flow-rate, will exhibit a maximum
(Figure 1 7).lgr70 The effect of flow-rate on heterogeneous chemistry when
the reaction is surface rather than reactant-supply rate-limited will be
discussed in Sec. 2.2.1.
2.7.4.2 Physical effects. The situation is very different when one con-
siders the effects of pressure and flow-rate on the physical properties of
the discharge and, in particular, the energy distributions of ions and
electrons. The effects of flow-rate on ion and electron energy distributions
are primarily an indirect consequence of the compositional changes
discussed above. However, pressure affects these distributions directly
FLOW RATE -
Figure 17: Generalized flow-rate dependence of etching rate for the case where
the surface reaction rate is fast and the reactant is generated in fheplasma. Two
limiting cases are obvious: (1) rate is reactant supply limited and (2) rate is re-
actant generation limited (after Reference 70).
Plasma Processing 217
by affecting collision rates. For example, the energy with which an ion
impactsasurfaceisnotonlydependentuponthesheathfield,rffrequency,
and ion mass but also the rate at which ions collide with neutrals as they
traverse the sheath.1g,56s71,72If the collision mean free path is greaterthan
the sheath thickness, then ions which enter the sheath from the plasma
boundary will be accelerated by the full sheath potential and will impact
the electrode with a narrow but highly energetic velocity distribution.
Alternatively, if the collision mean free path is much smaller than the
sheaththicknessand if theionsloseall theenergygainedfrom thefieldon
each collision, then
Ei = FX,
where Ei is the ion kinetic energy at the electrode and F is the sheath field,
assumed to be constant with position over one mean free path, X The truth
will generally lie somewhere between these two extremes.
Pressurealso affects the electron energy distribution. Manyelectron-
neutral collisions result in the formation of ion-electron pairs. If the electron-
impact ionization mean free path is small compared to the sheath thickness,
this can lead to substantial ionization in the sheath, which in turn will alter
the ion energy distribution at the electrode surface and the sheath
thickness.‘O In general, the overall charge density and plasma impedance
can be expected to change with pressure in a complex fashion. Since the
ion and electron collision mean free paths can be expected to scale
inversely with pressure, Equation 9 suggests that the natural variable with
which charged particle energy distributions and denities can be expected
to scale is neither pressure nor sheath field but rather the ratio F/P. This
has been long appreciated by scientists studying dc glow discharge
physics. As we will see in Sec. 2.1.7, it is also a useful parameter in
designing and understanding heterogeneous plasma-surface interactions.
2.1.5 Power Density. Many things can happen when the applied
power density is varied. Generally, ething rates, radical densities, charge
densities, and sheath fields increase initially and then saturate with in-
creasing power. One reason why saturation occurs may be that the plasma
volume often increases as the power is increased and electrons and ions
acquire greater energy. This expansion may result in the discharge“finding”
other grounds and discontunities can result in measured plasma para-
meters. If the plasma volume can be maintained at a constant volume, e.g. by
mechanical or magnetic73 confinement, then an increase in power corre-
sponds to an increase in power density. This in turn will result in higher
electron energies, sheath potentials, and ion energies.
a
A+Sol3 (10)
P
01)
Plasma Processing 219
i
I Bf
.i n
i
i
Bf
PV Ki-Kj
‘Bd
kT [
(Ki +l)(Kj+l)
1
Figure 18: Equivalent circuit for chemical vapor transport. The transport of
products from one surface to another is driven by an effective chemical potential
difference between the surfaces which may result in turn from differences in sur-
face temperature or ion bombardment. The effects of flow-rate are represented
by current generators so that in the absence of an effective chemical potential
difference, etching but not deposition may occur (after Reference 75). See text
for definitions of circuit elements.
where Ci is the area of the ith surface and ai=cri+Pi is the reaction conduc-
tance. Slow reactions correspond to large resistances. The diffusive resis-
tance between the two surfaces is given by,
pV (K’-Kj)
(13)
‘= kT(K’+l)(l<i+l)’
where KiJ = c$~‘//~~Jis the equilibrium constant at surface i,j. The term
chemical potential difference is used somewhat loosely here. Afinitevalue
for zeta implies that the equilibrium states for the two surfaces will not be
220 Semiconductor Materials
04)
where R = R,‘+R/+R, is the total circuit resistance. The reason for repre-
senting the surface reaction rates and diffusion rate as conductivites
should now be apparent. At zero flow the product current, or etching rate, is
limited by how fast etching occurs on one surface, deposition on the other,
and diffusion of reactants and products in between. If the heterogeneous
reaction rate happens to be limited by reactant generation from a homo-
geneous process, the product current will be limited instead by that
process. This can be seen in a formal fashion by supposing there is a
precursor to reactant A in Equation 10. In terms of the CVT equivalent
circuit (Figure 19), this means that there will be an additional resistance in
series with R,. In terms of the heterogeneous rate constant, ur-’ = uAn-’ +
UP ’
-l where up is the sum of the forward and backward rates for the
precursor reaction P +, A.75
It isimportanttonotethattheaboveequationfortheproductcurrentis
valid, within theframeworkof thesimplegenericequation,foranydeparture
from equilibrium. In equilibrium, Ki = Ki, there is no chemical potential
difference and no net transport of A or B from one surface to another. If a
chemical potential difference is maintained between surfaces i and j, for
example by application of a temperature or electrical potential difference
(see below and Reference 58), there will be net transport of A and B from
one surface to another. Whether B is deposited on surface i and etched
from surface j or vice versa depends on the sign of zeta, which in turn
depends not only on any temperature differences between the surfaces
but also on the reaction energetics. This can be seen most clearly by
expressing the equilibrium rate constants in terms of the surface tempera-
tures and the Gibbs reaction free energy, AG:
where~istheflowrateconstant,i.e.thereciprocalofthereactorresidence
time. Now, even in the absence of a chemical potential difference between
surfaces i and j, net transport of A and B can take place. Flow-rate can be
used as a process control variable in a fashion which is beyond the simple
variation of reactant concentrations by variation of residence time (Sec.
2.1.4).
i
SELECTIVITY = iJPDLy/~P,,,OXIDE
Iii = s+As
(17)
RT+ ART ’
,,v (ui,j~[Rd+Ri,j(~Rd+2)]
where As = k~
oL1
R’R j
and AR = ~R,(R’+R’+&R’R’+~~)
d
Both the driving potential difference and the reactive impedance are
modified by finite flow. At very large flow rates, the product current becomes
independent of both the chemical potential difference and the flow rate
and is limited only by the reactive resistance,
N,pC'I,
Z-B
M (194
where dz’/dt is the film thickness rate of change (i.e. the etching or
deposition rate), r is the gas-phase reactant lifetime, G is the gas-phase
reactantgeneration rate,pisthesubstratedensity,M isthesubstrategram
molecular weight, and N, is Avogadro’s number. The linear relationship ot
Equation 19 has been observed under a wide range of conditions. This
suggests that the assumptions inherent in Mogab’s theory may be too
restrictive.Thechemicalvaportransporttheorydescribedaboveallowsus
to examine loading effects when flow-rate effects cannot be ignored and
when reactant generation is not rate limiting.
Under conditions where the substrate transport rate is not limited by
reactant supply, the relationship between flow-rate, substrate surface
area and the transport rate is contained in Equation 17. At zero flow,
Note that Equation 19b also exhibits a linear relationship between the
reciprocal transport rate and substrate surface area. There is a direct
correspondence between the parameters in Equations 19a and 1 gb:
u.I *--). ,i
(20a)
transport current. In the case of zero flow, reactant A and product B are
transported in opposite directions from one surface to the other when the
chemical potential difference, <, is finite so that one surface acts as a
reactant generator for the other surface. Similarly, the reactant loss time
constant, r, corresponds to the sum of the reactive resistance at surface j
and the diffusive resistance (Equation 20~) both of which correspond to
reactant loss mechanisms.The important point is that despite the radically
different assumptions in the two theories, there is a correspondence of
sources and sinks such that the overall functionality remains the same.
Thus, observation of a linear loading effect is not sufficient to determine
the rate-limiting process.
At higher flow-rates, the relationship between reciprocal transport
rate and substrate area is no longer linear (see Equation 17). The loading
effect is predicted to go through a maximum as the flow-rate is varied.75
Thesensitivityofaprocesstothenumberofwaferswilldependontheflow-
rate in a non-linear fashion. Although some of the aspects of this theory
have been verified(e.g. that there isan effect of flow-rate on loading), more
experimental data are needed to assess the range of validity of Equation
17. One difficulty in obtaining such data is the change in discharge
composition which usually occurs when the residence time is changed.
(see Sec. 2.1.4). This could result in a change in the rate-limiting step from
heterogeneous to homogeneous. Another complication may arise when
more than one reactant is important, This effect alone can give rise toa non-
linear loading curve.81 If changes in flow-rate change the relative concen-
trations of these reactants further deviations from the simple theory above
can be expected. To test the range of validity for the CVT equivalent cir-
cuit model, studies of spontaneous reactions, i.e. without a plasma, would
be most appropriate.
2.2.2 Plasma Modified Chemical Vapor Transport. Until now we
have considered chemical vapor transport without considering the effects
of the plasma except to the extent that it modifies reactant concentrations
and diffusion coefficients. The major influence of the plasma, however, is to
modify the heterogeneous reaction rates by ion bombardment of the
surfaces. Although electron enhancement of heterogeneous rates has
been demonstrated82*83, the sheath fields ordinarily are such as to repel
electrons and negative ions from device surfaces (see Sec. 2.1.2). In this
section we will see how this effect can be treated in a formal fashion within
the framework of chemical vapor transport theory.28,58z74,75The following
section (2.2.3) will discuss the kinetic and microscopic origins of the
heterogeneous chemistry and plasma modifications.
In addition to the neutral reaction, Equation 10, we consider a parallel
ion-driven reaction:74t75
S + A+ f e- ++ B+ + e- ++ B
(21)
P+
This reaction is written in a formal fashion for simplicity. The ion need not
correspond to the neutral moiety nor must the product be an ion initially.
Electrons are included in Equation 21 merely to account for surface
Plasma Processing 225
cx = N*Zexp(-EA/kT) (22a)
,b = NBZexp(-EB/kT) (22b)
whereZisapreexponentialratefactor,fheisthedegreeofionizationofA,B
and “I@ = EAB - E*+ f3+is the difference between the neutral and ionic
activation energies. ’
Expressing the rate constants in this form allows us to explain a large
number of temperature studies where Arrhenius behavior has been ob-
served (Figure 20). For the case where the back reactions are negligible
and reactant supply is constant,74
For neutral dominated reactions, f, exp(U,/kT)< < 1 and the slope of 1,~s.
l/T gives the neutral activation energy, E,. For ion-dominated reactions,
f,exp(U,/kT)>>l and the slope gives -E,,. The fact that ion-dominated
reactions give smaller slopes (see Figure 20) indicates that the ionic
activation energy is smaller than the neutral activation energy. This is not
surprising considering that the ions have been accelerated to relatively
high energies by the sheath field. Another way to think of this activation
energy difference is to think in terms of effective temperatures. Because
the ion energy is superthermal, the effective temperature for the ion-
surface interaction is much larger than forthe neutral-surface interaction.
226 Semiconductor Materials
T('C)
300 280 260 240 220
I I I I I I I I
lOO(
g 100
1
Y
2
QL
3
t
10
P(W/cmV E(kcal/mole)
v 0.89 28.4 0
0 0.40 30.6 0
a 0.15 38.8
. 0.62 34.5 ,\ 0
1 0 -
0
I , ,
1.8 1.9 2.0
1ooorT
Figure 20: Arrhenius plots showing effects of temperature and power density
on etch rates of InP and relative In atom emission intensity in an 0.3 Torr dis-
charge through Clz. The emission intensities are proportional to etch rate. Note
the smaller slopes for the higher power densities suggesting that ion-bombard-
ment is reducing the overall activation energy (from Reference 181, courtesy of
V.M. Donnelly).
Plasma Processing 227
1+ fAexp(U&T)
IC, = I<
I + f,yxp(U,/kT) I’
- [AG-AU- kTln(fA/fB)
I<, -+ exp w
kT
where K, is the total (i.e. ionic and neutral) equilibrium rate constant, AU =
U, - U,. Thus, the effect of ion bombardment is a modification of the
effective free energy. Alternatively, we can think of the surface temperature
as being modified, T’ = T/(1-(kTInfJf,+AU)/AG]. If AlJ/kT>ln(f,/f& the
effectivetemperaturewill be hotterthantheactualsurface temperature. In
termsof the effective temperature, transport still goesfrom“hot” to”cold”
for endoergic neutral reactions and from “cold” to “hot” for exoergic
reactions. The difference is that “hot” and “cold” depend upon not only
surface temperature but also ion energies,which can be controlled by
frequency, dc bias, and pressure. Thus, we see the complementary nature
of ion-enhanced chemistry and thermally enhanced chemistry. Of course,
they are not the same thing. The product distributions are likely to be very
different since the thermal energy deposition will be statistically distributed
tothevariousdegreesoffreedom buttheionictranslationalenergymaybe
disposed in very specific ways.
2.2.2.7 Anisotropy. Now that weformallyunderstand theeffectsof ion
bombardment on heterogeneous reaction rates we can understand how
anisotropic patterning (see Sec. 3.4.2) is possible. In general we need to
consider transport between not only the two electrode surfaces but
among all surfaces. Specifically, when considering anisotropy, these
surfaces reside on the same electrode but are mutually orthogonal. The
surface which is perpendicular to the electric field lines will experience
more energetic ion bombardment than the surface which is parallel to the
electric field lines as long as the ion transport across the sheath is
anisotropic. This will occur at higher values of F/P (see Sec. 2.1.4 and
Figure 21).
The specific value of F/P which gives a particular anisotropy in the ion
transport directionality will depend upon the gas composition and operating
frequency(see Sec. 2.1). For example, the value of F/P, estimated from the
square root of the rf power density,** required for a given degree of
anisotropy is shown in Figure 21. Fora pure Cl, discharge, where resonant
charge exchange can be very effective in reducing the ion energy aniso-
Next Page
.A;, /p ,‘(’
. ,‘
I,.
,
.) ,, , 1,’
Figure 21: Anisotropy as a function of the ratio of the square root of rf power
density to pressure (approximately equal to F/P) for three discharges illustrating
the effect of gas composition on the degreeof anisotropy in the ion energy distribu-
tion (from Reference 28, reprinted by permission of the publisher, The Electro-
chemical Society Inc.) and the value of F/P required for a given anisotropy.
6
Physical Vapor Deposition
John A. Thornton
Department of Materials Science and Coordinated Science Laboratory
University of Illinois
Urbana, Illinois
1. INTRODUCTION
329
330 Semiconductor Materials
T-t
Ic, VANE PUMP -
VACUUM VACUUM 0
CHAMBER CHAMBER
1
SORPTION
ROUGHING
-29 GATE VALVE GATE VALVE GATE VALVE
tf3b
TRAP
/ I”I “I
DIFFUSION OR
TURBOMOLECULAR
PUrviP
TITANIUM
\ LN2 CRYOPUMP
CRYOPUMP MECHANICAL
PUMP
3!
SUBLIMATION 5-.
SPUTTER
PUMP c)
ION PUMP ”
MECHANICAL
PUMP c <
Figure 1: Schematic drawings showing vacuum pumping systems of the various types used in deposition technology: (a) dif- w
0
fusion or turbomolecular pump configuration, (b) getter pump configuration, and (c) cryogenic pump configuration.
332 Semiconductor Materials
discussed in References l-4. Our concern here is with the “vacuum state”
that is achieved and its implications on the deposition process.
Consider a cubic vacuum system with sides 1 m in length. The volume
is 1 m3. The internal surface area is 6 m2. An ideal metal surface contains
about 2~10’~ adsorption sites per cm2. When the chamber is exposed to
the atmosphere, an even largerdensityof molecules will become attached
onto the walls because of surface irregularities and multilayer adsorption.
The number of molecules per cubic centimeter in a room temperature gas
is about
where p is the pressure in Torr. Thus when pumping is initiated we have the
taskof removing(760 Torr)(3.3xl 016 molecules/cm3-Torr)(l 06cm3) = 2.5x
1025 moleculesfrom thevolumeandat least(2xl 0i5 molecules/cm2)(6xl O4
cm2) = 1.2x1 020 molecules from the walls.
Nowconsiderachamberfiiledtoan initial pressure P,with an idealgas,
which has no interactions with the walls other than reflections. When such
a chamber is evacuated by a pump of constant volumetric efficiency, the
pressure will decrease with time according to the equation.
wherethetimeconstant,r=V/S,isafunctionofthechambervolume(V)and
pumping speed(S). Suppose that the pumping system for our 1 m3 cham ber
has a speed of 1000 liters/set. Then z = 1 sec. Typical values are in the
range from 0.1 set to a few seconds.5 Thus Equation (2) predicts that the
pressure in our chamber will decrease exponentially once pumping has
commenced, dropping almost an order of magnitude during each 2 set
interval.
Now consider the removal of atmospheric gases from a chamber. The
relationship given by Equation (2) is obeyed, after a few seconds from the
start of evacuation and until a pressure of about 10 Torr is reached, as gas
is removed from thevolume of thechamber. Subsequently, theevacuation
will become rate-limited by outgassing from the chamber walls. Under
these conditions, the pressure will decrease much more slowly, obeying
an equation of the form
where Q(t) is the total outgassing rate from the surfaceswithin the chamber
atthetimet.Theoutgassingrate,Q(t),andthereforethechamberpressure,
decrease as a function of time, because internal diffusion and surface
desorption deplete the reservoirs of stored gas entrapped on thechamber
internal surfaces.
The dwell time of an atom or molecule on a surface under vacuum will
depend on the binding energy between the molecule and the surface, and
on the surface temperature. See Equation 13 in Section 6.1. Physisorbed
gases with binding energies of the order of 0.1 to 0.5 eV desorb quickly
Physical Vapor Deposition 333
whereAisthechambersurfacearea.Thetimet,isareferencepointwhere
q = q,. For approximate calculations t, can be taken as the point at which
the high vacuum valve is opened.
Consider the case of our 1 m3 chamber after three hours of pumping.
Assume that the chamber is constructed of stainless steel and that the
high vacuum pump has a speed of 600 liters/set. The specific outgassing
rate for stainless steel after 180 min. of pumping is seen in Figure 2 to be
about 4x1 O-*Torr-liters/set-cm”. From Equation 3 we estimate the chamber
pressure to be (4 X 10W8Torr-liters/set-cm2)(6 X 1O4 cm2)(600/liters/sec)
or 4 X 1 0e6Torr.
In practical deposition systems unwanted gases are a result of desorp-
tion from the deposition sources, substrates, and hot filaments as well as
the chamberwalls. Back-streaming gases from the pumps alsocontribute
contamination species. After prolonged pumping, the residual gases are
typically H,O, CO, CO,, 0, and N,.’ For a residual gas pressure p (given in
Torr), the impingement rate on a substrate surface is
\
\
\
\\
\
\
1o-9 7 \
\
Figure 2: Specific outgassing rate as function of time for various materials. From
Reference 6.
SUBSTRATES SUBSTRATES
/
COATING
FLUX
Pe
RESIDUAL
r-1
I _ \ GAS PRESSURE
SOURCE
1
VACUUM PUMPS PUMPS
CHAMBER ?
VACUUM Y
cn.
CHAMBER
c
5
B
TOTAL OUTGASSING FLUX Kp/ m
IMPURITY LEVEL - IMPURITY LEVEL - ;;
TOTAL COATING FLUX TOTAL COATING FLUX
8
8.
Z
a b 2
Figure 3: Schematic illustration showing influence of apparatus geometry on the way in which wall outgassing affects the coating
impurity level.
336 Semiconductor Materials
3. EVAPORATION
3.1 Introduction
In the evaporation process, vapors of the coating material are released
from a source because of heating. The source material may be in the liquid
or solid state, depending on its vapor pressure relative to its melting point.
Almost any conceivable method can be used to heat thesource. One of the
most common methods is resistive heating, either of the source material
itself or of a support containing the material. Other common heating
methods involve the use of an electron beam, a laser beam, or an arc
discharge to produce local surface heating of the source material.
The evaporation process is usually carried out at a sufficiently low
pressure(typically 1O-5to 1 0e6Torr) so that the evaporated atoms undergo
an essentially collisionless “line-of-sight” transport to the substrates. In
this connection it is useful to remember that the mean free path of gas
particles is about equal to
X= 5/p,, cm (6)
0 MELTING POINT
where r is the distance from the source to the substrate. The Co@ term
accounts for the fact that the substrate may not be perpendicular to the
lineofcentersconnectingthesourceandsubstrate.Thethicknessgrowth
rate of the film is given by
DC MW, (9)
PNA
where N, is Avogadro’s number, and M is the molecular weight and p is the
mass density of the deposit.
Depositionratesfortheevaporationprocessareclearlydependenton
the material being evaporated, the type of evaporation source, and the
position and orientation of the substrate surface. Consider the case of
aluminum evaporation from a 1 cm diameter source onto a substrate
located 15 cm directly above the source (COSQ = 1) and oriented to face
the source (Cosf3 = 1). Thus we have r = 15 cm (this is a typical distance),
and Equation 8 yields W, = W/707. From Table 1 we see that a source
temperature of 1220°C is required to provide p* = 1O-2Torr. We assume
that this temperature is used. Therefore, from Equation 7 with p* = lo-*
Physical Vapor Deposition 339
02 pressure dt
2000% 9
MgP2 >@W2)2
MgF2 1263 1130 0.91 MO
0W2)3
CaF2 CaF2,CnF 1418 1300 0.93 M0
Wd=~CosQCorcl
II J
Figure 5: Flux passing from small area source to elemental substrate area dAS
which inscribes solid angle dw.
prices to be discarded after one use if necessary. The wire or foil supports
must be fabricated from materials which have negligible vaporordissocia-
tion pressures at the operating temperatures. These temperatures are
typically in the range from 1000 to 2000°C. Wetting of the wire or foil
surface bytheevaporant isalsodesirableinordertoachievegood thermal
contact. Detailed recommendations pertaining to wire or foil support
materials for various evaporants are given in References 8 and 9. These
recommendations for a few evaporants are summarized in Tables 1 and 2.
The most commonly used support materials are tungsten, molybdenum,
and tantalum. Suitable wire or foil sources are available to evaporate small
charges of nearly all the elements except the refractory metals themselves.
Themaximumcapacityofwireandfoilsourcesistypicallyafewgrams.
The usual approach is to calculate the charge of source material that will
provide a given deposit thickness, using a relationship of the form of
Equation 8 for the apparatus geometry in question, and then to evaporate
the entire charge. The wire/foil approach is in general too time-consuming
for most production applications, but it is an effective method, for example,
for depositing test electrodes in laboratory studies.
3.3.2 Crucible Sources. Crucible sources are required to support
molten metals in quantities of a few grams or more. Since the melt is in
contact with the container for prolonged periods of time, the selection of a
noncontaminating and thermallystablecruciblematerial isveryimportant.
Detailed recommendations for crucible materials are given in References
8 and 9. The non-metallic support materials summerized in Tables 1 and 2
apply strictly to crucibles. Thus it is seen that graphite and the refractory
METAL FOIL SOURCE RESISTANCE HEATED CRUCIBLE
RESISTANCE HEATED WIRE
==I!@&\-
a
b
SUBLIMATION SOURCE
RF HEATED CRUCIBLE
RADIATION
SHIELDS
MOLTEN METAL
d e
Figure 6: Schematic illustrations showing several types of evaporation sources.
Physical Vapor Deposition 343
with flow orifices that are restricted but not necessarily small enough to
satisfy the free molecular flow conditions that are implicit in the particular
case of the Knudsen cell. The flow from effusion cells can be theoretically
predicted if the fluid dynamics of the flow through the orifice is properly
taken into account. However, as a general rule the emission characteristics
as a function of temperature from effusion cells, including Knudsen cells,
are determined experimentally. Figure 7 shows an array of effusion cell
sources used for multi-source deposition of semiconducting coatings of
CulnSe,. ” Effusion cells play a very important role in the process of
molecular beam epitaxy, which is discussed in Section 4.
3.3.6 Electron Beam Sources. Figure 8 showsaschematicdiagram
of an electron beam evaporation system. Since the beam is concentrated
on the evaporating surface, while other portions of the evaporant are
maintained at lower temperatures, the evaporant can form its own crucible.
Hence, interactions between evaporant and support materials are greatly
reduced. Electron energies are typically in the 3 to 10 KW range, with
power levels in the range2 to50 KW.Therefore, relatively highevaporation
rates can be achieved, even for the refractory metals. Rod-fed sources can
provide a large inventory of coating material. Therefore, electron beam
sources are the most commonly used evaporation sources for large scale
production applications. However, two difficulties that must be dealt with
are, first, that electron beam sources are vulnerable to spitting because of
, WART2 CRYSTAL
MONITOR
SUBSTRATE
VACUUM -
CHAMBER
/
SHIELD
SOURCE ’
HEATER
ASSEMBLY
/ PUMPING
THERMOCOUPLES
L-J PORT
SUBSTRATES
VAPOR FLUX
MAGNETICALLY FOCUSED
MOLTEN ELECTRON BEAM
POOL
EVAPORANT
SUPPORT
BEAM SOURCE
SOURCE MATERIAL
TO VACUUM
(MAY BE ROD WITH
PUMPS
AUTOMATIC FEED)
the high power densities at the point of beam impact, and second, that the
deposition flux is nonuniform as discussed in Section 3.4.
3.3.7 Other Types of Evaporation Sources. Several othertypes of
evaporation sources have been developed to deal with the particular
problems associated with forming stoichiometric coatings of alloys and
compounds. These are discussed in Section 3.5.
SUBSTRATES
SOURCE /’ -
by direct evaporation. This has led to the use of special methods such as
flash, two-source, and reactive evaporation. These processes are discussed
in Section 3.6.
3.5.2 Evaporation of Alloys. The constituents in alloys evaporate
independentlyof oneanother, mostlyassingleatoms. However, thevapor
pressures of the individual constituents are not equal to their pure metal
values at the temperature in question, because there is a contribution to
the chemical potential when one metal is dissolved in another. Most
metals evaporate incongruently, and this has led to the use of sputtering
where extreme composition control is necessary. An example is the
deposition of Nichrome (80%Ni - 2O%Cr) to form thin film resistors. How-
ever, in the important case of Permalloy(85%Ni - 15%Fe), the evaporation
issufficientlycongruent topermittheuseofsimplesinglesourceevapora-
tion for many applications8
Electron beam evaporation can significantly expand the range of
materials which can be evaporated with reasonable composition control.1°
This is possible because the electron beam source creates a small molten
region, as shown in Figure 10. During an incubation period the molten
region becomes deficient in the volatile species. The composition is then
rate-limited by the passage of material by diffusion from the solid into the
melt, across interface “A” in the figure. If the vapor pressure difference is
not too large for the constituent diffusion rates, a steady state is developed,
where the composition of the melt is just such as to produce a vapor
composition equal to that of the solid. It is reported that reproducible
compositions of Ni-20Cr, Ti-GAI, Ag-SCu, Ag-1 OCu, Ag-20Cu, Ag-30Cu,
and Ni-xCr-yAl-zY have been successfully achieved by electron beam
evaporation.1°
3.5.3 Evaporation of Compounds. In the evaporation of compounds
the transition to the vapor phase rarely occurs without changes to the
molecularspecies.Thusevaporation isusuallyaccompanied bymolecular
dissociation, association, or a combination of both processes. Dissociation
represents thermal decomposition and generally makes simple direct
evaporation impractical. The species formed in the direct evaporation of a
number of compounds are summarized in Table 2.
EVAPORATED
FLUX
RATE LIMITING
FLtiX
SOURCE -
MATERIAL
and compounds. In most cases, the vapors impinging on the substrate are
highly supersaturated, so that the film composition is not affected by the
condensation coefficients. (Condensation coefficients are discussed in
Section 6.1.) The most common problem is incomplete evaporation due to
particle ejection and deflection. Ni-Cr alloys, Cr/SiO cermets, GaAs, InP,
CL@, and BaTiO, are examples of materials that have been deposited by
flash evaporation.
3.6.2 Hot-WallEvaporation. Inthistechniquefilmsaregrownunder
conditions that areclose to thermodynamicequilibrium.15-16aAschematic
drawing of a hot wall evaporation apparatus is shown in Figure 11. The
evaporated flux is passed into an enclosure with walls held at a sufficiently
high temperature so that condensation is precluded. Accordingly, stoichio-
metric coatings can be deposited, even on substrates maintained at such
high temperatures that one or more of theconstituents has a lowcondensa-
tion coefficient. Since wall condensation is prohibited, the vapor pressures
of the volatile constituents simply build up until they deposit onto the
substrates at steady state rates that are equal to the rates at which they
enter the enclosure from the evaporation sources. For example, near-
stoichiometric CdTe films have been evaporated from a single CdTe
source maintained at 600°C with a hot-wall temperature of 5Oo”C, and a
SUBSTRATE
TEMPERATURE -T,
HOT WAL
T > T,
VACUUM
CHAMBER
\ EVAPORATION
SOURCES
PUMP
GAS
INLET
VACUUM -
CHAMBER I
HEATER
- SUBSTRATES
WALLT>T, -
NO CONDENSATION
- SUBLIMATION
PUMPING
PORT
“LI.““II
SELENIUM
(In & Cul I I CRUCIBLES
- HEATER
CONTROL
CONTROL
INDI IUM
BOA T
WATER COOLED
SHIELD PLATE
COPPER
BOAT
L
Figure 13: Multi-source evaporation system using open crucibles to deposit
CulnSe* semiconducting coatings. See Reference 19.
352 Semiconductor Materials
2
10-l 100 10’ 102 103
CLECTRON BEAM
4.1 Introduction
Molecularbeamepitxy(MBE)isamulti-sourceevaporationprocess,of
the type discussed in Section 3.6, which is done with extreme control over
the deposition parameters in order to exploit the kinetic processes of film
growth that are discussed in Section 6.1. MBE has been applied primarily
tothegrowthofsinglecrystalfilmsofcompoundsemiconductors.Thermal
molecular beams of each constituent of the film are directed to converge
onasinglecrystalsubstrateunderconditionssuitableforepitaxialgrowth.
Deposition rates are low (typically about 0.1 rim/see). The low deposition
rates reduce the temperature required to achieve epitaxialgrowth(Figure
47 in Section 6.2). The low growth rates are made feasible by the use of
ultra-high vacuum systems which have base pressures in the 1 O‘lOto lo-l1
Torr range and thereby reduce the residual gas contamination flux incident
on the substrates, as discussed in Section 2.
The slow growth rates permit very precise control of layerthicknesses
in the nm range. Shields are used to provide abrupt initiation or cessation
of the molecular beam fluxes and thereby to create sharp interfaces or
356 Semiconductor Materials
AUGER
ANALYZER
QUADRUPOLE MASS
SPECTROMETER ,, -) ION SPUTTER GUN
‘--*I-L
*_ FLUORESCENT
_.-_.-
3 NITROGEN
Ll.>nn, I_
EFFUSION CELLS
Figure 16: Schematic illustration of MBE deposition chamber with sources con-
figuredfor depositing GaAs type films.
Physical Vapor Deposition 357
istypicallycoveredwithacontaminationlayerconsistingofcarbonaceous
compounds and oxides which must be removed. The most detailed MBE
workhasinvolvedthedepositionofGaAsontosinglecrystaIGaAssubstrates.
The (100) plane is chosen for the technologically important reason that it
has orthogonal (1 10) type cleavage planes.27 A typical preparation pro-
cedure involves polishing with a diamond paste followed by etch-polishing
on an abrasive-free lens paper soaked with a sodium hypochlorite or
bromine/methanol solution.** The substrate is then rinsed in trichloro-
ethylene, methanol, and distilled water, following which it is boiled in
hydrochloric acid, free etched in sulfuric acid, again rinsed in distilled
water, and finally soldered with In to a MO backing plate. At this point the
substrates, which are generally contaminated with both oxygen and a
small amount of carbon, are introduced into the MBE chamber.
A typical in situ process within the MBE chamber involves the use of
ion bombardment toremoveforeign materialsfrom thesurface,andAuger
spectroscopy to verify that the surface is clean. The sample is then passed
into the deposition chamber and annealed to remove the surface damage
created by the sputter cleaning. In many cases thermal desorption rather
than ion bombardment may be used to avoid composition changes due to
preferential sputtering. It is well established that temperatures in the 525
535°C range evaporate the passivating oxide film on GaAs without causing
surface composition changes because of incongruent evaporation of the
GaAs.28 In any event, the HEED system is used to verify the crystalline
quality of the substrate surface prior to deposition.
The sources are then adjusted to tempeatures that provide the desired
deposition rates. The shutters are closed during this operation to prevent
depositionon thesubstrates.Thesubstratetemperatureisadjustedtothe
desired value and the appropriate shutters are opened to commence
deposition.The HEED system can be used to provide periodic verifications
that epitaxy is proceeding. The quadrupole mass spectrometer is used to
monitor the beam fluxes and the residual gas partial pressures. Finally the
coated substrates are withdrawn through the load-lock system.
ions appears to permit them to penetrate into the lattice of the growing
coating to a sufficient degree so that the incorporation probability is
increased.
The complexity of the doping problem is illustrated by the GaAs
technology. Commonly used dopants are Sn, Si, and Ge for n-type and Be
for p-type GaAs.27~*8Tin, which is the most commonly used n-type dopant,
illustrates the complexities that can occur. Surface segregation causes
the Sn to accumulate on the surface in a concentration which is several
ordersof magnitude largerthan that in the bulk.The rateof incorporation is
controlled bythe Sn surfaceconcentration and theGavacancyconcentra-
tion within the GaAs.27 The Sn segregated on the surface precludes the
formation of abrupt changes in doping concentration, since itsconcentra-
tion cannot be reduced to zero by simply closing the shutter at the Sn
source. Anothersource of complexity occurs because many of the dopants
are amphoteric. For example, under As-rich conditions Ge tends to be
incorporated as a donor, while under Ga-stabilized conditions (see Section
6.3) Ge is incorporated predominantly as an acceptor.27 These examples
illustrate how the dopant incorporation is dependent on relative substrate
arrival rates of As and Ge atoms as well as the doping flux itself.
The difficulties in forming p-type GaAs are even more severe. The
conventional acceptor dopants for GaAs such as Zn and Cd have high
vapor pressures and therefore low incorporation coefficients. Beryllium
has shown the most promising p-type doping properties, but isan extreme
toxicity hazard.27 The ion beam technique has been successfully used to
incorporate Zn+and provide carrier concentrations in the 101gcm-3 range.28”7
Manganese has been used but causesadverse surface degradation.27 Ion
beam deposition has also been proven effective in Si MBE.40,41
4.5 Applications
MBE is particularly effective when control overthickness, composition,
and doping profiles are critical to device performance. Such requirements
are often encountered in microwave and optoelectronics devices. These
needs have stimulated the development of MBE technology in general
and GaAs technology in particular. MBE has been used not only to produce
state-of-the-art performance in conventional structures, but also to produce
totally new types of thin film devices. Table 5 lists some devices which
contain epitaxial structures grown by MBE.
GaAsfield effect transistorsare typically used as low-noise microwave
signal detectors and microwave signal generators. Both low noise and
high-powerfield-effect transistors require n-type layers less than 1000 nm
thick Low-noise FET’s have been reported using 100 nm thick, heavily
doped, MBE GaAs layers.27 The linearity of power FETs can be improved
bytailoringthedopingprofile,arequirementthatcan beachievedbyMBE.
Microwavevaractors,mixerdiodes,andIMPATTdiodesareotherexamples
of devices in which controlled doping profiles are required and therefore
where MBE is useful.27
The formation of low-resistance contacts to n- and p-type GaAs, as well
as Schottky barrier diodes on GaAs, is also of great technological importance
Physical Vapor Deposition 363
Laser diodes
Waveguides
Integrated optics
Tnper couplers
Light-emitting diodes
Photodetectors
Other devices
Diodes
MIS capacitors
Superlattices
Tunnel triodes
Solar cells
cavity. The band energies of the active layers are such that the charge
carriers injected under forward bias are trapped in the active region. In
recent work, lasers with unique performance have been fabricated by
making the active layers have the form of quantum well superlattices
consisting of alternate layers of materials with different bandgapsand with
layerthicknesseslessthantheDebyelength.Thus, inoneexample,fourteen
GaAs quantum-well active layers, only -14 nm thick, were sandwiched
between AI,,,Ga,,,,As confinement layers ~13 nm thick.34j36 Injection
MQW laser diodes are of great importance in fiber optics communication
systems, because laser operation can be tailored to emit at frequencies
well above the standard lasing frequencies of the host material, and
therebytocontrol the losses in thefiberoptics. The deposition ofsuperlattice
structures with properties not found in homogeneous materials is an
active area of current research which can be expected to yield a host of
applications in the future.42t43
5. SPUTTERING
5.1 Introduction
Sputtering is a process whereby material is dislodged and ejected
from the surface of a solid or a liquid due to the momentum exchange
associated with surface bombardment by energetic particles. Asource of
coating material called the target is placed into a vacuum chamber along
with thesubstrates, and thechamberisevacuated toapressure typically in
therange5x10-4t05x10-7Torr.The bombardingspeciesaregenerallyions
of a heavy inert gas. Argon is most commonly used. The sputtered material
is ejected primarily in atomic form. The substrates are positioned in front of
the target so that they intercept the flux of sputtered atoms.
The most common method of providing the ion bombardment is to
backfill the evacuated chamber with a working gas in the 1 to 100 mTorr
pressure range and to ignite an electric discharge with the target serving
as the cathode or negative electrode. Such an apparatus configuration is
shown schematically in Figure 19. Applied potentials are typically between
500 and 5OOOV. Direct currents are generally used when the target
material is a good electrical conductor. Radio frequencies are used when
the target material is poorly conducting or an insulator. Deposits of poorly
conducting metallic compounds can also be formed by dc sputtering the
metallic component while injecting other constituents in the gas phase.
This is known as reactive sputtering. A voltage bias may be applied to the
substrates so that they are at a negative potential relative to the plasma
and therefore subject to an ion bombardment that can influence coating
properties. This is known as bias sputtering.
The most striking characteristic of the sputtering process is its uni-
versality. Since the coating material is passed into the vapor phase by a
mechanical (momentum exchange) rather than a chemical or thermal
process, virtually any material is a candidate coating. Films containing
almost every element in the periodic table have been prepared bysputter-
Next Page
CATHODE
WORKING
GAS FEED
ION FLUX
TO VACUUM
PUMPS
Richard B. Fair
Microelectronics Center of North Carolina
Research Triangle Park, NC
INTRODUCTION
455
456 Semiconductor Materials
CONTINUUM THEORY
J, = -D dC (1)
dx
= -*r a-l
JI-J, (2)
ax
s
-=- a D-ac
at ax I axI (4)
For the special case where D is constant and the surface concentration of
the impurity that is diffusing is fixed, then Equation 5 results:
a’c
D -=- i)c
(5)
ax2 6r
Fick’sSecondLawisacontinuityequationwhichdescribesthetimerateof
change of the impurity concentration. The diffusion constant D is in units of
cm2/sec. and the concentration C is usually in units of atoms/cm3.
Special Cases
Predeposition. UnderthespecialcasewhereDisconstant,thesurface
concentration of the diffusing impurity is fixed, the concentration of the
impurity at x=m is C( - ,t)=O for all time, and the concentration at any point
inthecrystalatt=OisC(x,O)=O, thenundertheseconditionsthesolutionto
Equation 5 is given as:*
Diffusion and Ion Implantation 457
Jl
J2
i i
I
I
I I
-I I- (cl
JI I I J2
Ax -t I
Ic
c (xJ.t)= coerfc -
I =I2JEJ’ (6)
I I
I
c Cx,t)= coelfc - * CR
z&it (7)
10-l
(9)
cB
I:] = 24% erfc-’
II
-
CCI (10)
a) c W = c, erfc (x/a& 1
c (I,&
b) -y-- = lo-’ ; xp& = 2.32 (from Fig. 2)
0
xj3
c) t=- : D = 22x10-” cm2/sec; t = 528.4 UC.
215D
w3+ 12 Si z 2B + (11)
The production of B,O, may come from either one of the reactions described
in Equations 12 or 13,
The source of the boron nitride in Equation 13 may be in the form of disks
460 Semiconductor Materials
about the size of a silicon wafer which are placed next to the wafers in the
diffusion furnace.
By varying the partial pressure of the gas phase of the dopant in the
furnace, it is possible to change the concentration of impurities in the
silicon. Henry’s Law relates the concentration of dopants that are introduced
in the furnace to the surface concentration:
Figure 4 shows Henry’s Law plotted. It can be seen that once the solid
solubility of boron in silicon is reached, Henry’s Law no longer applies.
Thus, most predeposition steps operate with a high enough partial pressure
in thedopantgasphasethatsolidsolubilityofthedopantisachievedinthe
silicon. This provides a natural control for reproducible diffusion results.
Diffusion and Ion Implantation 461
Solid rolubility of
2X1020 ----c-------_--_--
B in Si at 1100%
0’
/
/
I
/
1 /
/
0’
1 -/- Line correspondsto Henry’s law,
with H 12 X ld5 #m/cm3
.V 0 0.5 1
Pt (Torr)
1.5 2XlQ2
Ws + SiO, + 2P . (15)
Sources of P,O, vapor are solid P,O,, red phosphorus, POCI,, PBr,,
NH,H,PO, or PN.
For the predeposition step the goal is to deposit some number of
atoms/cm2 in the silicon substrate, Q(t). The way in which that number is
calculated is to integrate the totalconcentration percubiccentimeterfrom
0 to - as shown in Equation 16.
= co
Ipi.I 7r
(16)
c (x,t=o-l = Qs,,,
c (ms) = 0
then the distribution of impurities after diffusion for a time t is given by a
Gaussian function solution to Equation 5:
2
Q
c (x,t) = -jE$ “P
II
$ * (17)
lXIO1’ 410-4>2
c (1.t) = lxlols = - -
GE erp I 4Dt I
The solution to the equation above must be graphed and is shown where
the left-hand-side of the equation intersects the righthand side of the
equation in Figure 5. Note that if x = 0, the surface concentration C,
decreases as 1/A as shown in Equation 18. The junction depth can be
found by taking the log of both sides of Equation 17:
Q
co
(t)=J;rij; (18)
(19)
(20)
Diffusion and Ion Implantation 463
10-3
3
s
N.-
x
2
u
ii
I,oQ
10-4
l( 103 104
t (sac)
Figure 5: Graphed solution to the equation cited in the example in the text.
where
Diffusion Coefficients
Diffusion coefficients are based upon the atomic behavior of the atom
in a host lattice. Diffusion coefficients obey an analytical form as described
in Equation 21.
1x 105 n
0.6 .-
0.6 -.
0.4 _ Zdiii=O.lr
Figure 6: The Gaussian function: normalized concentration vs. distance for suc-
cessive times.
Impurity P As Sb B Al Ga In
0.6
~Dt)p,&(Dt)d,ive = 1.0
0.4
Exact
Uncoupled
10-3
3
Diffusion Mechanisms
The atomistic theory of diffusion is concerned with describing how an
atom gets from one part of a crystal to another. The lattice sites in a crystal
aregenerallytakenasthefixedlocationoftheatomsmakingupthecrystal.
It is known that the atoms oscillate around these lattice sites which are
their equilibrium positions. These oscillations lead to finite chances that
an atom will move from its lattice site to another position in the crystal.
There are several ways by which atoms can move from one site in the
crystal to another. These mechanisms are
00000 0000
00000 00&O
oo+ 00 0000
00000 (b) Thr interstitial diffusion mechanism.
(a) The vacancydiffusion mechanism.
0000
o@do 0
0000
(c) The interstitialcy mechanism.
change position along the x axis indicated in Figure 9. The atoms in this
simple case are taken to be located in planes at x, and x, + a, as shown in
the figure. The flux J is simply the concentration C times the velocity v:
Jx =cv (22)
The net flux is the difference between the flux to the right and from the left:
(24)
(29
x, %+a,
Figure 9: Flux in the x direction through the unit area in unit time. The planes
of the unit area are located at x = x0 and x = x0 + a,.
468 Semiconductor Materials
For motion by discrete jumps between planes a,, apart, the velocity is the
number of jumps per second, lY, times the distance a,, of each jump.
Equation 25 may now be written as
1 dC
J, = --&-,
2 dx (26)
1 *2
D=-aa,r
2 (27)
r = X,W
These terms are related to the Gibbs free energy change for vacancy
formation through the equation
where AS,,, and AH,,, are called the entropy and enthalpy of motion
respectively.
Now, the frequency, w, at which an atom and an adjacent neighboring
vacancy exchange can be written
Diffusion and Ion Implantation 469
Activated atom
0 0 0
88 0 0 0 88
(a) (c)
Figure 10: The sequence of (a), (b), and (cl show the movement of the atom
from a normal lattice site to an adjacent vacancy. Part (d) shows the variation of
free-energy as the atom moves from (a) to (c).
w=x,y (32)
where the frequency, y, is generally net known and is usually taken as the
lattice vibrational frequency of an atom about its equilibrium site, which is
of the order of 1013/sec-1. Now, from Equations 28,29 and 31 and 32 the
jump frequency for vacancy self-diffusion is
and
D = X, X, y (38)
Also the activation energy for vacancy diffusion depends upon the energy
necessary to form the vacancy and to move the lattice atom into an
adjacent vacancy.
1.2- + Si AT 0°K
CB
O.lleV “=
l.O- -f- t
0.44eV
E 0.8-
2 V-
g 0.6-
5
0.4 r
Figure 11: Estimated vacancy energy levels in the silicon band-gap at 0°K.
Diffusion and Ion Implantation 471
I- T=14003K
6- 6
I-
4-
(a) (b)
Figure 12: Calculated changes in the ratios of ionized to neutral vacancies at (a)
300°K and (b) 1400°K.
472 Semiconductor Materials
I l111111
n.I
LOG(n) LLECTRONS/cm3
can influence surface bonding which also affects surface quality. As these
processes produce point defects it is possible that extended structural
defects can grow in the silicon. Point defects can also influence the
precipitation of oxygen. Oxygen is incorporated into the crystal during
crystal growth, and precipitates during subsequent heat treatments. It is
known that these precipitates create good internal gettering sites for
metal impurities with subsequent impact on junction quality.
Point Defects
Point defects are defined as atomic defects. There are atomic defects
such as metal ions which can diffuse through the lattice as shown without
involving themselves with lattice atoms or vacancies. Another type of
atomic defect is the self-interstitial which in silicon is a silicon atom that is
bonded in a tetrahedral interstitial site. Examples of point defects are
shown in Figure 15.
One of the major controversies in solid state science currently is: what
is the dominant native point defect in silicon - the monovacancy or the
silicon self-intersitial? A brief review of the arguments for each species is
given below
The Monovacancy
From statistical thermodynamics, it is known that if avacancyisformed
by removing an atom from the crystal and depositing it on the surface, the
freeenergyofthecrystalwilldecreaseasthenumberofvacanciescreated
increases until a minimum in this free energy occurs. Because a minimum
inthefreeenergyoccursforacertainvacancyconcentrationinthecrystal,
the vacancy is a stable point defect. Other experimental observations
involving vacancies are listed below:‘O
Self-interstitial I
I
Substitution81 dopents
B,Al, 08 ; P, As, Sb
P-
After reviewing the balance sheet of pros and cons surrounding the
question of the native defect in Si, one isstill left with the question: what is
the native defect responsible for impurity diffusion and defect growth in
Si? So far we only have clues. However, the majority opinion currently is
that both types of point defects are important. Thermal equilibrium concen-
trations of point defects at the melting point are ordersof magnitude lower
in Si than in metals. Therefore, a direct determination of their nature by
Simmons-Balluffi type experiments24 has not been possible.The accuracy
of calculated formation and migration enthalpies appears to be within +
1 eV but do not help to distinguish whether vacancies or interstitials are
dominant in diffusion. The interpretation of low temperature experiments
on the migration of irradiation-induced point defects is complicated by the
occurrence of radiation-induced migration of self-interstitials.25*26 In addition,
thereareindicationsthatthestructureandpropertiesof pointdefectsmay
change from low to high temperatures. 27The observation of extrinsic type
dislocation loops in dislocation-free, float zcne Si showed that self-inter-
stitials must have been present in appreciable concentrations at high
temperature during or after crystal growth.28~2Q However, it is unclear
476 Semiconductor Materials
D, = Df + D; C D,= + D; (39)
where Di is the measured diffusivity and DiX, D!-, Di= and Di+
are the intrinsic diffusivities of the species through Interactions with
vacancies in the neutral, single acceptor, double acceptor ordonorcharge
states respectively. These individual contributions to the total measured
diffusivity were described in a previous section.
Analogous to the vacancy model, Si self-interstitials can be assumed
to be dominant such that C,>>C,. For such a model, dopant and self-diffusion
are assumed to occurviaan interstitialcy mechanism.30 Mobile complexes
consisting of self-interstitials in various charge states and impurities are
assumed to exist.
In principle, both vacancies and self-interstitials may occur simulta-
neously, and somewhat independently. Indeed, any relationship that may
exist between Cvand C, may be dominated by the Si surface which can act
asasourceorsinkforeitherspecies. If alocaldynamicalequilibriumexists
between recombination and spontaneous bulkgeneration, vacancies and
self-interstitials would react according to
V+IZ 0, (40)
where 0 denotes the undisturbed lattice. The law of mass action under
equilibrium for this reaction is
c, cv = cp cc. (41)
Di = D/ + Dt”e (42)
Diffusion and ton Implantation 477
Dsl = f, D, c + f, D, C; (44)
Experimental Observations
In order to understand whether vacancies or self-interstitials are
000000
op3 0
000000
Figure 16: A schematic diagram of the Watkin’s Replacement Mechanism.
478 Semiconductor Materials
Questions:
0 I+v~o
0 c,cv = cpc;
l Time to equilibrium?
In D
Figure 17: A diagram of the spectrum of the vacancy vs. self-interstitial debate.
Table 2
Oxidation-Enhanced Diffusion
As it was mentioned in the above discussion, oxidation generally
enhances the diffusion of Group III and Group V elements except for
antimony. These results are summarized in the Figure 19. In this figure
oxidation- enhanced diffusion is generally observed by depositing a silicon
nitride mask on the silicon surface which will prohibit oxidation in the
regions that it covers. Then oxidation is performed in a window opened to
the silicon surface, so that differential changes in junction depth can be
observed. In orderto explain these results, HUESproposed a model whose
essential points are:
PARTIAL
/+ DISLOCATION
Figure 18: A model of self-interstitial diffusion from the bulk to the partial dis-
location bounding a stacking fault. Under non-oxidizing conditions the concen-
tration of self-interstitials at the fault line, Cr L, is greater than the equilibrium
bulk interstitial concentration, Cr 0. Under oxidizing conditions, Cr is greater
than CrL until the retrogrowth temperature is reached.
fi _ D!/DY
I 1 (46)
we can write
Experiments that use the backside of the silicon wafer to inject self-
interstitials and thus observe diffusion on the frontside of the wafer are
illustrated in Figure 20.42j43On the wafersurface, films of Si,N, orSi,N, on
SiO, are deposited over previously diffused layers. On the backside a
windowisopenedwhosedistancefromthefrontsidesurfacecan bevaried
by etching. It can be seen in the figure that the backside oxidation can, in
fact, influence the diffusion of dopantson thefrontsidesurface.Theratioof
the junction depth under the oxidized portion to the non-oxidized portion
versus distance from the backside oxidizing surface is shown in Figure 21
for boron-phospohorus and antimony. These results were obtained in float
zone (FZ) silicon with no oxygen incorporated in the silicon. It can be seen
that self-interstitial diffusion lengths of the order of 200-300 microns were
obtained. These backside oxidation experiments show:
BN-arra BOsrea
‘JBN ‘JBO
Figure 20: Experiments illustrating the use of the backside of the silicon wafer
to inject self-interstitials in order to observe diffusion on the frontside of the
wafer (after Mizuo).
Figure 21: The ratio of the junction depth under the oxidized portion to the
non-oxidized portion of the wafer vs. distance from the backside oxidizing
surface.
’ +--A 00 CALCULAlEO
. .
DATA FROM TANIGUCHI tlal
hl-_
Figure 22: Measured and calculated values of boron and phosphorus diffusivity
as a function of total impurity doping. Data are divided into contributions to
substitutional impurity diffusion under non-oxidizing conditions, Dsl and the
enhanced contribution due to oxidation ADO.
Diffusion and Ion Implantation 485
1 +v -0. (50)
(HCIoxidation)
-
A {- OSF shrink
Figure 23: Diagram of SiCl formation during oxidation with the subsequent in-
jection of vacancies. The vacancy injection reduces the concentration of self-
interstitials in the bulk and causes oxidation stacking faults to shrink.
486 Semiconductor Materials
1 I I I I lllll I I I IllIll
10-l 100 10
PERCENT HCI IN 02
Figure 24: The effect of adding HCI to O2 on stacking fault length after oxida-
tion of silicon.
Silicon self-diffusion data over the range 850 to 1380°C are shown in
Figure 25.35,48,4gThe high temperature data show an activation energy of
5.02 eV while the lower temperature show a 4.25 eV energy. Watkins and
Corbett50reportedanactivationenergyforselfdiffusioninSiof3.9eV.This
result was obtained from low temperature annealing of E-centers, impurity-
vacancy pairs at 100” K. The cause for the continual decrease in activation
energy with decreasing temperture has been ascribed to different charge
state vacancies dominating self-diffusion in the various temperature
ranges.13~51~52ThusatverylowtemperaturesneutraIvacanciesmaydominate
self-diffusion. At high temperatures, both donor and acceptor vacancy
diffusion was considered important. An alternate view was expressed by
Seeger and Chik35 who suggested that in Si at low temperatures, self-
diffusion mainly occurs via vacancies, whereas at high temperatures it is
dominated by the interstitialcy mechanism. Their observations indicate a
change in theself diffusion mechanism and/or in theentropyand enthalpy
of self diffusion as a function of temperature. The change in entropy with
temperature can be accounted for by assuming that the form of the self-
interstitial changes with temperature. For example, the entropy would
increase due to a spreading out of the self-interstitial over several atomic
volumes.
Diffusion and Ion Implantation 487
TEMPERATURE
(‘C)
10-12 1300 1200 1100 1000 900
10-13
CALCULATED
(As-DOPED
TO8~1O~~,nr-~)
10-l’
10-15
2 1o-16 CALCULATED
(B-DOPEDTO
Z 2.5~10~~~11t’~
)-QSi=4.78CV
cy‘
E 10’17
.-
k?
10-18 . INTRINSICSi (HETTICH,
Ni IN INTRINSIC Si(SEEGER
6 CHIK)
B-DOPED Si (HETTICH, etaI)
10-20 A As-DOPED Si(FAIRFIELD 8 MASTERS)
. P-DOPED Si(SANDERS 8 DOBSON)
6 7 8 9 10 11
lO’/T (‘K-l)
Figure 25: Selfdiffusion data in intrinsic and heavily doped n and p-typesilicon.
[V-l n iv+&
-=-=-
[v-lj ni IV+] * (52)
Using the relation Dsi = 1/(2n,)D,[V), where the one-half term is the
correlation factor for the diamond lattice and nH is the number of lattice
sites, the Si self diffusion coefficient becomes
(53)
The expression for Dsi= cannot be obtained by analyzing the data in Figure
25. The values of the activation energies for each term are consistent with
theformationenthalpies, migration enthalpiesandaveragefreeionization
energies associated with each vacancy.53 Thus, at temperatures below
600°C Equation 53 predicts that C?siapproaches 3.89eV as the neutral
vacancy dominates self-diffusion. This agrees with the low temperature
value observed by Watkins and Corbett.50
into the lattice, and thiswould be the rate controlling mechanism as it is for
self-diffusion. Therefore, the difference between the activation energy for
self-diffusion and Group III or V impurity diffusion is less than the impurity-
vacancy pair binding energy, E,. This is illustrated in Figure 28 where a
particular long-range interaction potential is assumed. The potential energy
between a vacancy at a third coordination site and one infinitely removed
from the impurity atom is AQ. Thus the activation energy for impurity
diffusion, Q,, is proposed to be
\
1
012345678
VACANCY SITE ON COORDINATION
SHELL ABOUT THE IMPURITY
Figure 26: A schematic diagram showing a long range vacancy-impurity interac-
tion potential which could account for the lower activation energy of impurity
diffusion compared with selfdiffusion in silicon.
490 Semiconductor Materials
7, (INTERSTITIAL
--I--- fLUX)
-t_ ‘Si
I
C F’B
X
Figure 27: A schematic diagram showing the expected redistribution of a buried
layer under the influence of a flux of self-interstitials chemically pumped from
the surface by a phosphorus diffusion.
Diffusion and Ion Implantation 491
(c)
dopant atom-point defect pair.37 From our previous discussion the chemical
pumping effect is unlikely. Thus the question remains: which type of point
defect ismorelikelytopairupwiththedopantatoms?Theanswerappears
to depend in an unpredictable way on the type of dopant atom. For that
reason, the following sections will describe what is known about each of
the important dopants in Si.
Property Result
TEMPERATURE ("C)
lo-l2 ARSENIC
'\A EXTRINSICDATA
IO-l3
u
,”
2
E
”
z o GHOSTAGORE
10-17
;
10-18I
6 6.5 7 7.5 8 8.5 9
10-11,
E . 1ooo’c
0 02
{ ’ N2
A N2
10500 A 02
10-12 C0 WET 02
1lOOC 0 02
t
Figure 30: Arsenic diffusivity vs. total arsenic concentration. The solid curves
are calculated. The data represent diffusions in various furnace ambients.
where the formation of V-As+ and WAS+ pairs are included. By writing
down the equilibrium reactionsfortheformation of eachspecies in Equation
57 and setting C,, = n, then C, becomes
where K(T)‘s are the equilibrium constants Equation ,58 can be used to
describe the data in Figure 32.
Diffusion and ion Implantation 495
Dh e -J tac+X*, (59)
which is Fick’s First Law. If J is the flux of monatomic As, C,,, with a linear
concentration-dependent diffusivity, then
CA, ecA5
J = -Di - -
ni ax (60)
(61)
Di (nh$
DA =
4K&Th3
1+ (62)
1 + 3K&Tb2
102’ c
F
9oo"c
2020 MIN.
.-CT
DEPTH(Nm)
Figure 31: Total arsenic CT and free electron profiles in silicon following an
arsenic diffusion.
496 Semiconductor Materials
: IMPLANTED~DIFFUSED LlMIT(1050 cj
uE
-Z CHEMICAL SOURCE
10::
lx1016cm.3 (1050 Cj
I I I lllll I I111111
102” lo’!
CT (cm”)
Figure 32: Electron concentration vs. total arsenic for chemical source arsenic
diffusion and diffusion of an ion-implanted arsenic layer with the same inte-
grated concentration.
Y
v).
“5 10-14
xl o KENNEDY AND MURLEV
d
1 = 1000T.
CTO=8.5x102~ATOMS/cm3
10-15
10’9 1020 1021
C+ATOMS/cm3)
Figure 33: The effect of diffusivitv of arsenic vs. total concentration for dif-
fuiions into p-type silicon at 1000%.
CONCENTRATION (cnT3)
Figure 34: Arsenic and arsenic-vacancy pair concentrations vs. electron concen-
trations at two temperatures.
8
Microlithography for VLSI
R. Fabian Pease
Stanford University
Stanford, California
INTRODUCTION
( i) Spin coating the workpiece with a thin (0.1-2 urn) film of resist
541
542 Semiconductor Materials
Both negative and positive electron beam resist are widely used for
mask making. Because the image is built up by scanning a focused beam
sequentially over each address, the process is slow. Therefore, there is
considerable emphasis on resist sensitivity which is usually, unfortunately,
quantified as the dose in C/cm2 required to bring about the desired
chemical change.
Negative resists are long chain polymers containing groups which, on
electron bombardment, form crosslinks between adjacent chains.’ Such
groups include:
Those portions of the film that are crosslinked are insoluble in a suitably
chosen developer.Typical required doses range from lO.‘to 10-5C/cm2for
10 keV electrons.
Positive electron resists are also long-chain polymers. However, positive
resists contain groups which, on electron bombardment, cause chain
scission and hence locally reduced molecular weight. The exposed regions
are now selectively soluble in a suitably chosen developer. Examples of
such materialsarepoly(methyl methacrylate)(PMMA)*and poly(butene-1
sulfone) (PBS)3. The chain scission process also results in gas evolution
which may also promote the selective solubility of the exposed regions.
Microlithography for VLSI 543
Typical required doses are somewhat higher than those for negative
acting resist and range from 1O+ to 1OW4 C/cm2 for 10 keV electrons.
The resists as purchased are solutions of the above polymers in
volatile solvents. A few ml. of the solution are dispensed onto the center of
thewaferwhichisspunatapre-determinedratewhilethesolventevaporates
to leavea polymerfilm of predictable thickness.The theoryofthespinning
action is quite complicated because some material is lost through cen-
trifugal action and the viscosity of the remaining material changes during
evaporation of the solvent.4 Often the thickness, t, of the remaining film is
given by the expression:
t l klJ/&
SYSTEM
CHIP PATTERN
WRITING PATH
(Table mollon)
BEAM BLANKING
--------em
, LASER
BEAM DEFLECTION
I INTERFEROMETER _--__----- -_
I
WRITING TABLE I
, I
SERVO e .- J
MOTOR I
I
I
MAG 7
INPUT
Figure 1: (a) Writing strategy used in Bell Laboratories Electron Beam Exposure
System “EBES”: a single stripe of the circuit pattern is read out repeatedly from
memory then is written on to each chip site. (b) Schematic view of EBES. The
feedback system employs beam deflection to compensate for table position errors
(After D.R. Herriott et al.‘).
Microlithography for VLSI 545
ELECTRON LENS
BLANKING ELECTRODES
DEFLECTION COILS
- FINAL LENS
APERTURE
BEAM OF SEMI-ANGLE
OF CONVERGENCE a
WORKPIECE
Theelectronsenteringtheresistfilm notonlybringaboutthechemical
effectsalluded toearlier, butarealsoscattered bythefilmandsubstrateso
that the lateral extent of the interaction can exceed the diameter of the
impinging beam. This effect is known as the proximity effect and has been
the subject of large numbers of papers. The effect lends itself well to
546 Semiconductor Materials
ELECTRON-ELECTRON
INTERACTIONS
(I (rad)
Figure 3: Contribution to focused electron beam diameter as a function of semi
angle of convergence a. Often the total beam diameter is estimated by adding in
quadrature the contributing diameters. The electron energy, E, is IO keV with
energy spread AE = 3 eV, chromatic aberration coefficient 4 cm, spherical aber-
ration coefficient CS = IO ems, cathode current density 10 A/cm’, and column
length 66 ems.
r2
J(r) + Jb) exp (- 2)
='p
where Jr is the current density at radius r, and rP= 0.2 pm + 0.1 pm for the
case depicted in Figure 4a. At higher voltages two gaussian terms are
needed for an adequate description.
The chemical effect of electron bombardment has already been des-
cribed as generating crosslinks between adjacent polymer chains to
create an insoluble gel. We can invoke a simple model to describe quanti-
tatively the resist’s behavior in terms of the parameters, sensitivity and
contrast.
The starting point of the model is to assume the following:
X - brn X -- pm
T, l 1 - exp(-v, Ag Q/q)
Therefore, a plot of Tn versus log,, Q can be drawn (Figure 5). The general
form is clear and understandable. At doses below the critical dose Tn = 0.
Microlithography for VLSI 549
0.8 -
0.7x1 .6x10-lg
Q-, vm k3 C/cm:!
This looks quite reasonable since when the required dose is reduced for
largerV,,afilm madeof largermolecules requiresasmallerconcentration
of crosslinkstobondthemalltogether.Also,at IargevaluesofAtheenergy
dissipated per incident electron per unit depth, and g, the number of
crosslinks per eV dissipated will also tend to lower Q.
The other parameter is the contrast, or gamma, which is a measure of
the minimum ratio of exposures needed to bring about a required difference
in T, (often 0 to 50%). If we approximate the curve if Figure 5 as three
straight line segments(shown dotted) then theslopeof thecentral portion
is the contrast (y) defined as dT,/d(log,, Q). If this slope is equal to the
maximum slope of the solid curve then we can quantify y as
dTn
d(loglO)Q I l -
eloglOe +085 l
ElX
10,000 8
whence A + + 10 eV/cm/electron
3x0.4x10-4
We must now expose the resist on the wafer in accordance with the
maskpattern.Thesimplestwayofaccomplishingthisisbycontactprinting
or proximity printing. The formertechnique allows excellent resolution but
is prone to defects especially when many contacts are needed to assure
accurate alignment. Thus, contact printing is not in widespread use for
552 Semiconductor Materials
Microlithography for VLSI 553
I I
+=
P
f \ \r =l6
a
2 \
: \
\
\
zY 1 \
::
i
I-
-l 0.5 -
2
0
F I
s:
u
IL 10-8 10-7 10-6 10-5
EXPOSURE _ C/cm2
Figure 9: Two views of chromium masks made with electron beam exposure of
PBS resist (photographs courtesy of D.H. Dameron and J.P. Ballantyne).
Microlithography for VLSI 555
LOEJECT
Figure IO: Schematic view of doubly reflecting optics used in scanning pro-
jection printing.
556 Semiconductor Materials
because the image of the mask pattern occupies only a small portion (1 x 1
cm*) of the wafer area and the whole area is filled up by stepping the wafer
mechanically and repeating the exposures. Refracting focusing is usually
used for the optics. With this technique, VLSI chips with feature sizes of 1
l/4 pm are being manufactured. The specifications of two such steppers
are shown in Tables 3 and 4. With these techniques both positive and
negative photoresist are used although the former is becoming the most
popular.
In both the scanning projection and the stepping aligners the physics
of image formation is the same. An outline of the basic scheme is shown in
Figure 11. A mask is illuminated and the projection optics focuses an
image of the mask at the wafer surface. As with electron beam lithography
Throughput: 40 1% mm wafer/hour
Numerical Aperture: 0 JO
Magnification: l/5 x
FOCUSING SYSTEM
MASK (MIRROR OR LENS) WAFER
PATTERN --
/
/
/
/
/
/
/
5
RAYS APERTURE IMAGE OF MASK
DEFINING Q PATTERN OF
MAGNIFICATION M
Paul S. Ho
IBM Thomas J. Watson Research Center
Yorktown Heights, NY
INTRODUCTION
575
576 Semiconductor Materials
is
Figure 1: Progress in silicon chip technology since 1960 (Source: Siemens). Note
that the degree of integration quadruples about every 3 years. While the develop-
ment in microprocessors may show after the 32 bit level, the advance in dynamic
memories remains constant (Reference 1).
Protective
PbSn
solder pad Cu-Sn mrermetallic
Silicon 0. I5 pm Cr.Cr,O)
degree of flexibility for the functional design of the chip, to enhance the
performance and the level of integration of the whole system, new packaging
structuresare required which can utilize the high density and performance
of the device chips. This has brought forth significant improvements in the
performance and level of integration of the packaging system. As indicated
by the statistics in Figure 3, the wiring density in packaging has grown
exponentially with time, but at a rate less steep than that of the device
density.4 Although this has not been as well recognized in the past, it has
becomeclear recently that packaging isan important issue, particularlyfor
the high end computer systems where performance is the prevailing
factor.
Thiscan beillustratedbycomparingtheIBM3033and3081 computer
systems.5As shown in Figure4, the performance of the central processing
unit, as measured by system cycle time, can be divided into chip and
packaging portions. While the cycling time in the chip has been improved
about 20%, the improvement in the packaging portion is about threefold.
(The improvement in the circuit chips should not be measured by speed
alone since there issignificant enhancement in the circuit densitywhich is
not shown in Figure 4). The improvement is achieved primarily through a
500-
200-
loo-
20-
13 10-
z
c
$ 5-
9
B
P 2-
3
f?
1 I I I I I I
o-1 1 10 IO2 IO3 IO4 lo5
Figure 3: Progress in circuit density with time at the packaging module level.
The abbreviations associated with the data points represent different versions of
IBM packaging modules (Reference 4).
Metallization for VLSI Interconnect 579
Card
Module
Figure 4: Comparison of system performance of the IBM 3033 and 3081 cen-
tral processing units (Reference 14).
580 Semiconductor Materials
Module - 41.9
Card 38.8 -
100% 100%
3033 Technology
Module-to-card 22,560
670 (no card)
Card-to-hoard 4,000
WIRING STRUCTURE
P = aNb (1)
% I
I I
100
GATES/CHIP
I
1000
Figure 5: Results of a numerical analysis showing the “Rent’s rule” for a logic
chip. This rule correlates empirically the number of wire tracks required per cell
to the device density on a chip (Reference 10).
it
I
;i
2
E
k-
5
104-
0
dOTOTAl_
)///-
O
0
AVERAGE /
CIRCUITS/CHIP
Figure 6: The average length of interconnect and the total wire length on a logic
chip as a function of device density. The lengths are measured in circuit pitches,
i.e. the square root of the area per circuit (Reference 12).
packaging can become very complicated and usually requires a high level
of integration and optimization. This can be illustrated by the packaging of
the IBM 3081 system. This system contains two main levels of wiring
structure: the module and the board. The module is designated as the
thermal conduction module (TCM). Each TCM contains about 100 device
chips, each of which has about 2000 circuits. It is built into a ceramic
substrateofabout5.5 mmthickand 10cm bylOcmsquare(Figure8).The
wiring structure in the board is designed for each board to support 10
modules.15 It has a 20-layer structure containing about 6,000 connectins
with dimensions of 60 cm by 70 cm. The wiring complexity of this board is
comparable to that of the TCM.
It is interesting toestimatethewiring requirementsforthemoduleand
theboardonthebasisofRent’srule.Takingtheparameteratobe2.5andb
0.6, the number of interconnects required for each chip is about 250
connectionsforthe2000circuits.Tosupportthe100chipsononemodule
requires 16 times the I/O connections of an individual chip. This turns out
to be about 4,000 connections for each module. These requirements
evolve the module into a structure with more than 30 interconnect layers
Metallization for VLSI interconnect 585
B;IW phe
- ,
Figure 8: (a) Exploded view of the thermal conduction module assembly. (b)
Schematic drawing showing the wiring structure of the multilayer ceramic sub-
strate used in the thermal conduction module (Reference 14).
586 Semiconductor Materials
. _ _ _ _ _ _ _ _ _ _ _ _ _
--.
,------ 1
i .._...... 1
r 50
. _ ___
-.-.
--=.
-.-.
_
] :--+ : : : : i -- (
0
0 2 4 6 8 10
1 3 5 7 9
Channel position
Figure 9: Histogram of the maximum wire densities within a given channel cal-
culated using various methods of interconnect routing. The channel position
provides a measure of the length of the interconnect. The top dash line repre-
sents the results obtained by the simulated annealing technique (Reference 16).
Metallization for VLSI Interconnect 587
Scaling Factor
Doping concentration k k
Voltage I/k -1
Line resistance k k
ELECTRICAL CHARACTERISTICS
C = CWL
t OX
and (2)
wherep is the resistivity, L the wire length, W the wire width, t, thickness of
the wire, E the permittivity, and t,, the thickness of the oxide. Equation 2
showsthatift,andt,,canbescaledinthesamemannerasL,thelinedelay
would remain constant (seeTable 2 also). Therefore, the line delay becomes
an increasingly larger portion of the total circuit delay as the device
dimension decreases. Eventually it can become a substantial part or even
dominatethesystem response timeforverysmalldimensions. In practice,
linearscalingisdifficulttoaccomplishfordeviceswithsubmicrondimensions
sincet,andt,~areusuallylimitedtoabout0.5ymand1OOArespectively.In
addition, L of all the interconnects does not scale uniformly as we discussed
previously in the wiring placement section.
In general, the effect of the line delay will become a problem when the
minimum dimension reaches below about 2 pm. The impact is less for FET
memory circuits than for bipolar logic circuits because of the bipolar
circuits more complex wiring structure. For the 3-level bipolar chip with
about 2 pm minimum dimension shown in Figure 2, the wiring delay
constitutesasignificant portionofthesystem processing time. Inaddition,
for VLSI applications, the chip dimension is usually enlarged in order to
accommodate the high device density which, when combined with a more
complx wiring structure, results in broadening of the overall length distri-
bution of the interconnects. This widens the distribution in the RC time
constants with substantial increases for the portion of the wires with long
lengths.
To minimize the RC delay, it is usually important to use materials with
low resistivityand permittivityto build the interconnectstructure. InTables
4 and 5 are summarized some of the physical propertiesforthe commonly
used conductor and insulator materials. *‘a** (The resistivities for thin films
of these metals are not given here since they depend in general on various
parameters, e.g. method of deposition and grain structure. They are about
20-40% higher for pure and large grain films of Al and noble metals but can
be 2 to 3 times higher for impure refractory metal films.) Because of their
excellent conductivities, it isclearwhy Al, Au and Cu are the most commonly
used metals. However, in certain applications, because of processing
requirements (e.g. the annealing temperature) or device design (e.g., the
high-densityself-aligned polycidegate), materialsof lesserconductivities,
such as refractory metals, silicides and even highly doped polycrystalline
silicon are used.
Metallization for VLSI Interconnect 591
Electrical
5.5 4.4 13.0 0.004
Porcelain
Glass-
Ceramic [7]
For the insulators, SiO, is the universal material used to form the
dielectric layer on circuit chips. It has a low dielectric constant of 3.5 and
can be produced with extremely low defect density by oxidizing the Si
substrate to a thickness as small as l OO-2OOA. This makes it well suited for
gate insulator applications although forsubmicron devices, there is some
question regarding the integrity of SiO,asagate insulatorforthicknesses
below 100A.23 For interlevel insulation, silicon oxide up to 2 ,um thick
produced by sputtering or evaporation (often not of the exact SiO, stoichio-
metry) is often employed. Si,N, is frequently used in combination with SiO,
in spite of its high dielectric constants, Its excellent mechanical strength
makes it well suited to serve as a lithographic masking material. Examples
of the Si,N,/SiO, combined layer can be seen in Figure 2.
Forpackagingapplicationsceramicsformed byvariouscombinations
of oxides, particularly Al oxides and Si oxides, are common materials. For
example, the multilayerceramic module shown in Figure8 employs materials
of several oxide mixtures. The ceramic materials have relatively high
dielectric constants (about 7-8) and have to be processed at elevated
temperatures (above 1500°C). Special, and often complex, processes
have to be developed for the application of this type of material forforming
multilayer structures. For example, high temperature processing neces-
sitates the use of metals with high melting points, such as the refractory
metals. This increases the response time of the system due to the high
resistivity of these metals. To circumvent these difficulties, polymeric
materials are being considered to replace the ceramics in chips as well as
in packaging.24 The main advantages of this class of materials are the low
dielectric constant (about 3.5, similar to SiO,) and the low processing
temperature (usually below 400°C). High-temperature polymers such as
Polyimides are used to satisfy processing requirements where thermal
stability up to 400°C is required.
Some of the problems relating to wiring delays for VLSI applications
have been discussed by McGreivy. 25 He has considered the change of the
access time for a static NMOS (N channel) RAM with decreasing design
rules. His results are shown in Figure 10. The access time decreases
continuously with shrinking device dimensions down to about 1.5 to 2 pm.
Below that, the effect due to the RC delay of the interconnect becomes
observable and its magnitude depends on the resistivity of the material
used. For refractory metal gates with sheet resistivities of 1 ohm per
squarecm(sheet resistivityequalsp/t,),adecreasein theaccesstimecan
still be achieved below 1.5 pm although the gain is very small. The access
time is doubled what one would expect to achieve in an ideal scaling
model. For resistive polysilicon gates with 20 ohm per square cm sheet
resistivity, corresponding to a 1 pm thick gate with 200 ohm-cm resistivity.
The access time increases with decreasing geometry due to the RC delay.
Theaccesstimebecomesaboutanorderofmagnitudemorethantheideal
case in the submicron range.
He has also considered the problem of parasitic capacitance. His
results for the variation of the capacitance components of interconnect as
a function of line width are shown in Figure 1 1. Of the three capacitance
components, only the metal-to-substrate capacitance, C,,, decreases
Metallization for VLSI Interconnect 593
Figure 10: Variation in the access time of a 4k NMOS RAM with decreasing de-
sign rule for gate interconnects with different resistivity (Reference 25).
t-wm
+cw’ -+wm--.I
I
I
tfox
SiOz
Substrate
-C
7 me
-C
7 ms
2.0
I .5
C apacitance
(Relative
to c
ms
at 1 micron)
I .o
0.5
w m' ws
(microns)
MATERIAL REACTION
Fi = -vpi
Fi = yVci-v{&))-oiei-qiEi
where E, and Ei represent the deformation strain and the electric field
respectively.
Inthisform,thedrivingforceinamulticomponentsystemcanoriginate
596 Semiconductor Materials
fromthreetypesofsources.Thefirstsourceisrelatedtotheconcentration
gradient which is generally recognized as the diffusion term. The second
force comes from the internal chemical energy gradient representing the
driving force associated with the change in the chemical form of the ith
element, e.g. the compound phase or the composition in a concentrated
alloy. The last type relates to external constraints such as an applied stress
u, or an electric field E,.
Combining Equations3 and5, it isclearthat the existence of an atomic
flux requires not only that the atoms move, i.e., D,#O, but also that there is a
nonvanishingdrivingforce.Thustheroleofthedrivingforceisasimportant
as the diffusivity although the present discussion will address mainly the
diffusivity issue since the variety of driving forces makes a systematic
discussion of the role of the forces difficult. (One exception is the later
discussion on electromigration.)
The existence of an atomic flux by itself is not sufficient to cause
damageformation. Inorderforthattooccur,alocaldepletionoraccumulation
of materials is required. This condition can be expressed by the flux
continuity equation as
dci Ci.-,P
-V.Ji + +
-Z=
Tm/T
Figure 12: Summary of diffusivities via various types of structural defects. The
temperature scale is normalized to the absolute melting temperature T, (Ref-
erence 29).
6c T/T,,,=0.6,-----7 1 c T/Tm=O.5,__-_, 1
-6
0 2 4 6 8 IO 12 0 2 4 6 8 IO 12
Figure 13: Regimes of grain size (gs.) and dislocation density Pd over which lat-
tice diffusion (I), grain boundary diffusion (b), or dislocation diffusion (d) is
dominant. The calculation is based on steady-state diffusion through a thin film
specimen of an fee metal as a function of the homologous temperature (T/T,)
(Reference 27).
598 Semiconductor Materials
METALLIZATION RELIABILITY
(Ni-0.5Co)Au”(n)
\\
I I I I I
1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5
Tm 1 T(K)
Figure 14: Plot of the grain boundary diffusivity against the reciprocal normal-
ized temperature T,/T for (a) data obtained by sectioning techniques in thin
films and some bulk materials and (b) data obtained by permeation techniques
in thin films (Reference 30).
600 Semiconductor Materials
300 K
- THEORY
. PISl -s,
-!- ( ,O”Ocm”* )
A
Figure 15: Theoretical and experimental values of specific contact resistance.
Note that for doping level exceeding 1019 per cm3, R, is dominated by the tun-
neling process while for doping level of 1017 per cm3, thermionic emission domi-
nates and R, becomes constant (Reference 32).
Cosputtered Alloy 25
II 45-50
HISi,
VSiz 50-55
NbSiz 50
II 1000 -70
WSiz
8, >lOOO
FeSi, 700
SCHEMATIC PRESENTATIONOF
AI/Pd$i/Si REACTION
Pd+Si
a) INITIAL STAGE
AI/Pd
COMPOUND
IZFFU A13Pd4Si
Pd2 Si
‘:nSI ’ / l’ 7
::,
b) INTERMEDIATE STAGE
c) FINAL STAGE
BURIED LAYER
(b)
Figure 17: Schematic cross sections of silicide contacts to (a) bipolar and (b)
MOS devices. Note also in (b) the combined use of silicide and polysilicon for
gate metallurgy (from Reference 43).
widelyaccepted.47OnecancomparetheresistivitiesofsilicidesinTableVI
to those of metals in Table4. In Figure 17b, an example of a polyicide gate
used in a recent MOS device is shown.
Electromigration
Electromigration describes the movement of atoms in a metallic
conductor induced by the passage of a direct current. Its magnitude is
determinedbytheatomicdiffusivityandthecurrentdensity. Electromigration
induceddamageintheformofopensorshortsintheinterconnectlinesisa
result of a local divergence of the mass flux. This divergence can be
generated by various types of inhomogeneities, such as those from grain
size variation local heating and stress gradients.48
With regard to electromigration, the main impact of scaling is to
increase the current requirements of interconnecting lines. This problem
has two basic aspects, onefrom the increase in the current density and the
other from the reduction in device dimensions.4g From Table 2, the current
density is seen to increase linearly with the scaling of MOS devices and
Metallization for VLSI Interconnect 605
more than linearly for bipolar devices. This increases the driving force for
electromigration as well as the Joule heating generated in the conductor.
With the heating increase asj*p, the effect can be significantly higherthan
the driving force coming from the linear increase in j. The combination of
these factors can raise the conductor temperature giving rise to higher
atomic diffusivity and electromigration flux. This problem, together with
the increased powerdensity, necessitates an improvement in heat dissipa-
tion during device operation. The combined effect of these factors will
inevitably cause the electromigration flux to increase beyond thatcaused
by the increase in current density alone.
The increase in the current density can be estimated based on the
trend in device dimensions. For present devices with dimensions of about
3 pm or larger, the current density can reach 2 X 1O4and 5 X 1 04A/cm2for
MOS and bipolar devices, respectively. For the next generation of devices
with minimum line dimensions in the range of 1.5 to 3 pm, j increases to
about 5-l 0 X 1 O4 for MOS and 2-4 X 1O5 for bipolars devices. This trend
continues and can result in j exceeding lo5 and lo6 respectively, for
submicron MOSand bipolardevices.Thisisan increaseof 102-103timesin
the current density. Since an isolated metal wire can carry only about 1 O4
A/cm2 before melting, the Joule heating generated in a line carrying 1 O5
A/cm2 must be almost completely removed through the substrate and/or
the passivating overlayer. When the current reaches ~10~ A/cm2, any
imperfection in the substrate, such as processing defects or interfacial
barriers, can cause thermal runaway todestroythe line. Even without such
a catastrophe, the heating effects of such high current densities will
increase the rate of electromigration, resulting in a significant reduction in
the lifetime. In practice, this is reflected by an increase in the exponent n in
the lifetime equation of t,,=Aj-” exp (AH/IV)outside the normal range
between 1 and 2. This effect adds considerable difficulty in extrapolating
the lifetime under operating conditions for submicron lines from results
obtained in accelerated stress tests50
The otherelectromigration problem duetosize reduction isgeometry-
related and caused by scaling into the submicron range. For metal films, it
is generally observed that the grain size is about the same as the film
thickness. For a 1 p thick film, the common thickness of interconnecting
line, there will be only a few grains spanning across a l-3 pm wide line. At
the device operating temperature, the electromigration flux is confined to
grain boundaries. With a small number of the grains across the line, each
individual divergent site in the grain structure becomes potentially more
damaging since a line can fail without requiring a statistical linkage of
several divergent sites, as would be the case of a line many grains across.
This shortens the conductor lifetime while increasing the randomness of
the failure statistics, i.e. increasing the statistical deviation ain the lifetime.
Both trends have been observed in lifetest@’ as well as in computer
simulation for linewidths down to about 2 pm52. This effect is significant
since the extrapolated lifetime for device operation can be significantly
reduced by an increase in u.
Another effect which results from the increase in the grain size-to-line
width ratioisadecreaseintheroleofthegrain boundaryinmasstransport.
606 Semiconductor Materials
SUMMARY
REFERENCES
47. C.Y. Ting, SS. lyer, C.M. Osburn, G.J. Hu and AM. Schweighart in VLSI Science
and Technology/7982, ed. C.J. Dell’Oca and W.M. Bullis, The Electrochemical
Society, 1982, p. 224.
48. F.M. d’Heurle and P.S. Ho in Thin films-lnterdiftusion and Reactions. ed. by
J.M. Poate, K.N. Tu and J.W. Mayer, Wiley Interscience (1978). Chp. 8.
49. P.S. Ho. IEEE Proc. of 20th Sym. Reliab. Phys. San Diego, Ca. (1982). p. 288.
50. P.B. Ghate, IEEE Proc. of 20th Symp. Reliab. Phys. San Diego, Ca. (1982) p.292.
51. G.A. Scoggin. B.N. Agarwala, P. Peressini and A Browillard, Proc. 13th. IEEE
Symp. Reliab. Phys. (1975), p. 155.
52. J.M. Schoen. J. Appl. Phys. 51:513 (1980); K. Nikama. Proc. 19th IEEE Symp.
Reliab. Phys. (198 1). p. 175.
53. H.L. Huang, J. Vat. Sci, Technol. A3:705 (1985).
54. S. Vaidya and AK Sinha. IEEE Proc. 20th Sym. Reliab. Phys. San Diego, Ca.
(1982), p. 50.
55. S. Vaidya and A.K Sinha. Thin Solid Films 75:253 (1981).
56. J.K. Howard, J.F. White and P.S. Ho. J. Appl. Phys. 49:4083 (1978).
10
Characterization of Semiconductor Materials
Gary E. McGuire
Tektroni& incorporated
Beaverton, Oregon
610
Characterization of Semiconductor Materials 611
PHOTOELECTRON
AUGER ELECTRON
OR
IONIZING ELECTRON -(KQL2,3)
involved. For example, the kinetic energy of the Auger electron illustrated
in Figure 1 is typically described as
E = E, - E,, - E, - r$
KL1L2.3 213
where EK is the energy of the ionized core level, Et, .is the energy of the level
from which the electron originates to fill the inmal core hole, E~23 is the
energy level from which the Auger electron originates and C#Iis the work
function. Multiplecharacteristic Augertransitions may beobserved due to
the various core energy levels available for photoexcitation and the multiple
combination of energy levels available for de-excitation and Auger emission.
TheAugerelectronisusuallydescribed bythethreeenergylevelsinvolved
in its emission. Figure 1 depicts the KL, L,, Auger transition.
The kinetic energy of the AE is independent of the excitation source.
As a result, the tendency has been to use electron beams in the l-20 KV
potential energy range for excitation. Electron beams are the preferred
excitation source because they can be focused to a small spot size and
deflected to a region of interest on the sample.
The AE transition is characteristically a small feature sitting on a large
background of inelastically scattered electrons.2 The most prominent
feature in the electron spectrum is the contribution due to backscattered
electronsfrom the primary beam.Thedata has been presented historically
in the dN(E)/dE versus E format as a means to enhance the Auger signal.
More recently the data has been presented in the N(E) versus E due to the
availability of computers for background subtraction.
Figure 2 shows a schematic diagram of an Auger spectrometer. The
optics for the primary beam are coincident with the cylindrical mirror
1 10 100 1000
EnergyWI
Figure 3: Plot of electron inelastic mean free path versus energy which illustrates
the shallow sampling depth of the electron spectroscopies.
614 Semiconductor Materials
L Cl
-1 , I I I I I1 I I I I I I Il.1 1 I
200 400 600 800 1000 1200 1400 1600 1800 2000
SPUTTERING TIME
(min)
Figure 5: Auger depth profile of chemically vapor deposited Si3N4 on thermally
grown SiOz over Si.
I I I
L3%,5M4,5
%GaAs) Ga(Ga203)
WI/E
1525
(Ga2(
I I I I
055 1080 1105
KINETIC ENERGY (eV)
Figure 7: SEM image and 0 and S scanning Auger maps of zone refined Fe foil.
The numbers indicate the approximate orientation of the surface normal of the
various grains.
618 Semiconductor Materials
SURFACE
DIFFUSION
Photoelectron Spectroscopy
Photoelectron spectroscopy is a technique which has many similarities
toAES.ThesameenergyleveldiagramusedtodescribetheAugerprocess
may be used to describe the photoelectron process.12 Excitation of the
ionizing photoelectron may be accomplished through the use of a variety
of energetic photons orcharged particles. The primary focus in this text will
be on monochromatic x-ray excitation of photoelectrons (XPS). Use of a
monochromatic excitation source is essential to this spectroscopy, since
the photoelectron’s kinetic energy is directly dependent on the energy of
the excitation source. By knowing the energy of the x-ray with a high
degree of accuracy and measuring the kinetic energy (KE) of the emitted
photoelectron from the relationship:
one can determine the binding energy (BE) of any electron energy level
less than the photon energy.
A variety of electrostatic electron energy analyzers have been pro-
duced commercially. One of the most popular is the cylindrical mirror
analyzer, Figure 9, similar to that used for AES. A two stage, two cylindrical
mirror analyzers in tandem, device is employed to enhance the energy
resolution. An x-ray source, either an Al or Mg anode, mounted in proximity
to the sample is used for excitation. The x-rays flood a broad area of the
sample since they, unlike the electron source in AES, can not be easily
focused. The acceptance angle of the spectrometer determines the area
of analysis, which is typically a few millimeter diameter circle. With adjust-
able aperture slits the sampled area may be reduced to a few hundred
micrometers.
Figure 10 shows a schematic diagram of another common variety of
XPS spectrometer. It employs an x-ray monochrometer to enhance the
x-ray line width, eliminate satellite x-ray lines and focus the x-rays. The x-
Characterization of Semiconductor Materials 619
Precision Analyzer
Figure 9: Schematic diagram of a two stage cylindrical mirror analyzer and x-ray
source used for XPS.
NONOCHRO!CATOR
CRYSTAL
HEKISPHrRICkI.
ELECTRON GUI
Figure 10: Schematic diagram of a XPS system utilizing a bent quartz crystal x-
ray monochrometer in conjunction with an electrostatic lens and hemispherical
analyzer.
620 Semiconductor Materials
rays from an Al anode are allowed to diffract off of a bent quartz crystal
before interacting with the sample. The natural Al x-ray linewidth is approxi-
mately0.9 eVwhile that of the monochromizedsource is approximately0.4
eV.13Thefocusing propertiesofthemonochrometerproduceaspotsizeof
approximately 150 micrometers. l4 Due to the loss in x-ray intensity in
going through the monochrometer, most spectrometers of this type employ
an electrostatic lens to increase the collection efficiency of photoelectrons
going into the hemispherical analyzer and a position sensitive, multiple
array, detector to enhance the count rate.
Photoelectron spectrometers employ ion guns for in-depth profiling
as in AES. Since the area of analysis is much larger than in AES, the ion
beam isdefocused inordertogeneratea uniform ionflux.Thisreducesthe
ion etch rate but does not prevent one from monitoring signal intensity as a
function of ion sputtering time. In addition, many XPS systems have both x-
ray and electron beam sources for combined multi-technique analysis by
XPS and AES.
Figure 11 shows the Ag3d photoelectron spectrum. The trace illustrates
the relative simplicity of the spectra. The spectral features are Gaussian-
like sitting on a low background. The spin-orbit splitting of the energy
levels, in this case the 3d,,, and 3d,,,, is well characterized and easy to
recognize due to the predictable intensity ratios. The spectra are usually
plotted in the N(E) versus BE format even though the energy analysis is of
N(E) versus KE. Each element exhibits a unique set of photoelectron (PE)
transitions corresponding to its atomic energy levels. The PE transitions
are a function of atomic number so that the energy levels of adjacent
elements in the Periodic Table are all shifted to higher binding energy.
As suggested by this unique set of binding energies, XPS is a good
elemental surface analysis technique. Figure 12 shows a spectrum from a
0’
380 375 370 365 360
BE, eV
Figure 11: XPS spectrum of the Ag3d transition showing the spin-orbit splitting
into the 3d512 and 3d312 components.
Characterization of Semiconductor Materials 621
Oh
I I I I I
110 106 102 94
BINDING ENERGY%)
Figure 12: XPS spectrum of a Si surface with the native oxide showing the
chemical shift in the Si2p transition.
clean Si wafer to illustrate this point. Two Si2p transitions are observed,
one for elemental Si and one for the native oxide formed on the wafer as a
result of air exposure. One can get a feel for the surface sensitivity of XPS
since the native oxide thickness is typically less than 30A. The surface
sensitivity, as in AES, is controlled by the inelastic mean free path of the
electron as illustrated in Figure3, rather than the path length of thex-rays
used for excitation.
The ability to distinguish different oxidation states, as in thecase of Si
and SiO,, has been one of the recognized strengths of XPS. These chemical
shifts in the core level binding energies are due to changes in the valence
electrondensityduetocompoundformation.Although,thechemicalshifts
maybeaslargeas10-12eV,thereisfrequentoverlapformanycompounds
as illustrated in Table 1.15 However, there are many sources of chemical
information in the spectra. The sources include first the identification of
the elements present and their relativeconcentrations, then thechemical
shift of the cation to determine its approximate oxidation state and finally
the chemical shift of the anion to determine its oxidation state. Table 2
illustrates the magnitude of the chemical shifts observed for anionic
species X-*, X0,-* and X0,-* when X is S, Se and Te.16 Combination of this
information gives a detailed picture of the chemistry of the surface under
investigation in many cases. There are other spectral features which
provide additional chemical information but which are too detailed for the
scope of this text.
Figure 13 illustrates the use of XPS in the investigation of the anodiza-
tion of GaAs.7 Both the As3d and Ga3d transitions may be observed. By
combining XPS analysis with ion sputtering the composition of the ano-
622 Semiconductor Materials
Chromium, Cr 2% 24
Oxidation AE Relative to
Compound State Elemental State
S Se Te
x-2 -2 -1.4 - 0.8 - 0.6
x0 0 0 0 0
x03-2 +4 3.6 3.7 29
N(E)/E
40 20 0
BINDING ENERGYWI
Figure 13: XPS spectrum of anodized and annealed GaAs showing the chemi-
cally shifted Ga3d and As3d transitions at various depths in the oxide and at the
GaAs substrate.
peaks, one for the oxide and one for GaAs. The oxidation of many compound
semiconductors has been studied by XPS. The composition of the oxide,
especiallyasafunctionofdepth, hasbeenfoundtobestronglydependent
on the method and conditions of formation.
Investigation of metallization schemes used to contact semiconductor
devices is another key area where XPS has been applied. As an example
Figure 14 shows the Pt4f and Si2p spectra obtained from the silicides
which are formed when Pt is used to contact Si.17 Two different silicides
may be formed depending on the annealing conditions used in the process.
The PT 4f,,, and 4f,,, doublet exhibits a chemical shift of less than 1 eV
between the Pt,Si and PtSi phases which may be formed. This chemical
shift is easily detectable. The corresponding Si2p transitions have essen-
tially the same binding energy. The chemical information obtained from
XPScomplimentsthatobtainedfromavarietyofotherthinfilmandsurface
analysis techniques in the investigation of a variety of contact materials.
Since the photoelectron spectra of many elements exhibit only small
chemical shifts for a series of compounds in which the electronegativity
varies over a wide range, it is frequently necessary to examine the other
features of the spectrum. One of these features which frequently exhibits
useful chemical information even when the photoelectron spectra do not,
624 Semiconductor Materials
PI 4fdoublet
00 78 76 74 72 i
BINDING ENERGY,eV
Figure 14: XPS spectrum of the Pt4f,/, and Pt4fS/2 doublet and Si2p transition
of the two silicide phases, PtzSi and PtSi, of platinum.
Ag 374.0 1128.2
Ag,SO, 374.2 1132.2
Sn 492.9 1048.3
Sn02 494.8 1053.4
Primary
Ion Mass
Condenser Lens
Sample
SRJJ-TERING WELOS _
(atoms/ion) FOR ARGON -
Figure 17: Plot of the sputtering yield versus ion energy for argon on copper.
628 Semiconductor Materials
Re
10’ 5
Relative i -Th
iFnsity 10’ z
*
A&
@e I ’
Cd
Sb
ITI
I I I I I Auk I I I
30 40 so 60 70 80 90 100
1
100, . . . . . . . . . . , ,
High Purlly Silicon
Oxy5rn Bombardment
IO
,, I
t
5
s 0.1
.
.?
;; 0.01
z
cc
0.001
0.0001
Figure 19: Plot of the relative secondary ion intensity versus the mass-to-charge
ratio resulting from oxygen bombardment of high purity silicon.
6OkeV “P Implants
,..‘,‘.,,,,,,,
Normalized lo Doses Determined
By Neutron Acllvalion
Dose - I 09 x I O”at/cm
1’ 13x 10’sat/cm2
1'fOx10"at/cm7
: 5: 90~10'~at/crn~
Figure 20: Depth profiles of phosphorus implanted into silicon at 80 keV show-
ing the concentration versus depth.
630 Semiconductor Materials
Even with these limitations SIMS is a powerful tool for the analysis of
dopants in semiconductor materials. Figure 21 illustrates an example
where SIMS has been used successfully. Boron profiles in Si are shown
before and after laserannealing. 22Thedopantisredistributedasaresultof
laser annealing. The maximum B concentration decreases while the tail of
the profile diffuses to a greater depth. It is interesting to note that the
redistribution of B is concentration dependent. The tail of the implant
profile is redistributed to a much greater depth after laser annealing forthe
high dose implant than the low dose implant.
SIMS is one of a few analytical tools capable of distinguishing isotopes.
This has resulted in some well designed experiments that take advantage
of this feature. For example, Coleman et al. 23 utilized SIMS to investigate
“FYS
SECOND
RELATIVE
DEPTH
Figure 22: Depth profile of the isotopic distribution of oxygen in an anodic
oxide’grown on GaAs(001): curve (A) is for GaAs(001) anodized in HzO” then
Hz016 while curve (B) is for anodization first in Hz016 then Hz018.
632 Semiconductor Materials
75&t -_
‘;i
.E 181Tat+
II
= 6000 x25
(b) SALI
_t r- I I 1
Figure 23: The positive secondary ion SIMS spectrum in trace (a) is from a
GaAs wafer under static Ar+ bombardment. A SALI spectrum is shown in trace
(b) for the same surface under static Ar+ bombardment, plus an accompanying
248 nm laser for ionization.
Characterization of Semiconductor Materials 633
the sample by ion sputtering. A laser beam is positioned above the sample
so that it intersects with the vapor cloud of atoms as shown in Figure 24.
Thelaseristunedtothefrequencynecessarytoionizetheatomofinterest.
The ions are extracted from the sample region, then energy analyzed
before passing into a mass spectrometer. The spectrum like SIMS is a plot
of the intensity versus the mass-to-charge ratio.25
RIS reportedly is capable of achieving sensitivities reaching 1 part in
10’2 and aselectivity that eliminates ambiguity in the interpretation of the
results.26 Several factors contribute to the enhanced sensitivity of RIS
over similar techniques like SIMS. RIS has a lower background. The ions
that are generated during the sputtering process are extracted before the
laser pulse ionizes the remaining neutral species. As a result, only ions
generatedinthelaserpulsepassintothemassspectrometer.Additionally,
only selected ions are generated in the laser pulse. Figure 25 shows the
five basicschemes used in theresonance ionization process?‘j Typically, a
tunable dye laser is adjusted so that it emits precisely the correct wave-
length to excite an electron in an atom from its original state to a higher
state. Occasionally, a second photon from a second laser is used to excite
the atom further to an even higher state. After excitation, a second or third
photon is used to interact with the excited atom, causing the electron to be
released from the atom. Thereby, producing a positive ion and a free
electron. The key to RIS is choosing resonately excited states that can be
easily excited and that have large photoionization cross-sections, so that
they can be ionized with high efficiency. It is possible to saturate all of
these processes with commercially available lasers, so that an atom in the
initial state will be excited through the resonant intermediate states and
into the ionization continuum with unit probability during a single laser
pulse.
RIS is extremely selective, in that, only atoms of a given element are
ionized. The intermediate excited states through which the ionization
proceeds may be chosen such that they are uniquely characteristic of that
DEFLECTION
DOUBLE-FOCUSING
MASS SPECTROMETER
ENERGY MASS
ION SOURCE
1 2 3 4 5
Figure 25: Schematic diagram of the five basic schemes for RIS.
Figure 26: Correlation plot showing the boron concentration in silicon measured
by RIS versus the value determined by electrical measurements.
2
E, = W,, M,, @E,
1
where
M, cos0 + dMp2 - M,2sin2r3
K= (3)
MI + M2
Al REST
EO
AntR COLLISION
M2’“1
k=
I
M,cosl3 + ~4~4;-
4 + Mz
’
M:sin*O
I
From Acceleralor
u u Beam
L L
E, Et kE,
ENERGY OF BACKSCATTERED ION
Figure 30: A plot of the Rutherford backscattering yield versus the energy of
the backscattered ion with an accompanying illustration showing the scattering
location in the sampled depth.
As Deposlted I
I%, - NI
I -8
6- I
3OO’C’1 90’ Annealing i I
- N12SI
300 lC 90 mln
2 10 r ]
.O :
.O ;
-0 1
:o : i
t l * I
0 L-J
0.6 1.0 1.4 1.6
ENERGY (MeV)
Figure 31: RBS spectra of the phases of nickel silicide formed following the
deposition and annealing of nickel on silicon.
Next Page
1500 - 1
I I I I I I
RANDOM
I * n
<HO> ALIGNED:
VIRGIN
AS IMPLANTED a”,: MO
ANNEALED AT 550 ‘C FOR 50 min
iz ANNEALED AT 550 ‘C FOR 50 min, 850 ‘C FOR 10 min 3
F
600
0
0.4 8 0.40 0.32 0.24 0.16 0.06 0 0.08
DEPTH (prnl
Figure 32: FIBS spectra for Si(100) in the random and <I IO> aligned directions
before and after 80 keV 3oSi+ implant and subsequent anneal at 550” and 850°C.
INDEX
Acceptors 27
Adsorption
adatom 412 421
binding energy 411
Amphoteric 362
Atomic mobility 96 102
Auger electron spectroscopy (AES) 6 611
analysis 67 70 178 294
chemical shift 614
cylindrical mirror analyzer 613
detection limit 616
energy level diagram 611
profiling 66 614
Carrier
concentration 27
velocity 27
Channel 1 4
current 4
Charge 58 59
fixed oxide 58
mobile ionic 58
oxide trapped 59
Chemical vapor deposition (CVD) 6 80
boron nitride 115 122
polycrystalline silicon 106
Si 82 94
Silicon carbide 106 115 122
SiO2 75 106 115
Si3N4 75 105 106 115
Chemisorption 333 421
Child-Langmuir Law 209 397
CMOS 89
Collector 111
Condensation 410
Conductivity 308
Crucible 12 19 20
Crystal Growth
Czochralski 9
float-zone 9 20 475
seed crystal 12 82
Crystallographic orientation 84 133
Crystal puller 11
Diffusion (Cont.)
coefficient 69 463 495
continuum theory 456
diffusivity 469 476 482 487 495
501 596
dopant 67 69
Fick’s Laws 456 495
flux 466
grain boundary 596 606 616
Henry’s Law 460
interstitialcy 466 476 500
oxidation enhanced 70 479
oxygen 40 50 67 74
predeposition 459
surface 415
vacancy 465 491 496 501
Diode 96
Dislocation 6 13 23 24 32
475
misfit 94 178
Donor 27
oxygen 40
Dopant 15 18 19 20 21
63
MBE films 361
N-type 15 22 61
P-type 15 22 61
segregation 61
substitutional 477 488
Drain current 4
This page has been reformatted by Knovel to provide easier
navigation.
Index Terms Links
Etching (Cont.)
defect etch 91 92 128 172
planar 126
RCA clean 187
reactive ion 6
Schimmel 173
Secco 172 173
Sirtl 172 173
vapor 86 89 90 91
Wright-Jenkins 173 178
Yang 173
Evaporation 329 336
flash evaporation 348
hot-wall evaporation 349
rate 337
rate monitors 354
reactive evaporation 352
sources 340
sublimation 350
three-temperature 419
uniformity 345
of alloys 347
of compounds 347
G
This page has been reformatted by Knovel to provide easier
navigation.
Index Terms Links
Gate
delay 1
insulator 4
length 1
oxide 1 60 592
voltage 4
Gettering 34 41
backside 91 94 173
oxygen precipitates 173
vacuum 335
Gibbs free energy 468
Grain boundary, polycrystalline
silicon 73 91
Implantation 6
Impurity 6 15 18 43 90
gettering 65
N-type 84
P-type 84
substitutional 477
Inelastic mean free path 613
This page has been reformatted by Knovel to provide easier
navigation.
Index Terms Links
Infrared 36 38 654
epitaxial layer thickness measurement 658
Fourier transform (FTIR) 7
sensitivity 657
Ingot 11 12 15 18 20
22
Integrated circuits 9
Interface
SiOs/Si 58 59 61 62 65
66 67 68 70 71
traps 59
Interferometry 249
Interstitial 20
oxygen 37 38 41
self 455 472 499
silicon 69 70
Ion bombardment
cleaning 422
influence on microstructure 425
influence on stress 430
ion beam mixing 428
Ion implantation 512
critical energy 520
electronic stopping 519
implantation masking 522
ion channeling 527
Lindhard, Scharff & Schiott
(LSS) range theory 528
nuclear stopping 519
Junction 82
breakdown voltage 61
depth 1 67 506
Lifetime, minority-carrier 34
Lithography 541
Abbe principle 557
contact printing 551
electron beam lithography 6 543 570
ion beam 6
numerical aperture 555 558
optical 6 559
projection printing 555
wafer stepper 556
x-ray 6 571
M
This page has been reformatted by Knovel to provide easier
navigation.
Index Terms Links
Melting point 82
Metailization 575
multilevel 440
ohmic contact 438 598
polysilicon 592
Rent’s rule 581
silicides 439 601
step coverage 440
wiring structure 581
Microstructure 6 415
Mobility 27
atomic 90
Molecular beam epitaxy 330 355
apparatus 356
application 362
epitaxy growth 360
three-temperature evaporation 420
MOS device 58 65 73 87 89
91 109 133
gate 601 603
Multiquantum well laser (MQW) 363
Oxygen 15 19 20 31 32
35 37 91
diffusion 39 40
dissolved 9 38
Plasma 191
anisotropy 191 193 227 252
chemical effects 213
DC plasma 193 195 196
feedstock composition 210
fundamental aspects 192
mechanism 230
modeling 241
non-equilibrium 191
plasma body 193
plasma sheath 193 195 208 209 393
395
reactor loading 223
RF plasma 193 194 195 196 394
surface chemistry 217
Plasma etching 242
Al 263
compound semiconductors 264
defects 257
endpoint detection 249
equipment 244
oxide 260
patterning 242
radiation damage 257
selectivity 254
Si 259
silicides 259
silicon nitride 264
throughput 256
uniformity 253
This page has been reformatted by Knovel to provide easier
navigation.
Index Terms Links
Q
Quartz crucible 9 11 12 48
reactor 82
S
Scanning electron microscopy 96 158 162 640
backscattered electrons 641
energy dispersive x-ray spectrometer 642
secondary electrons 641
wavelength dispersive analyzer 645
Scanning transmission electron
microscopy (STEM) 7 475 646
bright field image 647
darkfield image 647
electron diffraction 648
electron energy loss spectroscopy 649
This page has been reformatted by Knovel to provide easier
navigation.
Index Terms Links
Silicon monoxide 12 19
Silicon nitride 55 105 111 122
Solubility, oxygen 39
Spreading resistance 18
Sputter assisted laser ionization
(SALI) 632
Sputtering 329 364
alloys 376
bias 424 442
cleaning 422
compounds 376
glow discharge 378
ion beam 392
magnetron 382
mechanisms 365
planar diodes 378
reactive 400
RF sputtering 394
sputtered species 369
triode 381
yield 371
Stacking faults 6 32 33 69 70
91 104 155 157 173
178 183 500 647
oxidation-induced 499
Stereographic projection 133
Sticking coefficient 421
Target poisoning
reactive sputtering 402
Thermal decomposition 86
Threshold voltage 4 61
Torr 330
Transistor 82
Trap, defect 91
Wafer 22 23 25 34
flatness 25
Work damage 22 23
Work function 60
X-ray diffraction 70
X-ray lithography 115
X-ray photoelectron spectroscopy (XPS) 6 70 303 618
Auger transition 624
chemical shift 621
detection limit 624
energy analyzers 618
monochrometer 618
profiling 620
X-ray orientation 22
X-ray topography (XRT) 7 653