Crystal Structure and Transport Properties of Oxygen-Deficient

Jan 29, 2018 - Department of Mechanical Engineering, University of South Carolina, Columbia, South Carolina 29208, United States. ‡ Department of ...
0 downloads 0 Views 4MB Size
Article www.acsaem.org

Cite This: ACS Appl. Energy Mater. 2018, 1, 822−832

Crystal Structure and Transport Properties of Oxygen-Deficient Perovskite Sr0.9Y0.1CoO3−δ Tianrang Yang,† Victoria F. Mattick,‡ Yan Chen,§ Ke An,§ Dong Ma,§ and Kevin Huang*,† †

Department of Mechanical Engineering, University of South Carolina, Columbia, South Carolina 29208, United States Department of Chemical Engineering, University of South Carolina, Columbia, South Carolina 29208, United States § Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States ‡

S Supporting Information *

ABSTRACT: The present work reports a systematic study on temperature-dependent local crystal structure, oxygen stoichiometry, and electrical/electrochemical properties of an oxygendeficient Sr0.9Y0.1CoO3−δ (SYC10) perovskite using variabletemperature neutron diffraction (VTND), thermal gravimetric analysis, and electrical/electrochemical methods, respectively. The VTND reveals that the crystal symmetry of SYC10 remains P4/mmm tetragonal up to 900 °C. The tetragonal symmetry reflects the net effects of temperature and oxygen stoichiometry on crystal symmetry. The observed p-type electronic conductivity behavior originates from the charge-ordering between the two distinctive Co-sites. The partial oxide-ion conductivity and diffusivity obtained from oxygen permeation measurements are 2.3 × 10−2 S cm−1 and 7.98 × 10−8 cm2/s at 800 °C in air, respectively. The electrochemical oxygen reduction reaction kinetics of the SYC10 cathode is primarily limited by the charge-transfer process at low temperatures (600−650 °C) and oxide-ion migration from the cathode into the electrolyte at high temperatures (700−800 °C). KEYWORDS: perovskite, oxygen stoichiometry, charge-ordering, electronic conductivity, oxide-ion conductivity

I. INTRODUCTION A key to the success of intermediate-temperature solid oxide fuel cells (IT-SOFCs) is to discover electrocatalytically active oxides to promote the limiting oxygen reduction reaction (ORR) at the cathode. It has been demonstrated in the past that perovskiterelated and oxygen-deficient complex oxides are a promising class of cathode materials for IT-SOFCs.1−8 Among them, BaCoO3−δ (BCO) and SrCoO3−δ (SCO) based simple perovskite oxides attract the most interest because of their high intrinsic catalytic activity with regard to ORR. However, due to their overly high Goldschmidt tolerance factor (>1.0), pure BCO and SCO are not stable under high-temperature and ambient-pressure conditions, making their use impractical. Previous studies have shown that a proper doping of the A- and/or B-cation in BCO and SCO can stabilize the perovskite structure when subjected to elevated temperatures, which has become a major driver for active studies on these two classes of materials in recent years.9−20 Another challenge for BCO- and SCO-based perovskites to be used as practical cathodes is their much higher thermal expansion coefficient (TEC; >20 ppm K−1) when compared to the electrolyte (∼10 ppm K−1). A direct use of these two materials in SOFC as a bulk cathode is, therefore, inadequate. To solve this issue, the ORR-active, but TEC-high, BCO- and SCO-based materials are usually utilized as a cathode in the form of nanoparticles (NPs) impregnated into a TEC-compatible scaffold.21−23 In this design, the TEC of the cathode is determined by the scaffold, not the attached NPs. Another critical issue © 2018 American Chemical Society

particularly associated with BCO-based perovskites is their high chemical propensity to react with CO2 in air, resulting in a gradual formation of BaCO3 and long-term performance decay. In comparison, SCO-based perovskites are less reactive with CO2; thus they have a better potential to be considered a practical IT-SOFC cathode.24,25 In recent years, there have been many reports on SCO-based oxygen-deficient perovskites in the literature.13−20,26−28 A general conclusion from these prior studies is that the thermal stability of these materials can be significantly enhanced through donor-doping (dopant with higher oxidation state than the host cation) on either A- and B-cation. Among them, Sr0.9Y0.1CoO3‑δ (SYC10)26,27 has been reported as an IT-SOFC cathode material with low polarization resistance (Rp).27 From the roomtemperature structural analysis using X-ray diffraction (XRD), both tetragonal (P4/mmm, No. 123)26 and cubic (Pm3̅m, No. 221) structures have been observed to exist in SYC10.27,28 Since SYC10 operates at elevated temperatures when employed as a SOFC cathode, the knowledge of the local structure and oxygen stoichiometry at elevated temperatures would be more meaningful than that under ambient temperature to interpret the ORR activity. However, the conventional XRD techniques, whether at room temperature or at higher temperatures, cannot provide Received: December 14, 2017 Accepted: January 29, 2018 Published: January 29, 2018 822

DOI: 10.1021/acsaem.7b00275 ACS Appl. Energy Mater. 2018, 1, 822−832

Article

ACS Applied Energy Materials

