Highly efficient and robust cathode materials for low-temperature solid oxide fuel cells: PrBa0.5Sr0.5Co(2-x)Fe(x)O(5+δ).

Solid oxide fuel cells (SOFC) are the cleanest, most efficient, and cost-effective option for direct conversion to electricity of a wide variety of fuels. While significant progress has been made in anode materials with enhanced tolerance to coking and contaminant poisoning, cathodic polarization still contributes considerably to energy loss, more so at lower operating temperatures. Here we report a synergistic effect of co-doping in a cation-ordered double-perovskite material, PrBa0.5Sr0.5Co(2-x)Fe(x)O(5+δ), which has created pore channels that dramatically enhance oxygen ion diffusion and surface oxygen exchange while maintaining excellent compatibility and stability under operating conditions. Test cells based on these cathode materials demonstrate peak power densities ~2.2 W cm(-2) at 600°C, representing an important step toward commercially viable SOFC technologies.

The demand for clean and sustainable energy has stimulated great interest in fuel cells, which allows direct conversion of chemical fuels to electricity. Among all types of fuel cells, solid oxide fuel cells (SOFCs) have the potential to offer the highest energy efficiency and excellent fuel flexibility123456789. To make SOFC technology affordable, however, the operating temperature must be further reduced so that much less expensive materials may be used for other cell components and balance of plant5. Unfortunately, SOFC performance decreases rapidly as the operating temperature is reduced, especially the cathode for oxygen reduction reaction (ORR)1011. While La1−SrMnO3 (LSM) is widely used as the cathode material for yttria-stabilized zirconia (YSZ)-based SOFCs because of its excellent compatibility with YSZ electrolyte and other cell components; the cathodic polarization loss at lower temperatures is unacceptable. Accordingly, extensive efforts have been devoted to the search for more active cathode materials toward ORR at lower temperatures.

In particular, mixed conducting oxides with simple perovskite structure have been studied as alternative cathode materials for intermediate-temperature (IT)-SOFCs, including Sm0.5Sr0.5CoO3− (SSC)12, Ba0.5Sr0.5Co0.8Fe0.2O3− (BSCF)13, and La0.6Sr0.4Co1−FeO3−(LSCF)14. Among them, BSCF was reported to have very low area-specific resistance (ASR), about 0.05 ~ 0.08 Ω cm2 at 600°C. However, its applicability to commercial fuel cells is hindered by its limited compatibility with other cell components and long-term stability.

Recently, cation-ordered double-perovskite structures such as LnBaCo2O5+ (Ln = Pr, Nd, Sm, and Gd) have attracted much attention due to their easier oxygen ion diffusion, faster surface oxygen exchange, and higher electrical conductivity at lower temperatures than simple perovskite cathode materials151617, as corroborated by electrical conductivity relaxation (ECR)15 and ion exchange depth profile (IEDP)1718 studies. This family of double perovskite compounds have a general formula of AA′B2O5+, where A = Y or a trivalent lanthanide ion, A′ = Ba or Sr, and B is a first row transition metal ion or a mixture of them. These compounds have a layered structure with a stacking sequence of …[A′O] [BO2] [AO] [BO2]…, similar to the structure of the cuprate superconductors. The vacant sites are arranged so that the coordination number for A cation is eight when δ = 0. All oxygen vacancies are confined only to the LnO plane, and so is oxygen migration18. The high concentration of mobile oxygen species may be responsible for the high diffusivity of oxide ions in the bulk and the enhanced surface activity toward ORR15.

Some interesting effects of ion substitution on the properties of LnBaCo2O5+ have been reported192021. For example, the replacement of Pr by Gd was reported to diminish the concentration of oxygen defects, the electrical conductivity, and the performance as a cathode in a fuel cell19. On the other hand, a reverse trend was found for the substitution of Ba by Sr in the same material20. Further, various transition metal ions have been introduced into the B-site, including Fe, Cu, and Ni. As Co is substituted by Fe, the thermal expansion coefficient (TEC) and electronic conductivity decreased; however, the oxygen ion diffusivity, ORR activity, and stability increased. Thus, the properties can be tailored by the type and the amount of proper ion substitution.

