# Oxide diffusion in innovative SOFC cathode materials

Y. Hu, V. Thoréton, C. Pirovano, E. Capoen, Christine Bogicevic, N. Nuns, A.-S. Mamede, G. Dezanneau, R.N. Vannier

► To cite this version:

Y. Hu, V. Thoréton, C. Pirovano, E. Capoen, Christine Bogicevic, et al.. Oxide diffusion in innovative SOFC cathode materials. Faraday Discussions, 2014, 176, pp.31 - 47. 10.1039/C4FD00129J. hal-01790588

HAL Id: hal-01790588

https://hal.science/hal-01790588v1

Submitted on 30 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

HAL
open science

# Oxide diffusion in innovative SOFC cathode materials

Y. Hu, V. Thoréton, C. Pirovano, E. Capoen, Christine Bogicevic, N. Nuns, A.-S. Mamede, G. Dezanneau, R.N. Vannier

<figure><img src="image_3.png" /></figure>

► To cite this version:

Y. Hu, V. Thoréton, C. Pirovano, E. Capoen, Christine Bogicevic, et al.. Oxide diffusion in innovative SOFC cathode materials. Faraday Discussions, 2014, 176, pp.31 - 47. 10.1039/C4FD00129J. hal-01790588

HAL Id: hal-01790588

https://hal.science/hal-01790588v1

Submitted on 30 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
PAPER

View Article Online
View Journal | View Issue

# Oxide diffusion in innovative SOFC cathode materials

Y. Hu,ª V. Thoréton,ᵇ C. Pirovano,ᵇ E. Capoen,ᵇ C. Bogicevic,ª N. Nuns,ᵇ A.-S. Mamede,ᵇ G. Dezanneauª and R. N. Vannier*ᵇ

Received 11th June 2014, Accepted 28th July 2014

DOI: 10.1039/c4fd00129j

Oxide diffusion was studied in two innovative SOFC cathode materials, Ba₂Co₉O₁₄ and Ca₃Co₄O₉+δ derivatives. Although oxygen diffusion was confirmed in the promising material Ba₂Co₉O₁₄, it was not possible to derive accurate transport parameters because of an oxidation process at the sample surface which has still to be clarified. In contrast, oxygen diffusion in the well-known Ca₃Co₄O₉+δ thermoelectric material was improved when calcium was partly substituted with strontium, likely due to an increase of the volume of the rock salt layers in which the conduction process takes place. Although the diffusion coefficient remains low, interestingly, fast kinetics towards the oxygen molecule dissociation reaction were shown with surface exchange coefficients higher than those reported for the best cathode materials in the field. They increased with the strontium content; the Sr atoms potentially play a key role in the mechanism of oxygen molecule dissociation at the solid surface.

## Introduction

With the possibility of cogeneration of electricity and heat, solid oxide fuel cells are promising. The first units were commercialised in Japan in 2009 and demonstrators are running in Europe. Their principles are based on a dense ceramic electrolyte which separates two compartments, one containing hydrogen, and the other containing air. Due to a difference in oxygen chemical potential between the two compartments, oxygen molecules are reduced at the cathode into oxide ions, which migrate through the electrolyte to the anode to react with hydrogen molecules and produce water. Hydrogen is usually produced by methane reforming, however, internal reforming of methane is also envisaged. The typical temperatures used range between 700 and 1000 °C. The main electrolyte materials which are currently studied are Yttria Stabilised Zirconia (YSZ),

<footnote>ªLaboratoire Structures, Propriétés et Modélisation des Solides (SPMS), Ecole Centrale Paris, Grande Voie des Vignes, 92295 Châtenay-Malabry, France</footnote>

<footnote>ᵇUniv Lille Nord de France, CNRS UMR 8181, Unité de Catalyse et de Chimie du Solide – UCCS, ENSCL, Université Lille 1, CS 90108, 59652 Villeneuve d'Ascq Cedex, France. E-mail: rose-noelle.vannier@ensc-lille.fr; Fax: +33(0)320436814; Tel: +33(0)320436583</footnote>
Gadolinia Doped Ceria (GDC) and lanthanum gallates (LSGM). At the anode, a Ni-YSZ cermet containing a 40% volume of nickel is usually used, but it is prone to carbon deposition when used with hydrocarbon fuels. It is also sensitive to sulphur poisoning, and suffers from nickel agglomeration and redox stability with prolonged usage. Alternative materials, such as chromo-manganites (La<sub>0.75</sub>Sr<sub>0.25</sub>)Cr<sub>0.5</sub>Mn<sub>0.5</sub>O<sub>3</sub> or titanates such as lanthanum-doped SrTiO<sub>3</sub> have been proposed. At the cathode, strontium-substituted lanthanum manganite (La<sub>1-x</sub>Sr<sub>x</sub>MnO<sub>3-δ</sub>, LSM) is conventionally used. However, it displays low oxide ion conductivity, and composites with YSZ have been developed in order to allow use at temperatures lower than 800 °C. To operate at lower temperatures, Mixed Ionic Electronic Conductors (MIEC) are usually preferred in order to allow the oxygen reduction reaction to be sprayed on the whole electrode volume, and not only at the electrode/electrolyte interface. As an alternative to LSM, the La<sub>0.6</sub>Sr<sub>0.4</sub>Co<sub>0.2</sub>Fe<sub>0.8</sub>O<sub>3-δ</sub> (LSCF) perovskite has been widely studied.¹ When used with a GDC electrolyte and after optimisation of its microstructure, an area specific resistance (ASR) of only 0.13 Ω cm² was obtained at 603 °C for this compound,² the target for the cathode being 0.15 Ω cm² at 700 °C. This high performance is likely due to the presence of Co(III) in the structure, which is known to be particularly active for the reduction of oxygen molecules into oxide ions.³ Indeed, cobaltites are known to display high oxide ion transport parameters. Among these, cobaltite with a composition of La<sub>0.6</sub>Sr<sub>0.4</sub>Co<sub>0.2</sub>O<sub>3-δ</sub> is one of the best catalysts for the oxygen molecule reduction reaction with a surface exchange coefficient (k*, cm s⁻¹) of 1.1 × 10⁻⁶ cm s⁻¹ associated with an oxygen self-diffusion tracer (D*, cm² s⁻¹) coefficient of 2.0 × 10⁻⁸ cm² s⁻¹ at 680 °C under 230 mbar.⁴ One drawback of these materials is their high thermal expansion coefficient (TEC). The aforementioned composition has a value of 23 × 10⁻⁶ K⁻¹, which is twice the thermal expansion of YSZ (10 × 10⁻⁶ K⁻¹)⁵ and GDC (12 × 10⁻⁶ K⁻¹).⁶ However, partial substitution of cobalt with iron in LSCF led to a good compromise with a TEC in the range of 14–15 × 10⁻⁶ K⁻¹ from 500 to 700 °C.⁷ More recently, promising performances as SOFC cathodes have been shown for two innovative cobaltites, Ba₂Co₉O₁₄ and Ca₃Co₄O₉+δ. An ASR of 0.5 Ω cm² was first measured on a symmetric cell made of a 70 wt% Ba₂Co₉O₁₄–30 wt% GDC composite electrode deposited on a GDC electrolyte.⁸ The good electrochemical performance of this compound was later confirmed by Li *et al.*, who obtained ASRs of only 0.133 Ω cm² and 0.068 Ω cm² at 750 and 850 °C, and powder densities of 450 mW cm⁻² and 770 mW cm⁻² on a cell made of a LSGM electrolyte with a Ba₂Co₉O₁₄–samarium-doped ceria composite cathode at 800 and 850 °C, respectively.³ This compound exhibits a three dimensional structure based on [Co₈O₈] blocks sandwiched in between [Ba₂CoO₆] layers (Fig. 1a).⁹,¹⁰ With cobalt atoms at different valences, it displays high electronic conductivity with a value of 240 S cm⁻¹ at 650 °C, but as with other cobaltites, its thermal expansion coefficient is high (22 × 10⁻⁶ K⁻¹), which necessitates the use of a composite, at least for mechanical compatibility, between the electrode and the electrolyte. After optimisation of its composition and thickness, an ASR of only 0.08 Ω cm² at 700 °C was recently measured on a symmetric cell made of a 50 wt% Ba₂Co₉O₁₄–50 wt% GDC composite electrode with a GDC electrolyte.¹¹