used for VTND runs. With the initial RT oxygen stoichiometry determined from iodine titration or ND refinement, the oxygen stoichiometry as a function of temperature can then be mapped according to weight changes. 2.4. Iodometric Titration. In addition to ND refinement, the initial oxygen stoichiometry (3 − δ) at RT was also determined by iodometric titration. The procedure is described as follows. Ultrapure N2 is bubbled into 30 mL of diluted 1 M HCl solution for 1 h to eliminate the dissolved O2 and Cl2. Then 50 mg of powder was dissolved in the above-mentioned solution under a N2 atmosphere, followed by adding 10 mL of 0.2 M KI solution to reduce all Co3+/4+ to Co2+. Next, 0.01 M Na2S2O3 is added dropwise to titrate the released elementary iodine by the reduction reaction. A few drops of saturated starch solution was used to indicate the end point of titration. With a known Na2S2O3 concentration, the average Co-ion oxidation state of the sample can then be calculated based on the following charge balance: n Co(2 + n) + + nI− = Co2 + + I 2 (1) 2

accurate information on the oxygen positions/occupancies due to a much weaker scattering contribution of oxygen to the total Bragg scattering than that of cations. In this study, we report a systematic study on temperaturedependent local crystal structure and oxygen occupancy (stoichiometry) in SYC10 determined by variable-temperature neutron diffraction (VTND). To complement VTND results, conventional thermogravimetric analysis (TGA) and iodometric titration for oxygen nonstoichiometry, oxygen permeation for partial oxideion conductivity, and impedance spectroscopy for polarization resistance are also carried out. With these fundamental data in hand, we attempt to correlate the crystal structure and oxygen nonstoichiometry with transport properties such as partial oxideion conductivity, oxygen-related diffusivity, and polarization resistance to further elucidate the ORR mechanism in SYC10.

II. EXPERIMENTAL PROCEDURE 2.1. Sample Synthesis. The SYC10 was prepared by a conventional solid-state reaction method. Stoichiometric amounts of SrCO3 (≥99.9%, Aldrich), Co3O4 (99.7%, Alfa Aesar), and Y2O3 (99.9%, Alfa Aesar) were first ball-milled in ethanol for 3 h. After drying, the powder was then pressed into pellets, followed by calcination at 1000 °C for 12 h with a 3 °C min−1 heating and cooling rate. The pellets were then broken up and ball-milled for 3 h, and calcined again at 1050 °C for another 12 h. Another 3 h ball-milling was followed before making the final pellets or rectangular bars. The final pellets, or rectangular bars, embedded in SYC10 powder were sintered at 1230 °C for 10 h to ensure phase purity and homogeneity. 2.2. Variable-Temperature Neutron Diffraction. Variabletemperature neutron diffraction experiments were conducted on the sintered pellets (∼17 mm in diameter and ∼10 mm in length) situated in an MgO crucible using VULCAN, the time-of-flight Engineering Diffractometer at the Spallation Neutron Source in Oak Ridge National Laboratory.29,30 The temperature was varied from the base furnace temperature (70 °C) to 900 °C, while the atmosphere remained in ambient air. The VTND data were collected continuously during both heating and cooling with a heating rate of 10 °C min−1 below 300 and 5 °C min−1 between 300 and 900 °C and a constant cooling rate of 10 °C min−1. During heating, the samples were held at 300 °C for 3 h and for 2 h at three other temperatures, i.e., 500, 700, and 900 °C. Only data collected during the last 1 h with total neutron counts of 4 × 106 at each temperature were used for the structural refinement. The beam gauge volume was determined to be 5 × 10 × 5 mm3 by the incident slits and receiving collimators. Room-temperture (RT) patterns were collected ex situ from a 5 g powder sample. The Rietveld refinements were performed with a GSAS program using the EXPGUI interface.31,32 Typically, a refinement starts with trials and errors using the model-free Le Bail method.33 The background and peak profiles were modeled with the shifted Chebyshev polynomial function and pseudo-Voigt function with Finger−Cox− Jephcoat asymmetry correction, respectively. The DIFA (one of the three parameters defining the relationship between measured time-offlight and d-spacing),31 unit cell, background, and profile parameters from Le Bail fitting were used as the initial inputs into the Rietveld refinement and kept fixed during the first stage. Scaling factors, atomic positions, isotropic thermal displacement parameters, and absorption coefficients were then refined separately. Following that, all the fixed parameters were released and refined together. Anisotropic thermal displacement factors were refined depending upon allowance by the data. Oxygen occupancies were the outcomes of the final run. The DIFA and absorption coefficients of in situ ND were considered temperature-independent. The crystal structure was constructed with Vesta 3.0.34 2.3. Thermogravimetric Analysis. TGA was performed with a NETZSCH STA 448 TGA/DSC thermal analyzer. Approximately 50 mg of powder was placed in an alumina crucible under a flow of 60 mL min−1 air. The sample was held at a reference state of 70 °C for 2 h to remove the absorbed H2O, followed by the same temperature profiles

I 2 + 2S2 O32 − = 2I− + S4 O6 2 −

(2)

To ensure accuracy, the concentration of Na2S2O3 solution needs to be calibrated. To do this, 25 mg of K2Cr2O7 powder is first dissolved in the 1 M HCl solution mixed with an excessive 0.2 M KI solution. Cr2O72− is expected to be fully reduced to Cr3+ by KI. Following the same procedure described above, the produced elementary iodine is titrated using the Na2S2O3 solution. The concentration of Na2S2O3 can then be determined from the charge balance between reactions 2 and 3. Cr2O7 2 − + 6I− + 14H+ = 2Cr 3 + + 3I 2 + 7H 2O

(3)

