Defects and Transport in Crystalline Solids
Defects and Transport in Crystalline Solids
OF OSLO
at
Department of Chemistry
Faculty of Mathematics and Natural Sciences
University of Oslo
2007
Version: September 2007
12
Preface
In 1972 Per Kofstad published his "Non-stoichiometry, diffusion and electrical
conductivity in binary metal oxides". It has been a popular textbook in defect chemistry
of oxides worldwide, not least because it contained a comprehensive review of defect
structure and defect-related properties of all binary metal oxides. It followed Kofstad's
equally well-recognized book "High temperature oxidation of metals" from 1967,
revised and published under the title "High temperature corrosion" in 1987.
In the 25 years that have passed since 1972 the field of defect chemistry in oxides
has grown and developed in many ways. New and improved measurement techniques,
new synthesis techniques, as well as skilled and devoted investigators have produced
better data for many oxides, partly encouraged by industrial developments in the use of
solid electrolytes (solid oxide fuel cells, gas separation membranes, ceramic electrodes,
catalysts, etc). Computer modelling enables detailed analysis of defect chemistry and
transport processes, and helps us interpret experimental data. Although the increasing
accuracy and complexity is welcome and unavoidable, the sound overview and ability
to analyse a situation may be obscured or lost. The newcomer may find the field
difficult to enter,
Along with the development in the application of theory and modelling, the field
has changed focus from binary oxides to ternary or higher oxides, much due to the
growing awareness that perovskite and related structures tolerate very large
substitutions and variations in valence and stoichiometry. Here we find lanthanum
aluminate insulators, titanate and zirconate electroceramics with ferroelectricity,
piezoelectricity, etc, the high-Tc cuprate and other superconductors, high-temperature
electronic conductors such as doped lanthanum cobaltate and chromite, high
temperature ionic conductors such as barium indate and doped lanthanum gallate,
mixed conductors, and alkaline earth cerates which dissolve protons and become proton
conductors at high temperatures.
A revision of the book from 1972 as such is thus utopic in many respects. In the
early 1990s Kofstad decided to author a new text that represented a development of the
general and theoretical first part of the former book, while he intended to include only a
few examples of data from important example oxide systems. After having sketched the
book he was struck by cancer in 1993. He fought the cancer for four years, with the
obligation to finish the book as one of the driving forces. However, deprived of much
of his usual strength, he was unable to finish the work. I became engaged as a co-author
at that time with the purpose of completing the book.
I used the new text as a curriculum when I took over the teaching of Kofstads
course in defect chemistry at the University of Oslo, thereby being able to adjust and
improve it step by step. However, the burdens of taking over the research group has
kept me from finishing the book. Furthermore, a textbook on defect chemistry by
Kofstads close friend and colleague, Donald M. Smyth was published that took away
some of the immediate need for our textbook. All in all I have felt that the new book
needs more work and refinement. With the present text I am still in the middle of that
process, and feedback is thus as welcome as ever.
Truls Norby, February 2006
13
14
Contents
Preface ......................................................................................................... 13
Copyrights and use of graphic material ..................................................... 14
Contents ....................................................................................................... 15
1.
Introduction............................................................................................... 110
Common crystal structures....................................................................... 113
Close-packing and voids.......................................................................... 113
Less close-packed arrangements .............................................................. 115
Some simple structures for oxides and other ionic compounds................. 116
Point defects............................................................................................... 121
Notation for description of point defects.................................................. 121
Notation for charged defects.................................................................... 122
Electroneutrality ....................................................................................... 124
Stoichiometry and nonstoichiometry........................................................ 125
Defects in stoichiometric compounds. ..................................................... 125
Defects in nonstoichiometric oxides. ....................................................... 127
Associations of point defects, point defect clusters, extended defects...... 128
Example cases of defect clusters.............................................................. 129
Extended defects, shear structures. .......................................................... 129
Modulated structures/infinitely adaptive structures .................................. 130
Summary ................................................................................................... 131
Literature................................................................................................. 131
Exercises ................................................................................................. 132
2.
DEFECT REACTIONS................................................................21
Introduction................................................................................................. 21
The three rules for writing defect reactions ............................................... 21
Mass balance ............................................................................................. 21
Electroneutrality........................................................................................ 22
Ratios of regular lattice sites...................................................................... 22
Examples of reaction equations for defects ................................................ 23
Stoichiometric compounds instrinsic disorders ....................................... 23
15
Introduction.................................................................................................... 3.1
Point defects and equilibrium in an elemental solid - a simple statistical
thermodynamic approach ....................................................................................... 3.2
Thermodynamics of chemical reactions ........................................................ 3.5
Thermodynamics and point defects............................................................... 3.8
Virtual chemical potentials of point defects .................................................. 3.8
Ideal and non-ideal models ........................................................................... 3.9
Examples of defect equilibria in pure, stoichiometric metal oxides ........... 3.10
Schottky defects ......................................................................................... 3.10
Frenkel defect pairs .................................................................................... 3.12
Intrinsic ionisation of electrons................................................................... 3.12
Defect equilibria in nonstoichiometric oxides ............................................. 3.14
Oxygen deficient oxides ............................................................................. 3.14
Oxide with excess metal ............................................................................. 3.17
Simultaneous presence of oxygen vacancies and interstitial metal ions ....... 3.18
Metal-deficient oxide.................................................................................. 3.20
Metal oxides with excess oxygen................................................................ 3.22
Examples of defect structures involving both oxygen deficiency and excess
................................................................................................................................ 3.22
Oxygen deficit............................................................................................ 3.23
Oxygen excess............................................................................................ 3.23
Stoichiometric condition............................................................................. 3.24
Summary ...................................................................................................... 3.26
Literature.................................................................................................... 3.27
Exercises .................................................................................................... 3.27
4.
Introduction.................................................................................................... 4.1
Constant concentration of impurity/dopant.................................................. 4.3
Schottky equilibria ....................................................................................... 4.3
Frenkel equilibria ......................................................................................... 4.4
Oxygen-deficient oxides doped with lower valent cations............................. 4.5
16
DIFFUSION................................................................................... 5.1
Introduction.................................................................................................... 5.1
Models of diffusion ......................................................................................... 5.1
Fick's first law .............................................................................................. 5.2
Potential gradients as the driving force ......................................................... 5.3
Example of chemical potential gradient acting on neutral particles ............... 5.4
Simplified model for one-dimensional diffusion ........................................... 5.5
Random diffusion......................................................................................... 5.8
Fick's second law .......................................................................................... 5.10
Measurements of the diffusion coefficients by tracer techniques................. 5.12
Diffusion mechanisms .................................................................................. 5.14
Vacancy mechanism................................................................................... 5.14
Interstitial mechanism ................................................................................ 5.15
Interstitialcy mechanism............................................................................. 5.15
Diffusion of protons in metal oxides. .......................................................... 5.17
Further considerations of factors that affect the diffusion coefficient in
crystalline solids .................................................................................................... 5.18
Free proton transport .................................................................................. 5.20
Temperature and oxygen pressure dependence of diffusion in metal oxides
................................................................................................................................ 5.21
Dependencies related to the concentration of defects .................................. 5.21
Temperature dependence of the attempt frequency .................................. 5.22
Resulting analysis of the diffusion coefficient Dr........................................ 5.24
17
Introduction.................................................................................................... 6.1
Transport in an electrical potential gradient ................................................ 6.1
Charge carriers in ionic compounds.............................................................. 6.3
The Nernst-Einstein relation between mobility and diffusion coefficient .... 6.4
Electronic conductivity in oxides ................................................................... 6.7
The band theory............................................................................................ 6.8
Effects of charged impurities or defects. Extrinsic semiconductors ............. 6.12
The Fermi level and chemical potential of electrons ................................... 6.17
Charge mobilities of electrons and electron holes ....................................... 6.18
Nonstoichiometric semiconductors ............................................................. 6.21
Ionic conductivity ......................................................................................... 6.25
Correlation effects: tracer diffusion and ionic conduction ........................... 6.26
Literature...................................................................................................... 6.28
Exercises ....................................................................................................... 6.29
7.
19
110
Figure 1-1. Schematic representation of some defects in an elemental solid. From A. AlmarNss: Metalliske materialer.
In addition to the structural defects, crystals also contain electronic defects, i.e.
electrons and electron holes that are relatively free to move in the crystal. The
electronic defects may either be formed through internal excitation of valence electrons
or they may be formed in association with point defects. If these electronic defects are
localised (trapped) at regular sites in the structure, the electronic defects are termed
polarons or - from a chemical point of view - valence defects. Defect electrons or
electron holes trapped at point defects often make otherwise transparent materials
coloured, and composite defects involving point defects and trapped electronic defects
are termed colour centres.
A complete description of the point and electronic defects in a compound and
their concentrations as a function of the activities (partial pressures) of the constituents
and the temperature is often termed the defect structure of the compound.
The reason for our interest in and concern with defects and imperfections is that
they determine a number of properties of crystalline solids. Thus diffusion or ionic
conductivity in crystalline solids takes place because of the presence of defects. Point
defects are responsible for lattice diffusion, dislocation diffusion takes place along line
defects and grain boundary and surface diffusion along planar defects. Solid state
diffusion, in turn, determines or strongly influences a number of properties or processes
such as mass transport in solids, solid state reactions, sintering, high temperature creep,
gas-metal reactions resulting in the formation of compact layers of reaction products,
etc. Electronic defects determine properties such as electronic conductivity, electrooptical properties, etc. Heterogeneous catalysis is closely related to the defect chemistry
of catalysts.
111
112
Figure 1-3. ABABAB(hcp) clos-packing (a) and ABCABC(fcc) close-packing (b) of spheres.
From Shriver and Atkins; Inorganic Chemistry.
The third layer (layer C) could alternatively be placed such that the spheres in this
layer will not be directly over either the atoms in the A or B layers, see (b) in the figure
above. The stacking sequence is now ABCABC... and this close packing is known as
cubic close packing (ccp) as this gives rise to a cubic (face-centred) crystal structure.
The following figure shows how the two close-packed structures appear when their
heaxonal and face-centred cubic unit cells are emphasized.
113
Figure 1-4. Schematic views of the hcp (a) and fcc (b) structure types. From Shriver and Atkins;
Inorganic Chemistry.
The close-packed structures contain two different types of empty space or voids
between the spheres, cf. the figures below. One type of void is termed a tetrahedral void
because each of these voids is surrounded by four spheres which centres are at the
corners of a tetrahedron. The other type of void is called an octahedral void. In this case
each void is surrounded by six spheres, three in one layer and three in the neighbouring
layer, and the centres of these six spheres lie at the corners of an octahedron. In the
close-packed structures there are two tetrahedral voids and one octahedral void per
atom (sphere). It may be noted that an octahedral void will fit a sphere with radius
0.414r where r is the radius of the bigger close-packed spheres. A tetrahedral void is
smaller and may fit a sphere with radius 0.225r. In different crystal structures of simple
metal oxides which can be considered to consist of a close-packing of oxygen atoms,
the metal ions often occupy all or part of the tetrahedral and/or octahedral interstices in
regular patterns.
Figure 1-5. Octahedral and tetrahedral hole in close-packed layers of speheres. From Shriver and
Atkins: Inorganic Chemistry.
114
Figure 1-6. Small spheres show locations of octahedral (left) and tetrahedral (right) voids in fcc
structure. From Shriver and Atkins: Inorganic chemistry.
Figure 1-7. Body-centered cubic (bcc) packing of spheres. From Shriver and Atkins: Inorganic
Chemistry.
Figure 1-8. Simple cubic (sc) packing of spheres. From Shriver and Atkins: Inorganic Chemistry.
115
Figure 1-9. The NaCl structure. Note that the centre anion is six-coordinated, as are also all
other ions. From Shriver and Atkins: Inorganic Chemistry.
The oxides MgO, CaO, SrO, BaO, CdO, CoO, NiO, FeO, TiO, NbO, and VO
possess this crystal structure. The last six of these oxides may be highly
nonstoichiometric and as such their NaCl structures highly defective.
In the sulfides FeS and NiS, the structure is similar, but the close-packing of
anions is hexagonal rather than cubic.
116
Figure 1-10. The CaF2 structure. (a) The cations form an fcc sublattice, while the anions form a
simple cubic arrangement. Every second anion cube has a cation in its centre. (b) The anions can be
taken to occupy tetrahedral voids in the cation sublattic.
117
Figure 1-11. The rutile TiO2 structure. From Shriver and Atkins: Inorganic chemistry.
In addition to rutile itself, the oxides SnO2, MnO2, VO2, MoO2, RuO2, GeO2 a.o.
possess regular or distorted rutile structures.
The corundum (-Al2O3) and ilmenite (FeTiO3) structures
The corundum structure is named after the high temperature modification of
alumina, -Al2O3. This structure may be described as a hexagonal close-packing (hcp)
of oxygen ions with the trivalent aluminium ions occupying 2/3 of the octahedral sites.
As the metal ions occupy octahedral sites, each metal ion is octahedrally coordinated
and surrounded by six oxygen atoms, while each oxygen atom is surrounded by four
metal ions. Two and two metal ions occupy neighbouring interstitial octahedral sites,
and the two corresponding AlO6 octahedra are linked by common faces. Other oxides
possessing the corundum structure are -Fe2O3, Cr2O3, Ti2O3, and V2O3.
A number of oxides with the nominal formula ABO3 have the corundum
structure when A and B have an average valence of 3 and are of approximately the
same size. The ilmenite structure with the nominal formula ABO3 and named after the
mineral FeIITiIVO3 is similar to the corundum structure in that the oxygen ions can be
described as hexagonal close-packed, the valences of the two cations have an average
valence of 3 and occupy 2/3 of the octahedral sites. However, in this case the two metal
ions (Fe2+ and Ti4+) are arranged in alternate layers.
118
Figure 1-12. Perovskite ABX3 structure. Left: Shown as a fcc packing of A (corners) and X (faces)
atoms, with the B in the center. Right: The A cation is 12-coordinated.( The A cation is depicted small,
but is in reality much bigger than the B cation). From Shriver and Atkins: Inorganic chemistry.
119
Table 1-1. Some structures of oxides based on close-packed oxygen ion sublattices.
Formula
MO
MO
4:4
M2O
8:4
M2O3,
ABO3
6:4
MO2
6:3
AB2O4
fcc of anions
hcp of anions
When the cations get somewhat larger they force the anions apart, and we get
non-close-packed structure types like the fluorite and its derivatives (pyrochlore, rareearth oxide structures, etc.)
When we involve even larger cations, being of the same size as the anions, the
large cations take part in the close-packing together with the anions. In this way we get
the perovskite structure and its derivatives.
120
Point defects
In order to understand and describe properties of inorganic compounds that are
dependent on the presence of point defects, it is necessary to be able to express the
concentrations of the defects under various conditions. This, first of all, requires a
system of notation to describe the point defects and the constituent atoms of the
inorganic compounds and their charges.
Figure 1-13. Schematic view of electronic structure of boron- and phosphorous-doped silicon.
From http://acre.murdoch.edu.au/refiles/pv/text.html
Defect electrons and electron holes that are free to move in the oxide have
effective negative and positive charges, respectively. They are written e / and h . If the
electron, for instance, is associated with a cation on a regular site - and may as such be
considered a valence defect - the defect may be written M M/ .
123
Electroneutrality
Crystals that we will deal with are considered to be electrically neutral. Then the
sum of all positive charges must equal the sum of all negative charges:
z[S
]=0
(1.1)
where the brackets denote concentration and z is the number and sign of charges on the
species Sz. For instance, in an oxide MO where the major defects are oxygen vacancies,
metal interstitials, and defect electrons, the electroneutrality condition can be written
2[v O ] + 2[M i ] - [e / ] = 0
or
2[v O ] + 2[M i ] = [e / ]
(1.2)
Despite its simplicity, and despite its being used in e.g. traditional aqueous
chemistry, the electroneutrality condition tends to cause both conceptual difficulties and
occasional errors (by misplacing the coefficients z). In order to help, we suggest that
one learns it by heart and uses it accordingly. Secondly, that one realises that it is not a
chemical reaction, but a mathematical relationship. It thus maintains no mass balance or
charge balance, it simply counts positive and negative charges from whatever
concentrations of charged species that are present and requires them to be equal.
124
125
Schottky disorder
A stoichiometric crystal with Schottky disorder contains equivalent
concentrations of cation and anion vacancies. A stoichiometric oxide MO contains
equal concentrations of metal and oxygen ion vacancies. In a stoichiometric oxide MO2
the concentration of oxygen vacancies is twice as large as that of the metal ion
vacancies. It may be noted that the formation of Schottky defects can only occur at
outer and inner surfaces or dislocations and will diffuse into the crystal until
equilibrium is reached.
Figure 1-14 Schottky defect pair in NaCl. From Shriver and Atkins: Inorganic Chemistry.
Frenkel disorder
A stoichiometric crystal with Frenkel disorder contains the same concentrations
of metal vacancies and metal interstitial ions. Contrary to the Schottky defects, Frenkel
defect pairs can be formed directly inside the crystal.
Figure 1-15. Cation Frenkel defect pair in AgCl. From Shriver and Atkins: Inorganic Chemistry.
Figure 1-16. Anti-site defect pair in CuAu intermetallic. From Shriver and Atkins: Inorganic
Chemistry.
127
128
Figure 1-17. Real and simulated image of shear plane in WO3-x. From Shriver and Atkins:
Inorganic Chemistry.
130
Summary
Due mainly to the gain in enthalpy of electrostatic attractions and electron orbital
overlap, solids usually organise their constituent ions, atoms, or molecules in a crystal
lattice. In such lattices the species take specific positions relative to each other, and this
pattern is repeated in 1, 2 or 3 dimensions. Different compounds take on different
geometrical arrangements depending on the properties of the constituents. This is a
complex function of the charge distribution of each atom or ion, in turn a function of
nuclear charge, ionisation energies, electron affinity and shape of the valence electron
orbitals. However, it is often useful to view the structures as close-packing or other
packings of rigid spheres, with systematic filling by smaller spheres of voids so created.
Because of the preferred occupation of specific sites in the lattice, together with the
tendency of distinct preferences in valence for the various elements, most compounds
are basically stoichiometric, that is, the ratio of the different atoms is a simple ratio of
small integers. In competition with this minimisation of enthalpy with perfect,
stoichiometric lattices, we have the tendency of increasing the entropy at T > 0 K and
this is done by the creation of defects. At any temperature there is thus an equilibrium
concentration of defects. The defects comprise point defects (vacancies, interstitials and
substitution), electronic defects (conduction electrons, holes, valence defects), and
associates, clusters, lines or planes of defects. The overall defect structure is a result of
energy minimisation by the enthalpy and entropy of formation and mixing, the
electroneutrality condition (charge balance), mass balance, and preservation of the
structure (site ratio balance). A defect is denoted by its constituent (a chemical element
symbol or a vacancy (v)), a subscript indicating lattice site or interstitial site (i) and a
superscript indicating effective positive or negative charge ( or ').
Literature
Frenkel, J. (1926) Z. Phys., 35, 652.
Krger, F.A. (1964) The Chemistry of Imperfect Crystals, North-Holland,
Amsterdam, and Wiley, New York.
Krger, F.A. and Vink, H.J. (1956) in Solid State Physics (eds F. Seitz and D.
Turnbell) 3, 307, Academic Press, New York.
Nray-Szab, I. (1969) Inorganic Crystal Chemistry
Budapest.
Akadmiai Kiad,
Rao, C.N.R. and Raveau, B. (1995) Transition Metal Oxides VCH, New York.
Schottky, W. (1935) Z. Phys. Chem., B29, 335.
Shriver, D.F., and Atkins, P.W. (1999) Inorganic Chemistry, 3rd ed. Oxford Univ.
Press., Oxford.
131
Exercises
1. Translate the defect names used in Fig. 1-1 and 1-2 into your first (or second)
language.
2. Calculate the number of atoms in the bcc and sc unit cells of Figs. 1-7 and 1-8.
(If you are not familiar with the procedures: Count all atoms, but divide by 2, 4,
or 8 for atoms that are shared with neighbouring cells by being in faces, edges
or corners.)
3. Calculate the number of atoms in the hcp and fcc unit cells of Figs. 1-4 (a) and
(b).
4. Calculate the number of ions of each kind in the NaCl and ABX3 cells in Figs.
1-9 and 1-12.
5. Calculate the number of atoms and voids in the cells in Fig. 1-6.
6. Calculate the volume filling ratio of balls in the structures of Figs. 1-4 (b), 1-7,
and 1-8. What would it be in 1-4 (a)?
7. Why would you say are metals generally heavier (have higher densities) than
ionic compounds (salts and ceramics)?
8. Sketch the periodic table of the elements, but include where possible the typical
binary oxides of the elements in their normal (most common) oxidation states,
in terms of formulae (stoichiometry) and structure type. Indicate also, based on
what you know about their tendency to take on neighbouring oxidation states,
whether the oxide is expected to be stoichiometric, have oxygen-deficiency, or
oxygen-excess.
9. List the main types of 0-, 1-, 2-, and 3-dimensional defects in crystalline solids.
10. Write the Krger-Vink notation for the following fully charged species in MgO:
Cation and anion on their normal sites, oxygen vacancy, magnesium vacancy,
interstitial magnesium ion.
11. Write the Krger-Vink notation for the following species in ZrO2: Cation and
anion on their normal sites, oxygen vacancy, zirconium vacancy, yttrium dopant
substituting Zr, interstitial carbon atom, cluster of yttrium dopant and oxygen
vacancy, nitrogen ion (N3-) substituting for oxygen ion.
12. Write the Krger-Vink notation for the following fully charged species in
CaTiO3: Calcium vacancies, titanium vacancies, oxygen vacancies, Ti ions on
Ca sites and vice versa, Ti interstitials.
13. Write the Krger-Vink notation and for the possible species and defects in a
material of your choice preferably a compound you are working with.
132
3 2
=0.74,
1-7:
3
8
=0.68,
1-8: ______,
7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. -
133
1-4a:
3 2
=0.74.
2. Defect reactions
2. Defect reactions
Introduction
As discussed in the previous chapter several different types of point defects
may be formed in metal oxides and other inorganic compounds. In principle, all
types of defects will be present, but in general, only a small number of different
defects will predominate.
When defect-dependent properties are to be interpreted, it is important that
the defect concentrations as a function of temperature and the activities of the
crystal components are known. For this purpose it is necessary to formulate
reactions and write equations for the formation (or annihilation) of the defects.
The defect formation may either occur internally in the solid or through
reactions with the environment. In the following, the rules for formulating defect
reactions will be described and applied to different defect structure systems, while
in the next chapter, conditions for equilibrium and equations relating equilibrium
defect concentrations to temperature, activities (partial pressures) of the
components in a compound, impurity concentrations, etc., will be discussed.
Mass balance
The defect reaction must balance with respect to the mass, i.e. the number
and types of atoms involved in the defect reaction must be the same before and
after the defect formation or annihilation. Vacancies, which only represent empty
21
2. Defect reactions
sites, have zero mass and do not count. Also electronic defects are considered not
to count in the mass balance.1
Electroneutrality
The compounds are and should remain electrically neutral. With the perfect
crystal as reference, this means that the total effective charge is the same before
and after the formation of the defects.
This means that the net charge on the left and right hand sides of a reaction
equation must be the same. This charge may be counted in terms of effective or
real charges. However, effective and real charges may in general not be balanced
against each other, and one must therefore avoid using both effective and real
charges in the same equation.
One may choose to count and balance also the mass of electrons, but then one must be
sure to include also non-defect electrons (electrons in the valence band).
2
One may choose to operate with sites or energy levels - for electronic defects, and one
must then take into account the emty levels of the conduction band and the occupied levels of the
valence band.
22
2. Defect reactions
crystallite sizes for which the number of surface atoms is insignificant to that of
the bulk atoms.
Defect reactions generally do not result in significant changes in the number
of surface atoms. In the formation of a vacancy, for instance, a bulk atom is
transferred to the surface, but in the same process a previous surface atom, in turn,
becomes a bulk atom. (It is worth mentioning, however, that the creation of the
new bulk atom has consequences for the energy involved, as we shall discuss in a
later chapter.)
Schottky disorder
As described in the previous chapter, the Schottky disorder involves the
presence of equivalent amounts of cation and anion vacancies. In an oxide MO
this means that the crystal contains equal concentrations of metal and oxygen
vacancies. The overall formation of such a defect pair within the crystal involves
the transfer of a pair of cations and anions on regular lattice sites from the bulk to
the surface. In reality the defects are formed at external and internal surfaces or
23
2. Defect reactions
dislocations and subsequently diffuse into the crystal until they are randomly
distributed. In writing the equation, one is only interested in the initial and final
states, and one disregards the kinetics of the defect reaction.
If one starts with a pair of cations and anions on regular lattice sites within
the crystal, MM and OO, one must also take into account that the formation of the
Schottky pair results in the formation of two new lattice sites, and the overall
equation may thus be written
M M + OO = v M// + vO + M M + OO
(2.1)
However, in this equation MM and OO on both sides may be cancelled, and the net
equation therefore becomes
0 = v M// + vO
(2.2)
Frenkel disorder
For the Frenkel disorder the predominant defects are either limited to the
cations and anions, and the disorder involves the presence of equal numbers of
vacancies and interstitial ions in a sublattice in a crystal. In the formation of a
Frenkel defect pair, a cation on a normal site is transferred to an interstitial site,
and no new lattice sites are created in the process. If, for the sake of illustration,
the interstitial ion and the resulting vacancy are assumed to be doubly charged, the
formation of a Frenkel defect pair may be written
M M = v M// + M i
(2.3)
24
2. Defect reactions
0 = e/ + h
(2.4)
e x = e / + h
(2.5)
We can also choose to keep some kind of track of the sites or energy levels for
electrons;
(2.6)
but both the two latter reactions represent levels of detail that we normally need
not consider.
As mentioned earlier, the electronic defects may be localized (valence
defects). In these cases the reactions are connected with individual atomic sites.
For instance the intrinsic ionization (disproprtionation) of Fe3+ ions into Fe2+ and
Fe4+ ions would be written
x
/
2 Fe Fe
= Fe Fe
+ Fe Fe
(2.7)
Similarly, intrinsic ionization may also take place by charge transfer. In ilmenite,
FeTiO3, we may for instance observe the reaction Fe2+ + Ti4+ = Fe3+ + Ti3+, which
in terms of defects is written:
x
Fe Fe
+ TiTix = Fe Fe
+ TiTi/
(2.8)
Nonstoichiometry
The ratio of cation to anion lattice sites is the same whether a compound is
stoichiometric or nonstoichiometric. But as nonstoichiometry means that there is
an excess or deficit of either cations or anions, nonstoichiometry also means that
there is an excess of a certain type or types of defects relative to that in the
stoichiometric condition. If the predominating type of defects are charged,
complementary electronic defects are created in order to conserve electrical
neutrality. The extent of nonstoichiometry and the defect concentrations in
25
2. Defect reactions
inorganic compounds are functions of temperature and activities (partial
pressures) of their components.
First we consider, as an example the formation of oxygen deficiency in an
oxide MO2. The overall reaction may be written
MO2 ( s) = MO2 x ( s) + 2x O2 ( g )
(2.9)
From this equation it is qualitatively seen using le Chatelier's principle that the
oxygen deficit increases with decreasing oxygen pressure. Conversely, for oxides
with excess oxygen the nonstoichiometry increases with increasing oxygen
pressure.
In the following various defect reactions which are encountered in oxides
with oxygen or metal excess or deficit will be considered. As the activity of the
metal component is usually negligibly small compared to that of the oxygen
activity under most experimental conditions, the nonstoichiometry in metal oxides
is correspondingly a result of the interaction and exchange of oxygen between the
metal oxide and the surrounding gas atmosphere. In the following examples of
formation of nonstoichiometric defects, only cases where the metal oxides interact
with gaseous oxygen are illustrated. However, it should be borne in mind that the
corresponding nonstoichiometry and defects may be formed by interaction with
metal if it is experimentally feasible to control the activity of the metal component
in the surroundings of the crystal.
Oxygen-deficient oxides
An oxygen vacancy is formed by the transfer of an oxygen atom on a
normal site to the gaseous state. No change in the number of lattice sites takes
place. This defect reaction may be written
(2.10)
OOx = vOx + 12 O2 ( g )
In this equation it is assumed that the oxygen vacancy is neutral, i.e. the two
electrons of the O2- ion are associated with the vacancy or its immediate
neighbourhood and the vacancy has zero effective charge. As described in Chapter
1, the two electrons trapped at or near the vacancy may, depending on the
temperature and vacancy concentration, be excited and transferred away from the
vacancy. Correspondingly, the oxygen vacancy acts as a donor and becomes
singly and doubly charged:
vOx = vO + e /
(2.11)
vO = vO + e /
(2.12)
26
2. Defect reactions
The formation of the doubly ionized oxygen vacancy can be written as the total
reaction:
OOx = vO + 2e / + 12 O2 ( g )
(2.13)
OOx + 2 M Mx = vO + 2 M M/ + 12 O2 ( g )
(2.14)
If this takes place in the oxide MO2 the valence of the M atoms then partly reduce
from +4 to +3 when the oxide becomes non-stoichiometric (MO1-x).
M Mx + 2OOx = M ix + O2 ( g )
(2.16)
From this we can state that the metal interstitial, like the oxygen vacancy, is an
electron donor. It may be noted that in both cases it is the effectively neutral or
incompletely ionized defects that is a donor, not the fully ionized defects.
We have until now written the formation of metal defects by exchange with
oxygen in the surrounding atmosphere. We may in principle choose to write the
same in terms of exchange with the metal component in the surroundings. For
27
2. Defect reactions
instance, metal insterstitials in MO can be formed in equilibrium with the metal
vapour:
M ( g ) =M i + 2e /
(2.17)
Metal-deficient oxides
In metal-deficient compounds metal vacancies are the predominating point
defects. In an oxide MO a metal vacancy is formed by reacting oxygen gas with
the oxide. This creates new oxygen lattice sites and therefore an equivalent
number of new metal lattice sites - which are vacant - is also formed:
1
2
(2.18)
O2 ( g )= v Mx + OOx
In this reaction it is assumed that the vacancy that is formed is effectively neutral,
that is, two electrons are taken from the surroundings of the vacancy in order to
form the O2- ion. The lacking electrons may be taken up from the valence band,
forming electron holes here;
v Mx = v M// + 2h
(2.19)
v Mx + 2 M Mx = v M// + 2 M M
(2.20)
This correspondingly means that the valence for some M atoms changes from +2
for to +3.
The metal vacancy, whether it is effectively neutral (electrons missing from
its neighbourhood) or ionized (electrons supplied from the rest of the crystal) it
illustrates that the effective size of a point defect may extend beyond the site
itself; the charge may be distributed over at least the closest neighbours of the
defect. Accordingly, defect chemistry works fine under the ionic model and with
the assumption of integer charges, even if the compound is far from ideally ionic
the point defect is just a little bigger than one site. It works because in defect
chemical as in other chemical reactions, the electrons have to choose to go or stay
no half electrons are involved. It then does not matter how close that electron
was of the point defect as long as it was associated with it. For most defect
chemical considerations, it then also suffices to consider the defect as a true point
defect, and that the electrons are associated with the point defect itself.
28
2. Defect reactions
(2.21)
O2 ( g ) = Oix
No new lattice sites are formed in this reaction. The neutral interstitial oxygen
atoms may in principle be ionised to yield electron holes and interstitial oxygen
ions with negative effective charges, so that the total reaction becomes, for
instance,
1
2
O2 ( g ) = Oi// + 2h
(2.22)
BSix = BSi/ + h
(2.23)
PSix = PSi + e /
(2.24)
29
2. Defect reactions
The neutral defects, the acceptor and donor, were here formed by dissolving the
foreign elements in the same oxidation state as the host, namely 0. It is the desire
of the crystal to fulfil the octet rule around each atom the covalent nature of the
bond that causes the ionization in this case.
We may do something similar with ionic compounds such as oxides. For
instance, hypothetical LiO, containing divalent lithium, Li2+, may be dissolved in
2+
NiO to form effectively neutral Li Ni
species; Li Nix . Since Li is unstable in
oxidation state +2, it easily accepts an electron from the valence band of the
crystal:
/
Li Nix = Li Ni
+ h
(2.25)
Thus, in this case, it is the formal oxidation state according to an ionic model that
is the argument for why the acceptor is ionised. Similar arguments hold for donors
in ionic compounds.
When dealing with aliovalent foreign elements it is common to assume
directly the ionisation of the donor or acceptor into an effectively charged state.
We often refer to the cases as donor- or acceptor-doped systems, but keep in mind
that it is the unionised species that constitutes the donor or the acceptor.
(2.26)
and thus the dissolution of Mh2O3 in M1-yO increases the concentration of metal
vacancies in the parent oxide.
210
2. Defect reactions
But as will be discussed in detail in the next chapter an increased
concentration of charged metal vacancies simultaneously results in a decrease in
the concentration of electron holes. The dissolution of Mh2O3 may thus
alternatively be written with a defect reaction involving annihilation of electronic
holes:
(2.27)
When the dissolution of Mh2O3 is written in this way, the equation emphasises
that the concentration of electron holes is reduced. However, as the concentration
of metal vacancies and electron holes are interrelated, it should be emphasised that
both processes take place. Actually, the dissolution of the higher valent cation
decreases the concentration of all effectively positive defects and increases the
concentration of all negatively charged defects. Qualitatively, this may simply be
seen as a consequence of the electroneutrality equation or, if one wishes, of Le
Chatelier's principle. Details of these aspects will be further discussed using defect
diagrams in Chapter 4.
If one has the choice of writing a doping reaction by creation or
annihilation of defects, it is usually more meaningful to choose the creation of
defects, since this can describe the doping reaction to levels beyond the defect
concentrations that existed in the undoped oxide, and since this describes the
defects that will dominate when the doping level gets high.
So far we have added a higher valent dopant oxide to M1-yO. Let us
alternatively consider what takes place when a lower valent dopant oxide, Ml2O,
is added to M1-yO. If the Ml+ ions dissolve substitutionally in M1-yO, the
dissolved Ml+ ions have one negative effective charge, Ml M/ . This will be
compensated by the formation of positive effective charge or the annihilation of
negative effective charge. In our example oxide this will mainly affect electron
holes and metal vacancies, and the two defect reactions may thus be written
Ml 2 O + 12 O2 ( g ) = 2 Ml M/ + 2h + 2OOx
(2.28)
Ml 2 O + v M// = 2 Ml M/ + OOx
(2.29)
Thus, the addition of Ml2O to M1-yO has the opposite effect of addition of
Mh2O3; it increases the concentration of electron holes and decreases the
concentration of metal vacancies. One may note that the first of the two reactions
is an oxidation and thus involves uptake of oxygen, while the latter is not. It is
typical of such dissolutions of foreign aliovalent elements in oxides that they can
be charge compensated by means of electronic defects in which case we are
dealing with reduction or oxidation and uptake or release of oxygen or by means
of point defects in which case we have no change in oxidation states and no
exchange of oxygen gas.
211
2. Defect reactions
.
This
charge must be compensated either by formation of negative effective
MhM
charges in the form of electrons or by annihilation of positive effective charges,
oxygen vacancies. These defect reactions may be written
Mh5+,
(2.30)
(2.31)
Thus when MO2-x is doped with a higher valent oxides Mh2O5 the
concentration of electrons is increased and the concentration of oxygen vacancies
is decreased.
When the same oxide, MO2-x, is doped with a lower valent oxide, Ml2O3,
the negative effective charge of the dissolved atoms Ml M/ is compensated by
annihilation of electrons or formation of oxygen vacancies:
Ml 2 O3 + 2e / + 12 O2 ( g ) = 2 Ml M/ + 4OOx
(2.32)
Ml 2 O3 = 2 Ml M/ + vO + 3OOx
(2.33)
(2.34)
2. Defect reactions
occupies. It may be noted that a foreign cation dissolving interstitially always will
have a positive effective charge and thus affect the defect structure in a similar
manner as a higher valent cation dissolved substitutionally. We shall return to this
and other examples in subsequent chapters and in connection with the review and
discussion of a few individual oxides systems.
(2.35)
3OOx + 2 NH 3 ( g ) = 2 N O/ + vO + 3H 2 O( g )
(2.36)
H 2 + 2OOx = 2OH O + 2e /
(2.37)
or
H 2 O( g ) + 2OOx = 2OH O + 2e / + 12 O2 ( g )
(2.38)
in which case the protons dissolved are charge compensated by the formation of
defect electrons. Of course, the electrons may interact with other defects in the
oxide so that the protons in effect are compensated by formation of other negative
defects or by the annihilation of positive defects. From the dissolution reaction
and through the interaction with native defects in the oxide it is clear that the
dissolution of hydrogen in metal oxides is dependent both on the partial pressure
of the hydrogen source (e.g. water vapour or hydrogen) and of oxygen. These
aspects will be described in more detail in a later chapter.
The dissolution of hydrogen from its oxide, H2O, is in principle similar to
dissolution of other foreign cations. However, the possibility of a controlled water
vapour pressure and the fast diffusion of protons makes it much easier to attain
213
2. Defect reactions
and vary (and more difficult to completely avoid) an equilibrium content of
protons in the oxide.
In terms of defect chemistry the dissolved proton, located on a normal
oxygen ion as hydroxide, may also be considered to constitute an interstitial
hydrogen ion, and as such it is also in the literature alternatively written H i . One
just has to bear in mind that the protons do not occupy regular interstitial positions
(voids).
Under reducing conditions, where hydrogen is stable in oxidation state 0 (as
H2 in the gas phase) it might be possible to foresee neutral hydrogen atoms
dissolved in oxides, probably interstitially, as H ix . Under even much more
reducing conditions could also hydride ions be expected to become stable, e.g. as
dissolved substitutionally for oxygen ions, as the defect H O .
(2.39)
//
O2 ( g ) = vCa
+ vTi//// + 3OOx + 6h
(2.40)
Similar principles should be applied also in cases where one and the same
element is distributed on different crystallographic sites. For instance, Y2O3 has a
structure where all oxygen ions are not strictly equal. Similarly, distorted
perovskites may have unequal oxygen sites. In the pyrochlore structure, A2B2O7,
there are 6 oxygen sites of one type and 1 of slightly different coordination and
energy (and one which is structurally empty and thus to be regarded as an
interstitial site). In principle the fomation or annihilation of crystal units has to
214
2. Defect reactions
maintain the ratio between those different sites in all such cases. However, this is
so far hardly ever practiced in defect chemistry.
Contrary to binary oxides, ternary and higher oxides can have nonstoichiometry not only in terms of the oxygen-to-metal ratio, but also internally
between the various cations. This is in practice often a result of synthesis. For
instance, it may be difficult to weigh in exactly equal numbers of moles of Ca and
Ti precursors when synthesizing CaTiO3, so that the synthesized material has a
permanent number of vacancies on one of the cation sites. Such non-stoichiometry
may also be a result of equilibria. For instance, if A-site deficiency is energetically
favourable over B-site deficiency in the compound ABO3, we may at very high
temperatures (e.g., during sintering) see a preferential evaporation of the A
component:
AAx + OOx = v A// + vO + AO ( g )
(2.41)
(2.42)
It may be noted that these reaction equations do not violate the site ratio
conservation requirement of the ternary oxide.
When we earlier doped elementary or binary compounds the reaction was
fairly straightforward. When we dope a ternary or higher compound, however, the
reaction may be less obvious we have some choices. It is quite common,
however, to do the synthesis and write the equation in such a way that one takes
out a corresponding amount of the host element that is substituted. If we, for
instance, want to dope LaScO3 with Ca substituting for La, we go for a
composition La1-xCaxScO3. In order to see how we write the doping reaction in
this case we first just look at the trivial normal synthesis:
1
2
(2.43)
x
La 2 O3 + 12 Sc 2 O2 = La La
+ Sc Scx + 3OOx
Accordingly, we then write the defect reaction for the doping in the way that we
let there be Sc2O3 reserved for the CaO:
/
CaO + 12 Sc 2 O3 = Ca La
+ Sc Scx + 52 OOx + 12 vO
215
(2.44)
2. Defect reactions
Summary
Defect reactions can be written in a similar manner as ordinary chemical
reactions. They must confer with the requirements for conservation of mass,
charge, and ratios of lattice sites, but are allowed to increase or decrease the total
number of lattice units.
We have treated defect reactions mostly in cases for binary ionic
compounds, notably oxides. Reactions for elemental crystalline solids follow the
same basic rules and principles, but without the site ratio conservation
requirement. Also higher compounds, e.g. ternary oxides, follow the same rules,
only often with slightly more complex site conservation considerations and more
defects.
Normally, one should seek to use reactions that create rather than annihilate
defects, and one should avoid using reactions that form more defects than
necessary (as such reactions are then a sum of simpler reactions).
We have formulated defect reactions which describe intrinsic ionic and
electronic disorder, nonstoichiometry, variable ionisation of point defects, and
substitutional dissolution of aliovalent cations and ions and anions. Aliovalent
elements may be compensated by electronic defects or by point defects, of which
the former involve red-ox-reactions.
We have treated hydrogen defects specially, as they arise from a special
source and in its most stable form, the proton, take on a special size and type of
defect.
Exercises
The excercises in this chapter mainly provide training in 1) formulating a
reaction in terms of the appropriate reactants and products, and 2) balancing and
checking them with respect to a) mass conservation, b) charge conservation, and
c) site ratio conservation.
1. Write a reaction for the formation of Schottky defects in NaCl.
2. Write a reaction for the formation of Schottky defects in MO2.
3. Write a reaction for the formation of Schottky defects in Cu2O.
4. Write a reaction for the formation of anion Frenkel defects in MO.
5. Write a reaction for the formation of anion Frenkel defects in CaF2.
6. Write a reaction for a charge transfer between the cation and anion in
CeO2, i.e. for reduction of the cerium ion and oxidation of the oxygen ion.
Write the same process as an intrinsic ionisation assuming delocalised
electronic defects.
216
2. Defect reactions
7. Write a reaction for the formation of fully ionised oxygen vacancies and
electrons when oxygen is lost in the reaction M2O3 = M2O3-x + x/2 O2(g).
8. Write reactions for the direct formation of fully ionised metal interstitials
in MO, M2O3, and MO2.
9. Write reactions for the direct formation of fully ionised metal vacancies in
MO, M2O3, and MO2.
10. Write reactions for formation of sodium vacancies in NaCl by exchange
with a) chlorine gas and b) sodium gas.
11. We have in the text claimed that charges can be associated directly with
ideal point defects, even if the compound is not ideally ionic, and that this
normally does not affect the results of the calculations. Under what
conditions will it start to matter that the charge of the defect actually
extends outside the defect itself?
12. Write reactions for dissolution of CaO into the oxide ZrO2-x
13. Write reactions for dissolution of MfO into M2-xO3 when we assume that
the foreign metal Mf is small and therefore dissolves interstitially.
14. Write reactions for dissolution of CaO substitutionally into the anionFrenkel dominated Y2O3.
15. Write reactions for dissolution of ZrO2 substitutionally into the anionFrenkel dominated Y2O3.
16. Write reactions for dissolution of CaF2 into CaO, assuming that Fdissolves substitutionally and that CaO is dominated by Schottky defects
in the pure, undoped state.
17. Write a reaction for dissolution of protons in an oxide with water vapour
as source and with oxygen interstitials as compensating defects.
18. Write reactions for dissolution of protons in an oxide M1-xO with water
vapour as source.
19. Consider an oxide M2O3 which is acceptor-doped by dissolving the oxide
MlO. The dopants are compensated by oxygen vacancies. Write a defect
reaction for the dissolution of MlO. Write a reaction for the dissolution of
protons from water vapour into this acceptor doped oxide, by annihilation
of the oxygen vacancies. Finally, write a reaction for the dissolution of
MlO in the presence of water vapour and without forming oxygen
vacancies. Can the latter reaction be constructed as a sum of the two
previous reactions?
20. Hydride ions have been suggested to dissolve in oxides accompanied by
protons by disproportionation of hydrogen gas. Write defect reactions for
this in the case that the hydride dissolves a) interstitially, and b)
substitutionally. (Assume in both cases that the oxide is perfectly
stoichiometric and has no defects as a starting point).
21. Write defect reactions for formation of A-site Frenkel, B-site Frenkel, and
anion Frenkel disorders in a perovskite ABO3 (assuming both cations are
trivalent).
217
2. Defect reactions
22. Write a defect reaction for the formation of Schottky disorder in spinel
MgAl2O4.
23. Suggest a defect reaction for formation of Schottky-type disorder in
LaMnO3 when it is assumed that cation vacancies are formed in the La
sublattice only.
24. Write a defect reaction for dissolution of Sc2O3 in CaTiO3 during
synthesis, assuming that Sc substitutes Ti forming acceptors compensated
by oxygen vacancies.
218
3. Defect equilibria
3. Defect equilibria
Introduction
After formulating the correct defect reactions in a compound, we turn our
interest towards the equilibrium relations for those reactions. These correlate the
equilibrium concentrations of the different defects with temperature, activities
(partial pressures) of the components in the compound, and other parameters
which affect the defect structure. The defect equilibria will here generally be
described by the law of mass action.
The change in the Gibbs (free) energy of a chemical reaction in a closed
system is given by3
(3.1)
G = H TS
where H is the change in enthalpy and S the change in entropy of the reaction
and T is the absolute temperature.
The reaction will proceed spontaneously as long as G < 0, until G has
reached its minimum. Here, G = 0, and this relation represents the condition for
chemical equilibrium.
The creation of single, unassociated point defects in an elemental,
crystalline solid increases the internal energy of the system and the enthalpy of the
defect formation is positive. But the configurational entropy of the system also
increases, and the equilibrium concentration of the defects will be reached when
the Gibbs energy of the system is at minimum. Thermodynamically, point defects
will thus always be present in a crystal above 0 K.
All types of defects will in principle be formed. However, the free energies
of formation of the different types and systems of defects can have widely
different values, and correspondingly it is often found that certain defects
predominate in a particular solid. The relative concentrations of the different types
of defects will be a function a temperature and other variables. Thus, defect
equilibria with large positive enthalpy of formation, for instance, which are not
favoured at low temperatures, may become important at high temperatures.
From a thermodynamic point of view a solid containing point defects
constitutes a solid solution where the point defects are dissolved in the solid. In
analogy with liquid solutions the solid may be considered to be the solvent and the
point defects the solute. Similarly the defect equilibria may be treated in terms of
the thermodynamics of chemical reactions and solutions.
A closed system is one where the mass (chemical entities) of the system is confined, but
where energy (heat) can be exchanged with the surroundings. The expression for the Gibbs energy
arises as a result of the definition of entropy and of the first law of thermodynamics which
preserves the total energy of the system and the surroundings.
3.1
3. Defect equilibria
(3.2)
G = nv ( H TS vib ) TS conf
S conf = k ln W
3.2
3. Defect equilibria
where k is Boltzmann's constant. The thermodynamic probability W represents in
this case the number of different distinguishable ways whereby nv vacancies may
be distributed on N+nv lattice sites and is given as
W=
( N + nv )!
N ! nv !
(3.5)
For large numbers of N and nV, which is typical for real crystals, Stirling's
approximation ( ln x!= x ln x x for x>>1) may be applied, whereby Eq. 3.4 takes
the form
S conf = k ( N ln
N + nv
N + nv
+ nv ln
)
N
nv
(3.6)
This entropy scales with the size of the system (the number N) but is not
proportional with the number of defects nv.
Intuitively, and by insertion in the equation above, we have Sconf = 0 at nv
= 0. Thus, Sconf = Sconf and by insertion into the expression for the change in free
energy we get
G = nv (H TS vib ) kT ( N ln
N + nv
N + nv
+ nv ln
)
N
nv
(3.7)
The changes in enthalpy, entropy and Gibbs energy are illustrated schematically in
Fig.3.1.
3.3
3. Defect equilibria
Figure 3-1. Example of the variation in enthalpy, vibrational and configurational entropy, and
Gibbs energy with the concentration of vacancies in an elemental solid.
G
nv
= H TSvib + kT ln
=0
n V
N + nv
(3.8)
nv
represents the fraction of the total number of sites which are
N + nv
vacant and is through rearrangement of Eq.3.8 given by
The term
nv
S
H
= exp( vib ) exp(
)
N + nv
k
kT
(3.9)
nv
S
H
= exp( vib ) exp(
)
N + nv
k
kT
3.4
(3.10)
3. Defect equilibria
nv
represent the site fraction of vacancies in the
N + nv
S vib
H
G 0
) exp(
) = exp(
)
R
RT
RT
(3.11)
We will see later that this standard Gibbs energy change is the one we have when
the reaction proceeds with all reactants and products in their standard reference
state (activity of unity). It is not the same as the general Gibbs energy change
calculated for the progression of the defect formation reaction earlier.
i = (
G
) T, p,n ,n ,...
n i
1
(3.12)
3.5
3. Defect equilibria
(3.13)
When a chemical reaction takes place, the numbers of various constituents change,
and dG changes correspondingly. For an open system (i.e. where we may add or
remove heat and volume) at constant temperature and pressure, dG is given by
dGT , p = 1 dn1 + 2 dn 2 + .... i dni
(3.14)
At equilibrium we have
dGT , p = i dni = 0
(3.15)
(3.16)
aA + bB = cC +dD
the change in Gibbs energy is given by the difference in the total free energy in
the final and initial state
G = c C + d D (a A + b B )
(3.18)
The chemical potential of the constituent "i" in the mixture can be written
i = i0 + RT ln ai
(3.19)
where ai is the activity of constituent "i" in the mixture and 0i is the chemical
potential of constituent "i" at a chosen standard state of unit activity. The activity
of pure solids or liquids is usually taken as unity at atmospheric pressure (1 bar or
105 Pa), while for gases the activity is unity when the partial pressure of the gas
constituent is 1 bar.
In solutions the activity is not unity, but may be equated with the mole
fraction, Xi, in ideal solutions:
(3.20)
ai = X i
3.6
3. Defect equilibria
In non-ideal solutions the activity is related to the mole fraction through the
activity coefficient, i,
(3.21)
ai = i X i
aCc a Dd
a Aa a Bb
(3.22)
(3.23)
equilibrium
(3.24)
K is termed the equilibrium constant and it relates the activities (or concentrations
under ideal conditions) of products and reactants when equilibrium has been
attained at a given temperature. G may as in Eq. 3.1 be expressed in terms of
the standard enthalpy change, H, and entropy change, S:
G 0 = H 0 TS 0 = RT ln K
(3.25)
H 0
S 0
H 0
) = exp(
) exp(
)
RT
R
RT
(3.26)
Thus, we have arrived at the same type of expression for the equilibrium constant
as we did in our intitial example from statistical thermodynamics. From that, we
know that the entropy term contained in K = exp(S/R) contains vibrational
terms only, and not the configurational entropy.
From the equation above we see that the temperature dependence of K is
given by
3.7
3. Defect equilibria
d ln K
H 0
=
d (1 / T )
R
(3.27)
and that a plot of lnK vs 1/T (vant Hoff plot) will give S/R as the intercept and
-H/R as the slope.
This general treatment is applicable to any ideal solutions whether these
are gaseous, liquid, or solid, and it applies also to defect chemistry as we have
examplified initially.
i = i0 + RT ln ai
(3.28)
The virtual potential differs from the true chemical potential by an undefined
constant which is incorporated in . For this reason the absolute value of the
virtual chemical potential can not be determined experimentally.
In real changes in crystals where the site ratio is maintained and
complementary structure elements (or building units) are created or annihilated,
the undefined constants in the virtual chemical potential of each separate structure
element cancel out, and the overall change is characterised by real changes in
3.8
3. Defect equilibria
Gibbs free energies. Real changes may thereby be described in terms of virtual
potentials. Krger et al. pointed out that a similar situation occurs in aqueous
electrolytes. In this case, arbitrary thermodynamic potentials are assigned to
separate ions, although it is impossible to determine experimentally the separate
potentials.
3. Defect equilibria
detailed understanding of defect interactions and how they influence defectdependent properties of inorganic compounds.
Despite its shortcomings the ideal solution approach will be used in this
book in treating defect equilibria. Defect interactions will generally be treated by
assuming the formation of associated defects as is, for instance, exemplified above
for the formation of singly charged metal vacancies through the association of an
electron hole with a doubly charged metal vacancy. Furthermore, the defect
equilibria will be treated using the law of mass action.
Concentrations of the structure elements may be expressed in many units. In
semiconductor physics it is common to express the concentrations in number per
cm3. This has the disadvantage, however, that when two compounds with different
molecular (packing) densities contain the same number of defects per cm3, the
number of defects per molecular unit is different in the two compounds, and vice
versa. Furthermore, the molecular unit size (packing density, or unit cell
parameters) of a compound changes as a function of its nonstoichiometry and
temperature, and in this case the unit of defects per cm3 does not unequivocally
reflect the relative defect concentration per molecular unit or atom site under
different conditions.
In the following, concentration will in most cases be expressed as the
number of defects or atoms per molecule or site, i.e. in molar or site fractions,
while we for some purposes need to use the number of defects per cm3. As
mentioned before, the value of the equilibrium constant depends on the units of
concentration that are employed, but it is a simple matter to convert the values of
the equilibrium constant from one system to another.
Examples of defect
stoichiometric metal oxides
equilibria
in
pure,
In this section we will see how we write defect reactions for selected simple
cases of combinations of defects, express the equilibrium constants, and combine
this with an electronutreality condition to express the concentration of dominating
defects as a function of temperature and activities. We will stay strictly with
statistical thermodynamics, such that entropy changes in the defect reaction can be
interpreted strictly as a change in the vibrational entropy. Approximations will be
implemented in order to visualise the simplicity f the method.
Schottky defects
As illustration let us consider the defect equilibrium for Schottky defects in
the oxide MO and where the metal and oxygen vacancies are doubly charged. The
defect equation for their formation is given in Chapter 2 as
0 = vO + v M//
(3.29)
3.10
3. Defect equilibria
The corresponding defect equilibrium may at low defect concentrations be written
K S = X v X v // =
O
[vO ] [v M// ]
[O] [M]
(3.30)
where KS is the equilibrium constant for a Schottky defect pair, X represent site
fractions, and square brackets express concentrations of species. [O] and [M]
represent the concentration of oxygen and metal sites, respectively.
If we let square brackets express concentration in mole fraction (i.e. the
number of the species per mole of MO), then [O] = [M] = 1, and then
K S = [vO ][v M// ]
(3.31)
G S
S
H S
) = exp( S ) exp()
RT
R
RT
(3.32)
It may be useful to repeat that the entropy change , SS, does not include the
configurational entropy changes associated with the formation of the defects, but
mainly involves the change in vibrational entropy.
When the Schottky pair dominates the defect structure, we may see from
the reaction equation or from the electroneutrality condition that in MO the
concentrations of the metal and oxygen vacancies are equal:
[vO ] = [v M// ]
(3.33)
SS
H S
) exp()
2R
2 RT
(3.34)
Under these conditions the defect concentrations are only dependent on the
temperature; they are independent of the activities of the components M and O2.
The slope in a vant Hoff plot of lnK vs 1/T would be -HS/2R. The factor 2
results from the square root of the equilibrium constant, in turn resulting from the
fact that there are two defects created, in comparison with the one we got in the
element used as example earlier.
It may be useful to note that the procedure we have applied is to write
down the defect reaction and the mass action law expression for its equilibrium
constant, and the prevailing electronuetrality condition. Thus, we have two
independent equations for the two unknowns (provided that either KS or all of
SS, HS, and T are known).
3.11
3. Defect equilibria
(3.35)
It may be noted that a vacant interstitial site has been included to keep track of the
sites where the interstitials are going. The corresponding defect equilibrium is
written as
KF =
X M X v //
i
X M x X vx
M
[ M i ] [v M// ]
//
[i] [ M ] [ M i ][v M ]
=
=
[ M Mx ] [vix ] [ M Mx ][vix ]
[ M ] [i ]
(3.36)
If there is one interstitial site per MO, and if the defect concentrations are small
such that [ M Mx ] = [vix ] = 1 , we obtain the simpler
K F = [ M i ][v M// ]
(3.37)
S F
H F
)
) exp(2R
2 RT
(3.38)
(3.39)
3.12
3. Defect equilibria
The equilibrium constant may then be expressed as site fractions, which in our
case immediately simplifies into an expression containing only volume or molar
concentrations:
Ki =
[ M M ][M M/ ]
[ M Mx ] 2
(3.40)
(3.41)
S i
H i
) exp(
)
2R
2RT
(3.42)
(3.43)
(3.44)
K i/ = exp(
Eg
RT
) = K 0,/ i exp(
H i
)
RT
(3.45)
Eg
2 RT
) = (K 0,/ i )1 / 2 exp(
H i
)
2 RT
(3.46)
3.13
3. Defect equilibria
X v X M2 /
O
X O x X M2 x
O
pO1 /22 =
[vO ][M M/ ]2
[OOx ][M Mx ] 2
(3.47)
p O1 /22
(3.48)
If these oxygen vacancies and the compensating electronic valence defects are the
predominating defects in the oxygen deficient oxide, the principle of
electroneutrality requires that
[ M M/ ] = 2[vO ]
(3.49)
S vO
H vO 1 / 6
) p O2
) exp(3R
3RT
(3.50)
This shows that the concentrations of the dominating defects (defect electrons and
oxygen vacancies) increase with decreasing oxygen partial pressure. It
furthermore shows the relation between defect concentrations and entropy and
enthalpy changes of the defect reaction. The factor 3 that enters results from the
formation of 3 defects in the defect reaction.
Now, let us do the same treatment, but for an oxide where defect electrons
are delocalised in the conduction band. The formation of oxygen vacancies is then
written:
OOx = vO + 2e / + 12 O2 ( g )
(3.51)
3.14
3. Defect equilibria
Since delocalised electrons in a band cannot be treated by classical statistics, we
abandon the site fraction principle, and express all activities simply as
concentrations, and denote the new equilibrium constant with a prime4:
K
/
vO
[vO ] 2 1 / 2
= x n pO2
[OO ]
(3.52)
If these oxygen vacancies and the compensating electrons are the predominating
defects in the oxygen deficient oxide, the principle of electroneutrality requires
that
n = 2[vO ]
(3.53)
H vO 1 / 6
) pO2
3RT
(3.54)
4
This text introduces the use of a prime in the equilibrium constant to denote cases where
statistical thermodynamics has not been used throughout, and where the entropy of the reaction
thus is not well defined. This notation is new and not in common use, as is neither the distinction it
notifies. One may therefore expect also that this text relaxes its stringency in this respect in later
chapters, where the distinction is less in focus.
3.15
3. Defect equilibria
If we want to investigate the more complex situation of involving neutral
and partially ionised oxygen vacancies and the excitation of the trapped electrons
(described in Chapter 2) we start by writing down the stepwise reactions and their
corresponding equilibria:
OOx = vOx + 12 O2 ( g )
K vOx
(3.55)
[vOx ] 1 / 2
= x p O2
[OO ]
(3.56)
vOx = vO + e /
/
K vO
1 =
(3.57)
[vO ]
n
[vOx ]
(3.58)
vO = vO + e /
K
/
vO 2
(3.59)
[vO ]
= n
[v O ]
(3.60)
/
/
where K vOx , K vO
1 , and K vO 2 are the respective equilibrium constants. Now the
principle of electroneutrality requires that
n = [vO ] + 2[vO ]
(3.61)
The concentrations of the electrons and the neutral, singly, and doubly charged
oxygen vacancies are related through the equations above, and by combination
expressions for each of the defects may be obtained. The electron concentration is
given by
/
/
1 / 2
n 3 = K vOx K vO
1 ( 2 K vO 2 + n) p O2
(3.62)
(3.63)
Under these conditions the concentrations of electrons and oxygen vacancies are
relatively small, and under these conditions the doubly charged vacancies are the
predominating oxygen vacancies, and the concentrations of electrons and oxygen
vacancies are proportional to pO12 / 6 . Of course, this situation corresponds to the
simple case presented initially for the doubly charged vacancies, and obviously
/
/
/
K vOx K vO
2 K vO 2 = K vO .
3.16
3. Defect equilibria
If, on the other hand, n >> 2 K vO/ 2 , then
/
1/ 2
n = [vO ] = ( K vOx K vO
pO12 / 4
1)
(3.64)
Under these conditions the concentration of electrons is relatively large, and the
predominant oxygen vacancies are then singly charged. Furthermore, the
concentrations of electrons and oxygen vacancies are then proportional to p O12 / 4 .
A general tendency similar to that of oxygen deficient oxides applies to
metal deficient oxides; in the oxide M1-yO the metal vacancies are doubly charged
at very small deviations from stoichiometry and tend to become singly charged
with increasing nonstoichiometry.
It may be noted that the neutral oxygen vacancies are not affected by
charged defects, nor do they affect the electroneutrality.
OOx + M Mx = M i + 3e / + 34 O2 ( g )
(3.65)
(3.66)
where KMi is the equilibrium constant. If these defects are the predominating
ones, and we as before assume small defect concentrations, the electroneutrality
condition and the oxygen pressure dependence of the interstitial metal ions and the
electrons is
/
3[ M i ] = n = (3K Mi
)1 / 4 pO32 / 16
(3.67)
Thus, also in this case the concentration of the defects increases with decreasing
oxygen activity and is proportional to p O-3/16
. This oxygen pressure dependence is
2
different from that for formation of singly as well as doubly charged oxygen
vacancies. In such a case it would thus in principle be possible to decide from
measurements of electron concentration (e.g. via electrical conductivity) or of
nonstoichiometry as a function of oxygen activity whether the predominating
defects are triply charged interstitial metal ions or oxygen vacancies with one or
two charges.
However, had the metal interstitials predominantly been doubly charged in
the oxide M2+yO3, a p O-1/4
dependence would have been the result, and additional
2
3.17
3. Defect equilibria
studies of defect-dependent properties (e.g. self-diffusion of metal or oxygen)
would be needed to distinguish this situation from that of singly charged oxygen
vacancies.
It may be noted that it is not only the absolute number of charges on a
defect that determines the p O 2 dependence, but also the difference between the
actual charge and the charge given by the nominal valence of the atoms involved.
Thus, in the oxide MO (or M1+yO) doubly charged metal interstitials will give a
p O-1/6
dependence of the defect concentrations (as opposed to the case above) and
2
in this case be indistinguishable from doubly charged oxygen vacancies.
Simultaneous presence
interstitial metal ions
of
oxygen
vacancies
and
In an oxide MaOb where the ratio of metal to oxygen is larger than the
stoichiometric ratio a:b it is a priori difficult to predict whether interstitial metal
ions or oxygen vacancies predominate. In principle both types of defects may be
important, at least in certain regions of nonstoichiometry. In the following, let us
consider a case where it may be necessary to take into account the simultaneous
presence of interstitial metal ions and oxygen vacancies.
Consider an oxide with a stoichiometric composition MO2. Let us further
assume for the sake of illustration that when the oxide is nonstoichiometric the
important point defects are doubly charged oxygen vacancies and doubly charged
interstitial metal ions. (It mat be noted that the metal interstitials are not fully
ionised in this case.) The composition of the nonstoichiometric oxide may
accordingly be written M1+yO2-x.
The defect equations for the formation of these two types of defects may be
written
OOx = vO + 2e / + 12 O2 ( g )
(3.68)
2OOx + M Mx = M i + 2e / + O2 ( g )
(3.69)
(3.70)
2
x -2
x -1
K Mi
2 = [ M i ]n p O2 [OO ] [ M M ]
(3.71)
n = 2[vO ] + 2[ M i ]
Two limiting conditions may be considered:
3.18
3. Defect equilibria
When [vO ] >> [ M i ] then
/
n = 2[vO ] = (2 K vO
)1 / 3 pO12 / 6
(3.73)
as obtained also earlier for the same conditions. By inserting this relationship for n
into the equilibrium constant for formation of the minority defects, metal
interstitials, we obtain
/
/
2 / 3 2 / 3
[ M i ] = K Mi
pO2
2 ( 2 K vO )
(3.74)
(3.75)
By inserting this relationship for n into the equilibrium constant for formation of
the minority defect, oxygen vacancies, we obtain
/
/
2 / 3 1 / 6
[vO ] = K vO
(2 K Mi
p O2
2)
(3.76)
3.19
3. Defect equilibria
It may be noted that under these conditions the concentration of oxygen
vacancies increases with the oxygen pressure, as shown in the figure.
The overall situation may be illustrated schematically as shown in Fig.3.3,
in which the concentrations of the three defects are plotted as a function of the
partial pressure of oxygen. In this diagram it is assumed that [vO ] is 105 times
larger than [ M i ] at 1 atm O2, and with this assumption it is seen that the
interstitial cations will predominate at a partial pressure below about 10-17 atm.
Many of the high melting oxides are stable over large oxygen pressure ranges and
it is not improbable that this hypothetical situation, with oxygen vacancies
predominating at high partial pressures of oxygen (and small deviations from
stoichiometry) and interstitial cations predominating at low partial pressures
(relatively large deviations from stoichiometry) may generally apply to oxygen
deficient oxides at elevated temperatures. In actual cases more than one ionisation
state of the defects may of course also have to be considered.
In this section we have employed and illustrated a very important sequence
of actions often used to solve ideal situations in defect chemistry, but also
implicitly used in traditional aqueous and other chemistry. We write the list down
in a numbered list for clarity:
1.
2.
3.
4.
5.
6.
Metal-deficient oxide
3.20
3. Defect equilibria
The reactions and defect equilibria for the formation of single, unassociated
neutral metal vacancies and the subsequent excitation of electron holes in MO
may for small defect concentrations be written
1
2
O2 ( g )= v Mx + OOx
(3.77)
(3.78)
v Mx = v M/ + h
(3.79)
[v M/ ]
p
[v Mx ]
/
K vM
1 =
(3.80)
v M/ = v M// + h
/
K vM
2 =
(3.81)
[v M// ]
p
[v M/ ]
(3.82)
/
/
where K vMx , K vM
1 , and K vM 2 are the equilibrium constants for the respective
defect equilibria.
When the metal vacancies and their complementary electron holes are the
predominating defects, the electroneutrality condition reads
p = [v M/ ] + 2[v M// ]
(3.83)
(3.84)
(3.85)
/
At high concentrations of electron holes ( p >> K vM
2 ) we obtain
3.21
3. Defect equilibria
/
1/ 2 1/ 4
p = ( K vMx K vM
p O2
1)
(3.86)
(3.87)
(3.88)
K i/ = np
(3.89)
K AF = [Oi// ][vO ]
(3.90)
3.22
3. Defect equilibria
In these equilibria we have assumed small defect concentrations such that the
concentrations of normal lattice sites and empty interstitial sites have been
assumed constant and equal to unity and thus omitted from the expressions.
It should be noted that the defect equilibria are interrelated, and through a
2
/
combination of the equations it may be shown that K i/ K AF
= K vO
K Oi/ . Thus, only
three out of the four equilibria are sufficient to describe the defect structure of the
oxide.
The full electroneutrality condition is given by
2[vO ] + p = 2[Oi// ] + n
(3.91)
Oxygen deficit.
At large oxygen deficit the following approximation may be made
n = 2[vO ] >> 2[Oi// ], p
(3.92)
We first insert this into the appropriate equilibrium to find the concentrations of
the dominating defects in the usual manner:
/
n = 2[vO ] = (2 K vO
)1 / 3 pO12 / 6
(3.93)
and then insert this into other equilibria to find the concentrations of the minority
defects. By inserting the expression for the concentration of vacancies into the
anion Frenkel equilibrium, we obtain for the concentration of interstitials:
/
[Oi ] = 2 2 / 3 K AF ( K vO
) 1 / 3 pO1 /26
(3.94)
By inserting the expression for the concentration of electrons into the intrinsic
electronic equilibrium, we obtain for the concentration of holes:
/
p = K i/ (2 K vO
) 1 / 3 p O1 /26
(3.95)
Oxygen excess
For relatively large excess oxygen, that is, when
p = 2[Oi// ] >> 2[vO ], n
(3.96)
3.23
3. Defect equilibria
we may in a manner analogous to the preceding case, derive the following
relations:
p = 2[Oi// ] = (2 K Oi/ )1 / 3 p O1 /26
(3.97)
(3.98)
n = K i/ (2 K Oi/ ) 1 / 3 p O12 / 6
(3.99)
Stoichiometric condition
At or close to stoichiometry two alternative limiting conditions must be
considered, namely dominance by intrinsic electronic ionisation or by anion
Frenkel disorder.
If intrinsic ionisation of electrons predominates, and thus
p = n = ( K i/ )1 / 2 >> 2[Oi// ],2[vO ]
(3.100)
The concentrations of electrons and electron holes (n and p) are then independent
of oxygen pressure. The point defect concentrations are obtained by insertion of
the expression for the concentrations of the electronic defects into the appropriate
equilibria, and we obtain:
[vO ] =
/
K vO
pO12 / 2
/
Ki
(3.101)
(3.102)
3.24
3. Defect equilibria
Figure 3-4. Schematic presentation of oxygen point defects and electronic defects as a
function of oxygen pressure in an oxide which depending on the partial pressure of oxygen may
have an excess or deficit of oxygen. Intrinsic electronic equilibrium is assumed to predominate at
stoichiometric composition.
[O ] = [v ] = K
1/2
AF
p n ( K i/ )1 / 2
>> , =
2 2
2
(3.103)
then [Oi// ] and [vO ] are independent of the partial pressure of oxygen, while the
concentrations of electronic defects are given by
/
1 / 4 1 / 4
n = ( K vO
)1 / 2 K AF
pO2
(3.104)
1 / 4 1 / 4
p = ( K Oi/ )1 / 2 K AF
p O2
(3.105)
3.25
3. Defect equilibria
Figure 3-5. Schematic presentation of the concentration of oxygen point defects and
electronic defects as a function of oxygen pressure. Oxygen defects are assumed to predominate at
stoichiometric composition
From the diagrams we can conclude that at low oxygen activities the oxide
has oxygen deficiency and will be an electronic n-type conductor, because the
mobility of the electrons is always much higher than that of oxygen vacancies. At
high oxygen partial pressures the oxide will correspondingly have oxygen excess
and be a p-type electronic conductor. At intermediate oxygen activities it may be
be a mixed n- and p-type electronic conductor in the case of intrinsic electronic
disorder, while it may exhibit ionic or mixed ionic/electronic conduction in the
case of anion Frenkel disorder.
Similar diagrams and analyses may be made for many other combinations of
non-stoichiometric and stoichiometric defect situations of pure (undoped) oxides
and other ionic compounds.
Summary
Chemical reactions for defects can be formulated and treated using the
mass-action law. As for other chemical reactions, equilibrium constants can be
defined in terms of the activities of the defects and other species. Under the
normal constrictions we can approximate activities with concentrations of defects
and partial pressures of gases. The equilibrium constants can also be expressed in
3.26
3. Defect equilibria
terms of the Gibbs free energy change of the reaction, which in turn is expressed
by the entropy and enthalpy changes.
We have seen that the standard entropy change of the standard Gibbs energy
change contains only vibrational terms provided that the configurational terms are
properly handled in expressing the equilibrium conditions. This is often possible
in defect idealised defect chemistry since we can apply classical statistical
thermodynamics. However, for non-classical defects such as delocalised electrons
and holes this approach is less meaningful, and the entropy change may then be
interpreted in a less straightforward manner. We have introduced the use of a
prime (as in K/) to denote equilibrium constants that are not expressed according
to statistical thermodynamics and where the entropy change in the Gibbs energy
change may have other than vibrationan contributions.
The mass action expressions can be combined with the full or limited cases
of the electroneutrality condition to obtain exact or approximate (limiting case)
expressions for the concentration of defects. Such concentrations are typically a
function of the oxygen partial pressure, temperature, and doping. It is common to
illustrate defect structures for oxides by plotting log defect concentrations vs log
p O2 (Brouwer diagrams) or ln or log defect concentrations vs 1/T. We have
shown these principles and techniques through examples of intrinsic ionic and
electronic disorder, various types of nonstoichiometry and variable ionisation of
point defects.
We have restricted our treatment to simple cases and made simplifications
where possible. Assumption of small defect concentrations have allowed the
assumptions that the concentrations of normal lattice atoms and empty interstitials
sites are constant, with activities equal to unity. Larger defect concentrations can
to a first approximation be taken into account by including mass and site balances
into the expressions used to solve the defect structure.)
Literature
Krger, F.A. (1964) The Chemistry of Imperfect Crystals, North-Holland,
Amsterdam, and Wiley, New York.
Krger, F.A., Stieltjes, F.H. and Vink, H.J. (1959), Philips Res. Rept. 14,
557.
Exercises
1. Do the insertion of Stirlings approximation and the derivation that leads
up to the relation between equilibrium constant and entropy and enthalpy
changes for the elemental reaction EE = vE + EE (Eqs. 3.3 and onwards).
2. Write the reaction for formation of Schottky defects in MO2 and find an
expression for the defect concentrations as a function of the equilibrium
constant and thermodynamic parameters and temperature. Note in
3.27
3. Defect equilibria
particular how the solution deviates from the one obtained for the oxide
MO treated in the text.
3. In the case of cation Frenkel defects in MO we assumed in the text that we
had one interstitial site per MO. Derive the expression for the defect
concentrations if the structure consists of fcc close-packed O ions, with M
ions on each octahedral hole, and with interstitial sites on all tetrahedral
holes (Hint: consult Chapter 1).
4. For a defect situation in an oxide dominated by doubly charged oxygen
vacancies and electrons, sketch the vant Hoff plot (Logarithm of defect
concentration vs 1/T) and a double-logarithmic plot of defect
concentrations vs pO2 (Brouwer diagram).
5. Consider a metal oxide MO1-x dominated by singly and doubly charged
oxygen vacancies. Suggest a condition that expresses the changeover from
dominance of one to dominance of the other. Next, find an expression for
the pO2 (as a function of equilibrium constants) of this transition point.
6. Consider a metal oxide M2-xO3 dominated by doubly and triply charged
metal vacancies. Suggest a condition that expresses the changeover from
dominance of the doubly charged to dominance of the triply charged.
Next, find an expression for the pO2 (as a function of equilibrium
constants) of this transition point.
7. Find a general expression for m in the [e/] pO21/m relationship for oxides
MOb with non-stoichiometry dominated by fully ionised metal interstitials.
(Hint: Find m as a function of b). Use this to find m in the case of M2O,
MO, M2O3, MO2, M2O5 and MO3.
8. Find a general expression for m (as a function of the valency z of the
cation) in the [h.] pO21/m relationship for binary oxides with nonstoichiometry dominated by fully ionised metal vacancies. Use this to find
m in the case of M2O3.
9. Sketch a full Brouwer diagram (log defect concentrations vs log pO2) for
an oxide M2O3 dominated by fully ionised oxygen and metal vacancies at
under- and overstoichiometry, respectively. Assume that Schottky defects
predominate close to stoichiometric conditions. (Hint: the main goal is to
obtain and illustrate the pO2-dependencies. Use the rules we listed for such
constructions of Brouwer diagrams.)
10. For the defect situation derived in the preceding exercise for M2O3 do the
following analyses:
i.
find a condition and the expression for the pO2 when the oxide is
fully stoichiometric.
ii.
iii.
3. Defect equilibria
iv.
11. Sketch a full Brouwer diagram (log defect concentrations vs log pO2) for
an oxide MO2 dominated by fully ionised oxygen and metal vacancies at
under- and overstoichiometry, respectively. Assume that intrinsic
electronic equilibrium predominates close to stoichiometric conditions.
(Hint: the main goal is to obtain and illustrate the pO2-dependencies. Use
the rules we listed for such constructions of Brouwer diagrams.)
12. For the defect situation derived in the preceding exercise for MO2 do the
following analyses:
i.
find a condition and the expression for the pO2 when the oxide is
fully stoichiometric.
ii.
iii.
13. Sketch a full Brouwer diagram (log defect concentrations vs log pO2) for
an oxide ABO3 dominated by fully ionised oxygen and metal vacancies at
under- and overstoichiometry, respectively. Assume that Schottky defects
predominate close to stoichiometric conditions. You may assume that both
cations are trivalent, but discuss also the effect it would have if A was
divalent and B tetravalent. (Hint: the main goal is in any case to obtain and
illustrate the pO2-dependencies. Use the rules we listed for such
constructions of Brouwer diagrams.)
14. Sketch a full Brouwer diagram (log defect concentrations vs log pO2) for
an oxide AB2O4 dominated by fully ionised metal interstitials and
vacancies at under- and overstoichiometry, respectively. Assume that
Frenkel defects predominate close to stoichiometric conditions. You may
assume that A is divalent and B trivalent. (Hint: the main goal is to obtain
and illustrate the pO2-dependencies. Use the rules we listed for such
constructions of Brouwer diagrams.)
15. In the text we considered an oxide MO2 that had oxygen deficiency by
both fully charged oxygen vacancies and doubly charged metal
interstitials. The figure illustrating the case was drawn assuming that the
ratio between the concentrations of the two defects was 105 at 1 atm. pO2.
Calculate the exact pO2 at which the concentrations are equal.
16. In the text we considered an oxide MO2 that had oxygen deficiency by
both fully charged oxygen vacancies and doubly charged metal
interstitials. The figure illustrating the case is not strictly a Brouwer
diagram, but rather a diagram showing the exact solution of the defect
3.29
3. Defect equilibria
structure, thus also correctly describing the transition region. The full
solution can be reached in several ways:
i.
ii.
iii.
iv.
A program that solves the mixture of different types of equations numerically has been
written by T. Norby. Compared to the other computerized approaches it executes very slowly.
3.30
4.1
4.2
(4.1)
The defect structure accordingly has two domains given by the limiting cases of
this expression:
If the foreign cation is a minority species, i.e., [ Mf M ] << [v X ] then
[v X ] = [v M/ ] , and from the Schottky equilibrium
K S = [v X ][v M/ ]
(4.2)
we obtain
[v X ] = [v M/ ] = K 1/2
S
(4.3)
(4.4)
,i.e., the metal vacancy concentration is fixed by the concentration of the foreign
cation. The concentration of the minority species, [v X ] , can be found by insertion
into KS:
[v X ] =
KS
[ Mf M ]
(4.5)
4.3
Frenkel equilibria
The effects of higher valent cation impurities on Frenkel equilibria are
analogous to those described for the Schottky equilibria. When divalent Mf M
impurities/dopants are dissolved in MX with predominating cation Frenkel defect
pairs ( v M/ and M i ) the electroneutrality condition is given by
[ Mf M ] + [ M i ] = [v M/ ]
(4.6)
4.4
(4.7)
[vO ] 2 1 / 2
n pO2
[OOx ]
(4.8)
/
=
K vO
(4.9)
(4.10)
4.5
Figure 4-2. Brouwer plot of the concentrations of defects as a function of oxygen partial
pressure in an oxygen deficient oxide predominantly containing doubly charged oxygen vacancies,
showing the effects of a constant concentration of lower valent cation dopants,
[ Mf M/ ] .
(4.11)
The effects of the foreign cations are evaluated by combining this with the
equations for the defect equilibria for the singly charged metal ion vacancies. As a
result we find that the lower valent cations will increase p and decrease [v M/ ] . The
variations in p and [v M/ ] as a function of the partial pressure of oxygen in such a
case are illustrated schematically in the Brouwer diagram below. At low oxygen
activities interstitial metal ions or oxygen vacancies may become important,
depending on the equilibrium constants for the Frenkel and Schottky defect
equilibria. In the figure it is assumed that singly charged metal interstitials become
the important native point defects at low oxygen activities.
Figure 4-3. Effect of a constant concentration of foreign lower valent cations on the
concentration of point and electronic defects as a function of oxygen pressure in a metal deficient
oxide predominantly containing singly charged cation vacancies.
4.7
(4.12)
In this case [v M/ ] will increase and p will decrease with increasing concentration
of Mf M . The variations in p and [v M/ ] as a function of the oxygen partial pressure
is illustrated schematically in the following figure, and the variation in the
concentration of electrons (n) as minority defects is also included.
Figure 4-4. Brouwer diagram of the effects of a higher valent cation impurity/dopant on the
concentration of point and electronic defects as a function of oxygen pressure in a metal deficient
oxide predominantly containing singly charged metal vacancies.
(4.13)
We assume that the amount of MlO added is well below the solubility so
that [ Ml M/ ] = constant. Fig. 4.5 shows a Brouwer diagram of the defect situation
as a function of oxygen partial pressure when the level of aliovalent dopant is
higher than the level of intrinsic disorder. At the lowest pO2 the oxide is oxygendeficient and oxygen vacancies and electrons predominate. As these defects
decrease with increasing p O2 we hit the level of the acceptor dopant. From here
on the acceptor will be compensated by a constant concentration of the positive
defect, i.e., the oxygen vacancies. This is analogous to the situation we obtained
when we earlier considered this case in particular.
In this situation the concentration of defect electrons decreases and that of
holes increases with increasing p . The two terms cross at n = p = K i1 / 2 , and
assuming that Ki is independent of the doping, this is necessarily the same level as
that of the intrinsic electronic disorder in the pure, undoped material. At this point
also the doped material may be said to be stoichiometric; it contains oxygen
vacancies exactly matching the presence of the lower-valent cation (acceptor): If
we consider the doped material to consist of M2O3 and MlO, both these
constituents are stoichiometric, and the additional presence of reduced and
oxidised states (n and p) is effectively zero at this point. The doped oxide can be
said to be stoichiometric with respect to the valence of its constituents.
By further increasing p O2 the concentration of holes will eventually become
dominating, taking over the dominance from the oxygen vacancies in terms of
compensating the charge of the acceptors. We hereby create higher oxidation
states (holes) and fill up the oxygen sublattice (removing oxygen vacancies).
During the increase in p O2 the concentration of metal vacancies increases and
eventually balances the concentration of oxygen vacancies. This level is given by
the Schottky equilibrium constant (KS) and represents another stoichiometric
point, this time with respect to the crystal structure of the host oxide (M2O3).
The situation encountered in this region is analogous to that depicted in Fig.
4.3, and with further increase in p O2 the charge due to the increasing
concentration of metal vacancies will become more important than that of the
acceptor level, and the oxide enters into a situation of dominance by metal
deficiency.
4.9
Figure 4-5. Brouwer diagram of the effects of lower valent dopant on the defect structure of
M2O3..
4.10
Figure 4-6. Brouwer diagram of the effects of higher valent dopant on the defect structure
of M2O3.
4.11
(4.14)
By taking the activity of the excess MlO(s) phase as well as that of oxygen ions as
unity, the equilibrium expression can be written
[ Ml M/ ]2 p 2 [OOx ]3
Kp =
= [ Ml M/ ]2 p 2 p O12/ 2
2
1/ 2
a MlO(s) p O 2
(4.15)
The important difference from the previous treatments is that now we treat this
dissolution as an equilibrium. When the compensating defects are electron holes,
the simplified electroneutrality condition can be expressed by [ Ml M/ ] = p , and by
insertion into the equilibrium expression we obtain
1/ 8
[ Ml M/ ] = p = K 1/4
p p O2
(4.16)
Thus, the solubility, i.e., [ Ml M/ ] , increases with p O2 until all MlO is eventually
dissolved. Thereafter the behaviour will become as with constant acceptor
concentrations (all dissolved or frozen-in), described in Fig. 4.5.
When taking into account the defect equilibrium for the formation of oxygen
vacancies and the intrinsic electronic equilibrium (np = Ki), the concentration of
oxygen vacancies is under this situation seen to be proportional to p O12/ 4 . At lower
oxygen activities the doubly charged oxygen vacancies will thus eventually
become the charge compensating defects at the expense of electron holes. In this
case it may be more relevant to describe the dissolution of MlO by the following
defect reaction and equilibrium expression:
2MlO(s) = 2 Ml M/ + vO + 2OOx
(4.17)
4.12
K v = [ Ml M/ ]2 [vO ]
(4.18)
The electroneutrality condition and its combination with the equilibrium then
yield:
[ Ml M/ ] = 2[vO ] = (2K v )1 / 3
(4.19)
4.13
Ke =
= [ MhM ]2 n 2 p1O/22
(4.20)
(4.21)
(4.22)
4.14
Figure 4-8. Brouwer plot of the effect of the dopant MhO2 in M2O3 in amounts exceeding
the solubility limit.
All in all, as a general rule the solubility of lower or higher valent oxides in
a parent oxide will be dependent on the oxygen activity when the compensating
defects are electronic defects, while they will be independent of the oxygen
activity when the compensating defects are fully ionised point defects. It may also
be noted that the temperature dependencies of the solubility will also be different
in the two cases.
Partially ionised point defects act as combinations of point and electronic
defects and generally gives oxygen activity dependent solubilities of aliovalent
dopants.
In real cases it may be difficult to reach equilibrium if the diffusion of the
cations to and from defect sinks (dislocations, grain boundaries a.o.) is slow.
Another point worth mentioning is that the excess second phase may well be a
compound of the parent and dopant oxide. For instance, MgO is normally present
in excess of the solubility in so-called high-purity alumina (Al2O3) and the
second phase in that case is spinel MgAl2O4. However, this does not alter the
oxygen activity dependencies; the ternary compound serves as a source/sink for
the dopant, with constant activity, and the defect chemical treatment remains
essentially the same.
(4.23)
(4.24)
K=
[OH O ] 2 n 2 p1O/22
(4.25)
[O Ox ]2 p H 2O
(4.26)
Let us keep the partial pressure of oxygen constant and vary the partial pressure of
the water vapour. Protons are then dissolved in the oxide following Eq. 4.24.
However, the concentration of electrons is given by Eq. 4.26 and is independent of
p H 2O . Insertion into Eq. 4.25 gives that the proton concentration is proportional to
p1H/22O (at constant p O 2 ). This is illustrated in the Brouwer diagram in the figure
below.
Figure 4-9. Brouwer plot of effects of water vapour on defect concentrations in oxygen
deficient M2O3-d.
4.17
At sufficiently high partial pressure of the water vapour, the dissolved protons
become the predominating point defects with positive effective charge and the
electroneutrality condition becomes
(4.27)
n = [OH O ]
By combination with Eq. 4.25 we obtain that the concentrations of protons and
electrons are then proportional to p1H/24O :
n = [OH O ] = K 1 / 4 p1H/24O pO12 / 8
(4.28)
This situation is part of the figure above. The figure also shows that the
concentration of oxygen vacancies decreases with the water vapour pressure,
illustrating that in general all charged native minority defects becoming dependent
on water vapour pressure when protons are dominating defects.
4.18
Figure 4-10. Brouwer plot of the effect of water vapour (at constant oxygen pressure) on
defect concentrations in acceptor-doped, oxygen deficient M2O3.
(4.29)
(4.30)
When, as in the last example, protons are dissolved at the cost of oxygen
vacancies, the reaction may be written
H 2 O(g) + v O + 2O Ox = 2OH O
(4.31)
i.e. water vapour dissolves to fill the oxygen vacancies and replace their positive
charge with protons. Reaction 4.31 has been studied for a number of oxides, and it
is found that it always involves a negative enthalpy change such that oxygen
vacancies and protons are dominant at high and low temperatures, respectively.
For some oxides, the enthalpy is sufficiently negative that protons remain the
dominant defects up above 1000C in practically dry atmospheres and to even
higher temperatures under wet conditions. Other oxides have only small negative
enthalpies and have dominant protons up to only moderately high temperatures,
sometimes sufficiently low to kinetically prevent the reaction with water (Eq.
4.31) to reach equilibrium
The relative dominance of protons vs oxygen vacancies also varies with the
acceptor doping level: The higher charge of the vacancy relative to the proton
makes the former relatively more dominant at higher acceptor levels, while
protons are relatively more dominant at moderate acceptor levels.
Let us also briefly recollect the variable solubility of aliovalent dopants
treated previously in this chapter. While the acceptors are compensated by oxygen
vacancies, their solubility does not vary with anything else than the temperature,
as shown before. However, the dissolution of acceptors by simultaneous
dissolution of protons,
2MlO(s) + H 2 O(g) = 2 Ml M/ + 2OH O + O Ox
(4.32)
4.20
(4.33)
Thus, water vapour may in some cases play a role in the synthesis of solid
solutions containing acceptors.
By inspection of Fig. 4.10 it is evident that at sufficiently high water vapour
levels, the concentration of metal vacancies will become significant and
eventually dominant; the oxide attains a metal deficiency (or, in other cases,
oxygen excess), compensated by protons. Such defect structures can be considered
to be the first step towards the phase limit of the oxyhydroxide or hydroxide.
(4.34)
If all the species are randomly distributed, the equilibrium can be written
Ka =
[(v O Ml M ) ]
[v O ][Ml /M ]
(4.35)
4.21
Summary
Impurities or dopants may dissolve interstitially or substitutionally, and the
defects formed are often easily ionised. Depending on whether they accept or
donate electrons they become acceptors and donors. In sufficient quantities these
may dominate the defect structure of the oxide. Thus, aliovalent cations
substituting host cations are often introduced in large quantities in oxides in order
to give high concentrations of ionic defects (for solid electrolytes) or electronic
defects (for electrodes etc.) Foreign species may be treated like other defects, but
are often assumed to be present in a fixed concentration due to the limited
availability of the species (full solubility) or because the concentration is frozenin. However, we have shown that the concentration of foreign species under
certain conditions may change as a function of temperature and for some defect
structures as a function also of p O 2 . Protons dissolved from water vapour is a
special case which often has a large influence on the properties of oxides. They
form hydroxide groups with the oxygen ions of the oxide. These are effectively
positively charged and through them the defect structure may become a function
of the water vapour partial pressure.
Exercises
1. Write the doping reaction for the doping of the Schottky-dominated metal
halide MX by the higher valent metal halide MhX2. Find the analytical
solution expressing the concentration of metal vacancies as a funtion of
doping concentration. Sketch also the temperature-dependency of all the
defects involved. (Hint: This is the initial illustrative case in the text use
the equilibria from that treatment. Remember that the doping reaction may
4.22
4.23
4.24
5 Diffusion
5. Diffusion
Introduction
Numerous chemical reactions or micro-structural changes in solids take
place through solid state diffusion, i.e. the movement and transport of atoms in
solid phases. In crystalline solids, the diffusion takes place because of the
presence of defects. Point defects, e.g. vacancies and interstitial ions, are
responsible for lattice diffusion. Diffusion also takes place along line and surface
defects which include grain boundaries, dislocations, inner and outer surfaces, etc.
As diffusion along linear, planar and surface defects is generally faster than in the
lattice, they are also termed high diffusivity or easy diffusion paths. Another
frequently used term is short circuit diffusion.
The relative contribution of the different types of diffusion in oxides and
other inorganic compounds are functions of the temperature, partial pressures or
activities of the constituents of the compounds, the microstructure, grain size,
porosity etc. Grain boundary and dislocation diffusion generally have smaller
activation energies than lattice diffusion and as a result become increasingly
important the lower the temperature in solids with a given microstructure.
In the literature on diffusion and diffusion-controlled reactions or processes
one encounters many different terms that describe the diffusional behaviour under
different experimental conditions: tracer and self-diffusion of atoms and ions,
diffusion of defects, chemical diffusion, ambipolar diffusion, a.o. Many of these
are used for treating diffusion in compounds, and in the following chapters these
phenomena and terms will be described in more detail. Here we will start out with
a few simple phenomenological descriptions, and for simplicity we look only at
diffusion of neutral, independent particles.
Models of diffusion
It is a well known phenomenon that heat flows from hot to cold regions.
Such a flow of heat in a one-dimensional temperature gradient is described by
Fourier's law
j q =
dT
dx
(5.1)
5.1
5 Diffusion
where jq is the heat flux density, i.e., the flow of heat per unit area of the plane
dT
through which the heat traverses per second,
is the temperature gradient, and
dx
(kappa) is the thermal conductivity. It may be noted that the minus sign reflects
that the heat flows from high to low temperatures (downhill).
j particles = D
dc
dx
(5.2)
dc
the concentration gradient of the
dx
particles and D the diffusion coefficient. As in the equation for the heat flux, the
minus sign reflects that the particles flow from high to low concentration of
particles (downhill). This relation is named after A. Fick who first formulated this
relation. The particle flux and gradient is illustrated schematically in Fig.5.1.
Figure 5-1. Schematic illustration of Fick's first law. The negative of the particle
concentration gradient is the "driving force" of the diffusion.
jparticles represents the number of particles (or moles) crossing a unit area
or m2) per unit time (seconds). If the concentration of particles is expressed
in number of particles (or moles) per cm3 and the distance x in cm, the diffusion
coefficient has the dimension cm2s-1. In SI units the concentration is expressed in
number per m3 and the diffusion coefficient has the dimension m2s-1.
(cm2
dc
, may be
dx
considered to be an expression of the "driving force" for the particle flux (Eq.5.2).
The larger the concentration gradient, the larger the particle flux. When Ficks
5.2
5 Diffusion
first law is applied to uncharged (neutral) and independently diffusing particles, it
is valid in the sense that the coefficient is a constant (independent of concentration
and gradient). This may apply e.g. to dilute solutions of neutral defects in solids.
Such defects may be vacancies, interstitials, and impurities in metals. In ionic
solids they may comprise neutral interstitial impurities, homovalent substituents,
or isotopic species. The applicability or inapplicability of Ficks first law will
become clearer later in this chapter and in the forthcoming chapters.
(5.3)
vi = BiFi
Fi =
dPi
dx
(5.5)
where P is a potential. The negative sign is, as above, due to the fact that the
transport takes place from higher to lower values of P, see Fig.5.2.
5.3
5 Diffusion
Figure 5-2. Schematic illustration of Fick's first law for transport in a potential gradient.
When one combines Eqs. 5.3 - 5.5, the particle flux becomes
ji = ci vi = ci Bi
dPi
dx
(5.6)
ji = ci vi = ci Bi
d i
dx
(5.7)
(5.8)
i = i + kT lnai
ai =
ci
ci0
(5.9)
5 Diffusion
d i
d ln ci kT dci
= kT
=
dx
dx
ci dx
(5.10)
When one substitutes the expression for the chemical potential gradient in Eq.5.7,
the particle flux becomes
ji = ci Bi
d i
dc
= Bi kT i
dx
dx
(5.11)
(5.12)
Di = BikT
By combination with Eq.5.12, Eq.5.11 takes the form of Fick's first law as given
in Eq.5.2:
ji = Di
dci
dx
(5.2)
It is thus important to realise that Eq. 5.2 is fully valid only for ideal cases of
diffusion of neutral particles in chemical potential gradients only. If the particles
are charged, we need to take into account both the electrical potential gradent and
requirements to the combination of fluxes to maintain a given total current and
electroneutrality. These cases will be treated later on.
In the following, ideal conditions will be assumed and concentrations will
be used for activities of atoms, ions and various types of defects. But it should be
recalled that this always represents an approximation.
5 Diffusion
Fig.5.3. The two neighbouring planes under consideration are termed plane 1 and
2. The number of particles per unit area in plane 1 and 2 is termed n1 and n2,
respectively, and let us further assume that n1>n2. Consider further that the
particles in plane 1 and 2 may jump from one plane to a neighbouring plane at a
jump frequency . The particles in plane 1 have an equal probability of jumping to
plane 2 and to the neighbouring plane in the opposite direction. The total number
of particles jumping out of a unit area of plane 1 per unit time is equal to the
product of number of particles per unit area times the jump frequency: n1. As the
particles may jump in opposite directions, the number of particles jumping from
plane 1 to 2 is given by n1. Correspondingly, the number of particles jumping
from unit area of plane 2 to plane 1 is given by n2. The difference in the jump
rates is equal to the net flux density of particles:
j particles = 12 (n1 n 2 )
(5.13)
As all particles have the same jump frequency, there is a net flow of
particles from plane 1 to 2 because there are more atoms per unit area in plane 1
than in plane 2.
(5.14)
j particles = 12 (c1 c 2 ) s
5.6
5 Diffusion
dc
. The relation
dx
between c1 and c2 can then be expressed by seeing that a difference is equal to the
gradient times the length;
c 2 c1 = s
dc
dx
(5.16)
j particles = 12 s 2
dc
dx
(5.17)
j particles = 16 s 2
dc
dx
(5.18)
By comparing with Fick's first law, Eq. 5.2, it is seen that the diffusion coefficient
D in the three-dimensional case is given by
D = 16 s 2
(5.19)
If we consider a large number of jumps, n, which occurs during the time t, then
n
t
(5.20)
ns2 = 6Dt
5.7
5 Diffusion
The expression for the diffusion coefficient in Eqs.5.19 and 5.21 provides an
interesting qualitative description of what is going on during the diffusion. By
way of example, the diffusion coefficient for interstitial diffusion of oxygen atoms
in niobium metal at 800 C is approximately D = 7.10-8 cm2s-1. The jump distance
can be assumed to be 1.65 (1.65.10-8 cm), and then from Eq. (5.19) the jump
frequency is about 1.54.109 s-1. Thus, each oxygen atom makes a tremendously
large number of jumps per second. But it should then also be recalled that the
atoms vibrate with a (Debye) frequency of 1012-1013 s-1, and thus only a small
fraction - about 1 in 103 or 104 - of the vibrations leads to a jump of the atoms.
We shall consider this in more detail when the different atomistic mechanisms are
treated below.
Although the number of jumps is tremendously large, the mean
displacement of each atom is relatively small most of the time it moves back and
forth. In the diffusion process it is not possible to observe the individual jumps of
the atoms, and it is necessary to find a relation between the individual atom jumps
for large number of atoms and the diffusion phenomena which may be observed
on a macroscopic scale. The problem is to find how far a large number of atoms
will move from their original sites after having made a large number of jumps.
Such relations may be derived statistically by means of the so-called random walk
method.
Random diffusion
Let us consider that the jumps of the atoms are random, i.e. that the jumps of
the atoms are independent of all the previous jumps and can occur in all
directions. In that case the displacement of a diffusing atom from the starting point
after n number of jumps, Rn, is given by the algebraic sum of the individual jump
vectors:
n
Rn = s1 + s 2 + ...s n = s j
(5.22)
j =1
If the individual jumps take place with equal probability in all directions and
the individual jump distances are equal, this algebraic sum equals zero. This does
not mean that the diffusing atom remains at its starting point after n jumps, but
only that jumps in "positive " and "negative" directions are equally probable. In
fact the total displacement may have any value between zero and ns.
In order to obtain a value for the magnitude (length) of the sum vector, one
squares Eq. 5.22:
5.8
5 Diffusion
n
n 1
Rn2 = Rn = s j + 2
j =1
(5.23)
sk
j =1 k = j +1
If, as above, the individual jump vectors are equal, i.e., |s1| = |s2| ... |sj| = s,
as, for instance, in crystals with cubic symmetry, and if they are random and
uncorrelated, then the second term on the right hand side of the equality sign in
Eq. 5.23 will approach zero for large numbers of jumps, as on an average sj and sk
have an equal chance of being positive and negative. The first term on the righ
hand side, on the other hand, is always non-zero and positive, and thus represents
the average displacement length for a large number of jumps:
n
Rn2 = s j = ns 2
(5.24)
j =1
The mean displacement is given by the square root of Rn2 and is termed the root
mean square displacement:
Rn = Rn2 = n s
(5.25)
From this it is seen that the mean displacement is proportional to the square root
of the number of jumps times the individual jump distance.
By combining Eqs. 5.21 and 5.24 we may express the random diffusion in
terms of the diffusion coefficient we dealt with in diffusion down a concentration
gradient in the 3-dimensional cubic case:
Rn2 = ns 2 = 6 Dr t
(5.26a)
R n = 6 Dr t
(5.26b)
5 Diffusion
displacement in each orthogonal direction. Thus, the mean diffusion length in one
direction in a three-dimensional cubic crystal is:
(5.27)
x = 2 Dr t
This length is thus shorter than the displacement radius, since displacements in the
y and z directions are wasted for x dimension displacement.
As described above and using oxygen diffusion in niobium as an example,
the oxygen atoms on an average exchange positions approximately 1.54.109 times
per second at 800 C. From the same considerations one may also estimate that an
oxygen atom has randomly covered total jump distances of 25.4 cm and 914.7 m
after 1 second and 1 hour, respectively. But what is the mean displacement? From
Eq 5.27 one may estimate that the one-dimensional root-mean-square
displacement after 1 hour only amounts to 0.022 cm. Thus, the mean displacement
is very small and on an average the oxygen atoms spend most of their time
jumping "back and forth".
What we did above was to consider random jumps, and then we related
jumps, times and distances to a term we recognised from earlier, namely the
diffusion coefficient, and we named it the random diffusion coefficient Dr. The
diffusion coefficient was in turn something we recognised while considering
diffusion down a one-dimensional concentration gradient from Ficks first law.
One should note, however, that random diffusion and the random diffusion
coefficient can be considered and expressed and quantified also in the absence of
a concentration gradient and also for charged particles.
As will be described below, diffusion is often measured by using tracer
atoms and one then obtains values of the diffusion coefficient of the tracer atoms.
Depending on the diffusion mechanism the tracer diffusion is in most cases not
completely random, but is to some extent correlated with previous jumps. This
will be further discussed later on.
5.10
5 Diffusion
larger than in plane 2, since Fick's first law applies in both planes (and at any
position) at any time.
The change in concentration per unit time at any position is proportional to
the gradient in flux at that position
c
j
=t
x
(5.28)
Although the concentration and concentration gradient change with time, Fick's 1.
law is valid at any one time and position and thus
c
j
c
==
(D )
t
x x x
(5.29)
c
j
2c
==D 2
t
x
x
(5.30)
The intuition of the qualitative arguments as well as the double derivative of the
concentration in Eq.5.30 tell us that particles will flow from convex to concave
regions in terms of distribution of particles.
Figure 5-4. Schematic illustration of Ficks 2nd law; the concentration gradient changes
with time.
Equations 5.29 and 5.30 are representations of Fick's 2nd law. It may be
solved explicitly under certain boundary conditions that may be closely
approximated experimentally (e.g. Crank (1975)). A couple of examples of this
are given in the following.
5.11
5 Diffusion
c=
co
x2
exp(
)
4 Dt t
2(Dt t )1/2
(5.31)
c is the activity (or concentration) of the tracer at a distance x from the surface, co
is the activity originally present on the surface, and t is time of the diffusion
anneal. Dt is the tracer diffusion coefficient. Following Eq.5.31 it is determined by
1
plotting lnc vs x2, in which case the resultant straight line has the slope .
4 Dt t
Plots of c vs x and of lnc vs x2 according to Eq.5.31 is illustrated in Fig.5.5a and
b, respectively. At the point where the activity is half of the activity at the surface,
then x = 2.77 Dt t . This distance corresponds approximately to the root-meansquare penetration distance (Eq. 5.26b).
5.12
5 Diffusion
(5.32)
Dt t .
It may be noted that Eqs.5.31 and 5.32 involve the dimensionless parameter
x/(2 Dt t ), and accordingly the penetration and the amount of the diffusing
species dissolving in the solid are proportional to the square root of time.
Figure 5-6. Graphical presentation of diffusion profile when the surface concentration of
the diffusing species remains constant with time. It is assumed that co=0 (Eq.5.32).
The penetration by the diffusing species may be measured by means of the socalled sectioning method, that is cutting, grinding or etching off thin sections of
layers of the sample parallel to the plane surface and subsequently determining the
concentration of the diffusing species in each section. By cutting the specimen
normal to or at an angle to the plane surface, the penetration of radioactive tracers
may also be measured by means of so-called autoradiography.
5.13
5 Diffusion
Electron microprobe analysis and secondary ion mass spectroscopy (SIMS)
combined with sputtering techniques also provide excellent tools for studying
penetration and diffusion of foreign ions. As oxygen does not have a radioactive
isotope suitable for tracer studies, SIMS is particularly useful for studying oxygen
diffusion employing the stable 18O isotope.
In the isotope exchange method the vapour of the diffusing component
surrounding the sample is enriched with either a radioactive or an inactive isotope,
and the diffusion is measured by following the exchange of the isotope with the
sample. The diffusion coefficient may be evaluated if diffusion is the slower
process and processes at the surface are rapid. Alternatively, the sample itself may
be isotopically enriched, and the increase in concentration of the isotope in the
vapour phase may be measured.
Diffusion rates may, in principle, also be determined from any property or
reaction which depends on atomic mobility. By way of illustration, ionic
conductivity of the anion is directly proportional to the anion diffusion coefficient
(see Electrical conductivity). From high temperature solid state reactions,
sintering, oxidation of metals etc. diffusion coefficients may be evaluated
provided the detailed mechanism of the processes are known. Examples of this
will be given in Chapter 7.
Diffusion mechanisms
Lattice diffusion takes place through the movement of point defects. The
presence of different types of defects gives rise to different mechanisms of
diffusion. These are illustrated schematically for elemental solids in the following
descriptions. But they also apply to metal oxides and other inorganic compounds
when the diffusion is considered to take place in the sublattices of the cations or
anions.
Vacancy mechanism
The diffusion is said to take place by the vacancy mechanism if an atom on
a normal site jumps into an adjacent unoccupied lattice site (vacancy). This is
illustrated schematically in Fig.5.7. It should be noted that the atoms move in the
direction opposite the vacancies.
5.14
5 Diffusion
Interstitial mechanism
If an atom on an interstitial site moves to one of the neighbouring interstitial
sites, the diffusion occurs by an interstitial mechanism. This is schematically
shown in Fig.5.8. Such a movement or jump of the interstitial atom involves a
considerable distortion of the lattice, and this mechanism is probable when the
interstitial atom is smaller than the atoms on the normal lattice positions.
Diffusion of interstitially dissolved light atoms, e.g. H, C, N, and O in metals
provides the best known examples of this mechanism.
Interstitialcy mechanism
If the distortion becomes too large to make the interstitial mechanism
probable, interstitial atoms may move by another type of mechanism. In the
5.15
5 Diffusion
interstitialcy mechanism an interstitial atom pushes one of its nearest neighbours
on a normal lattice site into another interstitial position and itself occupies the
lattice site of the displaced atom. This mechanism is illustrated schematically in
Fig.5.9
In the interstitialcy mechanism one may distinguish between two types of
movements. If the atom on the normal lattice site is pushed in the same direction
as that of the interstitial atom, the jump is termed collinear (Fig.5.9). If the atom is
pushed to one of the other neighbouring sites so that the jump direction is different
from that of the interstitial atom, the jump is termed non-collinear.
Other mechanisms.
In elemental solids also other mechanisms have been proposed. The
crowdion is a variant of the interstitialcy mechanism. In this case it is assumed
that an extra atom is crowded into a line of atoms, and that it thereby displaces
several atoms along the line from their equilibrium positions. The energy to move
such a defect may be small, but it can only move along the line or along
equivalent directions.
For metals it has also been proposed that diffusion may take place through a
so-called ring mechanism, but this mechanism is improbable in oxides or other
inorganic compounds.
5 Diffusion
is displaced to the neighbouring vacancy. The movement of the cluster is
completed when atom 2 is displaced from a normal site to the interstitial site of
the new location of the cluster.
Figure 5-10. Proposed mechanism for the translational movement for a 4:1 cluster. Open
circles represent oxygen ions, closed circles iron atoms, open circles with an "x" iron interstitials,
and cubes iron vacancies.
5 Diffusion
ion has a smaller radius and charge than the oxide ion and may as such be
expected to have a smaller activation energy for diffusion than the oxygen ion.
Also other species such as water molecules and hydronium ions, H3O+, may serve
as vehicles for protonic diffusion, notably in relatively open structures.
Figure 5-11. Schematic illustration of free transport of protons in metal oxides (Grotthuss
mechanism).
Dr = 16 s 2 = 16 s 2
n
t
(5.19)
Vacancy diffusion.
Let us consider vacancy diffusion in an elemental solid or a cation or an
anion sub-lattice. The number of jumps per unit time, , depends on several
factors. First, it depends on the jump frequency towards an adjacent site.
Furthermore, it is also proportional to the number of sites to which the atom may
jump, i.e., the number of nearest neighbour positions of the atom, Z. Finally, the
atom may only jump if a vacancy is located on an adjacent site, and this
probability is given by the fraction (concentration) of vacancies in the crystal, Nd.
Thus, is in this case given by
5.18
5 Diffusion
(5.33)
= ZN d
(5.34)
(5.35)
where is a geometrical factor involving the factor 1/6 (from Eq.5.21, the factor
Z from Eq.5.33 and the relation between the jump distance and the lattice
parameter. For vacancy diffusion in a bcc lattice is thus equal to unity, From the
same considerations it may also be shown that = 1 for vacancy diffusion in fcc
lattices.
Interstitialcy diffusion
Consider an atom on a normal lattice site of a cation or anion sub-lattice. If
this atom is to move by the interstitialcy mechanism, an atom on a nearest
neighbour interstitial site has to push the atom on the normal site to a
neighbouring interstitial site. Thus for this diffusion mechanism an atom may only
diffuse when it has an interstitial atom on a neighbouring site, and as for vacancy
diffusion the diffusion coefficient of the atoms is proportional to the fraction
(concentration) of interstitial atoms or ions in the sub-lattice.
Interstitial diffusion
When one considers interstitial diffusion of an interstitially dissolved
species in dilute solid solution, essentially all the nearest neighbour interstitial
sites of the same type are unoccupied and available for occupancy by the diffusing
interstitial atoms. Thus the interstitial atom may jump to any of the nearest
neighbour interstitial sites and in this case Nd is equal to unity. The interstitial
diffusion coefficient is then given by
5.19
5 Diffusion
Dr = 16 s 2 Z
(5.36)
(5.37)
5 Diffusion
must be treated quantum mechanically and to their dynamics in a dynamic lattice
of much heavier oxygen ions. For instance, the lighter isotope has a lower sticking
probability after an otherwise successful jump. In effect, this reduces the effective
diffusivity of all hydrogen isotopes and it can to a varying degree counteract or
even seemingly reverse the classical effect. The neglectance of the latter have in
general made many investigators ascribe higher diffusivities for protons compared
to deuterons to the classical effect, while in reality the non-classical zero-point
energy difference appears to be the main contributor to the observed effect.
The possibility of tunnelling as a major component of diffusion is not
expected to apply to protons except at very low temperatures, and the orders-ofmagnitude isotope effects that would be expected for proton vs deuteron or triton
diffusion have not been reported for oxidic materials at elevated temperatures.
For the sake of simplicity let us first consider the case of diffusion by a
vacancy mechanism in a pure elemental solid. The diffusing atoms may only
make a jump when a neighbouring site is vacant. Thus the jump frequency n/t is
proportional to the vacancy concentration and as described in Chapter 3 the
fraction of vacancies may be written
N d = exp(
G d
S
H d
) = exp( d ) exp(
)
RT
R
RT
(5.38)
where Gd, Sd, and Hd denote the Gibbs free energy, entropy, and enthalpy of
formation of the vacancies (defects). In elemental solids Hd is positive and the
vacancy concentration increases with increasing temperature.
Vacancies in an oxygen deficient oxide.
5 Diffusion
of illustration, let us consider an oxygen deficient oxide, MaOb-, in which
doubly charged oxygen vacancies are the predominating point defects. If intrinsic
ionisation and effects of impurities can be neglected, the concentration of the
oxygen vacancies is as described in previous chapters (Eq. 3.68) given by
(5.39)
S v
O
3R
) exp(
H v
O
3RT
formation, respectively.
In such an oxide the oxygen vacancy concentration may also be determined
by the presence of dopants or a sufficiently large level of impurities with negative
effective charge. If the dopant is doubly negatively charged (as, for instance,
Ca2+, an acceptor dopant, in ZrO2), then
[vO ] = [A //M ]
(5.40)
In this case [vO ] will be independent of oxygen pressure and most often
temperature (see Chapter 4).
In a similar manner one may obtain the temperature and oxygen pressure
dependencies of the concentration of any defect (majority or minority) when the
defect structure is known. Of course the concentrations of acceptors or donors will
enter when they dominate the defect situation and the water vapour partial
pressure may enter when protons are dominating defects.
5.22
5 Diffusion
Figure 5-12. Potential energy of atom diffusing in a solid. Hm is the activation energy for
the diffusion.
= exp
G m
S
H m
= exp m exp
RT
R
RT
(5.41)
where Gm, Sm, and Hm represent the free energy, entropy and enthalpy
change, respectively, connected with the movement of the atom from the
equilibrium position to the top of the potential barrier and (nu) represents the
vibration frequency. is often assumed to equal the Debye frequency, i.e. about
1013 Hz, as an order of magnitude approximation. (For protons in oxides the O-H
stretching frequency of about 1014 Hz is used.) Sm is often assumed to be a
small, positive term (below, say, 10 J/molK).
As a more detailed analysis of , Zener suggested that may be
approximated as = /a Hm / M , where is a structure- and mechanismdependent factor, a is a lattice parameter, and M is the reduced mass of the
oscillator. Intuitively, a particle vibrating in an energy valley will vibrate faster (
increases) when the walls become steeper (Hm increases). One may also view
the particle as vibrating on a spring; the frequency becomes higher when she
spring is shorter (/a decreases), when the spring is stiffer (Hm increases) or the
particle becomes lighter (M decreases). Thus, the temperature-independent (preexponential) term in the diffusion coefficient increases when Hm increases, so
that the two tend to counteract each other. Experimental observations of this is
sometimes referred to as the Meyer-Neldel effect (Meyer and Neldel (????)).
5.23
5 Diffusion
Dr = a 02 exp
S d + S m
(H d + H m )
exp
R
RT
(5.42)
D = D0 exp
Q
RT
(5.43)
where Q is termed the activation energy and Do the the pre-exponential factor.
By comparing Eqs.5.42 and 5.43 it is seen that the activation energy, Q, in
this case comprises
(5.44)
Q = Hd + Hm
Correspondingly, Do is given by
D0 = a 02 exp
S d + S m
R
(5.45)
If experiments are carried out under such conditions that the concentration
of defects, Nd, is constant and independent of temperature, e.g. at sufficiently low
temperatures that the defect concentration is frozen in, then Dr is given by
S m
H m
H m
exp
= D0 exp
R
RT
RT
(5.46)
and correspondingly the activation energy is under these conditions given simply
by
Q = Hm
(5.47)
5.24
5 Diffusion
Oxygen vacancy diffusion in oxygen-deficient oxides.
2
0
1 1/ 3
4
Dr = a ( )
1 / 6
O2
exp
S m
exp
H v
O
RT
H m )
(5.48)
Dr thus increases with decreasing oxygen pressure. The activation energy for the
diffusion is in this case given by
Q=
H v
O
(5.49)
H m
S m
H m
exp
R
RT
(5.50)
Thus in this case the activation energy is equal to that of the mobility of the
oxygen vacancies: Q = Hm. A similar situation would arise if the vacancy
concentration was frozen in rather than determined by acceptor doping.
Such a transition from intrinsic to extrinsic diffusion may take place when
the temperature is lowered from high temperatures, where the native point defects
predominate, to low temperatures, where the point defect concentrations are
determined by the impurity concentration or, in other cases, frozen in. The
temperature dependencies and the corresponding change in activation energy of
the random diffusion in such a case are illustrated in Fig. 5.13.
The above situations represent ideal cases. As the temperature is decreased
defect interactions may become increasingly important. This may be treated as
formation of associated defects. In the non-stoichiometric (intrinsic) case one may
for instance have association between charged vacancies and electrons forming
singly charged or neutral vacancies. This will change the oxygen pressure
dependency of Nd and probably the temperature dependencies of Nd and , but
probably not dramatically. Of larger effect, and more frequently observed, are the
associations between the mobile vacancies and the relatively stationary acceptors
in the extrinsic regime: The associated vacancies can be regarded as immobilised,
5.25
5 Diffusion
and the concentration of mobile vacancies (Nd) starts to decrease with decreasing
temperature. This is seen as an increasing activation energy of diffusion with
decreasing temperature in many heavily doped oxides. One may note that instead
of expressing the effect as a changing concentration of free vacancies one may
express it as a changing mobility in that the activation energy for diffusion is
increased by the trapping energy exerted by the acceptor. However, the simple
model of free and stationary (associated) point defects appears capable of
explaining most behaviours fairly well.
Figure 5-13. The diffusion coefficient for oxygen diffusion by the vacancy mechanism in an
oxygen deficient oxide in which oxygen vacancies are the predominant native point defects. At high
temperatures the oxide exhibits intrinsic behaviour and at reduced temperatures extrinsic
behaviour (i.e. the oxygen vacancy concentration is determined by the concentration of lower
valent cations).
D0 = 16 a 02 exp
S m
R
(5.51)
5.26
5 Diffusion
Assuming that ~ 1013 s-1 and ao = 1.5.10-8 cm, and as it is probable that Sm >
0, one may estimate a lower limit of Do of
Do > 5.10-4 cm2s-1
(5.52)
(5.53)
5.27
5 Diffusion
Dr N = DV Nd
(5.54)
where, as stated above, N denotes the fraction of sites occupied by atoms, often
approximated as ~1. This relation can be generalised to be a very important and
useful approximation for any point defect:
Dr N = Dd Nd
(5.55)
Excercises
1. Random (self) diffusion
a) The self-diffusion coefficient of a metal with cubic structure can be
expressed as
D=
1n 2
s
6t
where n/t represents the jump frequency (i.e. number of jumps n over a time t).
Close to the melting point most fcc and bcc metals have D 10-8 cm2/s.
i) If the jump distance is 3 , what is the jump frequency near the melting
point?
ii) What is the relation between this frequency and the vibrational
frequency?
iii) How far has one atom traveled after 1 hour?
iv) What is the root mean square displacement after one hour?
v) What is the root mean square displacement in one dimension after one
hour?
5.28
5 Diffusion
b) For a metal with cubic structure the diffusion coefficient can also be
expressed as
D = a 02N d
where is a geometric factor, a0 the lattice constant, is the jump frequency, and
Nd is the defect concentration. Derive the value for for vacancy diffusion in a
metal with fcc structure.
5.29
5 Diffusion
5.30
6. Electrical conductivity
6. Electrical conductivity
Introduction
In the preceding chapter we have described and discussed diffusion of
particles in solids and particularly of ions and defects in metal oxides. The driving
force for the diffusion has been taken to be the negative value of the particle
gradient or more precisely the negative value of the chemical potential gradient.
When using isotopes as tracers one may study self-diffusion, i.e. diffusion of the
components in the oxide (metal and oxygen ions) in a homogeneous oxide; in this
case the isotopic tracer gradient is the driving force for the diffusion.
In this chapter the transport of electrical charges will be described and
discussed. In metal oxides the electrically charged particles comprise ions and
electrons. The ionic charge carriers comprise the cations, anions, and foreign ions
(e.g. impurity ions, dopant ions and protons) and the electronic charge carriers are
the electrons and electron holes. The concentrations of the charge carriers are
directly related to the defect structure of the oxide and in this chapter we will
derive expressions for the temperature and oxygen pressure dependence of the
electrical conductivity. The discussion will be limited to transport of charges in
chemically homogeneous metal oxides (no chemical potential gradient) but with
an electrical potential gradient as the driving force. In the next chapter transport of
ionic and electronic charge carriers in metal oxides which are simultaneously
exposed to chemical and electrical potential gradients, i.e. electrochemical
potential gradients, will be discussed.
As the mobilities of electrons and electrons holes are normally much higher
than those of ions, most oxides are electronic conductors. One type of charge
carrier often predominates in an oxide under particular conditions of temperature
and oxygen pressure. An electronically conducting oxide is an n-conductor if
transport of electrons predominate and a p-conductor if electron holes prevail.
However, some oxides are or may become ionic conductors or mixed
ionic/electronic conductors depending on the temperature and oxygen pressure
often as a result of appropriate doping with aliovalent foreign ions. Some oxides
may also exhibit proton conductivity in hydrogen- or water vapour-containing
atmospheres; predominant proton conductivity in such oxides is in some cases
observed at reduced temperatures (< 600-700 C).
6.1
6. Electrical conductivity
d
F = -zie dx = zieE
(6.1)
where is the electrical potential and E = - dx is the electric field. The flux of
particles of type i is the product of the concentration ci, the particle mobility Bi,
and the force F:
ji = ciBi F = zie ciBiE
(6.2)
(6.3)
(6.4)
(6.5)
where i = zie ciui is the electrical conductivity due to the charge carriers of type i.
The electrical conductivity is determined by the product of the concentration ci of
the charged particles, the charge zie on the particles and the charge carrier
mobility, ui. It should be noted that Eq.6.4 is an expression of Ohm's law. The unit
for the electrical conductivity is Siemens per cm, Scm-1 (one Siemens is the
reciprocal of one ohm and in older literature the electrical conductivity is
expressed as ohm-1cm-1). The unit for the charge is coulomb, the concentration of
charge carriers is expressed as the number of charge carriers of type i per cm3, and
charge carrier mobility in units of cm2/Vs. (Although the SI unit for length is m,
cm is being used in the following as it is still by far the one most commonly used
in the literature).
It may be noted that in the above terminology, F, E, ii, zi, ui and ji may each
be positive or negative. ui and zi always have the same sign, and as long as no
other forces than the the electrical act, ii and Ei always have the same sign, and ji
and F always have the same sign. Bi and i are always positive, and it is common
also to neglect the sign when specifying charge mobilities ui.
6.2
6. Electrical conductivity
(6.6)
The ratio of the partial conductivity i to the total conductivity is termed the
transport (or transference) number of species i:
ti =
(6.7)
(6.8)
where c, a, n, and p are the cation, anion, electron and electron hole
conductivities, respectively.
Following Eq. 6.7 the individual conductivities may be written in terms of
their transport numbers: c = tc , a = ta, n = tn and p = tp. Using these
values Eq.6.7 takes the form
= (tc + ta + tn + tp)
(6.9)
It may be noted that the sum of the transport numbers of all the charge carriers
equals unity:
tc + ta + tn + tp = 1
(6.10)
The total electrical conductivity is often given by the sum of the ionic
conductivity, ion = c + a, and the electronic conductivity, el = n + p, and
the total conductivity can then be written
6.3
6. Electrical conductivity
= ion + el
(6.11)
Often only one type of charge carrier dominates the charge transport, and in
many cases and as an approximation contributions from minority carriers are
neglected. For oxides the mobilities of electrons and electron holes are usually
several orders of magnitude (~104 - 108) larger than those of the ions, and even
when the concentration of electron or electrons holes is smaller than that of the
ionic charge carriers (or, more precisely, than that of ionic charge carrier defects)
the oxide may still be a predominantly electronic conductor. The relative
importance of ionic and electronic conductivity will often vary greatly with
temperature and oxygen pressure. This will be illustrated in the following
chapters.
(6.12)
By combining this relation with Eqs. 6.4 and 6.5 one obtains the following
relation between the random diffusion coefficient and the charge carrier mobility
and the electrical conductivity:
kT
kT
Di = kTBi = ui z e = i 2
i
cizi e2
(6.13)
This relation is called the Nernst-Einstein relation. This relation and also the
effect of an applied electric field on migration of charged species in a
homogeneous crystal may also be derived from the following model, in which we
will understand also when and why conduction is termed a linear process.
Consider a one-dimensional system with a series of parallel planes separated
by a distance s (cf. Fick's first law in Chapter 5). It is assumed that the system is
homogeneous and that the volume concentration of the particles in the planes is ci.
The particles in neighbouring planes 1 and 2 have equal probability of jumping to
the neighbouring planes. In the absence of any external kinetic force, the number
of particles which jump from plane 1 to plane 2 and from 2 to 1 per unit time is
6.4
6. Electrical conductivity
1
equal and opposite and given by 2 cis. In a homogeneous system there will be no
net transport of particles.
When there is no applied electric field, the activation energy associated with
the jumps is Hm. When an electric field E is applied, the jump frequency in the
positive direction will be increased and that in the negative direction decreased in
that the activation energies are changed. In the forward direction the activation
1
energy is reduced to Hm- 2 ziesE and in the reverse direction increased to Hm+
1
2 ziesE.
The net particle flux is given by the difference in number of jumps in the
forward and reverse directions:
1
(6.14)
z esE
H m - i
S m
2
where forw = exp(
) exp(kT
k
and
z esE
H m + i
Sm
2 ).
rev = exp(
) exp(k
kT
Reverse jump
Forward jump
ciesE/2
Y
G
R
E
N
E
Rest position
-ci esE/2
Hm +ciesE/2
H m
Hm -ciesE/2
s/2
DISTANCE
6.5
6. Electrical conductivity
Equation 6.14 then becomes
ziesE
ziesE
1
ji = 2 cis{exp ( 2kT ) - exp (- 2kT )}
where = exp(
(6.15)
Hm
Sm
) exp(- kT ).
k
When ziesE << 2kT, which is valid for normal electrical measurements in
bulk materials (and when Ohm's law is applicable), the difference in exponentials
in Eq. 6.15 may be written ziesE/kT (since ex - e-x = 2x for x<<1). Eq. 6.15 then
takes the form
zieE
1
ji = 2 s2ci kT
(6.16)
(6.17)
We have by this shown that net flux density of a hopping charge carrier in an
electrical field in the small-signal (linear) range is proportional to the random
diffusion coefficient. As we have shown before, the flux of particles with a charge
zie may also be expressed in terms of charge carrier mobility or conductivity:
iE
ji = zie ciBiE = ci ui E = z e
i
(6.18)
and when one combines Eqs. 6.17 and 6.18 one obtains various forms of the
Nernst-Einstein relation (Eq. 6.13):
kT
kT
Di = Bi kT = ui z e = i 2
i
cizi e2
(6.13)
or, rearranged,
i = (zie)2ciDi/kT
6.6
6. Electrical conductivity
It is emphasised that the relation is derived assuming random diffusion and that
the mobilities and conductivity through this relation connects to the random (or
self) diffusion coefficient Dr. It is thus meaningful for relating electrical and
diffusional transport of atoms and ions. For electrons and holes this is only
meaningful when they migrate by an activated hopping mechanism.
From the Nernst-Einstein relation it is also seen that the temperature
dependence of the product iT is the same as that of Dr. Thus in evaluating the
activation energy associated with the diffusion coefficient from conductivity
measurements, it is necessary to plot (iT) vs 1/T.
It is also important to note that in the derivation it is implicitly assumed that
the ions and electrons move independently of each other, e.g. that there is no
interference between ionic and electronic flows.
(6.19)
where n and p are the electron and electron hole conductivities, n and p
the charge carrier concentrations of electrons and electron holes, respectively, and
un and up are the carrier (or drift) mobilities of electrons and electron holes. As
mentioned above, one type of charge carrier will often dominate; however, in
special cases where an oxide is close to stoichiometric both n- and p- conductivity
may contribute significantly to the electronic conductivity.
In defect-chemical equations the concentrations of electrons and electron
holes are often written in terms of the law of mass action without specifying
where the electronic defects are located. The concentrations of electronic defects
are often interpreted in terms of the band theory of solids and in the following is
6.7
6. Electrical conductivity
given a brief account of this theory and the relationship between this theory and
the law of mass action.
(6.20)
6.8
6. Electrical conductivity
In materials where the valence band is only partly occupied or where it
overlaps with the conduction band (no forbidden band gap) the electrons can
move freely in the available energy levels and we have metallic conduction.
In an insulator or semiconductor at 0 K there is a band gap and the valence
band is completely filled with electrons while the conduction band is completely
empty. In such cases no electronic conduction takes place when an electric field is
applied.
Conduction band
Electron energy
EC
Forbidden
energy gap
Eg=EC-EV
EV
Valence band
Distance through crystal
Figure 6-2. Schematic illustration of the energy band diagram for a pure
semiconductor
Intrinsic ionisation
When the temperature is increased, electrons in the valence band are excited
across the forbidden energy gap to the conduction band. This is the intrinsic
ionisation. The electrons in the conduction band and the unoccupied electron sites
in the valence band (electron holes) can move in an electric field. The electron
holes behave as though they were positively charged and move in the opposite
direction of the electrons. The intrinsic ionisation thus produces pairs of electron +
electron hole charge carriers. When the electronic conductivity is due to intrinsic
ionisation only, the semiconductor is called an intrinsic semiconductor. According
to this model electron and electron hole conductivities increase with increasing
concentrations of electrons in the conduction band and electron holes in the
valence band.
In many oxides, and particularly those with large percentage ionic bonding,
periodic fluctuations of the electric potential associated with each ion become too
large (and energy bands too narrow) so that the band model provides an
inadequate description or theory. In this case the electrons or holes may be
considered to be localised at the lattice atoms. Such localised electronic defects
are termed polarons and may from a chemical point of view be considered to
constitute valence defects.
6.9
6. Electrical conductivity
When classical statistics (the Boltzmann approximation) can be used it may
be shown that the concentration of electrons and electron holes (expressed in
number per cm3) is given by
EC-EF
n = NC exp(- kT )
(6.21)
EF-EV
p = NV exp (- kT )
(6.22)
NC and NV represent the number of available states per unit volume - density of
states - in the conduction and valence bands. EC is the energy of the lowest level
of the conduction band and EV the highest level in the valence band. The
parameter EF is termed the Fermi level and we will return to this later on.
When the electrons occupy a narrow band of energies close to EC, a
parabolic relation between NC and EC may be assumed and it may be shown that
NC is then given by
8m*e kT 3/2
NC = (
)
h2
(6.23)
where me is the effective mass of the electron. When the holes occupy a narrow
region of energies close to EV, NV is correspondingly given by
NV = (
8m*h kT 3/2
)
h2
(6.24)
where mh is the effective mass of the electron hole. The effective masses of
electrons and electron holes are seldom accurately known and as an
approximation the free electron mass is then used in these relations.
The defect reaction and corresponding equilibrium between electrons in the
conduction band and electron holes in the valence band can in terms of the law of
mass action be expressed as
0 = e' + h.
(6.25)
n.p = Ki
(6.26)
6.10
6. Electrical conductivity
where Ki is the equilibrium constant for the intrinsic ionisation. Combination with
Eqs. 6.21 and 6.22 gives
Eg
Ki = n.p = NC NV exp (- kT )
(6.27)
where Eg is the band gap, Eg = EC - EV. Eg may thus be considered the energy of
the intrinsic ionisation. It is emphasised that Eq. 6.27 presupposes that classical
statistics apply. It may also be noted that the Fermi level is eliminated in this
expression of the law of mass action. It may be noted again that in Eq. 6.27 the
concentrations of electrons and electron holes are expressed in terms of number
per volume unit (e.g. per cm3).
Intrinsic electronic semiconductor
(6.28)
(6.29)
It may be noted that NC and NV are temperature dependent and that also un and up
may have various dependencies on temperature. If the latter are not exponential
(as in diffusional hopping conduction processes) the exponential term of the
energy gap tend to dominate the temperature dependence and as an approximation
Eq. 6.29 is then often written
Eg
(6.30)
6.11
6. Electrical conductivity
Effects of charged
semiconductors
impurities
or
defects.
Extrinsic
6.12
6. Electrical conductivity
defects at 0 K, these point defects must themselves be charged at 0 K and thus
have levels outside the band gap.
In the following we consider imperfections introducing levels inside the forbidden
energy gap as illustrated in Fig.6.3. When the imperfection introduces an energy
level which is located below the lower edge of the conduction band (EC),
electrons can be thermally excited to the conduction band. As such the
imperfection donates an electron, it is called a donor and the corresponding energy
level a donor level. Under conditions where such imperfections dominate the
defect structure the oxide becomes an n-conductor. Correspondingly,
imperfections with energy levels just above the upper edge of the valence band
(EV) are termed acceptors; electrons in the valence band may be excited to the
energy level of the imperfection (acceptor level) and the oxide may become a pconductor.
Electron energy
Conduction band
ED
Donor level
EC
Ed
Acceptor level
EA
Ea
EV
Valence band
Distance through crystal
6.13
6. Electrical conductivity
Effects of donors
Let us again describe these processes in terms of the law of mass action.
Thus the ionisation of a donor Dx may be written
Dx = D. + e'
(6.31)
(6.32)
(6.33)
If the energy of the donor state is ED, and the donor state is mainly empty
([Dx]<<ND) then the concentration of electrons in the donor state, [Dx], may
following Boltzmann statistics (be expressed by
[Dx] = ND exp (-
(ED-EF)
kT
(6.34)
[Dx]
Since we assume N <<1, and that Boltzmann statistics also apply to conduction
D
band electrons and by further combining Eqs. 6.21 and 6.32-6.34, we get an
expression for KD:
KD = NC exp (-
EC-ED
kT
Ed
) = NC exp (- kT )
(6.35)
where Ed = EC-ED represents the ionisation energy of the donor (cf. Fig. 6.3). Ed
may in Eq. 6.35 be considered to be the enthalpy of the ionisation of the donor.
When no other imperfections are present, the situation described here can be
approximated by
n = [D.] ND
(6.36)
6.14
6. Electrical conductivity
and if the donors are present in an invariable amount as a fully soluble or frozenin impurity, then n = [D] ND = constant. From Eq. 6.32, we then get the
2 -1
minority concentration of neutral unionised donors: [Dx] = ND NC
exp(+Ed/kT).
If, on the other hand the donor level is mainly unionised, then [D]<<[Dx]
ND and Boltzmann statistics give
(EF-ED)
[D] = ND exp(- kT )
(6.37)
and then, by combination with Eqs. 6.21, 6.32, and 6.33 we again obtain KD = NC
exp(-Ed/kT) as in Eq. 6.35.
If the electroneutrality condition in the oxide is still given by
n = [D. ]
(6.38)
(but now [D]<<ND), then the concentration of electrons is, from insertion in Eq.
6.32:
E
d
n = [D] = (KD ND)1/2 = (NC ND)1/2 exp (- 2kT )
(6.39)
The reader may find it useful to note the similarity between the form of this
expression and Eq. 6.28. However, we also remind ourselves that ND may be a
constant concentration of impurities or it may be a constant or varying level of
other point defects such as anion vacancies or cation interstitials.
In the above we have seen that the same equilibrium constant KD for Eq.
6.31 applies whether the donor level is approximately empty (fully ionised) or
approximately fully occupied (practically unionised, or, in this case, neutral). But
our treatment does not hold for intermediate situations, i.e. what we may call
partly filled donor levels. In the latter case Boltzmann statistics do not represent a
sufficiently good approximation, and instead Fermi-Dirac statistics must be used.
Let us recapitulate briefly the treatment of the ionisation of a constant
concentration of a donor in the forbidden gap of an otherwise pure, stoichiometric,
ideal semiconductor: At low temperatures, the concentration of electrons is given
by a minor degree of ionisation of the donors, as given by Eq. 6.39. If we neglect
the temperature dependence of NC, the situation can be illustrated as in the righthand part of Fig. 6.4; the concentration of electrons increases with an apparent
6.15
6. Electrical conductivity
enthalpy of Ed/2. At a sufficiently high temperature (middle part of Fig. 6.4) all
donors are ionised, and the concentration of electrons becomes constant (Eq.
6.36). At even higher temperatures intrinsic semiconduction may predominate,
i.e., n = p, and in principle the temperature dependence becomes as illustrated in
the left hand part of Fig. 6.4, with an apparent enthalpy close to Eg/2 (cf. Eq.
6.25). Thus at low temperatures this oxide is an n-conductor due to the ionisation
of the donors and at high temperatures intrinsic ionisation predominates. The
behaviour over the entire temperature range could in principle be solved from the
full electroneutrality equation
n = [D]+ p
(6.40)
in combination with the constancy of the donor concentration and the expressions
for the equilibrium constants Ki (Eq. 6.27) and KD (Eq. 6.35). However, as noted
above, the expression for KD would not be valid between the intermediate and
low temperature (right hand side) domains in Fig. 6.4.
Log n
Eg/2
n=N D
Ed /2
1/T K
(6.41)
6.16
6. Electrical conductivity
The equilibrium constant for the defect reaction is given by
[A'] p
= KA
[Ax]
(6.42)
Following a similar treatment as for the donor, KA may under the assumption that
Boltzmann statistics apply to the valence band (small concentration of holes) and
the acceptor (nearly full or nearly empty acceptor levels) be expressed by
EA-EV
Ea
KA = NV exp (- kT ) = NV exp (-kT )
(6.43)
where EA is the energy level of the acceptor and Ea = EA-EV is called the
ionisation energy of the acceptor (cf. Fig. 6.3).
For an acceptor doped oxide the temperature dependence of the
concentration of holes will be qualitatively analogous to that of electrons in the
donor doped case (cf. Fig. 6.4) with the same constrictions.
(6.44)
The chemical potential of electrons, for instance, the conduction electrons, may be
written in terms of the chemical potential in a standard state, (cond. electrons)
and a term for the entropy of mixing:
6.17
6. Electrical conductivity
n
e = (cond. electrons) = (cond. electrons) + kT ln N
(6.45)
(6.46)
el = {2e(
Eg
2k 3/2 * . * 3/4 3/2
) (me mh ) T exp (- 2kT )}(un + up)
2
h
(6.47)
completely free electron is the mass equal to the real mass, me=me . As mentioned
above, the values of the effective masses are not accurately known, and the value
of me is often used as an approximation.
In order to obtain an accurate description of the temperature dependence of
the electronic conductivity it is necessary to consider the temperature
dependencies of the charge carrier mobilities.
6.18
6. Electrical conductivity
Non-polar solids itinerant electron model
(6.48)
Eg
(6.49)
If, on the other hand, the scattering is mainly due to irregularities caused by
impurities or other imperfections, the charge carrier mobility is proportional to
T3/2, e.g.
un,imp = const. T3/2
(6.50)
(6.51)
and from the temperature dependencies given above it is evident that impurity
scattering dominates at low temperature while lattice scattering takes over at
higher temperature.
Polar (ionic) oxides
When electrons and electron holes move through polar compounds, such as
ionic oxides, they polarise the neighbouring lattice and thereby cause a local
deformation of the structure. Such an electron or electron hole with the local
deformation is termed a polaron. The polaron is considered as a fictitious single
particle.
6.19
6. Electrical conductivity
When the interaction between the electron or electron hole and the lattice is
relatively weak, the polaron is referred to as a large polaron. Large polarons
behave much like free carriers except for an increased mass caused by the fact that
polarons carry their associate deformations. Large polarons still move in bands,
and the expressions for the effective density of states in the valence and
conduction bands are valid. The temperature dependence of the mobilities of large
polarons at high temperatures* is given by
ularge pol. = const. T-1/2
(6.52)
(6.53)
Thus in this case the value of the band gap can be deduced by plotting log(el/T)
vs 1/T.
The large polaron mechanism has been suggested for highly ionic nontransition metal oxides. Such oxides have large energy gaps (Eg 6 eV), where
the band structure is characterised by large band widths. Due to the large band gap
one expects electronic semiconductivity to be low in these oxides, and to
predominate possibly only at high temperatures.
For other oxides it has been suggested that the interactions between the
electronic defects and the surrounding lattice can be relatively strong and more
localised. If the dimension of the polaron is smaller than the lattice parameter, it is
called a small polaron or localised polaron, and the corresponding electronic
conduction mechanism is called a small polaron mechanism.
The transport of small polarons in an ionic solid may take place by two
different mechanisms. At low temperatures small polarons may tunnel between
localised sites in what is referred to as a narrow band. The temperature
dependence of the mobility is determined by lattice scattering and the polaron
mobility decreases with increasing temperature in a manner analogous to a broad
band semiconductor.
However, at high temperatures (for oxides above roughly 500 C) the band
theory provides an inadequate description of the electronic conduction
mechanism. The energy levels of electrons and electron holes do not form bands,
but are localised on specific atoms of the crystal structure (valence defects). It is
*
"High temperatures" are temperatures above the optical Debye temperature, . For oxides
~(h)/2k, where h is the Planck constant, k the Boltzmann constant and the longitudinal
optical frequency which for an oxide is ~1014 s-1.
6.20
6. Electrical conductivity
assumed that an electron or electron hole is self-trapped at a given lattice site, and
that the electron (or electron hole) can only move to an adjacent site by an
activated hopping process similar to that of ionic conduction. Consequently it has
been suggested that the mobility of a small polaron can be described by a classical
diffusion theory as described in a preceding chapter and that the Nernst -Einstein
can be used to relate the activation energy of hopping, Eu, with the temperature
dependence of the mobility, u, of an electron or electron hole:
e
Eu
(6.54)
Nonstoichiometric semiconductors
Corresponding expressions for el for nonstoichiometric electronic
semiconductors readily follows by considering the temperature and oxygen
pressure dependence of the concentration of the electronic defects.
For nonstoichiometric oxides the concentration of electronic defects is
determined by the deviation from stoichiometry, the presence of native charged
point defects, aliovalent impurities and/or dopants. The concentration of
electronic defects can be evaluated from proper defect structure models and
equilibria. Various defect structure situations have been described in previous
chapters and at this stage only one example - dealing with oxygen deficient oxides
with doubly charged oxygen vacancies as the prevalent point defects - will be
described to illustrate the electrical conductivity in nonstoichiometric oxides.
Oxygen deficient oxides
6.21
6. Electrical conductivity
2.
OO = VO + 2e' + 2 O2
(6.55)
(6.56)
n = 2[VO ]
(6.57)
(6.58)
The total electrical conductivity is given by the sum of the conductivity of the
electrons and of the oxygen vacancies:
2.
t = 2 e [VO ] uV2. + e n un
(6.59)
2.
2 e [VO ] uV2. represents the ionic conductivity due to the oxygen vacancies and
O
.
2
where uV is the mobility of the oxygen vacancies. However, if the electrons and
O
oxygen vacancies are the prevalent charge carriers, the contribution due to oxygen
vacancies can be ignored due to the much higher mobility of electrons than
oxygen vacancies, and the oxide is an n-conductor where the conductivity can
then be written
-1/6
t = n = e n un = e un (2KV2. )1/3 pO
O
6.22
(6.60)
6. Electrical conductivity
As described in previous chapters the equilibrium constant for the formation of
doubly charged oxygen vacancies and 2 electrons is given by
2.
KV 2 .
O
= exp (
SV O
k
) exp(-
H V 2.O
kT
(6.61)
When one combines Eqs. 6.61 and 6.62 the n-conductivity may be written:
t = n = e un exp (
SV 2.O
3k
) exp (-
H V 2.O
-1/6
) pO2
3kT
(6.62)
Let us further assume that the electrons are small polarons and thus that the
mobility of the electrons are given by Eq. 6.54. The conductivity can then be
expressed by
n = const.
1
T
exp (-
H V 2.O / 3 + E u
kT
-1/6
) pO2
(6.63)
-1/6
Thus following this equation the n-conductivity is proportional to pO , and if this
2
Log n-conductivity
defect structure situation prevails over a temperature range from T1 to T5, one will
obtain a set of isotherms of the n-conductivity as shown in Fig.6.5.
pO2 -1/6
T5
T4
T
T2 3
T1
6.23
6. Electrical conductivity
Furthermore, if it can be assumed that mobility of the charge carriers (defect
electrons) is independent of the defect concentration, then a plot of the values of
log10(T) at a constant oxygen pressure yields a straight-line relationship as
H V 2.O
1
illustrated in Fig. 6.6. The slope of the line is given by +Eu,
2.303k
3
where the factor 2.303 is the conversion factor in changing from lne to log10. The
activation energy is given by the term
E =
H V 2.O
3
+ Eu.
(6.64)
6. Electrical conductivity
described in the next chapter (Ch. 7) dealing with electrochemical transport in
metal oxides.
Ionic conductivity
Ionic conductivity follows the NernstEinstein relationship, for hopping
species, derived early in this chapter:
i = zie ci ui = (zie)2ciDi/kT
(6.65)
It may be noted that the species considered may be a defect (e.g. oxygen
vacancies) or a constituent (oxide ions). For defects the mobility and diffusivity
are large and roughly constant, while the concentration is small and variable. For
the constituent, the mobility and diffusivity are small and variable (with defect
concentration) while the concentration is large and roughly constant. The
conductivity (here oxide ion conductivity by the vacancy mechanism) is the same.
In many cases of utilizing ionic conduction, the concentration ci of ionic
defects is constant, given by a dopant. For instance, ionic conductivity in yttriadoped zirconia is determined by the concentration of oxide ion vacancies, in turn
given as charge compensating the yttrium acceptors, 2[vO..] = [YZr/]. If the
concentration of acceptors is given in mole-fraction, then it is necessary to
multiply the resulting mole fraction of vacancies by the formula density or molar
density of the compound in order to obtain the volume density required for
insertion in Eq. 6.65:
vO.. = 2e [YZr/]/2 * cZrO2 uvO.. = 2F [YZr/]/2 * CZrO2 uvO..
(6.66)
where cZrO2 and CZrO2 are, respectively, the molecular and molar densities (number
of formula units or moles per unit volume) of the oxide.
There are of course also cases also of practical interest where the
concentration of ionic defects vary, e.g. with temperature in intrinsically
disordered compounds, and with temperature and non-stoichiometry in nonstoichiometric compounds. In proton conducting oxides the proton conductivity
varies with proton concentration, typically a function of water vapour partial
pressure.
We leave further learning about ionic conduction to exercises and Chapter
7 (electrochemical transport), and here only briefly mention a couple of aspects of
ionic transport that relates it in more detail to diffusion.
6.25
6. Electrical conductivity
(6.67)
The value of f is governed by the crystal structure and the diffusion mechanism.
Dt ci (zie)2
kT
f =D =
r
i
(6.68)
This equation is applicable to any diffusion process for which the atom jump
distance is equal to the displacement of the effective charge, e.g. for vacancy and
interstitial diffusion.
However, in interstitialcy diffusion the charge displacement is larger than
the atom jump distance, and a displacement factor S must be included in the
Nernst-Einstein relation. In collinear interstitialcy diffusion (Fig. 5.9) the effective
charge is, for instance, moved a distance twice that of the tracer atom and Dt/Dr is
given by
Dt
Dr
Dt ci (zie)2
ikT
= S
=S
(collinear)
(6.69)
6. Electrical conductivity
ratios of Dt/Dr are in good agreement with the assumption that f = 0.78, i.e. that
the Na-ions diffuse by a vacancy mechanism.
However, such a simple relationship was not found for AgBr. AgBr is also a
cationic conductor and comparative values of Dt (diffusion of Ag in AgBr) and of
values of Dr evaluated from conductivity measurements are shown in Fig. 6.7.
From studies of the effect of Cd-dopants on the ionic conductivity it could
be concluded that cationic Frenkel defects predominate in AgBr. Thus the
diffusion was therefore expected to involve both vacancy diffusion and transport
of interstitial ions. The experimentally measured ratios of Dt/Dr varied from 0.46
at 150 C to 0.67 at 350C. For vacancy diffusion a constant ratio of 0.78 (=f)
would have been expected, and the diffusion mechanism could thus be ruled out.
For interstitial diffusion f=1, and this mechanism could also be excluded.
6.27
6. Electrical conductivity
(x - x) 2
exp ()
4D t t
(6.70)
N
O
I
T
A
R
T
N
E
C
N
O
C
X = uiEt
Literature
Friauf, R.T. (1957) Phys. Rev. 105, 843; (1962) J. Appl. Phys. 33 suppl.,
494.
Manning, J.R. (1962) J. Appl. Phys 33, 2145; Phys. Rev. 125, 103.
6.28
6. Electrical conductivity
Exercises
1. Cobalt oxide: The electronic conductivity of Co1-yO at 1350C and pO2 =
0.1 atm is 25 S/cm. Thermogravimetric measurements show that y = 0.008 under
the same conditions. It is assumed that singly charged cobalt vacancies are the
dominating point defects. Identify the charge carriers responsible for the
conductivity and calculate their charge mobility. (Assume that the density of CoO
at 1350C equals that at room temperature, 6.4 g/cm3. Atomic weights MCo =
58.93, MO = 16.00.)
2. Nickel oxide: Assume that doubly charged nickel vacancies and electron
holes are the dominating defects in Ni1-yO under oxidising conditions. At 1245C
and pO2 = 1 atm we know the following for the compound:
iv) Calculate the charge mobility of the nickel vacancies and the ionic
conductivity under the conditions referred to above.
3. Ca-stabilised ZrO2 (CSZ)
We shall here consider a densely sintered ZrO2 doped with 15 mol% CaO
(Zr0.85Ca0.15O1.85).
6.29
6. Electrical conductivity
i) Assume that the oxide contains doubly charged oxygen vacancies compensating
the Ca dopant. What are the site-fractions of dopants and of oxygen vacancies?
ii) Derive equations showing how the minority concentrations of defect
electrons and electron holes vary with pO2 in this oxide under the given conditions.
iii) The conductivity of this oxide is independent of pO2 from oxidising to
very reducing conditions. What can we deduce from this?
iv) Simpson and Carter (J. Am. Ceram. Soc. 49 (1966) 139) measured the
self diffusion coefficient for oxygen in Zr0.85Ca0.15O1.85 and found it to be DO =
2.0*10-7 cm2/s at 1100C. Calculate the electrical conductivity based on this.
v) Find also the diffusion coeffeicient and charge mobility for the oxygen
vacancies.
6.30
7. Electrochemical transport
7. Electrochemical transport
Electrochemical potential
In the chapter on diffusion we learned that random diffusion is driven by thermal
energy and that it in crystalline solids requires defects. Moreover, we learned that a net flux
of a species can result from a gradient in its chemical potential. In the chapter on electrical
conductivity we furthermore saw that a net flux of charged particles would result in a
gradient of electrical potential.
In ionic media (materials, liquids, and solutions) the chemical and electrical potentials
act simultaneously, and it is often convenient to combine them into an electrochemical
potential. For the species i we have
(7.1)
i= i+ z i e
This is true whether the species is real or a defect, but in general we shall onwards deal
mainly with real species and real charges zi.
The electrochemical potential gradient is accordingly, in the one-dimensional case,
d i
dx
d i
dx
+ zi e
d
dx
(7.2)
The combination of chemical and electrical potentials and potential gradients forms the basis
for the treatment of all mass transport processes involving charged species (ions) in ionic
solids, and is the theme of this chapter. The theory is termed Wagner-type after Carl Wagner,
who was the first to derive it, originally to describe oxidation of metals.
7.1
7. Electrochemical transport
ji =
- D i ci d Pi - D i ci d i - D i ci d i
d
=
=
[
+ zi e
]
kT dx
kT dx
kT
dx
dx
(7.3)
Via the Nernst-Einstein equation we may substitute conductivitiy for the random diffusivity
and obtain the following alternative expression for the flux density:
ji =
- i d i
d
[
+ zi e
]
2
dx
dx
( z i e)
(7.4)
If the species is charged, the flux density for i gives rise to a partial current density ii:
ii = z i e j i = -
i [ d i + e d ]
zi
zi e
dx
(7.5)
dx
The net current density in the sample is obtained by summing the partial current densities
over all the species k:
itot = k z k e j k = - k
k [ d k +
zk e
dx
zk e
d
]
dx
(7.6)
The sample is next assumed to be connected to an external electric circuit. This may be real,
with electrodes and wires and an electrical supply or load, or it may be absent, in which case
we say that we have an open circuit. Since the sample makes a series connection with the
external circuit, the currect in the sample must be the same as in the external circuit.
For a bare sample, a gas permeation membrane, or for an open-circuit fuel cell or other
electrochemical device, the total current is zero. For fuel cells or electrolysers in operation,
on the other hand, the current is non-zero.
By using the definition of total conductivity, tot = k k , and the definition of
transport number, t k =
k
k
=
, we obtain the following very important expression:
tot k k
d
i
t d k
= - tot - k k
dx
tot
z k e dx
(7.7)
7.2
7. Electrochemical transport
This relates the electrical potential gradient to the total (net, external) current density, the
total conductivity, and the transport number and chemical potential gradient of all charge
carriers.
The chemical potentials of charged species are not well-defined, and we need to
represent them instead by chemical potentials of neutral species. For this purpose we may
assume equilibria between neutral and charged species and electrons, i.e. in the
electrochemical red-ox- reaction
(7.8)
S = S z + ze
where S is a neutral chemical entity and z may be positive or negative. The equilibrium
condition for this is expressed in terms of the chemical potentials of products and reactants:
d S z + zd e- d S = 0
(7.9)
d S z = d S zd e-
We insert this for all ionic species n in the expression for the electrical potential gradient. The
entry (among k) for electrons is left unsubstituted. By using k t k = 1 we obtain
d
i
t d n 1 d e
= - tot - n n
+
dx
e dx
tot
z n e dx
(7.11)
for which the chemical potentials now refer to the neutral forms of each carrier.
7.3
7. Electrochemical transport
II
II
II
II
1
i tot
tn
I d = I tot dx - n I z n e d n + I e d e
(7.12)
1
II I = itot dx - n t n d n + ( e, II e, I )
e
e
I tot
I zn
(7.13)
II
II
We further assume that the voltage is measured on each side using the same inert metal, e.g.
Pt. This eliminates the difference between chemical potentials of electrons on the two sides,
and the voltage measured between the two sides is
II
II
U II I = II I = i tot dx - n t n d n
e
I tot
I zn
(7.14)
II
U II I = - n
I
tn
d n
zn e
(7.15)
We shall later see how this is used to calculate transport numbers based on open circuit
voltage measurements of cells exposed to a well-defined gradient in chemical activities.
Alternatively, if tn is known to be unity for ionic charge carriers this expression will yield the
open circuit voltage of a fuel cell or a galvanic sensor.
If current is drawn from the sample, as in a fuel cell or battery under load, we need to
know how itot varies with x. If we assume that it is constant, we get
U II I = itot
tot
II
II
- n t n d n = itot rtot - n t n d n
e
e
I zn
I zn
7.4
(7.16)
7. Electrochemical transport
where X is the thickness of the sample and rtot is the area-specific resistance of the sample.
Alternatively, itotrtot (current density and area specific resistance) may be replaced by RI,
(sample resistance and current):
II
U II I = IR - n
I
tn
d n
zn e
(7.17)
The voltage of a fuel cell or battery is thus composed of a thermodynamic part that
arises from the chemical gradient, and the well-known IR-term which arises from the
limiting kinetics (transport) in the sample. The letters IR may be taken to reflect current
times resistance, or to mean internal resistance. We will later see that the direction of the
current ends up such that the IR-term reduces the voltage over the fuel cell or battery.
ji =
d i
d k
t i itot
i
[
z i k t k
]
2
z i e ( z i e)
dx
z k dx
(7.18)
The first term in the right hand side simply says what we expect: If there is a net current, a
flux density of species i will be set up proportional to the total current density and the
transport number ti, divided by the species charge.
The equation above is a rather general expression that we can use to calculate flux
densities of one charged species in the company of many other species. However, it reflects
the flux density and gradients and properties at a particular point through the membrane. The
gradients will adjust according to the varying materials properties so as to maintain a constant
flux density everywhere what we call steady state. In order to implement this, we integrate
7.5
7. Electrochemical transport
the flux density expression over the thickness of the membrane and require that the flux
density remains constant:
II
II
j idx = j i X =
II
t i itot
i
t
dx
[d i z i k k d k ]
2
zi e
(
z
e
)
zk
i
I
(7.19)
If the transport number ti and the total current density itot can be taken as constant through the
membrane we further obtain:
ji X =
II
t i i tot
i
t
X
[d i z i k k d k ]
2
zi e
(
)
z
e
zk
i
I
(7.20)
jO 2 =
je =
- O 2
4e
- e
e
d O2
d
2e
]
dx
dx
(7.21)
d e
d
e
]
dx
dx
(7.22)
By summing up the currents from the two we obtain (as in Eq. 7.7) the expression for the
potential gradient:
7.6
7. Electrochemical transport
d
i
t 2 d O 2 t e d e
= - tot + O
+
dx
e dx
tot 2e dx
(7.23)
In order to relate the potantial gradient to the chemical potential of neutral species we
introduce electrochemical equilibria between neutral and charged species of oxygen:
(7.24)
2
O 2 (g) + 4 e = 2 O
O 2(g)
+ 4d
= 2d
(7.25)
(7.26)
and knowing that the sum of transport numbers for oxygen ions and electrons in the oxide we
have defined is unity, we further get
2 d O ( g )
d
1 d e
2
= - i tot + t O
+
dx
e dx
tot 4e dx
(7.27)
whereby we have exemplified Eq. 7.11. We now integrate over the thickness of the sample to
obtain the voltage over it, eliminating the chemical potential of electrons (as in Eq. 7.14 and
onwards):
1
= IR + t O 2 d O2 ( g )
4e I
II
U II I
(7.28)
Using
O ( g ) = O0 ( g ) + kT ln pO
2
(7.29)
7.7
7. Electrochemical transport
d O2 ( g ) = kTd ln p O2
(7.30)
we get
II
U II I
kT
= IR +
t O 2 d ln pO2
4e I
(7.31)
U II I = IR + t O 2
II
kT pO2
ln
= IR + t O 2 E N
4e pOI 2
(7.32)
t O 2 =
U II I
EN
(7.33)
U II Ref =
pOII
kT
ln Ref2
4e p O2
(7.34)
7.8
7. Electrochemical transport
Fuel cell
If we apply a solid electrolyte with t O 2 = 1 and expose it to a gradient in oxygen partial
pressure, we get
U II I
II
kT pO2
= IR +
ln
= IR + E N
4e pOI 2
(7.35)
This would be the voltage of a fuel cell. With I = 0, we obtain the OCV of the cell to be the
Nernst voltage, of course. We furthermore see that the cell voltage diminishes due to the IR
loss at increasing current, I. A plot of U vs I gives R as the slope.
The power delivered by the cell is the current multiplied with the voltage;
P = UI = I 2 R + IE N
(7.36)
and it goes through a maximum at U=EN/2. Under these conditions, the power is split halfhalf between external work and internal loss, so the electrical efficiency is 50%. The
efficiency increases as the power and current decrease, and commonly fuel cells are designed
to run at ca. U=2/3 EN.
jO 2 =
t O 2 itot O 2 t e d O 2
t 2 i tot O 2 t e d O2 ( g )
d e
[
2
]= O
2
2e
dx
dx
2e
dx
4e
8e 2
(7.37)
The steady state restriction is introduced by integrating over the thickness of the membrane:
7.9
7. Electrochemical transport
II
II
j O 2 dx = j O 2 X =
II
2 t e
t O 2 itot
dx O 2 d O2 ( g )
2e
8e
I
II
jO 2
(7.38)
II
kT
i
= tot t O 2 dx 2 O 2 t e d ln pO2
2eX I
8e X I
(7.39)
If itot is not zero, i.e., some external current is drawn, then the first term says that the oxygen
ion flux is proportional to that current and to the oxygen ion transport number. If the latter is
independent of x, then the integral is straightforward. The last term adds to the flux; it is
driven by an oxygen pressure gradient, and it is proportional to the oxygen ion conductivity
and the transport number of electrons. This term describes permeation of oxygen ions,
requiring transport of both oxygen ions and electrons. Moreover, it is inversely proportional
to the thickness of the sample, X. Under open circuit conditions we may simplify to get:
II
jO 2 =
II
kT
kT
O 2 t e d ln pO2 = 2 tot t O 2 t e d ln pO2
2
8e X I
8e X I
(7.40)
We now apply our derived expression for the example case of an oxygen-deficient oxide
MO1-x. This contains oxygen vacancies compensated by electrons, and we have earlier shown
that the defect concentrations are given by n = 2[vO ] = (4 K v/ )1 / 3 pO12 / 6 . Since the mobility of
O
electrons is usually much higher than that of ionic defect, and with similar concentrations of
electrons and oxygen vacancies, we may assume that the material is mainly an n-type
electronic conductor. Therefore t e 1 . Furthermore, O 2 = v = 2e[vO ]uv = 0 pO12 / 6 ,
O
jO 2 =
kT 0
8e 2 X
II
1 / 6
O2
(7.41)
d ln pO2
We substitute d ln pO2 =
1
dpO2 and obtain
pO2
7.10
7. Electrochemical transport
jO 2
II
kT
= 2 0
8e X
7 / 6
O2
(7.42)
d pO2
jO 2 =
6kT 0
( p OII2 ) 1 / 6 ( p OI 2 ) 1 / 6
2
8e X
(7.43)
One may note that we here have utilised the Wagner approach for electrochemical transport
of two defect species and combined it with the solution of the defect chemical equilibria to
obtain an accurate description of the flux of oxygen ions. Numerical values for the flux of
course require that one knows 0 , which contains many parameters relating to formation and
mobility of the defects. In the right-hand parenthesis the negative exponent 1/6 has the
consequence that the smaller of the two partial pressures of oxygen has the biggest influence.
The larger partial pressure makes a comparatively smaller difference and may be neglected to
a first approximation for large gradients. For instance, consider a gas separation membrane
exposed to air on the one side and a very reducing gas (e.g. methane) on the other. The flux
would then be essentially proportional to ( pOred2 ) 1 / 6 while it would largely be independent of
pOox2 . It would thus not help much to compress the air to 100 atm, while removal of oxidised
products like H2O on the reducing side would lower pOred2 and have a large effect.
We have seen above that the transport coefficient term that enters in the flux of oxygen
ions in a mixed oxygen ion conductor is a product which we may call the ambipolar
conductivity
amb = tot t O t e = O t e =
2-
2-
O e
O +e
2-
2-
(7.44)
-
This reflects the fact that so-called ambipolar conduction of two species is necessary (in this
case for charge compensation).
7.11
7. Electrochemical transport
The corresponding terms transformed to diffusivities are called the ambipolar diffusion
coefficients for oxygen:
(7.45)
D amb = DO 2- t e-
Such terms, like partial conductivity and self-diffusion, reflect thermal random mobility
as before, but are restricted by the slower of the two ambipolar diffusing species (ions or
electrons).
Many other cases of combined transport coefficients are in use, e.g. the combined
(additive) transport of oxygen and metal ions commonly that we shall address later (and
exemplify by the high temperature oxidation of metals), the combination of two diffusivities
involved in interdiffusion (mixing) processes, and the mass transport in creep being rate
limited by the smallest out of cation and anion diffusivities in a binary compound. As some
of these sometimes are referred to as ambipolar or chemical diffusivities, we want to stress
the above simple definition of ambipolar transport coefficients as relevant for membrane
applications using mixed conductors.
~ dci
ji = - D
i
dx
(7.46)
Unlike our previous treatment in this chapter (which uses Dr, B, u or ), Ficks 1st law and
the chemical diffucion coefficient says in general nothing directly about the self-diffusion
coefficient. Neither is, as we shall see, the chemical diffusion coefficient in general a constant
in any sense. Only when the species i is an effectively neutral particle that moves
independently of other species will the chemical diffusion coefficient be identical to the selfdiffusion coefficient. If the flux and the concentration gradient are in terms of an isotope
7.12
7. Electrochemical transport
(tracer), the coefficient is the tracer diffusion coefficient, which is almost equal to the selfdiffusion coefficient:
(7.47)
Dt = f Dr
where the correlation factor f is a number usually equal to or somewhat smaller than unity,
depending on structure and diffusion mechanism, reflecting the difference of transport
between distinguishable and indistinguishable diffusing species. For vacancy diffusion, f
usually falls around 0.75 depending on structure.
Chemical diffusion is sometimes referred to as diffusion under a chemical potential or
concentration gradient. However, we have seen that self-diffusion properly describes this, and
that all diffusion phenomena we consider here are merely minor perturbations of thermally
induced self-diffusion. Why then do we need to be concerned with chemical diffusivity? We
use it whenever we want to or have to relate our flux to the concentration gradient. Under
steady state processes, our flux vs force equivalents provides sufficient description, and can
be integrated to yield steady state fluxes through membranes, as we have shown; Thus, for
this purpose, the self-diffusion coefficients and conductivities are useful since they behave in
a simple manner vs defect concentrations, and since we in general know the chemical and/or
electrical potential gradents applied. However, if we know the concentration gradients and
how the chemical diffusivity varies with concentration, we can use chemical diffusivity.
In a transient experiment, however, parts of the material take up or give away matter
(leading to changes in the concentration of species with time). This is sometimes referred to
as a capacitive effect and is expressed in Fick's 2nd law:
dc
dj
=- =
dt
dx
dc
)
dx
dx
d( D
(7.48)
2
dc
dj
d c
=- =D 2
dt
dx
dx
(7.49)
7. Electrochemical transport
diffusion coefficient in case of a transient in isotope composition (whether measured
radiographically, mass spectrometrically, thermogravimetrically or electrically) or chemical
diffusion coefficients in case of a transient in chemical composition (measured
thermogravimetrically or electrically or by other means). The latter is usually executed by
stepping the oxygen activity by gas composition, total pressure, or electrochemical means.
It may be noted that certain tracer experiments actually yield an interdiffusion
coefficient between two tracers. The diffusion coefficient for oxygen in practice cannot and
need not distinguish between the different isotopes in use (usually 18O and 16O, more rarely
17
O). However, for protons, the difference to deuteron (or triton) transport coefficients is
considerable (e.g. a factor of 2) and should be taken into account when interpreting the
results.
The multitude of transport coefficients collected can thus be divided into self-diffusion
types (total or partial conductivities and mobilities obtained from equilibrium electrical
measurements, ambipolar or self-diffusion data from steady state flux measurements through
membranes), tracer-diffusivities, and chemical diffusivities from transient measurements. All
but the last are fairly easily interrelated through definitions, the Nernst-Einstein relation, and
the correlation factor. However, we need to look more closely at the chemical diffusion
coefficient. We will do this next by a specific example, namely within the framework of
oxygen ion and electron transport that we have restricted ourselves to at this stage.
In a system with oxygen vacancies as the only oxygen point defects, the concentration
of oxygen ions is the concentration of oxygen sites minus the concentration of vacancies:
(7.50)
cO 2 - =cO 2 -- cvO..
such that
dcv ..
dcO 2 =- O
dx
dx
(7.51)
7.14
7. Electrochemical transport
We manipulate this in the following manner:
dcvO.. d ln pO 2
dcv ..O d ln pO 2
dcO 2 ==dx
dx d ln pO 2
d ln pO 2 dx
(7.52)
and insert the result into the flux equation for oxygen ions, Eq. 7.37 (open circuit, itot = 0):
jO 2 =
O te d O (g )
8e
dx
O t e kT d ln pO
2
8e
dx
O t e kT d ln pO dcO
2
8e
dcv
dx
(7.53)
We have thus derived an expression for the flux that relates it to the concentration gradient of
oxygen ions. We may now therefore compare this expression with that of the oxygen flux in
terms of Fick's 1st law involving the concentration gradient in oxygen ions and the chemical
diffusion coefficient:
~ 2 - dcO 2 jO 2- = - D
O
dx
(7.54)
and we obtain
O 2 t e kT d ln pO2
~
DO 2 - =
dcv
8e 2
(7.55)
Further, by using the Nernst-Einstein-relation and inserting the relation DO 2 cO 2- = DV O.. cV O.. we
get
D c
D
Dv t e
O
~ 2 - = DO 2 cO 2 t e d ln pO2 = vO vO t e d ln p O2 = vO t e d ln pO2 =
DO
d
ln
c v
2
dc v
2
dcv
2 d ln cv
O
O
O
O
2
d ln p O2
(7.56)
7.15
7. Electrochemical transport
(7.57)
E = E N - I Ri - ( c + a )
7.16
7. Electrochemical transport
where I is the current, EN is the open circuit (Nernst) voltage, Ri is the ionic electrolyte
resistance, and c and a are the cathode and anode overpotentials, respectively. (The
conventions and practices on measuring and reporting the sign of various voltages may vary,
so that the equation may take various forms with respect to signs).
Similarly, the current of matter through an electrodeless gas separation, mixed
conducting membrane will be limited by ionic conduction, the surface kinetics, and by the
internal electronic resistance, the latter replacing the external circuit in the previous case.
However, we may need to express the loss of potential in terms of the chemical
potential of oxygen.
If we look at the system without any gradients or currents applied, there will be
equilibrium at both the electrodes and the surfaces. Depending on the application and the side
of the membrane, the electrochemical processes in operation on an oxygen ion conductor may
be:
(7.58)
2O 2 (g) + 4 e = 2 O
and
(7.59)
2H 2 (g) + O = H 2 O(g) + 2 e
In all of these the electrons can be supplied to or from an electrode or a mixed conductor
surface. Thus, basically, the reactions are the same in the electrolyte/electrode as in the mixed
conductor, although the physical presence of the electrode can promote (by catalysis) or
obstruct the reaction.
Most surface reactions are believed to be made up of a sequence of reactions steps.
Some possible step reactions in the case of the reduction of oxygen gas are given by the
following reaction sequence;
O 2 (g) = O 2,ads
O 2,ads + e = O 2,ads
O 2,ads + e = 2 O ads
2O ads + e = O ads
2-
2-
O ads = OO
7.17
7. Electrochemical transport
The last of these steps may involve bulk defects (e.g. a vacancy that the adsorbed ion may
occupy) but we will not attempt here an analysis of surface or electrode kinetics in terms of
defect chemistry. Determination of the nature of the rate limiting step of a reaction is a
difficult task. Studies of surface exchange kinetics using 18O or by relaxation methods can
give information about the temperature and pressure dependence of the rate determining step.
However, it should be pointed out that information from such methods do not give direct
information about the exact mechanism and which intermediate species that are involved.
Each step will have its own specific rate and resistance to the mass transport across the
membrane. In lack of detailed knowledge of the rate limiting step(s) we normally have to
relate to the overall reaction.
At equilibrium the overall reaction takes place forth and back at equal rates r0 as a
result of thermal energy. The activation process forth and back may not be symmetrical in
this case, but that will not affect us at present. Similarly to the random diffusion coefficient
connected with these thermal disorders in bulk (Dr = 1/6 r0 s2) we introduce an exchange
coefficient for interfaces; ki = r0si which says something about the thermal fluctuations of the
reaction across the interface of thickness si. We can alternatively express this as an exchange
flux density
(7.60)
j 0 = k i c = r0 si c
(7.61)
i0 = j0 ne = k i nec = r0 si nec
where n is the number of electrons involved in the reaction. If we now apply a small force
over the interface, we get a perturbation and a net flux, given by
j = -c k i
dP
kT
(7.62)
where c is the volume concentration of species in the interface, and dP is the potential step
over the interface. The net current density is
7.18
7. Electrochemical transport
i = jne = cne k i
dP
dP
= i0
kT
kT
(7.63)
kT i
ne i0
(7.64)
This is the linear or ohmic (small overpotential) form of the kinetics expression, as we expect
from our specification of small forces and perturbations. We note that the interface thickness
si has dropped out and is mostly not an interesting parameter at this level.
By measuring the overvoltage vs current density by voltammetric DC or impedance
spectroscopy AC methods, we may obtain the charge transfer resistance
Re =
kT
kT
= 2 2
ne i0 n e k i c
(7.65)
The overpotential enters directly into the sum of voltages over the fuel cell or electrolyser
circuit, and we can solve the system in terms of current and voltage. As the force and
perturbations of the reaction becomes larger, the full so-called Butler-Volmer exponential
dependence of current and flux on overpotential (and force) must be taken into account:
kT
i
ln( )
ne
i0
(7.66)
7.19
7. Electrochemical transport
We resort first to chemical diffusion and assume that the interface process is one which
can be described as species flowing from a high to a lower concentration. Then
- D(
dc
) = -k cinterface
dx bulk
(7.67)
, and D as well as k must refer to chemical coefficients. This requirement can be included in
the solution of Fick's 2nd law for transients, and D and k can in principle be found
independently by fitting transient data to this solution. The k's obtained in this way can be
used to predict steady state fluxes if they are measured as a function of e.g. pO 2 or at the pO 2
extremes which the membrane is going to be subjected to. Otherwise, the k values are not
immediately easy to relate to other data.
In a tracer experiment, we obtain Dt which is close to the self diffusivity (Dt = f Dr) and
then the exchange coefficient should also reflect the tracer exchange coefficient kt . ki and
this should be comparable to that calculated as above from measurements of Re; the
difference should say to what extent the presence of an electrode has changed the kinetics of
the rate limiting step of the surface reaction.
In order to analyse the effect of interface kinetics on fluxes through a membrane under
steady state conditions, we may take an approach similar to that for fuel cells and
electrolysers. Principally, no electrodes are involved, and voltages are therefore not directly
appropriate as in the other cases. Instead we may sum up chemical potential changes of
oxygen (or hydrogen) through the membrane and its interfaces. We start out with the
membrane interior, and from earlier we have for a mixed oxygen ion and electron conductor
of thickness L (earlier also referred to as X or x):
1 O2 t e
d O2 ( g )
L I 8e 2
II
jO 2 =
(7.68)
If we introduce ambipolar diffusivities and conductivities and assume these constant, the
equation may be integrated to obtain
jO 2 =
D amb cO 2- O 2 = amb O 2
2kT
L
8 e2 L
(7.69)
7.20
7. Electrochemical transport
and, by rearrangement, the chemical potential drop over the bulk of the membrane is
O2 =
- 2 jO 2 kTL
- jO 2 8 e 2 L
(7.70)
amb
D amb cO 2-
This drop over the membrane bulk is given by the external full gradient minus the losses at
the interfaces. In the present terms, the flux through an interface is, from Eqs. 7.62 and 7.65
given by
jO 2- =
- cO 2- k i Ointerface
2
2kT
- interface
O2
(7.71)
8 e2 R e
interface
=
O2
- 2 jO 2 kT
k i cO
2-
(7.72)
=- jO 2 8 e2 Re
Now, we sum up the total gradient in chemical potential of oxygen over the bulk of the
membrane and the two interfaces
(7.73)
Otot2 = Obulk
+ 2 interface
O2
2
and after insertion of the above expressions for bulk and interface, we can solve for the flux
of oxygen ions through a membrane with two equally contributing interfaces under the
conditions given:
(7.74)
7.21
7. Electrochemical transport
This result simply states that the flux is given by the driving force (which easily transforms to
Nernst voltage) divided by the total ambipolar resistance, in our case namely that of the bulk
and twice that of one interface.
In lack of full understanding of the relations involved in interface kinetics of
membranes, we may refer to a critical thickness; the thickness of the bulk membrane where
the interface and bulk impose equally large restrictions on the flux. In any membrane
application, going much below this thickness is not of much use. We see from the above that
this thickness is Lcrit=D/k for the thickness of the part of a membrane relating to one
interface, or Lcrit =2D/k for thickness of a membrane with two equally limiting interfaces,
typically the two surfaces.
Experimentally the value of Lcrit has been found to be of the order of 100 m for many
polished samples of different fluorite and perovskite materials, indicating that there is a close
relationship between e.g. vacancy concentration in bulk and the surface exchange. Roughing
the surface should increase relatively the surface flux, and this has been shown to work.
General expressions
We have in the previous section considered ionic transport by oxygen ions (anions)
only. This was done for simplicity, because it represents many real applications, and because
oxygen ions relate most directly to the oxygen (non-metal) activity which we may control
most conveniently.
Now we will include cation (metal ion) transport. Some systems have predominant
transport of cations, and we thus need to see how this relates this to the non-metal activity.
Some systems have both anion and cation transport so we need to take both into account. In
some cases we may not know which one is dominating and as we shall see it is not always
easy or necessary to distinguish them.
In an inorganic compound the total ionic current density, iion, is given by the sum of the
current densities of anions, ian, and the cations, icat. From our previously derived expressions
for partial current densities we thus obtain
an d an cat d cat
+
z an e dx
z cat e dx
7.22
(7.75)
7. Electrochemical transport
where zan and zcat represent the valences and thus ionic charges of the anions and cations,
respectively.
In order to relate the current density of the cations and anions let us consider the
equilibrium involving the formation of the compound MaOb from its ions:
(7.76)
aM zcat + bO z an = M a Ob
(7.77)
ad M z cat + bd O z an = d M aOb ( s )
The fact that d M a Ob ( s ) = 0 arises from MaOb(s) being a pure condensed phase. It may further be
noted that
(7.78)
az cat = bz an
d cat
b d an z cat d an
=
=
dx
a dx
z an dx
(7.79)
From this we note that the cation and anion chemical potential gradients are the negative of
each other in a binary compound, a relation that is used extensively, mathematically or by
intuition.
By adding z cat e
d
to both sides of the equation and rearranging on the right hand side
dx
we obtain
d cat z cat d an
=
dx
z an dx
(7.80)
7.23
7. Electrochemical transport
By combining Eqs 7.75. and 7.80 the total ionic conductivity may be expressed by
iion =
ion d cat
z cat e dx
ion d an
z an e dx
(7.81)
Thus, the ionic current density can be expressed in terms of the ionic conductivity (sum of
cationic and anionic conductivities) and the gradient in the chemical potential of cations or
anions.
From this it becomes clear that the derivations of oxygen ion fluxes and currents in
various situations done earlier could have been done for the case of ions (sum of cations and
anions) simply by inserting ion instead of O 2 , but with for instance the oxygen activity
gradient as the driving force.
Membrane walk-out
Metal ions are transported from the low oxygen pressure to the high oxygen pressure
side of the membrane while metal vacancies are transported in the opposite direction. During
the process oxygen is liberated and oxygen and metal sites are annihilated at the low oxygen
pressure side, e.g. for a metal deficient oxide M1-O with doubly charged metal vacancies:
2'
1
OO + VM + 2h. = 2 O2(g)
(7.82)
while equivalent number of lattice sites are formed at the high oxygen pressure side:
1
2
2'
O2(g) = OO + VM + 2h.
(7.83)
In this way the oxide membrane is actually moving in laboratory space in the direction of the
higher oxygen potential.
From this we see that even a very minor transport of metal ions may be detrimental to
the operation of an oxygen separation membrane or fuel cell: Over the many years of
operation the electrolyte or membrane will simply walk out of its housing, towards the high
oxygen pressure.
7.24
7. Electrochemical transport
Demixing of oxide solid solutions
An additional important aspect of diffusional transport of metal ions in a chemical
potential gradient is that in a homogeneous crystal of an oxide solution, e.g. (A,B)1-O such
as for instance (Co,Ni)1-O, a demixing process begins to take place. Both cations move by
vacancy diffusion. When one of the cations in (A,B)1-O, e.g. A2+, has a higher mobility than
the other, the A2+ ions will move faster to the side of the higher oxygen potential. The solid
solution is enriched in AO at this side and thus becomes kinetically demixed. After extended
time a steady state concentration profile is reached. This can be formulated through the use of
the appropriate transport equations and by taking into account the conditions of
electroneutrality and of local thermodynamic equilibrium. The process is schematically
illustrated below.
Figure 7-1. Schematic illustration of kinetic demixing of the oxide solid solution (A,B)1-O in a oxygen
potential gradient.
7.25
7. Electrochemical transport
dx
1
= k *p
dt
x
(7.84)
x 2 = 2k *p t + C 0 = k p t + C 0
(7.85)
where k p ( = 2k *p ) and k *p are expressions for the rate constants and Co is an integration
constant. The oxide thickness may alternatively be measured in terms of the oxygen uptake
by the metal specimen during oxidation by using, for instance, thermogravimetry. The oxide
thickness or oxygen uptake is parabolic with time, and the oxidation is termed parabolic.
Accordingly k p and k *p are termed parabolic rate constants. The figure below illustrates the
variation in the oxide thickness or oxygen uptake as a function of time.
Figure 7-2. The variation in oxide thickness as a function of time (x vs t and x2 vs t) during parabolic
oxidation of metals
7.26
7. Electrochemical transport
Carl Wagner developed the theory for parabolic oxidation assuming that the reaction is
governed by lattice diffusion of the reacting ions (metal and oxygen ions) or transport of
electrons. The important aspect of the Wagner oxidation theory is that the parabolic rate
constant is expressed in terms of independently measurable properties, i.e. the electronic and
ionic conductivity of the oxide or alternatively in terms of the self-diffusion coefficients of
the reacting ions. This thus provides an interpretation of the reaction mechanism and a
theoretical basis for changing and improving the oxidation resistance of metals and alloys.
The Wagner theory has been one of the most important contributions to our understanding of
the high temperature oxidation of metals.
The basic assumption of the theory is that the lattice diffusion of the reacting ions or the
transport of electrons through the dense scale is rate-determining for the overall reaction.
Lattice diffusion is assumed to take place because of the presence of point defects, and the
migrating species may alternatively be considered to constitute lattice and electronic defects.
Wagner further assumed that ions and electrons migrate independently of each other and that
local equilibria exist within the oxide. Such transport processes through a dense, single-phase
scale growing by lattice diffusion is illustrated in the figure below. One part (a) illustrates the
transport of metal and oxygen ions and electrons, while the other (b) illustrates the transport
processes when the predominant diffusion processes involve metal vacancies and interstitials.
Figure 7-3. Transport processes during growth of a dense, single phase scale growing by lattice
diffusion. a) transport processes illustrated by transport of reacting ions, b) transport processes illustrated by
transport of diffusion of important point defects (in this case assumed to be metal vacancies and interstitials).
7.27
7. Electrochemical transport
As diffusion through the scale is rate-determining, reactions at the phase boundaries are
considered to be rapid, and it is assumed that thermodynamic equilibria are established
between the oxide and oxygen gas at the oxide/oxygen interface and between the metal and
the oxide at the metal/oxide interface.
The overall "driving energy" of the reaction is the Gibbs free energy change associated
with the formation of the oxide, e.g. MaOb, from the metal M and the oxygen gas.
Correspondingly, a gradient in the partial pressure (activity) of oxygen exists across the scale
from the partial pressure of oxygen in the ambient atmosphere at the outer surface to the
partial pressure of oxygen at the metal/oxide interface. The latter partial pressure is the
decomposition pressure of the oxide in equilibrium with its metal.
The main driving force for the transport through the plane is the chemical potential
gradient. But when one considers the transport of the reacting ions and of the electrons
through the scale, it may also be noted that the mobilities of the cations, anions and electrons
are not equal. Because of this difference, a separation of charges takes place in a growing
scale. This creates a space charge (diffusion potential) that opposes a further separation of
charges, and a stationary state is reached for which no net electric current flows through the
scale. In describing the transport of ions and electrons through the scale, it is thus necessary
to take into account both the transport due to the gradient of the chemical potential and that
due to the gradient of the electrical potential, i.e. the electrochemical potential gradient. Thus,
the treatment of the transport through the oxide scale is identical to that done earlier for
electrochemical transport through a mixed conducting oxide.
The net ionic current is from earlier given by
iion =
ion d an
z an e dx
ion d O
2e
(7.86)
dx
iel =
el d el
1e dx
el d el
e
(7.87)
dx
As no net current flows through the scale during the scale growth, then itot = iion + iel = 0.
By solving this system in the usual manner (solving for the electrical potential gradient,
inserting into the expression for the ionic current density and inserting the equilibrium
condition between oxygen ions, oxygen molecules and electrons) we obtain
7.28
7. Electrochemical transport
iion =
ion t el d O ( g )
2
4e
dx
ion t el kT d ln p O
(7.88)
4e
dx
The rate of oxygen ions reacting per unit area is obtained by dividing the ionic current density
by the charge of the oxygen ions. Moreover, the rate dn/dt of oxide molecular units MaOb
formed per unit area is obtained by further division by b:
i
t d O2 ( g )
t kT d ln pO2
dn 1 dnO
=
= ion = ion2 el
= ion 2el
dt b dt
2eb
dx
dx
8e b
8e b
(7.89)
Similar to the integrations we have done earlier over the thickness of the oxide, we now
integrate over the instantaneous thickness x and from the ambient oxygen pressure, pO0 2 , to
dn
pO
1
dn kT 2
=
t
d
ln
p
ion
el
O
2
dt 8e 2 b pi
x
O
2
(7.90)
We have chosen to organize the equation and direction of integration such that the rate is
positive if p O0 2 > p Oi 2 . In general the directionalities and signs of the fluxes and processes that
go on pose a problem for us. One may most often neglect this issue because the output,
namely growth (increase in amount of oxide or thickness of scale) is intuitively positive.
However, it may not be just that easy, since there is nothing preventing scale reduction
(reverse of growth) if the outer oxygen partial pressure is smaller than the activity at the
metal/oxide.
Many readers will probably wonder about another small apparent problem: We
integrate over the scale to take into account the steady state condition, namely that the fluxes
are constant throughout the scale. Still, the scale grows, so how can it be steady state? The
answer is that it is steady state (constant fluxes) for that moment of time with that
instantaneous thickness of the scale. So our integration is for a given moment in time. At
another moment the integration will give a different flux.
7.29
7. Electrochemical transport
kT p O 2
one form of the parabolic rate constant and is in the following written kt:
o
dn
1
= kt
dt
x
(7.91)
kt thus represents a time-independent coefficient, while the flux and thickness of the scale
vary according to the parabolic relation above.
dn
As written in this form dt
cm2sec and x in cm. The derivation of Wagner's equation that we here have done for growth
of oxide scales, may be applied to many other gas-metal reactions.
It may be noted that although the total particle current is equal to the rate of growth of
the scale, the rate-determining process may either be diffusion of ions or transport of
electrons depending on the properties of the scales. If the scale is an electronic conductor, tel
1, the diffusion of the reacting ions through the scale is rate-determining, while the
transport of electrons through the scales is rate-determining if the scale is an ionic conductor,
tion 1. The oxygen pressure dependence of the reaction will depend on whether the scale is
an ionic or electronic conductor and in the latter case on the type of nonstoichiometry that
prevails in the oxide scale.
In our derivation, the growth rate is proportional to iontel = tottiontel = (an + cat) tel.
However, the conductivity of the metal and oxygen ions in MaOb can be expressed in terms
of the self-diffusion coefficients of the metal and oxygen ions in MaOb through the Nernst Einstein relation (zi2e2Dici = kT i). Let us insert this relationship in Eq.7.89 and assume that
the oxide is an electronic conductor (tel 1). Let us furthermore take into account that
concentrations of metal ions (cations), cM, and oxygen ions, cO, in an oxide MaOb with
relatively small deviation from stoichiometry are related through cM/cO = a/b = |zan|/zcat. The
expressions for the flux through the growing scale (expressed in terms of the number of
MaOb units per cm2sec.) and for kt then become
d ln pO2
d ln pO2
c z2 z
c z
dn
= O an ( cat D M + DO )
= O ( cat D M + DO )
dt
dx
dx
8b z an
2b 2
7.30
(7.92)
7. Electrochemical transport
pO
1
dn c 0 2 z cat
1
=
(
D M + DO )d ln pO2
= kt
dt 2b pi
2
x
x
O2
(7.93)
cO denotes the number of oxygen ions in MaOb per cm3 and DM and DO are respectively the
self-diffusion coefficients of the metal and oxygen ions in MaOb. It should be noted that DM
and DO are the self-diffusion coefficients for random diffusion of the respective ions.
In Eqs.7.92 and 7.93, kt represents the rate of formation of MaOb units per cm2 per
second for an MaOb scale of thickness 1 cm. One may as said before alternatively express the
parabolic rate constant in terms of the rate of growth of the oxide thickness
dx
1
= k *p
dt
x
(7.94)
The dimensions of kp are then cm2sec-1, and Eq.7.93 then takes the form
pO
1
dx b dn 1 2 z cat
1
(
D
D
)
d
ln
p
=
=
+
= k *p
M
O
O
2
dt c 0 dt 2 p i
2
x
x
O2
(7.95)
dx z cat
=
dt 4
1
1
D
d
ln
p
= k *p
O2
i M
x
x
pO 2
o
pO
2
(7.96)
7.31
7. Electrochemical transport
In order to integrate this we need to analyse the defect structure of the oxide. The formation
of the metal vacancies may be written as
b
b
O2 ( g ) = v M / + h + OOx
2a
a
(7.97)
'
If other defects may be neglected and the electroneutrality condition is given by [VM ] = p,
[v M / ] +1 = K V p Ob / 2 a
/
(7.98)
where K V / is the equilibrium constant for the formation of the metal vacancies (Eq. 7.97)
M
and [v M ] denotes the fraction of the metal ion sites that are vacant. On the basis of this
relation, d ln p O2 may be written
d ln pO2 =
2a ( + 1)
d ln[v M / ]
b
(7.99)
When this expression for d ln p O2 is introduced in Eq.7.96 , the expression for kp becomes
k *p = ( + 1) Dv / {[v M / ]o [v M / ]i }
(7.100)
If the concentration of metal vacancies at the oxide/oxygen interface is much larger than at
the metal/oxide interface, i.e. [v M / ] o >> [v M / ]i then
k *p = ( + 1) Dv / [v M / ]o = ( + 1) DMo
(7.101)
where DMo = Dv / [v M / ]o is the self-diffusion coefficient of the metal ions in Ma-yOb at the
M
oxide/oxygen interface. Thus under these conditions the parabolic rate constant and the selfdiffusion coefficient of M in Ma-yOb has the same temperature and oxygen pressure
dependences.
7.32
7. Electrochemical transport
Growth of oxygen-deficient MaOb-y on M
The same procedure may be used to derive the parabolic rate constant for the growth of
oxygen deficient MaOb-y on M. let us assume that the oxygen vacancies have charge , i.e.
.
OO = VO + e' + 2 O2
i
(7.102)
(7.103)
where Dv is the self-diffusion coefficient of the oxygen vacancies, [vO ]i the concentration
O
of oxygen vacancies and DOi the self-diffusion coefficient of oxygen in MaOb-y at the
i
metal/oxide interface. DO0 is the oxygen self-diffusion coefficient at 1 atm. O2, and pO is
2
the oxygen activity at the metal/oxide interface. For this case the parabolic rate constant is
independent of the ambient oxygen pressure, and the temperature dependence is given by that
of the oxygen self-diffusion coefficient in MaOb-y at the metal/oxide interface. In this respect
it may be further noted that the oxygen self-diffusion coefficient at the metal/oxide interface
is given by that the product of the self-diffusion coefficient of oxygen at constant oxygen
i
i
pressure (at 1 atm O2) and (pO )-1/2(+1) , i.e. DOi = DO0 (pO )-1/2(+1) . In these terms the
2
activation energy of the parabolic rate constant is under these conditions given by the
i
activation energy associated with DO0 and the enthalpy term associated with (p )-1/2(+1) .
O2
Similar treatments may be given for oxide scales growing by interstitial metals ions or
interstitial oxygen ions.
7. Electrochemical transport
feature of these reactions is that it is the electronic transport through the scale that is ratedetermining.
Sintering
Stages of the sintering process
Sintering is the process whereby powders and small particles agglomerate and grow
together to form a continuous polycrystalline body. As a rule of thumb the material must be
heated to 2/3 of the melting temperature to achieve considerable sintering.
The overall sintering process may conveniently be divided into different stages. The
principal stages involved is illustrated schematically below.
7.34
7. Electrochemical transport
Figure 7-4. Principle stages during the sintering process: a) particles in contact prior to sintering; b)
neck growth between particles during initial stage of sintering; c) further neck growth and formation of
continuous pore channels along grain boundaries; and d) further sintering in closing channels and formation of
closed pores at grain corners.
During an initial stage the surface roughness is decreased and the particles begin to
adhere and grow together. At the end of this stage grain growth begins to occur and grain
boundaries and metastable pore phase is established. Only modest shrinkage of a few percent
takes place during the initial stage.
During an intermediate stage grain growth continues, and the cross-sectional area of the
pore phase decreases. The final stage begins when the body achieves 90-95% theoretical
density. At this point the pore phase becomes discontinuous, and the closed pores usually
remain at grain boundaries. The final stage may involve complete removal of the remaining
pores, leading to a completely dense material. Alternatively, it may involve discontinuous
growth of the large grains at the expense of the small ones, and closed pores may as a
consequence be isolated inside the grains. If the latter process occurs, complete densification
becomes extremely difficult.
7.35
7. Electrochemical transport
Any surface atom or molecule is subjected to a resultant inward attraction because of
the unsaturated chemical bonds at the surface. The surface therefore tends to contract to the
smallest possible area. In the case of a spherical particle with radius r, the interior is subjected
to an excess force of r2p, where p is the pressure. This is counteracted by the surface tension
acting along the circumference: 2r. When equating these opposing forces, the pressure is
given by
p=
2r
r2
2
= r
(7.104)
From this it is seen that the vapour pressure of a sphere is larger the smaller the radius
and as a result a large sphere will grow at the expense of a neighbouring small sphere. This
type of grain growth in polycrystalline solids is termed Ostwald's law. Similarly the vapour
pressure of convex surfaces is larger than the vapour pressure of a concave surface and as a
result surface roughness of particles is reduced during the initial stage of sintering.
7. Electrochemical transport
2
O 2 (g) + 4 e = 2 O
(7.105)
H 2 (g) = 2 H + 2 e
(7.106)
O 2(g)
H 2(g)
+ 4d
= 2d
= 2d
(7.107)
+ 2d
(7.108)
We insert these into the flux equations for the species and utilise that the sum of all transport
numbers equals unity. We furthermore use that d = kTdlna kTdlnp, where a and p are
activity and partial pressure, respectively, and obtain
d
t 2 - kT d ln pO 2 t H + kT d ln p H 2 1 d e
i
= - tot + O
+
dx
4e
dx
2e
dx
e dx
tot
(7.109)
This can now be integrated to obtain the voltage of a cell, or inserted into a flux equation of a
species of interest in the usual manner. We now have a system with two chemical driving
forces for electrochemical transport; that of oxygen activity and that of hydrogen activity.
Other cases
We have in the present version of the treatment of electrochemical transport omitted
many cases. These include a more full coverage of proton transport in oxides and transport by
other hydrogen species and other foreign species. Moreover, we have left out the case of
solid-solid reactions. Also, creep and sintering were given more phenomenological than
defect-chemical treatments.
We have throughout assumed isothermal conditions. If we instead of (or in addition to)
chemical driving forces considered a temperature gradient, our integrations over membranes
would have been from one temperature to another. The voltage we would have obtained
7.37
7. Electrochemical transport
would be the thermoelectric voltage. It has several contributions, but in general it is related to
the transport numbers and the gradients in concentrations of defects. The variations of
concentration is in turn given by the entropy of the defects.
More on these themes will be included in future updates.
7.38
7. Electrochemical transport
Excercises
1. Electrochemical cell with stabilised zirconias (YSZ and CSZ)
a) Consider an yttria-stabilised zirconia (YSZ) material with 90 mol% ZrO2 and 10
mol% Y2O3. What is the site fraction of oxygen vacancies in the material? We often
use and refer to the YSZ as an electrolyte. What is an electrolyte? What requirements
does it put on a material?
b) We make ourselves a disk of YSZ, attach Pt electrodes to both faces, and place the
disk in a cell where we can expose the two faces to different gases. If the gases have
different partial pressures of oxygen, a voltage will be set up over the electrodes, and
we can draw current from the cell. Write equations for the two half-cell-reactions that
run when current is drawn. Write the expression for the voltage over the cell given
that the material is a purely ionic conductor. Which (high or low pO2) is the positive
pole?
c) Smith, Meszaros and Amata (J. Am. Ceram. Soc., 49 (1966) 240) found that the
permeabilities at high temperatures of Ar and N2 through calcium-stabilised zirconia
(CSZ) were immeasureable, while that of O2 was significant and proportional to
pO21/4. Suggest an explanation for these observations.
2. Electrochemical cell for determination of thermodynamic properties
a)
Indicate how you can use a YSZ cell with Pt electrodes to determine
the thermodynamic parameters for formation of a binary oxide (e.g. NiO).
b)
- (Mo)Th(s),ThF4(s)|CaF2(s)|ThF4(s),ThC2(s),C(s)(Mo) +
in order to measure the thermodynamic properties of ThC2. The Mo serves as
inert electrodes. The cell was kept under an inert atmosphere. The CaF2(s) is a solid
electrolyte a fluoride ion conductor in the temperature range of operation of the
cell (700-900C). The signs indicate the polarity of the cell. Write the reactions at
each electrode when current is drawn from the cell. Also, express the voltage of the
cell as a function of the Gibbs energy, enthalpy and entropy of formation of ThC2.
7.39
7. Electrochemical transport
3. Wagners oxidation theory
Wagners equation for the oxidation rate of a metal, forming the oxide MaOb, can in
one form be written
2
dn
RT
1
=
(
t
+
t
)
t
d
ln
p
cat
an
el
O
2
dt 8bF 2 1
x
a) Specify what the different symbols denote in this equation.
b) Examine the equality of units on both sides of the equation.
c) What are the model and requirements for the theory. Why do the different terms
enter the equation as they do?
d) This type of behaviour is called parabolic. What does it reflect? Identify the
parabolic rate constant.
e) Which deviations from the requirements of the theory may give rise to nonparabolic (e.g. linear) behaviour?
f) Assume that the oxide has formula MO and is a p-type semiconductor with doubly
charged metal vacancies as dominating point defects (DM >> DO). Assume that
pO2o >> pO2i and find how the oxidation rate varies with the outer oxygen
pressure pO2o.
g) Assume that the oxide has formula MaOb and that it is an n-type semiconductor
with Mi as dominating point defects and that DM >> DO. How does the oxidation
rate in this case vary with the outer oxygen partial pressure (assuming also now
that pO2o >> pO2i)?
h) The tracer self-diffusion coefficient for Co in CoO in air at 1200C is 9*10-9
cm2/s. Calculate the parabolic rate constant in cm2/s (which is kp in dx/dt = kp/x)
when cobalt oxidizes to CoO under the same conditions. The dominating point
defects in CoO are assumed to be singly charged metal vacancies vCo.
4. Sintering
a)
What happens during sintering? What are the driving forces? What are the
contributing mechanisms and what may be the rate limiting factors?
b)
7.40
7. Electrochemical transport
5. Creep
a)
What is creep? What are the driving forces, mechanisms, and rate limiting
factors?
b)
6. General considerations
b. The expression uphill diffusion is used for diffusion against a concentration
gradient, referring to Ficks 1st law. From what you know about
electrochemical transport, what might cause uphill diffusion?
c. Eq. 7.4 is written with conductivity as the transport coefficient. Flux of a
neutral species is not affected by the electrical potential gradient and can not
be expressed by conductivity. It is therefore meaningful to express the
equation using the random diffusion coefficient for the chemical potential
gradient and the conductivity for the electrical potential gradient. Do this
splitting/substitution. (This is mainly a simple exercise in the Nernst-Einstein
relation). Check that both parts of the resulting equation will have units that
confer with flux density.
7. Transport number and partial conductivity measurements
d. A membrane of a mixed oxygen ion and electron conducting oxide is equipped
with Pt electrodes and exposed to two different oxygen partial pressures: On
one side is dry air (pO2 = 0.21 atm) and on the other is pure oxygen (pO2 = 1.00
atm). What is the Nernst voltage of this cell at 1000C?
e. A voltage of 0.01 V is read over the cell at 1000C using a high impedance
voltmeter. What is the average ionic transport number of the oxide under these
conditions?
f. Make a sketch of the sample and gases and external circuitry during the
voltage measurement. Indicate internal and external fluxes and currents (with
directions).
g. The total conductivity of the sample is measured and the result is 0.01 S/cm.
What are the partial ionic and electronic conductivities?
h. Assuming that the sample is 0.5 mm thick, what is the flux of oxygen gas
through the sample during the voltage measurements?
i. If you had measured this flux in some way, and used this as your (only)
transport measurement, what would be the transport parameter that you find?
7.41
7. Electrochemical transport
8. Sensor
a) A sensor is made out of a membrane of yttria.stabilized zirconia (YSZ) with Pt electrodes
and is operated at 700C. As reference is used pure oxygen at 1 atm. A voltage of -950 mV is
read over the sensor (measuring on the unknown side vs the reference side). What is pO2 at
the unknown side?
b) If you assume that the unknown side is hydrogen gas at 1 atm, what does the sensor
voltage say about the water vapour content of that hydrogen gas?
9. Fuel cell
a) A fuel cell made of YSZ electrolyte with Pt electrodes is operated at 700C and runs on
hydrogen as fuel and air as oxidant, both at atmospheric pressure. Both gases contain 2%
water vapour. What is the open circuit voltage of the fuel cell?
b) Assuming that the electrolyte has a conductivity of 0.1 S/cm and is 100 m thick, calculate
the short-circuit current density, the maximum power density and the current and power
densities at a voltage of 2/3 of the open circuit voltage. (Assume in all cases that the gas
compositions remain unaffected by the current in the cell, and that electrode polarization is
insignificant).
c) Assuming the operation at 2/3 of the OCV voltage, what is the area of electrolyte needed
for a power plant of 300 MW capacity?
d) The assumption that the gas compositions are unaffected is not realistic for a real power
plant; what would it mean in terms of fuel utilization?
e) The system we have described can be operated at other temperatures. Still neglecting
electrode polarization, what are the factors and processes that enter into the temperature
dependency of the power density of the fuel cell? (Hint: refers to defects and transport.)
10. Gas separation membranes
a)
A gas separation membrane that can separate oxygen from air is constructed from
La2NiO4+d, a material where we for simplicity may assume that the defect structure is
dominated by doubly ionized oxygen interstitials compensated by electron holes. In the
following assume further that the cell is operated at 1000C and that the membrane is 100 m
thick and has no surface kinetics limitations. Assume also that the oxygen ion conductivity
(by interstitials) is 1 S/cm at 1 atm O2 and that the electronic (hole) conductivity is 200 S/cm
at 1 atm O2. The membrane is operated with atmospheric air on one side and pumped to
effectively 0.1 atm at the other side. What is the area-specific flux of oxygen?
7.42
7. Electrochemical transport
b)
What would increase the flux more: Compress the air to 10 atm or pump the vacuum side to
0.01 atm?
1. Find an expression for the chemical diffusion coefficient for oxygen ions.
2. During operation the membrane may possibly suffer from both membrane
walk-out and decomposition. Explain both phenomena qualitatively.
c)
d)
7.43
Tables
Greek letters
Constants
m/s2
alfa
Gravitational acceleration
9,80665
beta
1,6605*10-27 kg
gamma
Electron mass
me
9,110*10-31
delta
Elementary charge
kg
-19
-19
1,602*10
epsilon
Electron volt
eV
1,602*10
zeta
Proton mass
mp
1,673*10-27
kg
mn
-27
eta
Neuytron mass
1,675*10
kg
8
theta
2,99792*10
m/s
iota
Boltzmanns constant
1,381*10-23
J/K
-5
kappa
8,6174*10
-34
eV/K
lambda
Plancks constant
6,626*10
Js
my
Rydbergs constant
1,097*107
m-1
ny
Bohrs constant
2,18*10-18
23
ksi
Avogadros number
NA
6,022*10
mol-1
omikron
Gas constant
8,31451
J/molK
pi
0,0820578
Latm/molK
rho
Faradays constant
96485
C/mol
sigma
Gravitational constant
6,672*10-11
Nm2/kg2
tau
Permeability in vacuum
1,257*10-12
H/m
-12
8,854*10
F/m
ypsilon
Permittivity in vacuum
phi
9,0*109
Nm2/C2
khi
Magnetic constant
km
2*10-7
N/A2
psi
Stefan-Boltzmann constant
5,67*10-8
W/m2K4
omega
Wiens constant
2,90*10-3
mK
Vm,298K 24,4651