In contrast, a thermal expansion coefficient of only 10 × 10⁻⁶ K⁻¹ was measured for the well-known Ca₃Co₄O₉+δ thermoelectric compound,¹² which has a misfit structure composed of CdI₂-type [CoO₂] layers sandwiched in between
<figure><img src="image_4.png" /><figcaption>Fig. 1 Crystal structures of (a) Ba₂Co₉O₁₄ and (b) Ca₃Co₄O₉+δ.</figcaption></figure>

[Ca₂CoO₃-δ] rock salt slabs (Fig 1b).¹³,¹⁴ With an electronic conductivity of 100 S cm⁻¹ at 700 °C and oxygen vacancies in the rock salt layers, mixed ionic electronic conductivity was expected for this compound. When tested in a symmetrical cell with a GDC electrolyte, ASRs of 4.94 Ω cm² for the pure compound and 1 Ω cm² for a 70 wt% Ca₃Co₄O₉+δ-30 wt% GDC composite electrode were first measured.¹⁵,¹⁶ After optimisation of the composition and thickness, this was decreased to 0.5 Ω cm² at 700 °C for a 50 wt% Ca₃Co₄O₉+δ-50 wt% GDC composite electrode with a thickness of 21 µm, and further decreased to 0.35 Ω cm² when calcium was partly substituted with 10% of strontium.¹⁷ These good performances were confirmed by Samson *et al.*,¹⁸ who reported a polarization resistance of 0.64 Ω cm² at 600 °C for a 50 vol% Ca₃Co₄O₉+δ-50 vol% GDC composite electrode. Using a samarium-doped ceria electrolyte, Zou *et al.*¹⁹ measured a maximum power density of 430 mW cm⁻² at 700 °C for a 70 wt% Ca₂.₉Bi₀.₁Co₄O₉+δ-30 wt% Ce₁.₈Sm₀.₂O₁.₉₅| Ce₁.₈Sm₀.₂O₁.₉₅|Ce₁.₈Sm₀.₂O₁.₉₅-Ni button cell.

More recently, oxygen diffusion in this material was confirmed by combining ¹⁸O/¹⁶O isotope exchange, electrical conductivity relaxation and pulse isotopic exchange techniques.²⁰ Despite low diffusion with a coefficient of 2.7 × 10⁻¹⁰ cm² s⁻¹ at 700 °C, this study showed high kinetics toward oxygen dissociation with a surface exchange coefficient as high as 1.6 × 10⁻⁷ cm s⁻¹.

In the present paper, oxygen diffusion in these two innovative cathodes was investigated, with emphasis placed on the impact of the partial substitution of calcium with strontium on the oxygen transport properties of Ca₃Co₄O₉+δ.

## Experimental

Ba₂Co₉O₁₄ and (Ca₁-xSrx)₃Co₄O₉+δ (x = 0.1, 0.2) powders were prepared by solid-state reaction from the corresponding carbonates and oxides BaCO₃ (Prolabo – 99.5%, dried at 100 °C), CaCO₃ (Prolabo – 99.5%, dried at 100 °C), SrCO₃ (Aldrich – 99.8%, dried at 100 °C) and Co₃O₄ (Alfa Aesar, 99.7%). Reactants were carefully ground in an agate mortar and then calcined repeatedly at 900 °C with intermediate grinding. The purity of the samples was checked by X-ray diffraction (D8 Advance Bruker AXS diffractometer, Cu Kα radiation, 1D LynxEye PSD detector). Fullprof software [v2.05, July 2011] was used for refinement of the unit cell parameters.²¹ An Accupyc II 1340 micromeritic pychnometer was used to measure the powder density. The absolute oxygen stoichiometry of the sample was determined by iodometric titration.²²
Thermal analyses were conducted under air, using a TG/DTA Setaram Setsys apparatus. 40 mg of the sample powder was introduced into a Pt crucible, and the experiment was carried out at a heating rate of 5 °C min⁻¹ under flowing air (0.4 L h⁻¹) from room temperature to 950 °C. A blank experiment with an empty crucible was performed in order to correct the data for the effects of buoyancy. Thermogravimetry was also carried out for the evaluation of the oxygen stoichiometry of the sample as a function of pO₂. In this case, a Setaram 92-1750 apparatus was used. The experiment was carried out using 50 mg of powder in a Pt crucible under the flow of an O₂/N₂ mixture (5 L h⁻¹), with three successive oxygen partial pressures (100, 200 and 400 mbar), and the temperature ranging from 600 to 750 °C (50 °C step) for each pressure. For this set-up, the buoyancy was negligible and the data were directly used to derive the thermodynamic factor at each temperature.