2.5. Electrical Conductivity Measurement. Bar samples (5 mm × 2 mm × 15 mm) were used for total conductivity (σ) measurements by the standard four-probe method. Both T- and PO2-dependent σ were measured from 800 to 500 °C with an interval of 100 °C and PO2 of 1, 0.5, 0.2, 0.1, and 0.05 atm. Since it takes time for O2 to fully equilibrate with SYC10, we adopted the following wait-time for each temperature to ensure the accuracy of the conductivity: 3 h at 800 °C, 5 h at 700 °C, 10 h at 600 °C, and 18 h at 500 °C. Clearly, a lower temperature requires longer wait-time to reach a true gas/solid equilibrium, which also infers that this material has a low oxygen surface exchange rate at low temperatures. The resistance data were collected in the form of V−I curve with the E-I testing module available in the CorrWare program operated on a Solartron 1470 multichannel potentiostat. 2.6. Partial Oxide-Ion Conductivity Determination. The partial oxide-ion conductivity in SYC10 was determined by the classical Wagner equation using the measured oxygen permeation flux data:

JO = − 2

RT 16F 2L

ln POpermeate 2

∫ln P

feeding O2

t O2 −teσtotal d[ln PO2]

(4) −2 −1

where JO2 denotes the oxygen permeation flux (mol cm s ); R is gas constant 8.314 J mol−1 K−1; F is Faraday constant 96,485 C mol−1; L is membrane thickness (cm); and tO2− and te are the ionic and electronic transport numbers, respectively. For a MIEC with a dominant σe, te ≈ 1, eq 4 is then simplified into JO = − 2

RT 16F 2L

ln POpermeate 2

∫ln P

feeding O2

σ O2 − d[ln PO2]

(5)

Since σO2− varies with PO2 by a power law relationship:

σ O2 − = σ O°2 −POn 2

(6)

substitution of eq 6 into eq 5 followed by integration gives

JO = 2

823

A [(POpermeate )n − (POfeeding )n ] 2 2 n

(7) DOI: 10.1021/acsaem.7b00275 ACS Appl. Energy Mater. 2018, 1, 822−832

Article

ACS Applied Energy Materials A=−

σ Oo2 −RT 2

16F L

of 0.01−1 atm by electrochemical impedance spectroscopy (EIS). The EIS spectra were collected within a frequency range of 0.01 Hz to 1 MHz and AC stimulus amplitude of 10 mV using a Solartron 1470/1455B electrochemical station. The cell consisted of an electrolyte GDC (Gd0.2Ce0.8O2−δ, Fuel Cell Materials) sandwiched between two identical screen-printed SYC10 cathodes. The effective electrode area was 0.712 cm2. Silver paste and silver mesh were used as the current collector. The obtained EIS spectra were analyzed by an equivalent circuit model using ZSimpWin Demo software.

(8)

Here A and n values are constants that can be obtained by fitting at a fixed Pfeeding using eq 7. The partial experimental JO2 vs Ppermeate O2 O2 2− oxygen-ion conductivity σO can then be calculated from eq 6 with the σ°O2− value obtained from eq 8. The actual permeation cell setup schematic is given in Figure S1 of the Supporting Information (SI). Briefly, a dense SYC10 membrane with a thickness of 1.15 mm was first sealed onto a supporting alumina tube by a commercial silver paste (Shanghai Research Institute of Synthetic Resins). Air flowing at 100 cm3(STP) min−1 was fed to one surface of the membrane, while the other surface was swept with ultrapure Ar at variable flow rates. To vary PO2 at the permeate side, the sweeping Ar flow rate was changed from 5 to 200 cm3(STP) min−1. All gas flow rates were controlled by mass flow controllers (Smart-Trak 50 series). Since the cell needs to be first sealed at 800 °C, all measurements were made from 800 °C down to 650 at 25 °C interval. The equilibration time at a specific temperature varied from 2 h at 800 °C to 12 h at 650 °C. The concentrations of O2 and N2 in the effluent were analyzed by an online gas chromatography (Agilent 490). The final flux densities of oxygen (JO2) were calculated by the following equation:

JO = 2

CO2 Q 1 − CO2 − C N2 A

III. RESULTS AND DISCUSSION 3.1. Structural Determination. The temperature-dependent d-spacing evolution of SYC10 obtained from VTND is shown in Figure 1. All the diffraction peaks remain the same during the temperature excursion from RT to 900 °C, except shifts in peak positions which resulted from lattice expansion/shrinkage during heating/cooling cycles, respectively. No phase transition is found by VTND within the temperature range studied.

(9)

where CO2 and CN2 are the measured concentrations of O2 and N2, respectively; Q is the flow rate of sweeping Ar gas (cm3(STP) min−1); and A is the effective area of the sample (0.745 cm2). Any detectable N2 is corrected in eq 9. The detected N2 flux is very small, typically 500 °C do O2-sites show vacancies. In comparison, major oxygen vacancies (V•• o ) are populated at O3-site within the temperature range measured. A further plot of c-axis coordination of 2h Sr (ideally 0.25) and 2g O (ideally 0.75) shown in Figure 4b reveals variations of the structural distortion with temperature. Overall, the structural distortion is resulted from two competing effects: higher temperature promoting symmetry and oxygen stoichiometry variations promoting distortion. 3.2. Lattice Constants vs Temperature. The lattice constants a and c of the SYC10s tetragonal structure are plotted against temperature in Figure 5a,b, respectively. It appears that there is a deviation from the linear trending at 500 °C, above which a greater increase of a and c with temperature is observed. We believe that this increased lattice expansion reflects an additional contribution from chemical oxygen loss at the O2-site above 500 °C as shown in Figure 4a. Since O2 is the bridging oxygen between the two distinct CoO6 octahedra layers, its deficiency could create dimensional changes along both a and c

