(Norbert H. Nickel, Evgenii Terukov) Zinc Oxide
(Norbert H. Nickel, Evgenii Terukov) Zinc Oxide
(Norbert H. Nickel, Evgenii Terukov) Zinc Oxide
Sub-Series
I.
II.
III.
IV.
IOS Press
Springer (formerly Kluwer Academic Publishers)
IOS Press
Springer (formerly Kluwer Academic Publishers)
The NATO Science Series continues the series of books published formerly as the NATO ASI Series.
The NATO Science Programme offers support for collaboration in civil science between scientists of
countries of the Euro-Atlantic Partnership Council. The types of scientific meeting generally supported
are Advanced Study Institutes and Advanced Research Workshops, and the NATO Science Series
collects together the results of these meetings. The meetings are co-organized by scientists from
,
NATO countries and scientists from NATO s Partner countries countries of the CIS and Central and
Eastern Europe.
Advanced Study Institutes are high-level tutorial courses offering in-depth study of latest advances
in a field.
Advanced Research Workshops are expert meetings aimed at critical assessment of a field, and
identification of directions for future action.
As a consequence of the restructuring of the NATO Science Programme in 1999, the NATO Science
Series was re-organized to the four sub-series noted above. Please consult the following web sites for
information on previous volumes published in the Series.
http://www.nato.int/science
http://www.springeronline.com
http://www.iospress.nl
Norbert H. Nickel
Hahn-Meitner-Institut Berlin,
Germany
and
Evgenii Terukov
A.F. Ioffe Physico-Technical Institute,
St. Petersburg, Russia
ISBN-10
ISBN-13
ISBN-10
ISBN-10
ISBN-13
ISBN-13
1-4020-3474-1 (PB)
978-1-4020-3474-9 (PB)
1-4020-3473-3 (HB)
1-4020-3475-X (e-book)
978-1-4020-3473-2 (HB)
978-1-4020-3475-6 (e-book)
Published by Springer,
P.O. Box 17, 3300 AA Dordrecht, The Netherlands.
www.springeronline.com
Contents
Contributing Authors
ix
Preface
xv
15
25
35
37
vi
47
59
69
83
99
113
115
125
133
157
159
vii
183
195
197
211
217
225
Index
239
Contributing Authors
A. Kh. Abduev
Russian Academy of Sciences, Makhachkala, Russia
A. K. Akhmedov
Russian Academy of Sciences, Makhachkala, Russia
Ya. I. Alivov
Institute of Microel. Technology, Chernogolovka, Russia
H. R. Alves
Justus-Liebig-Universitt Giessen, Giessen, Germany
A. Sh.Asvarov
Russian Academy of Sciences, Makhachkala, Russia
B. M. Ataev
Daghestan Scientific Centre of RAS, Makhachkala, Russia
P. Barquinha
New University of Lisbon, Caparica, Portugal
V. A. Baturin
National Academy of Sciences of Ukraine, Kyiv, Ukraine
I. V. Blonsky
National Academy of Sciences of Ukraine, Kyiv, Ukraine
ix
xi
xii
xiii
xiv
Preface
xvi
the fact that most growth techniques either use hydrogen as a transport gas or
H containing precursor gases such as NH3. Therefore, in order to achieve ptype conductivity ZnO has to be subjected to a post-growth thermal treatment
at temperatures above 650 C. The purpose is two-fold: (i) interstitial
hydrogen donors are eliminated and (ii) H passivated nitrogen acceptors are
activated. The properties of H in ZnO expand well beyond the formation of
donors and the neutralization of impurities and deep defects. The presence of
hidden H has been suggested, which might be molecular H2 that can be
activated to form IR active O-H centers. Moreover, in sputtered ZnO the
total H concentration can be as high as 1021 cm-3. In addition, with increasing
H concentration the average H binding energy increases; H becomes more
strongly bound. Besides N, phosphorous and arsenic have been successfully
used to fabricate p-type ZnO. Since the fabrication of p-type ZnO becomes
more and more reliable the fabrication of p-n homojunction LEDs is not far
off.
Recently, a new and rapidly developing area emerged, namely the field of
transparent electronics based on ZnO transparent thin-film transistors. These
devices exhibit an on-off ratio of 6 orders of magnitude, a channel mobility
of 1 3 cm2/Vs and can be fabricated at room temperature. Hence, already in
the near future invisible large area electronics based on ZnO will
revolutionize the electronics market. Although first small steps towards ZnO
based optoelectronic applications are currently being made there is still a vast
lack of fundamental knowledge regarding doping mechanisms, influence of
impurities, formation of point and line defects, and their interactions with
hydrogen. Some controversies and questions that arose during the workshop
are based upon the fact that a materials basis for state-of-the-arte ZnO has
not been established yet. Depending on the fabrication process the materials
qualities such as impurities, their concentration, H content, number of point
defects and grain boundaries can vary significantly.
This book reviews the latest experimental and theoretical studies in the
field of ZnO. The chapters reflect the topics that were discussed during the
NATO advanced research workshop, which brought together researchers
from NATO countries, Russia, Ukraine, and Kazakhstan to encourage cross
fertilization and wider dissemination of the advances in this important
research field.
We would like to thank NATO for funding the Advanced Research
Workshop. Above all, we would like to thank all participants for their
interest in this workshop and active work in the vibrant field of ZnO.
Norbert H. Nickel
Evgenii Terukov
Chapter 1
THE SCOPE OF ZINC OXIDE BULK GROWTH
Robert Triboulet1, Vicente Munoz-Sanjos2, Ramon Tena-Zaera2, Mari
Carmen Martinez-Tomas2 and Sad Hassani1
1
Abstract:
The techniques classically used for the growth of bulk ZnO crystals, meltgrowth, vapour growth, solution-growth and hydrothermal growth, are
reviewed. In the case of vapor growth, numerical simulations suggest that
such species as H2 + H2O, CH4, Zn, C or Fe could act as sublimation
activators. The electrical and structural properties of crystals obtained by these
various techniques are compared.
Key words:
ZnO synthesis, melt-growth, vapor-growth, solution-growth, hydrothermalgrowth, electrical and structural properties
1.
INTRODUCTION
R. Triboulet et al.
That is why research related to the ZnO bulk growth has received a
considerable interest during these last years.
As a result of the small ionic radii of its constituting elements, mainly the
oxygen one, and of its subsequently very short and energetic chemical bond,
ZnO has a very high melting point of about 1900C. Furthermore, because
of the high electronegativity of its oxygen and zinc and because of its usual
off-stoichiometry, ZnO shows an extremely high reactivity with any
surrounding material at high temperature. This intrinsic property will be
shown to affect strongly its growth, not only its melt-growth but also its
vapour-growth.
In this chapter, the various techniques classically used for the growth
of bulk ZnO crystals will be reviewed.
2.
3.
4.
R. Triboulet et al.
4.1
CSVT experiments
the one found for ZnSe. Furthermore, the variation of the instantaneous
growth rate as a function of time has been found to be non-linear, as shown
in Figure 2.
14
Ts = 600 C
C
12
10
6
Ts = 550 C
1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
-1
1000//4 (K )
Figure 1. ZnO instantaneous growth rate as the function of the inverse source temperature for
two source temperatures, 550C and 600C.
R. Triboulet et al.
From this variation of the growth rate with time and from the small H
value, compared to the ZnSe one, we suggested that the vapor growth of ZnO
with hydrogen occurred according to a complex chemically assisted
sublimation mechanism involving a competition between two different
mechanisms, one associated with ZnO, the other with Zn.
4.2
Numerical simulation
1 .6
F(ZnO)-sub /F(ZnO)-crys
1 .4
1 .2
1 .0
0 .8
0 .6
1 e -5
1 e -4
1 e -3
1 e -2
1 e -1
T o ta l p re s s u re (a tm )
Figure 3. Calculated ratio between ZnO mass fluxes at both interfaces as a function of total
pressure under equilibrium conditions.
10
5.
R. Triboulet et al.
Such solvents as PbF227 and V2O5/P2O5 mixtures28 have been used for the
flux growth of ZnO crystals. Using PbF2 as the solvent in sealed Pt
crucibles, ZnO crystals have been grown by THM.29 PbCl2 has been found to
be a very good solvent of ZnO as well, but showing also great reactivity with
any surrounding material but platinum.30 The same authors have found the
Zn-In alloys to be good solvents of ZnO without reactivity with silica. NaCl
has been found recently to be an excellent solvent of ZnO allowing solutiongrowth in silica at temperatures of about 900C.31 Figure 4 pictures a ZnO
crystal grown using NaCl as a solvent.
6.
11
hydrothermal method has been shown suitable for the growth of large ZnO
crystals from (OH)- solutions at temperatures <500C under high pressure
(15 to 150 MPa) with a temperature difference 'T
T ~ 3 - 80K and a growth
rate in the range 0.05-0.3 mm/day.32-37 More recently, high quality ZnO
crystals, as demonstrated by X-ray rocking curves in the 40 arcsec range and
sharp PL peaks, have been grown hydrothermally at 355 C with a 'T of
10C from NaOH/KOH solutions as the solvent, pressures of 18,000 to
22,000 psi (maximum pressure 1500 atm).38 The crystals present a
dislocation density < 500 cm-2, a X-ray double diffraction FWHM ~130
arcsec and carrier mobilities ~175 cm2/Vs @ RT. The O and Zn surfaces of
(0001) planes have been found to behave differently. GaN layers deposited at
750 C on such substrates show DDX FWHM of 735 arcsec. Using a
KOH/LiOH solvent, ZnO bulk single crystals about 10 mm in dimension
have been grown at temperatures less than 400 C at pressures ranging from
830 to 1110 Mpa.39
Hydrothermal growth of ZnO remains very active in the Institute of
Crystallography of the Russian Academy of Sciences.40 Typical growth
conditions are40: temperature of the dissolution zone 270-350C, temperature
difference between the growth and the dissolution zones, 'T ~ 20-80C,
pressure 15-140 MPa, solvent KOH + LiOH, silver liner of capacity 500 cm3
put in autoclaves of capacity 1000 cm3. The growth rate is shown to depend
on the crystal orientation on the basis of the polar character of the crystal
structure along the c-axis and the presence of negatively charged species in
hydrothermal solutions.40 Typical crystals have a cross section equal to about
3 cm and a height along the hexagonal axis equal to 1-1.5 cm.41 They contain
concentrations 10-2-10-4 wt. % of Li and also concentrations 10-3-10-4 wt. %
of Cu, Mg, Si, Fe, Mn and Ag.
Very spectacular results in the hydrothermal growth of ZnO have been
recently reported by Oshshima et al.42 from Fukudas laboratory in Senda.
2-in ZnO single crystals (size up to 505015 mm3) have been grown using
KOH/LiOH @ T
T= 300-400C and P = 80-100 MPa. Rocking-curve FWHM
of 8 arc.sec for (002) reflections and 25 arc.sec for the (101) reflections are
reported (to be compared with 42 and 370 arc.sec respectively for a CVTgrown one). But incorporation of Li and K originated from LiOH and KOH
was observed at concentrations of 0.5-12 ppm for Li and 0.04-0.2 ppm for K.
In addition, incorporation of Fe and Al was also observed. Note that the
complex equipment used includes a 3 m length and 200 mm inner diameter
Pt liner preventing impurity incorporation from the autoclave inner surface.
12
7.
R. Triboulet et al.
KOH/LiOH
2-in
8
Pmax ~ 80-100 MPa
Growth rate ~ 8 m/h
Source 1000-1200C
SCVT
dist. source-substrate 3 inch
40
81016
res. atm. H2+H2O; seeds
40 mm
Growth rate: 40 m/h
CVT
C as an activator
cm-size
28
cm2/Vs
130
Ref.
150
3
3335
175
38
13
20
Whatever the method used, the electrical and structural properties of the
crystals are very close except the record structural quality recently reported
for hydrothermally grown crystals by Ohshima et al.38 Note the low growth
rates of the vapor grown crystals and more again of the hydrothermally
grown ones which could be a handicap in view of industrial production.
8.
CONCLUSIONS
In spite of the numerous efforts dedicated to the ZnO bulk growth for
several years, high quality large and cheap ZnO substrates are not yet
available. The best crystals obtained so far as far as the crystallographic
quality is concerned have been grown by the hydrothermal method, whose
possibilities are limited by an extremely weak growth rate and the
unavoidable incorporation of impurities coming from the solvent, in addition
to complex equipment. The better knowledge and understanding of the
mechanisms of vapor transport can let expect significant progresses in this
field in a near future.
13
REFERENCES
1. R. Triboulet, Proceedings SPIE 4412, 1 (2000).
2. G. Agarwal, J. E. Nause, and D. N. Hill, Mat. Res. Soc. Symp. Proc. 512, 41 (1998).
3. D. C. Reynolds, C. W. Litton, D. C. Look, J. E. Hoeschler, B. Claflin, T. C. Collins, J.
Nause, and B. Nemeth, J. Appl. Phys. 95, 4802 (2004).
4. M. Shiloh and J. Gutman, J. Crystal Growth 11, 105 (1971).
5. W. Piekarczyk, S. Gazda, and T. Niemyski, J. Crystal Growth 12, 272 (1972).
6. K. Matsumoto, K. Konemura, and G. Shimaoka, J. Crystal Growth 71, 99 (1985).
7. K. Matsumoto and G. Shimaoka, J. Crystal Growth 86, 410 (1988).
8. K. Matsumoto and K. Noda, J. Crystal Growth 102, 137 (1990).
9. M. Hirose and I. Kubo, Jpn. J. Appl. Phys. 8, 402 (1969).
10. Y. S. Park and D. C. Reynolds, J. Appl. Phys. 38, 756 (1967).
11. M. Hirose, Y. Furuya, and I. Kubo, Jpn. J. Appl. Phys. 9, 726 (1970).
12. E. Scharowsky, Z. Phys. 135, 318 (1953).
13. K. F. Nielsen, J. Crystal Growth 3/4, 141 (1968).
14. K. J. Fischer, J. Crystal Growth 34, 139 (1976).
15. I. Kubo, J. Phys. Soc. Jap. 16, 2358 (1961).
16. T. Takahasi, A. Ebina, and A. Kamigawa, Jpn. J. Appl. Phys. 5, 560 (1966).
17. E. A. Weaver, J. Crystal Growth 1, 320 (1967).
18. M. Hirose and Y. Furuya, Jpn. J. Appl. Phys. 9, 423 (1970).
19. D. C. Look, D. C. Reynolds, J. R. Sizelove, R. L. Jones, C. W. Litton, G. Cantwell, and W.
C. Harsch, Solid State Comm. 105, 399 (1998).
20. J.-M. Ntep, S. Said Hassani, A. Lusson, A. Tromson-Carli, D. Ballutaud, G. Didier, and R.
Triboulet, J. Crystal Growth, 207, 30 (1999).
21. S. Hassani, A. Tromson-Carli, A. Lusson, G. Didier, and R. Triboulet, Phys. Stat. Sol. 229,
835 (2002).
22. A. Mycielski, A. Szadkowski, L. Kowalczyk, B. Witkowska, W. Kaliszek, B. Chwalisz, A.
Wysmolek, R. Stepniewski, J. M. Baranowski, M. Potemski, A. Witowski, R. Jakiela, A.
Barsz, P. Aleshkevych, M. Jouanne, W. Szuszkiewicz, A. Suchocki, E. Lusakowska, E.
Kaminska, and W. Dobrowsolski, Proc. 11th International Conference on II-VI
Compounds, to be published.
23. J.-M. Ntep, M. Barb, G. Cohen-Solal, F. Bailly, A. Lusson, and R. Triboulet, J. Crystal
Growth 184/185, 1026 (1998).
24. R. Triboulet, J.-M. Ntep, M. Barb, P. Lemasson, I. Mora-Sero, and V. Munoz, J. Crystal
Growth 198/199, 968 (1999).
25. R.Tena-Zaera, M. C. Martnez-Toms, S. Hassani, R. Triboulet, and V. Muoz-Sanjos, J.
Crystal Growth, to be published.
26. N. Ramachandran, C. Su, and S. L. Lehoczky, J. Crystal Growth 208, 269 (2000).
27. K. F. Nielsen and E. F. Dearborn, J. Phys. Chem. 64, 1762, (1960).
28. B. M. Wanklyn, J. Crystal Growth 7, 107 (1970).
29. G. A. Wolf and H. E. LaBelle Jr., J. Am. Ceram. Soc., 48, 441 (1965).
30. J.-M. Ntep, PhD thesis, Paris 2000.
31. S. Hassani, G. Didier, P. Galtier, and R. Triboulet, to be published.
32. R. A. Laudise and A. Ballman, J. Phys. Chem. 64, 688 (1960).
33. E. D Kolb and R. A. Laudise, J. Am. Ceram. Soc.49, 302 (1966).
34. R. A. Laudise, E. D. Kolb, and A.J. Caporaso, J. Am. Ceram. Soc. 47, 9 (1964).
35. E. F. Venger, A. V. Melnichuk, L. Yu. Melnichuk, and Yu. A. Pasechnik, phys. stat. sol.
(b) 188, 823 (1995).
14
R. Triboulet et al.
Chapter 2
GROWTH MECHANISM OF ZnO LAYERS
A. Kh. Abduev1, A. Sh.Asvarov1, A. K. Akhmedov1, I. K. Kamilov1 and S.
N. Sulyanov2
1
Abstract:
Key words:
1.
INTRODUCTION
15
N.H. Nickel and E. Terukov (eds.),
Zinc Oxide A Material for Micro- and Optoelectronic Applications, 1524.
2005 Springer. Printed in the Netherlands.
16
2.
EXPERIMENT
3.
RESULTS
3.1
17
Figure 1. ZnO tetrahedral aggregate. The inset shows the model of Zn12O nucleus-cluster
18
3.2
19
Figure 2. SEM micrograph of cross-sections of ZnO layers prepared with various target-tosubstrate distances: a - 30 mm; b - 60 mm; c - 100 mm.
20
Figure 3. SEM micrographs of surfaces (a1, b1) and cross-sections (a2, b2) of ZnO layers
deposited by dc-magnetron co-sputtering (a - Tsub=3000C; b - Tsub=5000C).
21
Figure 4. XRD date of: curve 1 - without co-sputtering (only sputtering of ZnO:Ga target);
curve 2 - instrumental resolution function; curve 3 - co-sputtering (Tsubstrate=3000C); curve 4 co-sputtering (Tsubstrate=5000C)
22
Figure 5. SEM micrographs of the cross-section of ZnO layers, deposited at various discharge
currents for the Zn target: a - 0 mA; b - 25 mA; c - 50 mA; d - 90 mA.
4.
23
CONCLUSION
ACKNOWLEDGMENTS
This research was supported by Analytical Centre of DSC RAN.
REFERENCES
1.
2.
3.
4.
Sunglae Cho, Yunki Kim, Yi Sun, George K.L., Appl. Phys. Lett. 75, 2761 (1999).
M. H. Hiang, S. Mao, H. Feick et al, Science 292, 1897 (2001).
K. Minegishi, Y. Koiwai, Y. Kikuchi et al. , Jpn. J. Appl. Phys. 36, L1454 (1997).
Y. R. Ryu, W. J. Kim, H. W. White. , J. of Cryst. Growth, 219, 419 (2000).
24
5. P. F. Carsia, R. S. McLean, M. H. Reilly et al. , Appl. Phys. Lett. 82, 1117 (2003).
6. T. Minami, MRS BULLETIN AUGUST, 38-44 (2000).
7. S. Roy, S. Basu, Bull. Mater. Sci. 25, 513 (2002).
8. C. R. Gorla, N. W. Emanetoglu, S. Liang et al., J. Appl. Phys. 85, 2595 (1999).
9. B. S. Li, Y. C. Liu, Z. Z. Zhi et al., J. Vac. Sci. Technol. A 20, 1779 (2002).
10. H. Z. Wu, K. M. He, D. J. Qiu, and D. M. Huang, J. of Crystal Growth 217, 131 (2000).
11. T. Ohgaki, N. Ohashi, H. Kakemoto et al., J. Appl. Phys. 93, 1961 (2003).
12. M. Chen, Z. L. Pei, C. Sun et al., J. Crystal Growth 220, 254 (2000).
13. K. Tominaga, T. Murayama, I. Mori et al., Thin Solid Films 386, 267 (2001).
14. Y. Igasaki, T. Naito, K. Murakami et al., Appl. Surf. Science 169-170, 512 (2001).
15. F-J. Haug, Zs. Geller, H. Zogg et al, J. Vac. Sci. Technol. A 19, 171 (2001).
16. T. Minami, T. Miyata, T. Yamamoto et al., J. Vac. Sci. Technol. A 18, 1584 (2000).
17. B. M. Smirnov, Phys. Usp. 37, 621 (1994).
18. V. N. Tsytovich, Phys. Usp., 40, 53 (1997).
19. A. Kh. Abduev, A. K. Achmedov, V. G. Baryshnikov et al., Technical Physics Letters 26,
332 (2000).
20. A. Kh. Abduev, A. Sh. Asvarov, A. K. Achmedov et al., Technical Physics Letters, 28, 952
(2002).
21. G. Knuyt, C. D. Quaeyhaegens, J. Haen et al., Thin Solid Films 258, 159 (1995).
22. I. Sieber, N. Wanderka, I. Urban et al., Thin Solid Films 330, 108 (1998).
23. S. N. Sulyanov, A. N. Popov, and D. M. Kheiker, J.Appl.Cryst. 27, 934, (1994)
Chapter 3
KINETICS OF HIGH-TEMPERATURE DEFECT
FORMATION IN ZnO IN THE STREAM OF
OXYGEN RADICALS
M. B. Kotlyarevsky1, I. V. Rogozin2, and A. V. Marakhovsky2
1
Abstract:
Key words:
1.
INTRODUCTION
25
N.H. Nickel and E. Terukov (eds.),
Zinc Oxide A Material for Micro- and Optoelectronic Applications, 2534.
2005 Springer. Printed in the Netherlands.
26
M. B. Kotlyarevsky et al.
27
2.
KINETIC ANALYSIS
This method consists of increasing the atomic chalcogen flux JB to the above
equilibrium value for a given pressure and temperature via dissociation
of chalcogen molecules in response to external influences. As shown,10 the
fraction of dissociated oxygen molecules in an RF discharge attains 15%. At
pO2 1 Pa, the dissociation of 1% of the oxygen molecules produces atomic
oxygen pressure pO ~10-2 Pa. The same activated oxygen pressure can be
achieved via thermal dissociation of molecular oxygen at 1000 K and a
pressure pO2 pO2 K D 102 GPa where KD is the equilibrium constant of the
reaction O2 = 2O + D (D
( is the dissociation energy of molecular oxygen).
Thus, RF activation during heat treatment in the vapor phase produces an
atomic oxygen pressure that is essentially unattainable via thermal
dissociation under typical heat-treatment conditions.
28
M. B. Kotlyarevsky et al.
subsurface layer
Figure 2. (a) Schematic illustrating the ADC equilibrium on the surface of an AB crystal and
(b) the corresponding electronic transitions: (1) transition resulting from the interaction of a
surface B atom with a VA vacancy ( J ddes ), (2) transition resulting from the interaction of
a surface B atom with a hole ( J ddes ), (3) transition resulting from the recombination of a
).
B adatom with a VBx vacancy ( J
rec
29
the crystal surface, we took into account the interactions of neutral and
charged defects with the vapor phase.
System of rate equation describing processes of defect formation can be
written in the form:
dN B'
dW
JB
J dec
2J des 2
J cr
dVA
dW
J cr
dVA'
dW
dVB
dW
J ddes
dp
dW
J rec
g [VA ] h p[
p VB VA' .
J rec
'
A
J rec
] J ddes
J rec
(1)
(2)
(3)
(4)
J rec
(5)
(6)
30
M. B. Kotlyarevsky et al.
TC
GdesB
k
Eg
Ea1
Ed1
GA
Gddiifff B
(7)
Figure 4 shows the critical temperature as a function of the atomicoxygen flux. These data demonstrate that the critical temperature can only be
raised substantially by increasing JB. The calculations show that for ZnO Tc =
m2s1 in really reachable
800 K at chalcogen fluxes of orders 1016-1017 cm
experiments. Thus, it follows from the above thermodynamic and kinetic
analyses that one way of controlling the chemistry of native defects in II-VI
compounds is by increasing JB to above the equilibrium value at a given
pB2 and temperature via activation of chalcogen molecules.
3.
-2
log ( D ef. ) ( cm )
14
13
31
'
V Zn
VO
*
VO
'
V Zn
12
h
11
V Zn
10
VO
9
8
14
VO
V Zn
15
16
17
18
19
20
-2 -1
log JS ( cm s )
Figure 3. Defect densities in the surface layer of ZnO the flux of oxygen radicals (atoms) at
1200 (solid lines) and 1000 K (broken lines).
magnetic filter. The atomic oxygen flux, JO=1016 to 1017 cm-2s-1, was
determined as described. ZnO layers grown at an oxygen radical flux JO |
1.51017 cm-2s -1 and annealing temperatures in the range of 670870 K were
p-type, as determined by Hall effect measurements. The lowest resistivity ptype ZnO layers, with U | 6102 :cm and hole mobility p ~23 cm2/(Vs),
were obtained at 710 K. Raising the annealing temperature to 810 K was
found to increase U to about 3107 :cm. The layers produced at Tann t 910 K
were n-type, with the resistivity above 109 :cm. These data indicate that the
critical annealing temperature of ZnO is about 910 K, in accordance with
theoretical estimates.
In the UV range, the PL spectrum of oxygen-enriched ZnO ((p-type)
shows a strong line at 369.5 nm with a full width at half maximum (FWHM)
of 14 meV. In the visible range, we observe a weak band centered at 400 nm.
Comparison of the luminescence spectra of n- and p-type ZnO indicates that
the 369.0 and 369.5-nm lines are due to bound excitons. As shown by
Butkhuzi et al.,9 the 369.0 and 369.5-nm emissions arise from excitons
bound to neutral donors and acceptors, respectively. Increasing the annealing
time at 710 K increases the intensity of the 400-nm band. The spectra of
samples annealed for 4 h show only the 400-nm band, which spans the entire
32
M. B. Kotlyarevsky et al.
TC (K)
1000
900
800
700
600
12
14
16
18
20
22
-2 -1
log JS (cm s )
Figure 4. Critical temperature as a function of the atomic oxygen flux.
4.
The kinetics of defect formation in a heterophase system such as a singlecrystal compound of one chalcogen, B' (e.g.: selenium), exposed to a flow of
radicals of another chalcogen, B'' (e.g.: oxygen), are more difficult to
analyze. The general trends of defect formation in such systems must,
however, be the same as above. The formation of VA vacancies in an AB'
crystal is accompanied by the transition of A atoms to the surface layer,
where they combine with B'' atoms to form AB' dimers. In the system of rate
equations (1-5) one must take into account the formation of molecules
containing both chalcogens (e.g.: SeO2) and desorption of such molecules.
This process results in the formation of VB' defects in the crystal, whose
composition gradually changes to AB'' owing to B'' in-diffusion through
VB' vacancies.
33
JO
J des
d O
dN O
dW
>
@ JO
dN ZZn
dW
J cr
dNSSe
dW
J ddes
2J des
d O2
J rec
J rec
J des
(8)
(9)
(10)
(11)
34
M. B. Kotlyarevsky et al.
REFERENCES
1. A. N. Georgobiani, M. B. Kotlyarevsky, V. N. Mikhalenko, USSR Inventors Certificate
no. 684, 810 (1977).
2. N. K. Morozova and O. M. Morozova, Izv. Akad. Nauk SSSR, Neorg. Mater. 17, 1335
(1981).
3. M. T. Nyges, Cand. Sci. (Phys.Math.) Dissertation, Tallinn (1974).
4. V. A. Nikitenko, S. A. Stenli, and N. K. Morozova, Izv. Akad. Nauk SSSR, Neorg. Mater.
24, 1830 (1988).
5. V. N. Mikhalenko, B. P. Dementev, M. B. Kotlyarevsky, and A. N. Georgobiani, Izv.
Vyssh. Uchebn. Zaved., Fiz. 8, 150 (1978).
6. A. N. Georgobiani, M. B. Kotlyarevsky, and V. N. Mikhalenko, Tr. Fiz. Inst. im. P. N.
Lebedeva, Akad. Nauk. SSSR. 138, 79 (1983).
7. A. N. Georgobiani, M. B. Kotlyarevsky, and V. N. Mikhalenko, V.N., Izv. Akad. Nauk
SSSR, Neorg. Mater. 18, 12 (1982).
8. M. B. Kotlyarevsky, A. N. Georgobiani, I. V. Rogozin, and A. V. Marakhovskii, Zh. Prikl.
Spektrosk. 70, 86 (2003).
9. T. V. Butkhuzi, A. V. Bureyev, A. N. Georgobiani, N. P. Kekelidze, and T. G. Khulordava
J. Cryst. Growth 117, 366 (1992).
10. A. N. Georgobiani, M. B. Kotlyarevsky, V. V. Kidalov, and I. V. Rogozin, Neorg. Mater.
33, 232 (1997) [Inorg. Mater. (Engl. Transl.) 33, 185].
Chapter 4
ELECTRICAL PROPERTIES OF ZnO
David C. Look1, Bruce B. Claflin1, Gene Cantwell2, Seong-Ju Park3, and
Gary M. Renlund4
1
Abstract:
Key words:
1.
INTRODUCTION
ZnO has long been touted as a potential candidate for UV emitters, but
the lack of a good, reproducible p-type technology has hindered its
development in the past.1 Fortunately, in the last few years, much progress
has been made in this area and a few p-n homojunctions and heterojunctions
have even shown light emission.2 In spite of this recent success, much better
performance will be needed if ZnO is to compete with GaN-based lightemitting diodes (LEDs) in the blue and UV regions. One of the key required
elements will be the growth of bulk and epitaxial ZnO with low background
37
N.H. Nickel and E. Terukov (eds.),
Zinc Oxide A Material for Micro- and Optoelectronic Applications, 3746.
2005 Springer. Printed in the Netherlands.
38
2.
Few semiconductor laboratories in existence today are without a Halleffect system capable of room-temperature operation, at least. The only basic
requirements are a current source, a voltmeter, and a small magnet, of
typical field strength 0.1 1.0 T. If the current is imposed perpendicular to
the magnetic field, a voltage will be generated in the mutually perpendicular
direction, and this is the classical Hall voltage. The current, magnetic field
strength, and Hall voltage are all easily measurable quantities, and their
combined values immediately give the Hall coefficient R, where R is related
to the carrier concentration, n or p, by n = r/eR. Here e is the electronic
charge and r is the so-called Hall factor, which is usually close to unity but
also may be calculated by fitting the temperature-dependent Hall mobility
data. Indeed, the Hall mobility PH is simply defined as PH = RV, where V is
the conductivity, measured from the voltage drop parallel to the current. To
relate n and PH to physical parameters, such as the donor concentration, we
must solve the Boltzmann transport equation, which describes the relaxation
of carrier momentum via carrier scattering. If the various scattering
mechanisms are elastic, then a relaxation time W can be defined, and the
Boltzmann equation can be easily solved to give
PH
e W
m* W
(1)
39
W ( )
wff 0
dE
wE
E
o
f
3/ 2 wff 0
E
dE
0
wE
E
W n ( E ) E 3/2
W n ( ) E 3/ 2 e- E / kT dE
E 3/ 2 e- E / kT dE
(2)
Here, E is the carrier energy and the second expression in Eq. 2 is valid for
Boltzmann (nondegenerate) statistics. We must note that W is the total
relaxation time, made up of contributions from various phonon and impurity
scattering mechanisms. The additive quantity is the relaxation rate, W-1, so
that the total relaxation rate is given by
W-1(E
( ) = Wac-1(E) + Wpe-1(E) + Wpo-1(E
(E ) + Wii-1(E) + Wdis-1(E)
(3)
W ac
S ! 4 U d s 2 E -1/2
2
21/2 E1
3/ 2
kT
(4)
N Di
1 n / I
i
(5)
Di
40
g 0i
e
g 1i
D Di
k
N C 'T 3 / 2 e
ED 0i
kT
(6).
2000
NA = 2 x 10
15
5
cm
-3
3
P (cm /V s)
1500
1000
500
0
100
200
300
T (K)
Figure 1. Temperature-dependent mobility of bulk ZnO.
400
10
41
15
Theo
ory: ND1=9
9 x 10
17
17
-3
n (cm )
16
10
15
-3
ND22=1 x 10
cm , ED2=61 meV;
NA= 2 x 10
cm , fro
om mo
obility
15
10
-3
cm , ED1=31 meV;
10
3
-3
15
20
-1
10 /T (K )
Figure 2. Temperature-dependent carrier concentration of bulk ZnO.
3.
We first will discuss the identities of the acceptor and two donors found
in the VP-grown ZnO represented by Figs. 1 and 2. The acceptor was found
by positron annihilation to be the Zn vacancy VZn.5 This is somewhat
surprising in light of secondary-ion mass spectroscopic (SIMS) data on
similar material which shows that N is present at the 1017-cm-3 level.6
However, we believe that the N is passivated by H that is present during the
VP growth process,7,8 and that a further supply of H is responsible for the
shallow 31-meV donor, ND1. This assertion is supported by annealing
results, shown in Figs. 3 and 4. The annealing temperature, 715 qC, is high
enough to get rid of H, and indeed, the donor of shallow slope in Fig. 4 has
disappeared. Furthermore, the acceptor concentration has increased from
about 2u1015 cm-3 to about 1u1016 cm-3, as illustrated by the strong decrease
in low-temperature mobility shown in Fig. 3. Since H-N complexes have
42
2500
unannealed
anneal: 715 qC, 30 min
1500
P (cm /V s)
2000
1000
500
0
100
200
300
T (K)
Figure 3. Effects of annealing on mobility of bulk ZnO. The fits (not shown) give acceptor
concentrations of NA = 1.8u1015 cm-3 for the unannealed sample, and 1.2u1016 cm-3 for the
annealed sample.
43
10
16
10
15
10
14
-3
n (cm )
unannealed
annealed: 715 qC, 30 min
10
20
30
3
40
50
-1
10 /T (K )
Figure 4. Effect of annealing on carrier concentration of bulk ZnO. The fits (not shown) give
donor concentrations of 2.8u1016 and 3.7u1016 cm-3, and energies of 37 and 66 meV,
respectively, for the unannealed sample; and donor concentrations of 2.1u1016 and 2.8u1016
cm-3, and energies of 44 and 77 meV, respectively, for the annealed sample.
10
8
unannealed
annealed 715 qC
3.3628 eV
3.3329
6
4
x10
3.3596
3.3484
3.3296
3.3566
3.33
3.34
3.35
3.36
E (eV)
Figure 5. Effect of annealing on 4-K photoluminescence of bulk ZnO.
44
4.
P-TYPE ZnO
20
10
5
theory
17
-3
cm )
50
Pp
2
1
100
200
300
T (K)
Figure 6. Hole concentration p and mobility Pp for an As-doped ZnO sample. The solid line is
a theoretical fit to the mobility.
Here, the mobility Pp is fitted (solid line) with a very simple, singlevalence-band model, assuming an effective mass of 0.64 m0. The only fitting
parameter is the donor concentration, which turns out to be ND = 8u1019 cm3
. Although the temperature dependence of the hole concentration p cannot
be easily fitted, because it is partially degenerate, we can estimate that p |
NA ND | 1019 cm-3 at high temperatures, giving NA | 9u1019 cm-3. Indeed,
SIMS results show an As concentration of mid-1019, consistent with As
being the dominant acceptor in this material. The fact that As can be doped
to such a high concentration is quite surprising, because, unlike the NO case,
As is much larger than O, and should not substitute for it readily. Some
theoreticians have recently suggested that As does not enter the lattice
simply as AsO, but as a much more complex acceptor-type defect, AsZn-
45
2VZn11. Although this defect is consistent with our data, still it must be
confirmed by more direct means.
The most studied acceptor dopant up to now, by far, is N, and this
element has been shown by electron paramagnetic resonance data to go
predominantly onto the O site and have acceptor nature12,13. It has an
acceptor transition at ~ 90 meV, for [N] ~ 1019 cm-3, and this energy can be
estimated to be about 130 150 meV at infinite dilution14. Thus, it is perhaps
not surprising that p-type ZnO can be produced with N doping;
however, it is still somewhat of a mystery as to why As, P, and Sb also work.
In any case, the recent success in creating p-type ZnO bodes well for the
future of the ZnO LED industry, although much more research will be
necessary before the situation is well understood.
5.
SUMMARY
REFERENCES
1. D. C. Look, Mater. Sci. and Eng. B 80, 383 (2001).
2. D .C. Look and B. Claflin, phys. stat. sol. (b) 241, 624 (2004).
3. D. C. Look, Electrical Characterization of GaAs Materials and Devices (Wiley, New
York, 1989).
4. B. R. Nag, Electron Transport in Compound Semiconductors (Springer, Berlin, 1980).
5. F. Tuomisto, K. Saarinen, and D. C. Look, Phys. Rev. Lett. 91, 205502 (2003).
6. D. C. Look, D. C. Reynolds, C. W. Litton, R. L. Jones, D. B. Eason, and G. Cantwell,
Appl. Phys. Lett. 81, 1830 (2002).
46
Chapter 5
ELECTRICAL PROPERTIES OF ZnO THIN
FILMS AND SINGLE CRYSTALS
M. Grundmann1, H. von Wenckstern1, R. Pickenhain1, S. Weinhold1, B.
Chengnui1, and O. Breitenstein2
1
Abstract:
We have investigated the electrical properties of ZnO single crystals and thin
films. Low resistive ohmic contacts were realized by annealing thermally
evaporated Ti-Au contacts at 400C in vacuum. Hall effect measurements
provided information about dominant scattering mechanisms in the single
crystals and the epitaxial thin films. Dominant donor levels were found in an
energy range from 37 meV to 350 meV below the conduction band for the
different samples. Schottky contacts of good quality ((j(-5V) = 4 x 10-5 A/cm2, n
= 1.3) were produced on thin films after etching with HNO3. Deep level
transient spectroscopy yielded the activation energy of deep traps for single
crystals of 300 meV and 540 meV and for the thin films of 300 meV,
respectively.
Key words:
1.
INTRODUCTION
The wide band gap semiconductor ZnO has attracted much attention in
recent years. This is stimulated by the development of transparent
electronics, the growing need of emitters and detectors working in the ultraviolet range of the electromagnetic spectrum, and the progress in the
realization of ferromagnetic semiconductors with Curie temperatures close
to or above room temperature. ZnO is nowadays applied as transparent front
contacts of solar cells competing with the commonly used material indium
tin oxide (ITO). The low cost of ZnO as well as the radiation hardness1 make
47
N.H. Nickel and E. Terukov (eds.),
Zinc Oxide A Material for Micro- and Optoelectronic Applications, 4757.
2005 Springer. Printed in the Netherlands.
48
M. Grundmann et al.
2.
SAMPLES
3.
49
as-deposited
400C 3min vacuum
400C 3min N2
600C 10min N2
0.006
0.004
Current (A)
0.002
0.000
-0.002
-0.004
-0.006
-0.008
-1.5
-1.0
-0.5
0.0
0.5
1.0
1.5
Voltage (V)
4.
50
M. Grundmann et al.
Figure 2. Lock-in thermography images of PLD sample grown with a) rotating substrate
holder, b) fixed substrate holder, and c) fixed substrate holder but the outer 2 mm of the
sample where removed. In (a, c) the circular contacts are in the top and bottom corners. In part
(b) the contacts are stripe-like at the top and the bottom edge.
5.
5.1
The Hall mobility H of the thin films is highest within a narrow region of
free electron concentrations n ranging from about 2 to 5 x 1016 cm-3 as
depicted in Fig. 3. In general, the Hall mobility of the thin films is lower
51
5.2
To model the temperature dependence of the Hall data we first fitted the
Hall mobility accounting the intrinsic scattering mechanisms piezoelectric
scattering and deformation potential scattering on longitudinal acoustic
phonons and polar optical scattering (for them we used material constants, no
fitting parameters are needed for the intrinsic scattering mechanisms) as
52
M. Grundmann et al.
103
H (cm/Vs)
SCVT
HVPE II
HVPE I
thin film
102
10
20
40
60 80 100
200
400
T (K)
Figure 4. Hall mobility of single crystals and a thin film.
53
1017
1016
SCVT
thin film
1012
1014
n (cm 3)
n (cm 3)
1015
1013
HVPE I
1011
1012
1011
1010
HVPE II
1010
3.00
3.25
HVPE II HVPE I
2
3.50
3.75
1000/T (1/K)
10
12
4.00
14
1000/T (1/K)
Figure 5. Free electron concentration (markers) and corresponding fits (lines) of single
crystals and a thin film. The inset shows n(T) of HVPE I and HVPE II.
54
M. Grundmann et al.
energy of effective mass donors. The most probable cause of this level is Al
diffusing during growth from the sapphire substrate in to the thin film. The
lowest compensation ratio J of about 0.02 have the single crystal grown by
SCVT and the PLD thin film. This is a very low ratio indicating the high
quality of present state of the art grown ZnO. The parameters of the Hall
analysis are summarized in table 1.
6.
Figure 6. Electron beam induced current scans on a) a sample treated with N2O plasma and b)
a sample treated with HCl prior to the metal deposition.
The dependence of the Schottky barrier height on the metal work function
is weak10 and the ideality factors of the diodes are far from ideal indicating
the presence of surface states, recombination in the space charge region, or
diffusion limited currents. Etching the surface with HNO3 leads to Schottky
contacts of good quality on thin films as we found recently. The current
density at 5V is about 4 x 10-5 A/cm and the ideality factor determined for
voltages between 0.05 V and 0.22 V is 1.3. The current density vs. voltage
for this diode is depicted in Fig. 7.
55
10-3
jd (A/cm)
10-4
10-5
10-6
10-7
10-8
-5
-4
-3
-2
-1
Ud (V)
Figure 7. Current density in dependence on the voltage for a Schottky diode realized on a ZnO
PLD thin film with Pd after etching with HNO3.
6.1
6.2
56
M. Grundmann et al.
80
70
300 K
60
175 K
50
40
50 K
30
20
4K
10
3.2
3.3
3.4
3.5
3.6
energy (eV)
and PLD thin films. The series resistance of the samples HVPE I and II is too
large for measuring DLTS with a capacitance bridge working at 1 MHz
which we have at our disposal. We applied rate windows between 64 and
10725 Hz and found two traps in the single crystal and one trap in the thin
film. Arrhenius plots of the temperature shift of the peak maximums with
changing rate window are shown in Fig. 9. The slope in this plot corresponds
0
-1
SCVT
thin film
ln(en/T2)
-2
-3
-4
-5
-6
-7
-8
3.0
3.5
4.0
4.5
5.0
5.5
6.0
1000/T (1/K)
Figure 9. Arrhenius plot for deep levels found in the single crystal grown by SCVT and the
thin film
57
to the activation energy Et. A defect with Et = 300 meV is found in both
samples and was also found in HVPE I and II by the TDH experiments. The
existence of the defect in all samples suggests that it is an intrinsic defect
likely the oxygen vacancy. A deeper trap with an activation energy of 540
meV is contained in the SCVT sample.
7.
SUMMARY
REFERENCES
1. D. C. Look, Mater. Sci. Eng., B 80, 383 (2001).
2. R. Schmidt, B. Rheinlnder, M. Schubert, D. Spemann, T. Butz, J. Lenzner,E. M.
Kaidashev, M. Lorenz, A. Rahm, H. C. Semmelhack, and M. Grundmann, Appl. Phys.
Lett. 82, 2260 (2003).
3. E. M. Wong and P. C. Searson, Appl. Phys. Lett. 74, 2393 (1999).
4. M. J. Vellekoop, C. C. G. Visser, P .M . Sarro, A. Venema, Sensors and Actuators, A21A23, 1027 (1990).
5. M. Lorenz, E. M. Kaidashev, H. v. Wenckstern, V. Riede, C. Bundesmann, D. Spemann,
G. Benndorf, H. Hochmuth, A. Rahm, H.-C. Semmelhack und M. Grundmann, Solid-State
Electronics 47, 2205 (2003).
6. C. Bundesmann, M. Schubert, D. Spemann, T. Butz, M. Lorenz, E. M. Kaidashev, N.
Ashkenov, H. Neumann, and M. Grundmann, Appl. Phys. Lett. 81, 2376 (2002).
7. E. M. Kaidashev, M. Lorenz, H. v. Wenckstern, A. Rahm, H.-C. Semmelhack, K.-H.
Han,G. Benndorf, C. Bundesmann, H. Hochmuth und M. Grundmann, Appl. Phys. Lett. 82,
3901 (2003).
8. O. Breitenstein, M. Langenkamp: Lock-in Thermography - Basics and Use for Functional
Diagnostics of Electronic Components, Springer-Verlag, Berlin, 2003.
9. L. J. van der Pauw, "A Method of Measuring Specific Resistivity and Hall Effect of Discs
of Arbitrary Shape", Philips Research Reports 13, No. 1, pp. 1-9 (1958).
10. H. von Wenckstern, E. M. Kaidashev, M. Lorenz, H. Hochmuth, G. Biehne, J. Lenzner, V.
Gottschalch, R. Pickenhain, and M. Grundmann., Appl. Phys. Lett. 84, 79 (2004).
Chapter 6
STRUCTURE, MORPHOLOGY, AND
PHOTOLUMINESCENCE OF ZnO FILMS
V. A. Karpina1, V. D. Khranovskyy1, V. I. Lazorenko1, G. V. Lashkarev1, I.
V. Blonsky2 , and V. A. Baturin3
1
Abstract:
Key words:
1.
INTRODUCTION
59
N.H. Nickel and E. Terukov (eds.),
Zinc Oxide A Material for Micro- and Optoelectronic Applications, 5968.
2005 Springer. Printed in the Netherlands.
60
V. A. Karpina et al.
Besides that, ZnO has some interesting properties, which make it very
perspective for functional devices nowadays. ZnO is recognized as a
promising optoelectronic material. Due to a large exciton binding energy the
effective excitonic generation in ZnO could be achieved at and above room
temperature. Indeed, optically pumped stimulated emission and lasing have
been demonstrated at room temperature.4-6 There are papers devoted to the
deposition of p-type ZnO.7,8 Hence, we can expect that effective lightemitting devices based on ZnO will be available in the near future.
As semiconductor ZnO can be used for detection of UV-A (300-400 nm)
radiation. By doping with Mg the range can be tuned to the UV-B and UV-C
(200-280 nm) regions. Such a wide interval of sensing spectra makes
ZnxMg1-xO UV detectors to be used in many applications such as solar UV
radiation monitoring, ultra-high temperature flame detection, and airborne
missile warning systems. There are reports on the design of photoconductive
as well as Schottky types of ZnO UV detectors.2,3
Single crystal ZnO films have a strong optical anisotropy and can be used
for modulation of UV radiation. ZnO modulators with a contrast of 70:1 and
an operation speed of 100 ps have been designed already.9
ZnO appears to be a ferroelectric material upon doping with non-central
ions such as Li+, Mg2+.10,11 The Curie temperature can be increased to 340 K
when ZnO is simultaneously doped with Li and Mg. Recently, ferroelectric
thin films have been used for random access memory and other integrated
ferroelectric devices. So, ZnO can be one of the candidates for such systems.
ZnO is recognized as the best candidate for diluted magnetic
semiconductor in the recently emerged new field of electronics spintronics.
Spintronics (spin + electronics) focuses on the creation of functional devices
based on the control of ferromagnetic ordering by an electric field or light.
According to several theoretical predictions ZnO doped with transition
metals is expected to be ferromagnetic with a high Curie temperature.12,13
Therefore, a transparent ferromagnetic could be created. This shows that
ZnO is a promising material for the development of a variety of functional
devices.
2.
EXPERIMENTAL PROCEDURE
ZnO films were obtained by evaporation of zinc metal from a special cell
in oxygen atmosphere. The deposition was carried out on mica and glass
substrates. The deposited films were annealed at 400 C in air and at higher
temperatures. For RF magnetron sputtering we used pure Zn as a target and a
gas mixture of Ar+O2 for the working gas. These samples were deposited on
glass and sapphire substrates. Zinc acetylacetonate Zn(AA)2 was used as a
61
3.
3.1
62
V. A. Karpina et al.
Figure 1. XRD spectra of ZnO/glass films deposited by RTE (a), MS (b), PEMOCVD (c).
ZnO films deposited by RTE and PEMOCVD were annealed at 500 C for two hour in
ambient air.
a)
b)
Figure 2. SEM images of ZnO deposited by RTE (a) and MS (b) and annealed at 500 C.
63
Fig. 1
Figure 3. AFM images of ZnO films deposited by MS at Ts=150 C (a) and Ts=300 C (b).
64
V. A. Karpina et al.
Fig. 2
Figure 4. AFM images of ZnO deposited by PEMOCVD at Ts=200 C (a), Ts=250 C (b), and
Ts=300 C (c).
3.2
65
9.5 MW/cm
m
PL intensity,a.u.
0,8
0,6
0,4
0,2
0,0
1.2 MW/cm
m
-0,2
350
400
450
500
550
600
Wavelength,nm
a)
4 .4 M W /c m
3 ,5
PL intensity, a.u.
3 ,0
2 ,5
2 ,0
1 ,5
1 ,0
1 .2 M W /c m
0 ,5
0 ,0
-0 ,5
340
360
380
400
420
W a ve le n g th , n m
b)
Figure 5. RT PL spectra of ZnO films deposited by MS on sapphire (a) and RTE on mica (b).
ZnO films were annealed at 500 C.
66
V. A. Karpina et al.
77 K
2.5
10
77 K
8
Io
1.5
PL intensity, a.u.
PL Intensity, a.u.
2.0
0.5Io
0.1Io
1.0
0.5
0.0
6
4
2
0
360
380
400
420
Wavelength (nm)
Wavelength (nm)
(b)
(a)
Integrated PL intensity
35
30
25
20
15
10
5
0
10
2
Figure 6. PL spectra of ZnO films deposited on mica by RTE and annealed at 500 C (a), at
high excitation by an N2 laser (b), and at low excitation with a He-Cd laser (c). PL integrated
intensity at N2 excitation.
At low intensity there are two PL lines located at 369 and 374 nm (Fig.
6a). As the excitation intensity increases the intensity of the peaks also
increases. These two bands widen and overlap. Besides, the long-wave part
of the spectra strongly widens. At low excitation using a He-Cd laser we
67
observed three lines (Fig. 6b). Two lines are located at the same positions
and the third line is positioned at 382.5 nm. The energy difference between
the second and third line is 73 meV, what corresponds to LO phonon replica.
Whereas we observed a broadening of the line at increasing of pumping
intensity it could be supposed that stimulated emission with simultaneously
generation of acoustic phonons occurs. Besides that we also suggest that
twofold line is obliged to the recombination of excitons bound on two
different types of defects.
Finally, we have measured the time-dependent behavior of the PL
spectra at 77 K. At a delay time of 0-2 ns we observed strong UV PL (Fig.
7). At steady state conditions a deep-level yellow-red emission connected
with defect luminescence is observed. So it could be concluded that exciton
recombination is a faster process than recombination via deep levels.
normalized PL intensity
1.2
1.0
1 ns
steady state
0.8
0.6
0.4
0.2
0.0
300
400
500
600
700
800
Wavelength (nm)
4.
SUMMARY
ZnO films were deposited by RTE, MS, and PEMOCVD. In the case of
RTE the ZnO films have a large-grained structure with no preferred
crystallite orientation. On the contrary, ZnO films deposited by PEMOCVD
and MS have fine-grained structure with an average grain-size of 40-60 nm.
Under certain growth conditions ZnO films show a preferred orientation with
the c-axis normal to the substrate. However, ZnO films prepared by
68
V. A. Karpina et al.
ACKNOWLEDGEMENTS
We are thankful to A. Yu. Pyryatinsky, V. N. Kadan, A. N. Komarov, V.
M. Vereschyak, and P. M. Lytvyn for their interest and input to this paper.
We are indepted to the Institute of Single Crystals of the National Academy
of Sciences of the Ukraine for providing us with sapphire substrates. This
work is partially supported by the program "Nanosystems, nanomaterials,
nanotechnologies" (project H27-04).
REFERENCES
1. P. P. Gorbick, A. A. Dadykin, and I. V. Dubrovin, Reports of The National Academy of
Science of Ukraine, 1, 84 (2001).
2. W. Yang, R. D. Vispute, S. Choopun, R. P. Sharma, T. Venkates, and H. Shen,
Appl.Phys.Lett. 78, 2787 (2001).
3. S. Liang, H. Sheng, Y. Liu, Z. Huo, Y. Lu, and H. Shen, J. Cryst. Growth, 225, 110
(2001).
4. Z. K. Tang, G. K. L. Wong, P. Yu, M. Kawasaki, A. Ohtomo, H. Koinuma, and Y.
Segawa, Appl. Phys. Lett. 72, 3270 (1998).
5. H. Cao, Y. G. Zhao, H. C. Ong, S. T. Ho, J. Y. Dai, J. Y. Wu, and R. P. H. Chang, Appl.
Phys. Lett. 73, 3656 (1998).
6. Y. Chen, H. Ko, S. Hong, and T. Yao, Appl. Phys. Lett. 78, 1469 (2001).
7. Y. R. Ryu, S. Zhu, D. C. Look, J. M. Wrobel, H. M. Jeong, and H. W. White et al., J.
Cryst. Growth, 216, 330 (2000).
8. M. Joseph, H. Tabata, and T. Kawai, Jpn.J.Appl.Phys. 38, L1205 (1999).
9. M. Wraback, H. Shen, S. Liang, C. R. Gorla, and Y. Lu, Appl.Phys.Lett. 74, 507 (1999).
10. A. Onodera, N. Tamaki, Y. Kawamura, T. Sawada, and H. Yamashita, Jpn. J. Appl. Phys.
35, 5160 (1996).
11. M. Joseph, H. Tabata, and T. Kawai, Appl. Phys. Lett. 74, 2534 (1999).
12. T. Dietl, H. Ohno, F. Matsukura, J. Gibert, and D. Ferrand, Science 287, 1019 (2000).
13. K. Sato and H. Katayama-Yoshida, Jpn. J. Appl. Phys. 39, L555 (2000).
Chapter 7
OPTICS AND SPECTROSCOPY OF POINT
DEFECTS IN ZnO
Vladimir Nikitenko
Moscow State Transport University, Obraztsova str. 15, 103055 Moscow, Russia
Abstract:
Key words:
1.
INTRODUCTION
70
2.
Vladimir Nikitenko
EXPERIMENT
3.
71
400
1
J (arb. units)
JD(Zni )
(a)
300
2
JD(Zni)ohQ=E(JD)-Ed(1-1/n
200
Ed=0.05 eV
100
0
1
J (arb. units)
100
JD(Lii )
(b)
2
JD(Lii)ohQ=E(JD)-Ed(1-1/n
1
JD(Zni)
50
3700
3750
Ed=0.033 eV
3800
3850
O ()
Figure 1. Photoluminescence spectra of a hydrothermally grown ZnO single crystal (a) before
and (b) after Li-doping by heat treatment in Li2CO3 at 700 qC; T=4.2 K; laser excitation at
O=2950 .
72
Vladimir Nikitenko
moment of the oxygen ion O| with the magnetic moment of the nucleus of
isotope 7Li (j = 3/2).
In a non-excited crystal, the substituting lithium ion is surrounded by four
oxygen ions O|| (with two electrons). The center becomes paramagnetic LiZn
(in a more complete notation [Li+Zn O|] when a hole is trapped (or an
electron is released) by one of the nearest neighbor oxygen ions (or by one of
the four oxygen ligands).
As an example needed for further discussions we give Fig. 2. Optical
transitions (3) and (4) during excitation of the sample form paramagnetic
states of lithium at T = 30 K. The capture of free electrons by the traps keeps
N, I (arb. units)
100
80
60
40
20
0
1(LiZn-O')
2(TL)
0
T (K)
73
8000
7192
6989
6000
7209
7240
7307
3000
7121
7135
4000
7054
5000
6989
J (arb. units)
7000
2000
1000
0
7000
7200
7400
7600
7800
O ()
igure 3. Red photoluminescence spectrum of grown hydrothermally ZnO(Fe) single crystal
annealed in Li2CO3 at 973 K for seven days; T=4.2 K; excitation with a DKSSh-3000 lamp
(O=3500 ).
74
Vladimir Nikitenko
N, J (arb. units)
1(LiZn-O')
200
3+
3(FeZn)
2(LiZn-O')
150
5(TL)
4(TL)
100
6(TL)
x10
50
0
50
100
150
200
T (K)
Figure 4. Dependence of the spin density N of the paramagnetic centers and of the
thermoluminescence intensity (TL): (LiZn+-O|) at =30 K (1) and at =70 K (2), FeZn3+ at
=70 K (3) on the annealing temperature of the hydrothermal ZnO single crystal and TL
curves the temperature of excitation =5 K (4) and 80 K (5) TL at O=5900 ; at =70 K
and O=7200 (6). The original single crystal (1,4,5); ZnO:Fe single crystal annealed in
Li2CO3 at T = 973 K for seven days (2,3,6). Oexc. = 3800 .
75
76
Vladimir Nikitenko
luminescence, the optical depth of the lithium level is 1-1.3 eV above the
valence band.1-4 Thus, thermal and optical measurement yield different
defect ionization and neutralization energies for ZnO which is a natural
consequence of the high bond ionicity and may result from a charge-transferinduced polarization of the crystal lattice.
The spectrum of the red luminescence has a complex structure. We were
able to separate three components of red luminescence that each have their
own kinetics and that differ by the character of excitation. The radiation band
at 7300-7800 is predominantly excited by light with O1 = 3700-4800 and
has a three exponential time dependence of afterglow with time constants
of 3 msec, 15 msec, and 75 msec. The relatively strong emission line at 7175
has maxima in the excitation spectrum at 3700 and 5200 . It is
characterized by an afterglow time constant of W = 33 sec. Several narrow
lines in the 6900-7200 region are well excited with O=4000-4800 and
are characterized by W1 = 3 and W2 = 22 msec.
The red luminescence extends into the infra red region (Fig. 5). The
results obtained point to the requirement of further refinements of the
emission transition scheme in the case of the red luminescence of ZnO.
The model, according to which the red luminescence of ZnO represents
luminescence of ions Fe3+Zn (3d5 configuration) is the most encompassing
Omaxx ()
I (arb. units)
11722
10000
11863
11564
12393
12000
14000
16000
O ()
Figure 5. The infra-red photoluminescence spectrum of a ZnO-Fe,Li single crystal at T=4.2
K; excitation at 4880 (argon laser).
77
on our view now. The red luminescence of ZnO is the analogue of the
well-known orange-red luminescence of manganese ions Mn2+ (for example,
in ZnS), which have the same configuration 3d5 or the red luminescence of
ions Fe3+Zn in ZnS. This radiation corresponds to the optical transitions
4
T1(4G) o 6A1(6S).
Irradiation of the ZnO Li samples by protons leads to the appearance of
a photosensitive EPR signal (f signal) which corresponds to the center with
an axial symmetry and g-factor g|| = 1.9948 and gA = 1.9963. This signal is
very important for understanding the processes of defect formation, since it
leads to gives rise to two alternative interpretations for the microscopic
structure of the intrinsic defects that determine the deviation from the
stoichiometry in ZnO.6
According to the one viewpoint, the f-signal should be assigned to the
oxygen vacancies V+O (F+ center) that are formed by capture of one
electron by anionic vacancies with double positive charge. In this case the
detected hyperfine structure of the signal is explained by the interaction of
the unpaired spin of the electron with the nuclei of the 67Zn isotope (content
is 4.1%, spin I=5/2). The condition for photoexcitation of the f-signal of
EPR is such that the oxygen vacancy VO (F center) represents a deep
donor with a photoionization energy of 2.3 eV. In this case the low
resistivity in ZnO is determined by shallow donors of the interstitial zinc
Zni with an ionization energy Ed = 0.05 eV, which by forming the impurity
band or by recoiling the electrons to the empty band, are responsible for the
specific EPR signal with the factor g=1.96.
From another viewpoint the recorded f-signal of EPR should be
attributed to the interstitial oxygen in the paramagnetic state O|i or Oi, and
VO should be considered as the shallow donor specifying the EPR signal
with the factor g | 1.96 (there is an anisotropy g|| = 1.9557 and gA = 1.9573
at low temperature) and high electrical conductivity of ZnO.
The investigations of the processes of the charge exchange of the centers
by the EPR method in the mode of pulsed heating made it possible to
determine the temperature dependence of the relative concentration of the
paramagnetic centers that specify the f-signal of EPR and of the reference
LiZn-centers in the ZnO-Li single crystals irradiated with protons.
Concurrently, we studied thermoluminescence (Fig. 6). The spectra of the
photoluminescence of crystals and those of the excitation of the F+ and LiZn
paramagnetic centers are shown in Ref. 6.
An important feature of the results obtained is the recapture of the
electrons released from traps at T1 | 30-60 K and T2 | 110-150 K by the
paramagnetic centers that form the f-signal of EPR (Fig. 6, curves 1). Only
when these centers finish their work does the transformation of the lithium
paramagnetic centers into the diamagnetic state begin (T2 | 150-180 K,
78
Vladimir Nikitenko
79
2(LiZn-O')
100
N, I (arb. units)
80
+
1(F )
60
4(red TL)
40
20
0
3(y-oTL)
50
100
150
200
250
T (K)
Figure 6. Dependence of the concentration of paramagnetic centers N at T = 30 K (full
symbols) and of the intensity of thermoluminescence (curves 3,4) on the temperature of
annealing of the ZnO single crystal irradiated with protons (E = 7.7 MeV. Tdose 1015
protons/cm2): 1) F+-centers; 2) paramagnetic states LiZn; 3 and 4) the intensity of
thermoluminescence at = 6000 and 7200 (the excitation hv = 2.6 eV).
example, [FMg]+ and [FLi] in CaO or [FCu] in ZnS). The presence of the
impurity in the FA -center is usually confirmed by the hyperfine interaction of
the spin moment of the unpaired electron with the magnetic moment of the
impurity nucleus, which is fixed in the ZnO-Li crystals irradiated with
electrons.6,7
The values of the g-factor (g|| = 1.9951 and gA = 1.9956), hyperfine
interaction constants, and the electron wave function for FA-centers in ZnO
correlate well with similar parameters for the FA- and F+-centers in CaO and
ZnS. Here, LiZn seems to be the impurity in the FA-center in the crystals we
have investigated. This assumption is based on the fact that the investigated
hydrothermal single crystals were doped with lithium with a concentration of
NLi | 1017 cm-3 during the process of growing. The unpaired electron of the
FA-center interacts with the impurity having nuclear spin J = 3/2, which is
also well satisfied for lithium, which in the natural state contains 92.58% of
the isotope 7Li with J = 3/2.
In the temperature range 530-660 K in the investigated ZnO crystals,
thermally stimulated ionic processes were detected, at the first stage of
which (T = 530-630 K) "pure" anionic vacancies disappear and at the second
80
Vladimir Nikitenko
Figure 7. Energy scheme of basic point defects in zinc oxide. The figure shows the energy
level position relative to the edge of the conduction and valence band; the optical and thermal
depths of the corresponding level are given in parentheses and square brackets, respectively
(eV); T = 77 K. Luminescence: violet (V), yellow-range (YO) and green (g). The unbroken
line is the energy level occupied with one electron; the dotted line is the energy level of the
empty state.
optical transition on the impurity, the exciton - like state of two particles is
formed, 3d10(Cu+)h , the so-called acceptor exciton, whose formation and
breakup leads to the appearance of the characteristic structure of edge
absorption (the , and lines) and of green multi phonon luminescence.8-10
It should be noted that in recent years the concept of acceptor and donor
excitons bound with the impurities of transition metals in semiconductors has
found theoretical and experimental justification for a variety of compounds
A2 B 6 .
The given model is confirmed by the application of green emission
intensity in low voltage cathodoluminescent phosphors, single crystals and
ZnO films doped with Cu.11
4.
81
CONCLUSION
ACKNOWLEDGMENT
The author is grateful to . E. Tarkpea (investigations of EPR), A.
Freiberg, I. P. Kuz'mina, S. V. Mukhin and I. V. Pykanov for their help in
this work.
REFERENCES
1. W. Hirschwald, P. Banasewicz, L. Ernst, Curr. Top. Mater. Science. 7, 143 (1981).
2. J. P. Kuz`mina, V.A. Nikitenko, Zinc Oxide. Production and Optical Properties [in
Russian], Moscow (1984).
3. V. A. Nikitenko, J. Appl. Spectroscopy , 57, 783 (1992).
4. D. Zwingel, J. Luminescence. 5, 385 (1972).
5. V. A. Nikitenko, S. V. Mukhin, J. P. Kuz`mina, V. G. Galstyan, S. G. Stoyukhin,
Inorganic Materials, 31, 1243 (1995).
6. V. A. Nikitenko, K. E. Tarkpea, I. V. Pykanov, S. G. Stoyukhin, J. Appl. Spectroscopy, 68,
502 (2001).
7. K. Tarkpea, A. Ots, V. A. Nikitenko, J. Phys. Chem. Sol. 55, 1353 (1994).
8. R. Dingle, Phys. Rev. Lett. 23, 579 (1969).
9. D. J. Robbins, S. G. Bishop, J. A. Savage, P. Porteous, J. Phys. C. 14, 2847 (1981).
10. D. J. Robbins, D. C. Herbert, P. J. Dean, J. Phys. C. 14, 2859 (1981).
11. Ya.J. Alivov, M. V. Chukichev, V. A. Nikitenko, Phys. and Tech. Semiconductors, 38, 34
(2004).
Chapter 8
WHISPERING GALLERY MODES IN
HEXAGONAL ZINC OXIDE MICRO- AND
NANOCRYSTALS
Thomas Nobis, Evgeni M. Kaidashev, Andreas Rahm, Michael Lorenz, and
Marius Grundmann
Universitt Leipzig, Fakultt fr Physik und Geowissenschaften, Institut fr Experimentelle
Physik II, Linnstrae 5, D-04103 Leipzig, Germany
Abstract:
The resonator properties of zinc oxide (ZnO) micro- and nanocrystals grown
by a novel high pressure pulsed laser deposition process have been
investigated at room temperature by cathodoluminescence (CL), spatially
resolved CL-imaging and polarization resolved micro-photoluminescence (PL) within the visible spectral range. The spectra exhibit a series of
comparatively sharp and almost equidistant resonance lines. Using a simple
plane wave interference model and taking into account the spectral
characteristic n(Z) of the refractive index of ZnO, we can unambiguously
attribute those lines to whispering gallery modes (WGMs) of a two
dimensional hexagonal resonator. The predicted resonator diameters agree
well with the measured crystal sizes. Tapered, high aspect ratio ZnO
nanoneedles furthermore allow systematic investigations of the WGMs as a
function of cavity diameter D down to zero. Hence, the transition from a
multi-mode to a single mode cavity is directly observed. -PL experiments
demonstrate that the WGMs are mainly TM polarized.
Key words:
Whispering
polarization
1.
gallery
modes,
nanostructure,
microcrystal,
nanocavity,
INTRODUCTION
84
2.
2.1
Dielectric resonators
85
equation for the given geometry. Steady state solutions of this equation are
then usually called eigenmodes of the considered resonator.
In the case of a hollow metallic resonator as an example for a so-called
closed
d cavity, a discrete set of eigenmodes exists related to a discrete
spectrum of eigenmode energies.10 This is due to the boundary conditions for
the electric and magnetic field that lead to a complete suppression of wave
propagation outside the resonator. Hence, no energy is radiated from the
hollow resonator; the whole radiation field is restricted inside the metallic
cavity.
In contrary, a dielectric resonator is an open cavity. Because of the
continuity conditions for the electric and magnetic field at the cavitys
boundary, the radiation field is not exclusively restricted to the resonator
itself, but actually extended to infinity. This results in a permanent radiation
loss of energy and hence, every imaginable field distribution decays
exponentially with a lifetimeW, therefore not representing a steady state
solution.11 The actual solutions of this problem form a continuum of socalled scattering states consisting of an incoming plane wave and an
outgoing scattered wave. Therefore, the eigenmodes of a dielectric cavity
appear as maxima within the continuous spectrum of the total scattering
cross section in dependence of the energy of the incoming wave. Those
maxima, i.e. those incoming wave energies that exhibit the strongest
interaction with the given resonator, are usually called resonance energies or
resonances of the dielectric cavity.11 Due to the above mentioned finite
lifetime, all resonances exhibit a finite line width. The larger the losses of a
given resonance the shorter is the life time and the larger is the resonances
line width.
2.2
86
Figure 1. (a) The light confined within a hexagonal cross section can be assumed to circulate
inside the cavity as indicated by the white arrows. The geometry of the cavity can be
described by the radius of the incircle Ri, by the radius of the circumscribing circle R, or by its
related diameter D = 2R using the geometric relation R = 2Ri/3. Note that the circumference
of the inscribed white hexagon has a length of 6Ri. (b), (c) Two possible types of polarization
of two-dimensional whispering gallery modes, exemplary given for a hexagonal resonator.
For TE modes the electric field is transverse to the longitudinal axis of the resonator, for TM
modes the magnetic field is.
For large mode numbers N > 70 a simple plane wave model (PWM) has
been deduced.16 Its main idea is, that the light wave circulates around
exploiting the process of multiple total internal reflection (TIR, see Fig. 1a)
and finally interferes with itself when having completed one full circulation
within the resonator. To enforce constructive interference the total phase
87
shift of the wave along its path has to be an integer multiple of 2S, i.e. only
entire wave trains are allowed to perform multiple circulations generating a
standing wave. Taking into account the polarization-dependent negative
phase shift that occurs during the process of TIR,10 we obtain the following
equation:
h c
n E
a cta
arctan
t
(1)
EN v
1
Ri
(2)
and thus, WGMs generally shift to higher energies with decreasing cavity
diameter. Note that, although TIR suggests no way for the light to leave the
resonator, emission nevertheless occurs at the corners of the hexagon.16 After
all, radiation losses are always allowed due to the reasons given in part 2.1.
Therefore it is possible to optically detect WGMs. Since the geometrical
model leading to eq. (1) is a combination of basic ray and plane wave optics,
it is expected to be valid only for R >> O, i.e. for N >> 1.16 Anyhow, we show
in this work, that eq. (1) describes the experimentally found WGMs
even down to mode number N = 1.
88
3.
3.1
Sample properties
3.1.1
Growth
Luminescence
89
Figure 2. Different shapes of PLD grown ZnO micro- and nanocrystals. (a) Hexagonal
microcrystals with diameters of about 2 m. (b) ZnO nanowires with diameters around 100
nm. (c) Needle like nanostructures. (d) Microcrystals exhibiting a dodecagonal cross section.
3.2
3.2.1
90
Figure 3. Room temperature VIS emission CL spectra of a ZnO microcrystal (curve 1) and of
thin film material (curve 2) scaled to the same intensity and vertically shifted for clarity.
Arrows mark theoretical energetic positions of WGMs due to eq. (1) for TM polarization and
D = 2.86 m. The inset shows a SEM image of the investigated microcrystal. The scalebar has
a length of 2 m. The experimental cavity diameter is Dexp = (2.90 r 0.06) m.
3.2.2
For the investigated crystal, the correct mode numbering is shown also in
Fig. 3, the respective theoretical diameter yields to Dtheory = 2.86 m with a
small minimum-maximum spreading of less than 10 nm. Although theory
predicts both TM and TE polarization, all calculations have been performed
for TM type, since polarization-dependent micro-photoluminescence
experiments showed, that WGMs are preferentially TM polarized (see
Figure 7 below). The required data for n||(E) were obtained from
ellipsometry measurements on PLD grown ZnO thin film samples.20
91
Figure 4. Diagram Ri(E) to determine the correct mode numbering for the spectrum of
Figure 3. Equation (1) enables to calculate the radius of the cavity-incircle Ri=R
= i(E,N)
N out of a
given resonance number N and resonance energy E. The resulting set of curves is given here
for N = 1041 as a set of black lines approximately exhibiting a hyperbolic shape (as
expected from eq. 2). The index of the curves, i.e. the respective resonance number, is shown
on the left side of the scheme marked by black arrows. The experimentally detected resonance
energies are given as white circles. Since the resonances have to be numbered in ascending
order, every white line represents a possible peak numbering. As can be seen from the
diagram there exists only one peak numbering that results in an energy-independent incirlce
radius, denoted by the horizontal line. Hence, the starting number for the peak numbering of
the spectrum in Figure 3 yields to N0=21 with Ri = 1.24 m, this leads to a respective cavity
diameter D=2.86 m.
92
3.3
Although the simple PWM obviously is valid for WGMs with mode
numbers in the range of about N = 2030, i.e. for optical cavities in the
micrometer regime, it does not necessarily need to be valid for
nanostructures. Furthermore, prism-shaped microcrystals only enable
analysis of WGMs for fixed cavity diameter given by the cross section of the
particular microcolumn. Those restrictions can be vanquished by focusing on
nanostructures with a needle-like shape.
3.3.1
ZnO nanoneedles
Figure 5. SEM image of the investigated nanoneedle. All images have been obtained at 45
perspective. (a) SEM image of an array of ZnO nanoneedles containing the reported one
marked by a black rectangle. The scale bar has a length of 10 m. (b) Larger scale SEM
image of the investigated nanoneedle. The needles diameter is continuously tapered
approaching zero at the top. (c) High resolution SEM image indicating the hexagonal cross
section of the needle. The scale bar has a length of 300 nm. (d) Experimentally determined
shape of the needle. The obtained characteristic of the diameter D vs. scan position x can be
fitted using a potential law D v x0.523 r 0.007, leading to a square root like shape of the needle.
93
Detection
94
3.3.3
3.3.4
95
scan along the needles longitudinal axis is given in Fig. 8. In grey scales,
the resonances appear as a set of bright lines, that actually exhibit a
curvature as expected from eq. (2). The theoretical values are shown with
dots and cross symbols in Fig. 8. They are directly calculated without free
parameter using the experimentally determined shape of the needle, i.e. Ri(x
( ),
see Fig. 5, and the above mentioned thin-film data for n||(E) and
nA(E).20 The dots (TM) match the measured resonance energies with a very
good agreement. Additionally, the well fitting blue arrows within Fig. 6a
confirm, that the simple PWM gives a good description of the basic physics
of WGMs even in nano-sized crystal geometries within the limits of the
resonance line widths.
Two dimensional plot of spectra recorded along a line scan on the needles
longitudinal axis. The left vertical axis shows the line scan position x, the right one refers to
the respective needle diameter D. The grey scales refer to the spectral CL-intensity. The
spectral maxima, i.e. the measured WGM energies, appear as bright belts going from the
bottom left corner to the right upper one. With
h decreasing diameter all resonances shift
systematically to higher energies. The white dots give theoretical TM-resonance energy
positions obtained from eq. (1), white crosses give the same for TE-polarisation. Without
adjusting free parameters there is a very good agreement between experiment and theory.
96
3.3.5
Discussion
97
and 8.) These deviations probably reveal the limits of the simple PWM. To
discuss this fact more in detail we note, that the application of numerical
methods to predict resonance energies at the end always depends on a
precise determination of the considered geometry, e.g. Ri, and a good
knowledge of the refractive index n. SE microscopy measurements of the
cavity diameter D is performed with an error of at least 1%. The
ellipsometrically determined values for n|| for bulk material yield an error in
n of about 2%. The actual refractive index of the nanocrystals may possibly
differ additionally. These effects lead to an error in E of about 3%, i.e. 50 80 meV in the considered spectral range, which is indicated as error bars in
Fig. 6a. Compared to the line widths of the WGMs, eq. (1) for TM
polarization gives a very good description of the spectral position of the
resonant modes.
4.
CONCLUSION
ACKNOWLEDGMENTS
This work was supported by the Deutsche Forschungsgemeinschaft
within FOR 522 (Project Gr 1011/12-1).
REFERENCES
1. S.L. McCall, A.F.J. Levi, R.E. Slusher, S.J.Pearton, and R.A. Logan, Appl. Phys. Lett. 60,
289 (1992).
2. S. Chang, N.B. Rex, R.K. Chang, G. Chong, and L.J. Guido, Appl. Phys. Lett. 75, 166
(1999).
3. K.J. Luo, J.Y. Xu, H. Cao, Y. Ma, S.H. Chang, S.T. Ho, and G.S. Solomon, Appl. Phys.
Lett. 77, 2304 (2000).
4. X. Liu, W. Fang, Y. Huang, X.H. Wu, S.T. Ho, H. Cao, and R.P.H. Chang, Appl. Phys.
Lett. 84, 2488 (2004).
98
5. A.F.J. Levi, R.E. Slusher, S.L. McCall, S.J. Pearton, and W.S. Hobson, Appl. Phys. Lett.
62, 2021 (1993).
6. S.M. Spillane, T.J. Kippenberg, and K.J. Vahala, Nature 415, 621 (2002).
7. M. V. Artemyev, U. Woggon, R. Wannemacher, Appl. Phys. Lett. 78, 1032 (2001).
8. I. Braun, G. Ihlein, F. Laeri, J.U. Nckel, G. Schulz-Ekloff, F. Schth, U. Vietze, .
Weiss, and D. Whrle, Appl. Phys. B 70, 335 (2000).
9. U. Vietze, O. Krau, F. Laeri, G. Ihlein, F. Schth, B. Limburg, and M. Abraham, Phys.
Rev. Lett. 81, 4628 (1998).
10. J. D. Jackson, Classical electrodynamics (Wiley, New York, 1999).
11. J.U. Nckel, 2-d Microcavities: Theory and Experiments, in Cavity-Enhanced
Spectroscopies, edited by R.D. van Zee and J.P.Looney, (Academic Press, San Diego,
2002) pp. 185-226.
12. N.S. Kapany and J.J. Burke, Optical Waveguides, (Academic Press, New York, 1972)
13. S.-B. Lee, J.H. Lee, J.S. Chang, H.J. Moon, S.W. Kim, and K. An, Phys. Rev. Lett 88,
033903 (2002).
14. C. Gmachl, F. Capasso, E.E. Narimanov, J.U. Nckel, A.D. Stone, J. Faist, D.L. Sivco, and
A.Y. Cho, Science 280, 1556 (1998).
15. M.Lohmeyer, Opt. Quant. Electron. 34, 541 (2002).
16. J. Wiersig, Phys. Rev. A 67, 023807 (2003).
17. J. Wiersig, J. Opt. A. 5, 53 (2003).
18.A. M. Morales and C. M. Lieber, Science 279, 208 (1998).
19. M. Lorenz et al., unpublished.
20. R. Schmidt, B. Rheinlnder, M. Schubert, D. Spemann, T. Butz, J. Lenzner, E.M.
Kaidashev, M. Lorenz, A. Rahm, H.C. Semmelhack, and M. Grundmann, Appl. Phys. Lett.
82, 2260 (2003).
21. J. Christen, M. Grundmann, and D. Bimberg, J. Vac. Sci. Technol. B 4, 2358 (1991).
22. Th. Nobis, E.M. Kaidashev, A. Rahm, M. Lorenz, J. Lenzner, and Marius Grundmann,
Nano Lett. 4, 797 (2004).
Chapter 9
PROPERTIES OF DISLOCATIONS IN
EPITAXIAL ZnO LAYERS ANALYZED BY
TRANSMISSION ELECTRON MICROSCOPY
E. Mller1, D. Livinov1, D. Gerthsen1, C. Kirchner2, A. Waag3, N. Oleynik4,
A. Dadgar4, and A. Krost4
1
Abstract:
Key words:
1.
INTRODUCTION
100
E. Mller et al.
2.
EXPERIMENTAL TECHNIQUES
101
3.
3.1
102
E. Mller et al.
Figure 1. Cross-section weak-beam images of the same area of a MOVPE-grown ZnO layer
on Al2O3(0001) with a 40 nm GaN buffer layer taken with (a) g = {0002} under g,3g
conditions and (b) g = {11-20} under g,3g conditions.
Figure 2. Bright-field image of the interface between the GaN/Al2O3(0001) interface and the
ZnO layer taken with g = {11-20}. The black line marks a GaN island.
103
Figure 3. Cross-section weak-beam images of the same area of a MOPVE-grown ZnO layer
on Al2O3(0001) with a thick GaN buffer layer taken with (a) g = {0002} and (b) g = (11-20)
taken under g,3g conditions.
104
E. Mller et al.
are characterized by one {1-100} plane inserted from below as expected due
to the larger a lattice parameter of the ZnO.
Figure 4. (a) HRTEM image of the interface between the GaN buffer layer and ZnO along the
<11-20>-zone axis and (b) Fourier-filtered image showing only {1-100} lattice fringes. The
arrows indicate {1-100} fringes inserted from the substrate side.
105
Figure 5. Schematic representation of the orientation of a misfit dislocation (bold black line)
and possible Burgers vectors (grey vectors) in the (0001) interface plane
3.2
R2
r2
l 2 2 1
ln
4SHH 0 r R
Q
(1)
where Q=q e is the charge per dislocation length and H = 8.4 the dielectric
constant of the material (Viswanatha et al., 2004). The Read radius R =
q/
describes the radius of space charge region where Nsc is
density of compensating charges available for screening.
Electron holography is a TEM technique that allows the retrieval of the
phase of the electron wave. This provides the possibility to study the
interaction of the incident electrons with the electrostatic potential of a
charged dislocation with high spatial resolution. Fig.6 indicates the
106
E. Mller et al.
Figure 6. Geometry used for electron holography investigations: the electron beam is oriented
along the z direction and the dislocation along y perpendicular to the electron-beam direction.
CE (V0 V (r )) dz
0
(2)
107
dislocations and electron beam are parallel. The dislocation strain field can
locally tilt the lattice such that dynamical diffraction conditions apply in the
vicinity the dislocation, i.e. Bragg reflections are strongly excited, even if
kinematical conditions are chosen for the undisturbed crystal. The
contribution of dynamical phase shifts can be minimized in our approach as
will be shown later in this section.
On the basis of Eq.(2), a comparison of the phase shift measured by
electron holography and the calculated 'I yields the line charge q. Figure 7
shows the evaluated amplitude and phase of an electron hologram for the
transmitted beam. The dislocation line direction is indicated by the dashed
line in the amplitude image. Due to the noise, the position of the dislocation
is visible only by a weak dark contrast in the phase image.
Figure 7. Reconstructed amplitude and phase of the electron wave in an area containing a
dislocation indicated by the dashed line in the reconstructed amplitude.
108
E. Mller et al.
measured phase change with respect to the undisturbed crystal, the thickness
integration interval in Eq.(2) was chosen to be twice the Read radius.
Therefore, an accurate value for the sample thickness is not required.
Figure 8. Measured phase change perpendicular to the dislocation line and fit according to Eq.
(1) and Eq.(2).
The value for the line charge is reliable only if the effects of dynamical
diffraction can be neglected. To assess the influence of dynamical diffraction
on the phase shift, the worst case was considered by assuming that the
dislocation locally tilts the crystal from kinematical into {0002} or {11-20}
two-beam conditions. The resulting dynamical phase shift was calculated
using the Bloch wave formalism of the EMS program package (Stadelmann,
1987). Figure 9 shows the phase difference for the transmitted beam in the
undisturbed crystal as a function of the thickness comparing the phase under
kinematical conditions and for {0002} or {11-20} two-beam conditions. For
the tilt angle the experimental values (1.5 perpendicular and 5 along the
dislocation line with respect to the <1-100> zone axis) were taken into the
account. It follows from Fig. 9, that dynamical diffraction is negligible if the
distorted region in Bragg condition is smaller than 7 nm. With the
experimental geometry with a dislocation orientation perpendicular to the
electron-beam direction the assumption of a laterally limited distortion field
is more likely to be satisfied compared to a parallel orientation.
109
Figure 9. Phase difference for the transmitted beam in the undisturbed crystal comparing the
phase under kinematical conditions and for (lower curve) {0002} or (upper curve) {11-20}
two-beam conditions as a function of the thickness.
110
E. Mller et al.
4.
SUMMARY
REFERENCES
Alexander, H., and Teichler, H., Dislocations, 1991, in: Materials Science and Technology,
Vol. 4, chap.6, W. Schrter, ed., North-Holland Pub. Company, pp. 249.
Cherns, D., and Jiao, C.G., 2001, Phys. Rev. Lett. 87, 205504.
Dadgar, A., Oleynik, N., Forster, D., Deiter, S., Witek, H., Blsing, J., Bertram, F., Krtschil,
A., Diez, A., Christen, J., and Krost, A., 2004, J. Cryst. Growth 267, 140 .
Gerthsen, D., 1986, phys. stat. sol. (a) 97, 527.
Hong, Soon-Ku, Ko, Hang-Ju, Chen, Y., Yao, T, 2000, J. Cryst. Growth 209, 537.
Hull, D., Introduction to Dislocations, 1975, 2ndd edition, Pergamon Press, pp.122.
111
Kruse, P., Rosenauer, A., and Gerthsen, D., 2003, Ultramicroscopy 96, 11.
Landolt Brnstein, 1987, Springer-Verlag, Heidelberg, 165
Lehmann, M., and Lichte, H., 2002, Microscopy and Microanalysis 8, 447.
Narayan, J., Dovidenko, K., Sharma, A.K., and Oktyabrsky, 1998, S., J. Appl. Phys. 84, 2597.
Osipyan, Yu. A., Petrenko, V.F., Zaretski, A.V., and Whitworth R.W., 1986, Adv. Phys. 35
115.
Read, W.T., 1954, Phil. Mag. 45, 775.
Stadelmann, P.A., 1987, Ultramicroscopy 21, 131.
Strecker, A., Salzgeber, U., and Mayer, J., 1993, Prakt. Metallographie 30, 482.
Vigu, F., Venngus, P., Vzian, S., Lagt, M., and Faurie, J.-P., 2001, Appl. Phys. Lett. 79,
194.
Viswanatha, R., Sapra, S., Satpati, B., Satyam, P.V., Dev, B. N., and Sarma, D.D., 2004, J.
Mater. Chem. 14, 661.
Williams, D.B., and Carter, C.B., Transmission Electron Microscopy, 1996, Plenum Press,
New York.
Chapter 10
MUON SPIN ROTATION MEASUREMENTS ON
ZINC OXIDE
E. A. Davis
Department of Materials Science and Metallurgy, University of Cambridge, Pembroke Street,
Cambridge CB2 3QZ, United Kingdom
Abstract:
The theoretical suggestion that hydrogen might form shallow donor states in
zinc oxide and hence account for the n-type conductivity normally found in
undoped samples has been confirmed by experimental studies using
mounium, a light pseudo-isotope of hydrogen. Characteristic frequencies in
muon spin rotation experiments yield a hyperfine constant that is ~104 times
smaller than that of vacuum-state muonium, indicating an extended orbital and
a shallow centre. Temperature-dependence studies yield an ionization energy
of about 30 meV. Band-offset diagrams and measurements on other
semiconductors suggest that hydrogen forms shallow donor centres when the
electron affinity of the host material is more than about 4 eV; otherwise the
hydrogen level lies deep in the energy gap.
Key words:
1.
INTRODUCTION
116
E. A. Davis
hydrogen and so its Bohr radius and ionization energy are essentially the
same as those of hydrogen.
Experiments using muons are generally easier and more direct than
experiments1 using hydrogen itself. One reason for this is the high sensitivity
of the technique, by which one can essentially see individual muons.
Another and more fundamental aspect is that hydrogen frequently forms a
U system, making the neutral state inaccessible under thermal
negative-U
equilibrium conditions. Such conditions do not apply in the muon
experiments, permitting observation of both the ionized (charged) and
unionized (neutral) states.
The theoretical prediction2 that hydrogen forms a shallow donor state in
ZnO has been confirmed by muon implantation studies3,4. Below 40 K a
distinctive beating of the muon precession signal provides the required
signature. Fourier transformation reveals a central line at the muon Larmor
frequency accompanied by two symmetrically disposed satellites, the
separation of which yields a hyperfine constant about 10,000 times smaller
than that for vacuum-state muonium. Confirmation that this signal represents
the extended orbital of a shallow centre has been obtained by studies as a
function of temperature, which reveal the disappearance of the satellites with
an ionization energy of approximately 30 meV. More recent measurements5
in longitudinal magnetic fields provide further evidence for these findings.
2.
The experiments were undertaken at the ISIS pulsed muon facility at the
Rutherford Appleton Laboratory situated near Oxford in England. Positive
muons of energy 4 MeV and 100% spin polarized are brought to rest within a
few tenths of a millimetre in powder or single-crystal samples mounted
within a cryostat.
The technique of muon spin rotation involves applying a magnetic field
perpendicular to the direction of the incoming beam of muons (transverse to
their spin) and monitoring the resulting precession signal via the emission of
positrons that are emitted preferentially in the direction of the muon spin at
the moment of its radioactive decay. For the bare muons this is simply the
Larmor frequency but for muonium several frequencies are observed. In the
case of a small hyperfine constant, one can easily reach the so-called
Paschen-Back regime in moderate fields and then a triplet of lines is seen in
a Fourier transform of the raw data.
Asymmetry (%)
117
20
0
-20
0
Probability Density
0.06
0.05
10
Time (s)
15
20
ZnO powder
5K
0.04
0.03
0.02
0.01
0.00
2.0
2.2
2.4
2.6
2.8
3.0
3.2
3.4
Frequency (MHz)
Figure 1. Muonium spin precession signal for a ZnO powder sample at 5 K. The upper plot is
the raw time-domain spectrum (corrected for the muon decay) while the lower plot is the
corresponding frequency spectrum. The central line corresponds to the Larmor frequency of
the bare muon (ionized muonium) and the two symmetrically disposed satellites are associated
with muonium. The dotted curve is a theoretical fit using a powder-pattern lineshape.
118
E. A. Davis
100 60
40
20
10
Asymmetry (%)
20
26(7) meV
2.5(9) meV
15
Diamagnetic fraction
10
10
20
30
40
50
60
70
10
100
00
800 60
100
60
12
2
12.5
ZnO powder
TF 200G
25
Asymmetry (%)
T (K)
80
-1
1000 / T (K )
40
40
T (K)
20
20
12.5
ZnO powder
TF 200G
2.4(8) meV
33(12) meV
4
2
0
10
Paramagnetic fraction
20
30
40
50
60
70
80
-1
1000 / T (K )
Figure 2. Arrhenius plots of the temperature dependences of the amplitudes of the lines shown
in Figure 1. (Courtesy of J M Gil.)
3.
FURTHER ANALYSIS
Assuming for the sake of simplicity an isotropic centre, the effective Bohr
radius a
can be obtained directly from the hyperfine constant A
,
since this scales as the third power of the radius:
13
a0 A0 A
119
(1)
a0H me m
(2)
where H is the relative permittivity and m
the effective mass for electrons
in ZnO. Using the values H = 8 and m
me = 0.24, gives a
= 1.7 nm, in
fair agreement with the value deduced from the data.
In a similar vein we can estimate the ionization energy from the
hyperfine constant:
I
13
I 0 A A0 / H
(3)
where I 0 is the Rydberg = 13.6 eV. This yields I
= 51 meV, which can be
compared with the hydrogenic value, given by:
I
I 0 m me / H 2
(4)
4.
SINGLE-CRYSTAL STUDIES
D
E
2
(5)
120
E. A. Davis
Figure 3. Frequency spectrum derived from muonium spin rotation measurements on a singlecrystal sample of ZnO as a function of orientation (from reference 5).
Of the possible sites that hydrogen (muonium) can occupy in the wurtzite
lattice (see Figure 4) the BCparr site has been shown theoretically to have the
lowest energy, i.e. to be the most stable site9. Experimentally we cannot
distinguish between this site and the site antibonding to oxygen along the caxis.
BCpar
BCperp
ABApar
ABAperp
ABCpar
ABCperp
Figure 4. Six principal sites for hydrogen in the wurtzite structure. BC = bond centre and AB
= antibonding. Superscripts denote sites adjacent to the cation(C) and anion (A). Perp and
par refer to directions perpendicular (not exactly) and parallel to the c-axis (from ref. 11).
121
5.
1.2
SiO2
1.0
AlN
MgO
ZnSe
ZnS
0.8
GaAs
GaP
SiC
0.6
Ge
Si
0.4
0.2
CdTe/S/Se
ZnO
InN
0.0
Figure 5. Correlation of the normalized muonium hyperfine constant with electron affinity.
(Courtesy of S F J Cox reference 11.)
122
E. A. Davis
such as SiO2 and diamond on the left (for which A
A0 is close to unity) to
semiconductors in the middle of the plot. The five semiconductors with an
effectively zero (on this scale) hyperfine constant are those in which
muonium forms shallow centres. It is evident that the dilation of the
wavefunction from atomic-like to extended occurs rather suddenly on this
plot, corresponding to an electron affinity of about 3.7 eV.
The concept of the electron affinity of the host being the all important
factor influencing the deep to shallow transition is implicit in the band-offset
diagrams proposed by Van de Walle12. An example is shown in Figure 6 for
a few materials. In this diagram Ec represents the energy of the bottom of the
conduction band and Ev the top of the valence band, both plotted on an
absolute energy scale, i.e. with respectt to the vacuum level. The dashed line
marked +/- gives the energy at which the formation energies of the positively
and negatively charged states of hydrogen are equal. This represents the
position at which the Fermi level would be pinned in a negative U system.
If this level lies in the conduction band (as in ZnO and InN) hydrogen
forms a shallow donor level. If the level lies in the band gap (as for the other
three materials) then hydrogen forms a deep centre. The predictive nature of
this model is currently being tested for other semiconductors, in particular
oxides other than ZnO.
Peacock and Robertson14 have questioned whether the hydrogen level
really does lie at a constant depth below the vacuum level. A subtle point
relevant to this question is that the donor state actually lies at the level +/0,
i.e. the energy through which the Fermi level would pass when the centre
ionizes, changing from the neutral to the positively charged state. This level
differs by U
U/2 from the +/- level and so, even if the +/- level is invariant, the
123
+/0 level is not expected to be so. Under thermal equilibrium conditions the
Fermi level is pinned at +/- (for a ve U system). Such conditions do not
apply in the muonium experiments and so we are able to explore the higher
+/0 level directly, without the need to illuminate the sample to reveal the
neutral state, as required for, say, ENDOR experiments using hydrogen
itselff1.
ACKNOWLEDGEMENTS
The author wishes to acknowledge fruitful and enjoyable collaboration
with several colleagues during the experimental work reported here. These
are S F J Cox (Rutherford Appleton Laboratory and University College
London), P J C King, J S Lord and S P Cottrell (Rutherford Appleton
Laboratory), J M Gil, H Alberto, R Vilo, J Piroto Duarte and N Ayres de
Campos (Coimbra University, Portugal), and R Lichti (Texas Tech
University, USA).
REFERENCES
1. D. M. Hofmann, A. Hofstaetter, F. Leiter, H. Zhou, F. Henecker, B. K. Meyer, S. B.
Orlinskii, J. Schmidt, and P. G. Baranov, Phys Rev Lett. 88, 045504 (2002).
2. C. G. Van de Walle, Phys. Rev. Lett. 85, 1012 (2000).
3. S. F. J. Cox, E. A. Davis, S. P. Cottrell, P. J. C. King, J. S. Lord, J. M. Gil, H. V. Alberto,
R. C. Vilo, J. Piroto Duarte, N. Ayres de Campos, A. Weidinger, R. L. Lichti, and S. F.
C. Irvine, Phys. Rev. Lett. 86, 2601 (2001).
4. J. M. Gil, H. V. Alberto, R. C. Vilo, J. Piroto Duarte, P. J. Mendes, L. P. Ferreira, N.
Ayres de Campos, A. Weidinger, J. Krause, E. A. Davis, S. P. Cottrell, and S. F. J. Cox,
Phys. Rev. B 64, 075205 (2001).
5. H. V. Alberto, R. C. Vilo, J. Piroto Duarte, N. Ayres de Campos, R. L. Lichti, E. A. Davis,
S. P. Cottrell, and S. F. J. Cox, Hyperfine Interact. 136/137, 471 (2001).
6. D. C. Look, D. C. Reynolds, J. R. Sizelove, R. L. Jones, C. W. Litton, G. Cantwell, and W.
C. Harsch, Solid State Commun. 105, 399 (1998).
7. D. C. Reynolds and T. C. Collins, Phys. Rev. 185, 1099 (1969).
8. D. C. Look, C. Coskun, B. Clafin, and G. C. Farlow, Physica B 340-342, 32 (2003)
9. E V. Lavrov, J. Weber, F. Brrnert, C. G. Van de Walle, and R. Helbig, , Phys. Rev. B 66,
165205 (2002).
10. K. Shimomura, K. Nishiyama, and R. Kadono, Phys Rev Lett. 89, 255505 (2002).
11. S. F. J. Cox, J. Phys.: Condens. Matter 15, R1727 (2003).
12. C. G. Van de Walle and J. Neugebauer, Nature 423, 626 (2003).
13. E. A. Davis, S. F. J. Cox, R. L. Lichti, and C. G. Van de Walle, Appl. Phys. Lett. 82, 592
(2003).
14. P. W. Peacock and J. Robertson, Appl. Phys. Lett. 83, 2025 (2003).
Chapter 11
HYDROGEN DONORS IN ZINC OXIDE
M. D. McCluskey and S. J. Jokela
Department of Physics, Washington State University, Pullman, WA 99164, U.S.A.
Abstract:
Zinc oxide (ZnO) has emerged as a leading material for micro- and
optoelectronic applications.
Although the fabrication of ZnO, from
nanocrystals to bulk single crystals, is well established, a major roadblock for
fabricating optoelectronic devices is the lack of reliable p-type doping. The
presence of compensating donors inhibits the growth of p-type ZnO. In this
paper, studies pertaining to the microscopic structure and doping properties of
hydrogen in ZnO are described. Results from infrared (IR) spectroscopy are
consistent with a model where the hydrogen attaches to a host oxygen atom, in
an anti-bonding orientation, which is not aligned along the c axis. These
hydrogen complexes are unstable, however, perhaps due to the formation of
H2 molecules.
Key words:
1.
INTRODUCTION
126
2.
EXPERIMENTAL TECHNIQUES
Single crystal, c-cut ZnO samples were purchased from Cermet, Inc.
Samples were sealed in quartz ampoules along with hydrogen or deuterium
gas. In addition, ZnO powder was placed in the ampoule, in order to prevent
the decomposition of the ZnO crystal due to hydrogen reduction. Hydrogen
diffusion occurred at a temperature of 700-800C for 10-100 hr, followed by
a rapid quench to room temperature by immersion of the ampoule in water.
The sample was then retrieved by breaking the ampoule.
IR spectra were taken at room temperature (300 K) and liquid-helium
temperatures (5-15 K), using a Bomem DA8 Fourier transform infrared
(FTIR) spectrometer and an InSb detector. For the low-temperature
measurements, a Janis continuous-flow liquid-helium cryostat with wedged,
IR-transparent windows was utilized. Hall-effect measurements, in the Van
der Pauw geometry, were performed at room temperature using a system
from MMR Technologies. Wires were attached to the ZnO using silver paint,
which provided adequate Ohmic contacts for the electron concentrations
(1017 cm-3) in these samples.
3.
RESULTS
3.1
127
Absorbance
0.15
0.10
T = 300 K
0.05
T=8K
0.00
3300
3320
3340
3360
-1
1
1/ N 0
At
(1)
0.25
Absorbance
-3
4.4
4.0
17
0.20
3.6
3.2
0.15
2.8
0.10
2.4
0.05
128
2.0
0
200
400
600
800
1000
0.00
Time (hr)
Figure 2. Electron density (circles) and IR absorbance of the O-H peak (squares) as a function
of time, in hydrogenated ZnO.
3.2
Polarized IR spectroscopy
129
3.3
High-pressure measurements
[0001]
Zn
O
H
Figure 4. Ball-and-stick diagrams of the antibonding (left) and bond-centered (right) models
for hydrogen donors in ZnO.
130
1.1 GPa
0 GPa
3300
3320
3340
3360
-1
Wave numbers (cm )
4.
UNSOLVED MYSTERIES
131
132
5.
CONCLUSIONS
ACKNOWLDGMENTS
The authors would like to thank S. Limpijumnong, M. Stavola, C.G. Van
de Walle, and S.B. Zhang for helpful discussions. This work was supported
by NSF Grant No. DMR-0203832.
REFERENCES
1.
2.
3.
4.
5.
Chapter 12
HYDROGEN-RELATED DEFECTS IN ZnO
STUDIED BY IR ABSORPTION SPECTROSCOPY
E.V. Lavrov,1,2 F. Brrnert,1 and J. Weber1
1
Technical University Dresden, 01062 Dresden, Germany; 2Institute of Radioengineering and
Electronics of RAS, 101999 Moscow, Russia
Abstract:
Most of the ZnO growth techniques result in n-type conductivity of the crystal,
which hinders the progress of ZnO applications for electronic devices.
Hydrogen incorporated during the process of crystal growth is now considered
as a likely source of the n-type conductivity of ZnO. Infrared absorption
spectroscopy provides detailed insights into the physical properties of the light
impurities in the semiconductor matrix and is, therefore, an excellent tool to
explore the structure of the hydrogen-related defects embedded in ZnO. We
report a number of hydrogen-related defects observed in hydrothermally
grown ZnO and ZnO grown from the vapor phase studied by Fourier
transform infrared absorption spectroscopy. Three IR absorption lines at
3611.3, 3349.6, and 3312.2 cm1 at 10 K are observed after hydrogenation of
the vapor phase grown ZnO. The line at 3611.3 cm1 is tentatively assigned to
a bond-centered hydrogen, whereas the other two are identified as a Zn
vacancy decorated with two hydrogen atoms. The nature of a hydrogen center
responsible for the 3577.3 cm1 line observed at 10K in as-grown
hydrothermal ZnO is discussed.
Key words:
1.
INTRODUCTION
134
E. V. Lavrov et al.
However, a new cause for the n-type conductivity of ZnO has recently been
proposed: First-principles investigations based on density functional theory
suggest that hydrogen in ZnO occurs exclusively in the positive charge state,
i.e., it always acts as a donor.4 Recent muon spin rotation (PSR)5 and EPR6
studies confirmed the presence of hydrogen-related shallow donors.
However, the molecular structure of these defects remains unknown.
Local vibrational mode (LVM) spectroscopy is a tool suited to study this
problem. The knowledge about LVMs supplies a detailed information about
the physical properties of light impurities embedded in ZnO. The
frequencies of the LVMs reveal directly the chemical binding of hydrogen
with its neighbors due to the dependence on the atomic structure of the
hydrogen-related defects.
To our knowledge, the first infrared absorption study of hydrogen-related
LVMs in ZnO was made in the end of the 70th by Grtner and E. Mollwo.7,8
They investigated copper doped ZnO annealed in a hydrogen atmosphere
and identified a set of LVMs originating from different copper-hydrogen
defects.
Interest for hydrogen behavior in ZnO stimulated further studies of
LVMs hydrogen-related defects. McCluskey et al. reported an O-H stretch
mode at 3326.3 cm1 observed at 8 K in ZnO samples that were annealed at
700 oC in hydrogen atmosphere.9 The 3326.3 cm1 line was assigned to a
hydrogen atom in an antibonding configuration oriented at an angle of 100o
to the c axis.
Recently, Nickel and Fleischer reported a Raman scattering study of asgrown state-of-the-art ZnO. A number of Raman lines were detected, which
were ascribed to LVMs of N-H and C-Hn complexes with n=1,2,3.10
Here, we summarize our results of an infrared absorption study of
hydrogen-related centers in ZnO grown by different techniques.11,12 A
number of O-H LVMs has been detected. The lines at 3611.3, 3349.6, and
3312.2 cm1 are observed after hydrogen plasma treatment in nominally
undoped ZnO grown from the vapor phase. The 3611.3 cm1 line is
tentatively identified as a LVM of bond centered (BC) hydrogen, whereas
the other two modes are assigned to LVMs of a zinc vacancy decorated with
two hydrogen atoms. The origin of an IR absorption line at 3577.3 cm1
detected at 10 K in as-grown hydrothermal ZnO is discussed.
2.
EXPERIMENTAL DETAILS
Two types of ZnO crystals have been employed in the current study.
Nominally undoped n-type ZnO crystals grown from the vapor phase (VP) at
135
3.
3.1
Results
136
E. V. Lavrov et al.
where Hr and Dr are the reduced masses of the 16O-H and 16O-D species,
respectively. The LVM frequencies strongly suggest that the new lines are
stretch LVMs of the O-H species. When the sample is treated with a mixed
hydrogen and deuterium plasma, four additional lines at 3346.6, 3315.2,
2484.6, and 2463.0 cm1 are observed (bottom spectra in Fig. 1). No extra
lines are seen near the 3611.3 and 2668.0 cm1 lines. This implies, that two
different defects are responsible for the 3611.3, 3349.6, and 3312.2 cm1
lines. The LVM at 3611.3 cm1 belongs to a defect containing one hydrogen
atom, whereas the LVMs at 3349.6 and 3312.2 cm1 originate from a defect
containing two nonequivalentt hydrogen atoms. We label these defects H-I
and H-II.
Figure 2 shows the results of isochronal annealing experiments on ZnO
samples treated with a hydrogen plasma at 350 oC. It follows from the figure
that the H-I defect anneals out at 500 oC, whereas H-II is more stable and
disappears from the spectra at 600 oC. Parallel to the annealing of H-I and HII, another line at 3191.6 cm1 grows in the spectra. It represents a stretch
LVM of the well-known Cu-H defect.7,8 It reaches its maximum intensity at
600 oC, when both H-I and H-II are already gone from the spectra.
Figure 1. Absorption spectra measured for vapor phase grown ZnO samples at 10 K after: H-,
D-, and H+D-plasma treatments at 350 oC. Unpolarized light.
137
Figure 2. The intensity (integrated absorbance) of the 3191.6, 3312.2, 3349.6, and 3611.3
cm1 lines measured at 10 K for a vapor phase grown ZnO sample vs. the annealing
temperature. The sample was exposed to a hydrogen plasma at 350 C.
Bond angles of the O-H species comprising the H-I and H-II defects with
respect to the c axis of the crystal may be obtained from the polarized IR
absorption spectra. For the stretch modes, the direction of the dipole moment
should coincide with the bond direction, provided each bond vibrates
independently.
The H-I defect contains only one O-H bond, therefore the assumption that
its direction coincides with the induced dipole moment is fulfilled
automatically. The H-II defect contains two O-H species and according to
the isotope substitution data (see bottom spectrum in Fig. 1), the coupling
between the two LVMs of the defect is rather weak, i.e. substitution of one
hydrogen atom comprising the defect with deuterium does not change the
LVM frequency of the second O-H bond considerably. Thus, the O-H bonds
of H-II vibrate independently, and from the polarized absorption spectra we
get directly the bond angles.
The results of polarization dependent absorption measurements are
presented in Fig.3. It shows the normalized intensities of the LVMs
138
E. V. Lavrov et al.
belonging to the H-I and H-II defects as a function of the polarizer angle (-)
with respect to the c axis of the sample. As one can see from the figure, the
LVM of H-I and the LVM of H-II at 3312.2 cm1 are polarized along the c
axis, while the LVM of H-II at 3349.6 cm1 is polarized perpendicular to c.
Figure 3. Normalized LVMs intensities of the H-I (a) and H-II (b) defects measured at 10 K
as a function of the angle - between polarizer and the c axis of the sample.
The best fit to the experimental data shown in Fig. 3 by the solid lines is
achieved when the O-H bond of the H-I defect is parallel to the c axis with
the accuracy of 10o. The same result is obtained for the O-H bond
responsible for the 3312.2 cm1 line of H-II, whereas the O-H bond giving
rise to the 3349.6 cm1 line is nearly perpendicular to the c axis, i.e.:
- =100o.
3.2
Molecular models
3.2.1
H-I complex
The experimental observations indicate that this defect has one hydrogen
atom primarily bound to an oxygen atom with the O-H bond aligned close to
the c axis. H-I is observed in all our nominally undoped VP ZnO samples.
We propose that interstitial hydrogen at the BC|| site is the most likely
candidate for this center. Here, the symbol || stands for the configuration
139
with the O-H unit parallel to the c axis; the other bond directions are denoted
by the symbol A. Indeed, the first-principles calculations indicate that in
this configuration hydrogen is strongly bound to an oxygen atom, with the
neighboring Zn atom relaxing outward by almost 40 % of the bond length.4
The results of Ref. 4 showed that the BC||, BCA, AB|| and ABA configurations
all have very similar formation energies, although calculations with
increased convergence and additional relaxation indicate that BC|| may be
lower than the others by about 0.1 eV.11 The experimental results
conclusively indicate a configuration oriented parallel to c, but experiment
alone cannot distinguish between BC|| and AB||. Both the calculated
formation energies4 and the calculated vibrational frequencies11 strongly
favor BC||, indeed, the calculated LVM for BC|| is about 300 cm1 lower than
the LVM at BC||. All these findings point to interstitial hydrogen at BC||
(H+BC) as the most likely candidate for the H-I defect.
Nonetheless, we cannot exclude that this interstitial hydrogen might be
bound to an impurity. Indeed, the fact thatt the H-I signal is stable up to about
350 oC in the annealing experiments (see Fig. 2) may point to a
configuration that is more strongly bound than isolated interstitial hydrogen,
which has a diffusion barrier of about 0.9 eV.15 Interstitial hydrogen bound
to an isoelectronic or acceptor impurity could account for the H-I line,
provided the hydrogen would still be in the BC|| configuration.
It is known that the conductivity of ZnO crystals exposed to hydrogen
increases.16 With our samples we have found that the resistivity is reduced
from 50 to 5 :cm after hydrogen plasma treatment. This finding is definitely
consistent with hydrogen being introduced as an interstitial and acting as a
shallow donor. However, the resistivity data alone cannot help to distinguish
between hydrogen being present as an isolated donor (H+BC) or passivating
an acceptor impurity or native defect, since both scenarios increase the ntype conductivity of ZnO.
3.2.2
H-II complex
140
E. V. Lavrov et al.
with two hydrogen atoms may explain our data. Interstitial oxygen and
oxygen antisites are ruled out by the weak coupling between the two LVMs
of H-II. Therefore, we propose that a zinc vacancy decorated with two
hydrogen atoms, VZnH2, is the most probable model for H-II.11
This assignment is also supported by the results of cathodoluminescence
studies of ZnO after hydrogen plasma treatment. Sekiguchi et al. have shown
that such a treatment made at 400 oC fully suppresses the, so-called, green
luminescence band centered around 2.2 eV.18 We have observed the same
behavior by means of photoluminescence for our ZnO samples. On the other
hand, based on first-principles pseudopotential calculations, Kohan et al.
suggested that VZn is involved in the green luminescence in ZnO.19 Thus,
passivation of VZn with two hydrogen atoms leads to a reduction of the
green luminescence, and appearance of the H-II LVMs.
As mentioned above, the resistivity of our ZnO samples is reduced from
50 to 5 : cm after hydrogen plasma treatment, which is also in favor of the
VZnH2 model. Indeed, passivation of VZn acting as an acceptor in ZnO should
increase the n-type conductivity of the samples.
4.
141
142
E. V. Lavrov et al.
after thermal and hydrogen plasma treatment. The last finding strongly
indicates that the H-I* defect does not include a substitutional Ni atom, as
tentatively suggested in Ref. 12.
We could not establish a correlation/anticorrelation between the 3577.3
cm1 line and any other spectroscopic signal of HT ZnO. We take this as an
indication that the center giving rise to this line is related to the presence of
native defects. The LVM frequency of the center points to VZn, Oi, or OZn as
possible candidates. Of these three, we rule out VZn because, as we have seen
Fig. 3
Figure 5. Polarized absorption spectra measured at 10 K for a virgin hydrothermally grown
ZnO sample: top E || c, bottom E A c.
in the previous section, zinc vacancies tend to bind two hydrogen atoms
resulting in VZnH2 defects, which are not stable at annealing temperatures
above 600 oC.
Indirect indication that Oi or OZn might be responsible for H-I* is the
annealing temperature of this center. It is known that ZnO crystals, when
heated in air at temperatures above 1000 oC become nonstoichiometric in
composition through formation of excess zinc atoms resulting from the more
rapid loss of oxygen atoms.22 It means that of all native defects Oi and OZn
143
5.
CONCLUSIONS
ACKNOWLEDGMENTS
The authors greatly acknowledge C. G. Van de Walle, and B. Bech
Nielsen for helpful discussions, R. Helbig for supplying the ZnO samples,
and the Russian foundation for Basic Research (Grant No. 02-02-16030).
REFERENCES
1. D.C. Look, D.C. Reynolds, J.R. Sizelove, R.L. Jones, C.W. Litton, G. Cantwell, and W.C.
Harsch, Solid State Comm., 105, 399 (1998).
2. M. Joseph, H. Tabata, and T. Kawai, Jpn. J. Appl. Phys., 38, L1205 (1999).
144
E. V. Lavrov et al.
3. D. C. Look, J. W. Hemsky, nad J. R. Sizelove, Phys. Rev. Lett., 82, 2552 (1999) and
references therein.
4. C. G. Van de Walle, Phys. Rev. Lett., 85, 1012 (2000).
5. S.F.J. Cox, E.A. Davis, S.P. Cottrell, P.J.C. King, J.S. Lord, J.M. Gil, H.V. Alberto, R.C.
Vilo, J. Piroto Duarte, N. Ayres de Campos, A. Weidinger, R.L. Lichti, and S.J.C. Irvine,
Phys. Rev. Lett., 86, 2601 (2001).
6. D.M. Hofmann, A. Hofstaetter, F. Leitner, H. Zhou, F. Henecker, B.K. Meyer, S.B.
Orlinskii, J. Schmidt, and P.G. Baranov, Phys. Rev. Lett., 88, 45504 (2002).
7. F. G. Grtner and E. Mollwo, phys. status solidi B, 89, 381 (1978).
8. F. G. Grtner and E. Mollwo, phys. status solidi B, 90, 33 (1978).
9. M. D. McCluskey, S. J. Jokela, K. K. Zhuravlev, P. J. Simpson, and K. G. Lynn, J. Appl.
Phys., 81, 3807 (2002).
10. N. H. Nickel and K. Fleischer, Phys. Rev. Lett., 90, 197402 (2003).
11. E. V. Lavrov, J. Weber, F. Brrnert, C. G. Van de Walle, and R. Helbig, Phys. Rev. B, 66,
165205 (2002).
12. E. V. Lavrov, Physica B, 340-342, 195 (2003).
13. G. Mller and R. Helbig, J. Phys. Chem. Solids, 32, 1971 (1971).
14. R. Helbig, J. Cryst. Growth, 15, 25 (1972).
15. D. G. Thomas and J. J. Lander, J. Chem. Phys., 25, 1136 (1956).
16. E. Mollwo, Z. Phys., 138, 478 (1954).
17.CRC Handbook of Chemistry and Physics, 82ndd edition, David R. Lide (ed.), CRC Press
LLC, 2001-2002, p. 9-77.
18. T. Sekiguchi, N. Ohashi, and Y. Terada, Jpn. J. Appl. Phys., 36, L289 (1997).
19. A. F. Kohan, G. Ceder, D. Morgan, and C. G. Van de Walle, Phys. Rev. B, 61, 15019
(2000).
20. U. Kaufmann, P. Koidl, and O. F. Schirmer, J. Phys. C, 6, 310 (1973).
21. H. J. Schulz and M. Thiede, Phys. Rev. B, 35, 18 (1987).
22. H. E. Brown, Zinc Oxide: Properties and Applications, New York, Int. Lead Zinc
Research Organization, 1976.
Chapter 13
INFLUENCE OF THE HYDROGEN
CONCENTRATION ON H BONDING IN ZINC
OXIDE
N. H. Nickel
Hahn-Meitner-Institut Berlin, Kekulstr. 5, 12489 Berlin, Germany
Abstract:
Despite of the recent developments the properties of hydrogen in ZnO are not
fully understood. While in some ZnO single crystals the amount of H detected
by gas effusion experiments is similar in concentration to the number of free
electrons, a large number of samples exhibits much higher H concentrations.
Raman backscattering experiments reveal that a significant amount of H is
accommodated at sites that do not give rise to enhanced conductivity; a large
amount of H is bound to carbon and nitrogen impurities forming CHX and
NH complexes. Depending on the method of the sample preparation H
concentrations as high as 3.0u1021 cm-3 have been observed. To gain further
insight into the properties of H in ZnO data obtained from H effusion
measurements of single crystal and sputter deposited ZnO have been analyzed
to deduce the hydrogen binding energy as a function of the H chemical
potential. In samples with a low H content six peaks were observed in the H
density-of-states distribution. With increasing H concentration the peaks
broaden. In addition, an increase of the H concentration results in a
pronounced increase of the average H binding energy. This observation
indicates that the properties of H expand well beyond the formation of donors
and the neutralization of impurities and deep defects.
Key words:
1.
INTRODUCTION
146
N. H. Nickel
2.
147
EXPERIMENTAL
Two sets of samples were used in this investigation. The first set consists
of commercially available ZnO single crystals with a (0001) orientation and
an O face. These crystals were fabricated by Eagle-Picher. The second set
of samples were Al-doped ZnO thin-films. These specimens were grown on
quartz substrates by dc magnetron sputtering at a power density of 22.6
W/cm2 and a substrate temperature of 150 C. The ZnO films had a thickness
ranging from 0.11 m to 3.83 m. From a structural point of view the
sputtered ZnO samples were composed of columnar grains with an average
grain diameter of about 0.16 m. The length of the grains is only limited by
the thickness of the film. The samples were characterized with Raman
backscattering experiments in a conventional macro-Raman setup using the
488 nm laser line of an Ar+ ion laser. Information on hydrogen content and H
bonding was obtained from H effusion measurements. For this purpose the
samples were placed in an ultrahigh vacuum system and heated to a
temperature of up to 950 C with a constant heating rate of 20 K/min. During
the temperature ramp the flux of molecular hydrogen was measured with a
quadrupole mass spectrometer. Prior to each measurement background
spectra were taken and subsequently subtracted from the sample data. To
obtain absolute values for the H2 flux the effusion spectra were calibrated by
measuring the know neon flux through a capillary.
3.
148
N. H. Nickel
(a)
(b)
after H effusion
2800
2900
3000
3100
3200
-1
Z (cm )
Figure 1. Raman backscattering spectra of single crystal ZnO before (a) and after hydrogen
effusion (b). The measurements were performed using the 488 nm line of an Ar+ laser and a
laser power of 190 mW. Background data were subtracted.
149
13
10
-2 -1
H2 flux (cm s )
12
10
11
10
crystalline ZnO
10
10
10
200
400
600
800
1000
T (C)
Figure 2. Molecular hydrogen flux of single crystal ZnO as a function of temperature. The
sample was heated with a heating rate of 20 K/min.
150
4.
N. H. Nickel
10
d (m)
0.11
13
2.06
-2 -1
H2 flux (cm s )
10
3.83
12
10
11
10
sputtered ZnO:Al
10
10
200
400
600
800
T (C)
Figure 3. Molecular hydrogen flux vs temperature for sputter deposited ZnO thin films. The
experiments were performed with a heating rate of 20 K/min. For clarity only three effusion
spectra are plotted.
151
Table 1. Hydrogen concentration in single crystal and sputter deposited ZnO. The asterisk
marks single crystal ZnO.
H concentration (cm-3)
3.31021
7.91020
1.01020
1.01020
5.21016
Thickness (m)
0.11
0.5
2.06
3.83
500 (*)
CH
N
H
dE
(1)
f
wCH
.
wP H
(2)
152
N. H. Nickel
700
high T peak
T (C)
600
500
400
low T peak
F
kT ln
F0
(3)
-3
0.59
1.35
(b)
0.92
1
0
0.0
-1
-3
single crystal
y
ZnO
1.55
0.5
1.0
E H (eV)
10
1.5
2.0
(c)
153
0.9
1.1
1
1.3
5
0
4
-1
0
2
1.03
1.18
1.4
20
(a)
21
20
NH (10 eV cm ) NH (10 eV cm )
-1
-3
16
-1
-3
NH (10 eV cm ) NH (10 eV cm )
0
0.96
(d)
1.44
2
0
0.0
0.5
1.0
1.5
2.0
E H (eV)
Figure 5. Hydrogen density-of-states distribution in single crystal ZnO (a) and Al doped ZnO
thin films (b) (d). Zero energy corresponds to the H transport site.
For sputter deposited ZnO the H content and thus NH increases up to 81021
cm-3eV-1 (see Fig. 5d). The data plotted in Fig. 5(b) show an important
result. The H effusion spectra were obtained on 2 different samples with
different thickness and the analysis yielded identical density-of-states
distributions indicating that the method to derive the H DOS is reproducible
and reliable.
Further analyzing the data in Fig. 5 reveals a striking observation. While
in single crystal ZnO about 45 % of all H atoms are accommodated with
binding energies of less than 1.0 eV with respect to the H transport site an
increase of the H concentration is accompanied with an increase of NH in the
energy range 1.0 eV < E H < 1.8 eV. With increasing H content the
amount of H atoms accommodated with binding energies of less than 1.0 eV
decreases to about 25 %. This suggests that H creates larger clusters
presumable in a similar way as it has been observed in other
semiconductors.17
5.
SUMMARY
154
N. H. Nickel
ACKNOWLEDGEMENT
The author is grateful to C. Klimm for the preparation of the ZnO thin
films and to K. Brendel for stimulating discussions and his constant interest
in the topic.
REFERENCES:
1. C. G. Van de Walle, Phys. Rev. Lett. 85, 1012 (2000).
2. S. F. Cox, E. A. Davis, S. P. Cottrell, P. J. C. King, J. S. Lord, J. M. Gil, H. V. Alberto, R.
C. Vilo, J. Pironto Duarte, N. Ares de Campos, A. Weidinger, R. L. Lichti, and S. J. C.
Irving, Phys. Rev. Lett. 86, 2601 (2001).
3. D. M. Hofmann, A. Hofstaetter, F. Leiter, H. Zhou, F. Henecker, B. K. Meyer, S. B.
Orlinskii, J. Schmidt, and P. G. Baranov, Phys. Rev. Lett. 88, 045504 (2002).
4. B. K. Meyer, H. Alves. D. M. Hofmann, W. Kriegseis, D. Forster, F. Bertram, J. Christen,
A. Hoffmann, M. Straburg, M. Dworzak, U. Haboeck, and A. V. Rodina, phys. stat. sol
(b) 241, 231 (2004), and references therein.
5. N. Y. Garces, N. C. Giles, L. E. Halliburton, G. Cantwell, D. B. Eason, D. C. Reynolds,
and D. C. Look, Appl. Phys. Lett. 80, 1334 (2002).
6. D. C. Look, J. W. Hemsky, and J. R. Sizelove, Phys. Rev. Lett. 82, 2552 (1999).
7. M. D. McCluskey, S. J. Jokela, K. K. Zhuravlev, P. J. Simpson, and K. G. Lynn, Appl.
Phys. Lett. 81, 3807 (2002).
8. E. V. Lavrov, J. Weber, F. Brrnert, C. G. Van de Walle, R. Helbig, Phys. Rev. B 66,
165205 (2002).
9. S.-Y. Lin, J. Appl. Phys. 80, 1399 (1996), and references therein.
10. D. M. Joseph, R. Balagopal, R. F. Hicks, L. P. Sadwick, and K. L. Wang, Appl. Phys. Lett.
53, 2203 (1988).
11. M. O. Manasreh, J. M. Baranowski, K. Pakula, H. X. Jiang, and J. Lin, Appl. Phys. Lett.
75, 659 (1999).
12. C. G. Van de Walle, presentation at the MRS Spring Meeting in San Francisco 2004.
13. J. W. Bozzelli and R. B. Barat, Plasma Chem. Plasma Proc. 8, 293 (1983).
155
Chapter 14
VALENCE BAND ORDERING AND MAGNETOOPTICAL PROPERTIES OF FREE AND BOUND
EXCITONS IN ZnO
A.V. Rodina1, M. Strassburg2, M. Dworzak2, U. Haboeck2, A. Hoffmann2,
H. R. Alves3, A. Zeuner3, D. M. Hofmann3, and B. K. Meyer3
1
Abstract:
The Zeeman splitting of free and bound excitons is analysed theoretically for
both * and *symmetries of the upper valence band. The Zeeman splitting of
the ** holes in magnetic field parallel to the hexagonal c axes of the
crystal is found to be smaller (larger) than the Zeeman splitting of the
electrons. The magnetic-field dependence of the bound exciton optical
transition energies measured in the Faraday geometry with right and left
circular polarized light can be well described with * as well as with *holes.
However, the fitting of the experimental dependence of the Zeeman splitting
on the angle between magnetic field and c axis proves the * symmetry of the
holes involved in all low temperature transitions and thus of the upper valence
band in bulk ZnO. The determined parameters are in excellent agreement with
known experimental data on the fine structure of the free exciton ground state
in ZnO. The hole g-factors determined for excitons bound to an ionized and
neutral impurity centers are close to each other and to the g-factor calculated
for the *7 hole in the free exciton ground state. This points to a donor character
of the neutral impurity centers. This conclusion is confirmed by temperature
dependent magneto-transmission and magneto-luminescence measurements.
Key words:
1.
INTRODUCTION
Bulk ZnO is a direct band gap wurtzite semiconductor with the valence
band maximum split by the crystal field and spin-orbit interaction. The
159
N.H. Nickel and E. Terukov (eds.),
Zinc Oxide A Material for Micro- and Optoelectronic Applications, 159170.
2005 Springer. Printed in the Netherlands.
160
A. V. Rodina et al.
symmetry of the upper valence subband (A-subband) has been the subject of
controversy (9 or 7 character) for more than 40 years.1-15 Based on the
polarization properties of the free exciton transitions, the inverted ordering of
the valence subbands as shown in Figure 1(a) was first suggested by
Thomas1 and supported later in many articles.3-7,11-18,21 However, different
argumentation led some researches2,8-10,19,20 to the conclusion that the
ordering of valence subbands in ZnO is the same as in other wurtzite
semiconductors like GaN, CdSe, CdS (see Figure 1(b)).
(a)
(b)
E
*7
*7
Egap=3.4376 eV
T=4.2 K
*7
A
A
*9
AB = 4.9 meV
B
*7
43 7 meV
V
EBC = 43.7
B
*7
C
Figure 1. Energy band structure and the symmetry of the free exciton ground state in wurtzitetype semiconductors for (a) ZnO and (b) GaN, CdSe, CdS.
2.
161
THEORY
P B Bg
B exc
P B B(
//
h
),
(1)
162
A. V. Rodina et al.
(a) excitons with *7 holes: *7 x *7 =*
* 5+ * 1+ * 2
B
B = 0
B ||c
*5
* 1, * 2
E||c
V V-
E||c
*1
B = 0
B ||c
*5
*6
E||c
V V-
E||c
*1
Figure 2. Level schemes and selection rules for the ionized bound exciton transitions.
//
h
|2 c o s 2
163
(2)
|2 s in 2 4
Therefore, for the * holes the Zeeman splitting in any magnetic field
remains always smaller than the Zeeman splitting of the electrons. In
contrary, for the * holes some critical angle 4crr should exist so that
gh(4crr)=ge, and gh(4)>ge for 0q d 4 < 4cr and
d gh(4 )< ge for 90q t 4 >4cr.
Examination of the selection rules for optical transitions in the B // c // k
configuration shows that the upper-energy transition is expected to be active
for the right-circular polarized light V+ for all cases presented in Fig. 3 cases.
(a) (D
D0,X((*7))
D0
g |0
A
h
| g h|||| | 1.4
(c) (A
A0,X)
, )
h
B
+
+1/2
B ||c
-1/2
1/22
A0(*9)
+1/22
g e | 1.95
| gh|||| | 0.75
ghA |0
+1/22
d)) ((A
A0,,X)
X))
((d)
-1/2
1/2
-3/2
3/22
V-
ge | 1.95
-1/2
1/22
g e | 1.95
| gh|||| | 3.0
V
+1/22
A0((*
*7)
D0 +3/22
ge | 1.95
g hA
+1/2
V
e
-1/2
1/22
D0,X((*9))
(b) (D
-1/2
1/22
+3/22
ghA 0
-1/2
+1/2
||
h
| g || 3.3
B ||c
-3/2
3/22
Figure 3. Level schemes and selection rules for neutral bound exciton transitions.
3.
We used nominally undoped (as grown from Eagle-Picher, see for details
Ref. 15) and nitrogen doped bulk ZnO crystals. A nitrogen doped ZnO:N
sample was prepared from the as grown one by ion implantation (2 MeV,1013
cm2) followed by a thermal annealing process at 900 C for 15 min (see also
Ref. 15).
164
A. V. Rodina et al.
4.
I4
(b) ZnO:N
I2
k//c
(a) ZnO
as grown
3,356
3,360
165
k//c
3,364
,
Energy (eV)
3,368
Figure 4. Photoluminescence spectra of undoped as grown ZnO (a) and nitrogen doped
ZnO:N (b) at 4.2 K.
166
A. V. Rodina et al.
(a) B ||| c
V-
PL Intensity
y [norm.
[
u.]
u]
(c) B A c
( B=5T
(b)
I3/I2
90
5T
80
4T
60
3T
40
2T
4T
3T
2T
20
0T
I2
3,3675
3,3680
Energy [eV]
0T
0
3,366
3,367
3,368
Energy [eV]
3,366
3,367
3,368
Energy [eV]
Figure 5. Photoluminescence spectra of ZnO:N in the energy region of the I2/I3 lines at 4:2 K
for different magnetic fields in Faraday configuration (case (a): B // k // c, E A c) and Voigt
configuration (case (c): B A c, k // c, E A c), and for different angles between the c axis and
the direction of the magnetic field B at 5 T (case (b)).
The magnetic field and angular dependencies of the emission lines I2 and
I3 allow us to determine the symmetry and the g values of the valence band
hole involved in the corresponding bound excitons.14,15 The magnetic field B
/// c linearly splits the I2 transition (allowed in E A c configuration) as well as
the I3 transition (allowed in E // c configuration). The upper-energy
component of the I2 line is active for the right-circular polarized light (see
dashed line in Fig. 5 (a)) as predicted in Fig. 2. The linear splitting of the I2
line in B // c can be described with the positive exciton effective g- factor
gexc = 0.71 that is close to the calculated values in Fig.2 (a) and (b). However,
the experimental value gexc = 2.66 determined for the I3 line is in better
agreement with the expected value for the * holes (case (a) in Fig. 2) and
differs significantly from the respective value in case (b). Furthermore, from
the fitting of the angular dependence of the transition energies14,15 we
obtained gh// = -1.24 in the excellent agreement with the calculated value for
the * hole in the exciton ground state. This demonstrates that the hole in the
lowest exciton state and thus the upper A valence subband in ZnO also have
the * symmetry.
The Zeeman behavior of the emission lines from I4 to I9 is different from
that of the I2/II3 lines. Their linear splitting in B A c indicates transitions
originating from excitons bound to neutral impurities. We discuss the
emission lines I9 and I4 in detail. These lines representt the dominant
(c) B A c
(b) B = 5 T
V+
PL Intensity
y [norm.
[
u.]
V-
167
5T
90
5T
80
4T
60
3T
4T
3T
40
2T
2T
20
0T
I9
3,3565
0T
3,3570
I9
3,3565
Energy [eV]
3,3570
Energy [eV]
3,3565
3,3570
Energy [eV]
Figure 6. Photoluminescence spectra of ZnO:N in the energy region of the I9 line at 4:2 K.
Cases (a), (b) and (c) are the same as in Fig. 5.
E ( B, )
E0 1/ 2
( ))) ,
(3)
where E0 is the zero field transition energy and gh(T) is given by Eq. (2). For
the *9 holes the ghA =0 by the symmetry and the Eq. (3) leads to the
linear dependence off E(B,4) on cos((4) (see dashed lines in Fig. 7 (b)). For
the *7 holes in ZnO the value of ghA | 0 and the nonlinear behavior of
E(B,4) on cos((4) can not be detected as well (see solid lines in Fig. 7 (b)).
However, the symmetry of the holes can be unambiguously established by
the analysis of the angular dependence in Fig. 7 (b). Indeed, no critical
168
A. V. Rodina et al.
(a) k||B || c
(b) B = 5 T
(c) B A c, k || c
Energy (eV)
3,3635
I4
3,3630
3,3625
V + ZnO
V - ZnO
3,3620
B (Tesla)
EA c ZnO
1,0
0,5
Cos(T)
0,0
B (Tesla)
Figure 7. Zeeman splitting of the bound exciton line I4. Case (a): B // k // c, E A c. Case (c):
B A c, k // c, E A c. Case (b): the angular dependence for B = 5 T. Symbols are experimental
data measured with circular polarized light (V + and V -) and with unpolarized light (E A c) as
denoted in the Figure, and lines are fits. For the fits the hole is assumed to be of *7
symmetry with | gh|||< ge (solid lines) and of *9 symmetry with | gh|||> ge (dashed lines).
169
5.
CONCLUSION
ACKNOWLEDGEMENTS
The work of A. V. Rodina was carried out during the stay at the Institute
for Solid State Physics, TU Berlin, and supported by the Deutsche
Forschungsgemeinschaft (DFG).
REFERENCES
1.
2.
3.
4.
5.
6.
7.
8.
170
A. V. Rodina et al.
Chapter 15
FUNDAMENTAL OPTICAL SPECTRA AND
ELECTRONIC STRUCTURE OF ZnO CRYSTALS
V. Val. Sobolev and V. V. Sobolev
Udmurt State University, 426034 Izhevsk, Russia
Abstract:
The analysis of ZnO exciton spectra is given. In the band and exciton theory,
necessity of taking into account hybridization of O 2p and Zn 3d states is
confirmed. Full sets of ZnO exciton spectra are obtained for 1.6, 4.2, and 90 K
for the first time. The main parameters of the three exciton series are
determined. Five variants of full sets of the fundamental polarized optical
functions of ZnO crystal were determined in the energy range 026 eV for
polarizations E A c and E || c. Three variants of dielectric permittivity and
volume electron energy loss polarized spectra of ZnO crystal were
decomposed into the transverse and longitudinal transition band components.
The main features and parameters of the five variants of the optical functions
spectra were established. A theoretical analysis of the obtained components
was performed, based on the known theoretical band schemes.
Key words:
1.
INTRODUCTION
The polarization of the exciton spectra of ZnO does not follow the
classical model. All the three exciton series are strongly polarized and can be
observed only with E A c (A, B) and E || c (C). This has led to contradictory
position schemes of the two highest valence bands (*9, *7) and values of 'so,
'cr. There are also certain discrepancies in the results of different papers on
the nature of the lines of complex discrete reflectivity and absorption spectra
of ZnO at 4.2 K, as well as on the bands and excitons parameters.113
Polarized reflectivity spectra of ZnO have long been measured in the
wide energy region of fundamental absorption 026 eV, and some other
optical functions have been calculated from them.4,7,14,15 However, full sets
171
N.H. Nickel and E. Terukov (eds.),
Zinc Oxide A Material for Micro- and Optoelectronic Applications, 171182.
2005 Springer. Printed in the Netherlands.
172
2.
CALCULATION TECHNIQUES
3.
173
the positions of the exciton lines. This agrees well with the simple classical
exciton model and contradicts the polariton model.
The problem has long been discussed of a connection between the
differences in the exciton spectra and the triplet band structure of ZnO and
other hexagonal crystals.111 In all these papers, the simplest model is
accepted on the assumption of invariability of the symmetry of the three
bands: in some papers, inversion of two upper bands is assumed,1,12 i.e.: the
highest band is *7 instead of *9, while other papers propose to retain the
conventional interleaving of bands *9, *7, and *7.2,3,10 This problem is
theoretically considered.9
It was shown in the quasicubic approximation of the group theory that the
relative energies of three upper valence bands with energies E1, E2, E3
are determined by the values of 'so and 'cr:
E2
A2
2 B 3, E3
A2 2 B 3
(1)
where 2 A ' so ' cr ; B ' so ' cr ; E1 = 0; 'so, 'crr << Eg. The values of 'so
for ZnO obtained with Eq. (1) vary in their absolute values and even in their
sign: 0.009,1 0.015,12 0.01610 eV. The calculations of 'so and 'crr are
known that use the classical model of band interleaving (*9, *7, *7) and
inversion (*7, *9, *7). The values of 'crr proved to be almost identical in all
the papers, regardless of the band model: 'crr = 0.041,1 0.035512 eV for *7,
*9, *7; 0.041,3 0.04310 eV for *9, *7, *7. However, the values of the second
parameter ('so) turned out to be mostly negative: 'so 0.009,1 0.01512 eV for
*7, *9, *7; 0.005,3 0.01610 eV for *9, *7, *7. Large divergence of the
calculation results for 'so seem to be due to too small a distance between two
( , B) exciton series and larger relative errors of ((BA
(A
)/A
/ .
Two experimental results have been known, which show a possibility of
considerable change of the parameters of all the three valence bands of ZnO.
First, polarization of all the three exciton series of ZnO conceptually differs
from their polarization in other A2B6 crystals that is in good agreement with
the conventional band model. This fact undoubtedly shows that in ZnO, not
only the symmetry of the second band changes, but the symmetry of the
third band also. Second, according to the photoemission spectra (PES), the
valence band of zinc d-states (Zn 3d) is located below the top of the valence
band of ZnO by ~ 8 eV, while below its bottom only by ~ 3 eV.19 Therefore,
Zn 3d states can intermix with O 2p states and change the symmetry of the
upper valence bands.
According to the theoretical calculations,20 O 2p states form ~ 67 74%
of the highest valence band, while Zn 3d states form the remaining almost
174
' so ( p O) a
so
(d Zn)
(2)
According to the estimations, 'so(d Zn) | 0.37 eV, a = 0.07. Oxides BeO
and MgO, lacking d-bands, have 'so(p O) | 0.034 eV. Hence, 'so(ZnO) |
(0.0340.026) eV < 0.01 eV.
Lambrecht and his colleagues have theoretically shown11 that the
energies 'so and 'crr depend very strongly and almost linearly on the location
of the d-levels of zinc ((E
Ed). According to the authors, the most optimal
variant is that with Ed | 6.25 eV together with the band inversion model.
However, the experimental results give Ed | 8.5 r 0.4 eV.19 Unfortunately, it
is in the energy region E > 6.7 eV that the calculations of this paper are not
reliable. Thus, both the theoretical and experimental results indicate that zinc
d-states make a substantial contribution to the formation of the highest
valence band of ZnO. It is the hybridization of O 2p and Zn 3d states that
strongly changes the symmetry of the three valence bands and, consequently,
the polarization of three exciton series A, B, C, as well as the splitting
energies of the three bands.
4.
175
Table 1. Energy values E (eV) of H2 and Im H1 maxima, bands areas S, and oscillator
strengths f u 105 of ZnO excitons, calculated from the R(E) spectra of papers Lagois, et al.4
(L), Park, et al.2 (P), and Klucker, et al.5 (K)
N
B1
C1
A1
H2
ImH1
H2
ImH1
H2
ImH1
L,
E
3.3766
3.3792
3.3830
3.3948
3.421 3.439
0.16
0.0005
0.20
0.008
0.32
0.007
1.6 K S
f 105 (107)
fu
7.5
2.4
9.4
38
15
33
P,
E
3.3746
3.3772
3.3796
3.3904
3.419 3.433
0.20
0.0006
0.16
0.006
0.50
0.006
4.2 K S
f 105 (107)
fu
9.4
2.8
7.5
2.8
23.5
28
K,
E
3.349
3.38
3.349
3.38
3.387 3.424
90 K
S
0.50
0.008
0.50
0.008
0.50
0.008
f 105 (107)
fu
24
38
24
38
24
38
N
A2
B2
C2
ImH1
H2
ImH1
H2
ImH1
H2
L,
E
3.43
3.43
3.44
3.43
3.47
3.47
1.6 K S
0.02
0.0005
0.03
0.0005 0.05
0.0007
f 105 (107)
fu
0.94
2.35
14
2.4
2.35
3
P,
E
3.42
3.42
3.43
3.43
3.46
3.47
4.2 K S
0.04
0.04
0.02
0.001
0.1
0.003
f 105 (107)
fu
1.9
1.9
0.94
4.7
4.7
1.4
K,
E
3.41
3.41
3.41
3.43
3.453
3.49
S
0.10
0.10
(0.10)
0.008
0.05
0.001
90 K
f 105 (107)
fu
4.4
4.4
4.7
38
2.4
4.7
5.
176
0.25
0.02
0.15
0.00
2.5
3
1 4
0.01
-1
P (10 cm )
0.20
2
5
3.0
(a)
3.5
4.0
0.10
0.05
0.00
10
15
E (eV)
20
0.25
(b)
2
0.02
-1
P (10 cm )
0.20
0.15
0.00
2.5
25
5
3.0
0.10
3.5
4.0
0.05
0.00
10
15
E (eV)
20
25
Figure 1. The P (a, b) and E2H2 (c, d) spectra of ZnO crystal for EAc (a, c) and E||c (b, d) in
the range of 025 eV (in the inset in the range of 2.54.0 eV), calculated from the
experimental R spectra (1, 2, 3) and from the theoretical H2 and H1 spectra (4, 5)
1000
2
30
H (E )
800
20
600
(c)
3
4
10
0
2.5
5
3.0
177
3.5
4.0
400
200
0
40
800
20
0
2.5
600
15
E (eV)
4
5
3.0
3.5
4.0
200
25
(d)
400
20
H2 (E )
1000
10
10
15
E (eV)
20
25
Figure 1. Continued: The P (a, b) and E2H2 (c, d) spectra of ZnO crystal for EAc (a, c) and
E||c (b, d) in the range of 025 eV (in the inset in the range of 2.54.0 eV), calculated from
the experimental R spectra (1, 2, 3) and from the theoretical H2 and H1 spectra (4, 5)
178
The maxima of the theoretical R(E) spectra of papers ref. 24 and 25 differ
one from the other by ~ 0.21 eV, while the differences between these and
the maxima from papers5,14,15 come up to 1.5 eV, with the intensity
distribution of the maxima remaining mostly similar. A noticeable general
similarity was also discovered as a result of direct comparison of the
remaining experimental and calculated, and theoretical spectra of the optical
functions (Fig. 1). However, the detailed comparison of the theoretical and
experimental spectra of the optical functions showed, naturally, substantial
differences in the intensities and locations of the maxima.
In papers 3336, the photoemission spectra of ZnO crystal were
experimentally studied. According to the presented results Zn 3d electron
energy levels of ZnO ((E
Ed) are lower than the maximum of the valence bands
for ~ 8.5 eV,19 7.5 eV,34 and 7.5 eV,21 while the total width of the highest
valence band ('E
'Ev) is equal to ~ 5.0 eV,34 8.0 eV,19 and 6.0 eV.21 From the
theoretical band calculations, 'E
'Ev is equal to ~ 6 eV24 and 4.2 eV,25 while Ed
25
| 4.2 eV; zinc d-levels being ignored in paper ref. 21. This means that both
important band parameters of ZnO ('E
'Ev and Ed) are undervalued in paper
ref. 24 by ~ 2 eV and 4 eV for 'E
'Ev and Ed, respectively. This must have led,
in the first place, to underestimation of the energy values of the theoretical
transitions from the lower sub bands of the highest valence band and
alteration of the transitions intensities.
6.
Spectra of H2 and ImH-1 of the three variants of optical functions sets ((A,
B, C
C) were deconvoluted by us into components. The main parameters of the
obtained components are given in Tables 1, 2 ((A), 3, 4 ((B), and 5 (C),
C as well
as in Fig. 1.
In Table 2, the calculated values of energy maxima Ei and band areas Si
of transverse transitions were compared for the three variants of calculations
of optical function sets. These variants differ comparatively little in energy
Ei. The principal inconsistencies are observed in the values of Si. This is
caused by the differences in the experimental R(E) spectra, i.e. by using of
ZnO samples of different quality and by different registration techniques for
polarized reflectivity spectra. The determined components of H2 and ImH-1
spectra are caused by direct interband transitions or metastable excitons,
except for the most long-wavelength of them, which are associated with free
excitons. The theoretical band calculations of ZnO20,2325 strongly differ in
the bands dispersion and positions. This makes it difficult to propose a
179
0.36
1.59
3.4
0.26
2.56
3
3.58
0.62
1.01
3.7
3.7
0.84
1.9
4
4.34
4.42
1.11
0.9
3.8, *
5
4.7
4.9
1.11
1.01
5.0, *
6
6.2
6.1
2.67
1.21
7.5, '
7
6.7
6.7
2.36
1.2
8
8.4
8
0.57
0.21
9.5, ', 6
9
8.8
9.3
9
9.3
4.22
1.39
1.05
0.2
10, U
10 11.1
11.4
11.2
11.4
0.67
1.66
0.12
0.35
11.0, S
11 12.5
12.7
12.9
12.9
6.85
2.2
0.84
0.32
12.6, U
12 13.7
13.7
13.7
13.8
1.78
6.83
0.17
0.77
14.0, S
13 14.7
14.7
15
15
3.8
3.51
0.31
0.3
14.6, U
14 16
16
15.5
1.87
1.69
0.14
0.13
16.7, '
15'
16.4
16.3
0.07
15 17.2
17.3
17.5
17.5
0.77
0.45
0.06
0.03
17.8, U
16 18.7
18.4
18.7
18.3
0.48
0.89
0.03
0.06
16
19.1
19.1
17'
20.1
19.7
20.3, U
17 20.7
21.5
20.9
20.9
1.47
1.89
0.11
0.28
18 22.9
23.4
23.4
23.7
0.34
1.08
0.03
0.08
22.0, U
19 25.1
25.1
2.46
0.19
20
26.8
26.8
2.31
0.17
7.
CONCLUSIONS
180
3.10
0.36
0.44
0.15
0.35
2
3.18
3.19
3.20
0.08
0.22
0.09
2
3.46
3.48
3.40
3.35
0.68
0.05
0.26
0.95
3
3.58
0.62
3.69
3.70
3.70
0.91
0.84
0.05
4
4.34
4.10
4.25
4.15
1.11
0.75
0.14
2.16
5
4.60
4.70
4.80
0.63
0.13
6
6.20
5.38
5.13
5.25
2.67
0.39
0.96
1.19
7
7.20
7.3
6.70
7.40
1.10
0.22
2.36
1.91
8
8.2
8.40
0.06
0.57
9
8.80
8.8
9.00
9.1
4.22
0.53
1.39
0.65
10 11.1
11.4
11.2
11.1
0.67
0.21
1.66
1.08
11 12.5
12.8
12.3
12.9
12.0
6.85
2.74
0.35
2.20
0.67
12 13.7
13.7
12.7
13.7
13.7
1.78
0.38
1.68
6.83
6.77
13 14.7
15.1
14.0
15.0
14.8
3.80
7.58
1.20
3.57
0.38
14 16.0
15.5
15.5
16.0
1.87
4.32
1.69
15 17.2
17.3
17.6
17.5
17.9
0.77
3.75
3.88
0.45
1.46
16 18.7
19.5
18.8
18.7
0.48
1.42
0.16
0.89
17 20.7
21.9
20.4
20.9
21.9
1.47
2.60
2.26
1.89
5.05
18 22.9
23.3
23.4
23.3
0.34
0.24
1.08
19 25.1
24.7
24.3
24,7
2.46
1.09
2.50
0.05
20
26.8
2.31
181
ACKNOWLEDGEMENTS
The paper was supported by the Fundamental Natural Sciences Center,
State University of Saint-Petersburg, Russia.
REFERENCES
1.
2.
3.
4.
182
5. R. Klucker, H. Nelkowski, Y. S. Park, M. Skibowski, and T. S. Wagner, phys. stat. sol. (b)
45, 265 (1971).
6. V. V. Sobolev, V. I. Donetskich, and E. F. Zagainov, Fizika i Tekhnika Poluprovodnikov
12, 1089 (1978).
7. V. Sobolev, (Stiintza, Kishinev, 1980).
8. P. Kuzmina and V. A. Nikitenko, Poluchenie i opticheskie svoystva, (Nauka, Moscow,
1984).
9. J. Hopfield, J. Phys. Chem. Sol. 15, 97 (1960).
10. Reynolds, D. C. Look, B. Jogai, and C. W. Litton, Phys. Rev. B 60, 2340 (1999).
11. L. Lambrecht, A. V. Rodina, and B. Segal, Phys. Rev. B 65, 75207 (2002).
12. H. Yoshikawa and S. Adachi, Jpn. J. Appl. Phys. 36, 6237 (1997).
13. G. E. Jellison and L. A. Boatner, Phys. Rev. B 58, 3586 (1998).
14. J. L. Freeouf, Phys. Rev. 7, 3810 (1973).
15. R. L. Hengehold, R. J. Almassy, and F. L. Pedrotti, Phys. Rev. B 1, 4784 (1970).
16. V. V. Sobolev and V. V. Nemoshkalenko, Elektronnaya struktura poluprovodnikov
(Naukova dumka, Kiev, 1988).
17. V. V. Sobolev, A. I. Kalugin, V. Val. Sobolev, and V. I. Kormilets, Wide Bandgap
Materials 8, 87 (2001).
18. V. Val. Sobolev and V. V. Sobolev, Phys. Low Dim. Struct. 5/6, 65 (2003).
19. C. J. Vesely, R. L. Hengehold, and D. W. Langer, Phys. Rev. B 5, 2296 (1972).
20. D. Vogel, P. Kruger, and J. Pollmann, Phys. Rev. B 54, 5495 (1996).
21. S. A. Leontyev, S. V. Koshcheyev, V. G. Devyatov, A. E. Cherkashin, and E. P.
Mikheyeva, Zh. Strukturnoy Khimii 38, 867 (1997).
22. S. Bloom, and J. Ortenburger, phys. stat. sol. (b) 58, 561 (1973).
23. S. B. Zhang, S.-H. Wei, and A. Zunger, Phys. Rev. B 37, 8958 (1995).
24. J. R. Chelikowsky, Solid State Commun. 22, 351 (1977).
25. J. Xu and W. J. Ching, Phys. Rev. B 48, 4335 (1993).
Chapter 16
PHOTO-INDUCED LOCALIZED LATTICE
VIBRATIONS IN ZnO DOPED WITH 3D
TRANSITION METAL IMPURITIES
Alexey Kislov
Ural State Technical University, Mira Str. 19, 620002 Ekaterinburg, Russia
Abstract:
The purpose of this work is twofold: to simulate vibrational spectra of ZnO crystals
with nickel impurities in different charge states and to analyze vibrational
structures of electroabsorption spectra of ZnO:Ni based on calculations of
localized vibrations induced by the impurities.
Key words:
1.
INTRODUCTION
184
Alexey Kislov
Figure 1. Diagrams of energy levels of a 3d impurity [(a) and (b)]: transitions involving an
ionization of the impurity (D, A) and transitions to hydrogen-like states (d, a) are shown. Fig.
1(c) is an illustration of the transitions (d, a) in a space near a Ni impurity.
185
The formation of the impurity exciton leads to an appearance of photoinduced lattice vibrations. Sokolov1 has offered an experimental
electroabsorption (EA) method, which has a very high sensitivity for
detecting such vibrations. However, the analysis of experimental data
supposes a theoretical study of the vibrational spectrum of the defective crystal.
The purpose of our work was to simulate vibrational spectra of ZnO crystals
containing nickel impurities in different charge states with respect to the
lattice and also to analyze vibrational structures of EA spectra of ZnO:Ni based
on calculations of localized vibrations induced by the impurities.
2.
METHODS OF CALCULATION
The lattice relaxation near the impurities was calculated on the basis of
the molecular static method2. The crystal with the impurity was dividing into
regions: an inner region 1 and an outer region 2+3 (see Fig. 2). The energy of
the static lattice of this system can be written as2
&
Ests E
E(( ,[ )
( )
&
( , )
&
([ )
(1)
186
Alexey Kislov
&
where& E1( ) is the energy of the region 1, E3 ([ ) is the energy of region 3, and
&
E2 ( , ) is the interaction energy of regions 1 and &3. x is the coordinate
describing the lattice configuration of region 1 while [ is the displacement in
&
region 3 and is formally distinguished from x .
The ion displacements in region 1 can be determined by direct
&
minimization of the defect energy with respect to x , i.e.: by solving the
equation
i,
i j
wwE
E st
wx
wx
x(
i, j, ,
i j
w 2 E st
w xi w x j E
&
x(
u stjE
(2)
g iD ( )
Im
uiD
( ) .
(3)
iD
| (
|[
0))
] 1 | uiD
3N
s 1
| isD |2
.
Z i00 Z s2
2
(4)
u n 1
(D
)u n
bn 1u n 1
(5)
1
D 0
...
0
b1
0 0 ... 0
b2 0 ... 0
2
1 ... ... ... ... .
(6)
187
In this case the diagonal element of the Green's function are described by the
easily calculated continued fraction
Gu1 ( )
Z 2 a1
b1
(7)
b2
Z 2 a2
...
bn1
2
Z an bnt ( )
1
2 i, j
1
(r )
k P2
ij 2 i i
i
j i
(8)
1
[
2 i, j
1
( )]
k P2
ij
2i i i
( )
ij
j i
V Buck ( ij )
X X
XY
YX
YY
i j
i j
i j
i j
c
s
c
c
s
s
r
|r
r | |r
r | |r
rs |
ij
j
i
j
i
j
i
ij
exp( |
s
j
s
i
|/
ij
)
Ciij
|
s
j
s 6
i
(9)
(10)
188
3.
Alexey Kislov
890.0
1265.0
88760.0
0.300
0.300
0.149
Cij (eV/ )
(Ref. 6)
0.0
0.0
0.0
ki (eV/ 2)
116.4
51.8
Observed
3.25 (Ref. 6)
5.21 (Ref. 6)
42.2 (Ref. 7)
209.7 (Ref. 8)
121.1 (Ref. 9)
105.1 (Ref. 9)
210.9 (Ref. 8)
42.47 (Ref. 8)
44.29 (Ref. 8)
8.33 (Ref. 6)
8.84 (Ref. 6)
3.94 (Ref. 8)
4.00 (Ref. 8)
8.0 11.5 (Ref. 8)
17.7 (Ref. 10)
As an example, Fig. 3 plots the phonon dispersion curves for three highly
symmetric directions in the Brillouin zone of the perfect ZnO crystal.
Comparison of the theoretical and experimental frequencies8 shows good
agreement for the acoustic branches. The densities of phonon states of the
perfect ZnO crystal calculated by integrating over the Brillouin zone are
displayed in Fig. 4. Comparison of the results of our calculation and a calcu-
189
Figure 3. Phonon dispersion curves for ZnO crystals. The points denote experimental values.8
lation based on the valence force model11 shows clearly that the approximation
of pairwise central interactions is satisfactory. The difference in the position of the
phonon forbidden region and of some maxima is due to limitations associated
with the pairwise interaction approximation used in the calculations.
The formation of the nickel excitons results in the lattice distortion near
them and the induced lattice vibrations. Conditions of their occurrence are
defined by the charged impurity because the removed hydrogen-like type
carrier practically does not influence the deformation of the lattice near the
charged impurity. Thus, our analysis of the vibrational background of the
zero phonon line of the EA spectrum of the nickel exciton is based on results of
a simulation of the lattice dynamics of the ZnO crystal with Ni+1 or Ni+3 ions.
An important question is how to correctly describe the interaction between
the defect and the remaining atoms of the crystal. The parameters of the shortrange interaction potentials between the Ni impurities in various charge
states and its surrounding ions being taken as the corresponding parameters
190
Alexey Kislov
Figure 4. Total density of phonon states in ZnO crystals. Dashed line - result of the
calculation from Ref. 11; solid line - our calculation.
for Zn ions. Of course, we took into account the change of mass and charge
associated with the substitution of Ni+1 or Ni+3 for Zn.
Table 3 lists the displacements of nearest neighbors around the charged
nickel. The calculations show that the nearest neighbors of the Ni+1 impurity are
quite removed from it. For Ni+3 they shift toward the impurity.
The effect of charged nickel ions on the vibration spectrum of ZnO crystals
was studied by calculating the symmetrized LVDS in perfect and defective
crystals using a cluster of 1000 ions (region 1). For the calculation of the
Coulombic component of the diagonal elements of the force constant
matrixes of the ion-ion interactions by the Ewald method about 4000 ions
were considered for region 2.
We have calculated the symmetrized LVDS projected onto a region
containing two coordination spheres around a zinc ion or its nickel
replacement. Nickel ions in ZnO crystals are in the positions with the
symmetry point group C3v. The 36-dimensional displacement space splits
into the irreducible representations of the C3vv point group: 3(36u36)= 91 +
32 + 12.
191
Table 4. The defect formation energy and the displacements around the defect.
Defect
Ni+1
'Ed (eV)
18.4
Ni+3
-30.3
Displacement ()
0.251
-0.017
-0.240
0.027
Figure 5. LVDS projected onto displacements with A1-symmetry of the center ion (a) and the
O-2(z) ion (b).
192
Alexey Kislov
For ZnO:Ni the EA spectrum depicted in Fig. 6 (a) was measured at liquidhelium temperature T=4.2 K;12 therefore, its vibrational part is shifted into the
Stokes region of the spectrum and is associated with phonon creation.
Table 5. The frequencies of localized vibrations in ZnO with nickel impurities.
Defect
Type of mode
Type of symmetry
Frequency (THz)
Resonance
A1
3.0, 13.8
Gap
A1
10.0
Ni+1
Local
A1
16.8
Resonance
E
3.0, 13.8
Local
E
16.8
Ni+3
Resonance
A1
Gap
Local
Resonance
A1
A1
E
Gap
Local
E
E
4.
CONCLUSIONS
193
Figure 6. Spectral dependence of the EA second harmonic D2 for ZnO:Ni (a). Symmetrized
LVDS projected onto displacements with A1- and E-symmetries of the ions of two
coordination spheres around a Ni+1 impurity (b).
194
Alexey Kislov
REFERENCES
1.
2.
3.
4.
5.
6.
7.
8.
Chapter 17
ZnO WINDOW LAYERS FOR SOLAR CELLS
Walther Fuhs
Hahn-Meitner-Institut Berlin, Kekulstrae 5, D-12489 Berlin, Germany
Abstract:
ZnO layers are very frequently used as window layers in photovoltaic solar cell
structures. In this application high electrical conductivity should be
combined with a low value of the optical absorption constant in the visible
range. Thus an appropriate measure of the performance is the ratio r = / and
the optimization consists in realizing high by high carrier concentration and
mobility with a minimum optical loss due to free carrier absorption at long
wavelength. Such layers can be prepared by numerous techniques most
frequently by magnetron sputtering either from compound targets or reactively
from metallic targets. Optimized doped ZnO-layers have resistivities of 1.52u10-4 cm with Hall mobilities of up to 60 cm2/Vs, refractive indices of
approximately 2.0 and an average transmittance of 85 % in the visible range. It
appears that for the value of the physical limit as defined by ionized
impurity scattering has been reached for homogeneously doped thick layers.
This paper will discuss various aspects of the application of ZnO window
layers in photovoltaic solar cells.
Key words:
1.
INTRODUCTION
198
Walther Fuhs
concentration and mobility with a minimum optical loss due to free carrier
absorption at long wavelength. However, for a particular application this
figure of merit is not the only parameter that has to be considered. Further
aspects may become important such as physical, chemical and thermal
durability, stability, homogeneity, chemical etchability, workfunction,
method of preparation, toxicity, availability of the materials or simply cost
factors. So far tin oxide, SnO2, prepared by chemical vapor deposition
(CVD) is most commonly used in low-energy windows and indium tin oxide
(ITO) made by spray pyrolysis or sputter deposition is the preferred material
for flat panel displays.1 ZnO films are used in particular for transparent
contacts in photovoltaic solar cells (window layers). However, due to
potentially lower cost, material availability and easy chemical etchability
ZnO may become attractive also for display applications in the future.
This article deals with various aspects of ZnO window layers in
photovoltaic cells. The copper indium diselenide CuInSe2 (CIS) thin-film
solar cell uses a ZnO window with a thickness of about one m. Here ZnO
forms a heterojunction with a CdS buffer layer.2 In the superstrate
configuration of the amorphous silicon a-Si:H pin solar cell light enters
through the glass substrate which is covered by a thick ZnO layer. This layer
is chemically structured to enhance light scattering to enable the efficient use
of sun light.3 A third example is a solar cell on basis of a silicon wafer which
uses a heterojunction formed between a highly doped thin a-Si:H emitter
layer an the crystalline silicon wafer. In this case the task of the window
layer is twofold: enhancement of the lateral conductivity and reduction of the
reflection losses. The latter limits the thickness to 80 nm only.4 These
examples demonstrate the various aspects in the application: high lateral
conductivity, formation of heterojunctions (ZnO/CdS, ZnO/a-Si:H) and the
generation of textures for enhanced light trapping by chemical etching.
2.
199
10
9
21
6
20
19
10
18
10
II
10
P H (cm / Vs)
-3
n (cm )
10
10
15
20
25
30
III
f(O2) (sccm)
Figure 1. Electron concentration n and electron mobility measured by Hall effect of ZnO
films prepared by reactive magnetron sputtering as a function of the oxygen flow rate f(O2).
SEM micrographs taken on samples in the three regions of low, medium and high f(O2).8,9
In a reactive sputtering process the oxygen flow rate f(O2) is the most
relevant parameter. Fig.1 displays a typical example of the influence of f(O2)
on physical properties and structure. Hall effect measurements show that the
free carrier concentration n decreases continuously with f(O2) whereas the
electron mobility attains a maximum at medium values of f(O2). This
variation of the n and clearly reflects the change from metallic behavior at
low f(O2) (region I) to oxide formation (region III) at high f(O2) which is
related with an increase of the optical transmission T. These changes are
accompanied by structural variations in the ZnO layers. The SEM
200
Walther Fuhs
3.
TRANSPORT PROPERTIES
201
Intensity (a.u.)
(a)
2
ZnO
Zn(OH) x
1
as grown
536
534
532
530
528
-3
C H (cm )
10
22
10
21
10
20
10
19
(b)
50
100
150
200
250
depth (nm )
Figure 2. (a) XPS-Photoelectron spectra of ZnO deposited by magnetron sputtering taken with
an excitation energy of Eex = 1486 eV (O1s-spectra): (1) as deposited, (2) after annealing at
400 C; (b) Depth profile of the hydrogen concentration CH determined by nuclear reaction
technique (15N). Same ZnO samples as in (a).10
arguments it has been emphasized that the physical advantage lies in the
incorporation of particle energy from the plasma into the growing film. Two
directions have been pursued to obtain the low resistivity values at room
temperature:
202
Walther Fuhs
203
U :
:cm
cm
101
10-1
10-3
(a)
(cm/Vs)
10-5
100
10
1015
(b)
1017
1019
n (cm-3)
1021
Figure 3. Resistivity U (a) and Hall mobility (b) as a function of the carrier
concentration n.6 Data from various laboratories: ZnO single crystals, (, ) magnetron
sputtering impurity doped, # magnetron sputtering without impurity doping. The dashed line
is the result of a theoretical estimate for impurity scattering.14
Using glass substrates it has been observed that the resistivity increases
appreciably with decreasing layer thickness at values below of 100-300 nm
mainly due to a decrease of the electron mobility. A possible reason for such
behavior is the different growth and morphology in the early stages of film
growth.15
204
Walther Fuhs
4.
ZnO/SEMICONDUCTOR INTERFACE
In many applications not only the bulk properties of the ZnO layers are
important but also the morphology and electronic structure of the
ZnO/semiconductor interface which may have influence on the band profiles
of the solar cells, losses by contact resistance and instability.
t D = 500 s
Intensity (a.u.)
90s
SiO x
30s
Si(111):H
N an
104
102
100
98
205
nc-Si
a-SiOx
c-Si
Figure 5. Cross sectional TEM micrographs of the ZnO/c-Si interface.17
The complex structure of the ZnO/Si interface has also directly been
observed in cross-sectional TEM micrographs. An example is given in
Fig. 5. This micrograph clearly reveals the formation of a thin amorphous
film on the Si-surface (suboxide phase), the formation of nano-crystalline
phases followed by columnar ZnO growth.17
The above results showed that the ZnO/Si interface prepared by
magnetron sputtering is of high complexity and difficult to control. This
explains why all attempts to develop photovoltaic devices with reasonable
performance on basis of sputtered ZnO/Si heterojunctions have failed. This
may be a particular problem of the deposition method used. Recently it was
206
Walther Fuhs
shown that the formation of silicon suboxides at the interface does not occur
when the ZnO film is deposited by metal organic chemical vapor deposition
(MOCVD) from diethylzinc/H2O vapor on a heated substrate.10 The different
result obtained from these two methods may be explained by the fact that in
case of MOCVD only thermal energy is involved while magnetron
sputtering is characterized by additional energy input from the plasma into
the growing film. The quality of the interface may therefore strongly depend
on the energy of the species, which are responsible for film growth, and on
the ion-to-neutral ratio in the plasma.
5.
LIGHT TRAPPING
207
smooth (as deposited) and textured ZnO (short HCl-dip, optimized) window
layers. The cells had been prepared in the same run with an absorber layer
(a)
1.0
(b)
T
T, R
0.8
0.6
0.4
0.2
0.0
R
400
600
800
1000
wavelength (nm)
Figure 6. Optimized surface texture of ZnO-layer (a) and optical properties: transmission T
and reflectivity R (b).3
208
Walther Fuhs
1.0
1.0
0.9 1 m -Si
0.9
optimised
0.8
0.8
short dip
smooth
0.7
0.7
0.6
0.6
0.5
0.5
0.4
0.4
0.3
0.3
0.2
0.2
0.1
0.1
0
0
300 400 500 600 700 800 900 1000 1100
wavelength (nm)
absorbance
6.
CONCLUSION
Polycrystalline ZnO for window layers in solar cells are most frequently
prepared by magnetron sputtering. Optimized layers have resistivities of
1.5-2u10-4 cm with Hall mobilities of up to 60 cm2/Vs, refractive indices
of approximately 2.0 and an average transmittance of 85 % in the visible
range. The best values of the mobility in optimized films have reached the
209
intrinsic limit set by ionized impurity scattering. Optimized films grow with
columnar structure with the c-axis oriented perpendicular to the substrate.
Depending on the type of solar cell various aspects have to be considered:
high lateral conductivity, formation of heterojunctions (ZnO/CdS, ZnO/aSi:H) and the generation of textures for enhanced light trapping by chemical
etching.
ACKNOWLEDGMENT
The author gratefully to C. Pettenkofer for allowing the use of so far
unpublished results (Figs. 2 and 4) and to Ina Sieber and Frank Fenske for
fruitful discussion.
REFERENCES:
1. D. S. Ginley and C. Bright (guest editors), MRS Bulletin, 15 (2000).
2. J. Hedstrm, H. Ohlsen, M. Bodegard, A. Kylner, L. Stolt, D. Hariskos, M. Ruckh, and H.
W. Schock, Proc. 23rdd IEEE PVSC
C (1993).
3. B. Rech, J. Mller, T. Repmann, O. Kluth, T. Roschek, J. Hpkes, H. Stiebig, and W.
Appenzeller, Mat. Res. Soc. Symp. Proc. 762, A 3.1. (2003).
4. M. Schmidt, A. Schpke, O. Milch, T. Lussky, and W. Fuhs, Mat. Res. Soc. Symp. Proc.
762, A19.11.1. (2003).
5. K. Ellmer, J. Phys. D: Appl. Phys. 33, R1 (2000).
6. K. Ellmer, J. Phys. D: Appl. Phys. 34, 3097 (2001).
7. O. Kluth, G. Schpe, J. Hpkes, C. Agashe, J. Mller, and B. Rech, Thin Solid Films 442,
80 (2003).
8. S. Brehme, F. Fenske, W. Fuhs, E. Neubauer, M. Poschenrieder, B. Selle, and I. Sieber,
Thin Solid Films 342, 167 (1999).
9. F. Fenske, W. Fuhs, E. Neubauer, A. Schpke, B. Selle, and I. Sieber, Thin Solid Films
343, 130 (1999).
10. U. Meyer, Dissertation Universitt Cottbus (2004), U. Meyer and C. Pettenkofer,
unpublished.
11. S. F. J. Cox, E. A. Davis, S. P. Cottrell, P. J. C. King, J. S. Lord, J. M. Gil, H. V. Alberto,
R. C. Vilao, J. Pironto Duarte, N. Ayres de Campos, A. Weidinger, R. L. Lichti, and S. J.
C. Irvine, Phys. Rev. Lett. 86, 2601 (2001).
12. T. Minami, MRS Bulletin, 38 (2000).
13. T. Minami, H. Nanto, and S. Takata, Appl.Phys.Lett. 41, 958 (1982).
14. J. R. Bellingham, W. A. Phillips, and C. J. Adkins, J. Mater. Sci. Lett. 11, 263 (1992).
15. M. Birkholz, B. Selle, F. Fenske, and W. Fuhs, Phys. Rev. B 68, 205414 (2003).
16. J. Wennerberg, J. Kessler, and L. Stolt, Solar Energy Materials and Solar Cells 75, 47
(2003).
17. I. Sieber, N. Wanderka, I. Urban, I. Drfel, E. Schierhorn, F. Fenske, and W. Fuhs, Thin
Films 330, 108 (1998).
Chapter 18
ZnO/AlGaN ULTRAVIOLET LIGHT EMITTING
DIODES
E. V. Kalinina1, A. E. Cherenkov1, G. A. Onushkin1, Ya. I. Alivov2, D. C.
Look3, B. M. Ataev4, A. K. Omaev4, and C. M. Chukichev5
1
Abstract:
Key words:
1.
INTRODUCTION
212
E. V. Kalinina et al.
2.
EXPERIMENTAL PROCEDURE
3.
213
20
10
Intensity, a.u.
15
x7
ZnO
AlGaN
0
350
400
450
500
550
Wavelength, nm
Figure 1. Cathodoluminescence spectra of the ZnO and Al0.12Ga0.88N layers measured at room
temperature. Imax (ZnO) = 387 nm, FWHM (ZnO) = 21 nm.
214
E. V. Kalinina et al.
60
Current, mA
50
40
30
20
10
0
-1
-10
10
Bias Voltage, V
-20
10
EL intensity, a.u.
12
300 K
500 K
0
355
365
375
385
395
405
415
Wavelength, nm
425
435
445
Figure 3. Electroluminescence spectra of the n-ZnO/p-Al0.12Ga0.88N heterostructure of a lightemitting diode at temperatures of 300 and 500 K (Ic = 20 mA).
215
4.
CONCLUSION
REFERENCES
1. D. C. Look, Mater. Sci. Eng. B 80, 383 (2001).
2. D. M. Bagnall, Y. F. Chen, Z. Zhu, T. Yao, M. Y. Shen, and T. Gata, H, Appl. Phys. Lett.
73, 1038 (1998).
3. V. I. Zinenko and Yu. A. Agafonov, Solid State Electronics, 48, 2343 (2004).
4. A. E. Tsurkan, N. D. Fedotova, L. V. Kicherman, and P. G. Pasko, Semiconductors 6,
1183 (1975).
5. I. T. Drapak, Semiconductors 2, 624 (1968).
6. Ya. I. Alivov, J. E. Van Nostrand, D. C. Look, M. V. Chukichev, and B. M. Ataev, Appl.
Phys. Letters 83, 1 (2003).
7. B. M. Ataev, Ya. I. Alivov, E. Kalinina, V. V. Mamedov, G. A. Onushkin, S. Sh.
Makhmudov, and A. K. Omaev, ICPS2004, July 26-30, 2004, Flagstaff, Arizona, USA.
216
E. V. Kalinina et al.
Chapter 19
ZnO TRANSPARENT THIN-FILM TRANSISTOR
DEVICE PHYSICS
J. F. Wager
School of EECS, Oregon State University, Corvallis, Oregon, 97331, U.S.A.
Abstract:
Key words:
1.
INTRODUCTION
2.
The simplest approach for specifying the drain current-drain voltage (IDSVDS) electrical characteristics of a TTFT involves use of the square-law
model,9 as expressed by
I DS
V2
Z
Ci VGS VT VDSS DDS
L
2
217
(1)
218
J. F. Wager
Z
Ci
2L
(2)
IDS
pre-pinch-off
Z
Ci
L
IDS
VDS
VDS
2
x
post-pinch-off
Z
Ci
2L
IDS
2
VDS
VDSAT = VGS-VT
Figure 1. Drain current-drain voltage (IDS-VDS) curve at a given gate voltage (VGS),
illustrating the pre-pinch-off and post-pinch-off regimes of operation, which are distinguished
by VDSAT = VGS - VT.
PVDS
2S 2
GS
(3)
2 L
219
3.
There are three common mobility assessment methods,9 giving rise to (i)
effective mobility, Pefff, (ii) field-effect mobility, PFE, and (iii) saturation
mobility, Psat. Each type of mobility is now individually considered.
3.1
Effective Mobility
Z
Ci
L
VGS
VT VDDS ,
(4)
wII DS
D
wVDS V
GS
G
Z
Ci VGS VT
L
wP
VDDS ,
VDS
D
(5)
and solving for the mobility arising from this procedure, while neglecting the
term involving the change in mobility with respect to drain voltage yields,
Peff
gd
Z
Ci
L
(6)
VDDS o 0
Neglect of the term involving the change in mobility with respect to drain
voltage is appropriate for long-channel TTFTs for which velocity saturation
is not a relevant consideration.
Of the three common mobilities used for TTFT assessment, the effective
mobility is probably the best channel mobility figure-of merit since (i) it is
evaluated in the linear regime of operation where the channel is most
uniform and the concept of mobility as a proportionality constant between
the electric field and carrier velocity is most appropriate, and (ii) the term
220
J. F. Wager
3.2
Field-effect Mobility
wII DS
wVGS V
DS
D
(
Ci VDDS 1
L
GS
eff
) wP
.
wVGS
(7)
If the term involving the change in mobility with respect to a change in the
gate voltage is neglected, the defining equation for the field-effect emerges,
P FE
gm
Z
CiVDDS
L
VDDS
.
o
(8)
The most serious liability of PFE is the neglect in the derivation of PFE of
the term involving the change in mobility with respect to a change in the
gate voltage. This approximation is clearly not accurate at large gate
voltages where the accumulation channel thickness decreases and hence
electrons drifting in the channel experience more intense scattering at the
interface, resulting in a reduction of the channel mobility. Another PFE
demerit is that its estimation is sensitive to the magnitude of the source-drain
resistance, as is Peff. On the other hand, the use of PFE offers several
advantages, (i) it is evaluated in the linear regime of operation where the
channel is most uniform and the concept of mobility as a proportionality
constant between the electric field and carrier velocity is most appropriate,
(ii) VT is not required, and (iii) PFE is not singular when VGS = VT. The fieldeffect mobility is probably the most popular method used for channel
mobility estimation.
3.3
221
Saturation Mobility
Unlike Pefff and PFE, the saturation mobility, Psat, is evaluated from an
assessment of IDS in the post-pinch-off or saturation regime by taking the
square root of the post pinch-off expression for IDSAT, differentiating with
respect to VGS, and then solving the resulting equation for mobility to obtain,
P sat
w I DDSAT
wVGGS
Z
Ci
2L
,
VDS
GS
(9)
where the condition VDS > (VGS VT) ensures that this procedure is
accomplished in the post-pinch-off or saturation regime of the IDS-VDS curve.
Note that the numerator term is just the square of the slope of the linear
portion of a (IDSAT)1/2 VGS plot.
Of the three most common channel mobilities, Psat is probably used the
least. The most serious liability of Psatt is that it is evaluated in the postpinch-off regime of operation where the channel thickness and electric field
profile are highly non-uniform, the gradual channel approximation is not
valid, and the concept of mobility as a proportionality constant between the
electric field and carrier velocity is of questionable viability. However, Psat
offers several advantages, (i) VT is not required, and (ii) PFE is not singular
when VGS = VT, and (iii) its estimation is relatively insensitive to the
magnitude of the source-drain resistance, unlike Pefff and PFE.
3.4
P avg (
GS
gd (
Z
Ci (
L
GS
GS
G
,
ON
VDDS o 0
(10)
222
J. F. Wager
Pinc (
GS
wg
g d ( GGS )
wVGGS
Z
Ci
L
V
DS
D
(11)
o0
4.
OUTPUT RESISTANCE
IGS
CGD
+
VGS C
GS
S
223
IDS(total)
IDS
Rout
+
VDS
-
The important feature of Fig. 1 is that the output current, IDS(total), can
arise from either or both of the two output branches shown, and identified by
IDS and Rout. In the context of the present discussion, IDS represents the
square-law theory current-voltage characteristics, which can be written in
functional form as IDS(VDS,VGS), where VDS and VGS act as independent
variables. Thus, this IDS branch of the equivalent circuit represents transistor
behavior. In contrast, Routt is a simple resistor that accounts for the
conductivity of the bulk channel layer. If the channel layer is conductive,
Routt is small, most of the output current flows through Rout, and the TTFT
operates as a resistor. In contrast, if the channel is highly resistive, Routt is
large, and most of the output current flows through IDS so that the TTFT
functions as a transistor is supposed to function. Intermediate cases in which
the output current flows through both branches give rise to voltagecontrolled resistor operation. It is our experience that early-stage TTFT
development more often results in resistor or voltage-controlled resistor
behavior than transistor behavior, at least until the fabrication process is in
control.
ACKNOWLEGMENTS
This work was funded by the U.S. National Science Foundation under
Grant No. DMR-0071727 and by the Army Research Office under Contract
No. MURI E-18-667-G3.
224
J. F. Wager
REFERENCES
1. S. Masuda, K. Kitamura, Y. Okumura, S. Miyatake, and T. Kawai, J. Appl. Phys. 93, 1624
(2003).
2. R. L. Hoffman, B. J. Norris, and J. F. Wager, Appl. Phys. Lett. 82, 733 (2003).
t 82, 1117 (2003).
3. P. F. Carcia, R. S. McLean, M. H. Reilly, and G. Nunes, Appl. Phys. Lett.
4. J. Nishii, F. M. Hossain, S. Takagi, T. Aita, K. Saikusa, Y. Ohmaki, I. Ohkubo, S.
Kishimoto, A. Ohtomo, T. Fukumura, F. Matsukura, Y. Ohno, H. Koinuma, H. Ohno, and
M. Kawasaki, Jpn. J. Appl. Phys. 42, L347 (2003).
5. J. F. Wager, Science 300, 1245 (2003).
6. B. J. Norris, J. Anderson, J. F. Wager and D. A. Keszler, J. Phys. D 36, L105 (2003).
7. Y. Kwon, Y. Li, Y. W. Heo, M. Jones, P. H. Holloway, D. P. Norton, and S. Li, Proc.
Electrochem. Soc. 20003, 68 (2003).
8. R. L. Hoffman, J. Appl. Phys 95, 5813 (2003).
9. D. K. Schroder, Semiconductor Material and Device Characterization, 2nd ed., (Wiley,
New York, 1998).
Chapter 20
ZINC OXIDE THIN-FILM TRANSISTORS
E. Fortunato, P. Barquinha, A. Pimentel, A. Gonalves, A. Marques, L.
Pereira, and R. Martins
Materials Science Department, Faculty of Sciences and Technology of New University oof
Lisbon and CEMOP-UNINOVA, Campus da Caparica, 2829-516 Caparica, Portugal
Abstract:
Key words:
1.
Even though the field-effect transistor did not come into widespread use
until the 1960s, its invention predated both the junction and point-contact
transistors by many years. As it is normal with many innovations, its
practical realization was delayed until adequate materials and technologies
were available for its fabrication. We can even say that the Thin Film
Transistor (TFT) was the first solid-state amplifier ever patented. The basic
principle of the field-effect transistor (what we call now JFET) was proposed
by the first time by Julius Edgar Lilienfeld as early as 1925 and patented in
1930 Ref. (1) (see Fig. 1a) where an adaptation of the cross section of the
225
N.H. Nickel and E. Terukov (eds.),
Zinc Oxide A Material for Micro- and Optoelectronic Applications, 225238.
2005 Springer. Printed in the Netherlands.
226
E. Fortunato et al.
227
a)
Cu2S
Al
gate
drain
source
control gate
b)
insulator
CdS
source
drain
insulating substrate
a-Si:H
c)
source
drain
silicon nitride
gate
Fig. 1. Schematic cross sectional views of TFTs: a) Adapted from Lilienfeld,1 b) after
Weimer et al.,2 and c) after LeComber et al.6
228
2.
E. Fortunato et al.
Sol Gel
700 C (ref. 15)
Pulsed Laser
450 C (ref. 10)
0.2
107
not mentioned
enhancement
14
0.031-0.97
106
-1-2.5
enhancement
80
229
3.
EXPERIMENTAL DETAILS
3.1
ZnO optimization
The ZnO films (doped and undoped) were deposited onto soda lime glass
substrates by rf (13.56 MHz) magnetron sputtering using a ceramic oxide
ZnO target from Super Conductor Materials, Inc. with a purity of 99.99%
and 2 diameter. The sputtering was carried out at room temperature and the
argon deposition pressure was 0.15 Pa. The distance between the substrate
and the target was 10 cm and the rf power was varied between 50 and 175
W. The deposition rate was varied between 15 and 30 nm/min. The film
thickness was measured using a surface profilometer (Dektak 3 from Sloan
Tech.). The electrical resistivity (U) was measured as a function of
temperature in the range of 300 to 500 K using thermally evaporated
aluminum electrodes in a coplanar configuration. X-ray diffraction
measurements were performed at RT in air, using the Cu-KD line of a Rigaku
DMAX III-C diffractometer. The surface morphologies were analyzed using
a Field Effect Scanning Electron Microscope (FE-SEM, S-1400 Hitachi).
The optical transmittance measurements were performed with a Shimadzu
UV/VIS 3100 PC double beam spectrophotometer in the wavelength from
200 nm to 2500 nm.
3.2
ZnO-TFT fabrication
230
E. Fortunato et al.
Fig. 2. Schematic illustration of the ZnO based TFT structure. The ZnO channel layer and the
ATO gate insulator are 100 nm and 220 nm in thickness, respectively, while the ITO gate and
the GZO source and drain are 200 nm and 150 nm in thickness, respectively. The channel
length and gate width ratio varied between 1 and 65.
4.
231
RESULTS
Fig. 3 shows the dependence of the electrical resistivity and the average
optical transmittance at the visible spectra (between 400 and 700 nm), as a
function of rf power density (P). The highest resistivity (| 108 :cm) was
obtained for P = 5 W/cm2. For P around 5 W/cm2 the films became close to
stoichiometry with low structural defects and consequently higher resistive.
As we decrease or increase the rf power density from 5 W/cm2, a deviation
from stoichiometry is obtained accompanied by a decrease on the electrical
resistivity, due to a lower carrier concentration or/and electron mobility.
100
10
10
10
10
10
60
40
10
-1
10
-3
20
Transmittance (%)
Resistivity (:cm)
80
0
10
8
2
232
E. Fortunato et al.
0.35
8
8.9
0
0.30
7
7.7
0
0.25
6.4
0
0.20
0
0.15
0.10
0
5.
5.1
3.8
0.05
0
2.
2.6
31
32
33
34
35
36
37
(degrees)
(b)
Intensity (a.u.)
500
(002)
600
400
4
300
3
200
2
100
20
30
40
50
60
70
2T (degrees)
Fig. 4. (a) X-ray diffractograms of ZnO films corrected to the thickness of each sample
(ranging from 190 nm to 360 nm), deposited on soda lime substrates, using the Cu-KD line.
The large width of the peaks reveals a limited crystalline domain. For the film produced at 5
W/cm2, using the full-width at half-maximum of the peak at an angle of 34.1 , we calculated
an average crystallite size of the order of 10 r 1 nm using the Scherrer equation. The low
crystallite size is also associated with the fact that the films were processed at room
temperature without any type of post thermal annealing; (b) X-ray diffractogram of a ZnO
film produced with an rf power density of 5 W/cm2 for a wider 2T range. The film thickness is
130 nm.
233
Fig. 5. SEM micrograph with an apparent viewing angle of 40 of the same ZnO film
presented in Fig. 4(b).
Ci sattW
22L
L
ffor VDS
DS
VGS
G - VT
TH
(1)
234
E. Fortunato et al.
where Ci is the capacitance per unit area of the gate insulator. The obtained
Psat is around 20 cm2/Vs and the VTH is 21 V, showing that the ZnO-TFT
operates in the enhancement mode. Enhancement mode is preferable to
depletion mode since it is not necessary to apply a gate voltage to switch off
the transistor, because the circuit design is simpler and the power dissipation
is lower. The high value of Psatt agrees with the effective mobility (Pefff)18
determined for VDS d 2 V (in order to guarantee a uniform distribution of
carriers in the channel along the source to drain region) that has a value of
about 19 cm2/Vs. Here we also have to notice that the average field-effect
mobility (PFE) obtained from the transconductance18 data in the same VDS
region and for 30 V d VGS d 40 V is around 21 cm2/Vs. These data confirm
the high quality of the undoped ZnO films (improved crystallinity fraction of
the nanocrystals and low oxygen vacancies and/or Zn interstitials working as
donors) as well as the good channel-insulator interface obtained. The
magnitude obtained for VTH is directly proportional to the gate insulator
thickness. Reducing the gate insulator thickness could easily reduce these
values. The off-current is low, on the order of 1.310-9 A, and the on/off
ratio is 2105. That is, the TFT is able to sustain currents as high as 2.310-4
A. The gate voltage swing, S defined as the voltage required to increase the
drain current by a factor of 10 was 1.24 V/decade for the ZnO-TFT under
analysis. S is given by the maximum slope in the transfer curve (Fig. 6a).
Fig. 6(b) shows the IDS current curves as a function of the drain current
(VDS) for different VGS. The saturation was about 230 PA under a gate bias of
40 V. Besides the high value obtained for the saturation current, the device
exhibits hard saturation, evidenced by the flatness of slope of each IDS
curve, for large VDS, without the need of post heating treatments. This
indicates that the entire thickness of the ZnO channel layer is depleted and so
the device exhibits a high resistance (defined as the inverse of the slope of
the I-V curve for a given gate voltage, in the off saturation regime), which is
larger than 20 M:, as required for most circuit applications. The on
resistance (RON = 1/gd) defined as
RON
gd
1
L
Peff CiW ( GS -
TH
(2)
is of about 45 K: (VGS = 30 V). The relation between the ROFF and RON is
within what is expected for fast switching device behavior. Besides that, the
device response time, defined as being proportional to RONCi amounts to
values of about 25 ps/Pm2, which is highly promising for device
miniaturization.
-4
10
-5
10
-6
10
-7
10
-8
10
-9
(a)
0.015
1/2
10
(A )
-3
0.010
1/2
10
235
IDS
IDS (A)
0.005
10
20
30
0.000
40
VGS (V)
-4
2.5x10
VGS (V)
40
(b)
-4
2.0x10
IDS (A)
-4
1.5x10
35
-4
1.0x10
-5
5.0x10
30
25
20
10
20
30
40
VDS (V)
Figure 6. Typical ZnO-TFT characteristics with the channel layer deposited at room
temperature by rf magnetron sputtering, for a TFT with a width-to-length ratio of 1.4. (a)
Transfer characteristics for VDS = 20 V. The on/off ratio is 2105. The ZnO-TFT operates in
the enhancement mode with a threshold voltage of 21 V and a saturation mobility of 20
cm2/Vs. (b) Output characteristics for a ZnO-TFT. The saturation was about 230 PA under a
gate bias of 40 V. The ZnO-TFT exhibits hard saturation, evidenced by the flatness of slope
of each IDS curve, for large VDS. The dashed line represents the saturation drain current that
follows an exponential dependence on the voltage.
236
E. Fortunato et al.
100
Transmittance (%)
glass substrate
ITO film
80
60
40
ITO/ATO film
20
ZnO-TFT
0
500
1000
1500
2000
2500
Wavelength (nm)
Figure 7. Optical transmission spectra for the entire ZnO-TFT structure including the glass
substrate, were it is possible to observe that the ZnO-TFT is fully transparent to visible light.
For comparison the optical transmission spectra for the glass substrate (1.1 mm thickness), the
ITO, and the ITO/ATO films are also shown. The baseline and reference were measured in air
to quantify the total amount of transmitted light. The lower part shows a photograph of a glass
substrate (2.5 2.5 cm2) with ZnO based TFTs.
237
(7%). The absorbance at the near infrared region is associated with the
increase of the carrier concentration of the ITO and GZO films that leads to a
lower value associated to the plasma frequency.19
5.
CONCLUSIONS
ACKNOWLEDGEMENTS
We thank A. Pakkala from Planar Systems, Inc., Espoo, Finland, for
supplying the ITO/ATO substrates. The authors acknowledge A. Lopes for
the SEM analysis.
REFERENCES
1.
2.
2
3.
238
E. Fortunato et al.
Index
Acceptor 37, 69
Atomic force microscopy
Cathodoluminescence 83
Cluster 15
Columnar structure 15
Output conductance 217
dc-magnetron sputtering 15
Defects 133
Device physics 217
Dislocations 99
electrical activity 99
DLTS 47
Donor 37, 69, 115
hydrogen 115
Doping 197
Electrical properties 3
Electrical transport 197
Electroabsorption spectra 183
Electroluminescence 211
Electron affinity 115
Electron holography 99
Epitaxial layers 211
Electron paramagnetic resonance 69
photo-EPR 69
Exciton 69, 159, 171, 183
240
Magneto-luminescence 159
Magnetron sputtering 59
Metal organic vapor phase epitaxy 99
Microcrystals 83
Mobility 217
average 217
effective 217
field-effect 217
incremental 217
saturation 217
Molecular static method 183
MOVPE See Metal organic vapor phase
epitaxy
Muon spin rotation 115
Muonium 115
Muons 115
Nanocavity 83
Nanostructure 83
n-type 25, 37
O-H complexes 125
Optical spectra 171
Optoelectronic devices 225
Oxygen radicals 25
PEMOCVD
plasma enhanced MOCVD 59
Permitivity 171
Photoluminescence 59, 69, 83, 159
thermoluminescence 69
Puldes laser deposition 47
Point defects 69
Polarization 83, 171
Polycrystalline films 197
p-type 25, 37
Index
Raman spectroscopy 145
Recursion method 183
Reflectance 171
rf Magnetron sputtering 225
Schottky contacts 47
Structural properties 3
Thermal evaporation 59
Thermography 47
Thin-film transistors 217, 225
Transmission electron microscopy 99
Transparent electronics 225
Transparent thin-film transistor 217, 225
Ultraviolet emission 211
Valence band symmetry 159
Vapor-liquid-solid mechanism 15
Vibrational density of states 183
Vibrational modes 125
WGM See whispering gallery modes
Whisker 15
Whispering gallery modes 83
Window layers 197
Zeeman effect 159
Zn vacancy 133
ZnO synthesis 3
hydrothermal growth 3
melt growth 3
solution growth 3
vapor growth 3