Dense ceramics, 3 mm thick, were prepared by Spark Plasma Sintering at the IdF Platform in Thiais (France). A 20 mm (diameter) carbon die was used. The samples were sintered within 15 minutes with a 2 minute dwell at 850 °C, under a 50 MPa (8.8 kN) pressure. To remove carbon residues, the pellets were annealed for 12 h at 800 °C. A Hitachi S-3400N instrument was used to characterise the surface of the samples.

For isotopic exchange and electrochemical conductivity relaxation, sample surfaces were roughly polished in the first step, before they were polished down to a 1/4th of a micron with successive SiC papers and diamond pastes (Ra ~ 10~20 nm) in the second step. In the case of calcium cobaltites, to avoid the evolution of the surface roughness, samples were annealed for 72 h at 900 °C before the second polishing in order to release residual strains in the material.

The ¹⁸O/¹⁶O Isotope Exchange Depth Profiling technique was applied to all samples in order to derive the oxygen self-diffusion tracer (D*, cm² s⁻¹) and surface exchange coefficients (k*, cm s⁻¹). The diffusion time was chosen to get diffusion lengths significantly smaller than the sample dimensions to consider a 1D diffusion model perpendicular to the surface of the dense ceramic in a semi-infinite medium. In these conditions, assuming a first order surface exchange reaction and constant ¹⁸O concentration during the exchanges, the solution for a semi-infinite solid to the Fick's 2nd law of diffusion is given by the following equation:²³

$$C'(x, t) = \frac{C(x, t) - C_{bg}}{C_g - C_{bg}} = \text{erfc}\left[\left(\frac{x}{2 \times \sqrt{D^* \times t}}\right)\right] - \left[\exp\left(\frac{k^*x}{D^*} + \frac{k^{*2}t}{D^*}\right) \times \text{erfc}\left(\frac{x}{2 \times \sqrt{D^* \times t}} + k^* \sqrt{\frac{t}{D^*}}\right)\right] \quad (1)$$

where Cbg is the natural ¹⁸O ratio, Cg is the ¹⁸O ratio in the enriched exchange atmosphere, x is the distance from the surface of the specimen, and t is the time of the isotope exchange.

The exchange conditions are reported in Table 1. Experiments were carried out in between 650 and 750 °C under a pressure close to 210 mbar with pre-annealing
<table><caption>Table 1 Conditions for the ¹⁸O/¹⁶O isotope exchange</caption><thead><tr><th>Temperature (°C)</th><th>Pressure (mbar)</th><th>Pre-annealing time (hours)</th><th>Annealing time (min)</th><th>¹⁸O ratio in the gas phase</th></tr></thead><tbody><tr><td colspan="5"><strong>Ba<sub>2</sub>Co<sub>9</sub>O<sub>14</sub></strong></td></tr><tr><td>648</td><td>210</td><td>240</td><td>2860</td><td>0.79</td></tr><tr><td>700</td><td>210</td><td>96</td><td>2765</td><td>0.80</td></tr><tr><td>750</td><td>218</td><td>96</td><td>2790</td><td>0.80</td></tr><tr><td colspan="5"><strong>(Ca<sub>1-x</sub>Sr<sub>x</sub>)<sub>3</sub>Co<sub>4</sub>O<sub>9+δ</sub> x = 0.10, 0.20</strong></td></tr><tr><td>604</td><td>207</td><td>160</td><td>725</td><td>0.88</td></tr><tr><td>649</td><td>211</td><td>100</td><td>567</td><td>0.67</td></tr><tr><td>694</td><td>219</td><td>61</td><td>392</td><td>0.80</td></tr><tr><td>741</td><td>210</td><td>20</td><td>55</td><td>0.9</td></tr></tbody></table>

under normal oxygen (search grade oxygen 99.995% ALPHAGAZ) to reach thermodynamic equilibrium prior to the exchange under labelled oxygen. For both annealings, under the ¹⁸O enriched atmosphere or pre-annealing, the samples were annealed for the required time, and quenched. To determine the ¹⁸O content in the oxygen gas used for the exchange, a silicon wafer was oxidised overnight at 1000 °C under 200 mbar.

A ToF-SIMS⁵ machine (ION-TOF GmbH, 25 keV $Bi^+$ LMIG analysis gun, DSCS sputter gun) was used to measure the ¹⁸O diffusion profiles using the same experimental protocol as described in ref. 20. Both depth analysis from the sample surface (depth profile mode)²⁴ and edge analysis of the cut sample (line scan mode)²⁵ were performed and the profiles merged. The depth of the craters, after SIMS depth profiling, were measured with a KLA Tencor Alpha-step IQ profilometer. The same apparatus was used to measure the samples' roughness.

Electrical conductivity relaxation was also carried out on strontium-doped $Ca_3Co_4O_{9+δ}$. For this purpose, plane-parallel samples (17 mm in length and 1 to 2.4 mm in width) were cut from the strontium substituted dense ceramics using a diamond saw. The electrical conductivity was measured using a 4 point configuration while changing the oxygen concentration of the surrounding atmosphere, as shown in ref. 20. A gold paste was used to paint gold contacts which were sintered at 750 °C for 1 h at a rate of 100 °C h⁻¹ before the experiment. The distance between the two internal electrodes was about 14 mm.

Once the sample was at equilibrium at the given temperature, the $pO_2$ was changed step-wise and the conductivity was monitored as a function of time. The variation in $pO_2$ was performed by an abrupt switch between air and mixtures of $O_2$/Ar with oxygen partial pressures of 0.05 atm and 0.01 atm, respectively, at a 5 L h⁻¹ flow rate. In addition, the temperature-dependent electrical conductivity was measured under air with a 5 L h⁻¹ flow rate at a very slow speed of 25 °C h⁻¹.