The Rietveld refinement profiles for ND patterns collected at RT and a representative 500 °C are shown in Figure 2 as examples. The collected profiles show marginal differences with those calculated from crystal structure model based on space group P4/mmm (a × a × 2a unit cell). In addition to nuclear Bragg reflections, magnetic peaks are also observed at RT as shown in Figure 2a. Co-ions in SYC10 are G-type magnetically ordered with a propagation vector k = (1/2, 1/2, 0),35 which is the same as SrCo0.9Nb0.1O3‑δ. No magnetic peaks were observed in VTND patterns above 300 °C, suggesting that the Néel temperature of SYC10 is lower than 300 °C. Indeed, TN = 264 °C has been previously reported for this class of materials.36 The crystallographic space group P4/mmm (a × a × 2a unit cell) best describes the nuclear structure of SYC10; it contains atomic positions of Sr at 2h (1/2, 1/2, ∼0.25), Co1/Nb1 at 1a (0, 0, 0), Co2/Nb2 at 1b (0, 0, 1/2), O1 at 2f (1/2, 0, 0), O2 at 2g (0, 0, ∼0.75), and O3 at 2e (1/2, 0, 1/2). Structurally speaking, Co1 forms octahedral coordination with O1 and O2, while Co2 coordinates with O2 and O3. A unit cell structure is depicted in Figure 3, and the refined parameters are given in Table 1 with reliability factors of the refinement. The tetragonal structure of SYC10 can be viewed as a distorted 3C-perovskite cubic structure (Pm3̅m) with oxygen-ordering. Figure 4a shows that the O1-site occupancies increase from 825

DOI: 10.1021/acsaem.7b00275 ACS Appl. Energy Mater. 2018, 1, 822−832

Article

ACS Applied Energy Materials

Figure 4. (a) Oxygen occupancies variation with temperature; (b) c-axis coordination of Sr and O2 determined from VTND. The deviation of Sr and O2 from their ideal position (Sr, 1/2, 1/2, 1/4; O2, 0, 0, 3/4) indicates the distortion level.

Figure 5. Variations of lattice parameters and thermal expansion coefficients with temperature in SYC10: (a and c) along a-axis; (b and d) along c-axis.

VTND. The two curves generally follow the same trend, but the deviation between the two methods becomes larger above 700 °C. A possible reason for the enlarged deviation is the increased uncertainties due to the increased background noise in VTND refinement. It is interesting to note that both TGA and VTND confirm an oxygen gain occurring from RT to 300 °C. We believe that this oxygen gain is resulted from gaseous O2 occupying the “frozen oxygen vacancies” trapped by the slow SYC10/O2 exchange kinetics during cooling; the compacted samples used for enhancing cation diffusion and phase formation could be the reason for the latter. 3.4. Electrical Conductivity vs Temperature and Partial Pressure of Oxygen. The variations of electrical conductivity of SYC10 with T and PO2 are shown in Figure 7. In general, at a given temperature, the conductivity increases with PO2, suggesting

axes as observed in Figure 5a,b. The TGA curve shown in Figure 6a indeed indicates a greater weight loss above this temperature. In addition, the decrease of the average Co-ions oxidation state as shown in Figure 6b also contributes to the curvature of lattice constants a and c vs T. The fitted linear thermal expansion coefficients (LTECs, α) along a-axis (Figure 5c) and c-axis (Figure 5d) are listed in Table 2 for both low- and high-temperature regions. The average α is calculated by (2αa + αc)/3. 3.3. Oxygen Stoichiometry vs Temperature. The oxygen stoichiometry (3 − δ) at RT determined by iodometric titration is 2.580, a value which matches well with 2.543 obtained from VTND’s Rietveld refinement. The temperature dependence of 3 − δ with the initial value from iodometric titration can be further calculated from the TGA results of Figure 6a; the results are shown in Figure 6b and compared with those obtained from 826

DOI: 10.1021/acsaem.7b00275 ACS Appl. Energy Mater. 2018, 1, 822−832

Article

ACS Applied Energy Materials

Figure 6. (a) TGA curves measured in air; (b) oxygen stoichiometry (3 − δ) and Co-ion oxidation state from both ND and TGA. The initial value for TGA-derived (3 − δ) was taken from iodometric titration.

Table 2. Fitted Linear Thermal Expansion Coefficients α along a-Axis and c-Axis α×10−6/K−1

RT−500 °C

500−900 °C

a-axis c-axis av

19.5 19.5 19.5

34.8 22.1 30.6

where j is the surrounding cation; rij is the observed bond lengths between i and j; r0 is the ideal bond length (1.70 Å for Co3+−O38); and b is an empirical constant (0.37 Å). The rij values are taken from VTND refinement results listed in Table S1. The calculated BVS results are shown in Table 3, where a remarkable Table 3. Bond Valence Sums of Co-Ions

a p-type electron−hole conduction mechanism. However, at a given PO2, the conductivity decreases with temperature. This infers that increasing temperature substantially decreases hole concentration (or the oxidation state of Co-ions). Clearly, oxygen loss is one major reason for the lowering of hole concentration, and thus electronic conductivity. The average Co-ion oxidation state obtained from TGA and VTND as shown in Figure 6b is generally lower than +3.0, particularly between 500 and 900 °C. What causes the observed p-type electronic conduction behavior at such a low oxidation state of Co-ion is an intriguing scientific question that should be addressed. With the VTND data, we are able to perform bond valence sums (BVS) analysis for the oxidation state of Co-ions within their local coordinating environment. In theory, the BVS of a cation i is the sum of the individual bond valences vij. defined by37 ⎛ ro − rij ⎞ ⎟ vij = exp⎜ ⎝ b ⎠ (10)

Co-ion

RT

300 °C

500 °C

700 °C

900 °C

Co1 Co2

3.82 2.65

3.75 2.55

3.62 2.55

3.57 2.41

3.37 2.44

