Subscriber access provided by University of Birmingham
Communication
Defective interfaces in Yttrium-doped Barium Zirconate films and consequences on proton conduction Nan Yang, Claudia Cantoni, Vittorio Foglietti, Antonello Tebano, Alex Belianinov, Evgheni Strelcov, Stephen Jesse, Daniele Di Castro, Elisabetta di Bartolomeo, Silvia Licoccia, Sergei V. Kalinin, Giuseppe Balestrino, and Carmela Aruta Nano Lett., Just Accepted Manuscript • DOI: 10.1021/acs.nanolett.5b00698 • Publication Date (Web): 19 Mar 2015 Downloaded from http://pubs.acs.org on March 23, 2015
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Nano Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
Defective interfaces in Yttrium-doped Barium Zirconate films and consequences on proton conduction Nan Yang1,2,3, Claudia Cantoni4, Vittorio Foglietti1,5, Antonello Tebano1,2,5, Alex Belianinov6, Evgheni Strelcov6, Stephen Jesse6, Daniele Di Castro1,5, Elisabetta Di Bartolomeo2,7, Silvia Licoccia2,7, Sergei V. Kalinin6, Giuseppe Balestrino1,5 and Carmela Aruta1,2,5*
1.
National Research Council CNR-SPIN, University of Roma “Tor Vergata”,
Rome I-00133, Italy. 2.
NAST Center, University of Roma “Tor Vergata”, Rome I-00133, Italy.
3.
Engineering Faculty, Università degli studi Niccolò Cusano, Rome I-00166,
Italy. 4.
Materials Science and Technology Division, Oak Ridge National Laboratory,
Oak Ridge, TN 37831, USA. 5.
Department DICII, University of Roma Tor Vergata, Rome I-00133, Italy.
6.
Center for Nanophase Materials Sciences, Oak Ridge National Laboratory,
Oak Ridge, TN 37831, USA. 7.
Department of Chemical Sciences and Technologies, University of Roma
“Tor Vergata”, Rome 00133, Italy.
1 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Abstract Yttrium doped Barium Zirconate (BZY) thin films recently showed surprising electric transport properties. Experimental investigations conducted mainly by electrochemical impedance spectroscopy suggested that a consistent part of this BZY conductivity is of protonic nature. These results have stimulated further investigations by local unconventional techniques. Here we use Electrochemical Strain Microscopy (ESM) to detect electrochemical activity in BZY films with nanoscale resolution. ESM in a novel cross sectional measuring set-up allows the direct visualization of the interfacial activity. The local electrochemical investigation is compared with the structural studies performed by state of art scanning transmission electron microscopy (STEM). The ESM and STEM results show a clear correlation between the conductivity and the interface structural defects. We propose a physical model based on a misfit dislocation network which introduces a novel 2D transport phenomenon, whose fingerprint is the low activation energy measured.
Keywords: ionic conduction, interface defects, doped barium zirconate, electrolytes, SPM, STEM, perovskite oxide thin films
Doped perovskite oxides have been widely investigated in recent years as proton conducting solid electrolytes for a variety of electrochemical devices,
2 ACS Paragon Plus Environment
Page 2 of 25
Page 3 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
such as fuel cells, hydrogen sensors, electrolyzers and hydrogen pumps1, 2. Excellent chemical stability of BaZr0.8Y0.2O3-δ (BZY) makes it one of the most promising electrolyte materials for protonic fuel cells.3,
4
However, the
conductivity of BZY, reported by many groups varies widely introducing major uncertainties for its implementation in real devices5. Theoretical and experimental studies on doped barium zirconate proved that this material exhibits high proton conductivity when the grain boundary contribution to the transport properties is suppressed. 6, 7. This effect has been later investigated in epitaxial BZY, 1 µm thick films, grown on (100)-oriented MgO substrates confirming the theoretical predictions8 and showing a protonic conductivity, measured by Electrochemical Impedance Spectroscopy (EIS), higher than oxygen ion conductors, in the temperature range of 300 to 600 °C. These results stimulated further studies of BZY properties in epitaxial thin films especially focused on the interplay between structural and transport properties. Surprisingly, very high values of conductivity, at temperatures of 550 °C - 600 °C, have been recently measured by EIS in few tens of nanometers thick BZY film grown on (110) NdGaO3 (NGO)9, with a 10% lattice mismatch between the film and substrate. These results suggest that heavily strained interfaces may represent a key parameter for tailoring defect densities in thin epitaxial films, thus affecting their conductivity.
10
Here we
report nanoscale evidence of higher electrochemical activity at the interface
3 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
by Electrochemical Strain Microscopy (ESM), in both in-plane and crosssection geometry. The ESM study, coupled with the scanning transmission electron microscopy (STEM) analysis, shows a clear correlation between the conductivity of our perovskite thin films and the defective interfaces. Crosssectional measurements with different scanning probe microscopy techniques already demonstrated to be a powerful tool to study the interface properties of complex oxide thin films11,
12, 13
. We extend such an approach to ESM to
directly visualize the interface reactivity and compare it with the crosssectional structural properties obtained by STEM. We demonstrate that the microscopic origin of such a high conductivity in BZY films is the strongly defective interface between film and substrate. Local electrochemical activity in electrolyte materials can be investigated at high spatial resolution by ESM14,15. This technique utilizes the SPM tip as a moving electrode as it scans over the sample surface. When bias is applied to the sample surface, at a sufficiently high probe bias, the electrochemical reaction at the tip-surface junction can be activated, resulting in generation or annihilation of mobile ionic species below the tip. The movement of those mobile species will result in a change in the molar volume (electrochemical strain). The associated dynamic surface deformation can be detected at the 2–5 pm level16,
17
. Protonic defects in BZY are formed from water gas phase
dissociation, where the hydroxide ion fills an oxide ion vacancy and the proton
4 ACS Paragon Plus Environment
Page 4 of 25
Page 5 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
forms a covalent bond with lattice oxygen6. By ESM, further protons can be injected or extracted by electrochemical reactions (see below) giving rise to the surface deformation. The initial state of hydration of our samples can be estimated from the thermodynamic data available for bulk BZY (see Supporting Information). Our films are expected to be fully hydrated even at room temperature and in the range of experimental partial pressures used. Considering that the films in this study are much thinner than 1µm, the water uptake can approach equilibrium in less than a few hours, for the entire range of experimental conditions. In ESM the conductive tip, with a nanometer sized apex, confines the electric field to a small volume of material (comparable with tip contact area) in presence of a water meniscus at the point of contact between the hydrophilic tip and the film. When the tip is positively biased with respect to the grounded electrode, with a potential sufficiently higher than the redox potential V ~ 1.23 V, the water oxidation reaction occurs, generating protons, oxygen gas and losing electrons to the anode (see Supporting Information). Protons will migrate inside the film due to the high local electric field generated by the tip, resulting in the following chemical reaction:
4× + 2 → 4∙ + + 4
(1)
Under negative tip bias the water evolution reaction occurs:
4∙ + + 4 → 4× + 2
5 ACS Paragon Plus Environment
(2)
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
These reactions have low redox potential and will be prevalent at low voltage tip bias. Further oxygen evolution and reduction reactions, involving creation and annihilation of oxygen vacancies inside the film, can appear at higher bias voltages (see Supporting Information and discussion below). Our measurements aimed at finding a correlation between strongly defective interface and electrochemical activity, and elucidating the microscopic origin of the high conductivity in BZY thin films9.
Figure 1. (a) Topography of a 1 × 1 µm2 area of 500 nm thick BZY film on NGO. The 600 × 600 nm2 region selected for the ESM measurement is highlighted by the black square. (b) Spatial distribution of the cantilever’s contact resonance frequency over the region indicated in (a). (c) Schematic of the cross-sectional ESM measurement. The second electrode not shown in the figure is grounding the film surface. (d) Hysteresis loops averaged over different 20 nm × 600 nm areas of the cross section map reported in (h): (A) is the film region close to the BZY film surface, (B) is close to interface with NGO substrate and (C) is within the NGO
6 ACS Paragon Plus Environment
Page 6 of 25
Page 7 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
substrate region. (e)-(h), ESM hysteresis loop area maps at different maximal bias voltage: 24 V (e), 45 V (f), 52 V (g), 60 V (h). In panel (h) the vertical dashed black line indicates the interface between BZY and NGO. The colored dashed squares indicate the regions where the ESM loops of panel (d) were acquired.
The BZY film on NGO substrate was fractured to perform a transverse measurement across the interface. The surface of the fractured BZY/NGO is not atomically flat, as it can be seen in the AFM topography of Fig. 1a. The height variation in the averaged topography line profile is found to be of the order of 10 nm. The main reason for the presence of a hill-valley morphology in the cross-sectional topography lies in the fracturing dynamics of two different materials18. ESM measurements on the cross section allowed direct visualization of the electrochemical interface activity, as shown in Fig.1. The schematics of the cross sectional ESM geometry is reported Fig. 1c. By comparing the AFM topography and the ESM contrast of Fig.1e-h, the interface can be precisely located and the left and right sides can be directly attributed to the BZY layer and NGO substrate. The ESM signal shows a hysteretic behavior, which arises from the positive and negative nucleation biases (PNB and NNB) required to activate the chemical processes of ionic injection or annihilation. The shape of the loops from different regions is shown in Fig.1d. The PNB and NNB are identified in the hysteresis loops as the voltage corresponding to the inflection points. In addition to the ESM signal, we have measured the resonance frequency shift
7 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
of the tip, as shown in Fig. 1b. The tip resonance frequency gives direct information about the changes in mechanical properties of the tip-surface junction, which can depend on friction, viscous flow, and plastic deformation. The enhanced surface interaction experienced by the tip is equivalent to an increase in elastic modulus of the material, thus increasing the resonance frequency19. To rule out the topography effect, we can compare the resonance frequency map with the ESM response map. The spatial distribution of the resonance frequencies on the sample surface (Fig. 1b) shows a large variation that reflects changes in the topography. On the contrary, the ESM signal shows a bias dependent response which depends on the electrode interaction with the underneath material. At voltages just above the nucleation bias, a slight but well-defined contrast signal can be measured, as shown in Figure 1(e). A large variation of the ESM signal at the interface is clearly visible increasing the bias voltages, as shown in Fig. 1 f-h. We would like to stress that the ESM response depends on the distance from the interface and that the intensity at the interface exhibits a strong bias dependence. The hysteresis loop data of the cross-section map reported in Fig. 1d have been averaged in three regions: close to the surface (A), close to the interface (B) and within the substrate (C). Both ESM loop area and loop height increase when the cross-sectional measurement is performed in the region close to the interface. By increasing the applied bias the contrast can be improved
8 ACS Paragon Plus Environment
Page 8 of 25
Page 9 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
suggesting that electrochemical reactions are occurring close to the interface. This is in agreement with the conductivity results obtained from previous electrochemical impedance spectroscopy measurements in BZY films on the NGO substrates9. Our cross-sectional ESM measurements confirm that BZY films can be viewed as consisting of two different regions, with the most electrochemically active located near the interface.
2
Figure 2. (a) Hysteresis loops averaged on a 1 × 1 µm surface area of two films with
different thickness, 300 nm and 20 nm. (b) Schematic of the measurements on the film surface where the probing depth of about 20 nm is highlighted. G indicates the grounding electrode.
Fig. 2 shows the ESM signal recorded in the standard planar geometry for two BZY films that are 300 nm and 20 nm thick, respectively. The probing depth of the ESM measurement is about 20 nm14, thus in the case of the thinner film
9 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 25
we are more sensitive to the interface region. The ESM signal shown in fig 2a is clearly stronger for the thinner film, thus confirming that the electrochemical activity
is
enhanced
close
to
the
interface,
consistently
with
the
electrochemical behavior obtained in the cross-sectional geometry (Fig 1d). The similarity between the ESM measurements in the cross-section and planar configuration allows us to use the results obtained on the two films to study the electrochemical process close and far away from the interface. The ESM loop curves shown in fig 2 are relative to a bias tip ranging from -35 up to 35 V. In this region the shape of the two curves looks very different. In the 20 nm thick film, the shape of the hysteresis curves does not show any sign of saturation. These observations support the hypothesis that the measurement falls inside of the ohmic transport region16,
20,
where the
transport of charged defects, i.e. protons, is linearly related to the local electrical conductivity of the film. In this regime larger bias results in increased loop area. In the case of the 300 nm thick film, Fig. 2a shows that the ESM response saturates at around ±15V. There can be different mechanisms responsible for the saturation of the hysteresis loop. One possibility is related to a mass transport limited reaction; we rule this out, however, due to the fact that the thinner film does not show any saturation in the loop curve, under the same conditions as the thicker film. We believe that in our set-up a lower mobility of charged protons in thicker films, as compared to the thinner film,
10 ACS Paragon Plus Environment
Page 11 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
may lead to the formation of a space charge region in the vicinity of the tip, preventing further injection of mobile species and thus explaining the saturation effects observed in the loop curve. Once the electric field is removed, there is a drift-diffusion of protons inside the BZY film, thus relaxing the strain signal. Analysis of the relaxation time reported in Supporting Information shows a different behavior for BZY films with different thickness. However, it does not allow for a full understanding of the proton dynamics, because it requires an accurate knowledge of the number of carriers injected and of the possible presence of compensating carriers for protons which would eventually enable ambipolar diffusion. To further elucidate the transport mechanism, we conducted humidity and temperature dependent ESM. Furthermore, to avoid the uncertainty related to the choice of the bias window in ESM and in order to separate the reaction and transport processes, we utilize the first order reversal curve (FORC) method21. FORC is a powerful tool to study the on-site transport process for a typical electrochemical reaction-diffusion system. It consists in progressively increasing the maximum bias voltage Vmax of the ESM measurements thus allowing systematic tracking of the evolution of loops. The envelope of probe bias is a bipolar excitation waveform with increasing amplitude. This allows multiple hysteresis loops with increasing bias to be collected at each spatial location.
11 ACS Paragon Plus Environment
Nano Letters
A demonstration of the proton injection by reactions (1) and (2) is given by the humidity dependence measurements reported in Figure 3. It can be observed that the loop area of the 20 nm thick BZY film increases with Vmax and this effect is enhanced by increasing the humidity RH from 10% to 80%. The inset of Fig. 3 shows the comparison between the thicker (300nm) and the thinner (20nm) films. Above the bias tip voltage of ~30-35 V, the ESM loop area of the
Loop area (a.u.)
300 nm thick film also increases.
3 3
Loop area (a.u.)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 12 of 25
2
3
300 nm 20 nm
2 1 40% RH 0 0
2
20
40
60
Bias (V)
1 10% RH 20% RH 40% RH 80% RH
1 0
0
10
20
30
40
50
60
Bias (V) Figure 3. Loop-area vs peak bias voltage for 20 nm BZY film with increasing Relative Humidity. Inset : comparison of the Loop-area vs peak bias voltage for the 20 nm and the 300 nm films, made at 40% RH (similar with ambient condition).
This is a consequence of oxygen evolution/reduction reaction:
× ⇆ ∙∙ + + 2
(3)
In fact, in case of 300 nm thick film, loop area does not depend on RH (see Supporting Information). Larger threshold voltage will be necessary to activate
12 ACS Paragon Plus Environment
Page 13 of 25
oxygen evolution/reduction reaction than the proton injection/water evolution reaction. This may explain the larger threshold voltage observed in the case of the 300nm thick film. (see Supporting Information).22 Temperature-dependent ESM allows us to describe the electrochemical activity in terms of the activation energy Ea, which is one of the key parameters of interest in describing ionic transport. In order to extract the activation energy of the ionic motion, we analyze the averaged loop area across a 1 µm2 area from room temperature to 110 °C. The temperature dependence of averaged loop area for BZY films under different bias voltages are given as ln(loop area) vs. 1/T in Figure 4.
Ea= 0.1 eV
-6
31 V 37 V 42 V
-7 300 nm 20 nm
-6
-8 -9
ln(loop area)
ln (loop area)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
-7 -8 -9
Vmax=31 V -3
-3
-3
2.7x10 3.0x10 3.3x10 -1
1/T(K ) -3
2.7x10
-3
3.0x10
-3
3.3x10
-1
1/T(K ) Figure 4. Temperature dependence of averaged loop area for the 20 nm BZY film under different bias voltages. Inset : comparison of the measurements for the 20 nm and the 300 nm films at 31 V peak voltage bias.
The behavior was fitted using an Arrhenius equation, shown by a solid line. 13 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Although the used narrow temperature range of 100 °C limits the accuracy of the measurement, the activation energy (Ea) of the 20 nm film was clearly around 0.1 eV. The inset of figure 4 shows a comparison of the measured activation energies between the 20 nm and the 300 nm thick films, at Vmax =31 V and identical environmental conditions. While the value of Ea measured for the 20 nm thick film is surprisingly low, the film of 300 nm shows an activation energy of 0.3 eV which is much closer to the values reported in literature for bulk protonic conductors, namely about 0.4 eV6. The activation energy measured by ESM is also lower than the previous values obtained by EIS in the case of thin films8,9. However, the range of temperatures of ESM measurements, 20-100 °C, is different than that used to perform conventional EIS measurements, 400-600 °C. A consistent experimental scenario is that the transport phenomena prevailing at lower temperature, whose fingerprint is in the low activation energy, may be obscured at higher temperatures by a further conductive path with higher activation energy. Very recently it has been theoretically demonstrated that the expected value of the activation energy in the BZY bulk should be as low as 0.1 eV, but the actual higher value can be a consequence of the proton trapping23. Our experimental findings, as discussed below, suggest that the presence of defects at the substrate-film interface, mainly dislocations, may create new channels for protonic conduction where the trapping mechanism is possibly
14 ACS Paragon Plus Environment
Page 14 of 25
Page 15 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
inhibited. We used state-of-art aberration-corrected STEM/EELS analysis to correlate the local information on the electrochemical activity obtained by ESM with the local structural information. The STEM/EELS results are shown in Fig. 5.
Figure 5. Structural and composition defects near the substrate interface. a) High resolution Z-contrast STEM images of a 10-nm-thick BZY with the NGO substrate, showing numerous dislocations at the interface and associated tilting of BZY film grains. b), c) and d) electron energy loss (EELS) maps for Ba, Zr and O, respectively. The maps were obtained by plotting the integrated intensity of Ba-N4,5, Zr-M4,5, and O-K in a spectrum image acquired in a region indicated by the yellow box in a).
The majority of defects observed in Fig. 5a close to the interface originate from the large mismatch with the NGO substrate. The cubic lattice parameter of BZY perovskite cell is a = 4.22 Å, to be compared with the pseudocubic lattice parameter of NGO which is 3.86 Å. There is a network of dislocations
15 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
at the interface spaced ~ 10 u.c. apart . The amount of strain relieved by the interface dislocations can be estimated using the equation D = bZ/(fO-εO), where D is the average dislocation spacing, fO defines the original misfit between lattices (fO = 0.088), bZ is the Burgers vector component along the interface plane (bZ = 1/2 a), and εO is the residual misfit strain. The obtained value for εO indicates that a large part of the original misfit strain is relieved in the ~ 5 nm-thick interfacial layer containing the dislocations. The EELS spectra of Fig.5b-d show a Ba depletion and Zr accumulation close to the interface. It has been theoretically demonstrated that defects concentrations affect the proton incorporation in proton conductors.24 In particular, proton incorporation can be enhanced with increasing the cation vacancy concentration. In addition, the increased cation distances may reduce the trapping effects leading to a lower proton activation energy.6 However, we cannot discuss the exact interface stoichiometry variation because EELS spectra measurements at the Y edges were not feasible. In the following we concentrate our discussion on the structural defects and their effect on the proton conduction. To investigate the evolution of the strain above the 5 nm interfacial layer, we generated maps of the in-plane and out-of-plane BZY lattice spacing by calculating positions (xij, yij) of the center of mass of the Ba columns in the Zcontrast image of the region indicated in Fig. 6b by a dashed square. The c-
16 ACS Paragon Plus Environment
Page 16 of 25
Page 17 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
parameter, given by yij-yi(j-1) is mapped in Fig. 6a and its average along the x direction is plotted in the graph of Fig. 6c. The figure shows that the c-lattice parameter close to the interface is smaller than expected on the basis of an in-plane compressive strain. After several tens of unit cells the bulk c-lattice value is recovered.
Figure 6. (a), c-spacing map of the BZY thicker film over the region indicated by the dashed box in the Z-STEM images of panel (b). (c), BZY c-spacing average along the x direction in a, as a function of y (distance from interface). Dashed line indicates the BZY bulk value.
ESM results suggest that the thinner BZY film possesses a larger electrochemical activity and lower transport activation energy. The structural characterization reveals a large density of the dislocations near the interface. The defects are clearly identified as misfit dislocations spaced ~ 10 unit cells apart. These dislocations reduce the film misfit within a highly defective 17 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
interface layer of less than 5 nm in thickness. The remaining film strain, generated by lattice tilts and antiphase boundaries is relieved within a thickness of 50-60 unit cells or ~ 20 nm from the interface. There are many examples in literature, showing that dislocations can give a major contribution to modification of the transport properties of ionic and electronic conductors.25,
26
Dislocations have been reported to have an effect
of increased oxygen mobility in oxides,27 while a dislocation mediated mechanism reported recently has been shown to enhance the surface exchange coefficient, decreasing the relaxation time28. However, further experimental and theoretical studies do not confirm any evidence of fast diffusion of oxygen ions along the dislocation array29. In the case of proton conductors the situation is even more controversial because the theoretical studies reported in literature are mainly devoted to the space charge effects on the segregation energy of oxygen vacancies and protons in grain boundaries, and do not consider local spatial correlations between the defects within dislocation cores.30 Nevertheless, it is well known that proton transfer between octahedral sites determines the activation energy of the proton diffusion. One of the main factors affecting the activation energy is the distance between the cations and oxygen ions. However, density functional theory calculations of atomic-scale proton movements31 demonstrate that a structural variation of 5% in the octahedral oxygen distance of BZY should not
18 ACS Paragon Plus Environment
Page 18 of 25
Page 19 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
substantially vary the energy barriers of proton migration, which remain in the range 0.4-0.5 eV. Therefore, the c-axis compression alone cannot explain the very low activation energy obtained in the case of the thinnest BZY film. Very recently, a strong dependence of the dislocations structure at the SrTiO3/MgO heterointerface on the ionic transport properties has been demonstrated32. The correlation between our ESM and STEM results suggests that similar mechanisms can be active is this case. A quantitative estimate of the contribution that such an array of dislocations can provide to proton conductivity must rely on ad hoc atomistic simulations. Nevertheless, at the film-substrate interface, the atomic arrangement is highly perturbed along every dislocation line, where chemical defects are also observed. Both cation vacancies and structural distortion may play a relevant role: charge unbalance due to cation vacancy can be compensated by proton incorporation,24 while the local distortions of oxygen octahedra, which coordinate Zr or Y, may reduce the trapping effect.23 Therefore, the presence of a two-dimensional (2D) dislocation network at the interface with the substrate, as shown in Fig. 7, provides easily accessible sites for anchoring H+ and allows for the low activation energy to be measured in our thinner BZY films.
19 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Figure 7. Schematic drawing of the interface between BZY film and NGO substrate with the dislocation network as deduced from the STEM structural characterization.
In conclusion, many experimental results reported in literature show unusual high values of ion conductivity in very thin films which are explained in terms of defective interfaces.33 In this framework, our study demonstrates a clear correlation between the local electrochemical investigation by ESM with the structural studies performed by STEM. Our results suggest a 2D transport mechanism occurring across the misfit dislocation network with a very low activation energy which gives to a contribution to the high values of proton conductivity obtained in our thin BZY films. However, we stress that our experimental results and the qualitative explanation proposed are related to the in-plane conductivity. Conventional solid oxide fuel cells (SOFCs) consist of a cathode and anode on the opposite surfaces of an electrolyte membrane. There are two separated and sealed gas chambers which supply the fuel and
20 ACS Paragon Plus Environment
Page 20 of 25
Page 21 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
the oxidant gases. In the BZY-NGO system we discuss here, the nature of the high conductivity values requires the use of a planar technology without a membrane. Indeed, when the membrane is present, the transverse conductivity rules the transport properties. However, the use of free standing membranes is technologically very challenging considering the thickness of our BZY films. In perspective, to exploit the extraordinary properties of these strongly defective interfaces, dual chamber SOFCs using a planar technology should be fabricated using the state of the art microfluidic and microelectromechanical system technology.
Author information Corresponding Author *E-mail:
[email protected] Notes The authors declare no competing financial interest.
Acknowledgements The authors acknowledge META—Materials Enhancement for Technological Applications
Project
(FP7-PEOPLE-2010-IRSES—Marie
Curie
Actions,
PIRSES-GA-2010-269182. Italian MIUR is acknowledged for support through the FIRB Project RBAP115AYN “Oxides at the nanoscale: multifunctionality and applications” and PRIN Project 2010-2011 OXIDE, “OXide Interfaces: emerging new properties, multifunctionality, and Devices for Electronics and Energy. The research at ORNL was conducted at the Center for Nanophase Materials Sciences, which is sponsored at Oak Ridge National Laboratory by the Division of Scientific User Facilities, U.S. Department of Energy, in the 21 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
project CNMS2013-032 “Local electrochemical characterization of epitaxial thin films of proton conductor perovskite oxides”. CC acknowledges support by the U. S. Department of Energy, Basic Energy Sciences, Materials Science and Engineering Division.
References 1
Iwahara, H.; Yajima, T.; Hibino, T.; Ozaki, K.; Suzuki, H.; Protonic conduction
in calcium, strontium and barium zirconates, Solid State Ionics 1993, 61, 6569. 2
Ito, N.; Iijima, M.; Kimura, K.; Iguchi, S.; New intermediate temperature fuel
cell with ultra-thin proton conductor electrolyte, J. Power Sources 2005, 152, 200. 3
Ryu, K. H.; Haile, S. M.; Chemical stability and proton conductivity of doped
BaCeO3–BaZrO3 solid solutions, Solid State Ionics 1999, 125, 355. 4
Malagoli, M.; Liu, M. L.; Park, H. C.; Bongiorno, A. Protons crossing triple
phase boundaries based on a metal catalyst, Pd or Ni, and barium zirconate, Phys. Chem. Chem. Phys., 2013 15, 12525. 5
Babilo, P.; Uda, T.; Haile, S. M.; Processing of yttrium-doped barium
zirconate for high proton conductivity, Journal of Materials Research, 2007, 22, 1322-1330. 6
Kreuer, K. D.; Proton-conducting oxides, Annual Review of Material
Research, 2003, 33, 333-359. 7
Kreuer, K. D.; Aspects of the formation and mobility of protonic charge
carriers and the stability of perovskite-type oxides, Solid State Ionics 1999, 125, 285–302. 8
Pergolesi, D.; Fabbri, E.; D’Epifanio, A.; Di Bartolomeo, E.; Tebano, A.;
Sanna, S.; Licoccia, S.; Balestrino, G.; Traversa, E.; High proton conduction in grain-boundary-free yttrium-doped barium zirconate films grown by pulsed
22 ACS Paragon Plus Environment
Page 22 of 25
Page 23 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
laser deposition, Nat. Mater., 2010, 9, 846–852. 9
Foglietti, V., Yang, N., Tebano, A., Aruta, C., Di Bartolomeo, E., Licoccia, S.,
Cantoni, C.,& Balestrino, G., Heavily strained BaZr0.8Y0.2O3-x interfaces with enhanced transport properties, Applied Physics Letters, 104, 081612 (2014). 10
Garcia-Barriocanal, J., Rivera-Calzada, A., Varela, M., Sefrioui, Z., Iborra,
E., Leon, C., Pennycook, S. J., & Santamaria, J., Colossal Ionic Conductivity at Interfaces of Epitaxial ZrO2:Y2O3/SrTiO3 Heterostructures, Science, 321, 676-680 (2008). 11
Basletic, M., Maurice, J.-L., Carrétéro, C., Herranz, G., Copie, O., Bibes,
M., Jacquet, É., Bouzehouane, K., Fusil, S. & Barthélémy, A., Mapping the spatial distribution of charge carriers in LaAlO3/SrTiO3 heterostructures, Nature Materials 7, 621 - 625 (2008) 12
Chien, T. Y., Kourkoutis, L. F., Chakhalian, J., Gray, B., Kareev, M.,
Guisinger, N. P., Muller, D. A., & Freeland, J. W., Visualizing short-range charge transfer at the interfaces between ferromagnetic and superconducting oxides, Nature Comm., 4, 2336 (2013) 13
Chien, T. Y., Chakhalian, J., Freeland, J. W., & Guisinger, N. P., Cross-
Sectional Scanning Tunneling Microscopy Applied to Complex Oxide Interfaces, Adv. Funct. Mater. 23, 2565–2575 (2013) 14
Morozovska, A. N., Eliseev, E. A., Balke, N., & Kalinin, S. V., Local probing
of ionic diffusion by electrochemical strain microscopy: spatial resolution and signal formation mechanisms, J. Appl. Phys. 108, 053712 (2010). 15
Doria, S., Yang, N., Kumar, A., Jesse, S., Tebano, A., Aruta, C., Di
Bartolomeo, E., Arruda, T. M., Kalinin, S. V., Licoccia, S., Balestrino, G., Nanoscale mapping of oxygen vacancy kinetics in nanocrystalline Samarium doped ceria thin films, Appl. Phys. Lett. 103, 171605 (2013). 16
Kumar, A., Ciucci, F., Morozovska, A. N., Kalinin, S. V., & Jesse, S.,
Measuring oxygen reduction/evolution reactions on the nanoscale, Nature 23 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Chemistry, 3, 707-713 (2011). See also Supplementary Informations. 17
Kumar, A., Ciucci, F., Leonard, D., Jesse, S., Biegalski, M., Christen, H.,
Mutoro, E., Crumlin, E., Shao-Horn, Y., Borisevich, A. & Kalinin, S. V., Probing Bias-Dependent Electrochemical Gas–Solid Reactions in (LaxSr1–x)CoO3–δ Cathode Materials, Adv. Func. Mat. 23, 5027–5036 (2013). 18
Chien, T. Y., Liu, J., Chakhalian, J., Guisinger, N. P., & Freeland, J. W.,
Visualizing nanoscale electronic band alignment at the LaCaMnO3:SrTiO3 interface, Phys. Rev. B, 82, 041101(R) (2010). 19
Lucas, M., Mai, W., Yang, R., Wang, Z. L., & Riedo, E., Aspect ratio
dependence of the elastic properties of ZnO nanobelts, Nano Letters, 7, 13141317 (2007). 20
Yang, N., Belianinov, A., Strelcov, E., Tebano, A., Foglietti, V., Di Castro, D.,
Schlueter, C., Lee, T.-L., Baddorf, A. P., Balke, N., Jesse, S., Kalinin, S. V., Balestrino, G., & Aruta, C., Effect of Doping on Surface Reactivity and Conduction Mechanism in Samarium-Doped Ceria Thin Films, ACS Nano, 8, 12494–12501 (2014) 21
Balke, N., Jesse, S., Kim, Y., Adamczyk, L., Ivanov, I. N., Dudney, N. J. &
Kalinin, S. V., Decoupling electrochemical reaction and diffusion processes in ionically-conductive solids on the nanometer scale, ACS Nano, 4, 7349 (2010). 22
Bard, A. J.; Faulkner, L. R. Electrochemical methods: fundamentals and
applications. NewYork: Wiley, 1980. 23
Yamazaky, Y.; Blanc, F.; Okuyama, Y.; Buannic, L.; Lucio-Vega, J. C.;
Proton Trapping in yttrium-doped barium zirconate, Nat. Mater. 2013, 12, 647651. 24
Song, S.-J., Wachsman, E. D., Dorris, S. E., Balachandran, U., Defect
chemistry modeling of high-temperature proton-conducting cerates, Solid State Ionics 149, 1 – 10 (2002) 24 ACS Paragon Plus Environment
Page 24 of 25
Page 25 of 25
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
25
Otsuka, K.; Kuwabara, A.; Nakamura, A.; Yamamoto, T.; Matsunaga, K.;
Ikuhara, Y.; Dislocation-enhanced ionic conductivity of yttria-stabilized zirconia, Appl. Phys. Lett. 2003, 82, 877. 26
Chang, C. P.; Chu, M.-W.; Jeng, H. T.; Cheng, S.-L.; Lin, J. G.; Yang, J.-R.;
Chen, C. H. Condensation of two-dimensional oxide-interfacial charges into one-dimensional electron chains by the misfit-dislocation strain field, Nat. Commun. 2014, 5, 3522. 27
Szot, K., Speier, W., Bihlmayer, G. & Waser, R., Switching the electrical
resistance of individual dislocations in single-crystalline SrTiO3, Nature Materials, 5, 312 - 320 (2006). 28
Yan, L., & Salvador, P. A., Substrate and Thickness Effects on the Oxygen
Surface Exchange of La0.7Sr0.3MnO3 Thin Films, Appl. Mat. Int. 4, 2541−2550 (2012). 29
Metlenko, V., Ramadan, A. H. H., Gunkel, F., Du, H., Schraknepper, H.,
Hoffmann-Eifert, S., Dittmann, R., Waser, R. & De Souza, R. A., Do dislocations act as atomic autobahns for oxygen in the perovskite oxide SrTiO3?, Nanoscale 2014, 6, 12864. 30
Helgee, E.E.; Lindman, A; Wahnstrom, G. Origin of Space Charge in Grain
Boundaries of Proton-Conducting BaZrO3, Fuel Cells 2013, 13, 19-28. 31
Merinov, B.; Goddard III, W. Proton diffusion pathways and rates in Y-doped
BaZrO3 solid oxide electrolyte from quantum mechanics, J. Chem. Phys.
2009, 130, 194707. 32
Dholabhai, P. P.; Pilania, G.; Aguiar, J. A.; Misra, A.; Uberuaga, B. P.
Termination
chemistry-driven
dislocation
structure
at
SrTiO3/MgO
heterointerfaces, Nat. Commun. 2014, 5, 5043. 33
Leon, C.; Santamaria, J.; Boukamp, B. A. Oxide interfaces with enhanced
ion conductivity, MRS BULLETIN 2013, 38, 1056-1063.
25 ACS Paragon Plus Environment