To characterise the nature of the atoms at the uppermost surface of the sample, Low-Energy Ion Scattering (LEIS) was also applied to a dense ceramic ($Ca_{1-x}Sr_x)_{3}Co_{4}O_{9+δ}$ with $x = 0.1$ composition which had been heated at 800 °C in air during the 12 hours prior to the analysis. The analysis was performed using a Qtac¹⁰⁰ (ION-TOF GmbH) apparatus. Before analysis, the sample was subjected in situ to an atomic oxygen treatment ($2-3 \times 10^{-4}$ mbar) to remove possible organic
surface contamination. The study was carried out with a $^4$He$^+$ (3 keV) beam on a 1 mm$^2$ analysed area.

# Results and discussion

## Characterization of powder and ceramics

The purity of the synthesised powders and dense ceramics was confirmed by X-ray diffraction. As shown in Fig. 2, good pattern matching was obtained for both type of ceramics.

Corresponding unit-cell parameters are given in Table 2. For Ba$_2$Co$_9$O$_{14}$, they are in good agreement with the parameters reported by Ehora *et al.* (5.6963(8) Å, 28.924(6) Å)$^9$ and Sun *et al.* (5.6958(4) Å, 28.909(4) Å).$^{10}$ The small differences in parameters we observed between the powders and the ceramics, on one hand, and these authors, on the other hand, are likely due to small differences in the oxygen content due to different thermal histories. For the refinement of (Ca$_{1-x}$Sr$_x$)$_3$Co$_4$O$_{9+δ}$, the phases were considered as a mixture of two phases, one corresponding to the [CoO$_2$] layers and the second one corresponding to the [Ca$_2$CoO$_{3-δ}$] rock salt layers, with common *a*, *c* and *β* cell parameters constrained

<figure><img src="image_5.png" /><figcaption>Fig. 2 The XRD calculated profile compared to experimental data collected on dense pellets of Ba$_2$Co$_9$O$_{14}$ (a), and (Ca$_{1-x}$Sr$_x$)$_3$Co$_4$O$_{9+δ}$ with x = 0.10 (b).</figcaption></figure>

<figure><img src="image_6.png" /></figure>
<table><caption>**Table 2** Ba₂Co₉O₁₄ and (Ca₁⁻ₓSrₓ)₃Co₄O₉+δ x = 0.1, 0.2 unit cell parameters from XRD on the powder and dense ceramics</caption><thead><tr><th></th><th>Powder</th><th>Dense ceramics</th></tr></thead><tbody><tr><td><strong>Ba<sub>2</sub>Co<sub>9</sub>O<sub>14</sub> (R&#x0333;3m)</strong></td><td></td><td></td></tr><tr><td>a (Å)</td><td>5.6945(3)</td><td>5.700(1)</td></tr><tr><td>b (Å)</td><td>28.905(2)</td><td>28.926(6)</td></tr><tr><td><strong>(Ca<sub>1-x</sub>Sr<sub>x</sub>)<sub>3</sub>Co<sub>4</sub>O<sub>9+δ</sub> x = 0.10</strong></td><td></td><td></td></tr><tr><td>a (Å)</td><td>4.8429(3)</td><td>4.8433(9)</td></tr><tr><td>b<sub>H</sub> (Å)</td><td>2.8255(3)</td><td>2.825(4)</td></tr><tr><td>b<sub>RS</sub> (Å)</td><td>4.5968(6)</td><td>4.563(2)</td></tr><tr><td>c (Å)</td><td>10.8730(4)</td><td>10.881(2)</td></tr><tr><td>β (°)</td><td>98.137(9)</td><td>98.14(3)</td></tr><tr><td><strong>(Ca<sub>1-x</sub>Sr<sub>x</sub>)<sub>3</sub>Co<sub>4</sub>O<sub>9+δ</sub> x = 0.20</strong></td><td></td><td></td></tr><tr><td>A (Å)</td><td>4.8518(3)</td><td>4.8547(6)</td></tr><tr><td>b<sub>H</sub> (Å)</td><td>2.8225(5)</td><td>2.821(1)</td></tr><tr><td>b<sub>RS</sub> (Å)</td><td>4.6304(8)</td><td>4.6362(8)</td></tr><tr><td>c (Å)</td><td>10.8995(7)</td><td>10.904(2)</td></tr><tr><td>β (°)</td><td>98.123(7)</td><td>98.12(2)</td></tr></tbody></table>

to be the same, and different $b_H$ and $b_{RS}$, refined independently. The evolution of the unit cell parameters with the strontium content is in good agreement with the results reported by Wang *et al.*,²⁶ who showed almost no variation of the $b_H$ unit cell characterising the [CoO₂] layer periodicity, and an increase of the $b_{RS}$ unit cell with strontium content, in good agreement with the higher radius of Sr²⁺ (1.18 Å) compared to Ca²⁺ (1.00 Å).²⁷

As shown by Ehora *et al.*,⁹ Ba₂Co₉O₁₄ decomposed into CoO and a cubic BaCoO₄~2~ form at 1030 °C. To check the stability of the calcium cobaltites, TGA and DTA analyses were carried out on the two strontium-doped compositions. Their thermograms are compared to the signals measured for the non-doped compound in Fig. 3. They exhibit similar behavior. The first transition, which is observed at 130 °C for the pure compound, disappears when strontium is introduced into the structure. The second transition at 550 °C, corresponding to the beginning of a mass loss, is maintained at the same temperature, but the intensity of the associated endothermic peak on heating decreases with the increase in strontium content. All compounds decomposed before 950 °C, and the decomposition temperature was slightly lowered when the strontium content was increased. This decomposition process is reversible, but it is worth noting that recombination happens at a lower temperature when the substitution rate increases.

Due to the lack of stability at high temperature, it was not possible to sinter these compounds using conventional sintering, but dense ceramics with relative densities higher than 96% were obtained using spark plasma sintering, as confirmed by SEM (see Fig. 4).

¹⁸O/¹⁶O isotope exchange in Ba₂Co₉O₁₄

The normalised ¹⁸O/¹⁶O + ¹⁸O) isotope concentrations as a function of the penetration depth are given in Fig. 5 for the three studied temperatures. Both data
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>Fig. 3 Comparison of TGA and DTA analyses corresponding to the Ca₃Co₄O₉+δ parent compound (CCO, black solid line), (Ca₁-xSrx)₃Co₄O₉+δ with x = 0.10 (CSCO10, red line) and (Ca₁-xSrx)₃Co₄O₉+δ with x = 0.20 (CSCO20, black dotted line).</figcaption></figure>