Here we report a synergistic effect of co-doping (Sr on A-site and Fe on B-site) in a cation-ordered double-perovskite, LnBaCo2O5+, to create crystalline channels for fast oxygen ion diffusion and rapid surface oxygen exchange while maintaining the compatibility with the electrolytes for IT-SOFCs and the durability under operating conditions.

Results

Characterization of PrBa0.5Sr0.5Co2−FeO5+ (x = 0, 0.5, and 1.0)

X-ray diffraction (XRD) analysis of Fe-doped PrBa0.5Sr0.5Co2−FeO5+ (x = 0, 0.5, and 1.0) samples suggests that they had a single phase without detectible amount of impurity after firing at 1100 ~ 1150°C in air for 12 hours (Supplementary Fig. S1a). Rietveld analysis of the XRD data indicates that the PrBa0.5Sr0.5Co1.5Fe0.5O5+ (PBSCF05) compound has a tetragonal structure (P4/mmm) with a = 3.871 Å and c = 7.757 Å (Supplementary Fig. S1c). The lattice expands with increasing Fe content as the smaller Co3+(r = 0.545 Å) or Co4+(r = 0.530 Å) ions are replaced by the larger Fe3+(r = 0.645 Å) or Fe4+(r = 0.585 Å) ions (Supplementary Table S1). As temperature increased from room temperature to operating temperature, the cell volume was expanded, but the structure (P4/mmm) remains the same (Supplementary Fig. S2 and Supplementary Table S2).

Transmission electron microscopy (TEM) analyses (Fig. 1) of the as-synthesized and the annealed (700°C for 600 hours) PBSCF05 sample suggest that the tetragonal structure has lattice constants of a = 3.871 Å and c = 7.757 Å, similar to the XRD analysis (Supplementary Fig. S1a and Supplementary Table S1). Further, the high-resolution TEM image (Fig. 1b) also confirms the ordered stacking sequence [Ba(Sr)O]-[Co(Fe)O2]-[PrO]-[Co[Fe]O2]-[Ba(Sr)O2] of the structure. In particular, this crystal structure demonstrated remarkable stability at high temperatures; the atomic structure of the surface remains unchanged after annealing at 700°C for 600 hours (Fig. 1c and Supplementary Fig. S4), in stark contrast to simple perovskite cathode materials such as BSCF or LSCF. The latter is much less stable; surface segregation of Sr-enriched phases or SrO has been observed under similar conditions14, leading to gradual degradation in performance. As a result, a surface coating must be applied to suppress surface phase segregation22.

Figure 1

Transmission electron microscopy (TEM) analysis (a) A bright-field (BF) TEM image and an electron diffraction (ED) pattern obtained from an as-synthesized PrBa0.5Sr0.5Co2-FeO5+ sample. (b) A high-resolution TEM image of a grain in (a). (c) A BF-TEM image and an ED pattern of the sample annealed at 700°C for 600 hours. (d) A high-resolution TEM image of a grain after annealing.

Electrical and redox properties

The electrical conductivities of each sample at a given temperature (Fig. 2a) increased with increasing partial pressure of oxygen, p(O2), an indication of a typical p-type conductor. This is ideally suited for cathode application where electron holes are generated (or electrons are consumed) by the ORR. In general, the cathode experiences a lower p(O2) under fuel cell operating conditions due to cathodic polarization23. Therefore, sufficient electrical conductivity at relatively low p(O2) is important to ensure efficient current collection and long-term stability. For all samples examined in this study, the electrical conductivities vary from 10 to 103 S cm−1 under typical operating conditions, which are sufficient for the cathode application. As p(O2) was decreased below ~10−6 atm at 700°C, the conductivity dropped dramatically due to partial reduction of the oxide lattice2425. However, this stability limit can be shifted toward lower p(O2) by increasing the content of Fe doping.

