Implications of Cation Interdiffusion between Double Perovskite Cathode and Proton-Conducting Electrolyte for Performance of Solid Oxide Fuel Cells

: Chemical compatibility and cation interdiffusion between the double perovskite cobaltites R BaCo 2 O 6 − δ ( R = Gd, Pr) and proton-conducting electrolyte BaZr 0.8 Y 0.2 O 3 − δ were studied. Chemical interaction was found to occur already at 1100 ◦ C as a result of the partial dissolution of R BaCo 2 O 6 − δ ( R = Gd, Pr) in BaZr 0.8 Y 0.2 O 3 − δ . Analysis of the element distribution along the cross sections of diffusion couples R BaCo 2 O 6 − δ ( R = Gd, Pr)|BaZr 0.8 Y 0.2 O 3 − δ showed strong interdiffusion of cations, with cobalt being the most mobile one. Its diffusion depth in the electrolyte reaches up to several hundreds of micrometers. The addition of NiO as a sintering aid to BaZr 0.8 Y 0.2 O 3 − δ promotes cation diffusion especially through the grain boundary mechanism, increasing the diffusion depth of Co. The possible implications of cation interdiffusion on the performance of proton-conducting SOFCs are discussed based on the results obtained.


Introduction
Solid oxide fuel cells (SOFCs) are solid-state electrochemical devices that convert the chemical energy of the fuel (hydrogen, natural gas and other hydrocarbons) directly into electrical energy and possess high energy conversion efficiency and low environmental impact [1][2][3].At the same time, a number of problems impede large-scale applications of SOFCs, among them high-temperature degradation processes, chemical compatibility of cell components and the high cost of materials [1][2][3][4].The lowering of the operating temperature is the general way to promote the commercial application of SOFCs because it broadens the range of compatible construction materials, prolongs their lifetime and, consequently, reduces the cost of energy production.In this respect, proton-conducting solid oxide fuel cells (H-SOFCs) are generally considered to be better suited for practical applications due to lower operating temperatures and the absence of fuel contamination.In turn, conclusions on the chemical compatibility of materials are usually made based only on a costly and time-consuming trial-and-error method: heating a physical mixture of different powders at high temperature for several hours followed by analyzing its phase composition using X-ray diffraction (XRD).Such an approach may lead to contradictory results because of the kinetic limitations of solid-state interaction and non-zero detection limits of conventional XRD.In this respect, thermodynamic analysis of materials' compatibility is completely free from such limitations and, probably, represents the only way of unambiguous interpretation of compounds' stability under particular conditions.
The choice of thermodynamically stable and chemically compatible materials is a necessary but not sufficient condition for successful operation of SOFCs, because solidstate materials exposed to a thermodynamic potential gradient (gradient of temperature, chemical potential of elements, electrical potential) are subjected to kinetic degradation governed by diffusion processes [5].The applied gradients induce directed fluxes of the mobile components and lead to three basic degradation phenomena of materials: (i) kinetic demixing, (ii) kinetic decomposition and (iii) morphological instability [5].Thus, in SOFCs the different chemical natures of components and high operating temperatures (600-1000 • C) inevitably lead to cation interdiffusion and continuous kinetic degradation of materials and, therefore, determine the lifetime of the devices [4].
Obviously, understanding the diffusion processes at the electrode/electrolyte interfaces in SOFCs is of critical importance for their successful large-scale commercial application [4].Unfortunately, this important information is mostly lacking, and only a limited number of studies of cation diffusion phenomena in perovskite oxides may be found in the literature.Thus, Ca-and Sr-doped LaCrO 3 , and Sr-and Mg-doped LaGaO 3 were investigated by tracer diffusion method [6][7][8].Cation interdiffusion in the systems LaFeO 3 , LaCoO 3 , La 0.9 Sr 0.1 NiO 4 , La 2 Ni 0.8 Cu 0.2 O 4 and La 27 W 5 O 55. 5 was studied by diffusion couple method [9][10][11][12][13].The kinetic demixing process was studied only for La 0.5 Sr 0.5 Fe 0.5 Co 0.5 O 3−δ and La 0.3 Sr 0.7 CoO 3−δ [14,15].Comprehensive investigation was carried out for La 2 NiO 4+δ , a Ruddlesden-Popper phase, including determination of the self-diffusion coefficient of Ni 2+ [16], chemical tracer diffusion [17], interdiffusion coefficients for the Aand B-site cations [18], as well as kinetic decomposition of La 2 NiO 4+d under an oxygen chemical potential gradient [19].As for the interdiffusion between the components of SOFCs, to the best of our knowledge this has been studied only for LSCF|GDC|YSZ interfaces, where Gd-doped ceria (GDC) acts as a buffer layer between the La 1−x Sr x Co 1−y Fe y O 3−δ (LSCF) cathode and yttria-stabilized zirconia (YSZ) electrolyte [20,21].In turn, for protonconducting SOFCs, such studies have not yet been performed.
The thermodynamic analysis of the stability of the double perovskite PrBaCo 2 O 6−δ under different conditions has recently been carried out by us [22].As a result, its instability against chemical interaction with CeO 2 was found.The diffusion couple (PrBaCo 2 O 6−δ (PBC)-SDC) experiment showed interdiffusion of Pr and Sm and chemical interaction between PBC and SDC with the formation of barium cerate as a product, in full agreement with thermodynamic calculations.Both phenomena were shown to affect significantly both thermal expansion and total conductivity of the composite electrodes and lead to higher than expected TECs of the composites.Thus, the potential problems with cation interdiffusion have already been outlined [22].Formation of barium cerate on the border between PBC and SDC was also shown to cause a significant drop in the composites' total conductivity and should be detrimental for the long-term stability of the SOFC cell.One can expect that formation of BaCO 3 due to the chemical reaction of PBC with CO 2 from the ambient air, as predicted by thermodynamic calculations, will lead to lowering of the SOFC's performance as well.The PBC double perovskite was shown to be more suitable as an electrode for proton-conducting solid oxide fuel cells because of its thermodynamic stability against chemical interaction with barium-cerateand zirconate-based electrolytes [22].
In addition, the influence of the sintering aids such as, for example, widely employed NiO on the electrode-electrolyte compatibility and interdiffusion is completely unknown.

