Strain-Induced Tailoring of Oxygen-Ion Transport in Highly Doped

Dec 4, 2017 - Strain-Induced Tailoring of Oxygen-Ion Transport in Highly Doped CeO2 Electrolyte: Effects of Biaxial Extrinsic and Local Lattice Strain...
0 downloads 0 Views 2MB Size
Letter www.acsami.org

Cite This: ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Strain-Induced Tailoring of Oxygen-Ion Transport in Highly Doped CeO2 Electrolyte: Effects of Biaxial Extrinsic and Local Lattice Strain Junsung Ahn,†,‡ Sungjun Choi,† Kyung Joong Yoon,† Ji-Won Son,† Byung-Kook Kim,† Jong-Ho Lee,† Ho Won Jang,‡ and Hyoungchul Kim*,† †

High-Temperature Energy Materials Research Center, Korea Institute of Science and Technology, 5 Hwarang-ro 14-gil, Seongbuk-gu, Seoul 02792, Republic of Korea ‡ Department of Materials Science and Engineering, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Republic of Korea S Supporting Information *

ABSTRACT: We explored oxygen-ion transport in highly doped CeO2 through density-functional theory calculations. By applying biaxial strain to 18.75 mol % CeO2:Gd, we predicted the average migration-barrier energy with six different pathways, with results in good agreement with those of experiments. Additionally, we found that the migration-barrier energy could be lowered by increasing the tetrahedron volume, including the space occupied by the oxygen vacancy. Our results indicate that the tetrahedron volume can be expanded by larger codopants, as well as biaxial tensile strain. Thus, the combination of thin-film structure and codoping could offer a new approach to accelerate oxygen-ion transport. KEYWORDS: oxygen-ion transport, strain effects, solid-state ionics, doped CeO2, SOFC

I

bilayers was increased by around 1 order of magnitude, along with enhanced chemical stability of the δ-Bi2O3. Recently, Shi et al.5 demonstrated reliable manipulation of oxygen-ion conduction with strain energy in a free-standing nanoscale membrane of Ce0.8Gd0.2O1.9‑δ, utilizing an optimized Pt microelectrode design. Despite comprehensive studies on the development of fast oxygen-ion conductors using thin-filmbased nanostructures, it has remained unclear how much fast ionic conduction can be achieved by lattice strain, and how to control lattice strain in a material system. These two subjects are linked to the question of reproducibility of many reported experimental studies, and the inconsistency between computational predictions and experimental results.3 The conductivity enhancement reported by Garcia-Barriocanal et al.,2 in particular, has never been achieved in other studies, even those using similar material systems.9 Therefore, it is necessary to develop detailed computational predictions of oxygen-ion migration-barrier energy (Em) that take into consideration the effects of biaxial extrinsic strain under a wide range of strain conditions, and local lattice strain induced by dopants that have different ionic radii from host ions. Such predictions should substantially contribute to the development of superior and fast oxygen-ion conductors.

on-conducting nanostructures have attracted much attention for their applications in a variety of energy-conversion and energy-storage systems such as batteries,1 fuel cells,2−6 and electrolysis cells,7 owing to their improved ion-transport properties. In particular, with the development of vacuumbased deposition processes such as molecular beam epitaxy,8 atomic layer deposition,1 and pulsed laser deposition,3−7 the fabrication of thin-film ionic nanostructures has become an important technology that has led to great innovations. Over the past decade, attempts have been made to enhance ionic conductivity through combinations of such thin-film nanostructures (e.g., ultrathin,2 multicoated,4 and free-standing layers5); these efforts have mainly utilized the effects of strain (ε) resulting from external force or internal stress of thin-film structured materials. A report by Garcia-Barriocanal et al.2 on the colossal oxygen-ion conductivity at the interfaces of ZrO2:Y2O3/SrTiO3 heterostructures, which showed an enhancement of more than 8 orders of magnitude, was the initiator of research in this field. Since then, many researchers have observed various strain-induced enhancements of ion conductivity in nanostructures such as ultrathin, multicoated, and free-standing layers. Fluri et al.6 reported the suppression of oxygen-ion migration energy (−0.05 eV under ε = 0.35%) in an ultrathin Sm-doped CeO2 (SDC) film with thickness of 20 nm, whose tensile strain was induced by a SrTiO3 buffer layer. Sanna et al.4 also reported that the ion conductivity of a multilayer stack consisting of 20 CeO2:Gd/δ-Bi2O3:Er2O3 © XXXX American Chemical Society