Figure 2

Electrical and chemical properties of PrBa0.5Sr0.5Co2−FeO5+ (a) Electrical conductivities of PrBa0.5Sr0.5Co2−FeO5+ (x = 0, 0.5, and 1.0) measured at 700°C in different p(O2). (b) Oxygen non-stoichiometry of PrBa0.5Sr0.5Co1.5Fe0.5O5+ and Ba0.5Sr0.5Co0.8Fe0.2O3−δ as a function of p(O2) at 700°C.

To gain more insight into the redox properties, oxygen nonstoichiometry of the PBSCF05 samples, together with that of BSCF (for comparison), was characterized as a function of p(O2) at 700°C using coulometric titration (Fig. 2b). The isotherm of ordered PBSCF05 has a lower critical p(O2) for decomposition than that of BSCF, implying that PBSCF05 has higher redox stability or better durability under cathodic polarization, a favourable property for practical SOFC applications.

Power output and durability of fuel cells

The electrocatalytic activity for ORR of these materials was determined in symmetrical cells using impedance spectroscopy. The cathodic polarization resistance normalized by the geometric electrode area, the area specific resistance (ASR), can be readily calculated from impedance spectra acquired under open circuit conditions (inset of Fig. 3a). When the electronic conduction in the electrolyte is negligible, the cathodic polarization resistance can be adequately approximated by the diameter of the impedance loop (or the difference between the intersects with the real axis at high and low frequencies); otherwise, it has to be calculated with proper corrections to take into consideration of the partial shorting effect of electronic conduction through the electrolyte26. Summarized in Fig. 3a are the ASRs of the PrBa0.5Sr0.5Co2−FeO5+ (x = 0, 0.5, and 1.0) electrodes as determined from impedance spectra. It is noted that, when x = 0.5, the PBSCF05 showed the lowest ASR in the temperature range studied, ~0.33 Ω cm2 at 500°C and ~0.056 Ω cm2 at 600°C, which are lower than those reported for BSCF under similar conditions (e.g., 0.7 Ω cm2 at 500°C). In general, the population of mobile oxygen defects may contribute to the enhanced oxygen kinetics associated with oxygen bulk diffusion and surface exchange15. Higher concentration of mobile oxygen defects in the Ln-O layer due to larger amount of Fe doping (up to ~50%) may lead to faster oxygen kinetics and better electrochemical performance.

Figure 3

Electrochemical properties of PrBa0.5Sr0.5Co2−FeO5+δ in fuel cells (a) Arrhenius plot of reciprocal ASR for PrBa0.5Sr0.5Co2−FeO5+-GDC (x = 0, 0.5, and 1.0) cathode. The inset shows the ASR of PBSCF05-GDC cathode measured at 600°C in air under open-circuit conditions. (b) I–V curves and the corresponding power densities of test cells with PBSCF05-GDC cathode using humidified H2 (3% H2O) as the fuel and ambient air as the oxidant to the cathode at 500 ~ 650°C. (c) Peak power densities of cells with LnBa0.5Sr0.5Co1.5Fe0.5O5+-GDC (Ln = Pr and Nd) cathode. (d) Short term stability measurement for a test cell, Ni-GDC | GDC | PBSCF05-GDC, at a constant cell voltage of 0.6 V at 550°C.