difference in BVS between Co1 and Co2 is observed from RT to 900 °C, suggesting a charge-ordering at these two Co-sites, i.e., a mixed oxidation state of Co3+/4+ on the Co1-site and Co2+/3+ on the Co2-site. The determined oxidation state of Coions from oxygen stoichiometry reflects the average values of Co1 andCo2. With the charge-ordering mechanism, it is possible to postulate that the SYC10 can exhibit either p-type electron−hole conduction via Co1 or n-type excess electron conduction via Co2, depending on actual T and PO2 conditions. In general, hopping of electrons between two neighboring Co-ions is mediated by the conduction band formed by the orbital overlap between Co 3d and O 2p commonly known as a double exchange mechanism.39 The stronger the overlap, the broader the band and therefore the easier for electrons to hop. However, the presence

Figure 7. (a) Temperature dependence of σ at different Po2; (b) Arrhenius plots of σ. 827

DOI: 10.1021/acsaem.7b00275 ACS Appl. Energy Mater. 2018, 1, 822−832

Article

ACS Applied Energy Materials of oxygen vacancies at O3-sites in coordination with Co2 (see Figure 3) distorts the crystal field and weakens the overlap between Co2 3d and O 2p, resulting in a narrower conduction band formation and promoting n-type electronic conduction with a higher energy barrier.40 Meanwhile, the positively charged donordopant YSr• can also trap negatively charged n-type electrons, making the total electronic conductivity in SYC10 dominated by p-type electron holes. It is to be noted that the portion of n-type electronic conduction will increase with a narrowed band, increased oxygen vacancy concentration, and activation energy as T increases and PO2 decreases. Figure 7b indeed shows such a trend. Finally, the level of activation energy shown in Figure 7b seems to suggest that electron−holes conduction in SYC10 gradually transitions from a delocalized to a localized state as PO2 is decreased at a fixed T. This should also be the case as T is increased at a given PO2. 3.5. Partial Oxide-Ion Conductivity and Oxygen SelfDiffusion Coefficient. Partial oxide-ion conductivity (σO2−) of a mixed electron and oxide-ion conductor (MIEC) plays a critical role in oxygen electrocatalysis for reversible SOFCs. Therefore, σO2− is often used to evaluate whether an oxide is a suitable cathode material for SOFCs. However, determination of σO2− is a challenging task due to the domination of electronic conductivity, σe, in these materials. To separate σO2− from σe, a number of experimental methods have been proposed in the past,41−44 among which oxygen permeation is deemed the most reliable and reproducible method.45 Therefore, in this study we adopt the oxygen permeation method to determine σO2−. Figure 8 shows the measured oxygen flux density of SYC10 vs PO2 at different temperatures, along with the modeled curves

Figure 9. Partial oxide-ion conductivity, σO2−, vs temperature at PO2 = 0.21 atm.

σ O2 − =

⎛ − 1.16 ± 0.05 eV ⎞ 6.28 × 106 ⎟ exp⎜ ⎝ ⎠ T RT

(11)

The activation energy Ea = 1.16 ± 0.05 eV is close to that reported for oxide-ion conductors with an oxygen vacancyhopping mechanism.55 The oxygen vacancy diffusion coefficient, Dv, can be further calculated using the Nernst−Einstein equation,56 in which the oxide-ion conductivity from oxygen permeation, oxygen nonstoichiometry from TGA, and unit cell volume from VTND refinement are used: σ O2 − = Cv =

4F 2CvDv RT

(12)

[V •• O] Vm

(13)

where Cv is the oxygen vacancy concentration; Vm is the molar volume of the oxide, and [V•• O ] equals oxygen nonstoichiometry δ assuming that all oxygen vacancies are mobile at elevated temperatures. The oxide-ion transport in perovskite oxides is generally considered to occur via a vacancy mechanism. Therefore, combining eqs 12 and 13 gives Dv =

σ O2 −RTVm 4F 2δ

(14)

The oxygen-ion self-diffusion coefficient, DO, is related to Dv by DO = Figure 8. Oxygen flux density vs Ppermeate from 650 to 800 °C. O2

Cv δ Dv = Dv CO 3−δ

(15)

where CO is the concentration of oxide-ion. Combining eqs 14 and 15 leads to

using eq 7. The detailed model values are listed in Table S2. The oxygen permeation flux is nearly 1−2 orders of magnitude lower than those of similar systems SrCo0.9Nb0.1O3−δ13,46−52,64 and SrCo0.9Ta0.1O3−δ19,53,54 under a similar testing condition. For example, at 800 °C, a 1 mm thick SrCo0.9Nb0.1O3−δ48 and 1.36 mm thick SrCo0.9Ta0.1O3−δ53 can yield oxygen flux densities of 1.5 × 10−6 and 6 × 10−753 mol cm−2 s−1, respectively, vs 8.9 × 10−8 mol cm−2 s−1 for a 1.15 mm thick SYC10 of this study. Clearly, the cubic structure of B-sitedoped SCO35,54 favors a higher partial oxide-ion conductivity than the tetragonal structure A-site-doped counterpart. The σO2− of SYC10 at PO2 = 0.21 atm obtained from the oxygen flux is shown as a function of temperature in Figure 9 and analytically expressed by

DO =

σ O2 −RTVm 4F 2

(3 − δ)

(16)

The calculations of Vm and δ can be found in the SI. The temperature-dependent Dv and DO, along with Cv and CO in a temperature range of 650 to 800 °C is shown in Figure 10. The decreases of CO with T arise from both oxygen loss and lattice expansion (larger molar volume). However, oxygen loss plays a more pronounced role since Cv increases monotonically with T. Both Dv and DO shown in Figure 11 follow the Arrhenius relationship with the following analytical expressions: ⎛ − 1.13 ± 0.05 eV ⎞ ⎟ Dv = 7.14 × 10−2 exp⎜ ⎝ ⎠ RT 828