Received: September 5, 2017 Accepted: December 1, 2017

A

DOI: 10.1021/acsami.7b13440 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Letter

ACS Applied Materials & Interfaces

listed in Table 1, the predicted elastic properties of a SQS GDC structure were in general agreement with the experimental

In the present study, we used density-functional theory (DFT) to investigate the effects of biaxial extrinsic and local intrinsic lattice strain on oxygen-ion transport in a highly doped CeO2 structure. Taking the thin-film layout into consideration, we examined the effects of strain on oxygen-ion transport under biaxial strain conditions. The effects of local lattice strain induced by the highly doped structure were also analyzed using various atomic configurations and dopant ions with different radii. On the basis of our calculation results, we will propose a strategy to accelerate oxygen-ion transport using the effects of both biaxial extrinsic and local lattice strain. To describe the complex configuration of a highly doped CeO2 structure, we used a special quasi-random structure (SQS) model of 18.75 mol % Gd-doped CeO2 (GDC), which is close to the GDC composition with 20 mol % of Gd that is commonly used in practical applications. This approach generates the atomic structure of a random alloy in a computationally efficient way by reproducing the nearestneighbor (NN) pair and multisite correlation functions of an infinite random alloy.12,13 Even though the SQS model is not the lowest-energy configuration, it adequately preserves and shows the local environment of random alloys. Applying the SQS approach to highly doped CeO2, a SQS unit cell for GDC with 18.75 mol % Gd, as shown in Figure 1, was obtained; the

Table 1. Calculated Thermo-Mechanical Properties, Elastic Constant (Cij), Shear Modulus (G), Young’s Modulus (E), and Poisson’s Ratio (ν) of our SQS Model of CeO2 Doped with 18.75 mol % Gd (18.75 mol % GDC)a Property

This work

Reference

C11 (GPa) C12 (GPa) C44 (GPa) G (GPa) E (GPa) ν

285.6 84.3 48.1 65.0 171.3 0.31

− − − 71.710 187.011 0.3310

a

The materials reported in refs 10 and 11 were 10 mol % GDC and 20 mol % GDC, respectively.

results, indicating that the SQS cell and its DFT treatment effectively reproduced the structural features of GDC. (See the Supporting Information for computational details.) Our quantitative analysis of the oxygen-ion transport phenomenon in highly doped CeO2, i.e., SQS 18.75 mol % GDC, was carried out using the nudged−elastic band (NEB) method for different oxygen sites and migration pathways. Considering the number of Gd ions at the first NN cation sites, four different oxygen sites and six different oxygen−oxygen migration pathways in the out-of-plane ⟨100⟩ direction were observed. Note that only the oxygen ions hopping along the ⟨100⟩ direction were considered in this work, because the migration path in the ⟨100⟩ direction is usually reported as the lowest-energy path, rather than the paths in the ⟨110⟩ and ⟨111⟩ directions.16 Figure 2 shows various atomic configurations and migration pathways, along with the four designated oxygen positions.

Figure 1. Schematic illustration of strain-induced tailoring of oxygenion transport in highly doped CeO2. The image in the inset shows a SQS unit cell for CeO2 doped with 18.75 mol % Gd (18.75 mol % GDC) as used in this study.

Gd ions substitute for cations in six octahedral cation sites among the 32 cation sites, and thus three oxygen vacancies are placed under the electroneutrality condition as follows: × Gd 2O3 → 2Gd′Ce + 3OO + V″O

(1)