To characterize the performance of these new PBSCF cathode materials in a practical fuel cell, we used anode-supported cells based on a ~15 μm thick GDC electrolyte. Fig. 3b shows some typical I–V curves and the corresponding power densities of the test cells with PBSCF05-GDC as the cathodes at 500 ~ 650°C. The peak power densities of the cells were 0.71, 1.31, 2.16, and 2.90 W cm−2 at 500, 550, 600, and 650°C, respectively, much higher than those reported for the cells using a BSCF cathode (e.g., 1.01 W cm−2 at 600°C) which are prepared and tested under similar conditions. Fig. 3c shows the average peak power densities of cells with a cathode of LnBa0.5Sr0.5Co1.5Fe0.5O5+ (Ln = Pr and Nd)-GDC measured at different operating temperatures. The performance of a Ni-GDC | GDC | PBSCF05-GDC test cell was very stable under a cell voltage of 0.6 V at 550°C for 150 hours (Fig. 3d), demonstrating stable power output without observable degradation.

DFT prediction of possible elementary pathway for the ORR

To gain some insight into the structural features of these cathode materials, we performed density functional theory (DFT) calculations2728. Based on our XRD and TEM analyses (Supplementary Table S1 and Fig. 1), we constructed a model for PrBa0.5Sr0.5Co1.5Fe0.5O6.0 (PBSCFO6.0) (Supplementary Fig. S10). Calculations of oxygen-vacancy formation energy (EOV) in this model suggest that it is energetically more favorable to form oxygen vacancies in the vicinity of Co ions than that of Fe ions (Supplementary Fig. S11 and Table S3; EOV = 1.18 eV versus 1.36 eV). When two oxygen ions connected to the Co ion are removed, it reduces to PrBa0.5Sr0.5Co1.5Fe0.5O5.875 (PBSCFO5.875), forming a pore channel (Supplementary Fig. S12). However, we found that the activation energy for oxygen diffusion through the channel in PBSCFO5.875 is ~2.6 eV, which is much higher than the experimental values (~0.55 eV). To lower the reaction barrier, more oxygen ions connected to the Co ions were removed (Supplementary Fig. S12). Based on our DFT modeling results, several types of pore channels (Fig. 4a) could be proposed for the defective structures; one of the oxygen ion diffusion paths may follow a zig-zag type trajectory through the CoO plane perpendicular to the PrO plane (Fig. 4b), because the reaction barrier (~0.46 eV) for PrBa0.5Sr0.5Co1.5Fe0.5O5.75 (PBSCFO5.75) and PrBa0.5Sr0.5Co1.5Fe0.5O5.625 (PBSCFO5.625) is closer to the experimental value (~0.55 eV).

Figure 4

DFT calculations for elucidating the most probable elementary pathway for the ORR on the PBSCFO (a) Schematic illustration of PBSCFO with a pore channel labeled with a round box on the (010) plane. To illustrate the pore channel through the PrO plane and the CoO plane (perpendicular to the PrO plane), only the PrO and the CoO layers were shown. (b) A proposed mechanism for the surface ORR and the bulk diffusion via the pore channels in PrBa0.5Sr0.5Co2−FeO5+ (010). VE is an oxygen vacancy in the GDC electrolyte.

Further, a 2-D surface model was also constructed to simulate the interactions between O2 and the PBSCO(010) surface (Supplementary Fig. S13 and S14). DFT calculations suggest that the most probable elementary pathway for the ORR on the PBSCFO cathode can be described as follows (Fig. 4b and Supplementary Table S4): (1) adsorption and dissociation of molecular O2 on an oxygen vacancy near a Co ion, (2) incorporation of the dissociated oxygen ions into the pore channels; (3) diffusion of the oxygen ions through the pore channels (bulk diffusion), and (4) combination of the oxygen ions with oxygen vacancies in the GDC electrolyte.

Discussion

In conclusion, a class of cation-ordered, double-perovskite compounds display fast oxygen ion diffusion through pore channels and high catalytic activity toward ORR at low temperatures while maintaining excellent compatibility with electrolyte and good stability under typical fuel cell operating conditions. DFT analysis using simplified models suggests that the most attractive properties of these materials are the pore channels in the [PrO] and [CoO] planes that could provide fast paths for oxygen transport, which in turn accelerates the kinetics of surface oxygen exchange. More detailed understanding of the mechanistic details may help to rationally design better double-perovskite cathode materials for a new generation of high-performance SOFCs with enhanced durability.