Experimental Procedure
Powder samples of the double perovskites PrBaCo 2 O 6−δ (PBC) and GdBaCo 2 O 6−δ (GBC) and the solid electrolyte BaZr 0.8 Y 0.2 O 3−δ (BZY0.2) were prepared by means of a glycerol-nitrate method using Gd 2 O 3 , Pr 6 O 11 , Y 2 O 3 , BaCO 3 , Co and Zr(OH) 2 CO 3 •xH 2 O as starting materials.The purity of all the materials used was 99.99% (wt.).Zr(OH) 2 CO 3 •xH 2 O was preliminarily dissolved in concentrated HNO 3 (purity 99.99%), and the Zr concentration in the resulting solution was determined by the gravimetric method.The as-obtained solution with known concentration of Zr was used as a precursor for the solid electrolyte synthesis.
For preparation of BaZr 0.8 Y 0.2 O 3−δ , the stoichiometric quantities of BaCO 3 and Y 2 O 3 were dissolved in concentrated HNO 3 .This solution was mixed with the Zr(OH) 2 CO 3 •xH 2 O solution described above, and a stoichiometric amount of glycerol (99% purity) was added as a complexing agent and a fuel.The required weight of glycerol was calculated according to full reduction of the nitrates to molecular N 2 .The as-obtained solution was heated continuously at 100 • C till complete water evaporation and pyrolysis of the dried precursor.The resulting ash was subsequently calcined at 1100 • C for 24 h in air and then pressed into pellets (10 mm in diameter at 150 MPa) and sintered at 1550 • C for 12 h in air.A composite material of BaZr 0.8 Y 0.2 O 3−δ with 1 wt.% of NiO as a sintering aid was made by physical mixing of the components in an agate mortar in an ethanol medium.
For synthesis of the RBaCo 2 O 6−δ (R = Gd, Pr) powders, a stoichiometric mixture of the appropriate starting materials was dissolved in concentrated HNO 3 , and glycerol (99% purity) was added to the solution.The weight of glycerol was calculated as described above in the case of BaZr 0.8 Y 0.2 O 3−δ .As-prepared solutions were heated continuously at 100 • C till complete water evaporation and pyrolysis of the dried precursor.To yield the RBaCo 2 O 6−δ (R = Gd, Pr) powders, the product of pyrolysis was subsequently calcined at 900, 1000 and 1100 • C for 12 h in air followed by thorough regrinding in an ethanol medium in an agate mortar after each annealing step.
The phase composition of the as-synthesized samples was studied by means of X-ray diffraction (XRD) with an XRD 7000 diffractometer (Shimadzu, Kyoto, Japan) using Cu Kα radiation (λ = 1.5418Å).The phase identification was performed using the PDF-2 2021 database [23] and Match! software [24].The XRD patterns demonstrating phase purity of all the as-prepared oxides are presented in Figure S1.
The refinement of the lattice parameters was performed using the Le Bail method as implemented in Rietica 4.0 software [25].The refined unit cell parameters of the starting double perovskite and electrolyte powders are given in Table S1.
The chemical compatibility of RBaCo NiO (1 wt.%).Before annealing, each diffusion couple was pressed in a special device and fixed with a clamping mechanism.For each system, two independent experiments were performed, both at 1100 • C for 20 h and at 1200 • C for 48 h in air.After annealing, diffusion couples were placed in a silicon form, fixed with epoxy resin and cut perpendicular to the reaction interface.Then cross section of the reaction zone was polished down to 1 µm using diamond suspension, purified ultrasonically in distilled water and dried at 70 • C in a drying oven in air atmosphere.The microstructure and element distribution across the couples' sections were examined by means of scanning electron microscopy (SEM) using a VEGA 3 microscope (Tescan, Brno, Czech Republic) equipped with an Ultim Max 40 (Oxford Instruments, Abingdon, UK) detector for energy dispersive X-ray spectroscopy (EDX).
Additionally, diffusion couples with porous layers of RBaCo 2 O 6−δ (R = Gd, Pr) on top of the dense ceramic pellet of BaZr 0.8 Y 0.2 O 3−δ were studied to mimic the behavior of SOFCs' cathodes.These diffusion couples were prepared by screen-printing the cathode paste on the polished surface of the ceramic electrolyte pellet.The cathode paste was obtained by thoroughly mixing the powder of double perovskite cobaltite with a solution of polyvinyl butyral (5 wt.%) in ethanol.The as-prepared diffusion couple was annealed at 1100 • C for 6 h in air.The treatment of the reaction zone and subsequent investigation by SEM were performed as described above for the other diffusion couples.