Hence, it seems reasonable to investigate the effect of strain on the SQSs for 18.75 mol % GDC to account for the effect of a complex configuration of random alloys that has yet to be covered in DFT calculations. To investigate the effect of strain on oxygen-ion migration, we performed DFT calculations, applying anisotropic strain ranging from −3.78% to +2.14% along the lateral (x-y) and longitudinal (z) directions, with Poisson’s ratio (ν) defined as ν = εzz/(εzz − 2εxx) = εzz/(εzz − 2εyy), where εxx, εyy, and εzz are the strains in the x, y, and z directions, respectively. Because the epitaxial strain is experimentally limited within 1% by introducing dislocation above the critical thickness of a thin film,14,15 we believe it is reasonable to carry out extensive analysis in such a wide range of biaxial extrinsic strain. In addition, based on the predicted thermo-mechanical properties

Figure 2. Variety of oxygen-ion migration paths in the out-of-plane direction of SQS 18.75 mol % GDC. Pathi‑j indicates the migration path between Oi and the adjacent Oj site.

To derive the averaged migration-barrier energy (Em,a) from different values of Em,ij, which is the migration-barrier energy of Pathi‑j, we assumed that the SQS GDC cell had N parallel migration paths with the same migration length (L). In a consecutive migration step, we also assumed that the migration step with the highest Em,ij determined the migration-barrier energy of the path (see Figure S1). Then, the total conductivity (σ) of the parallel circuit can be expressed as B

DOI: 10.1021/acsami.7b13440 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Letter

ACS Applied Materials & Interfaces

Figure 3. (a) Em as a function of the biaxial extrinsic strain. The solid lines represent computational and experimental results obtained from the literature.6,18,19 The shaded area corresponds to the range of Em for bulk 20 mol % GDC reported previously.20−22 The dashed and dash-dotted lines are the guidelines of the symmetric and asymmetric paths, respectively. (b) Em as a function of tetrahedron volume of the initial state, including the space occupied by the oxygen vacancy. Same-color symbols indicate migration paths with the same number of dopant atoms at their first NN sites.

σ=

1 R total

L = A



1 L = R ij A

∑ σij

Aij A

18.75 mol % GDC indeed describes the nature of a highly doped material with a complex configuration. From these results of migration-energy calculations in a SQS, we found two novel tendencies for oxygen-ion transport in a strained GDC structure. First, the change in slope of the plot of Em versus strain according to the structural symmetry was confirmed. As shown in Figure 2, the structural symmetry indicates whether the dopant configuration of the tetrahedron (including an oxygen vacancy) that makes up the start and end points of the path is the same. In our SQS GDC structure, one group has a symmetric dopant configuration on the initial−final state (Path0‑0, PathI‑I, PathII‑II), whereas the other has an asymmetric dopant configuration (Path0‑I, PathI‑II, PathII‑III). The slope of Em versus strain was higher in the asymmetric state, indicating that the Em of the asymmetric path was more affected by the biaxial extrinsic strain. Second, Em and the tetrahedron volume (Vi) consisting of an oxygen vacancy and the first NN cations are inversely proportional to each other (see Figure 3b). This tendency can be easily understood, because oxygen-ion transport can be allowed by distorting the three-dimensional cation sublattice.16,18,19,23 Thus, the larger tetrahedron volume leads to a smaller lattice distortion. As is well-known, the critical-radius (rcrit) models reported in literature16,19,23 describe the migration process of oxygen ions to be in close correlation with the oxygen-ion radius and the two-dimensional aperture size formed by the edge cations. To consider the threedimensional factors (e.g., atomic configuration including nonedge NN cations and oxygen-ion transport pathway) governing oxygen-ion migration; however, we suggest a new descriptor: a “tetrahedron-volume” model. This model is an extension of the existing two-dimensional critical-radius model (see the linear correlation of Vi and rcrit in Figure S2), but enables more accurate prediction of the oxygen-ion migrationbarrier energy including complex biaxial extrinsic and local multidopant strain conditions. Consequently, it seems reasonable to suppose that variations in tetrahedron volume would be a key parameter in the modification of the barrier energy of oxygen-ion migration. In general, the tetrahedron volume seems to be modified by biaxial extrinsic tensile or compressive strain, which is generally induced by lattice mismatch at an interface between heterophases. Here, we would like to emphasize that the tetrahedron volume also varies with the effects of local lattice