Methods

Synthesis of electrode and electrolyte powders

Cathode materials LnBa0.5Sr0.5Co2−FeO5+ (LnBSCF) (Ln = Pr and Nd; x = 0, 0.25, 0.5, 0.75, and 1.0) were synthesized using a glycine-nitrate process (GNP). Stoichiometric amounts of Ln(NO3)3·6H2O(Aldrich, 99.9%, metal basis), Ba(NO3)2 (Aldrich, 99+%), Sr(NO3)2 (Aldrich, 99+%), Co(NO3)2·6H2O (Aldrich, 98+%), and Fe(NO3)3·6H2O (Aldrich, 98%) were dissolved in distilled water with proper amount of glycine. The solutions were heated up to 350°C in air and followed by combustion to form fine powders, which were calcined at 600°C for 4 hours. The resulting powders were then grinded and calcined again at 900°C for 4 hours. For the measurement of electrical conductivity and coulometric titration, the powders were pressed into pellets at 5 MPa and sintered in air at 1100 ~ 1150°C for 12 hours (to achieve relative density > 97%). For transmission electron microscopy (TEM), the PrBa0.5Sr0.5Co1.5Fe0.5O5+ (PBSCF05) sample was annealed at 700°C for 600 hours. The Ce0.9Gd0.1O2−(GDC) powders for electrolyte and NiO-GDC powders for anode were also synthesized using the GNP method. For preparation of cathode slurries, pre-calcined cathode and GDC powders (at a weight ratio of 60:40) were mixed using ball milling, together with an organic binder (Heraeus V006).

Fabrication of fuel cells and TEM sample preparation

Symmetrical cells with a configuration of electrode | GDC | electrode were used for impedance spectroscopy. GDC powder was pressed into pellets of ~1 mm thick and sintered at 1350°C for 4 hours in air to obtain a dense electrolyte membrane. Cathode slurries were then painted onto both surfaces of dense GDC electrolyte and fired at 950°C for 4 hours. Ni-GDC anode-supported cells with a configuration of Ni-GDC | GDC | cathode were fabricated using a drop-coating method. NiO powder, GDC powder, and starch (weight ratio of 6:4:1.5) were mixed using ball milling in ethanol for 24 hours. After drying, the NiO-GDC mixture was pressed into a pellet (~0.6 mm thick and 15 mm diameter). Thin GDC electrolyte membranes were prepared by a refined particle suspension coating technique. A GDC suspension was prepared by dispersing GDC powders (Aldrich) in ethanol with a small amount of binder (polyvinyl Butyral, B-98) and dispersant (Triethanolamine, Alfa Aesar) at a ratio of 1:10. The GDC suspension was applied to a NiO-GDC anode support by drop-coating, followed by drying in air and subsequent co-sintering at 1400°C for 5 hours. Cross-sectional samples for TEM analysis were prepared using a focused ion beam (FIB, Quanta 3D, FEI).

Electrical and electrochemical testing

Electrical conductivities of the LnBSCF cathode materials were determined in air using a four-electrode measurement. All four electrodes were made of Ag wire and Ag paste. The current and voltage were controlled/measured using a potentiostat (BioLogic) in the temperature range of 100 to 750°C with an interval of 50°C.

For symmetrical cells, two Ag wires were attached to each of the two electrodes using Ag paste. Each cell was mounted on an alumina tube using a ceramic adhesive (Aremco, Ceramabond 552). Impedance spectra were recorded under OCV in a frequency range of 1 mHz to 500 kHz with ac perturbation of 10 mA in the temperature range of 500 ~ 650°C. For the single cell tests, each cell was mounted on an alumina tube using a ceramic adhesive. Humidified (with 3 v% H2O) H2 was used as the fuel at a flow rate of 100 mL min−1 (passing through a water bubbler at 25°C), whereas ambient air was supplied to cathode as the oxidant. Impedance spectra and I–V polarization curves were obtained with a BioLogic Potentiostat. The I–V polarization curves were recorded between 500°C and 650°C.