Chemical Compatibility
The results of chemical compatibility experiments are shown in Figures 1 and S2.At a first glance, XRD patterns of the sample mixtures look very similar, at least, at temperatures lower than 1100 • C. Such XRD patterns can be readily found in the current literature along with a typical conclusion that they indicate the absence of chemical interaction between the components of the mixture investigated.However, the XRD pattern of the GdBaCo 2 O 6−δ + BaZr 0.8 Y 0.2 O 3−δ mixture calcined at the highest studied temperature (see Figure 1) clearly indicates the formation of a new compound isostructural to Gd 2 O 3 (possibly doped by Y 2 O 3 ).Furthermore, the intensities of the peaks corresponding to the components of the mixtures as well as peaks' positions (and, hence, lattice parameters) depend on the calcination temperature as seen in Figure 1c,d Therefore, the similarity of the XRD patterns of the starting mixture (denoted as 25 • C in Figure 1) and the one calcined at 1100 • C cannot be used as evidence of the absence of chemical interaction.
As seen in Figure 2a,b, lattice parameters of RBaCo 2 O 6−δ (R = Gd, Pr) oxides, first slightly increase but then decrease with increasing calcination temperature.This can be explained by the combined effect of two factors: (i) variation of the RBaCo 2 O 6−δ oxygen content, 6−δ, in a wide range [26,27] with temperature and the sample's prehistory and (ii) partial substitution of the smaller Y 3+ cation for larger Gd 3+ and Pr 3+ [28] in RBaCo 2 O 6−δ due to cation interdiffusion between components of the sample mixture.
As for BaZr 0.8 Y 0.2 O 3−δ , its unit cell parameter, as seen in Figure 2c, decreases with calcination temperature and more strongly in contact with GdBaCo 2 O 6−δ oxide.This can be, at least qualitatively, explained assuming cobalt dissolution in BaZr 0.8 Y 0.2 O 3−δ leads to lattice contraction as shown in [29], whereas dissolution of Gd and Pr, on the contrary, facilitates the lattice expansion of barium zirconate due to differences in the cation radii [28].

