Chenhuan
Zhao†
a,
Yifeng
Li†
a,
Wenqiang
Zhang
a,
Yun
Zheng
ab,
Xiaoming
Lou
a,
Bo
Yu
*a,
Jing
Chen
*a,
Yan
Chen
*c,
Meilin
Liu
*d and
Jianchen
Wang
a
aInstitute of Nuclear and New Energy Technology (INET), Collaborative Innovation Center of Advanced Nuclear Energy Technology, Tsinghua University, 30 Shuang’qing Road, Beijing 100084, P. R. China. E-mail: cassy_yu@mail.tsinghua.edu.cn; jingxia@mail.tsinghua.edu.cn
bDepartment of Chemical Engineering, University of Waterloo, Waterloo, ON, Canada
cGuangzhou Key Laboratory for Surface Chemistry of Energy Materials, New Energy Institute, School of Environment and Energy, South China University of Technology, Guangzhou 510006, China. E-mail: escheny@scut.edu.cn
dSchool of Materials Science and Engineering, Georgia Institute of Technology, Atlanta, GA 30332-0245, USA. E-mail: meilin.liu@mse.gatech.edu; Tel: +1-404-894-6114
First published on 13th September 2019
Solid oxide cells (SOCs) have the potential to be the most efficient energy storage and conversion systems. To minimize energy loss due to charge and mass transport associated with the operation of SOC systems at intermediate temperatures, electrodes and electrolytes containing different types of heterointerfaces have been designed, fabricated, and tested under various conditions. While heterointerfaces can significantly enhance not only the ionic and/or electronic conductivity but also the electrocatalytic activity and stability of SOC components, as predicted by theoretical calculations and demonstrated by experimental results, the mechanisms of these enhancements are yet to be fully understood. In this review, we start with an overview of the techniques for fabrication of heterointerfaces with controlled composition, structure, and morphology. Then, the latest developments in performance enhancement of SOCs with heterointerfaces are summarized, including boosting the ionic conductivity of heterostructured electrolytes (oxygen ion conductors and proton conductors) and increasing the electrocatalytic activity and durability of heterostructured electrodes (oxygen electrodes and fuel electrodes). Subsequently, we will highlight the unique attributes of heterointerfaces in the enhancement of the SOC performance and provide important insights into the mechanisms of performance enhancement in order to establish the scientific basis for rational design of better electrolyte and electrode materials. Finally, the remaining challenges in design and fabrication of novel materials for advanced solid-state electrochemical systems will be discussed, together with possible strategies to overcome these critical issues, new research directions, and future perspectives.
Broader contextLow-cost and efficient energy storage technologies are the enabler for the implementation of intermittent renewable energy technologies. Solid oxide cells (SOCs), when operated in reversible mode, have the potential to be the most efficient option for large-scale energy storage and power generation. However, broader commercialization of SOC technology hinges on breakthroughs in material development for high performance and long operational life. Heterostructured materials have demonstrated significantly enhanced charge and mass transport and electro-catalytic activity, especially at relatively low temperatures. This review provides a timely and comprehensive update on the latest advancements in the development of heterostructured SOC components. The methodology and engineering strategies are valuable not only to SOC technology but also to a wide range of energy applications for a clean and sustainable energy future, including photo/electrochemical water splitting and solar thermal CO2 reduction. |
A SOC consists of a dense electrolyte, an oxygen electrode, and a fuel electrode. During operation, the oxygen electrode provides active sites for the oxygen reduction or evolution reaction (ORR/OER)12–14 whereas the fuel electrode provides active sites for oxidation of fuel (e.g., H2, CH4, and CO) or production of fuel (e.g., electrolysis of H2O and CO2). Both electrodes must provide not only a sufficient number of active sites for the electrode reactions but also proper pathways for rapid transport of the species involved in the electrode reactions (such as ions, electrons, and gas molecules). The electrolyte in an SOC facilitates ion migration (oxygen ions, protons, or both) between the oxygen and the fuel electrodes while blocking the transport of electrons and gas molecules to avoid the short circuiting of the cell. SOCs can operate as a solid oxide fuel cell (SOFC) or as a solid oxide electrolysis cell (SOEC) as needed.15–18 SOFCs can convert chemical fuels directly to electricity with high energy efficiency, high energy and power densities, and fuel flexibility.3,4,7,19,20 SOECs can store “excess” electric energy in the form of chemical fuel (H2, CO, syngas, or other renewable fuels) through the electrolysis of H2O and/or CO221–24 while maintaining similar advantages such as high energy conversion efficiencies, high production rates, fast power cycling, and low manufacturing costs.
Despite these numerous advantages and broad application prospects, the large-scale application of SOC systems is severely inhibited by the high operating temperature (e.g., 800 to 1000 °C).18,25 Although such high operating temperatures can significantly accelerate the reaction kinetics and ionic conductivity, they can lead to detrimental interactions between cell components, leading to poor durability and high costs of utilization.3,26,27 Therefore, to develop more efficient, durable, and economically feasible SOC systems, the SOC operating temperature must be reduced to about 500 °C while maintaining sufficiently high performance.28–30 As the operating temperature is reduced, however, the cell performance will decrease exponentially due to slower oxygen reduction/evolution (ORR/OER) reaction kinetics on the oxygen electrode,31–35 lower ionic conductivity of the electrolyte,36 and increased polarization resistance of the fuel electrode.37 Therefore, new materials or novel structures for SOC components need to be designed to enable high ionic conductivity in the electrolytes, high electrocatalytic activity at the electrodes, and good compatibility with other cell components under operating conditions.38,39
Overall, these results demonstrate that heterointerface engineering is a promising approach to enhancing the SOC component performance at intermediate temperatures. Such enhancements can be attributed to many factors, including lattice strain, dislocation, cation segregation, changes in electronic structure, and expansion of triple phase boundaries (TPBs) (Fig. 2).4,50,51 For example, research into electrolytes52,53 and oxygen electrodes54,55 suggested that tensile strain at heterointerfaces can accelerate ionic transport by decreasing the migration barrier of charge carriers. In addition, results showed that the electronic structure near the heterointerface of LSC113/LSC214 oxygen electrodes can be different from their single-phase counterpart, leading to strongly enhanced ORR kinetics.56 The presence of cation segregation near the interfacial region was reported to provide more active sites for the ORR/OER.57,58 Furthermore, Ni/metal oxide interfaces of fuel electrodes have been reported to have higher coking resistance by facilitating water adsorption and water-mediated carbon removal.59,60 To guide the rational design of high performance SOCs, it is imperative to fully understand the role of heterointerfaces in each of the key cell components: the electrolyte, the oxygen electrode, and the fuel electrode.
It is noted that the factors that affect the properties of heterostructured SOC components (Fig. 2) are often inter-related. For examples, small lattice mismatches at a heterointerface may be readily accommodated by local lattice strains;52,53 when the lattice mismatches are sufficiently large, however, the formation of dislocations may become energetically favorable to relieve local strain and stabilize interfacial structures. Furthermore, it has been proven that lattice strain may influence defect formation energies and migration barriers,54,55 thus altering the distribution of charge carriers near heterointerfaces or in the space charge region (e.g., cation segregation at the interfacial region57,58). The presence of or a change in lattice strain, defects (e.g., oxygen vacancies or dislocations), and redistribution (accumulation or depletion) of charge carriers in space-charge regions can induce electron transfer and influence the electronic structure56 (or the energy landscape) of the interface, which determines the fundamental properties of the heterointerfaces. However, the specific correlations among these factors are yet to be fully understood.
This review aims to provide a comprehensive and in-depth overview of the impacts of heterointerfaces on the ionic transport, catalytic activity, and stability of heterostructured materials, providing important insights into the mechanisms of performance enhancement and establishing a scientific basis for the rational design of better electrolyte and electrode materials. The fundamentals of heterointerfaces are reviewed in Section 2, including basic concepts, structural characteristics, and controlled fabrication processes. Subsequently, applications of heterostructures in electrolytes, oxygen electrodes and fuel electrodes of SOCs are highlighted, along with possible mechanisms for the performance enhancement in those heterostructured SOC components, including lattice strain, dislocation, cation segregation, modified electronic structure, and extended TBPs at heterointerfaces. In addition, new research directions, future perspectives, existing challenges and possible solutions for the rational design of novel materials for advanced solid-state electrochemical systems involving heterostructures will be summarized. Overall, this review will provide a useful guide for the rational design of novel materials with unique functionalities for high-performance SOCs. The knowledge and methodology illustrated in this review are also suitable for the development of novel catalysts for other energy storage and conversation systems.
With technological developments and in-depth studies of heterojunctions, the definition and research of heterostructures have largely expanded. The extraordinary properties associated with heterointerfaces include high electrical conductivity,36,62 enhanced thermoelectric effects, boosted catalytic activity,63 super-wettability,64,65 and induced magnetization.66–69 The first study on the enhancement of ionic conductivity in composite materials was based on inorganic salt/Al2O3 and AgI/AgBr heterointerfaces in the 1980s.70 Subsequently, significantly enhanced ion conduction was also found in multi-layered CaF2/BaF236,71 and YSZ/STO epitaxial heterostructures.72 In addition, high-mobility electron gases at the heterointerface of LaAlO3/SrTiO3 were investigated by Ohtomo et al.73,74 and Nakagawa et al.75, which demonstrated extremely high carrier mobility exceeding 10000 cm2 V−1 s−1 as a result of the built-in polarity discontinuity at the interface. These remarkable enhancements in conductivity and activity inspired the application of heterostructured oxides in SOCs.
(a) The multilayer structure is a classical model for theoretical investigations in which heterointerfaces between two adjacent layers are parallel to the surface. Here, well-ordered interface structure enables the precise control of interfacial properties such as local compositions, lattice strains, electronic structures, etc. And in particular, the interfacial properties can play a dominate role in determining the functionality of the multilayer when the thickness of each layer approaches several nanometers.
(b) The VAN structure is another widely used model system for the investigation of the role of interfaces in the functionality of materials. Compared with the multilayer structure model, the unique advantage of the VAN structure is the maximum exposure of the active heterointerface. Furthermore, the interfaces in VAN structures are either directly exposed to gas (as electrodes) or along the direction of ion flow (as electrolytes), enabling full utilization of the heterointerfaces in SOCs. In addition, a VAN structure can sustain higher interfacial strain compared with a multilayer structure.
(c) Composite structures are often used in practical solid oxide cells; the heterointerface in this structure is random and the actual contact of the two phases is complicated. It is difficult to quantify the interfacial region between the two phases.
(d) The skeleton/surface coating structure is composed of a porous backbone decorated with a thin-film coating, which is also widely used in practical SOCs. In this structure, the backbones (or the materials serving as the mechanical support) generally possess good ionic or mixed conductivity, whereas the coating materials are generally highly active catalysts.
To be specific, PLD can be used to fabricate complex heterostructured thin films with desired stoichiometric ratios, in which a high-power pulsed laser beam is produced by an excimer laser (e.g. KrF) and focused inside a vacuum chamber through a lens to a target (Fig. 4a).103 Here, the target material is evaporated as a plasma plume in a high vacuum or in the presence of a background gas (e.g. oxygen of various concentrations) and subsequently deposited onto a substrate. And by optimizing operating conditions such as temperature, oxygen partial pressure, laser energy, etc., thin films with well-ordered crystal lattice units and specific orientations can be prepared. In addition, the use of PLD can produce various heterostructured oxide morphologies including bilayer,106 ML107 and VAN77 structures. For example, Schichtel et al.108 fabricated YSZ|Sc2O3 multilayers using PLD to study strain effects on the ionic conductivity.
Fig. 4 Schematics of a (a) PLD system,103 (b) MBE system104 and (c) up: tubular ALD reactor, down: layout of a plasma-enhanced ALD cycle for Cu deposition.105 Reproduced from ref. 103–105 with permission from RSC, copyright 2016; Baiutti et al; license Beilstein-Institute (https://www.beilstein-journals.org/bjnano/articles/5/70), copyright 2014 and American Chemical Society, 2015, respectively. |
The main advantage of the MBE method (Fig. 4b) is the excellent control of oxide growth as a result of the single-element fluxes, minimized kinetic energies and slow kinetics.104,109,110 The difficulty of the MBE method, however, is the inability to maintain desired cationic stoichiometry during growth, which requires cooperation with an in situ monitoring auxiliary system, including a quartz crystal microbalance (QCM), real-time absorption spectroscopy (AAS), reflection high-energy electron diffraction (RHEED) etc. Furthermore, applications of MBE and PLD are also limited by their expensive deposition systems.
In addition, ALD is applied in layer-by-layer structure fabrication, which is based on the gas–solid reaction between the atmosphere and a substrate (Fig. 4c).105,111 ALD can produce conformal and pinhole-free coatings of porous structures and is usually used in the fabrication of binary or ternary oxides.112,113 For example, high quality films of many single oxides (e.g., MgO,114 NiO115 and CuO116) can be readily coated onto perovskite LSCF backbones using ALD.
Electrolytes are key components of fuel cells, which largely determine SOC electrode compositions.3,88,121,122 To improve the performance of SOCs at intermediate temperatures, the electrolyte materials should meet the following requirements: (1) high ionic conductivity and negligible electronic conductivity under operating conditions; (2) good stability in various atmospheres over a wide range of temperature; (3) a small mismatch in the thermal expansion coefficient with electrodes and other cell components, (4) negligible interactions with electrodes and interconnect materials under operating and processing conditions, and (5) adequate mechanical strength or integrity.
High ionic conductivity of SOC electrolytes requires a sufficient concentration of mobile charge carriers with high mobility. Based on the types of ionic charge carriers, SOC electrolytes can be classified into oxygen ion conductors, proton conductors, and mixed ion conductors. To date, tremendous efforts have been contributed to the development of high-performance electrolyte materials with oxygen ion conductivity,123,124 including doped zirconia,125,126 doped ceria,125,127,128 lanthanum gallate-based oxides,129,130 and bioxide based materials.131,132 The activation energy for O2− transport in oxygen ion conductors is usually in the range of 50–60 kJ mol−1,133,134 indicating significantly lower ionic conductivity at reduced temperatures.
In addition, since the discovery of alkali earth doped cerates (ACeO3, A = Sr and Ba) with perovskite structures as promising proton conductors,135,136 researchers have devoted much attention to these materials.137–141 The general formula of perovskite type proton conductors can be expressed as AB1−xMxO3−δ, in which M denotes a trivalent dopant and δ represents the oxygen deficiency per unit cell.139,142 Here, the mechanisms of proton conduction include proton transfer, structural reorganization, diffusion motion of extended moieties, etc.135,136,139,143
Fig. 5 Schematic of the functional principles for heterostructured electrolytes. Insets: (left) Schematic view illustrating the distortion of the unit cell under lattice strain;145 activation energies for vacancy migration in CeO2 biaxially strained along [100] and [010];157 oxygen vacancy formation energies at the MgO layer in a TiO2-terminated SrTiO3 interface;158 lateral ionic conductivity versus inverse temperature for YSZ films deposited on different substrates and with film thicknesses in the range 420 to 58 nm;159 charge distribution in the STO/LAO/STO heterostructure;160 normalized vacancy concentration as a function of distance from the GB interface at 440 °C161,162 (right). Reproduced from ref. 145 and 157–161 with permission from RSC, copyright 2017; RSC, copyright 2012; Springer Nature, copyright 2014; WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim, copyright 2010; Springer Nature, copyright 2018 and RSC, copyright 2015, respectively. |
Material | Geometry | T (°C) | Relative increase in conductivity | Mechanism | Ref. |
---|---|---|---|---|---|
GdCeO/MgO | Bilayer | 500–800 | ×3 | Heat treatment that releases strain energy and removes dislocation | 148 |
YSZ/MgO | Bilayer | 400–800 | ×102.5 | Lattice mismatch strain | 43 |
YSZ/MgO YSZ/STO | Bilayer | 150–500 | ×103 | Combination of mismatched dislocation density and elastic strain | 159 |
Ce0.9Gd0.1O2−δ/STO | Bilayer | 450–900 | ×10–102 | Compressive strain | 163 |
YSZ/Al2O3 | Bilayer | 27–377 | ×100.5–1.5 | Tensile strain and edge dislocation | 164 |
BZY/NGO | Bilayer | 550–600 | ×102 | Dislocation | 144 |
Ce0.8Gd0.2O2−δ/Si3N4 | Bilayer | 300–500 | ×10 | Strain | 165 |
GDC/Si3N4/Si | Bilayer | 327 | ×102 | Strain | 145 |
YSZ/LAO, NGO, MgO, Al2O3 | Bilayer | 300–500 | ×10−2 | Film texture and grain boundary density | 166 |
SDC/STO/BZO/MgO | Bilayer | 600 | ×3 | Strain | 146 |
YSZ/Y2O3, Sc2O3, Lu2O3 | Multilayer | 610 | ×2–10 | Strain | 167 |
YSZ/CZO | Multilayer | 650 | Almost 0 | Strain | 168 |
GDC/ZrO2 | Multilayer | 377 | 18 (compared to YSZ single crystal) | Strain and dislocation | 169 |
CSZ/Al2O3 | Multilayer | 350–700 | ×102 | Structural disorder/mismatch | 150 |
YSZ/STO | Multilayer | 84–258 | ×108 | Atomic reconstruction | 72 |
YSZ/Y2O3 | Multilayer | 350–700 | ×10 | Dislocation | 170 |
YSZ/STO | Multilayer | 357–531 | — | p-Type conductivity | 171 |
YSZ/STO | Multilayer | 100–500 | ×105 compared to YSZ/Al2O3 | Cation interdiffusion | 172 |
(all electronic cond.) | |||||
Ce0.8Sm0.2O2−δ (SDC)/YSZ | Multilayer | 400–800 | ×10 compared to SDC | Tensile strain | 173 |
×10 compared to YSZ | |||||
YSZ/STO | Multilayer | Room T–500 | ×102 | Cation interdiffusion | 174 |
Sc2O3/YSZ | Multilayer | 420–780 | ×105 total | Strain | 108 |
YSZ/CeO2 | Multilayer | 400–700 | Almost 0 | Strain | 125 |
ZrO2/STO | Multilayer | — | Fluorite not stable for ε > 0.05 | Strain | 151 |
SrTiO3/YSZ/SrTiO3 | Multilayer | 227 | Almost ×103 | Defect redistribution | 175 |
YSZ/Gd2Zr2O7 | Multilayer | 250–475 | ×102 | Tensile strain | 176 |
YSZ/Y2O3 | Multilayer | 520 | ×102 | Strain | 177 |
Sc2O3/YSZ | Multilayer | 420–780 | ×105 total | Strain | 108 |
CeO2/Y2O3, La2O3, Gd2O3 | Multilayer | 610–1000 | ×10 | Lattice mismatch strain | 149 |
GDC/YSZ | VAN | 600 | ×2 | Strain and fast ionic transport rates | 46 |
YSZ/SDC, STO | VAN | 360 | ×102 | Strain and structurally mismatched vertical interface | 152 |
SrTiO3/SDC | VAN | 400 | ×10 | High crystallinity nanopillars | 153 |
YSZ/STO | Composite | 127–727 | ×101.5 @ 1000 K ×103.5 @ 400 K | Lattice strain | 53 |
SDC/LiNaCO3 | Composite | 300–650 | ×10 | Surface properties and electrolyte thickness | 156 |
CeO2/M2O3 | Composite | 227 | ×104 | Strain | 157 |
YSZ/STO | Composite | 200–600 | ×10−2 | Defect formation and transport, p-type conductivity | 178 |
YSZ | Calculation | 137–927 | ×10 | Strain | 179 |
BZY/MgO | Calculation | 300–700 | Almost 0 | Segregation tendency, but saturated | 180 |
ZrO2/Cr2O3 | Calculation | 527–927 | — | Defect redistribution | 181 |
Many multilayer structured oxides have been reported to show high conductivity and low activation energies due to the presence of heterointerfaces. For example, Azad et al.169 synthesized a multilayer structure consisting of Gd2O3-doped CeO2 and ZrO2 films grown on Al2O3(0001) using MBE and reported that the conductivity of the 10-layer GDC/ZrO2 heterostructure was 2.9 × 10−4 S cm−1 at 650 K in air, which was an ∼18 times enhancement as compared with YSZ single crystals. In another example, Peters et al.150 found that in CSZ (ZrO2 + 8.7 mol% CaO)/Al2O3 multilayer systems with coherent interfaces, the activation energy for ionic conduction was decreased from 1.47 eV to 0.96 eV as the number of CSZ layers increased from 1 to 5. There was a two orders of magnitude enhancement in the total conductivity when the individual CSZ layers got thinner from 0.78 μm to 40 nm (Fig. 6a and b). Sillassen et al.159 deposited YSZ thin films onto MgO(110)(111) substrates and observed remarkable conductivity enhancements; the ionic conductivity of the resulting YSZ thin film was increased to 1 S cm−1 at 500 °C, nearly 3.5 orders of magnitude higher than that of bulk YSZ.
Fig. 6 (a) TEM micrograph of a CSZ/Al2O3 multilayer system after heat treatment at 800 °C for ∼100 h. (b) Arrhenius plot of the total conductivities (σtot) measured on different CSZ/Al2O3 multilayer systems. Reprinted with permission from ref. 150. Copyright 2007 Elsevier. (c) High-resolution TEM image of a vertical SDC–STO interface in a cross-sectional view. Scale bar, 10 nm. (d) Temperature dependence of the real part of the ac conductivity (σac′) measured in the nanoscaffold SDC–STO with comparisons with the single phase thin film and literature data.153 Reprinted with permission from ref. 153. Copyright 2015 Springer Nature. (e) Z-STEM image of the BZY thin film and its interface with the NGO substrate, inset: Z-dependent intensity profile of the cations taken along the dashed line in the (−1 0 1) direction. (f) Arrhenius plot of the conductance for the whole measured set of film thicknesses.144 Reprinted with permission from ref. 144. Copyright 2014 AIP publishing. |
Compared with multilayer structures, vertically aligned nanocomposite (VAN) heterostructure films can potentially accommodate larger strain near the interface, creating fast ionic diffusion paths along the vertical interface.78,152 For example, Yang et al.153 prepared highly crystalline micrometer-thick vertical nanocolumns on a Nb:STO substrate using PLD (Fig. 6c) and reported that this SDC–STO VAN structure exhibited significantly enhanced oxygen ion conductivities, approaching 3 × 10−2 S cm−1 at 400 °C as recorded using scanning probe microscopy measurements. In contrast, the conductivity of SDC, YSZ, and STO single phase thin films is on the order of 10−3, 10−4, and 10−6 S cm−1, respectively, under the same conditions (Fig. 6d). Here, the activation energy for oxygen migration in the SDC–STO VAN structure (0.65 ± 0.02 eV) was also slightly lower than that of the SDC film (0.70 ± 0.04 eV). Lee et al.152 also demonstrated that a double-layered vertical YSZ–STO/SDC–STO heterostructured film prepared by PLD can provide a high ionic conductivity of ∼10−2 S cm−1 at 635 K, which was over 2 orders of magnitude higher than that of YSZ films (∼6 × 10−5 S cm−1). In addition, Su et al.46 reported that in GDC/YSZ VAN systems, two-phase strain coupling can improve the ionic conductivity along the vertical heterointerface. The ionic conductivity at 600 °C of the GDC/YSZ VAN thin film (0.96 S cm−1) was about 2 times higher than that of single phase YSZ thin films (0.41 S cm−1) and slightly larger than that of the GDC/YSZ multilayers (0.82 S cm−1), as determined from impedance spectroscopy. As a result, the corresponding single cell (Ni–YSZ|YSZ/GDC/VAN|LSC) exhibited peak power densities (PPD) of 0.488, 0.694 and 0.883 W cm−2 at 700, 750 and 800 °C, respectively.
Similarly, heterostructured proton conductors have been reported to exhibit higher proton conductivity than those of the corresponding single phases. For example, Foglietti et al.144 reported that the conductivity of a 10 nm thick film was about 2 orders of magnitude larger than that of a 500 nm thick film (Fig. 6e and f) and the average in-plane conductivity of the 10 nm thick film was 20 S cm−1 at 550–600 °C, as demonstrated by a series of strained BaZr0.8Y0.2O3−δ (BYZ) thin films of different thicknesses deposited on NdGaO3 (NGO) substrates using PLD. The enhancement in proton conductivity was attributed to fast ionic transport induced by dislocations along the heterointerfaces between dissimilar materials. In another example, Yang et al.182 characterized the electrochemical activities and transport properties of BaZr0.8Y0.2O3−δ (BZY)/NdGaO3 (NGO) heterointerfaces using electrochemical strain microscopy (ESM), demonstrating that the mismatch dislocation network reduced the activation energy (0.1 eV for the 20 nm thick film) and induced a 2D transport phenomenon, which will be discussed in detail in Section 3.2.2.
Tensile elastic strain is generally believed to accelerate ionic transport by decreasing the energy barrier of migration or by reducing the activation enthalpy for charge-carrier migration.52,53,157,187 For example, Kushima et al.53 investigated changes in oxygen ion migration barriers in YSZ under biaxial lattice strain using DFT and the nudged elastic band (NEB) method. They found that the oxygen migration pathways may diversify and coexist depending on the local distribution of defects such as vacancies and dopant cations along the migration path. Here, the researchers reported that at lower stress levels, tensile strain can expand the migration space, increase the distance between cations (Δ(Zr–Zr)) among the migration path in this configuration and weaken O–C bonds (Fig. 7), all of which can decrease vacancy migration barriers and exponentially increase the oxygen diffusivity. Specifically, under 4% biaxial tensile strain, a maximum 6.8 × 103 times diffusivity enhancement was obtained at 400 K for YSZ layers.
Fig. 7 (a) Change in the cation–cation distance Δ(Zr–Zr) as a function of tensile strain from 0.00 to 0.08; and (b) simulated bond strengths for O–C in the migration process under biaxial lattice strain. Reprinted with permission from ref. 53. Copyright 2010 the RSC. |
Furthermore, Souza et al.157 used quantitative equations to describe the impacts of stress on the activation enthalpy and ionic conductivity in fluorite structured solid solutions such as ZrO2–M2O3 and CeO2–M2O3. Here, the researchers used static atomistic simulations to calculate the energy change along oxygen vacancy migrations in strained CeO2 based on empirical pair-potentials (EPP) and evaluated the correlation between conductivity and strain using the equations:
(1) |
(2) |
Based on the model described by eqn (1) and (2), Korte et al.167 further established a computational model by simulating elastically deformable crystallites. And based on grain boundary dependent strain relaxation in solid electrolyte films, their equation can adequately approximate the relationship between the total conductivity and structural parameters when the lattice mismatches are below 10%:
(3) |
Subsequent experimental measurements were also consistent with these theoretical models and highlighted the dependence of ionic conductivity on lattice strain.108,159,189 Schichtel et al.108,189 systematically investigated the influence of interfacial stain effects at semi-coherent and coherent heterointerfaces such as YSZ(111)/RE2O3(111)/Al2O3(0001) and Al2O3(0001)/Sc2O3(111)/YSZ(111) on ionic transport along the interface. In these studies, the researchers reported that the conductivity decreased with increasing interfacial density and attributed this to the compressive strain in YSZ at the interface. This was also in agreement with model (3). Furthermore, Sillassen et al.159 suggested that eqn (3) more accurately explained strain effects on continuous, flat and coherent interfaces (or with slight dislocations).
Overall, interfacial tensile strain generally increases the ionic conductivity and decreases ionic migration barriers whereas compressive strain has the opposite effect. However, lattice mismatch for heterointerfaces can lead to tensile strain on one side and compressive strain on the other, in which heterointerfaces coupled with two opposite effects simultaneously can face uncertainty and complexity. Therefore, the aim of heterointerface engineering is to maximize positive effects through the optimal design of the morphology and selection of material systems.
Fig. 8 (a) HRTEM micrograph of a Sc2O3/YSZ heterointerface with a regular network of dislocations (white marks).108 Reprinted from ref. 108 with permission from RSC, copyright 2013. (b) HRTEM micrograph of a YSZ/MgO interface with mismatch dislocations marked by T symbols. The dotted rectangle exhibits a Burgers circuit in which the difference between atomic steps on each side represents the quantity of dislocations in each of the two orientations.159 Reprinted from ref. 159 with permission from WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim, copyright 2010. |
Recent studies also suggested that dislocations can critically impact the ionic transport properties.158,159,182,191–193 For example, Sillassen et al.159 found that the conductivity of YSZ/MgO systems with 58 nm-thick YSZ films can be enhanced by almost 3.5 orders of magnitude with a large lattice mismatch of −18.01% between the YSZ and MgO lattices. They assumed that the dislocation network along the semi-coherent YSZ/MgO interface (Fig. 8b) may play a role of fast transport pathways for oxygen diffusion.
In another study on the MgO/SiTiO3 system, Dholabhai et al.158 used atomistic simulations with LAMMPS and proposed that dislocation networks in the TiO2-terminated interface can reduce oxygen vacancy formation energies and accelerate ionic transport. Here, the researchers reported that oxygen vacancies preferred to form along dislocations lines with lower formation energies, which was similar to those in metallic nanocomposites.194 They also suggested that TiO2-terminated interfaces can exhibit analogous pipe diffusion mechanisms to metallic nanocomposites,195 which were responsible for remarkable ionic conductivities. In comparison, SrO-terminated interfaces possessed oxygen vacancies that were concentrated on terraces and were inhibited from migrating.
Based on TEM analysis and calculation, Foglietti et al.144 suggested that the enhanced in-plane proton conductivity was due likely to the high densities of defects, particularly dislocations, at the interface of the film. Similarly, Yang et al.182 proposed that the 2D dislocation network near the interfacial region can facilitate the proton conductivity by lowering the activation energies (Fig. 9a). STEM characterization identified mismatch dislocations spaced −10 unit cells apart at the BZY–NGO interface. Here, the atomic arrangements were highly disturbed with more oxygen vacancies along the dislocation network, allowing the redistribution of defects to induce charge imbalance and attract proton incorporation near the interface. In addition, they reported that the local distortion of oxygen octahedra can reduce the trapping effect, allowing the 2D dislocation network at the BZY/NGO interface to provide easily accessible sites to anchor H+ and decrease activation energies.
Fig. 9 (a) Schematic of the dislocation network at the interface of a BZY thin film and NGO substrate based on STEM characterization results.182 (b) Map of oxygen vacancy formation energies around one dislocation core.197 Reprinted with permission from ref. 182 and 197, both with permission from American Chemical Society, 2015. |
Despite many studies reporting positive impacts of dislocations on the ionic conductivity, some have also reported negative impacts. For example, Sun et al.196 used atomic simulations and found that strain fields can enrich the concentration of trivalent dopant cations and oxygen defects near the 1/2 〈110〉 {100} edge dislocations in CeO2 with depletion at the other edge. In the enriched zone, strong dopant–vacancy and vacancy–vacancy interactions can diminish oxide ion diffusion, whereas in the depleted zone, the lack of charge carriers can decrease the conductivity. In another study, Marrocchelli et al.197 found that the redistribution of defects resulted in overlapping electrostatic fields surrounding the dislocation cores with positive charge, which further influenced the electrical and ionic conductivity. And although the oxygen vacancy formation energies were lower than 2 eV at sites close to the dislocation cores of 〈100〉 {011} as compared with the SrTiO3 bulk (Fig. 9b), these researchers did not find any significant enhancement in the oxide ion mobility. Similarly, Schichtel et al.108 reported the negligible role of dislocations in the ionic conductivity along the heterointerface of Al2O3(0001)/Sc2O3(111)/YSZ(111). Yep et al.198 also reported that edge dislocations induced by kinking on the (220) slip plane near the YSZ/quartz interface can inhibit oxygen ion migration along the interfacial region. Overall, the complex role of dislocations in the ionic diffusion kinetics is still not well understood and requires further systematic research.
(4) |
The positive impact of space-charge effects on the ionic conductivity at interfacial regions has been highlighted in several studies.43,204,205 For example, Kosacki et al.43 proposed that due to space-charge effects, the ionic conductivity of 1.6 nm-thick YSZ/MgO interface layers can reach ∼2 S cm−1 at 700 °C, which was 3 orders of magnitude larger than that of the lattice. In another example, Karthikeyan et al.205 reported that in YSZ/MgO, the ionic conductivity increased as the film thickness decreased. They also found that the high activation energy indicated that the space-charge effect at the interface and grain boundary dominated the conductivity enhancements. Despite these promising results, however, a straightforward relationship between the space-charge effect and ionic enhancement still needs to be further investigated. Thus far, there were only a few reports about tuning the space-charge potential and defining the vacancy accumulation regions clearly.206–208 For example, Gäbel et al.207 reported that a nanocrystalline GDC thin film on a MgO substrate can exhibit relatively high ionic conductivities (2.7 × 10−4 S m−1 S at 700 °C in 10−5 bar O2) with a low space-charge potential of 0.19 ± 0.05 V, in which the researchers attributed this to the low depletion of oxygen vacancies at the grain boundaries.
Alternatively, there were studies demonstrating that the impacts of space-charge effects on conductivities are complicated. Some researchers found that the positively charged grain boundary cores can block the migration of oxygen vacancies, significantly influencing the electronic and ionic conductivities.209,210 For example, Adepalli et al.211 reported that dislocation-induced space-charge fields can enhance the n-type conductivity of SrTiO3 by 50 times at 10−5 atm p(O2) at 650 °C and decrease the p-type conductivity by 50 times and reduce the oxygen diffusion coefficient by three orders of magnitude at high oxygen pressures. In another study, Lupetin et al.212 reported that in SrTiO3, a reduction of the grain size to 30 nm can bring out a 3 orders of magnitude enhancement of the n-type conductivity and a significant decrease in the p-type conductivity and oxygen vacancy conductivity. The researchers attributed such decreases to increased space-charge effects near the interface as the particle size decreases. Similar effects have also been found for nanocrystalline ceria.213
Furthermore, Mebane et al.161 combined the Poisson–Boltzmann with the Cahn–Hilliard model to estimate space-charge effects in solid solutions with wide dopant levels. The ‘Poisson–Cahn’ theory is based on local and non-local chemical interactions and can yield activity coefficients for point defects. And in the case of heterostructured CeO2–Gd2O3, the model successfully predicted defect behaviors near the grain boundary over the entire concentration range. In addition, the model, included existing conductivity considering local dopant concentrations, covered various conditions of solid solutions, and finally gave a comprehensive prediction of both bulk and grain boundary conductivities.
In this section, we have discussed the role of heterointerfaces in electrolyte materials for SOCs. In summary, heterointerfaces can greatly influence the transfer properties of charge carriers (such as oxygen ions and protons) via various mechanisms including lattice strain, dislocation and space-charge effects. Actually, the evolution of the local lattice structure, cation arrangement, and electronic structure induced by heterointerfaces may not only affect the mass transfer, but also have a more profound influence on the primitive steps of electrode reactions, which will be discussed in detail in the next sections.
ORR/OER processes on mixed ionic and electronic conductor (MIEC) electrode materials involve several elementary steps. 215 The detailed ORR process includes the following steps (Fig. 10):
Fig. 10 Schematic of the elementary steps of an oxygen reduction reaction process on a heterostructured oxygen electrode. |
① Gas diffusion, O2(gas) → O2(near surface);
② Adsorption, O2(near surface) → O2(ads);
③ Dissociation, O2(ads) → 2O(ads);
④ Electron transport, O2(ads) + e− → O2−(ads), O2−(ads) + e− → O22−(ads), O22−(ads) → 2O−(ads), O(ads) + e− → O−(ads), O−(ads) + e− → O2−(ads);
⑥ Diffusion of O2− to the electrode/electrolyte interface;
⑦ O2− transport across the electrode/electrolyte interface.
Many of the elementary reaction steps above involve the transportation of ionic and electronic species. Therefore, optimal ionic and electronic conductivities as well as surface exchange coefficients are desirable for oxygen electrodes. Furthermore, oxygen electrode materials should meet additional requirements such as high chemical stability, compatibility with electrolyte materials, thermodynamic properties matching other materials, etc.
Material | Geometry | T (°C) | Performance | Mechanism | Ref. |
---|---|---|---|---|---|
La0.8Sr0.2CoO3−δ/LaSrCoO4+δ | Multilayer | 550 | k q (from EIS) ∼103–104 enhancement | Strain, space-charge effects and an increase in electronic structure or oxygen vacancies | 48 |
La0.8Sr0.2CoO3−δ/YSZ | Multilayer | 510 | Surface exchange coefficient ∼10 times | Sr decoration | 226 |
(La0.5Sr0.5)2CoO4±δ/La0.6Sr0.4CoO3−δ | Multilayer | 550 | k q ∼ 3 orders of magnitude enhancement | Heterointerface | 222 |
SrTi1−xFeO3/YSZ | Multilayer | 345 | Band gap increased from 2.5 ± 0.5 eV to 3.6 ± 0.6 eV | Surface electronic structure | 227 |
La1−xSrxCoO3−δ/SrTiO3 | Multilayer | 280–475 | Tensile strain with faster surface exchange (∼4 times) and diffusion (∼10 times) | Strain | 228 |
La1−xSrxCoO3−δ/LaAlO3 | |||||
La0.8Sr0.2CoO3−δ/La0.8Sr0.2MnO3−δ | Multilayer | 550 | k q (cm s−1) (from EIS) ∼10–102 enhancement | Mn substitution (cation diffusion) | 79 |
La0.8Sr0.2CoO3/(La0.5Sr0.5)2CoO4 | Multilayer | 250 | Energy gap decreased from 1.4 eV (LSC113) and 2.6 eV (LSC214) to 1.0 eV (LSC113/214) | Electron donation and transfer of oxygen vacancies across the heterointerface | 229 |
(La,Sr)2CoO4±δ/La1−xSrxMO3−δ | Multilayer | 550 | Surface exchange coefficients enhanced by 102 times | Sr segregation | 230 |
SrCoOx/(LaAlO3–(SrAl0.5Ta0.5O3)0.7) or SrTiO3 or DyScO3 or GdScO3 or KTaO3) | Multilayer computation | 300 | Oxygen activation energy barriers decreased by ∼30% | Strain | 231 |
La0.6Sr0.4CoO3−δ/Ce0.9Gd0.1O1.95 | Multilayer | 700 | Slight increase in 18O ratio | Extended TPB | 232 |
La0.8Sr0.2MnO3/SrTiO3 | Multilayer | 600 | Oxygen ion diffusion enhanced by 103 | Dislocation and strain | 192 |
La0.8Sr0.2MnO3/LaAlO3 | |||||
La0.8Sr0.2CoO3−δ/Nd2NiO4+δ | Multilayer | 300–415 | Oxygen exchange kinetics are much slower | Charge transfer | 233 |
La0.6Sr0.4Co0.2Fe0.8O3/PrxCe1−xO2 | Multilayer | 600 | PrO2/LSCF: nearly 6 times faster | Inducing higher oxygen vacancy concentrations | 234 |
CeO2/LSC: slight enhancement | |||||
La0.6Sr0.4CoO3−δ/LaSrCoO4±δ | Multilayer | 600–850 | k q ∼ 5–10 times that of the single phase | Sr enrichment at the interfacial region and stabilized against detrimental Sr segregation | 106 |
Nd0.5Sr0.5CoO3−δ/Nd0.8Sr1.2CoO4±δ | Multilayer | 500 | k q ∼ 102–103 times that of the single phase | Enriched Sr concentration and decrease in the valence state of Co | 107 |
La0.65Sr0.35MnO3/SrTi0.2Fe0.8O3 | Multilayer computation | — | Bandgap values (1 and 2 eV for LSM and STF) decreased to 0.3 and 1 eV for LSM and STF | Electronic structure, strain, defects | 235 |
La0.8Sr0.2CoO3/(La0.5Sr0.5)2CoO4 | Multilayer | 300 | 103 times increase in charge transfer rates to oxygen | Electronic structure | 56 |
La0.8Sr0.2CoO3/HfO2 decoration | Multilayer | 530 | k q ∼ 30 times enhancement | Optimum surface oxygen vacancy concentration | 236 |
La0.6Sr0.4CoO3−δ/Co2O3 decoration | Multilayer | 450 | Surface exchange resistance reduction by 13% per laser pulse | Reactivation of LSC surface | |
La0.8Sr0.2CoO3−δ/LaSrCoO4+δ | VAN | 320–400 | k q ∼ 10 times enhancement | Electronic activation and more stable cation composition at the interface | 77 |
La0.6Sr0.4CoO3−δ/LaSrCoO4+δ | Composite | 500 | Exchange coefficient k* (cm s−1) (from SIMS) enhanced by 103 | Fast oxygen incorporation paths along the heterointerface boundary | 216 |
La0.6Sr0.4CoO3−δ/LaSrCoO4+δ | Composite | 500 | Exchange coefficient σe (S cm−2) (from EIS) 10–101.5 | Heterointerface enhanced the oxygen catalytic reaction | 237 |
Fe2O3–LaCePrOx | Composite | 400–600 | Ionic conductivity increased by 10 times | Heterogeneity and interfacial conduction effects | 238 |
SP Ba0.5Sr0.5(Co0.7Fe0.3)0.6875W0.3125/DP Ba0.5Sr0.5(Co0.7Fe0.3)0.6875W0.3125 | Composite | 550–700 | ASR ∼3 times lower | Cooperative DP/SP effects | 239 |
Er0.4Bi1.6O3−δB/La0.8Sr0.2MnO3−δ | Composite | 550–650 | R p/5 | Particle sizes | 219 |
La0.6Sr0.4Fe0.8Co0.2O3−δ/Gd0.2Ce0.8O1.9−δ–ZrO2 | Skeleton/surface coating | 800 | ∼25% degradation rates after overcoating in 1100 hours of operation | Synergistic function of porosity, mixed conductivity and suppressed Sr enrichment | 240 |
La0.6Sr0.4CoO3−δ/ZrO2 | Skeleton/surface coating | 700 | Lower polarization resistances and 19 times slower degradation rates over 4000 h | Porosity, mixed conductivity and suppressed Sr enrichment | 241 |
La0.58Sr0.4Co0.2Fe0.8O3−δ/GDC | Skeleton/surface coating | 900 | k chem ∼ 6 times higher than pure LSCF | Facilitating dissociative adsorption of O2 | 242 |
La0.6Sr0.4CoO3−δ/Al2O3 | Skeleton/surface coating | 800 | Power output reduced from 1.62 W cm2 to 0.742 W cm2 | Blocking effect of excessive Al2O3 | 243 |
La0.6Sr0.4Co0.2Fe0.8O3−δ/MgO | Skeleton/surface coating | 650 | ASR reducing from 0.49 Ω cm2 to 0.31 Ω cm2 | Promotion of charge transfer process | 114 |
La0.6Sr0.4Co0.2Fe0.8O3−δ/CuO | Skeleton/surface coating | 750 | k chem ∼ 9.3 × 10−5 cm s−1, 3.5 times enhancement | Additional reactions on both the CuO surface and the LSCF–CuO-gas boundaries (3PBs) | 116 |
La0.6Sr0.4Co0.2Fe0.8O3−δ/NiO | Skeleton/surface coating | 800 | k q ∼ 6.9 × 10−4 cm s−1, 20 times enhancement | Promotion of oxygen incorporation process | 115 |
(La,Sr)2FeO4−δ/La0.8Sr0.2FeO3−δ | Skeleton/surface coating | 650–800 | ORR activity enhanced by 10 times | Higher catalytic activity of LSF214 and mismatch between LSF214 and LSF113 | 244 |
La0.8Sr0.2MnO3/Ba0.5Sr0.5Co0.8Fe0.2O3−δ | Skeleton/surface coating | 650–750 | Electrode resistances of ∼0.101, 0.049 and 0.023 Ω cm2 at 650, 700 and 750 °C respectively | Material intrinsic properties and fabrication method | 245 |
La0.84Sr0.16MnO3−δ/Bi1.4Er0.6O3 | Skeleton/surface coating | 750 | 10–102 ionic conductivity | Effective TPB extension | 246 |
PrSrCoMnO6−δ/La0.6Sr0.4Co0.2Fe0.8O3−δ | Skeleton/surface coating | 750 | R p of 0.093 Ω cm2 after running for 500 h with a performance enhancement of 61.7% as compared with blank LSCF | Material intrinsic properties and fabrication method | 247 |
La0.6Sr0.4Co0.2Fe0.8O3−δ/La0.85Sr0.15MnO3±δ | Skeleton/surface coating | 700 | Power density of 655 mW cm−2 at 700 °C | Formation of LSM(C) phase | 248 |
(La0.6Sr0.4)0.995Co0.2Fe0.8O3−δ/La2NiO4+δ | Skeleton/surface coating | 750 | 67% increase in peak power density, a low degradation rate of 0.39% for ∼500 h | Cation segregation of LSCF and favorable acceptance by LNO of Sr/Co doping. | 249 |
La2NiO4+δ/Ba0.5Sr0.5Co0.8Fe0.2O3−δ | Skeleton/surface coating | 600 | ASR decreased to 0.078 Ω cm2 | Shell increased TPB | 250 |
Complete resistivity to CO2 attacks at IT | |||||
Er0.4Bi1.6O3/La0.76Sr0.19MnO3+δ | Skeleton/surface coating | 550–750 | Power density of 162 W cm−2; excellent stability in SOFC, SOEC and reversible SOC operating models for over 200 h | Fast oxygen ion migration at the heterointerface | 6 |
(La,Sr)CoO3−δ/(La,Sr)2CoO4+δ | Computation | 500 | 400 times faster oxygen incorporation kinetics | Anisotropy and strain | 51 |
La0.8Sr0.2CoO3−δ/(La0.5Sr0.5)2CoO4−δ | Computation | 500–600 | Oxygen vacancy concentrations enhanced by 102–102.5 | Sr segregation | 223 |
(La1−ySry)2CoO4±δ/La1−xSrxCoO3−δ | Computation | — | Anomalous Sr segregation | Cation diffusion | 57 |
LaCoO3/LaAlO3 | Computation | — | (100) oriented film possessed better OER performance | Spin-state transition of Co | 40 |
Based on the compositions, heterostructured oxygen electrodes can be sorted into four categories: (1) perovskite/Ruddlesden–Popper (ABO3/A2BO4), (2) perovskites with different cations (ABO3/A′B′O3), (3) perovskite/simple metal oxides, and (4) others, which will be discussed respectively as follows.
Fig. 11 (a) Mapping profile of 18O signals by SIMS at the surface above LSC113/LSC214/CGO after oxygen isotope exchange.47 (b) HAAADF STEM micrograph of the LSC113/LSC214 interface;48 (c) EIS results of the microelectrodes (∼200 μm) for LSC113 films with ∼0.1, ∼0.8, ∼5 and ∼15 nm LSC214 surface decoration and the LSC113 reference at 550 °C at 1% p(O2).48 (d) TEM image and selected area electron diffraction of a LSC113/LSF214 particle.244 Reproduced from ref. 47, 48 and 244 with permission from the Electrochemical Society, copyright 2008; and American Chemical Society publications, copyright 2010 and 2017, respectively. |
The ultra-fast oxygen exchange kinetics near the interface of LSC113/214 was confirmed by Crumlin et al.48 using epitaxial La0.8Sr0.2CoO3−δ (LSC113) thin films with partial decoration of (La0.5Sr0.5)2CoO4±δ (LSC214) (Fig. 11b). Here, the researchers reported that the LSC214 decorated LSC113 presented a 3–4 orders of magnitude enhancement in the surface exchange coefficient (kq) compared with the LSC113 single phase reference, reaching as high as ∼3 × 10−6 cm s−1 in 1 atm O2 at 550 °C (Fig. 11c). Furthermore, Lee et al.230 compared the exchange coefficients for LSC113/LSC214 and LSCF113/LSC214 heterostructures and reported that in 1 atm O2 at 550 °C, both LSCF113/LSC214 (∼7 × 10−8 cm s−1) and LSC113/LSC214 (∼2 × 10−7 cm s−1) exhibited significantly enhanced kq values relative to the LSCF113 reference (∼3 × 10−8 cm s−1) and the LSC113 reference (∼5 × 10−9 cm s−1), respectively.
Recently, Zhao et al.214 also reported that bilayered LSC113/LSC214 was able to maintain superior ORR activities at 873–1173 K relative to single-phase LSC113. At 973 K in air, kq of LSC113/LSC214 was 2.5 × 10−5 cm s−1, which was 10 times higher than that of the LSC113 reference. Moreover, in another study conducted by the same group, Zheng et al.11 fabricated a series of Nd0.5Sr0.5CoO3−δ/Nd0.8Sr1.2CoO4±δ multilayers with the same total thickness and different numbers of interfaces. By conducting 18O isotope exchange experiments (97.1% 18O isotope enriched, 773 K, 15 s) and utilizing SIMS, the researchers found that the 18O− signal strongly concentrated along the NSC214/NSC113 interface. With the increased layers of the heterostructure, kq showed a significant enhancement. Specifically, the bilayer NSC214/113 film (∼5.5 × 10−8 cm s−1) exhibited a ∼41 times higher kq value than that of the NSC113 reference (∼1.3 × 10−9 cm s−1), whereas the 18-layered NSC214/113 film (∼2.9 × 10−7 cm s−1) exhibited a 165 times enhancement.
Hong et al.244 used the infiltration method to fabricate a heterostructured (La,Sr)2FeO4−δ (LSF214)–La0.8Sr0.2FeO3−δ (LSF113) electrode (Fig. 11d). After repeated infiltration of LSF214 onto a LSF113 skeleton (4 times), numerous LSF214 nano-islands were generated, and its chemical surface exchange coefficient (kchem) increased from 1 × 10−5 cm s−1 to 4 × 10−4 cm s−1 accordingly, as calculated by ECR experiments at 750 °C.
Aside from bilayer or multilayer structures, Ma et al.77 also fabricated VAN structured LSC113/LSC214 oxygen electrodes that were coherently grown onto STO(001) and YSZ/GDC(001) substrates using PLD. Results from high resolution STEM-EDX verified the separation of LSC113 and LSC214 columns in near-vertical orientation with a ∼300 nm average size of each column. Here, the ASR of the VAN thin-film electrode (∼2 × 104 Ω cm2) was only ∼10% of LSC113 (∼1.6 × 105 Ω cm2) and LSC214 (∼3.5 × 105 Ω cm2) single phase electrodes in air at 320 °C.
ABO3/A2BO4 electrodes with different elemental composition of the ABO3 and A2BO4 components have also been reported to exhibit enhanced performance. For example, Zhang et al.249 demonstrated that the infiltration of La2NiO4+δ (LNO) into (La0.6Sr0.4)0.995Co0.2Fe0.8O3−δ (LSCF) can decrease the ASR value from 1.34 Ω cm2 to 0.042 Ω cm2 at 750 °C with the activation energy changing from 1.38 eV to 1.06 eV. Here, the resulting cell with the LSCF–LNO cathode demonstrated a 67% increase in peak power density (697 mW cm−2) relative to the cell with LSCF at 750 °C. And after long-term operation at a constant current density of 250 mA cm−2, the degradation rate was reported to be as low as 0.39% after 500 hours’ testing.
Fig. 12 (a) kq of LSC82 with a thin LSM82 coating (0.9 nm) and the LSC82 reference as a function of annealing time.79 (b) Comparison of the electrochemical activities of LSCF decorated with CeO2, PCO and PrO2.234 (c) Schematics of a LSCF backbone and conformal, dense PNM decoration as well as exsoluted PrOx nanoparticles.256 (d) Typical I–V–P curves and (e) long-term stability performance (at a constant cell voltage of 0.7 V) for Ni–YSZ anode supported cells with different cathodes: bare LSCF, PNM-, PrOx- and hybrid PNM–PrOx catalyst-decorated LSCFs with 3% humidified H2 as the fuel and air as the oxidant at 750 °C.256 Reproduced from ref. 79, 234 and 256 with permission from American Chemical Society, copyright 2014; American Chemical Society, copyright 2018 and RSC, copyright 2017, respectively. |
Lynch et al.248 also reported that the electrochemical properties of porous La0.6Sr0.4Co0.2Fe0.8O3−δ (LSCF) cathodes can be enhanced through a uniform La0.85Sr0.15MnO3±δ (LSM) decorating layer using the infiltration method. The current densities improved from 0.7 A cm−2 to 1.1 A cm−2, resulting in enhanced peak power densities (655 mW cm−2 at 700 °C) in 200 hours’ operation. Results from TEM and convergent beam electron diffraction (CBED) for the aged oxygen electrodes showed that the LSCF skeleton retained its perovskite structure after operation, and the LSM surface coatings exhibited a loss in crystallinity. In another example, Zhu et al.245 fabricated a (Ba,Sr)(Co,Fe)O3 porous cathode and decorated it with a (La,Sr)MnO3 nanofilm through solution impregnation. The polarization resistance of the LSM-decorated BSCF was found to be 0.023 Ω cm2 at 750 °C in air, which was ∼20% lower than that of the undecorated cathode. The researchers also reported that the LSMO layer on the BSCF backbone prevented inherent electrode poisoning by CO2, H2O and SO2 in air, resulting in considerable stability enhancements. Moreover, Wang et al.252 reported that the infiltration of perovskite SrCo0.9Nb0.1O3−δ (SCN10) onto LSCF can lead to an 89.4% reduction in ASR, indicating enhanced ORR/OER kinetics.
Tsvetkov et al.236 conducted systematic research into the decoration of less reducible cations on perovskite LSC, including V2O5, Nb2O5, TiO2, ZrO2, HfO2 and Al2O3, and concluded that related metal oxide decoration can improve the surface exchange kinetics with the exception of V2O5. A ‘volcano’ relation between kq and the oxygen vacancy formation enthalpy (ΔHVf(MeOx)) was proposed in which optimal kq enhancement (∼30 times) was achieved through HfO2 addition after 54 h in air at 530 °C.
Furthermore, Rupp et al.253 used in situ impedance spectroscopy during pulsed laser deposition (IPLD) to investigate the electrochemical activity of LSC during the deposition of well-defined monolayer fractions of Sr-, Co- and La-oxides. Here, the researchers reported that Co2O3 decoration using one laser pulse reduced the surface exchange resistance (Rsurfexch) by 13%, whereas SrO decoration using one laser pulse (∼13 pmol Sr) increased Rsurfexch by 42 ± 6%. The researchers in this study attributed this difference to the fact that Co decoration can induce more active sites whereas Sr decoration can deactivate LSC surfaces.
The decoration of Zr-based and Ce-based oxides onto perovskite surfaces has also been reported to improve the reaction activity and operation stability of oxygen electrodes.240–242 For example, Gong et al.241 decorated LSC surfaces with less electrocatalytically active ZrO2 using atomic layer deposition (ALD) and reported that the ZrO2 overcoating gradually became porous and mixed conducting after thermal exposure, resulting in elevated ORR properties. In addition, the researchers reported that the ZrO2 coating prevented LSC nanoparticles from agglomerating and suppressed Sr enrichment at the surface during long-term operation. As a result, the nanostructured LSC–ZrO2 electrode exhibited significantly lower polarization resistance and a 19 times slower degradation rate over 4000 h at 700 °C in air as compared with pristine LSC. In another study, Gong et al.240 reported that ZrO2 coatings can reduce the degradation rate of porous LSCF–GDC to ∼25% as compared with undecorated samples during 1100 hours of operation at 800 °C. Saher et al.242 also reported that LSCF skeletons with a porous nanoparticulate GDC layer can exhibit a high surface exchange coefficient (kchem) of ∼3 × 10−4 m s−1 at 900 °C, which was 6 times higher than that of pure LSCF. Here, the researchers suggested that the decoration of GDC particles facilitated the dissociative adsorption of O2, leading to enhanced activities for the heterostructure.
Dispersed MgO,114 NiO115 and CuO116 particles have also been highlighted as promising synergistic catalysts of LSCF electrodes. For example, Yang et al.114 reported that by infiltrating MgO particles onto LSCF skeletons, the ASR of the oxygen electrode can be reduced from 0.49 Ω cm2 to 0.31 Ω cm2 at 650 °C. In addition, researchers reported that the introduction of CuO can greatly improve the surface exchange kinetics of LSCF with kchem rising from 2.6 × 10−5 cm s−1 to 9.3 × 10−5 cm s−1 at 750 °C.116 Furthermore, coatings on cathodes can also enhance the cell performance, in which cells with Ni–YSZ|YSZ/SDC|LSCF–NiO can exhibit a high peak power density of 1.031 W cm−2 at 800 °C, which is ∼1.5 times higher than that of single cells with bare LSCF as the cathode.115
Despite these promising results, however, there are also studies highlighting the fact that surface decorations with simple metal oxides may reduce the electrochemical performance. For example, Kim et al.243 modified La0.6Sr0.4CoO3-based electrodes with a Al2O3 layer using ALD and reported that excessive Al2O3 coating on LSC-based electrodes can negatively affect the power output (reduced from 1.62 to 0.742 W cm2 at 800 °C). In addition, there are also controversial evaluations concerning CoOx and SrO coatings,254,255 in which the results suggest that the surface modification process needs to be rationally controlled. Here, proper amounts of decoration can activate electrode surfaces whereas excessive coating may further block active sites and degrade the performance.
Furthermore, Chen et al.256 demonstrated that a hybrid catalyst coating consisting of a conformal PrNi0.5Mn0.5O3 (PNM) thin film and exsoluted PrOx nanoparticles on LSCF obtained using a one-step infiltration process can significantly improve the ORR kinetics and durability (Fig. 12c). The heterostructured cathode reduced the polarization resistance to ∼0.022 Ω cm2 at 750 °C, which was 1/6 that of bare LSCF cathodes. Here, the researchers reported that in full cell operation, anode-supported cells with the hybrid catalyst-coated LSCF cathode exhibited an extraordinary peak power density of ∼1.21 W cm−2 (Fig. 12d) and excellent durability at 0.7 V for 500 h (Fig. 12e). And if applied in BaZr0.1Ce0.7Y0.1Yb0.1O3-based protonic fuel cells, electrodes combining PNM with exsoluted PrOx particles can exhibit excellent electrocatalytic activity with ASRs as low as 0.052 Ω cm2 at 700 °C and operation stability under 0.7 V for ∼500 h.257 In another study, a PrBa0.5Sr0.5Co1.5Fe0.5O5+δ (PBSCF) nanofibre-based cathode with PrOx nanoparticles exhibited high activity for oxygen reduction. The PPD of PBSCF–PrOx|SDC|Ni–BZCYYb reached as high as 0.37 W cm−2 at 500 °C with nearly dry methane as the fuel and air as the oxidant.258
Researchers have proposed that tensile strain can generally increase the free spacing in crystal lattices and decrease bond energies and weaken interactions between ions,53 thus facilitating the formation of oxygen vacancies.231,267 For example, Petrie et al.231 reported that the amount of oxygen deficiency in strained SrCoO3−x can increase with x changing from 0 to 0.25 as the tensile strain increases from +1% to +4.2% (Fig. 13a). Here, the researchers used DFT and proposed that tensile strain can decrease the electronic hybridization between Co 3d and O 2p and further weaken Co–O bonds, leading to a greater intercalation formation enthalpy (Hi) and more thermodynamic oxygen instability in the tensile strained SrCoOx. The researchers also reported that a 2% tensile strain can reduce the activation energy barrier (Ea) by ∼30% and facilitate the process of oxygen diffusion. Such an increase in Hi and decrease in Ea will both result in an oxygen deficient state and further influence the electrochemical properties.
Fig. 13 (a) Oxygen deficiency in tensile strained P-SCO determined by XRD and O-K and Co-L2 edge XAS measurements.231 (b) Vacancy formation energies of LSM as a function of strain using an oxygen chemical potential for 1173 K and 0.1 atm pO2.267 (c) Comparison of calculated Ovacancy formation energies for LSC/STO and LSC/LAO with different oxygen vacancy concentrations. The insets exhibit the smaller oxygen vacancy formation energies for LSC/STO.268 (d) Depth profiles of 18O− distributions for tensile strained (100) LSC82 thin films on STO and compressive strained LSC82 on LAO after oxygen isotope exchange at 400 °C for 5 min. Reproduced from ref. 228, 231, 267 and 268 with permission from WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim, copyright 2016; Elsevier, copyright 2017; WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim, copyright 2019 and American Chemical Society, copyright 2013 (https://doi.org/10.1021/nn305987x), respectively. |
In another example, Morgan et al.267 used ab initio calculations to develop a quantitative model for the evaluation of strain impacts on vacancy formation energies (Fig. 13b) and proposed that tensile strain can increase lattice volumes, which caused decreased vacancy formation energies. By systemically studying strained LaMnO3, La0.75Sr0.25MnO3, LaFeO3 and La0.75Sr0.25FeO3, the researchers inferred that the vacancy formation energies for most oxides generally declined under tensile strain; a 1% increase in the biaxial strain was estimated to correspond to a ∼30–100 meV decrease in the vacancy formation energy.
Recently, Chen et al.268,269 investigated the critical impact of lattice strain on the oxygen defect chemistry and OER activity of LSC in liquid solutions through a combination of DFT calculations and experimental measurements. Here, the DFT results showed that the oxygen vacancy formation energy of tensile LSC/STO was smaller than that of LSC/LAO (Fig. 13c). As a result, LSC/STO contained more lattice oxygen vacancies than LSC/LAO. These calculations were also consistent with the larger δ values observed in HRXRD and the increased Co2+ signals in XPS characterization of LSC/STO.
Han et al.51 suggested that strain effects can influence not only the formation of oxygen vacancies, but also the incorporation and transportation of oxygen defects. By using DFT+U calculations, the researchers proved that lattice strain near the LSC113/LSC214 heterointerfacial region can enhance the concentration and mobility of oxygen vacancies, leading to highly accelerated oxygen dissociation and incorporation kinetics.265,270 Here, the researchers reported that the tensile strain in LSC113(001) (+1.9%) reduced elastic interactions along the oxygen migration paths and facilitated both oxygen adsorption and oxygen vacancy formation in neighboring regions. In addition, the researchers also reported that the tensile strain on LSC113 in the (110) direction can lead to more available spacing among crystal lattices for oxygen migration, minimizing energy barriers. Furthermore, calculations in this study demonstrated that lattice strain can contribute to 3 × 102 times of conductivity enhancement at 500 °C (considered both in LSC113 and LSC214 together).
There were also some studies indicating that strain effects on electronic structures can greatly influence the charge transfer kinetics in ORR/OER processes. For example, Cai et al.259 studied the dependence of surface electronic structures on strain states in epitaxial La1−xSrxCoO3−δ (LSC) thin films with STO and LAO substrates. Results from in situ STM/STS revealed an electronic structural transformation from a semiconducting state with an energy gap of 0.8–1.5 eV at room temperature to a metallic-like state with no energy gaps at 200–300 °C. The researchers also found that in tensile strained LSC, the electronic density around the Fermi level was strengthened following this transition, facilitating charge transfer in the ORR. And by collaborating with Kubicek et al.,228 these researchers further demonstrated that the surface exchange coefficient (k*) and diffusion coefficient (D*) of the tensile strain of LSC–STO were ∼4 and ∼10 times larger than that of the compressive strain of LSC–LAO, respectively (Fig. 13d).
In a recent review article, Hwang et al.271 discussed the impacts of lattice strain on the electronic structure and electrocatalytic activity of perovskite oxides. They proposed that lattice strain can change the level of octahedral rotation and structural deviation as well as the bond length of metal-oxides, thus modifying the oxide electronic structures through 3d orbital occupancy and M 3d–O 2p orbital overlap.272,273 As the energy gap between M 3d and O 2p centers was considered as a quantitative descriptor of the electronic changes, the lattice strain may therefore have broad effects on not only the surface reactivity such as the oxygen adsorption and dissociation energy, but also the bulk defect energetics such as the oxygen vacancy formation and migration energy.
Zhang et al.249 used TEM/EDAX to observe Sr, Fe and Co signal peaks on LNO nanoparticles in LNO-infiltrated LSCF cathodes after long term testing, revealing the occurrence of the cation segregation phenomenon (Fig. 14a). Here, the researchers reported that the segregation of highly-active cobaltic oxides accelerated oxygen surface exchange and improved the electrocatalytic activity. Furthermore, the Sr diffusion increased the concentration of electron holes and enhanced the charge transport process. These resulted in a 67% increase in power output over long-term operation.
Fig. 14 (a) TEM image and EDAX spectra of LNO-infiltrated LSCF particles after 1000 h operation.249 (b) HADDF-STEM micrograph of the NSC113/NSC214 interface of the NSC214/113-17L sample and ratios of Sr/(Sr + Nd) and (Sr + Nd)/Co along the (001) direction.107 Reproduced from ref. 107 and 249 with permission from Elsevier, copyright 2014 and 2018, respectively. |
Researchers have attempted to uncover the driving forces of cation segregation at heterointerfaces. For example, Zheng et al.107 used HADDF-STEM to obtain direct evidence of cation segregation in the interfacial region of NSC113/214. According to their results, the stacking sequence of atoms along the c-axis was in the order of perovskite (Nd0.5Sr0.5CoO3−δ) and Ruddlesden–Popper (Nd0.8Sr1.2CoO4±δ), and the (Sr + Nd)/Co ratio conformed to the stoichiometry, respectively. However, the Sr/(Sr + Nd) ratio was found to be enhanced obviously at the interfacial region, revealing that Sr strongly segregated near the NSC113/214 interface (Fig. 14b). The researchers also used DFT calculations to reveal that E(SrNd) (the energy required to replace a surface Nd3+ cation with a lattice Sr2+ cation) decreased significantly at the interfacial region (−0.53 eV), indicating a strong thermodynamic trend that can drive Sr from the perovskite bulk to the interface.
Gadre et al.223 used DFT simulations to calculate cation diffusion at the LSC113/LSC214 heterointerface and reported that Sr tended to migrate to the interfacial region thermodynamically. Moreover, the heterointerface can stabilize high dopant Sr concentrations at the interface and suppress the formation of secondary passive particles (such as SrO). In addition, the researchers reported that the increased concentration of Sr can induce the formation of oxygen vacancies, leading to enhanced electrocatalytic activities. In addition, combining COBRA with ab initio calculations, Feng et al.57 provided atomic-scale structural information on the LSC113/LSC214 heterointerface, in which the signals at (0,0,0) and (0.5,0.5,0) (Fig. 15a and b) exhibited apparent Sr segregation at the interface and near the surface of LSC214. Here, the low SrLa substitution energy in the heterogeneous region was believed to be responsible for Sr segregation. Using ab initio calculations, Lee et al.230 found that the SrLa substitution energy in the LSCF113/LSC214 heterogeneous region was much lower than that of LSC113/214 (Fig. 15c), which was consistent with the weaker Sr interdiffusion trend of LSCF113/LSC214 as observed by AES.278
Fig. 15 Electron densities of as-deposited LSCO214/LSCO113/STO thin films along: (a) the (0 0 Z) and (b) (0.5, 0.5, Z) lines using COBRA.57 (c) Up: SrLa substitution energies in bulk LSC214, LSC113 and LSCF113 calculated by DFT simulations. Middle: heterostructure model used in DFT simulations. La/Sr (dark blue), Fe/Co (light blue, center of the octahedra) and O (red). Bottom: SrLa substitution energies of each layer near the heterointerfacial region.230 Reproduced from ref. 57 and 230 with permission from American Chemical Society, copyright 2014 (https://pubs.acs.org/doi/10.1021/jz500293d) and RSC, copyright 2015, respectively. |
In addition, Tsvetkov et al.236 proposed that the phase separation of Sr-rich particles was detrimental to the ORR kinetics and that the decoration of less reducible cations (, Me = V5+, Nb5+, Ti4+, Zr4+, Hf4+ and Al3+) onto perovskite LSC can balance Sr segregation with the ORR kinetics. Here, the researchers suggested that decoration can reduce the oxygen vacancy concentration acting as reactive sites, while the coatings can increase the Co oxidation state and decrease the electrostatic attraction of Sr enrichment onto the surface, further suppressing harmful phase separation and stabilizing the surface composition. Based on the results of APXPS/XAS, the researchers summarized that these two factors resulted in a ‘volcano’ relationship between kq and ΔHVf(MeOx). In another study, Rupp et al.253 reported that SrO decoration onto LSC severely deactivated the LSC thin film, whereas deposited Co3O4 re-activated the surface. Here, the researchers used IPLD and found that deposited Sr cations can migrate and attach to highly active sites of inhomogeneous surfaces, hindering ionic and electronic transport in the ORR process,279 leading to the deactivation of active sites. Alternatively, Co3O4 decoration can induce highly active sites on surfaces.280
Surface Sr segregation and separation are major causes of the degradation of Sr containing electrodes. It has been demonstrated that the deposition of a secondary phase to form a heterostructure can effectively suppress the phase separation of Sr. For example, Wen et al.281 proposed that ZrO2 coatings fabricated using ALD can suppress excess Sr segregation in LSC epitaxial films and improve the stability of oxygen electrodes. Here, the researchers suggested that the existence of cation exchange between Zr and Co near the interface can reduce the electrostatic force of surface Sr segregation, thus improving the thermal and electrochemical stability. In another study, Gong et al.241 confirmed that ZrO2 overcoats can optimize nanostructured LSC surfaces by preventing agglomeration and suppressing Sr segregation in long term operation. Similarly, the suppression of surface Sr segregation using surface coatings was also reported for LSM-decorated LSCF,41 Pr0.2Ce0.8O2/CeO2-decorated LSCF,234,282 PSM/PSCM-decorated LSCF,247,283etc.
By combining focused ion beam (FIB) milling235 with in situ scanning tunneling microscopy and spectroscopy (STM/STS), Chen et al.56 discovered that LSC214 in a LSC113/214 multilayer (ML) can be electronically activated. This electronic activation was highlighted by the smaller bandgap and a reversible transition from a semiconductor state to a metallic-like state with the temperature increasing to ∼250 °C (Fig. 16a). Here, the researchers reported that in the LSC113/214 interface, electron injection from LSC113 to LSC214 can result in band alignment and a decreased distance between the Fermi level and the conducting band of LSC214. The reduced distance can exponentially increase the quantity of available electrons and accelerate the kinetics of electron transfer from the electrode surface to adsorptive oxygen (Fig. 16b). Furthermore, Ma et al.77 reported that VAN-structured LSC113/214 exhibited similar electronic structure changes to ML LSC113/214. The energy bandgap of VAN was 1.5 ± 0.2 eV at room temperature, while metallic-like behavior with no energy gap was found above 250 °C. Here, the researchers proposed that electronic activation and stable surface composition can lead to enhanced ORR activity. In a recent study, Chen et al.80 used HRXRD (high resolution X-ray diffraction) and HAXPES (hard X-ray photoelectron spectroscopy) to probe the distribution of oxygen defects across the heterointerfaces of LSC113/LSC214 and LSC113/LNO214. Here, larger concentrations of oxygen vacancies were found in the LSC113 of LSC113/LSC214 than in the LSC113 of LSC113/LNO214 in a reducing environment, while both the LSC214 and LNO214 layers in the heterostructured samples were found to be more reduced compared with the respective single counterparts (Fig. 16c).
Fig. 16 (a) Transformation of apparent energy gaps in single phase LSC113 and LSC214 as well as multilayered LSC113/LSC214 at different temperatures.56 (b) Origin of the enhanced ORR activity at heterostructured LSC113/LSC214 layers.56 (c) Illustrations of how the LSC113/LSC214 and LSC113/LNO214 interfaces induce the formation of local electric fields and influence local oxygen defect equilibria; (left) the misalignment of Fermi levels in LSC113 and LSC214 (LNO214) single phases after the loss of oxygen through annealing at high temperatures, (right) the alignment of Fermi levels in LSC113 and LSC214 (LNO214) heterostructures after annealing through the creation of local electric fields.80 Reproduced from ref. 56 and 80 with permission from WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim, copyright 2013 and American Chemical Society, copyright 2018, respectively. |
Fig. 17 (a) Illustration of an idealized model of an integrated miniature half-cell with sufficient TPB; (b) TEM image and schematic of a LSC113 thin film with LSC214 particle decoration.48 Reproduced from ref. 48 with permission from American Chemical Society, copyright 2010. |
Crumlin et al.48 reported high ORR activities on the heterointerface of nano-island LSC214 (La0.5Sr0.5)2CoO4±δ decorated LSC113 (La0.8Sr0.2CoO3−δ) films, in which they proposed that the dispersive nano-island LSC214 can extend the TPB, leading to excellent ORR activity (Fig. 17b). Wan et al.288 also proposed that the oxygen reduction kinetics of PrBaCo2O5+δ–Sm0.2Ce0.8O1.9 (PBC–SDC) composites can be enhanced through the introduction of 3PB (PBC–SDC-gas), in which the polarization resistances of the oxygen electrode decreased from 0.281 Ω cm2 (for bare PBC) to 0.107 Ω cm2 (for the composite with 60 wt% SDC) at 700 °C. Here, the researchers suggested that the 3PB accounted for more than 80% of active sites in the oxygen incorporation processes.
Overall, the development of high-performance heterointerfaces is not only a promising strategy to enhance the electrode performance but can also greatly extend the scope of selection for new electrode materials. Therefore, recent progress in the utilization of heterostructured oxides as high performance fuel electrodes will be summarized in this section to provide direction for the research of more efficient and economical fuel electrode materials to further improve the SOC performance and stability.
In one example, Yang et al.293 synthesized a novel BaO–Ni/YSZ anode with a finely distributed Ni/BaO heterointerface using vapor deposition and subsequent reduction and tested a single cell (BaO–Ni/YSZ|YSZ/SDC|LSCF) in CO. Here the researchers achieved a typical PPD of ∼0.70 W cm−2 at 750 °C, which was more than twice that of conventional Ni–YSZ anodes.301,302 Furthermore, the cell in this study provided stable power outputs for ∼100 h with 500 mA cm−2 in wet CO, demonstrating superior coking resistance. In another example, Tao et al.60 found that the BaO/Ni interface can stabilize cathode surfaces against coking in hydrocarbon fuels, in which the peak power density of their oxide cell (Ni–BaZr0.1Ce0.7Y0.2O3−δ (BZCY)–SDC|SDC|LSCF–SDC) reached 570 mA cm−2 at 700 °C in methane in 150 minutes of operation with a low electrode resistance of 0.136 Ω cm2. Likewise, Islam et al.303 reported that 5 wt% BaO loading on Ni/YSZ through infiltration increased the peak power density by ∼40% at 1073 K in 100% CH4.
Researchers have studied the mechanisms of coking resistance properties provided by extrinsic metal oxide decoration.59,60,293 For example, Yang et al.293 proposed that under operating conditions, carbon can be reversibly generated and removed at Ni surfaces through 5 primitive steps, including: (1) CO adsorption onto Ni due to a high adsorption energy (−2.6 eV) and the water-gas shift reaction; (2) carbon deposition on the surface of Ni catalysts due to the dehydrogenation of hydrocarbon fuels or CO disproportionation; (3) strong H2O adsorption on the BaO sites of BaO/Ni(111), leading to barrierless O–H bond cleavage surrounding BaO; (4) the reaction between dissociated OH and adsorbed C on Ni surfaces to form an intermediate COH, which dissociates to CO and H; and (5) CO and H diffusion to the TPB and electrochemical oxidation to CO2 and H2O.
Alternatively, Li et al.59 have studied the origin of enhanced coking resistance through BaO decoration on Ni surfaces using in situ Raman spectroscopy. They found the formation of a thin Ba(OH)2 layer on Ni surfaces during operation in wet hydrocarbon fuels, which further enabled the conversion of deposited coke to CO or CO2 rapidly. Here, the possible chemical reactions occurring near the interface can be summarized as follows:
BaO + H2O ⇔ Ba(OH)2; Cad + Ba(OH)2 → BaO + CO + H2; CO + Ba(OH)2 → BaO + CO2 + H2; BaO + CO2 → BaCO3.
This mechanism was supported by the rapid decrease in the –OH intensity accompanied by the gradually increasing –CO3 signal as measured by in situ Raman spectroscopy (Fig. 19).59 Likewise, Tao et al.60 proposed that the BaO/Ni interface can accelerate water adsorption and enhance the reaction of water-mediated carbon removal, whereas fast oxygen ion transportation and generated steam at the anode can promote hydrocarbon reformation to enhance the performance and stability.
Fig. 19 (a) Integrated peak intensities of the carbon G-bands for bare Ni and BaO–Ni as functions of exposure time in wet C3H8 at 450 °C. (b) Collection of time-resolved Raman spectra on Ni powders and (c) on BaO modified Ni powders. Reproduced from ref. 59 with permission from American Chemical Society, copyright 2015. |
Despite enhanced surface stability, however, surface decoration needs to be precisely controlled to avoid negative effects such as the formation of insulating alkaline oxides that can block active sites. For instance, the covering of TPB active regions with BaO particles can increase electrode polarization resistances and partly lower power densities.59,294,304,305 To address this, Islam et al.303 recently developed a novel fabrication technique combining microwave irradiation with vapor deposition to control the distribution and quantity of deposited nanosized BaO particles to ensure that additive components were only deposited separately on Ni surfaces.
In addition to BaO, decoration with alkaline earth metal oxides such as MgO and CaO has also been reported to enhance the electrode stability through similar mechanisms. For example, Qu et al.294 modified Ni–SDC anodes using CaO nanoparticles and achieved only slight degradations in performance during operation in a wet CH4 atmosphere for 70 h at 650 °C whereas an unmodified Ni–SDC-based cell failed to generate power within 1.2 h in the same conditions. Hua et al.299 also infiltrated Ni–YSZ cermet anodes with MgO additives and reported an Rp of 0.93 Ω cm2 after 48 h operation, which was 31% lower than that of a NiO infiltrated anode. Furthermore, Lay et al.300 proposed that a moderate amount of MgO decoration on the Ni–YSZ surface can lead to accelerated steam gasification reactions and enable less generation of surface carbon fibers in CH4 atmospheres.
Vincent et al.295 fabricated a YSZ/La0.4Sr0.5Ba0.1TiO3 anode through infiltration and achieved a peak power density in the single cell of 200 mW cm−2 at 850 °C in H2/H2S (5000 ppm). Moreover, Xu et al.296 fabricated a CeO2/La0.3Sr0.7Ti0.3Fe0.7O3−δ composite anode and reported polarization resistances as low as 0.151 Ω cm2 at 850 °C in a CO atmosphere. More recently, Chen et al.258 designed a Ni–BaZr0.1Ce0.7Y0.1Yb0.1O3−δ (BZCYYb)/SDC hybrid SOFC anode (Fig. 20) and reported that the corresponding cell was able to operate in nearly dry CH4 fuel (with ∼3.5% H2O) for more than 550 h at only 500 °C without coke formation.
Fig. 20 (a) Schematic of a single cell with yellow and grey grains representing BZCYYb and Ni phases. (b) SEM of the anode after long term operation for ∼500 h in nearly dry CH4, indicating excellent coking tolerance. Reproduced from ref. 258 with permission from Springer Nature, copyright 2018. |
In another example, Hua et al.307 reported that decorating Ni/YSZ anodes with PrBaMn2O5+δ nanoparticles can significantly enhance the catalytic activity and durability, leading to a 1.35 W cm−2 peak power density in a biogas atmosphere (containing CH4, CO2 and small amounts of H2S) at 800 °C as well as negligible performance degradation after 100 h of testing. Here, the researchers suggested that through infiltrating the Ni/YSZ surfaces with perovskite oxides, the finely distributed perovskites with high catalytic activity can act as micro-anodes that can increase the TPB area largely and enhance the processing capacity of the anode for converting hydrocarbon fuels (Fig. 21). Furthermore, Irvine et al.4 suggested that the heterogeneous and nanosized decoration of perovskite oxides such as SrTiO3 and LaCrO3 on fuel electrode surfaces can improve the long-term stability, resulting from the superior stability over a wide range of environments and moderate electronic conductivities. This mechanism appears to also be applicable in perovskite-modified Ni/YSZ electrode systems as well.
Fig. 21 Comparison of Ni/YSZ anodes before (a) and after (b) PrBaMn2O5+δ infiltration. The infiltrated perovskite nanoparticles resisted carbon species deposition on the Ni/YSZ surfaces and enhanced carbon removal. Reproduced from ref. 307 with permission from Elsevier, copyright 2017. |
Self-assembled metallic nanoparticles can normally exhibit higher conductivities and catalytic activities than perovskite bulk phases, leading to significant enhancements in fuel electrode electrochemical activity.4,309–312 For example, Yang et al.311 proposed that exsolved NiO particles on SrFe0.85Ti0.1Ni0.05O3−δ (STFNi) electrode surfaces can significantly reduce the ASR from ∼0.4 to ∼0.15 Ω cm2 at 700 °C, in which the anode-supported cell with NiO–STFNi as the fuel electrode can exhibit high peak power densities of 1650, 1420 and 1100 mW cm−2 at 750, 700 and 650 °C respectively. Tsekouras et al.312 further compared La0.4Sr0.4Ni0.06Ti0.94O2.94 (LSNT) with Ni0 exsolution, La0.4Sr0.4Fe0.06Ti0.94O2.97 (LSFT) with Fe0 exsolution and La0.4Sr0.4TiO3 (LST) and found that the exsolution trend of Ni0 was generally larger than Fe0 exsolution. Here, electrolysis cells operating in 47% H2O/53% N2 at 900 °C at −1.5 V with the above cathodes showed reduced HTSE onset potentials (EHTSE) from −1.19 V for the LST cathode to 0.98 V and 0.63 V for the Fe0 and Ni0 exsoluted cathodes, respectively, with Rp decreasing from 1.7 to 1.5 and 1.2 Ω cm2. Likewise, Li et al.309 fabricated a perovskite-structured (La0.2Sr0.8)(Ti0.9Mn0.1)O3−δ (LSTM) skeleton with exsolved Ni nanocrystals on the surface and reported a ∼50% enhancement in catalytic activity for CO2 electrolysis as compared with a bare electrode. The high catalytic activity for CO2 electrolysis is likely associated with the high oxygen vacancy concentration.298 The exsolution of metal particles is often accompanied by formation and aggregation of oxygen vacancies near the interface. Theoretical calculations indicated that the exsolved Ni clusters and oxygen vacancies could synergistically weaken the C–O1 bond, thus activating the adsorbed CO2 molecule and enhancing the CO2 reduction. In addition, the increase in the oxygen vacancy concentration enhances the ionic conductivity of the fuel electrode through a vacancy-mediated mechanism.
Many studies have also confirmed that the surface composition and electrochemical performance of these exsolved particle systems experience no significant decline after many redox cycles, demonstrating considerable operational stability.7,314–317 For example, Liu et al.318 used a double-layered perovskite ((Pr0.4Sr0.6)3(Fe0.85Mo0.15)2O7) with exsolved Co/Fe metallic nanoparticles and found that no graphite was generated during polarization for 100 h of CO2 electrolysis, indicating excellent anti-coking properties. Furthermore, Song et al.319 fabricated a self-assembled dual-phase structure (R-LS0.6FCN) consisting of both SP and RP (SLF214 + LFO113) layers with exsolved Fe–Cu bimetal nanofibers and applied it as the fuel electrode of a SOFC. Here, the researchers reported that the peak power density of the solid oxide cell using R-LS0.6FCN as the anode yielded 1.03 W cm−2 at 850 °C in H2/3% H2O, which was higher than that of a similar single cell with Ni–Fe/SDC cermet as the anode. These researchers also reported that there was no obvious performance degradation after operation in H2/50 ppm H2S or syngas at 700 °C for 100 h, demonstrating high durability against coking/sulfur deposition. Here, such high stabilities were attributed to the formation of nanofiber networks without agglomeration and to the construction of an SP/RP heterointerface at the anode. More studies about the coking resistance capacity of perovskite-based fuel electrodes with exsolved metals are summarized in Table 3.
Exsolved metal | Perovskite precursor | Temperature (°C) | Atmosphere | Cell configuration | Performance | Ref. |
---|---|---|---|---|---|---|
Fe | La0.6Sr0.3Cr0.85Fe0.15O3−δ (LSCFe) | 800 | CH4 | LSCFe|YSZ–YSZ–LSM|YSZ | PPD ∼ 330 mW cm−2 | 320 |
Co | Pr0.5Ba0.5Mn0.9Co0.1O2.89 (PBMCo) | 900 | Syngas | PBMCo|GDC–GDC–YSZ–GDC–LSM|GDC | PPD ∼ 900 mW cm−2 | 321 |
Ni | La0.52Sr0.28Ni0.06Ti0.94O3 (LSNT) | 800 | 20% CH4/H2 | Slight carbon fiber growth after 4 h exposure | 313 | |
La0.9Mn0.8Ni0.2O3 (LMN) | 700 | CH4/CO2 | Negligible change for 24 h | 322 | ||
Cu | La0.5Sr0.5Fe0.8Cu0.15Nb0.05O3−δ (LSFCNb) | 800 | Wet syngas | LSFCNb–SDC–ScSZ–SDC–LSFCNb | PPD ∼ 640 mW cm−2; 9.4% degradation rate for 100 h at 0.7 V | 323 |
Co–Fe alloys | (Pr0.4Sr0.6)3(Fe0.85Nb0.15)2O7 (RP-PSFN) | 800 | Wet C3H8 | RP-PSFN–CFA|LSGM|BSCF | Stable operation for 100 h under 400 mA cm−2 | 324 |
Pr0.8Sr1.2(Co,Fe)0.8Nb0.2O4+δ (K-PSCFN) | 850 | C3H8 | K-PSCFN–CFA|LSGM|P-PSCFN | PPD ∼ 940 mW cm−2; negligible degradation during 26 cycles of redox testing | 325 | |
La0.5Ba0.5Mn0.8Fe0.1Co0.1O3−δ | 850 | Wet CH4 | CoFe–LBMF|LDC–LSGM|BSCF | PPD ∼ 503 mW cm−2; stable power output for 200 h | 326 | |
La0.3Sr0.7Cr0.3Fe0.6Co0.1O3−δ (LSCrFeCo10) | 850 | C3H8/H2O | LSCrFeCo10–GDC|LDC–LSGM|LSCF–GDC | PPD ∼ 700 mW cm−2 | 327 | |
(Pr0.4Sr0.6)0.3(Fe0.85Mo0.15)2O7 (DLP-PSFM) | 750 | C2H6 | Stable operation for 100 h under 650 mA cm−2 | 328 | ||
Fe–Ni alloys | Sr2FeMo0.65Ni0.35O6−δ (SFMNi) | 850 | Wet CH4 | SFMNi–LDC|LSGM|LSCF | PPD ∼ 500 mW cm−2, 24 h of stability under 200 mA cm−2 | 329 |
(La0.7Sr0.3)(Cr0.85Ni0.1125Fe0.0375)O3−δ (LSCNi-Fe) | 850 | 5000 ppm H2S-syngas | LSCNi-Fe|YSZ|YSZ–LSM/YSZ | PPD ∼ 400 mW cm−2 and a 1300 mA cm−2 maximum current density | 330 | |
La0.5Sr0.5Fe0.8Ni0.2O3−δ (LSFN) | 725 | Wet syngas | LSFN–SDC|ScSZ|SDC–LSFCNb | Stable operation for 20 h at 0.7 V | 331 | |
Fe–Pb alloys | La0.8Sr0.2Fe0.9Nb0.1Pd0.04O3−δ (LSFNP) | 800 | Dry C3H8 | LSFNP|SDC|NBCCo–SDC | 25 h operation under 80 mA cm−2 | 332 |
Fe–Re alloys | LaRexFe0.8O-La3ReO8 (x < 0.2) (LRF) | 860 | CH4/CO2 | No significant carbon deposition after 70 h operation | 333 | |
Fe–Cu alloys | (La0.75Sr0.25)0.9(Cr0.5Mn0.5)0.9(Cu1−xFex)0.1O3−δ | 850 | Dry CH4 | Cu0.5Fe0.5–LSCM|LSGM |LSM | PPD ∼ 40 mW cm−2, negligible degradation after 100 h operation and 10 redox cycles | 334 |
Ni–Ce alloys | La0.3Sr0.6Ce0.1Ni0.1Ti0.9O3−δ (LSCNT) | 900 | CH4 | LSCNT–YSZ|YSZ|LSM–YSZ | PPD ∼ 600 mW cm−2, no degradation for 80 h under 0.5 V | 49 |
Ni–Cu alloys | (La0.75Sr0.25)0.9(Cr0.5Mn0.5)0.9(Ni1−xCux)0.1O3−δ (LSCM-Ni1−xCux) | 800 | CH4/CO2 | Symmetric cell | 300 h stability at 0.9 V | 335 |
Alternatively, Neagu et al.313 proposed that the high stability of perovskite–metal heterosystems resulted from a combination of size effects and formation of robust interfaces, in which for a specific La0.52Sr0.28Ni0.06Ti0.94O3 system with in situ exsolved Ni metal (Fig. 22a), after operation in a 20% CH4/H2 environment at 800 °C for 4 h, there was nearly no carbon fiber on the exsolved Ni particles with sizes of 25–60 nm, while carbon growth was found on exsolved Ni particles of ∼80 nm in size. For comparison, there was obvious carbon fiber accumulation on La0.4Sr0.4TiO3 anodes infiltrated by 20 nm Ni particles. TEM micrographs indicated that exsolved nanoparticles were embedded at a considerable depth in the parent perovskite oxide support (Fig. 22b), resulting in a robust interface between the perovskite oxide and the exsolved metal particles. The TEM micrographs and electron diffraction revealed that the exsolved metallic particles could grow epitaxially from the parent perovskite oxide in a specific orientation (perovskite(111), Ni(111)), forming a well-ordered structure across the interface with good stability (Fig. 22c and d). In addition, the cation diffusion across the interface between the metal particles and the parent oxide greatly increases the adhesion of the metal particles to the oxide phases. The strong bonding between the exsolved metal particles and the parent perovskite oxide plays a vital role in preventing particle uplifting and subsequent carbon fiber growth, resulting in high coking resistance and stable catalytic activity in a carbon-containing atmosphere.
Fig. 22 (a) Schematic of the possible mechanism of B-site exsolution in a perovskite-structured La0.52Sr0.28Ni0.06Ti0.94O3 model system. (b) HRTEM image of the interfacial region between the perovskite matrix and exsolved Ni nanoparticles. (c) HRTEM micrographs at the perovskite–metal interface, demonstrating the epitaxic atomic planes and orientations. (d) Atomic arrangement at the Ni–La0.52Sr0.28Ni0.06Ti0.94O3 interface according to structure characterization. Reproduced with permission from ref. 313, copyright 2015 Macmillan Publishers Limited. |
The driving force for B-site cation exsolution was also investigated by several groups.4,310,336–339 According to Irvine et al.,4,310,313 the exsolution of B-site cations is dynamically driven through a combination of chemical and electrochemical reduction and generally initiates at surface defect sites. Specifically, in lattices experiencing reducing conditions, the reduction of B-site transition metals and the formation of surface defects (i.e. oxygen vacancies) become more favorable, collectively leading to metal nucleation at the surface and the lowering of nucleation barriers through crystal defects.30 More recently, Han et al.336 successfully controlled B-site exsolution in La0.2Sr0.7Ni0.1Ti0.9O3−δ (LSNT) thin films through lattice strain control, resulting in a high degree of Ni exsolution with more than 1100 particles μm−2 and tiny particle sizes of ∼5 nm in diameter. Here, thermodynamic analysis indicated that the small atomic diffusion coefficients of the exsoluted Ni under both tensile and compressive strain were much lower compared to conventionally deposited Ni, indicating a small thermodynamic potential of particle coarsening and agglomeration. Furthermore, the researchers suggested that due to the small size of the exsoluted Ni and strong interactions between the embedded Ni and perovskites, carbon fibers grew in “base growth” without particle growth. The key factor was deduced to be the mismatch-strain relaxation energy (ΔGmr), which can reduce nuclei sizes and nucleation barriers and simultaneously quicken nucleation rates, leading to smaller exsolution particles with higher densities.
Alternatively, Kwon et al.337,338 explained B-site exsolution in terms of thermodynamics, in which, based on density functional simulations, the researchers proposed that in a PrBa(Mn,T)2O5+δ (T = Fe, Co and/or Ni) double-perovskite model system, substituted cations such as Co2+ and Ni2+ possess lower exsolution energies than the host cation Mn2+ (−0.55 and −0.50 eV vs. −0.47 eV) and demonstrate more favorable thermodynamic trends. In addition, Co/Ni doped PrBa(Mn,T)2O5+δ exhibits a lower surface oxygen vacancy-formation energy than undoped samples (2.46 and 2.85 eV vs. 2.97 eV), which can also assist in stabilizing the self-grown metallic phase.
First, precise control of the composition and microstructure of heterointerfaces is necessary for the rational design of novel electrode or electrolyte materials via interface engineering. It has been reported that the atomic structure of the heterointerface such as local strain and dislocations can profoundly impact the ionic conductivity and/or the electrocatalytic activity of the heterostructure. The use of thin film techniques such as PLD, ALD and MBE presents unique advantages in controlling the atomic structure of the heterointerfaces, which makes them idea tools for synthesis of novel heterostructured SOC components. The challenges that need to be overcome are the high price of the equipment and the slow growth rates of these thin film techniques.
Secondly, advanced characterization techniques, particular the ones with in situ measurement capability, are needed to further reveal the role of heterointerfaces in determining the ionic conductivity and electrocatalytic activity of the heterostructure. Especially, it is necessary to probe how the interface evolves under the working conditions of SOCs, such as high temperature, a reactive gas environment and electrical bias. There have been several successful attempts to use in situ STM, XPS, and SIMS to investigate the evolution of the electronic structure, cation valence and ion distribution at the interfacial region at elevated temperature or under bias.11,56,233
Thirdly, although several factors have been proposed to explain the impact of heterointerfaces, an overall theoretical model that describes the role of heterointerfaces in determining the transport properties and catalytic activity of SOC components is still lacking. Advanced theoretical modelling and simulation need to be carried out to be combined with experimental results to uncover the effects of heterointerfaces, hence providing guidance for the rational design of novel materials.
Finally, besides improving the performance of each component of SOCs, system level optimization is also necessary. The compatibility of SOC materials with other parts of the system such as interconnection materials needs to be tested. The solid oxide cell and stack design and the thermal management of the whole system are also critical for the long stability of devices.
While SOC components are taken as examples in this review to introduce the influence of heterointerfaces on the conductivity, electrochemical activity and stability, the mechanism of the performance improvement and methodology for controllable fabrication of interfaces are applicable to other energy conversion systems for a wide range of applications such as photo/electrochemical water splitting, solar thermal CO2 reduction, and lithium ion batteries.
ALD | Atomic layer deposition |
ASR | Area surface resistance |
BSCF | Ba1−xSrxCo1−yFeyO3−δ |
BZCY | BaZr0.4Ce0.4Y0.2O3 |
BZO | BaZrO3 |
BZY | BaZrxY1−xO3−y |
CFA | Co–Fe bi-metallic alloy |
COBRA | Coherent Bragg rod analysis |
CSZ | ZrO2 + CaO |
DFT | Density functional theory |
DP-BSCFW | Double perovskite Ba0.5Sr0.5(Co0.7Fe0.3)0.687W0.312O3−δ |
EC | Electron conduction |
EELS | Electron energy loss spectroscopy |
EPD | Electrophoretic deposition |
HAXPES | Hard X-ray photoelectron spectroscopy |
HRXRD | High resolution X-ray diffraction |
HRTEM | High resolution transmission electron microscope |
HTPC | High-temperature proton conductor |
IC | Ionic conduction |
LAO | LaAlO3 |
LCP | LaCePrOx |
LNO | LaNiO3±δ |
LSC113 | La1−xSrxCoO3±δ |
LSC214 | La2−xSrxCoO4±δ |
LSCF | La1−xSrxCoyFe1−yO3±δ |
LSCM | La1−xSrxCoyMn1−yO3±δ |
LSCMC | (La0.75Sr0.25)0.9(Cr0.5Mn0.5)0.9Cu0.1O3−δ |
LSCMN | (La0.75Sr0.25)0.9(Cr0.5Mn0.5)0.9Ni0.1O3−δ |
LSGM | LaxSr1−xGayMg1−yO3−δ |
LSF | La1−xSrxFeO3±δ |
LSM | La1−xSrxMnO3±δ |
LSZG | Li13.9Sr0.1Zn(GeO4)4 |
MBE | Molecular beam epitaxy |
MIEC | Mixed ionic and electronic conductor |
ML | Multilayer |
NEXAFS | Near edge X-ray absorption fine structure |
NGO | NdGaO3 |
NSC113 | Nd1−xSrxCoO3±δ |
OER | Oxygen evolution reaction |
ORR | Oxygen reduction reaction |
PLD | Pulsed laser deposition |
PPD | Peak power density |
PSCM | PrSrCoMnO6−δ |
RP | Ruddlesden–Popper |
RP-PSFN | RP type layered (Pr0.4Sr0.6)3(Fe0.85Nb0.15)2O7 |
SCDC | Sm/Ca co-doped ceria |
SCFN | SrCo0.2Fe0.6Ni0.2O3 |
SDC | Samarium doped ceria |
SIMS | Secondary-ion mass spectrometry |
SOC | Solid oxide cell |
SOEC | Solid oxide electrolysis cell |
SOFC | Solid oxide fuel cell |
SP-BSCFW | Single perovskite Ba0.5Sr0.5(Co0.7Fe0.3)0.6875W0.3125O3−δ |
SSC–GDC | Sm0.5Sr0.5CoO3–Gd0.1Ce0.9O2−δ |
STEM | Scanning transmission electron microscope |
STF | SrTi0.2Fe0.8O3 |
STM/STS | Scanning tunnelling microscopy/spectroscopy |
STO | SrTiO3 |
TEM | Transmission electron microscope |
TOF-SIMS | Time of flight secondary ion mass spectroscopy |
TPB/3PB | Triple phase boundaries |
VAN | Vertical aligned nanocomposite |
XANES | X-ray absorption near-edge structure |
XAS | X-ray absorption spectroscopy |
XPS | X-ray photoemission spectroscopy |
YBCO | (Y,Ca)Ba2Cu3O7 |
YSZ | Yttira-stabilized zirconia |
Footnote |
† These two authors contributed equally to this article. |
This journal is © The Royal Society of Chemistry 2020 |