NORTHWESTERN UNIVERSITY
AStudyofOxidesforSolidOxideCells
ADISSERTATION
SUBMITTED TO THE GRADUATE SCHOOL
IN PARTIAL FULFILLMENT OF THE REQUIREMENTS
for the degree
DOCTOR OF PHILOSOPHY
Field of Materials Science and Engineering
By
Olivier Comets
EVANSTON, ILLINOIS
December 2013 UMI Number: 3605699
All rights reserved
INFORMATION TO ALL USERS The quality of this reproduction is dependent upon the quality of the copy submitted.
In the unlikely event that the author did not send a complete manuscript and there are missing pages, these will be noted. Also, if material had to be removed, a note will indicate the deletion.
UMI 3605699 Published by ProQuest LLC (2013). Copyright in the Dissertation held by the Author. Microform Edition © ProQuest LLC. All rights reserved. This work is protected against unauthorized copying under Title 17, United States Code
ProQuest LLC. 789 East Eisenhower Parkway P.O. Box 1346 Ann Arbor, MI 48106 - 1346 2
c Copyright by Olivier Comets 2013 All Rights Reserved 3
ABSTRACT
AStudyofOxidesforSolidOxideCells
Olivier Comets
As the world energy consumption increases, it is a question of global health to increase energy production e ciency and to reduce CO2 emissions. In that respect, solid oxide cells are solid state devices that convert directly fuel into electricity, or vice versa. In fact, when run in fuel cell mode, such devices produce electricity with e ciency up to twice that of current natural gas power plants. However, systems equipped with them have only seen limited commercialization owing to issues of cost, durability, and performance.
In this thesis, three di↵erent aspects of solid oxide cells are studied. First, the e↵ects of stress on the properties of mixed ionic electronic conducting oxides are considered.
Such oxides can be used as electrode materials, where they are often subject to large stresses, which can, in turn, a↵ect their performance. Hence, understanding the rela- tionship between stress and properties in such materials is crucial. Non-stoichiometry in strontium substituted lanthanum cobaltite is found to increase under tension and to decrease under compression. 4
Then, degradation taking place when the cell is run in electrolysis mode is discussed.
A high current allows for a high production rate of hydrogen gas. However, this can also lead to oxygen bubble nucleating in the electrolyte and subsequent degradation of the cell. The analysis conducted here shows that such nucleation phenomenon can be avoided by keeping the overpotential at the oxygen electrode below a critical value.
Finally, the growth and coarsening of catalyst nanoparticles at the surface of an oxide is studied. Scientists have developed new oxides for anodes in which a catalyst material is dissolved and exsolves under operating conditions. As the performance of the cell is controlled by the surface area of the catalyst phase, understanding the kinetics of the growth is critical to predict the performance of the cell. An approach is developed to study the growth of one particle, in the limiting case where only bulk transport is allowed. 5
Acknowledgements
As I reflect back at my time in the Department of Materials Science and Engineering at Northwestern, I realize how much I have learned, how many great people I have met, and how many amazing experiences I have lived. Undeniably, this department and the people I have met through it have played a major role in my scientific development, personal fulfillment, and my integration in the US.
First and foremost, I would like to thank my advisor, Peter Voorhees, for his knowl- edge and guidance while confronting me with such exciting and stimulating projects.
Knowing that a graduate school experience is both of academic and human nature, he encouraged me to develop my soft skills through various projects unrelated to work.
Peter, thank you for everything!
I would like to thank Scott Barnett, who played the role of a second advisor given the overlap of my research and his expertise, for insightful discussions and thrilling collaboration work. I would also like to thank my committee members Thomas Mason,
Kenneth Poeppelmeier, and Chris Wolverton for thoughts, suggestions, and insight.
I’m very fortunate to be part of such an amazing and complementary group as the Voorhees Research Group and I would like to individually thank each one of you:
Kuo-An, Tony, Thomas, Begum, Larry, Ian, Megna, Alanna, Eddie, Anthony, John T.,
Tom, John G., Kevin, Quentin, and Ashwin. I leave the group with memories of great scientific discussions, help in dire situations, and with great friendships. I am also very 6 grateful to the Barnett group for their thoughts and insights on Solid Oxide Cells, and namely to: David B., Scott, Kyle, Gareth, Ann, Beth, and David K.
This work would not have been possible without the many challenging and en- lightening discussions with our collaborators: professors Jason Nicholas, Stuart Adler,
Katsuyo Thornton, Dr. Hui-Chia Yu, and T. J. McDonald as well as with Prof. Anil
Virkar and Prof. Junichiro Mizusaki. The English in this thesis wouldn’t have been as good without the help of John, Alex, Ahmed and Kyle.
Finally, I would like to thank my loving family and friends for all their support during the process. I am grateful to my parents for the education they provided me with and to my parents, Aude, and Antoine for their constant encouragements. I am very glad to Dave Herman, Ahmed Issa, Carlos Alvarez and Begum Gulsoy for valuable friendships, great advice and much fun Ive had during grad school. Last but not least,
IwouldliketoacknowledgemyfriendPierreGarreauforhisconstantsupport,an infallible friendship, a lot of fun during grad school and many essential conversations weve had together.
This work was financially supported by the US Department of Energy (DOE) and the National Science Foundation (NSF). 7
Contents
ABSTRACT 3
Acknowledgements 5
List of Figures 11
List of Tables 14
Chapter 1. Introduction 15
Chapter 2. Background 17
2.1. ElectricityproductionintheUS 18
2.2. Solid Oxide Cells 20
2.2.1. Fuel cell mode 21
2.2.2. Electrolysis mode 22
2.2.3. Materials 23
2.2.4. Features 25
Chapter 3. The E↵ects of Stress on the Defect and Electronic Properties of
MixedIonicElectronicConductors 27
3.1. Introduction 27
3.2. Thermodynamics 29
3.2.1. Thermodynamic description of the system 29 8
3.2.2. Equilibrium conditions 31
3.2.3. NewfreeenergyfunctionandMaxwell’sequation 39
3.2.4. Chemical potential of oxygen under stress 41
3.3. E↵ects of stress on the non-stoichiometry 43
3.4. E↵ects of stress on the vacancy formation energy 45
3.5. E↵ectsofstressonthechemicalcapacitance 48
3.6.Comparisonsandpredictions 51
3.6.1. E↵ects of a hydrostatic stress on the properties of La0.8Sr0.2CoO3 52 3.6.2. Thin Films 56
3.7. Discussion 69
3.7.1. LSC thin films 69
3.7.2. Generalization to other mixed conductors 72
3.8. Conclusion and future work 73
Chapter 4. Oxygen Bubble Formation in Solid Oxide Electrolysis Cells 76
4.1. Introduction 76
4.2. Thermodynamics of nucleation 78
4.2.1. Thermodynamic model 82
4.2.2. Internal energies 84
4.2.3. Constraints 86
4.2.4. Equilibrium conditions 92
4.3. Driving force 94
4.3.1. Value of the oxygen potential 94
4.3.2. Expression of the oxygen potential 95 9
4.3.3. Expressions of the grand potentials 99
4.3.4. Change in the grand potential 103
4.3.5. Free energy change of nucleation 107
4.4. Results and discussion 112
4.4.1. Critical radius 112
4.4.2. Homogeneousandheterogeneousnucleation 114
4.4.3. E↵ects of parameters on the nucleation polarization 117
4.4.4. Critical current 118
4.4.5. Vacancy concentration 121
4.5. Conclusion and future work 122
Chapter 5. Growth and Coarsening of Nanoparticles on the Surface of an Oxide 125
5.1. Introduction 125
5.2. Background 127
5.2.1. Coarsening in 3D 128
5.2.2. Coarsening in 2D 129
5.3. Modeling considerations 129
5.4. Mathematicalformulationofthesystem 131
5.4.1.Governingequation 131
5.4.2. Boundary conditions 131
5.4.3. Particle growth rate 132
5.4.4. Undimensionalizingtheequations 133
5.5. Approach 135
5.5.1. Green’s function 135 10
5.5.2. Green’s theorem 140
5.5.3. Solving the equations 143
5.6. Extension of the model and future work 143
5.7. Conclusion 144
Chapter 6. Conclusion 146
References 148 11
List of Figures
2.1 World energy consumption in the world as predicted by the United
StatesEnergyInformationAdministrationin2011. 18
2.2 CompositionoftheelectricityproducedintheUSbyresources. 19
2.3 Projections for added electricity generation capacity as a function of
sources through 2040. 20
2.4 Schematic of a solid oxide cell running in fuel cell mode on hydrogen
gas. 22
2.5 Schematic of a solid oxide cell running in electrolysis mode on water. 23
3.1 System under consideration for the derivation of the equilibrium
conditions: oxide and gas phase delimited by an arbitrary interface
@ .33 V
3.2 Thought experiment to understand the e↵ect of stress on the
non-stoichiometry. 44
3.3 Non-stoichiometry as a function of the trace of the stress in LSC-82. 56
3.4 Schematic of the change in non-stoichiometry in a coherent and
dislocation-free thin film due to lattice mismatch with the substrate. 61 12
3.5 Schematic of the change in non-stoichiometry in a thin film grown on
asubstrateunderthermalstress. 62
3.6 Chemical capacitance versus oxygen partial pressure at T =873K
as estimated for bulk La0.6Sr0.4CoO3 ,asreportedina1.5µm-thick LSC film on GDC and according to the model. 64
3.7 Chemical capacitance versus oxygen partial pressure at T =793K
evaluated for bulk La0.8Sr0.2CoO3 ,asreportedfora45nm-thick LSCfilmonYSZandaccordingtothemodel. 67
4.1 Sketch of a SOEC under operation. If the current is high enough,
bubblescanformintheelectrolyte. 77
4.2 Schematic of oxygen bubble formation in the dense YSZ electrolyte
of a SOEC. 79
4.3 Sketch of the system under study for the derivation of the equilibrium
conditions: perfect YSZ lattice with a spherical bubble of oxygen. 83
4.4 Sketch of the oxygen potential near the oxygen electrode of a SOEC
underanappliedcurrent. 95
4.5 Driving force for nucleation explained in the perspective of the mole
fraction of oxygen vacancy. 109
g 4.6 Plot of the grand potential of the gas bubble, !v ,thehomogeneous
part of the grand potential of the oxide, !ox,andthenegative
of the elastic energy, W ,asafunctionoftheoxygenelectrode e polarization. 111 13
4.7 Critical radius of the nucleus versus electrode polarization for the
exact case as given by Eq. (4.85) and the approximation given by
(4.86). 113
4.8 Reversible work for the formation of a critical nucleus as a function
of the oxygen electrode polarization for the homogeneous nucleation
case (within a grain) and heterogeneous case (at a grain boundary). 117
4.9 Nucleation overpotential as a function of the surface energy of the
oxide. 119
4.10 Critical overpotential as a function of temperature, T ,andoxygen
O partial pressure at the oxygen electrode, PO2 .120
5.1 Schematic of catalyst nanoparticles precipitating at the surface of the
anode. 126
5.2 Schematic of the mechanism for the formation of catalyst particles at
the surface of the oxide. 127
5.3 Schematic of the configuration for the coarsening problem. 130 14
List of Tables
3.1 ParametersforLSC-82usedinestablishingFig.3.3. 57
3.2 Parameters for GDC-91, LSC-64, LSC-55 and LSC-73 used to
compute Cchem as a function of PO2 toestablishFig.3.6. 65
3.3 Parameters for LSC-82 used to compute Cchem as a function of PO2 toestablishFig.3.7. 68
4.1 Values of the parameters for nucleation of oxygen bubbles in 8-mol
% YSZ electrolyte. 110 15
CHAPTER 1
Introduction
The growing needs in energy and the depletion of the oil resources have made man consider new, cleaner and sustainable ways to produce energy. In that respect, asolidoxidecell(SOC)isasolidstatedevicethatconvertsdirectlychemicalenergy into electricity, or vice versa.Runinfuelcellmode,aSOCproduceselectricitybya direct oxidation of the fuel, skipping the conversion steps into thermal and mechanical energy present in the standard fossil fuel power plants. Production of electricity by this process is up to twice as e cient as in standard power plants. Such cells can also be run in electrolysis mode, to regenerate gas. SOCs are one of several di↵erent fuel and electrolysis cells, but are of interest because of higher e ciency, low emissions, fuel flexibility and potential long-term stability. Thanks to these advantages, solid oxide fuel cells have a wide range of applications from auxiliary power units in big-rig vehicles to dispersed stationary power generation. However, for such systems to be mass produced the issues of performance, durability, and cost must be addressed.
The goal of this thesis is to use thermodynamics and other mathematical tools to study di↵erent aspects of solid oxide cells. In doing so, we hope to gain better under- standing of the processes taking place in the cell and of the cell intrinsic limitations.
Such basic understanding is the cornerstone of SOC systems commercialization.
Chapter 2 o↵ers some background for the current thesis work. First, SOC Research is motivated by the increase in global energy demand, the need for e cient processes 16
to produce electricity and the increase in CO2 emissions. Then, an overview of the fundamentals of operations and the di↵erent components of a cell are discussed.
Chapter 3 discusses the e↵ects of stress in mixed ionic electronic conducting oxides.
Such oxides are used in a variety of di↵erent applications (e.g. sensors, SOCs) and often are in a state of stress (e.g. thermal, mismatch). As stress can a↵ect their performance, understanding how such oxide behave under stress will allow to better predict their performance in applications.
Chapter 4 presents the degradation of solid oxide electrolysis cells (SOECs) by oxygen bubble formation in their electrolyte. SOECs are used to regenerate gas using electricity. High rates of production are achieved with high currents. However, when the voltage applied to the cell is above a critical value, bubbles start forming in the electrolyte of the cells, leading sometimes to deleterious consequences.
Chapter 5 tackles the growth and coarsening of catalyst nanoparticles at the surface of an oxide. Scientists have developed novel anode materials, where catalyst is dissolved within the oxide and exsolves under operating conditions. Because the performance of the cell is controlled by the surface area of those catalyst particles, understanding the kinetics of the process is crucial to predict the performance of the cell.
Finally, Chapter 6 summarizes the main results of this thesis. For each of the projects, only the most important elements of the future work are recalled. 17
CHAPTER 2
Background
It is no surprise that the world energy demand is growing quickly. In 2008, the
United States Energy Information Administration (EIA) evaluated the world energy demand to grow by 2% yearly [1]. In other words, every 10 years the world adds capacities equivalent to the entire annual energy production of the US. The projections for the energy demand, reported in Fig. 2.1, show that growth is mainly driven by non-OECD countries. This fact can be understood when considering that non-OECD countries represented 80% of the world population but 50% of the energy consumed worldwide in 2008. The development of those countries is synonymous with a dramatic increase in energy demand. This then results in an unprecedented increase in the production of greenhouse gases (e.g. CO2)andotherbyproducts(e.g.heat).Developing e cient and cleaner ways to produce energy is thus crucial.
After considering the energy landscape in the US, it will be shown that electricity production represents a large share of the energy mix and is a highly ine cient process.
Solid oxide cells are devices that address that problem, capable of e ciently convert- ing chemical energy to electricity and vice versa. Various aspects of those devices: electrochemistry, materials and features will be presented in the second section. 18
Figure 2.1. World energy consumption in the world as predicted by the United States Energy Information Administration in 2011 [1]. The grow- ing demand is mainly driven by the non-OECD countries, who represent more than 80% of the global population.
2.1. Electricity production in the US
Today, electricity generation in the United States accounts for approximately 40% of the energy consumed in the US [2]. Fig. 2.2 represents how the electricity is split among the di↵erent resources. Fossil fuels (coal, gas, oil), which have by far the biggest carbon footprint, represent 71% of all the electricity produced in the US. The e ciency for a fossil fuel-based power plant is currently about 33% [3], i.e. two third of the energy used to produce electricity is wasted. Furthermore, electricity generation is responsible for 40% of the total US carbon dioxide production [4], almost exclusively due to fossil fuels. All this proves that electricity generation in the US remains a highly ine cient process and responsible for much greenhouse gas emissions. 19
Figure 2.2. Composition of the electricity produced in the US by re- sources [2]. ”Renewable” includes hydro, geothermal, solar, wind and tide. Fossil fuels (coal, gas and oil) contribute to 71% of the electricity produced in the US.
Petit et al. have shown a positive correlation between atmospheric levels of CO2 and the earth temperature [5]. Given the recent rise in CO2 concentration in the atmosphere, the Nobel prize-winning Intergovernmental Panel on Climate Change has predicted a temperature increase of 2 to 6 Cbytheendofthecentury[6]. Such an increase in temperature can in turn lead to a disruption of the earth’s fragile climate.
However, addressing the issue of electricity generation could result in vital progress in reducing CO2 emissions, controlling the atmosphere temperature and protecting the planet’s fragile equilibrium.
With a steady increase in the demand for electricity, the reduction of these emissions must result from the development of a more e cient electricity generation process. The
first step is to use cleaner fuels, e.g. natural gas rather than coal. For that matter, the
US EIA projects a drastic increase in the number of natural gas plants as additional capacities, along with renewable resources, which have little carbon footprint, as shown 20
Figure 2.3. Projections for added electricity generation capacity as a function of sources through 2040. Source: United States EIA [7].
in Fig. 2.3. Solid oxide cells can be operated on natural gas and have nearly twice the e ciency of current plants running on such fuels. However, for such systems to be mass produced and integrated into power generators, research is needed to lower the costs and limit degradation.
2.2. Solid Oxide Cells
ASolidOxideCell(SOC)isasolidstateelectrochemicaldevicecapableofe ciently converting the chemical energy of a fuel gas to electricity and vice versa. The cell can be run both in the fuel cell mode to produce electricity, and in the electrolysis mode to regenerate the gas. 21
2.2.1. Fuel cell mode
In conventional power plants, gas and oxygen are mixed together and combusted. The heat produced by this reaction is transferred to a fluid which then drives a turbine, activating an alternator to finally generate electricity. Given the number of di↵erent steps in the process and that its e ciency is limited by the Carnot cycle, the maximum theoretical e ciency of a traditional power plant is 47%. In a Solid Oxide Fuel Cell
(SOFC), no direct combustion take place. The reactants (fuel and air) are, in fact, spatially separated and involved in electrochemical reactions at electrodes, separated by an electrolyte, much like in a battery. However, unlike a battery, a fuel cell does not need to be recharged and will run as long as the reactants are supplied. Fig. 2.4 is a schematic of a SOC running in fuel cell mode on hydrogen gas. Oxygen is reduced on
2 the cathode to form oxygen ions (O ). Those ions are then transferred to the anode via the electrolyte, where they react with hydrogen gas to form water and regenerate electrons. Electrons are thus produced at the anode and consumed at the cathode generating a current. The reactions taking place are:
1 2 at the cathode: O (g) + 2e0 O (ox)(2.1) 2 2 !
2 (2.2) at the anode: H +O (ox) H O(g) + 2e0 2 ! 2 1 and overall: O (g) + H H O(g)(2.3) 2 2 2 ! 2
The overall reaction is a reaction of combustion. Finally, in a SOFC, the reaction of combustion has essentially been split up into it’s reduction and oxidation reactions in order to directly use the flow of electrons. 22
Figure 2.4. Schematic of a solid oxide cell running in fuel cell mode on hydrogen gas. Air and fuel are fed to the cell. At the cathode, air is reduced to oxygen ions. Those ions are then transported to the anode through the electrolyte, where they react with the fuel gas (here H2), forming water, and regenerating electrons. The flow of electrons from the anode to the cathode is then used outside of the cell to power appliances.
2.2.2. Electrolysis mode
The operating principle of a Solid Oxide Electrolysis Cell (SOEC) is the very opposite of that of a SOFC. Fig. 2.5 is a schematic of a solid oxide cell running on water in electrolysis mode. Water and an electric current are fed to the cell. Water molecules 23
Figure 2.5. Schematic of a solid oxide cell running in electrolysis mode. Water and electric power are fed to the cell. At the cathode, water is reduced to hydrogen gas and oxygen ions. Those ions are then trans- ported to the cathode through the electrolyte, where they recombine, regenerating oxygen gas and electrons. react with electrons at the cathode producing oxygen ions and hydrogen gas. Oxy- gen ions are then transported via the electrolyte to the anode where they recombine, regenerating oxygen gas and electrons.
2.2.3. Materials
Because the electrolyte and electrodes serve di↵erent purposes, di↵erent materials and geometries are used. 24
The role of the electrolyte is threefold: to physically separate the fuel and the oxidant, to transport the oxygen ions from one electrode to the other, while preventing the passing of electrons. Thus an electrolyte must:
(1) be fully dense,
(2) exhibit high ionic conductivity,
(3) have low electronic conductivity,
(4) be stable in oxidizing and reducing environment,
(5) be chemically compatible with electrode materials,
(6) have a thermal expansion relatively similar to that of the electrodes.
Typical electrolyte materials are 8 mol% yttria-stabilized zirconium (YSZ), Y2O3-ZrO2, and gadolinium doped ceria, Gd2O3-CeO2 [8,9]. As electrolyte resistance is a function of thickness, electrolytes are made very thin, typically < 10 µm.
The anode of a SOFC provides reaction sites for the oxidation of the fuel. It also supports the transport of the various species to and from those reaction sites: gases, oxygen ions and electrons. As a result, it should be porous to allow for the gases to di↵use, capable to conduct oxygen ions and electrons. Finally, anodes must have the following features:
(1) high porosity
(2) electronic conductivity
(3) ionic conductivity
(4) stability at high temperatures and in reducing environments
(5) mechanical compatibility with electrolyte
(6) chemical compatibility with electrolyte and interconnect 25
(7) catalyst activity toward the oxidation of the fuel
(8) fuel flexibility and resistance to impurities and carbon deposition
Typical anode comprise a mixture of nickel and YSZ. The metal phase (nickel) provides a path for the electrons, while the oxide phase (YSZ) provides the transport of oxygen ions and the pores allows for gas di↵usion. The pores are created from the reduction of nickel oxide to nickel after exposure to the fuel. In this case, the active sites are at the junction of the three phases, known as the triple phase boundaries. Often the anode is fabricated using a dual-layer geometry, where a first layer of thickness .5 1mm provides the support for the cell and a second layer of thickness 10 50 µmisthe functional layer.
The cathode of a SOFC is exposed to air. Like the anode, the cathode needs to be porous, capable of transporting both oxygen ions and electrons, compatible (me- chanically and chemically) with the electrolyte and interconnect and stable at high temperature. However, the anode is required to be stable under an oxidizing environ- ment while having a catalytic activity on the reduction of molecular oxygen. Typical cathodes are made of a mixture of strontium substituted lanthanum manganite and
YSZ.
2.2.4. Features
Various aspects of SOCs have caused them to emerge recently as a serious solution to address the problem of growing electricity demand. Solid oxide cells operate at very high temperature (between 400 and 1000 C), which make expensive catalysts unnec- essary. Furthermore, the high quality of the by-products can be used for cogeneration 26
(in fuel cell mode), boosting the e ciency of the process even more. Because the cell is entirely solid state, there are no moving parts, making it silent and easier to run. The electrolyte does not require any management, unlike sulfuric acid fuel cells for exam- ple. As a result, SOCs present a potential long life expectancy of 40,000-80,000 hours of operation. SOFCs can achieve e ciencies of 45 to 60%, and up to 90% with heat recovery [9]. Last, a tubular geometry has recently emerged promising much shorter start up times than the classic planar geometry, typically on the order of minutes.
Although this technology has been known for over 160 years, much more research is necessary to optimize the cell, control its long term degradation, and eventually bring the cost down, making it competitive with conventional less e cient technologies. 27
CHAPTER 3
The E↵ects of Stress on the Defect and Electronic Properties
of Mixed Ionic Electronic Conductors
3.1. Introduction
Oxides transporting both ionic and electronic species are called mixed ionic elec- tronic conductors (MIECs). Such materials are of particular interest in applications where simultaneous ion and electron conduction is required, such as at the electrodes of
Solid Oxide Cells [8]. Recent studies have shown that such oxides under stress (e.g. in thin film form) display very di↵erent properties —e.g. non-stoichiometry and kinetics— from that of bulk materials [10,11]. This di↵erence in behavior between the thin film configuration and bulk is often attributed to stresses developing in the film. Such stress can be thermal, originating from di↵erent operating and firing temperatures, or due to a misfit between the lattice parameters of the MIEC and the substrate. As oxides in most applications are in the form of thin films, stress is present and it is critical to understand its e↵ects to predict their behavior under operating conditions. Finally, because oxide lattice parameter, oxygen non-stoichiometry and electrical properties are closely related [12], stress will influence all of these simultaneously.
Two types of MIECs have been reported and classified according to their electronic conduction mechanism. The first has a metallic-like electronic conduction mechanism, 28 mediated by holes present in a partially filled delocalized band. Such behavior is de- scribed by the itinerant electron model developed by Mizusaki et al. [13]andLankhorst et al. [14–17]. La1 xSrxCoO3 (LSC) is a such example [18]. The second type has an electronic conductivity described by an activated electron hopping mechanism, also known as the small polaron model [19]. A small polaron is a defect created when an electronic carrier is trapped at a given site fostering a local distortion of the lattice.
The hopping of such defects (the carrier plus its polarization field) is responsible for the electronic conductivity in such materials as La1xSrxMnO3 (LSM) for x ¡ .2 [18]. Un- fortunately, the defect structure is also much more complex in such material, evidenced by extensive work on the topic, e.g. [20–29].
The purpose of this chapter is to illustrate the e↵ects of stress on the properties of mixed conductors with a perovskite structure in equilibrium with an atmosphere, through the example of LSC. In the first section, results from previous studies on elas- tically stressed crystals [30–34]areusedtoderivetheequilibriumconditionsandthe expression of the chemical potential of oxygen in the oxide under stress. In the follow- ing three sections, the expressions for the non-stoichiometry, the chemical capacitance and the vacancy formation energy are respectively derived for an oxide under stress.
The fifth section is dedicated to predictions of the model. First, changes in the non- stoichiometry of an oxide under a hydrostatic stress are considered. The e↵ects on the hole concentration and on the ionic and electronic conductivities are also discussed.
Second, predictions of the chemical capacitance as a function of oxygen pressure are compared to experimental data for a thin film under mismatch strain and thermal build up [10,11]. Third, calculations of the change in vacancy formation energy due to 29 stress are compared to ab initio results [35]. The qualitative agreement resulting from those comparisons show that stress is not the sole controlling factor for the di↵erence in behavior between thin film and bulk. Possible sources of the quantitative discrepancy are discussed in the last section. Finally, the treatment developed here is shown to extend to other mixed conducting oxides, with more complex defect equilibria.
3.2. Thermodynamics
The approach employed here was initially introduced by Cahn and Larch´e[30], and developed by Johnson and Schmalzried [31, 32], see [33] for a review. Swaminathan et al. used a similar approach to study the di↵usion of charged defects in ionic solids
[34, 36]. First, the equilibrium conditions are derived for an oxide under stress, in equilibrium with oxygen gas. Introducing a new free energy function allows us to derive a Maxwell’s equation. Integrating this equation between a stress free state and astateunderstressfinallyyieldstheexpressionofthechemicalpotentialofoxygenin the oxide.
3.2.1. Thermodynamic description of the system
Under consideration is a dislocation-free slab of oxide perovskite structure of general chemical formula ABO3 in equilibrium with a gas containing molecular oxygen, O2.
The oxide has three distinct sublattices: two for the cations (A and B), one for the oxygen ions (O) and the interstitial sites are all vacant. The species assumed to be present in the gas are molecular oxygen, O2,andothergasesthatdonotreactwith the oxide (e.g. N2). Elements from sublattices A and B are not soluble in the gas and 30 no other phases are assumed to form. The various thermodynamic densities relative to the crystal are referred to a reference or stress-free state, while those relative to the gas are referred to the actual state. Thermodynamic densities expressed per-unit-volume in the reference state are designated with a superscript 0.
The oxide used to derive the results in the rest of the chapter is La1 xSrxCoO3 , where x is the strontium substitution level and is the non-stoichiometry. Considering absolute charges, oxygen with a 2oxidationstateandneutralvacanciesarefoundon the O sublattice. The A sublattice is populated with lanthanum ions (III), strontium ions (II), and neutral vacancies, while the B sublattice is populated with cobalt ions
(III) and neutral vacancies. Considering relative charges and using the Kr¨oger-Vink
x x x notation, these elements are noted OO,VO·· ,LaA,SrA0 ,VA000,CoB,VB000. Any given ion must occupy a site on one of the subblattices. Dislocations, exchange of atoms between the anionic and cationic sublattices and interstitial atoms are not included in the model.
The internal energy density per unit volume of the oxide in the stress-free state, eox,istakentobeafunctionoftheentropysox,thedeformationgradienttensor ,the v0 v0 F electric displacement field D and the number densities of the di↵erent elements ⇢0 x , LaA
⇢0 , ⇢0 , ⇢0 x , ⇢0 , ⇢0 x , ⇢0 and ⇢0 : SrA0 VA000 CoB VB000 OO VO·· h·
ox ox (3.1) e s , , D,⇢0 x ,⇢0 ,⇢0 ,⇢0 x ,⇢0 ,⇢0 x ,⇢0 ,⇢0 v0 v0 F LaA SrA0 VA000 CoB VB000 OO VO·· h· ⇣ ⌘ Avariationofanyofthesevariablesinducesachangeintheinternalenergy:
ox ox ox x x ev =T sv + T : F + JE D + µLa ⇢La0 + µSr0 ⇢Sr0 + µV000 ⇢V0 0 0 · A A A A0 A A000
(3.2) + µCox ⇢0 x + µV ⇢0 + µOx ⇢0 x + µV ⇢0 + µh ⇢0 B CoB B000 VB000 O OO O·· VO·· · h· 31
where T ox is the absolute temperature of the oxide, T is the first Piola-Kirchho↵stress tensor, J =detF is the Jacobian of the transformation (also equal to the ratio of the volume of a cell in its deformed state to that in its non-deformed state J = dv/dv0),
@eox v x E is the electric field and µi = 0 is the chemical potential of specie i (i =LaA,SrA0 , @⇢i0 x x V000,Co ,V000,O ,V·· and h·). The symbol ” ”representstheclassicalscalarproduct A B B O O · while ”:” represents the tensorial scalar product.
g The internal energy density of the gas phase in the actual state ev is a function of
g b the entropy sv,thepressureinthebubbleP and the number density species: ⇢O2 and
⇢N2 —assuming nitrogen is the only other nonreactive gas. The internal energy of the gas phase is of the form:
g g b (3.3) ev sv,P ,⇢O2 ,⇢N2 Avariationofanyofthesevariablesinducesachangeintheinternalenergyof:
g g g b gas gas (3.4) e = T s P + µ ⇢O + µ ⇢N v v O2 2 N2 2
g gas gas where T is the temperature of the gas phase and µO2 and µN2 are the chemical potential of oxygen and nitrogen in the gas phase.
3.2.2. Equilibrium conditions
The thermodynamic equilibrium conditions are obtained using a Gibbsian variational approach, stating that the energy of an isolated system is at a minimum. An arbi- trary volume of the system containing both phases is first identified, as depicted V in Fig. 3.1. This volume is then isolated from the rest of the system and subject to 32 virtual perturbations. In order to do so, global constraints must be taken into account.
The condition of no heat flow translates into constant entropy, the absence of atomic
flux across the interface @ translates into constant number of atoms (O, La, Sr, Co V and N) and constant charge in the system [31]. Other constraints that need to be included are local constraints: electrostatics, lattice site conservation and mechanical considerations.
The total energy of the thermodynamic system defined by = + is the sum V Vox Vg of the internal energies of both phases
ox g (3.5) " = e dv0 + e dv +[surfaceterms] v0 v g ZZZVox0 ZZZV where ”[surface terms]” group all the integrals on the surfaces: @ and ⌃. Note that V for this problem, only the bulk equilibrium conditions are important, thus the surface integrals will not be explicitly treated.
As mentioned above, global thermodynamic constraints on the system must first be taken into account
(1) Constant entropy:
ox g (3.6) = s dv0 + s dv v0 v S g ZZZVox0 ZZZV
(2) Constant charge:
(3.7) = 3⇢La0 x +2⇢Sr0 +3⇢Co0 x + ⇢h0 2⇢O0 x dv Q A A0 B · O ZZZVox0 ⇣ ⌘ 33
(a) (b)
Figure 3.1. System under consideration for the derivation of the equilib- rium conditions: oxide and gas phase delimited by an arbitrary interface @ .Thethermodynamicdensitiesrelativetotheoxidearereferredtoa stressV free state, while those relative to the gas are referred to the actual state.
(3) Constant number of atoms:
(3.8) O = ⇢0 x dv0 +2 ⇢O dv OO 2 N g ZZZVox0 ZZZV
(3.9) La = ⇢La0 x dv0 N A ZZZVox0
(3.10) Sr = ⇢Sr0 dv0 N A0 ZZZVox0
(3.11) Co = ⇢Co0 x dv0 N B ZZZVox0
(3.12) N =2 ⇢N2 dv N g ZZZV
Nitrogen atoms remain in the gas and lanthanum, strontium and cobalt are
not soluble in the gas. 34
Those constraints are accounted for in the Lagrangian of the system:
(3.13) "⇤ = " T cS ONO LaNLa SrNSr CoNCo NNN
where " is the total energy of the system, defined by Eq. (4.10), Tc, o,andthe is—for i =O, La, Sr, Co and N— are the Lagrange multipliers associated with the constraints aforementioned. The first variation of this energy is given by
(3.14) "⇤ = " T c S o Q O NO La NLa Sr NSr Co NCo N NN
Substituting the expression of the internal energies, Eq. (4.5) and (4.9), and using the global constraints, (3.6)-(3.12), in that equation yields
ox ox "⇤ = [T Tc] sv + T : F + JE D +[µLax La 3eo c] ⇢La0 x 0 · A A ZZZVox0 ⇢
+[µSr0 Sr 2eo c] ⇢Sr0 + µV000 ⇢V0 A A0 A A000
x x +[µCo Co 3eo c] ⇢Co0 + µV000 ⇢V0 B B B B000
x x +[µO O +2eo c] ⇢O0 + µV·· ⇢V0 +[µh eo c] ⇢h0 dv0 O O O O·· · · g g b + [T Tc] sv P g ZZZV ⇢ gas gas +[µ 2 O] ⇢O +[µ 2 N] ⇢N dv O2 2 N2 2 +[surfaceterms](3.15)
All the variations appearing in Eq. (3.15) are not independent, they are linked via local constraints. 35
First, every site of each sublattice must be occupied either by an atom or a vacancy, ie.
A (3.16) ⇢0 x + ⇢0 + ⇢0 = ⇢ LaA SrA0 VA000
B (3.17) ⇢0 x + ⇢0 = ⇢ CoB VB000
O (3.18) ⇢0 x + ⇢0 ,⇢0 = ⇢ OO VO·· h· where ⇢A, ⇢B and ⇢O are the number densities of sites on each of the sublattices. The perovskite structure further requires ⇢A = ⇢B = ⇢O/3. This imposes a relationship between the concentrations of the di↵erent elements.
Furthermore, the electric displacement must satisfy Gauss law in the oxide:
(3.19) D =3⇢La0 x +2⇢Sr0 +3⇢Co0 x 2⇢O0 x r· A A0 B O
Noting the electric potential, we can rewrite:
(3.20) E D = D = [ ( D) ( D)] = ( D)+ ( D) · r · r· r· r · r·
Using this decomposition, the integral involving the electric displacement in the ex- pression of " simplifies to [31]:
E Ddv0 = eo 3⇢La0 x +2⇢Sr0 +3⇢Co0 x 2⇢O0 x dv · A A0 B O ZZZVox0 ZZZVox0 n o (3.21) + D nda @ ox+⌃ · Z V 0 36
The last transformation involves the elastic term T : F = Tji Fij using the Einstein notation. Using the divergence theorem, the integral on the elastic strain energy can be rewritten as
ox0 (3.22) Tji Fijdv = Tjinj uida Tji,j uidv ⌃ ZZZVox0 ZZ 0 ZZZVox0 where the index after the comma in Tji,j denotes a derivative with respect to the i-th component, nox0 is the normal to the interface pointing outward and u is the displacement vector.
Finally, using those local constraints, Eq. (3.16), (3.17), (3.18), (3.21) and (3.22) in (3.15), the first variation of internal energy of the system under the constraints is rewritten as
ox ox "⇤ = [T Tc] sv +(T ) u +[⌘Lax La 3eo c] ⇢La0 x 0 ·r A A ZZZVox0 ⇢
x x +[⌘Sr0 Sr 2eo c] ⇢Sr0 +[⌘Co Co 3eo c] ⇢Co0 A A0 B B
+[⌘Ox O +2eo c] ⇢O0 x +[⌘h eo c] ⇢h0 dv0 O O · · g g b + [T Tc] sv P g ZZZV ⇢ gas gas + µ 2 O ⇢O + µ 2 N ⇢N dv O2 2 N2 2 ⇥ ⇤ ⇥ ⇤ +[surfaceterms](3.23) 37 where more surface integrals have been added to the last term and the electrochemical potentials are defined as
(3.24) ⌘Lax = µLax µV +3eo A A A000
(3.25) ⌘Sr = µSr µV 2eo A0 A0 A000
(3.26) ⌘Cox = µCox µV 3eo B B B000
(3.27) ⌘Ox = µOx µV +2eo O O O··
(3.28) ⌘h· = µh· + eo
As all the variations in Eq. (3.23) are now independent, the bulk equilibrium con- ditions are read by setting the terms in brackets to 0:
the thermal equilibrium conditions imposes a uniform and constant tempera- • ture throughout the system:
ox g (3.29) T = T = Tc
the mechanical equilibrium condition imposes that •
(3.30) T =0 ·r
the chemical equilibrium condition states that the chemical potential of each • species is constant in the system, and specifically
x (3.31) ⌘LaA = La +3eo c 38
(3.32) ⌘ = +2e SrA0 Sr o c
x (3.33) ⌘CoB = Co +3eo c
(3.34) ⌘Ox = O 2eo c O
(3.35) ⌘h· = eo c
gas (3.36) µO2 =2 O
gas (3.37) µN2 =2 N
Because the crystal under study is a mixed ionic electronic conductor, we
can make the assumption that it is locally charge neutral, that is:
(3.38) 3⇢0 x +2⇢0 +3⇢0 x =2⇢0 x LaA SrA0 CoB OO
Using this new conditions in Eq. (3.23) simplifies the chemical equilibrium
condition to:
1 gas (3.39) µO = µ 2 O2
where µO is the chemical potential of oxygen in the oxide
(3.40) µO = µOx µV +2⌘h O O·· ·
Although such assumptions make the electrostatic term disappear from the
expressions, this is not in contradiction with the development of an electric 39
field [16]. This equilibrium is consistent with the reaction [37]
x . 1 (3.41) O +2h V·· + O (gas) O ! O 2 2
3.2.3. New free energy function and Maxwell’s equation
Applying a stress to an oxide changes its energy, which in turn a↵ects its chemical potential. This dependence can be determined by using a Maxwell equation for a free energy function [30–33]. This Maxwell equation is integrated from the initial (stress- free) state to the final (stressed) state yielding the oxygen chemical potential under stress.
We introduce a new free energy function for the oxide
ox ox ox (3.42) g = e Ts ij✏ij v0 v0 v0
where T is the temperature of the system, ✏ij and ij are the Eulerian strain and stress tensors that follow from standard linear elasticity, i, j =1, 2, 3andimplicitsummation over repeated indices from 1 to 3 is assumed. Thus, ✏ij ij represents the scalar product of those two tensors.
Using the same description as above, the change in the internal energy of the oxide in the limit of small strain is [31,33]:
ox ox e =T s + ✏ + ⌘ x ⇢ x + ⌘ ⇢ + ⌘ ⇢ v0 v0 ij ij OO OO VO·· VO·· h· h·
(3.43) + ⌘ x ⇢ x + ⌘ ⇢ + ⌘ ⇢ + ⌘ x ⇢ x + ⌘ ⇢ LaA LaA SrA0 SrA0 VA000 VA000 CoB CoB VB000 VB000 40
where ⌘i = µi + zi is the electrochemical potential of specie i, zi its charge (e.g. -2
x for i =OO)and is the electric potential. Note that the e↵ect of the electric energy has been factored into the electrochemical potentials. As mentioned in the previous paragraph, the variations present in this equation are not all independent, but are coupled via the conservation of sublattice sites, Eq. (3.16), (3.17) and (3.18), and the local charge neutrality, (3.38). Using those conditions in Eq. (3.43), the total derivative of the internal energy of the oxides simplifies to
ox ox (3.44) de = Tds + d✏ + µ d⇢ x + µ d⇢ x + µ d⇢ + µ d⇢ x v0 v0 ij ij O OO La LaA Sr SrA0 Co CoB
where µO is defined by Eq. (3.27). The variations in (3.44) are now all independent.
Using Eq. (3.44), one can evaluate the total derivative of the free energy gox,defined v0 by (3.42),
ox ox (3.45) dg = s dT ✏ijd ij + µOd⇢Ox + µLad⇢Lax + µSrd⇢Sr + µCod⇢Cox v0 v0 O A A0 B
Noting that the number density can be linked to the oxygen non-stoichiometry, ,in
O O La1 xSrxCoO3 : ⇢Ox = ⇢ ⇢V = ⇢ (1 /3), the total derivative of the new free O O·· energy rewritten as
ox ox 1 O (3.46) dg = s dT ✏ijd ij + ⇢ µOd + µLad⇢Lax + µSrd⇢Sr + µCod⇢Cox v0 v0 3 A A0 B 41
Finally taking the cross derivatives of the second and third terms yields a Maxwell relation involving the oxygen chemical potential:
@µ 3 @✏ (3.47) O = ij @ ⇢O @ ij T, kl=ij ,⇢i T, kl,⇢i ✓ ◆ 6 ✓ ◆
Knowing the constitutive equation for the strain, this equation can be integrated to provide the chemical potential of oxygen as a function of stress.
3.2.4. Chemical potential of oxygen under stress
Strain can result from numerous sources. Here we consider three such sources. One is a change in lattice parameter with temperature, in materials that have a nonzero thermal expansion coe cient. Similarly, a change in the non-stoichiometry can induce stress. There can also be mismatch strain that is a result of placing a thin film of one lattice parameter coherently (continuous lattice planes) on a substrate with another lattice parameter. Accounting for all of these sources of strain, the relationship between strain and stress is,
c a T (3.48) ✏ij = Sijkl kl + ✏ ( ) ij + ✏ ij + ✏ (T ) ij where S is the compliance tensor, ✏c( )=ec[ o]/3, ✏a and ✏T (T )arethemagnitude ijkl of the isotropic compositional, mismatch and thermal strains respectively. ec is the compositional strain coe cient and is defined by Chen et al. [12]as
ˆ c 1 1 @ ln V (3.49) e = C = 3 3 @cV O·· !T,P 42 where Vˆ is the specific volume of the oxide and o is the non-stoichiometry at which the mismatch strain is computed. Using Eq. (3.48) in (3.47) yields
@µ ec O = @ ⇢O ij ij T, kl=ij ✓ ◆ 6 assuming the compliance tensor and compositional strain coe cient are independent of non-stoichiometry. Integrating this equation between the oxide under a state of
o hydrostatic pressure P and stress, ij,yields
ec (3.50) µ ( , )=µ ( ,P o)+ ( + P o ) O ij O ⇢O kk kk
o where summation over repeated indices from 1 to 3 is assumed, µO ( ,P )isthe chemical potential of oxygen at non-stoichiometry and under hydrostatic pressure of P o =1atm,whichismeasuredexperimentally. P o =1atmisassumedto | kk| be the case in the rest of the chapter, Eq. (3.50) simplifies to
ec (3.51) µ ( , )=µ ( , 0) + ( ) O ij O ⇢O kk
where again µO ( , 0) designates the bulk chemical potential of oxygen at under no stress. Thus, the chemical potential of oxygen changes linearly with the trace of the stress, with a direction that depends on the change in the volume of the oxide with vacancy concentration. Since the compositional coe cient in mixed conducting oxides can be as large as ec 0.10 [12], stress may induce large changes in composition. Note ⇡ that the stress, kk( ), is a function of the non-stoichiometry. Solving Eq. (3.48) for 43 the stress,
(3.52) ( )=C ✏ ✏c ( ) ✏a ✏T (T ) kl klmn mn mn mn mn highlights that dependence on non-stoichiometry. Cklmn is the sti↵ness tensor, the inverse of the compliance tensor Sijkl.
In order to determine the non-stoichiometry under stress, the dependence of µO on