Analysis of Diffusion Couples
The results of SEM and EDX analysis of the cross sections of diffusion couples RBaCo 2 O 6−δ (R = Gd, Pr)-BaZr 0.8 Y 0.2 O 3−δ +NiO (0, 1 wt.%) after annealing are presented in Figures 3-5.They support the conclusions drawn above and based on the results of chemical compatibility study by XRD.As seen in Figures 3-5, cation interdiffusion can be observed, with Co being the most mobile cation demonstrating pronounced diffusion into the electrolyte pellet for all the couples studied.The largest diffusion depth, reaching 1000 µm, was found in the PrBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ +NiO (1 wt.%) couple annealed at 1200 • C for 48 h in air.Again, the diffusion of Pr and Co occurs through both grain boundaries and bulk of the electrolyte phase with the former mechanism being the dominating one.Interestingly, it seems in this case the cation diffusion in the electrolyte pellet so strongly depends on the presence of NiO that near the interface, where the concentration of Ni is relatively lower than that in the bulk of the electrolyte pellet due to nickel diffusion into the double perovskite pellet, a layer depleted in Co and Pr is formed.The maximum of Co concentration is, thus, shifted deep into the electrolyte pellet at a distance of 400-540 µm from the interface with double perovskite.In this area, two types of Co-containing compounds were observed.The one with an approximate composition Y 0.70 Pr 0.24 Ba 0.94 Co 1.76 Ni 0.12 O 6−δ was detected on the grain boundaries (Figure 4b).The other is a mixed cobalt-nickel oxide with Co to Ni ratio 2.04:0.9,which forms particles with a size of about 2-5 µm, as seen in Figure 4b.The composition of the bulk of the electrolyte's grains in the area of maximum cobalt content corresponds on average to Ba 2 O 6−δ (R = Gd, Pr) and BaZr 0.8 Y 0.2 O 3−δ was studied by annealing the corresponding powder mixtures (50:50 wt.%) at different temperatures in the range 1100-1250 • C for 12 h in air with subsequent phase identification by the XRD.Ceramic pellets (Ø10 mm × 3 mm) for cation interdiffusion experiments were prepared by dry pressing at 150 MPa and sintering for 12 h in air at 1240, 1250 and 1550 • C for PrBaCo 2 O 6−δ , GdBaCo 2 O 6−δ and BaZr 0.8 Y 0.2 O 3−δ +NiO (0, 1 wt.%), respectively.Density of the ceramic samples was estimated by the Archimedes method using water as a medium.The relative densities of the as-prepared pellets were 95% for GdBaCo 2 O 6−δ , 97% for PrBaCo 2 O 6−δ , 92% for BaZr 0.8 Y 0.2 O 3−δ and 95% for BaZr 0.8 Y 0.2 O 3−δ +NiO (1 wt.%).One surface of each pellet was polished down to 1 µm using diamond suspension.Then it was purified ultrasonically by means of an IL 10 ultrasonic generator (INLAB, Saint Petersburg, Russia) in an ethanol medium and dried at 70 • C in a drying oven in air atmosphere.Cation interdiffusion was studied by the diffusion couple method.The following diffusion couples were studied: GdBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ , GdBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ +NiO (1 wt.%), PrBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ and PrBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ + . For example, the intensities of the peaks corresponding to double perovskites RBaCo 2 O 6−δ (R = Gd, Pr) strongly decrease with calcination temperature.At the same time, XRD peaks of BaZr 0.8 Y 0.2 O 3−δ slightly increase in intensity and shift towards high 2Θ • angles.Similar changes can be observed in the XRD patterns of the mixtures containing NiO; see Figure S2.All these evidences undoubtedly indicate in favor of chemical interaction between double perovskites and a solid electrolyte which occurs already at 1100 • C and likely at lower temperatures.The XRD results seem to suggest at least partial dissolution of RBaCo 2 O 6−δ (R = Gd, Pr) in BaZr 0.8 Y 0.2 O 3−δ oxide.