recorded in the depth-profile and line-scan modes are reported. A good agreement between the two sets of data was observed. However, whatever the temperature of exchange, a high oxygen content was measured at the surface of the samples within the first micron of the depth; the thickness of this high oxygen content layer increased with temperature. It was therefore not possible to fit these data to the diffusion equation (eqn (1)).

Taking into account only the data obtained far from the surface, the transport parameters given in Table 3 were derived.

Surprisingly, despite very good electrochemical performances, rather low oxygen diffusion and kinetics towards oxygen molecule dissociation were shown. However, one must be careful with these results. The increase in ¹⁸O at the sample surface is not clearly understood. It is likely due to an oxidation
<figure><img src="image_9.png" /><figcaption>Fig. 4 SEM images of Ba₂Co₉O₁₄ (a), and (Ca₁-xSrx)₃Co₄O₉+δ with x = 0.10 (b) and 0.20 (c), showing the dense ceramics.</figcaption></figure>

process which occurs when the sample is cooled, but this hypothesis still has to
be confirmed.

Oxygen transport in strontium-doped Ca₃Co₄O₉

The oxygen transport properties in the two strontium-doped Ca₃Co₄O₉+δ
materials were characterised by both electrical conductivity relaxation and the
¹⁸O/¹⁶O isotope depth profile technique. The electrical conductivity at the
thermodynamic equilibrium corresponding to the x = 0.10 composition under
different oxygen pressures is given in Fig. 6a. An increase in the electrical
conductivity was observed when the oxygen partial pressure increased, con-
firming the p-type behaviour of these materials. As shown in Fig. 6b, the
electrical conductivity slightly decreased when strontium was introduced into
the structure. The evolution of the electrical conductivity corresponding to the