(2)

where R and A indicate resistance and area, respectively, and the subscript ij denotes the properties of the migration path with Em,ij. Because we only consider the ⟨001⟩ path in the cubic oxygen-ion sublattice, the area of each migration path can be assumed to be same for all paths; thus, the area fraction (Aij/A) can be determined by the fraction of a migration path which has Em,ij in the total migration paths, nij/N (N = 16 in this work). From the Arrhenius relation, the total conductivity can be expressed as σT =

∑ fij ·σijT = ∑ fij γijexp(−Em,ij /kT )

(3)

where k is the Boltzmann constant, T is the temperature, and γij is a pre-exponential factor related to the average jump distance, the total charge-carrier concentration, and the jump frequency prefactor.17 Because the jump frequency prefactor only depends on lattice vibration frequency and the migration entropy term, γij can be assumed to be equal for all migration paths. Finally, expressing the total conductivity in Arrhenius form with one single Em, Em,a can be derived as Em,a = −kT ln ∑ fi exp( −Em, i /kT )

(4)

The value of Em predicted by the DFT-NEB method is plotted as a function of the biaxial extrinsic strain, εxx (= εyy) in Figure 3a. We observed general monotonic behavior for all migration paths: the biaxial tensile strain decreased Em, whereas the biaxial compressive strain increased Em. This result is highly consistent with those of previous computational and experimental studies.6,18,19 In particular, Fluri et al.6 reported that a tensile strain of 0.35% in the thin-film structure of 15 mol % SDC decreased Em by ∼0.05 eV (∼5.9%), which is consistent with the result of this work (in which Em decreased by ∼0.03 eV (∼5.3%) at εxx = 0.35%). Using eq 4, we also calculated the value of Em,a, which is shown in Figure 3a. Interestingly, Em,a seemed to converge on the two most frequent migration-path energies (Em,0−0 and Em,I−I), i.e., Em,a = 0.65 eV at εxx = 0%; this was predictable, as oxygen ions might detour to avoid crossing a high Em. Furthermore, for the unstrained state, the Em,a value of 0.65 eV lay within the range reported in the literature for 20 mol % bulk GDC (0.55−0.76 eV).20−22 Therefore, the SQS of C

DOI: 10.1021/acsami.7b13440 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Letter

ACS Applied Materials & Interfaces

Figure 4. (a) Schematic illustration of local lattice strain of GDC (in PathII‑III) with large La3+ or small Sc3+ codopant under biaxial extrinsic tensile stress. (b) Em as a function of ionic radius of dopant in doped CeO2; the solid symbols represent Em of singly doped CeO2, and the half-filled symbols represent Em of codoped CeO2. For [28] (olive symbol), Em of the codoping system is plotted as a function of average ionic radius, because the two dopant atoms are substituted by the same amount. (c) Predicted change in Em in PathII‑III of GDC after application of biaxial extrinsic tensile strain (1%; orange arrow) and codopant substitution (with La3+ or Sc3+ ion; blue arrow).

agreement with those of the experiments. More specifically, Kahlaoui et al.28 reported that Ce0.8Sm0.1La0.1O1.9 has higher oxygen-ion conductivity and lower activation energy than singly doped CeO2, which is similar to the findings of the current work. Finally, the extent to which the Em can be reduced by applying both biaxial extrinsic and local lattice tensile strain is shown in Figure 4c. Given that the limitation of biaxial extrinsic tensile strain is 1%, which reduces Em by 9%, a further 26% reduction can be attained by codoping of Gd and La ions. In summary, we investigated the effect of strain on various migration paths in highly doped CeO2 produced by the SQS approach. It is worth noting that the calculated thermomechanical properties and oxygen-ion Em of highly doped CeO2 are comparable with values reported in the literature, which indicate that the SQS structure adequately describes the nature of highly doped materials. With a computationally efficient description of the complex atomic configuration, we found that the Em of oxygen-ion can be controlled by modifying not only the biaxial extrinsic strain but also the local lattice strain, for example, by changing the dopant configuration and number of dopants. Consequently, we propose that the limitation of the biaxial extrinsic strain below ε < 1% can be overcome by applying local tensile strain with large codopants, which would expand the tetrahedron space of the oxygen ions without changing the mobile-carrier concentration. Consequently, this work sheds light on a new approach to control strain in nanostructured oxygen-ion-conducting materials by considering both the external strain at the heterophase interface and the local strain induced by codoping.