(17)

DOI: 10.1021/acsaem.7b00275 ACS Appl. Energy Mater. 2018, 1, 822−832

Article

ACS Applied Energy Materials

Figure 10. (a) Dv, Cv and (b) DO, CO dependence on temperature.

between Rp and PO2 in an oxygen reduction reaction generally follows63,64 R p ∝ (PO2)−m

with m values giving the information on the rate-limiting steps: (1) m = 1, O2(g) → O2,adsorption, oxygen gas adsorbs on the cathode surface (2) m = 0.5, O2,adsorption → 2Oadsorption, adsorbed oxygen gas dissociates into oxygen atoms × (3) m = 0.25, Oadsorption + 2e′ + V•• O → OO, adsorbed oxygen atoms are reduced into oxide-ion by electrons and vacancies into lattice oxygen, which is also known as the charge-transfer process (4) m = 0, Oelectrode2− → Oelectrolyte2−, oxide-ion migrates from cathode into electrolyte The Nyquist plots measured at 650 °C from PO2 = 0.01 to 1 atm are shown in Figure 12 as an example; the plots at other

Figure 11. Arrhenius plot of Dv and DO.

⎛ −1.18 ± 0.05 eV ⎞ ⎟ DO = 2.42 × 10−2 exp⎜ ⎝ ⎠ RT

(19)

(18)

The activation energy Ea for DO is 1.18 eV, which follows the trend suggested by previous reports that Ea varies within 3.2−1.5 eV with x varying from 0 to 0.5, respectively, for La1−xSrxCoO3−δ systems.57,58 The DO values at 800 °C obtained in this work are also compared with others in Table 4. For a similar A-siteTable 4. Comparison of DO of This Work with Other Perovskite Cathodes in the Literature composition

DO (cm2 s−1) at 800 °C

ref

La0.35Sr0.65Co0.7Fe0.3O3‑δ (LSCF) SrFe0.33Co0.67O3‑δ (SFC) Ba0.5Sr0.5Co0.8Fe0.2O3‑δ (BSCF) Sr0.9Ce0.1CoO3‑δ (SCC) Sr0.9Y0.1CoO3−δ (SYC10)

9.0 × 10−7 5.2 × 10−6 2.98 × 10−6 6.0 × 10−8 7.98 × 10−8

59 60 61 62 this work

Figure 12. Impedance spectra of SYC10 symmetrical cell under different PO2 at 650 °C.

doping sample Sr0.9Ce0.1CoO3−δ, the obtained DO by this work is comparable, e.g., 7.98 × 10−8 vs 6 × 10−8 cm2 s−1. However, it is 1−2 orders of magnitude lower than other cubic structured perovskites listed. This observation also supports the later conclusion (section 3.6) that the rate-limiting step for ORR at SYC10 is oxide-ion migration from cathode into electrolyte at 700−800 °C because of relatively lower oxide-ion conductivity and diffusivity. 3.6. Cathode Polarization Behavior. The PO2 dependence of the area-specific polarization resistance (Rp) was measured to understand the polarization mechanisms. The relation

temperatures are similar but with a different resistance range. The plots consist of two semicircles, one at high frequency (Rp,high) and another at low frequency (Rp,low). Figure 13a shows a very weak dependence of Rp,high on PO2, i.e., m = 0.0008−0.03, suggesting the high-frequency arcs correspond well to step 4. For the low-frequency arc, Figure 13b, the m values increase with temperature from 0.36 at 600 °C to 0.97 at 800 °C, suggesting that the rate-limiting step changes from chargetransfer (step 3) to surface dissociation/adsorption (steps 1 and 2). To calculate the activation energy of SYC10 as a cathode, the EIS of symmetrical cells was measured as a function of tem829

DOI: 10.1021/acsaem.7b00275 ACS Appl. Energy Mater. 2018, 1, 822−832

Article

ACS Applied Energy Materials

Figure 13. PO2 dependence of polarization resistances: (a) high-frequency resistance and (b) low-frequency resistance.

The average oxidation state of Co-ions in SYC10 is very close to or smaller than +3.0, while p-type conduction is observed in a temperature range of 500−800 °C and PO2 = 0.05−1 atm. The BVS calculations indicate a mixed oxidation state of Co3+/4+ on Co1-site and Co2+/3+ on Co2-site. Hence, the charge-ordering could render SYC10 to exhibit either p- or n-type electronic conduction, depending on the actual T and Po2 condition. The partial oxide-ion conductivity and diffusivity of the tetragonal structured SYC10 measured from oxygen permeation are 2.3 × 10−2 S cm−1 and 7.98 × 10−8 cm2/s at 800 °C, respectively, which is 1−2 orders of magnitude lower than cubic structured LSCF, SFC, and BSCF. The lower oxide-ion conductivity appears to be the reason for the higher Rp exhibited by SYC10 in comparison to Nb- and Ta-doped SCO. The PO2 dependence of polarization resistances shows that ORR kinetics is primarily limited by surface-controlled charge-transfer process at low temperatures (600−650 °C) and oxide-ion migration from cathode into electrolyte at high temperatures (700−800 °C).

Figure 14. Arrhenius plots of Rp, Rp,high, and Rp,low at PO2 = 0.21 atm from 600 to 800 °C.