Redox property and oxygen nonstoichiometric

The redox properties and oxygen nonstoichiometry of LnBSCF cathodes and Ba0.5Sr0.5Co0.8Fe0.2O3− (BSCF) were measured using coulometric titration (CT) as a function of the oxygen partial pressure, p(O2). The detailed CT procedure is as described elsewhere2425. After purging 5% O2-Ar gas over the sample in the tube for 24 hours, p(O2) was determined from the OCV. The sample was allowed to equilibrate with the surrounding atmosphere until the change in potential was less than 1 mV h−1. Oxygen nonstoichiometry was determined through this procedure at 700°C over a wide range of oxygen partial pressure. Electrical conductivity was measured using a four-electrode configuration and a BioLogic Potentiostat.

DFT calculation

The Vienna ab initio simulation package (VASP)2729 was used for periodic density functional theory (DFT) calculations. As reported28, it is difficult to apply for the DFT + U approach since the concurrent optimization of two Ueff parameters of two B-site cations of PrBa0.5Sr0.5Co2−FeO5+ (PBSCFO) is problematic. Thus in this study, we carried out the spin-polarization method with the generalized gradient approximation (GGA) with the projector-augmented-wave method (PAW)30 using the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. To simulate the oxygen-deficient double perovskite PBSCFO, we constructed a tetragonal P4/mmm structure (a = b ≠ c). The kinetic energy cutoff for a plane wave basis set was 415 eV. The Brillouin zone was sampled with the (3 × 3 × 3) and (3 × 3 × 1) Monkhorst-Pack mesh k-points31 for bulk and surface calculations. For the 2-D surface calculations, the slabs were separated by a vacuum space of 15 Å in the direction perpendicular to the surface. In this study, the adsorption energy (Ead) on a CoFeO(010)-terminated surface (1/16 ML) was calculated by Ead = E[O-PBSCFO] − E[PBSCFO] − 1/2E[O2], where E[O- PBSCFO], E[PBSCFO], and E[O2] are the predicted electronic energies for an adsorbed O species on a PBSCFO surface, a bare PBSCFO surface, and a gas-phase triplet O2, respectively. The oxygen-vacancy formation energy (EOV)2832 was calculated by EOV = E[defective PBSCFO] + 1/2E[triplet O2] − E[perfect PBSCFO], where E[defective PBSCFO] and E[perfect PBSCFO] are the predicted electronic energies for defective and perfect bulk PBSCFO structures, respectively.

Other Characterization

X-ray powder diffraction (XRD) (Rigaku diffractometer, Cu Ka radiation) analysis was used to confirm the crystalline structures of samples. In situ XRD was obtained from room to operating temperature (Bruker, D8 Advance). The microstructures and morphologies of LnBSCF cathode samples were observed using a field emission scanning electron microscope (SEM) (Nova SEM). The TEM images were acquired with JEOL JEM 2100F with a probe forming (STEM) Cs corrector at 200 kV. A thermogravimetric analysis (TGA) was carried out using a SDT-Q600 (TA instrument, USA). TGA experiments were performed from 100°C to 900°C with a heating/cooling rate of 2°C min−1 in air. The room-temperature oxygen content values were determined by iodometric titration.

Author Contributions

S.C., S.Y., J.K., S.P., A.J., S.S. and J.K. fabricated samples and carried out analytical and electrochemical characterizations. H.Y.J. collected and analyzed the TEM data. Y.C. performed the DFT calculations. All authors contributed to writing the paper. J.S., G.K. and M.L. conceived and coordinated the project.