strain, e.g., the number of dopants (ndopant). This indicates that the barrier energy of oxygen-ion migration can be further modified by the local lattice strain induced by ndopant. This effect of dopants on the Em has mostly been attributed to two main factors: first, a lattice distortion caused by the difference between the ionic radii (r) of host and dopant cations; second, attractive columbic interaction between a negatively charged dopant (e.g., Gd′Ce) and a positively charged oxygen vacancy (VO·· ).16 Notably, Nakayama et al.16 reported that the Coulombic interaction between charged point defects became more pronounced as the coordination number of dopants around the oxygen vacancy increased, which is consistent with our work (e.g., PathII‑II). Thus, it has been suggested that lattice distortion, which is demonstrated by coupling of different ionic radii and the Coulombic interaction, affects the Em.24,25 Consequently, the tetrahedron shrinkage and increase of Em with increasing number of Gd ions at first NN sites (Figure 3b) could be attributed to the more prominent attractive columbic interactions between Gd′Ce and V··O. On the basis of the above conclusions, we propose a new strategy to tailor the barrier energy of oxygen-ion migration by employing biaxial extrinsic strain as well as local lattice strain. As shown in Figure 4a, combining cosubstitution with a dopant that is larger than a Gd ion and biaxial extrinsic strain in an epitaxial thin-film structure offers a good solution to expand the tetrahedron volume and lower the Em without any deterioration in oxygen-ion concentration. For example, we considered the unstrained PathII‑III, which led to the highest Em for all strain states, and the remaining Ce site (rCe = 0.94 Å) at the initial state of PathII‑III was replaced with a large La ion (rLa = 1.16 Å) or a small Sc ion (rSc = 0.87 Å).26 The predicted Em under the unstrained state (εxx = 0%) is plotted in Figure 4b as a function of the ionic radius of the dopant to show the effect of codoping on the Em in the SQS GDC structure; experimental results of codopants reported in the literature are also included for comparison.27−30 Interestingly, we found that cosubstitution of larger dopants effectively lowers the Em of oxygen ions. The value of Em in PathII‑III was reduced by ∼26% when a La ion substituted the Ce ion, whereas Em was increased when a Sc ion substituted the Ce ion. This result is consistent with our tetrahedron-volume model (see Figure S3), and is in general



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsami.7b13440. Computational details and further analysis on oxygen-ion migration in a highly doped CeO2 structure (PDF)



AUTHOR INFORMATION

Corresponding Author

*H. Kim. E-mail: [email protected]. D

DOI: 10.1021/acsami.7b13440 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Letter

ACS Applied Materials & Interfaces ORCID