<figure><img src="image_10.png" /><figcaption>Fig. 5 Normalised ¹⁸O diffusion profile measured on Ba₂Co₉O₁₄ dense ceramics after annealing under labelled oxygen at 650 °C, 700 °C and 750 °C for about 46 hours and corresponding fit (in red) taking into account the data obtained far from the surface.</figcaption></figure>
<table><caption>**Table 3** Tentative Ba₂Co₅O₁₄ transport parameters derived from normalised ¹⁸O isotope profile far from the sample edges</caption><thead><tr><th>T (°C)</th><th>k* (cm s<sup>-1</sup>)</th><th>D* (cm<sup>2</sup> s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>648</td><td>1.3 × 10<sup>-10</sup></td><td>2.5 × 10<sup>-12</sup></td></tr><tr><td>700</td><td>4.0 × 10<sup>-10</sup></td><td>5.0 × 10<sup>-12</sup></td></tr><tr><td>750</td><td>2.5 × 10<sup>-9</sup></td><td>6.0 × 10<sup>-12</sup></td></tr></tbody></table>

(Ca₁⁻ₓSrₓ)₃Co₄O₉+δ (x = 0.1) composition at 750 °C with the changes in surrounding atmosphere is shown Fig. 7. It is worth noting that oxidation and reduction are reversible processes.

Data obtained for both compositions in the 650–750 °C range were fitted to the theoretical model of diffusion. In these conditions, only the diffusion coefficient could be derived with precision, reflecting a D<sub>chem</sub>/k<sub>chem</sub> characteristic length that was much smaller than the size of the sample, and suggesting much higher surface exchange kinetics for these materials than for chemical diffusion.²⁸ Corresponding D<sub>chem</sub> values are given in Table 4 and shown in Fig. 8. The D<sub>chem</sub> values obtained when the air ↔ 0.05 atm and air ↔ 0.01 atm are close, and are similar to what was obtained for the pure material.²⁰ Moreover, the chemical diffusion coefficients are independent of whether the relaxation is a reduction or an oxidation process.

As shown in Fig. 8, as expected, the diffusion process is thermally activated, and activation energies of 1.03 ± 0.03 eV and 0.98 ± 0.04 eV were measured for the compositions with x = 0.1 and x = 0.2, respectively, which were slightly lower than the activation energy displayed by the pure compound (1.2 ± 0.02 eV), indicating that diffusion was easier for the strontium-doped compound, which also exhibits higher oxygen diffusion coefficients.

<figure><img src="image_11.png" /><figcaption>**Fig. 6** Evolution of the electrical conductivity of (Ca₁⁻ₓSrₓ)₃Co₄O₉+δ with an x = 0.10 composition at the thermodynamic equilibrium under different oxygen pressures (a), and comparison of the electrical conductivities of the x = 0 (CCO in black), x = 0.10 (CSCO10 in red) and x = 0.20 (CSCO 20 in blue) compositions as a function of temperature under air and 0.05 atm (b).</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>Fig. 7 Evolution of the electrical conductivity of (Ca$_{1-x}$Sr$_x$)$_3$Co$_4$O$_{9+δ}$ with an $x = 0.10$ composition at 750 °C with surrounding atmosphere changes, showing the reversibility of the oxidation and reduction processes.</figcaption></figure>

This increase in the oxygen diffusion coefficient with the substitution ratio was confirmed by $^{18}$O/$^{16}$O isotope exchange. The normalized $^{18}$O diffusion profile corresponding to the (Ca$_{1-x}$Sr$_x$)$_3$Co$_4$O$_{9+δ}$ with $x = 0.1$ composition, exchanged at 649 °C over 567 minutes under 211 mbar dry oxygen, is shown in Fig. 9. Again, a good agreement was obtained between data collected in depth profile mode and line scan mode.

However, in this case, as shown for the parent compound,²⁰ a small decrease in the $^{18}$O content was noticed at the sample surface, which is attributed to a possible carbonation or hydroxylation when back in the ambient atmosphere. In the inset, the same plot has been drawn on a decimal logarithmic scale to compare the experimental and theoretical profiles in depth. As previously noticed for the pure compound, a “tail” can be seen, most likely suggesting anisotropic

<table><caption>Table 4 $D_{chem}$ (cm$^2$ s$^{-1}$) values corresponding to (Ca$_{1-x}$Sr$_x$)$_3$Co$_4$O$_{9+δ}$ with $x = 0.1$ and $x = 0.2$</caption><thead><tr><th>T (°C)</th><th>Air-0.05 atm</th><th>0.05 atm-air</th><th>Air-0.01 atm</th><th>0.01 atm-air</th></tr></thead><tbody><tr><td colspan="5">(Ca<sub>1-x</sub>Sr<sub>x</sub>)<sub>3</sub>Co<sub>4</sub>O<sub>9+δ</sub> x = 0.10</td></tr><tr><td>650</td><td>1.37 × 10<sup>-7</sup></td><td>1.35 × 10<sup>-7</sup></td><td>1.34 × 10<sup>-7</sup></td><td>1.31 × 10<sup>-7</sup></td></tr><tr><td>675</td><td>1.92 × 10<sup>-7</sup></td><td>1.91 × 10<sup>-7</sup></td><td>1.85 × 10<sup>-7</sup></td><td>1.95 × 10<sup>-7</sup></td></tr><tr><td>700</td><td>2.58 × 10<sup>-7</sup></td><td>2.51 × 10<sup>-7</sup></td><td>2.52 × 10<sup>-7</sup></td><td>2.57 × 10<sup>-7</sup></td></tr><tr><td>725</td><td>3.47 × 10<sup>-7</sup></td><td>3.61 × 10<sup>-7</sup></td><td>3.44 × 10<sup>-7</sup></td><td>3.53 × 10<sup>-7</sup></td></tr><tr><td>750</td><td>4.68 × 10<sup>-7</sup></td><td>4.95 × 10<sup>-7</sup></td><td>4.73 × 10<sup>-7</sup></td><td>4.76 × 10<sup>-7</sup></td></tr><tr><td colspan="5">(Ca<sub>1-x</sub>Sr<sub>x</sub>)<sub>3</sub>Co<sub>4</sub>O<sub>9+δ</sub> x = 0.20</td></tr><tr><td>650</td><td>3.09 × 10<sup>-7</sup></td><td>3.17 × 10<sup>-7</sup></td><td>3.04 × 10<sup>-7</sup></td><td>3.25 × 10<sup>-7</sup></td></tr><tr><td>675</td><td>4.62 × 10<sup>-7</sup></td><td>4.54 × 10<sup>-7</sup></td><td>4.57 × 10<sup>-7</sup></td><td>4.58 × 10<sup>-7</sup></td></tr><tr><td>700</td><td>5.99 × 10<sup>-7</sup></td><td>6.40 × 10<sup>-7</sup></td><td>5.87 × 10<sup>-7</sup></td><td>6.37 × 10<sup>-7</sup></td></tr><tr><td>725</td><td>8.41 × 10<sup>-7</sup></td><td>8.50 × 10<sup>-7</sup></td><td>8.00 × 10<sup>-7</sup></td><td>8.24 × 10<sup>-7</sup></td></tr><tr><td>750</td><td>1.04 × 10<sup>-6</sup></td><td>1.14 × 10<sup>-6</sup></td><td>0.96 × 10<sup>-6</sup></td><td>1.08 × 10<sup>-6</sup></td></tr></tbody></table>
<figure><img src="image_13.png" /><figcaption>Fig. 8 Evolution of $D_{chem}$ values as a function of temperature, measured on the parent compound $Ca_3Co_4O_{9+δ}$ (CCO, green square), and $(Ca_{1-x}Sr_x)_3Co_4O_{9+δ}$ with $x = 0.10$ (CSCO10, red dot) and $x = 0.20$ (CSCO20 blue triangle).</figcaption></figure>

diffusion.²⁰ Nevertheless, as shown in Fig. 9, a good fit between the calculated and experimental data was obtained. The derived $k^*$ and $D^*$ parameters are given in Table 5. Their evolution as a function of temperature is shown on an Arrhenius plot in Fig. 10.

The tracer diffusion can be correlated to the chemical diffusion via the thermodynamic factor $γ$, according to the following equation:

$$D_{chem} = \gamma D^*$$

<figure><img src="image_14.png" /><figcaption>Fig. 9 Normalised ¹⁸O diffusion profile measured on the $(Ca_{1-x}Sr_x)_3Co_4O_{9+δ}$ with $x = 0.1$ composition, exchanged at 649 °C over 567 minutes under 211 mbar dry oxygen. The calculated profile is given in red.</figcaption></figure>
$$ \gamma = 1/2 \frac{d \ln pO_2}{d \ln C_O} \approx 1/2 \frac{\Delta \ln pO_2}{\Delta \ln C_O} \quad \text{when } \sigma_e \gg \sigma_i \quad (2) $$

where $\Delta \ln pO_2$ and $\Delta \ln C_O$ are the corresponding variations in the oxygen partial pressure and oxygen content, respectively, and $\sigma_e$ and $\sigma_i$ are the electronic and ionic conductivities, respectively.

For the sake of comparison, the evolution of $D_{\text{ECR}}^* = D_{\text{chem}}/\gamma$ is also reported in Fig. 10. The values are in good agreement. The slightly higher values observed for $D_{\text{ECR}}^*$ compared to the tracer diffusion coefficient values can be explained by a lack of accuracy in the determination of the thermodynamic coefficient because of small oxygen stoichiometry variations (see the values of the thermodynamic coefficient in Table 6).

An increase in the diffusion coefficient with strontium content is confirmed, and may be explained by the increase in the volume of the rock salt layers, which may facilitate oxygen mobility within these layers. However, with a value in the range of $10^{-10} \text{ cm}^2 \text{ s}^{-1}$, it is still low, which justifies the addition of gadolinia-doped ceria to improve the electrochemical performance.$^{15-19}$

With surface exchange coefficients of $3.0 \times 10^{-6} \text{ cm s}^{-1}$ and $4.0 \times 10^{-6} \text{ cm s}^{-1}$ at 700 °C for both the $x = 0.1$ and $x = 0.2$ compositions, these materials exhibit very fast kinetics for the oxygen reduction reaction compared to the best cathode materials in the field, since values of only $1.3 \times 10^{-7}$, $1.1 \times 10^{-6}$ and $1.3 \times 10^{-7} \text{ cm s}^{-1}$ were reported for $La_2NiO_{4+\delta}$ at 700 °C,$^{29}$ $La_{0.6}Sr_{0.4}CoO_{3-\delta}$ at $680 \text{ °C}$ and $GdBaCoO_{2.5+\delta}$ at $686 \text{ °C}$,$^{30}$ respectively.

To further understand oxygen transport in these materials, the $x = 0.1$ composition was characterised by LEIS. The initial spectrum, corresponding to the estimated removal of 0.2 monolayers, is compared to the one measured after the estimated removal of 3 monolayers in Fig. 11. Despite bulk contamination with sodium, which was also observed at the surface of the powders and ceramics of the parent compound,$^{20}$ mainly strontium and calcium were initially observed at the uppermost sample surface. No cobalt was initially observed, just a sub-surface signal was shown on the low energy side from the surface peak$^{31}$ in the background, indicating that cobalt atoms are beneath the first measured atomic layer. After ionic sputtering, cobalt was evidenced at the outermost surface, and an increase in calcium was also noticed. However, one

<table><caption>Table 5 $D^*$ (cm$^2$ s$^{-1}$) and $k^*$ (cm s$^{-1}$) values corresponding to $(Ca_{1-x}Sr_x)_3Co_4O_{9+\delta}$ with $x = 0.1$ and $x = 0.2$</caption><thead><tr><th>T (°C)</th><th>604</th><th>649</th><th>694</th><th>741</th></tr></thead><tbody><tr><td colspan="5">(Ca<sub>1-x</sub>Sr<sub>x</sub>)<sub>3</sub>Co<sub>4</sub>O<sub>9+δ</sub> x = 0.10</td></tr><tr><td>k*</td><td>1.4 × 10<sup>-8</sup></td><td>1.9 × 10<sup>-7</sup></td><td>3.0 × 10<sup>-6</sup></td><td>3.0 × 10<sup>-6</sup></td></tr><tr><td>D*</td><td>1.1 × 10<sup>-10</sup></td><td>2.5 × 10<sup>-10</sup></td><td>4.1 × 10<sup>-10</sup></td><td>1.0 × 10<sup>-9</sup></td></tr><tr><td colspan="5">(Ca<sub>1-x</sub>Sr<sub>x</sub>)<sub>3</sub>Co<sub>4</sub>O<sub>9+δ</sub> x = 0.20</td></tr><tr><td>k*</td><td>6.7 × 10<sup>-8</sup></td><td>3.0 × 10<sup>-7</sup></td><td>4.0 × 10<sup>-6</sup></td><td>2.0 × 10<sup>-6</sup></td></tr><tr><td>D*</td><td>3.0 × 10<sup>-10</sup></td><td>2.8 × 10<sup>-10</sup></td><td>6.1 × 10<sup>-10</sup></td><td>6.7 × 10<sup>-10</sup></td></tr></tbody></table>
<figure><img src="image_15.png" /><figcaption>Fig. 10 Arrhenius plot of k* and D* values compared to D<sub>chem</sub>/γ for the composition (Ca<sub>1-x</sub>Sr<sub>x</sub>)<sub>3</sub>Co<sub>4</sub>O<sub>9+δ</sub> with x = 0.1 (a) and x = 0.2 (b).</figcaption></figure>

<table><caption>Table 6 (Ca<sub>1-x</sub>Sr<sub>x</sub>)<sub>3</sub>Co<sub>4</sub>O<sub>9+δ</sub> (x = 0.1, x = 0.2) thermodynamic coefficients derived from TGA analysis as a function of temperature and oxygen partial pressure</caption><thead><tr><th>Temp. (°C)</th><th>x = 0.1</th><th>x = 0.2</th></tr></thead><tbody><tr><td>650</td><td>211</td><td>313</td></tr><tr><td>675</td><td>200</td><td>305</td></tr><tr><td>700</td><td>189</td><td>297</td></tr><tr><td>725</td><td>183</td><td>299</td></tr><tr><td>750</td><td>178</td><td>302</td></tr></tbody></table>

must be careful with this technique, since strontium may also segregate at the sample surface. The same behaviour was observed for other compounds. For SmCoO<sub>3</sub>, Fullarton *et al.* also showed that samarium was preferentially at the surface with only 5% of cobalt at the outermost surface, and segregation of strontium was evident at the surface for strontium-doped compounds.<sup>32</sup> Similarly, Viitanen *et al.* showed the absence of cobalt and iron in the outermost atomic layer of a La<sub>0.6</sub>Sr<sub>0.4</sub>Co<sub>0.2</sub>Fe<sub>0.8</sub>O<sub>3</sub> membrane.<sup>33</sup> We also observed calcium oxide at the uppermost surface of the Ca<sub>3</sub>Co<sub>4</sub>O<sub>9</sub> parent compound. Here, strontium seems to be mainly at the sample surface, and its presence may partly explain the increase in the kinetics towards the oxygen dissociation reaction when the strontium content increases.
<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /><figcaption>Fig. 11 $^4$He$^+$ (3 keV) LEIS spectra recorded on dense ceramics of $(Ca_{1-x}Sr_x)_3Co_4O_{9+δ}$ with $x = 0.1$ after annealing at 800 $^\circ$C for 12 hours prior to the analysis. The first spectrum (in black) corresponds to an estimated removal of 0.2 monolayers, and the last spectrum (in red) corresponds to an expected removal of around 3 monolayers.</figcaption></figure>

## Conclusions

Although oxygen diffusion was confirmed in the promising material $Ba_2Co_9O_{14}$, it was not possible to derive accurate transport parameters because of an oxidation process at the sample surface which still has to be clarified. In contrast, for $Ca_3Co_4O_{9+δ}$ derivatives, it was shown that oxygen diffusion was improved when calcium was partly substituted with strontium, likely due to an increase in the volume of the rock salt layers of the structure, which may help oxygen mobility. Although the diffusion coefficient remains low, interestingly, fast kinetics towards the oxygen dissociation reaction were shown, with a surface exchange coefficient higher than those reported for the best cathode materials in the field. The surface exchange coefficient is increased when the strontium content increases, which may be correlated to the presence of a higher amount of strontium at the sample surface, which may play a key role in the mechanism of oxygen molecule dissociation at the surface of the solid.

## Acknowledgements

The "Fonds Européen de Développement Régional (FEDER)", "CNRS", "Région Nord Pas-de-Calais" and French "Ministère de l'Éducation Nationale de l'Enseignement Supérieur et de la Recherche" are acknowledged for funding of X-ray and surface analysis facilities. V.T. is grateful to the University Lille 1 for funding his PhD grant. The authors are also grateful to the Agence Nationale de la Recherche for funding the ANR EVEREST (ANR-07-PANH-005), NanOxyDesign (ANR-10-BLAN-0814) and IDEA-MAT (ANR-12-BS08-0012-01) projects. Maxence Vandewalle, Nora Djelal and Laurence Burylo from UCCS are acknowledged for SEM, TG-DTA analyses and conventional XRD experiments, respectively.
References

1 *Functional materials for sustainable energy applications*, Woodhead publishing series in energy, ed. J. A. Kilner, S. J. Skinner, S. J. C. Irvine and P. P. Edwards, Woodhead Publishing Limited, 2012.

2 D. Marinha, J. Hayd, L. Dessemond, E. Evers-Tiffée and E. Djurado, *J. Power Sources*, 2011, **196**, 5084.

3 Y. Li, M. W. Xu and J. B. Goodenough, *J. Power Sources*, 2012, **209**, 40.

4 A. V. Berenov, A. Atkinson, J. A. Kilner, E. Bucher and W. Sitte, *Solid State Ionics*, 2010, **181**, 819.

5 O. Yamamoto, *Electrochim. Acta*, 2000, **45**, 2423.

6 K. Kordesch and G. Simader, *Fuel Cells and their Applications*, VCH Wiley, NY, Weinheim, 1996.

7 B. Fan, J. Yan and X. Yan, *Solid State Sci.*, 2011, **13**, 1835.

8 A. Rolle, N. Preux, G. Ehora, O. Mentré and S. Daviero-Minaud, *Solid State Ionics*, 2011, **184**, 31.

9 G. Ehora, S. Daviero-Minaud, M. Colmont, G. André and O. Mentré, *Chem. Mater.*, 2007, **19**, 2180.

10 J. Sun, M. Yang, G. Li, T. Yang, F. Liao, Y. Wang, M. Xiong and J. Lin, *Inorg. Chem.*, 2006, **45**, 9151.

11 I. Kehal, to be published.

12 K. Nagasawa, S. Daviero-Minaud, N. Preux, A. Rolle, P. Roussel, H. Nakatsugawa and O. Mentré, *Chem. Mater.*, 2009, **21**, 4738.

13 H. Muguerra, D. Grebille and F. Bourée, *Acta Crystallogr., Sect. B: Struct. Sci.*, 2008, **64**, 676.

14 C. D. Ling, K. Aivazian, S. Schmid and P. Jensen, *J. Solid State Chem.*, 2007, **180**, 1446.

15 A. Rolle, S. Boulfrad, K. Nagasawa, H. Nakatsugawa, O. Mentré, J. Irvine and S. Daviero-Minaud, *J. Power Sources*, 2011, **196**, 7328.

16 A. Rolle, V. Thoréton, P. Rozier, E. Capoen, O. Mentré, B. Boukamp and S. Daviero-Minaud, *Fuel Cells*, 2012, **12**(2), 288.

17 A. Rolle, to be published.

18 A. J. Samson, M. Søgaard, N. Van Nong, N. Pryds and N. J. Bonanos, *J. Power Sources*, 2011, **196**, 10606.

19 J. Zou, J. Park, H. Yoon, T. Kim and J. Chung, *Int. J. Hydrogen Energy*, 2012, **37**, 8592.

20 V. Thoréton, Y. Hu, C. Pirovano, E. Capoen, N. Nuns, A.-S. Mamede, G. Dezanneau, C. Y. Yoo, H. J. M. Bouwmeester and R. N. Vannier, *J. Mater. Chem. A*, 2014, **2**, 19717.

21 J. Rodriguez-Carvajal, *Phys. B*, 1993, **192**, 55.

22 L. Karvonen, H. Yamauchi and M. Karppinen, *Chem. Mater.*, 2008, **20**, 7143.

23 J. Crank, *The Mathematics of diffusion*, Clarendon Press, Oxford, 2nd edn, 1975.

24 J. A. Kilner, B. C. H. Steele and L. Ilkov, *Solid State Ionics*, 1984, **12**, 89.

25 R. J. Chater, S. Carter, J. A. Kilner and B. C. H. Steele, *Solid State Ionics*, 1992, **53–56**, 859.

26 L. B. Wang, A. Maignan, D. Pelloquin, S. Hébert and B. Raveau, *J. Appl. Phys.*, 2002, **92**, 124.
27 R. D. Shannon, *Acta Crystallogr., Sect. A: Cryst. Phys., Diffr., Theor. Gen. Crystallogr.*, 1976, **32**, 751.

28 M. W. den Otter, H. J. M. Bouwmeester, B. A. Boukamp and H. Verweij, *J. Electrochem. Soc.*, 2001, **148**(2), J1.

29 R. Sayers, R. A. De Souza, J. A. Kilner and S. J. Skinner, *Solid State Ionics*, 2010, **181**, 386.

30 A. Tarancón, S. J. Skinner, R. J. Chater, F. Hernández-Ramirez and J. A. Kilner, *J. Mater. Chem.*, 2007, **17**, 3175.

31 H. R. J. ter Veen, T. Kim, I. E. Wachs and H. H. Brongersma, *Catal. Today*, 2009, **140**, 197.

32 I. C. Fullarton, J. P. Jacobs, H. E. Benthem, J. A. Kilner, H. H. Brongersma, P. J. Scanlon and B. C. H. Steele, *Ionics*, 1995, **1**, 51.

33 M. M. Viitanen, R. G. Welzenis, H. H. Brongersma and F. P. F. van Berkel, *Solid State Ionics*, 2002, **150**, 223.