perature in air; the results are shown in Figure 14. The Rp is 0.24 Ω cm2 at 700 °C, which is almost twice the reported values of SrCo0.9Nb0.1O3‑δ and SrCo0.9Ta0.1O3‑δ.54 This observation appears to be consistent with the lower oxygen permeation flux of SYC10 than those B-site-doped SCO counterparts mentioned above. It also suggests that partial oxide-ion conductivity indeed plays an important role in determining RP. The low-frequency process has a higher activation energy than the high-frequency one (2.33 eV vs 1.16 eV), which suggests that the former may be related to a surface-limiting process. The overall Rp is mostly dominated by Rp,low at 600− 650 °C, but by Rp,high at 700−800 °C. Based on the mechanisms outlined above, it appears that, at low temperatures, ORR kinetics is limited by a slower charge-transfer (dissociative adsorption and electron transfer) process and at high temperatures the limiting ORR step transitions to oxide-ion transfer across the cathode/electrolyte interface.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsaem.7b00275. Oxygen permeation flux test setup schematic and additional data (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Kevin Huang: 0000-0002-1232-4593 Notes

The authors declare no competing financial interest.



IV. CONCLUSIONS In summary, the crystal structure of Sr0.9Y0.1CoO3−δ (SYC10) is characterized by variable-temperature neutron diffraction to be P4/mmm tetragonal from RT to 900 °C. The degree of tetragonal distortion is affected by two competing processes: high temperatures promoting a more symmetric structure and oxygen intake/loss promoting lattice distortion. The SYC10 structure contains an alternately arranged oxygen-deficient Co2 layers coordinating with O2 and O3, where V•• O exists and randomly distributes, and an almost oxygen-saturated Co1 layer coordinating with O1 and O2.

ACKNOWLEDGMENTS This work was funded by the Advanced Research Projects Agency-Energy (ARPA-E), U.S. Department of Energy, under Award No. DE-AR0000492; the Office of Fossil Energy, U.S. Department of Energy, under Award No. DE-FE-0023317; and the National Science Foundation, under Award No. CBET1464112. Neutron scattering was carried out at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory, which is one of the user facilities sponsored by the Scientific User Facilities Division, Office of Basic Energy Sciences, U.S. 830

DOI: 10.1021/acsaem.7b00275 ACS Appl. Energy Mater. 2018, 1, 822−832

Article

ACS Applied Energy Materials

(15) Cascos, V.; Troncoso, L.; Alonso, J. New Families of Mn+-doped SrCo1‑xMxO3‑δ Perovskites Performing as Cathodes in Solid-oxide Fuel Cells. Int. J. Hydrogen Energy 2015, 40, 11333−11341. (16) Baszczuk, A.; Dabrowski, B.; Kolesnik, S.; Chmaissem, O.; Avdeev, M. Structural and Physical Properties of Re-Substituted B-site Ordered and Disordered SrCo1‑xRexO3‑δ (x= 0.1, 0.25, 0.5). J. Solid State Chem. 2012, 186, 240−246. (17) Yang, W.; Zhang, H.; Sun, C.; Liu, L.; Alonso, J.; FernándezDíaz, M.; Chen, L. Insight into the Structure and Functional Application of the Sr0.95Ce0.05CoO3‑δ Cathode for Solid Oxide Fuel Cells. Inorg. Chem. 2015, 54, 3477−3484. (18) Aguadero, A.; Pérez-Coll, D.; Alonso, J.; Skinner, S.; Kilner, J. A New Family of Mo-doped SrCoO3‑δ Perovskites for Application in Reversible Solid State Electrochemical Cells. Chem. Mater. 2012, 24, 2655−2663. (19) Li, M.; Zhao, M.; Li, F.; Zhou, W.; Peterson, V. K.; Xu, X.; Shao, Z.; Gentle, I.; Zhu, Z. A Niobium and Tantalum Co-Doped Perovskite Cathode For Solid Oxide Fuel Cells Operating Below 500 °C. Nat. Commun. 2017, 8, 13990. (20) Gwon, O.; Yoo, S.; Shin, J.; Kim, G. Optimization of La1‑xSrxCoO3‑δ Perovskite Cathodes for Intermediate Temperature Solid Oxide Fuel Cells Through the Analysis of Crystal Structure and Electrical Properties. Int. J. Hydrogen Energy 2014, 39, 20806−20811. (21) Ai, N.; Jiang, S. P.; Lü, Z.; Chen, K.; Su, W. Nanostructured (Ba,Sr)(Co,Fe)O3‑δ Impregnated (La,Sr)MnO3 Cathode for Intermediate-temperature Solid Oxide Fuel Cells. J. Electrochem. Soc. 2010, 157, B1033−B1039. (22) Ren, Y.; Cheng, Y.; Gorte, R. J.; Huang, K. Toward Stabilizing Co3O4 Nanoparticles as an Oxygen Reduction Reaction Catalyst for Intermediate-Temperature SOFCs. J. Electrochem. Soc. 2017, 164, F3001−F3007. (23) Ding, D.; Li, X.; Lai, S. Y.; Gerdes, K.; Liu, M. Enhancing SOFC Cathode Performance by Surface Modification Through Infiltration. Energy Environ. Sci. 2014, 7, 552−575. (24) Zhang, Y.; Yang, G.; Chen, G.; Ran, R.; Zhou, W.; Shao, Z. Evaluation of the CO2 Poisoning Effect on a Highly Active Cathode SrSc0.175Nb0.025Co0.8O3‑δ in The Oxygen Reduction Reaction. ACS Appl. Mater. Interfaces 2016, 8, 3003−3011. (25) Yi, J.; Schroeder, M.; Martin, M. CO2-tolerant and Cobalt-free SrFe0.8Nb0.2O3‑δ Perovskite Membrane for Oxygen Separation. Chem. Mater. 2013, 25, 815−817. (26) Porras-Vazquez, J. M.; Slater, P. R. Synthesis and Characterization of Oxyanion-Doped Cobalt Containing Perovskites. Fuel Cells 2012, 12, 1056−1063. (27) Jiang, L.; Wang, J.; Xiong, X.; Jin, X.; Pei, Q.; Huang, K. Thermal and Electrical Stability of Sr0.9Y0.1CoO2.5+δ as a Promising Cathode for Intermediate-Temperature Solid Oxide Fuel Cells. J. Electrochem. Soc. 2016, 163, F330−F335. (28) Istomin, S. Y.; Drozhzhin, O.; Svensson, G.; Antipov, E. Synthesis and Characterization of Sr1‑xLnxCoO3‑δ, Ln= Y, Sm-Tm, 0.1 ≤ x ≤ 0.5. Solid State Sci. 2004, 6, 539−546. (29) An, K. VDRIVE: Data Reduction and Interactive Visualization Software for Event Mode Neutron Diffraction. ORNL Report No. ORNL-TM-2012-621; Oak Ridge National Laboratory: Oak Ridge, TN, USA, 2012. (30) Chen, Y.; Rangasamy, E.; dela Cruz, C. R.; Liang, C.; An, K. A Study of Suppressed Formation of Low-Conductivity Phases in Doped Li7La3Zr2O12 Garnets by In Situ Neutron Diffraction. J. Mater. Chem. A 2015, 3, 22868−22876. (31) Larson, A. C.; Von Dreele, R. B. GSAS. General Structure Analysis System; LANSCE, Los Alamos National Laboratory: Los Alamos, NM, USA, 1994. (32) Toby, B. H. EXPGUI, A Graphical User Interface for GSAS. J. Appl. Crystallogr. 2001, 34, 210−213. (33) Le Bail, A. Whole Powder Pattern Decomposition Methods and Applications: a retrospection. Powder Diffr. 2005, 20, 316−326. (34) Momma, K.; Izumi, F. VESTA 3 for Three-dimensional Visualization of Crystal, Volumetric and Morphology Data. J. Appl. Crystallogr. 2011, 44, 1272−1276.

Department of Energy. We thank Mr. M. J. Frost from SNS for the technical support of the neutron experiment. We also thank Libin Lei for helping with TGA measurements. This manuscript has been authored by UT-Battelle, LLC under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a nonexclusive, paidup, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan).