Structures of Epitaxial Colossal Magnetoresistive La0.8Ca0.2MnO3 Thin Films. Appl. Phys. Lett. 1998, 73, 3294−3296. (15) Wang, T.; Ganguly, K.; Marshall, P.; Xu, P.; Jalan, B. Critical Thickness and Strain Relaxation in Molecular Beam Epitaxy-Grown SrTiO3 Films. Appl. Phys. Lett. 2013, 103, 212904. (16) Nakayama, M.; Martin, M. First-Principles Study on Defect Chemistry and Migration of Oxide Ions in Ceria Doped with RareEarth Cations. Phys. Chem. Chem. Phys. 2009, 11, 3241−3249. (17) Ahamer, C.; Opitz, A. K.; Rupp, G. M.; Fleig, J. Revisiting the Temperature Dependent Ionic Conductivity of Yittria Stabilized Zirconia (YSZ). J. Electrochem. Soc. 2017, 164, F790−F803. (18) Alaydrus, M.; Sakaue, M.; Aspera, S. M.; Wungu, T. D. K.; Linh, N. H.; Linh, T. P. T.; Kasai, H.; Ishihara, T.; Mohri, T. A DFT Plus U Study of Strain-Dependent Ionic Migration in Sm-Doped Ceria. J. Phys. Soc. Jpn. 2014, 83, 094707. (19) De Souza, R. A.; Ramadan, A.; Horner, S. Modifying The Barriers For Oxygen-Vacancy Migration in Fluorite-Structured CeO2 Electrolytes through Strain: a Computer Simulation Study. Energy Environ. Sci. 2012, 5, 5445−5453. (20) Avila-Paredes, H. J.; Choi, K.; Chen, C. T.; Kim, S. DopantConcentration Dependence of Grain-Boundary Conductivity in Ceria: A Space-Charge Analysis. J. Mater. Chem. 2009, 19, 4837−4842. (21) Chen, L.; Chen, C. L.; Donner, W.; Liu, S. W.; Lin, Y.; Huang, D. X.; Jacobson, A. J.; Chen, X. Electrical Properties of a Highly Oriented, Textured Thin Film of The Ionic Conductor Gd:CeO(2‑δ) on (001) MgO. Appl. Phys. Lett. 2003, 83, 4737−4739. (22) Joo, J. H.; Choi, G. M. Electrical Conductivity of Thin Film Ceria Grown by Pulsed Laser Deposition. J. Eur. Ceram. Soc. 2007, 27, 4273−4277. (23) Kilner, J. A.; Brook, R. J. A Study of Oxygen Ion Conductivity in Doped Non-Stoichiometric Oxides. Solid State Ionics 1982, 6, 237− 252. (24) Kim, N.; Kim, B. H.; Lee, D. Effect of Co-Dopant Addition on Properties of Gadolinia-Doped Ceria Electrolyte. J. Power Sources 2000, 90, 139−143. (25) Omar, S.; Wachsman, E. D.; Jones, J. L.; Nino, J. C. Crystal Structure-Ionic Conductivity Relationships in Doped Ceria Systems. J. Am. Ceram. Soc. 2009, 92, 2674−2681. (26) Shannon, R. D. Revised Effective Ionic-Radii and Systematic Studies of Interatomic Distances in Halides and Chalcogenides. Acta Crystallogr., Sect. A: Cryst. Phys., Diffr., Theor. Gen. Crystallogr. 1976, 32, 751−767. (27) Jaiswal, N.; Upadhyay, S.; Kumar, D.; Parkash, O. Ionic Conductivity Investigation in Lanthanum (La) and Strontium (Sr) Co-Doped Ceria System. J. Power Sources 2013, 222, 230−236. (28) Kahlaoui, M.; Chefi, S.; Inoubli, A.; Madani, A.; Chefi, C. Synthesis and Electrical Properties of Co-Doping with La3+, Nd3+, Y3+, and Eu3+ Citric Acid-Nitrate Prepared Samarium-Doped Ceria Ceramics. Ceram. Int. 2013, 39, 3873−3879. (29) Kobi, S.; Jaiswal, N.; Kumar, D.; Parkash, O. Ionic Conductivity of Nd3+ and Y3+ Co-Doped Ceria Solid Electrolytes for Intermediate Temperature Solid Oxide Fuel Cells. J. Alloys Compd. 2016, 658, 513− 519. (30) Yeh, T. H.; Chou, C. C. Ionic Conductivity Investigation in Samarium and Strontium Co-Doped Ceria System. Phys. Scr. 2007, T129, 303−307.