Figure 1 .
Figure 1.XRD patterns of the powder mixtures RBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ (50:50 wt.%) annealed at different temperatures for 12 h in air: (a,c) R = Gd; (b,d) R = Pr.The analysis of the cross section of the GdBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ diffusion couple annealed at 1200 • C for 48 h in air, as seen in Figure 3a, reveals the formation of Gd 2 O 3 on the surface of the GdBaCo 2 O 6−δ pellet along the interface of the diffusion couple.This is obviously in agreement with the XRD results discussed above.In addition, slight interdiffusion (depth ~15 µm) of Gd 3+ ions in the BaZr 0.8 Y 0.2 O 3−δ electrolyte and Y 3+ into the GdBaCo 2 O 6−δ phase was detected.The diffusion of Co ions occurs mainly through the bulk of the electrolyte.The diffusion depth reaches 300 µm.Cobalt concentration in the electrolyte pellet is practically constant at a distance of up to 30 µm from the interface with GdBaCo 2 O 6−δ .The average chemical composition of this electrolyte layer was estimated as Ba 1.02 Zr 0.76 Y 0.16 Co 0.06 O 3−δ .In the case of the PrBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ couple annealed at 1200 • C for 48 h in air, significant diffusion of both Co and Pr into the BaZr 0.8 Y 0.2 O 3−δ pellet was found, as shown in Figure 4a.Their diffusion occurs both along the grain boundaries and through the bulk of BaZr 0.8 Y 0.2 O 3−δ .The former mechanism is dominating.The diffusion depth of Co reaches 560 µm.The analysis of the grain boundary phase's chemical composition gives near the interface PrBaCo 2 O 6−δ |BaZr 0.8 Y 0.2 O 3−δ the following averaged composition: Ba 1.02 Co 0.72 Y 0.21 Pr 0.09 O 3−δ .However, the composition of the grains' interior may be expressed as Ba 0.99 Zr 0.69 Y 0.16 Co 0.13 Pr 0.03 O 3−δ .

Figure 2 .
Figure 2. The unit cell parameters vs. annealing temperature of powder mixtures RBaCo 2 O 6−δ (R = Gd, Pr)-BaZr 0.8 Y 0.2 O 3−δ (50:50 wt.%):(a) GdBaCo 2 O 6−δ , (b) PrBaCo 2 O 6−δ and (c) BaZr 0.8 Y 0.2 O 3−δ .The addition of NiO as a sintering aid to the BaZr 0.8 Y 0.2 O 3−δ solid electrolyte strongly promotes interdiffusion especially along the grain boundaries, as can be observed in Figures 3b, 4b and S3.As a result, for example, the depth of Co diffusion from PrBaCo 2 O 6−δ to BaZr 0.8 Y 0.2 O 3−δ is almost doubled, up to 1000 µm when 1 wt.% of NiO is added to the electrolyte.The diffusion couple GdBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ +NiO (1 wt.%) annealed at 1200 • C for 48 h in air does not show the formation of Gd 2 O 3 as a reaction product, contrary to the couple without NiO.Instead, diffusion of Gd 3+ into the electrolyte phase occurs mainly through the grain boundary mechanism, as seen in Figure 3b.Similar behavior is evident for diffusion of Co ions.As a result, the formation of Ba 1.03 Co 0.81 Ni 0.03 Gd 0.13 O 3−δ oxide was detected on the grain boundaries of the electrolyte pellet near the interface with double perovskite, whereas the composition of the electrolyte grains' bulk corresponds to Ba 1.02 Zr 0.78 Y 0.16 Co 0.04 O 3−δ .In addition, the diffusion of Ni 2+ ions from the electrolyte pellet to the GdBaCo 2 O 6−δ phase was observed.

Figure 3 .
Figure 3. SEM images and concentration maps of reaction zones for diffusion couples after annealing at 1200 • C for 48 h in air: (a) GdBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ and (b) GdBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ +NiO (1 wt.%); BSE (back scattered electrons) image and element maps were obtained near the phase boundary of the diffusion couple.The red arrow marks the interface between the pellets.

Figure 4 .
Figure 4. SEM images and concentration maps of reaction zones for diffusion couples after annealing at 1200 • C for 48 h in air: (a) PrBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ ; BSE image and element maps obtained near the phase boundary of diffusion couple and (b) PrBaCo 2 O 6−δ -BaZr 0.8 Y 0.2 O 3−δ +NiO (1 wt.%); BSE image and element maps obtained in the area of Co diffusion maximum in the depth of about 500 µm from the phase boundary.The corresponding images near the phase boundary of the diffusion couple may be observed in Figure S3.The red arrow marks the interface between the pellets.