REFERENCES

(1) Duan, C.; Hook, D.; Chen, Y.; Tong, J.; O’Hayre, R. Zr and Y Co-Doped Perovskite as a Stable, High Performance Cathode for Solid Oxide Fuel Cells Operating Below 500 °C. Energy Environ. Sci. 2017, 10, 176−182. (2) Chen, Y.; Zhou, W.; Ding, D.; Liu, M.; Ciucci, F.; Tade, M.; Shao, Z. Advances in Cathode Materials for Solid Oxide Fuel Cells: Complex Oxides Without Alkaline Earth Metal Elements. Adv. Energy Mater. 2015, 5, 1500537. (3) Burriel, M. n.; Peña-Martínez, J.; Chater, R. J.; Fearn, S.; Berenov, A. V.; Skinner, S. J.; Kilner, J. A. Anisotropic Oxygen Ion Diffusion in Layered PrBaCo2O5+δ. Chem. Mater. 2012, 24, 613−621. (4) Hildenbrand, N.; Nammensma, P.; Blank, D. H.; Bouwmeester, H. J.; Boukamp, B. A. Influence of Configuration and Microstructure on Performance of La2NiO4+δ Intermediate-Temperature Solid Oxide Fuel Cells Cathodes. J. Power Sources 2013, 238, 442−453. (5) Oh, D.; Gostovic, D.; Wachsman, E. D. Mechanism of La0.6Sr0.4Co0.2Fe0.8O3 Cathode Degradation. J. Mater. Res. 2012, 27, 1992−1999. (6) Kim, Y.-M.; Chen, X.; Jiang, S. P.; Bae, J. Effect of Strontium Content on Chromium Deposition and Poisoning in Ba1‑xSrxCo0.8Fe0.2O3‑δ (0.3 ≤ x ≤ 0.7) Cathodes of Solid Oxide Fuel Cells. J. Electrochem. Soc. 2012, 159, B185−B194. (7) Yang, W.; Hong, T.; Li, S.; Ma, Z.; Sun, C.; Xia, C.; Chen, L. Perovskite Sr1‑xCexCoO3‑δ (0.05≤ x ≤ 0.15) as Superior Cathodes for Intermediate Temperature Solid Oxide Fuel Cells. ACS Appl. Mater. Interfaces 2013, 5, 1143−1148. (8) Zhu, Y.; Sunarso, J.; Zhou, W.; Jiang, S.; Shao, Z. Highperformance SrNb0.1Co0.9‑xFexO3‑δ Perovskite Cathodes for LowTemperature Solid Oxide Fuel Cells. J. Mater. Chem. A 2014, 2, 15454−15462. (9) Shao, Z.; Haile, S. M. A High-Performance Cathode for The Next Generation of Solid-Oxide Fuel Cells. Nature 2004, 431, 170−173. (10) Magnone, E. A Systematic Literature Review On BSCF-based Cathodes for Solid Oxide Fuel Cell Applications. J. Fuel Cell Sci. Technol. 2010, 7, 064001. (11) Hu, S.; Wang, Y.; Cazorla, C.; Seidel, J. Strain-enhanced Oxygen Dynamics and Redox Reversibility in Topotactic SrCoO3‑δ (0