Kyung Joong Yoon: 0000-0002-4161-5111 Ji-Won Son: 0000-0002-5310-0633 Jong-Ho Lee: 0000-0003-4481-6258 Ho Won Jang: 0000-0003-4074-4277 Hyoungchul Kim: 0000-0003-3109-660X Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This research was financially supported by the New & Renewable Energy Core Technology Program of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) funded by the Ministry of Trade, Industry & Energy, Republic of Korea (No. 20143030031430). This work also supported in part by the Energy Efficiency & Resources Core Technology Program of the KETEP granted financial resource from the Ministry of Trade, Industry & Energy, Republic of Korea (No. 20152020106100).



REFERENCES

(1) Elam, J. W.; Dasgupta, N. P.; Prinz, F. B. ALD for Clean Energy Conversion, Utilization, and Storage. MRS Bull. 2011, 36, 899−906. (2) 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 2008, 321, 676−680. (3) Yildiz, B. "Stretching" The Energy Landscape of Oxides − Effects on Electrocatalysis and Diffusion. MRS Bull. 2014, 39, 147−156. (4) Sanna, S.; Esposito, V.; Andreasen, J. W.; Hjelm, J.; Zhang, W.; Kasama, T.; Simonsen, S. B.; Christensen, M.; Linderoth, S.; Pryds, N. Enhancement of The Chemical Stability in Confined δ-Bi2O3. Nat. Mater. 2015, 14, 500−504. (5) Shi, Y.; Bork, A. H.; Schweiger, S.; Rupp, J. L. M. The Effect of Mechanical Twisting on Oxygen Ionic Transport in Solid-State Energy Conversion Membranes. Nat. Mater. 2015, 14, 721−727. (6) Fluri, A.; Pergolesi, D.; Roddatis, V.; Wokaun, A.; Lippert, T. In Situ Stress Observation in Oxide Films and How Tensile Stress Influences Oxygen Ion Conduction. Nat. Commun. 2016, 7, 10692. (7) Thieu, C. A.; Hong, J.; Kim, H.; Yoon, K. J.; Lee, J. H.; Kim, B. K.; Son, J. W. Incorporation of a Pd Catalyst at The Fuel Electrode of a Thin-Film-Based Solid Oxide Cell by Multi-Layer Deposition and Its Impact on Low-Temperature Co-Electrolysis. J. Mater. Chem. A 2017, 5, 7433−7444. (8) Cho, A. Y.; Ballamy, W. C. GaAs Planar Technology by Molecular-Beam Epitaxy (MBE). J. Appl. Phys. 1975, 46, 783−785. (9) Cavallaro, A.; Burriel, M.; Roqueta, J.; Apostolidis, A.; Bernardi, A.; Tarancón, A.; Srinivasan, R.; Cook, S. N.; Fraser, H. L.; Kilner, J. A.; McComb, D. W.; Santiso, J.; et al. Electronic Nature of The Enhanced Conductivity in YSZ-STO Multilayers Deposited by PLD. Solid State Ionics 2010, 181, 592−601. (10) Amezawa, K.; Kushi, T.; Sato, K.; Unemoto, A.; Hashimoto, S.; Kawada, T. Elastic Moduli of Ce0.9Gd0.1O2‑δ at High Temperatures under Controlled Atmospheres. Solid State Ionics 2011, 198, 32−38. (11) Atkinson, A.; Selcuk, A. Mechanical Behaviour of Ceramic Oxygen Ion-Conducting Membranes. Solid State Ionics 2000, 134, 59− 66. (12) Jiang, C.; Wolverton, C.; Sofo, J.; Chen, L. Q.; Liu, Z. K. FirstPrinciples Study of Binary BCC Alloys Using Special Quasirandom Structures. Phys. Rev. B: Condens. Matter Mater. Phys. 2004, 69, 214202. (13) Zunger, A.; Wei, S. H.; Ferreira, L. G.; Bernard, J. E. Special Quasirandom Structures. Phys. Rev. Lett. 1990, 65, 353−356. (14) Rao, R. A.; Lavric, D.; Nath, T. K.; Eom, C. B.; Wu, L.; Tsui, F. Three-Dimensional Strain States and Crystallographic Domain E

DOI: 10.1021/acsami.7b13440 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX