Forums
New posts
Search forums
What's new
New posts
Latest activity
Articles
Members
Current visitors
Log in
Register
What's new
Search
Search
Search titles only
New posts
Search forums
Menu
Log in
Register
Install the app
Install
Forums
Lab talk, the good, the bad, and the ugly
Zirconium
not easy to explain. or understand.. or accept...
JavaScript is disabled. For a better experience, please enable JavaScript in your browser before proceeding.
You are using an out of date browser. It may not display this or other websites correctly.
You should upgrade or use an
alternative browser
.
Reply to thread
Message
<blockquote data-quote="disturbed" data-source="post: 81213" data-attributes="member: 4763"><p>This ones going to be difficult to explain. I have some more information to regurgitate(love that)..</p><p></p><p>I am now in contact with chemists at the local college and talking to Ivoclar to get the specifics on the molecular exchange between Zirliner and Zirconium, as I have proof that IF there is one it is VERY weak. more likely there is not.</p><p></p><p>FROM THE JOURNAL OF MATERIALS CHEMISTRY accepted march 2010 sooo... recent</p><p>this is complicated info..start with Section IV (fig 4) Migration barrier as a function of biaxial lattice strain. than read conclusion.</p><p></p><p> this is above a lot of peoples understanding so I will be producing studies from colleges across the globe and doing my best to explain the chemistry over the next few months...</p><p></p><p></p><p>and yes.. its another attempt to show you all that a bond other than mechanical retention is IMPOSSIBLE to zirconium and that their will be more chipping and failures to come. if you only care about your restorations lasting a max of 5 years than go ahead and ignore this.. it will probably only confuse and uspet you..</p><p></p><p>If I find that materials are added to zirconium oxide that allow for a chemical bond that gives an adequate bond strength I will post.. but as the vitrification of the material happens around 2700 degrees C I doubt it is possible.</p><p></p><p>[ATTACH]3849[/ATTACH]</p><p></p><p></p><p></p><p>Oxygen ion diffusivity in strained yttria stabilized zirconia: where is the</p><p>fastest strain?</p><p>Akihiro Kushima and Bilge Yildiz*</p><p>Received 8th January 2010, Accepted 9th March 2010</p><p>First published as an Advance Article on the web 29th April 2010</p><p>DOI: 10.1039/c000259c</p><p>We present the mechanism and the extent of increase in the oxygen anion diffusivity in Y2O3 stabilized</p><p>ZrO2 (YSZ) under biaxial lattice strain. The oxygen vacancy migration paths and barriers in YSZ as</p><p>a function of lattice strain was assessed computationally using density functional theory (DFT) and</p><p>nudged elastic band (NEB) method. Two competing and non-linear processes acting in parallel were</p><p>identified to alter the migration barrier upon applied strain: (1) the change in the space, or electronic</p><p>density, along the migration path, and (2) the change in the strength of the interatomic bonds between the</p><p>migrating oxygen and the nearest neighbor cations that keep the oxygen from migrating. The increase of</p><p>the migration space and the weakening of the local oxygen–cation bonds correspond to a decrease of the</p><p>migration barrier, and vice versa. The contribution of the bond strength to the changes in the migration</p><p>barrier is more significant than that of the opening of migration space in strained YSZ. A database of</p><p>migration barrier energies as a function of lattice strain for a set of representative defect distributions in</p><p>the vicinity of the migration path inYSZwas constructed. This database was used in kinetic Monte Carlo</p><p>(KMC) simulations to estimate the effective oxygen diffusivity in strained YSZ. The oxygen diffusivity</p><p>exhibits an exponential increase up to a critical value of tensile strain, or the fastest strain. This increase is</p><p>more significant at the lower temperatures. At the strain states higher than the critical strain, the</p><p>diffusivity decreases. This is attributed to the local relaxations at large strain states beyond a limit of</p><p>elastic bond strain, resulting in the strengthening of the local oxygen–cation bonds that increases the</p><p>migration barrier. The highest enhancement of diffusivity in 9%-YSZ compared to its unstrained state is</p><p>6.8 103 times at 4% strain and at 400 K. The results indicate that inducing an optimal strain state by</p><p>direct mechanical load or by creating a coherent hetero-interface with lattice mismatch can enable</p><p>desirably high ionic conductivity in YSZ at reduced temperatures. The insights gained here particularly</p><p>on the nonlinear and competing consequences of lattice strain on the local bonding structure and charge</p><p>transport process are of importance for tuning the ionic transport properties in a variety of solid-state</p><p>conducting material applications, including but not limited to fuel cells.</p><p>Introduction</p><p>There is much interest in decreasing the Solid Oxide Fuel Cell</p><p>(SOFC) operating temperature to an intermediate-to-low</p><p>temperature range to improve chemical and mechanical stability</p><p>and reduce the cost. Enhancing the ionic transport in the electrolyte</p><p>and cathode materials is a key for this purpose. Over the</p><p>past several decades, developing new compositions with high</p><p>ionic mobility at low temperatures has been the primary activity.1</p><p>In recent years, however, nanoscale structures were shown as</p><p>means for potentially increasing the transport properties in</p><p>oxides of relevance to SOFCs.2–9 The importance of utilizing the</p><p>favorable interface properties of nanoscale structures shifted the</p><p>paradigm from developing new material compositions to</p><p>tailoring nano-structures using existing materials for the above</p><p>purpose.5,10 Among these, nanoscale hetero-layered oxides were</p><p>suggested to enhance the ionic conductivity along the interfaces</p><p>by orders of magnitude.5–9 The presence of an interface between</p><p>dissimilar oxides is commonly the source of charged mobile</p><p>defects and local anisotropic distortions and strained interatomic</p><p>bonds. These contribute to changes in the local electronic</p><p>structure; formation, coordination and redistribution of charged</p><p>defects; and mobility of defects in the vicinity of interfaces. Most</p><p>recently, an eight orders of magnitude increase in ionic conductance</p><p>was reported by Barriocanal et al.8 for 1–30 nm-thick yttria</p><p>stabilized zirconia (YSZ) layers coherently ‘‘strained’’ between</p><p>dielectric SrTiO3 (STO) layers. The exact nature of the ionic vs.</p><p>electronic conductance induced in the vicinity of the YSZ/STO</p><p>interface remains debatable.11 While these observations imply</p><p>that the hetero-interfaces with controlled lattice-strain and</p><p>defect-chemistry can play an important role in improving the</p><p>ionic mobility, results to date remain phenomenological without</p><p>systematic theoretical efforts. Beyond the SOFC materials</p><p>domain, similar questions for the role of hetero-interfaces on the</p><p>local bonding structure and the charge transport process are</p><p>open for tuning the transport properties in a variety of solid-state</p><p>conducting material applications.12–15 Systematic experimental</p><p>and theoretical investigation of the role of local strain states on</p><p>tuning the mechanism and magnitude of charge transfer properties</p><p>in these materials is prone to pursuit.</p><p>Department of Nuclear Science and Engineering, Massachusetts Institute</p><p>of Technology, 77 Massachusetts Avenue, Cambridge, MA, 02139, USA.</p><p>E-mail: <a href="mailto:byildiz@mit.edu">byildiz@mit.edu</a>; Tel: +617-324-4009</p><p>This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4809</p><p>PAPER <a href="http://www.rsc.org/materials" target="_blank">Journal of Materials Chemistry Home-High impact applications, properties and synthesis of exciting new materials</a> | Journal of Materials Chemistry</p><p>Two main mechanisms contribute to increase the ionic transport</p><p>at the hetero-interface: (1) a favorable strain state at the</p><p>interface to shift and/or change the symmetry of electron energy</p><p>levels to provide for improved charge transfer16 and ion</p><p>mobility,17 and (2) the alteration of the defect chemistry7,18 near</p><p>the interface to enhance the density and distribution of desired</p><p>charge carriers. To understand the underlying mechanism of</p><p>oxygen anion transport in anion conducting ceramics, a significant</p><p>number of theoretical studies using atomistic simulations</p><p>based on both density functional theory (DFT) and empirical</p><p>potentials have been conducted (examples related to this work in</p><p>ref. 19–25). However, most of them are focused on the ‘‘bulk’’ or</p><p>the ‘‘surface’’ ionic transport. Only few theoretical studies exist</p><p>for the anion transport at hetero-interfaces, on specific cases,</p><p>i.e. CaF2/BaF2, CeO2/YSZ,20 with no generalization of results.</p><p>Of particular interest for this paper is the effect of lattice strain</p><p>on the ionic conductivity in the YSZ, as may be encountered due</p><p>to a lattice mismatch near the YSZ/STO interface. In probing the</p><p>effect of strain on conductivity in YSZ, a number of experimental</p><p>and simulation studies have been conducted in the past.21–24 Even</p><p>though the past results imply a potential increase in the</p><p>conductivity, the magnitude of increase was far lower than what</p><p>was reported recently in ref. 8 for the YSZ/STO interface.</p><p>Furthermore, the reported increases in magnitude of the</p><p>enhancement are inconsistent among these past studies because</p><p>of the differences in conditions, structures and methods that were</p><p>considered. For example, results were based on examining either</p><p>incoherent hetero-interfaces which cannot sustain large lattice</p><p>strain,21 or low strain conditions,24 or molecular dynamics</p><p>simulation with empirical potentials23,24 which are not appropriate</p><p>for interfaces, or phenomenological modeling using</p><p>macroscopic elastic properties of the material.22 In essence, the</p><p>mechanism by which a biaxial lattice strain alters the anion</p><p>transport in the YSZ fluorite structure was not uncovered. Given</p><p>these past inconsistencies and open questions, we adopt a firstprinciples-</p><p>based approach in probing the mechanistic picture of</p><p>the ionic conduction in strained YSZ. Our specific objectives are</p><p>to elucidate the mechanism and to identify the possible extent of</p><p>increase in the oxygen ion conductivity in YSZ due to a biaxial</p><p>lattice strain.</p><p>In order to attain reliable results in probing the effect of strain</p><p>on ionic transport in YSZ, it is important to consider a realistic</p><p>representation of the defect structure, in particular the distribution</p><p>of vacancies and dopant cations in YSZ. This is because the</p><p>oxygen vacancy migration barrier, which can be altered by strain,</p><p>depends strongly on the local defect distributions.25 On the other</p><p>hand, it is difficult to characterize experimentally the vacancy</p><p>and cation distribution in YSZ at atomic resolution, therefore</p><p>this type of information is not available from literature</p><p>straightforwardly. In spite of the difficulties, past work on this</p><p>topic provided useful insights on the vacancy–vacancy, vacancy–</p><p>cation and cation–cation interactions, enabling a reasonable</p><p>description of the defect structure in YSZ.</p><p>Here we provide a summary of defect–defect interactions that</p><p>we considered for constructing the YSZ model. The final YSZ</p><p>model used in our work should reasonably capture these reported</p><p>defect–defect interactions. First of these interactions</p><p>considered, namely the vacancy–cation interaction, determines</p><p>the relative position of the vacancy with respect to the Zr and Y</p><p>in YSZ. In earlier studies for this interaction, there were contradicting</p><p>experimental results; whether the vacancy is first</p><p>nearest neighbor26–28 (1NN) or second nearest neighbor29–31</p><p>(2NN) to the Y cations. However, later experimental30–33 and</p><p>theoretical34–40 studies indicated that the vacancy favors to bind</p><p>to the host Zr cations, resulting in a configuration such that the</p><p>vacancy is 1NN to the Zr and 2NN to the Y cations. Using DFT</p><p>calculations, Bogicevic and Wolverton37 showed that the local</p><p>relaxation of atoms was responsible for the vacancy to favor the</p><p>2NN position to the Y cation; opposite to what is expected from</p><p>the electrostatic interactions alone. The same trend was reported</p><p>in doped CeO2 when the dopant radius is larger than the host</p><p>cation,41 as in YSZ. For the vacancy–vacancy and the cation–</p><p>cation interactions, they concluded that the vacancies locate</p><p>themselves as far away as possible from each other, and that the</p><p>1NN Y–Y pair configurations are slightly more favorable than</p><p>the 2NN Y–Y pairs.36,38 More recent experimental work using</p><p>selected area diffraction also indicated a random distribution of</p><p>oxygen vacancies, maximizing the distance between them, in up</p><p>to 10% yttria doped YSZ.42 Consistent with the findings of</p><p>Bogicevic and Wolverton, theoretical results by Predith et al.43</p><p>suggest even the existence of the ordering of the dopant cations</p><p>(Y) in YSZ at high doping concentrations, and this cation</p><p>ordering could exist locally in YSZ at the lower Y concentration</p><p>range. Among these three types of defect–defect interactions, the</p><p>vacancy–vacancy interaction is the strongest followed by the</p><p>vacancy–cation interaction, and the cation–cation interaction is</p><p>the weakest of the three.36 Supporting the importance of the</p><p>vacancy–vacancy interactions, it was shown recently that the</p><p>vacancy distribution around the migration path could significantly</p><p>impact the value of the migration barriers in YSZ.44</p><p>In this study, we considered only the strain state in the bulk</p><p>to provide for improved ionic mobility in YSZ and did not</p><p>include an interface explicitly as in Barriocanal’s work.8</p><p>Therefore, the model in this paper excludes the effect of altered</p><p>defect chemistry near the interfaces of YSZ as a source of</p><p>increase in oxygen diffusivity. In the following sections, we</p><p>describe the simulation approach to calculate the oxygen</p><p>vacancy migration path and energy barriers. We report the</p><p>energetically favorable distributions of oxygen vacancies and Y</p><p>cations in the YSZ model structure to enable a realistic representation</p><p>of YSZ in this study. The mechanism by which the</p><p>lattice strain increases or decreases the anion transport in YSZ</p><p>are discussed. Because the vacancy migration barrier depends</p><p>on the local defect configurations, a database of migration</p><p>barrier energies as a function of lattice strain for a set of</p><p>representative defect distributions in the vicinity of the migration</p><p>path are constructed and presented. Finally, the effective</p><p>tracer diffusivity of oxygen in YSZ at different strain states is</p><p>presented and interpreted in comparison to the recent results</p><p>reported for the YSZ/STO interface.8</p><p>Simulation approach</p><p>The YSZ model was established to capture the charged defect</p><p>interactions and the defect distribution in the bulk YSZ, as</p><p>summarized in the Introduction. Upon identifying the stable</p><p>vacancy–cation distribution, the effect of the biaxial lattice</p><p>strain in the xy-plane on the migration barrier and on the</p><p>4810 | J. Mater. Chem., 2010, 20, 4809–4819 This journal is ª The Royal Society of Chemistry 2010</p><p>effective vacancy diffusivity were calculated and interpreted.</p><p>We used first-principles-based simulations to directly and</p><p>accurately probe the atomistic-scale nature of the vacancy</p><p>migration in the presence of lattice strain. The kinetic Monte</p><p>Carlo (KMC) method was used for calculating the effective</p><p>diffusivities.</p><p>The first principles calculations for identifying the energies</p><p>were performed by the Vienna ab initio simulation package</p><p>(VASP),45,46 employing density functional theory (DFT) using</p><p>a plane-wave basis set. Projector-augmented wave (PAW)</p><p>method47,48 with plane waves up to the energy cutoff at 400 eV</p><p>was used and the exchange–correlation energy was evaluated</p><p>by the generalized gradient approximation (GGA) using the</p><p>Perdew–Burke–Ernzerhof (PBE) function.49 Y, Zr and O atoms</p><p>were described by 11 (4s24p65s24d1),12 (4s24p65s24d2),and 6</p><p>(2s22p4) valence electrons, respectively. The G point in the Brillouin</p><p>zone was selected in the calculation. We confirmed the</p><p>energy convergence by comparing the energy obtained from each</p><p>of the k-point meshes of 3 3 3, 5 5 5 and 7 7 7. The</p><p>calculated energies are within 0.6% error compared to the energy</p><p>obtained using only the G point. The same condition was applied</p><p>to all of the simulations in this study unless stated otherwise.</p><p>YSZ has the fluorite crystal structure with the cations</p><p>(Y3+, Zr4+) occupying the FCC lattice sites and the O2 ions</p><p>occupying the tetragonal sites. For every two Y atoms, one</p><p>oxygen vacancy is created to satisfy charge neutrality. The</p><p>oxygen vacancy migrates by hopping to the adjacent tetragonal</p><p>site as schematically illustrated in Fig. 1. The migration barrier</p><p>depends on the neighboring atoms 1–6 in the figure.25 The</p><p>surrounding cations within a larger cutoff radius to the migration</p><p>path must also be accounted for in order to accurately quantify</p><p>the effect of the strain on the migration barrier.50 However, this</p><p>leads to an impractical number of patterns (well above 10 000</p><p>configurations just by including the next nearest neighbor cations</p><p>to the migration path) to be handled by the DFT simulation.</p><p>Therefore, we selected a subset of defect structure patterns near</p><p>the migration path, forming a database of migration paths and</p><p>barriers. This database serves as input for the KMC calculations</p><p>in estimating the effective diffusivity as a function of lattice strain</p><p>and temperature. In constructing this database, first, we</p><p>considered only two combinations of atoms at the sites 3 and 4,</p><p>as Zr–Zr and Zr–Y, in the migration path shown in Fig. 1. The</p><p>migration barrier through the Y–Y pair was found to be significantly</p><p>higher than these two combinations, and this makes the</p><p>Y–Y pair unlikely on a vacancy migration path. This selection</p><p>was also validated by the work of Car et al.,50 which implied the</p><p>Zr–Zr and Zr–Y subset at sites 3 and 4 to be sufficient to capture</p><p>the quantitative trends in ionic conductivity in YSZ as a function</p><p>of Y2O3 doping concentration. Second, we considered the defect</p><p>configurations outside of the first neighbor cations to the</p><p>migration path to affect the barrier. The particularly important</p><p>impact of the vacancy interactions on the migration barrier was</p><p>recently reported by Pietrucci et al.44 Therefore, in creating the</p><p>migration path and barrier database, we varied the relative</p><p>positions of the Y cations and the vacancies with respect to the</p><p>initial and final positions of the vacancy on the migration path.</p><p>The details of the migration path configurations in this database</p><p>are provided in the following section. For each migration path</p><p>configuration, nudged elastic band (NEB) method51 was used to</p><p>find the vacancy migration barrier as a function of biaxial strain</p><p>on the xy-plane. Given that the vacancy migration path in YSZ is</p><p>rather simple with single hops between the adjacent tetragonal</p><p>sites, the saddle point is expected near the midst of the migration</p><p>path.25,50 Three images between the initial and the final states</p><p>were taken in the simulations. Here we used the climbing image</p><p>NEB method52 to calculate not only the minimum energy path</p><p>but also the transition state configuration at the saddle point.</p><p>The method drives the image with the highest energy up to the</p><p>saddle point, by maximizing its energy along the band connected</p><p>to the image and minimizing in all other directions. This way, the</p><p>method allows us to accurately calculate the saddle point energy</p><p>and configuration with a small number of images needed in the</p><p>calculation. The NEB simulations were conducted with fixed</p><p>lattice vectors.</p><p>The database created for the migration barrier values at</p><p>different lattice strains and the defect configurations was</p><p>provided as input for kinetic Monte Carlo (KMC) simulations.</p><p>KMC calculations were performed to quantify the effective</p><p>diffusivity of oxygen in YSZ as a function of strain and</p><p>temperature. The diffusivity calculated in this work is the oxygen</p><p>tracer diffusion constant. We expect that the lattice strain has</p><p>a similar order of magnitude impact on the mobility of an oxygen</p><p>ion in both the chemical and the tracer diffusion constants,</p><p>especially in the range of low vacancy concentrations.53 In the</p><p>simulation model, oxygen atoms were assumed to hop to the</p><p>nearest neighbor site on the oxygen sublattice through the path</p><p>shown in Fig. 1. The rate of hops in the migration paths were</p><p>obtained by the Boltzman relationship as,</p><p>vAB ¼ n0 expEAB</p><p>kBT (1)</p><p>where EAB, found by the DFT-NEB simulations, denotes the</p><p>migration barrier across the A–B cation pair at sites 3 and 4 in</p><p>Fig. 1. The attempt frequency was fixed at n0 ¼ 1013 s1, which is</p><p>appropriate for most metal oxide systems.50 The periodic</p><p>boundary conditions were applied to all three dimensions of the</p><p>simulation cell. The vacancy diffusivity, Dv was extracted from</p><p>mean square displacement hR2i calculated in the KMC simulation</p><p>using the Einstein relation,</p><p>hR2i ¼ 6Dvt. (2)</p><p>The oxygen diffusivity DO was then obtained using Dv</p><p>considering the balance between the fluxes of oxygen atoms and</p><p>vacancies as,</p><p>Fig. 1 Migration path of the oxygen vacancy. Positions 1–6 are occupied</p><p>by the cations Zr or Y.</p><p>This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4811</p><p>DO ¼</p><p>cv</p><p>1 cv</p><p>Dv (3)</p><p>where cv is the vacancy concentration fraction. Dv exponentially</p><p>depends on the effective barrier for oxygen vacancy migration,</p><p>Eeff, as,</p><p>Dv ¼ dv</p><p>0 exp </p><p>Eeff</p><p>kBT (4)</p><p>where, d0v</p><p>is a constant. At small values of cv, D0 scales almost</p><p>linearly with cv. On the other hand, the change in the migration</p><p>barrier, Eeff, contributes exponentially to Dv and D0.</p><p>Results and discussion</p><p>This section consists of five sub-sections, I–V. Sub-sections I–II</p><p>demonstrate that the YSZ model in this work captures the</p><p>reported vacancy and cation interactions in stabilizing the cubic</p><p>structure in YSZ. The results in I and II are important for</p><p>obtaining reliable values of the energy barriers for vacancy</p><p>migration in YSZ. Sub-section III reports the effect of defect</p><p>configurations in the vicinity of the migrating vacancy on the</p><p>migration energy barrier. This proves that, in order to assess the</p><p>effect of strain on migration barriers, one must consider a realistic</p><p>subset of defect configurations in YSZ, and we report these</p><p>in sub-section IV. Sub-section IV in particular illustrates the</p><p>mechanism by which strain alters the diffusion barriers for</p><p>oxygen vacancies in YSZ. Finally sub-section V quantifies the</p><p>effective increase in oxygen tracer diffusivity as a function of</p><p>biaxial lattice strain, based on the subset of defect configurations</p><p>considered in this work for YSZ.</p><p>I. Defect structure in YSZ: oxygen vacancy and yttrium</p><p>distribution</p><p>To determine the favorable location of vacancies with respect to</p><p>the Y cations in the YSZ model, we tested a number of defect</p><p>configurations. These models constitute 30 Zr, 2 Y and 63 O</p><p>atoms, corresponding to a low doping ratio, 3% Y2O3, allowing</p><p>only one vacancy in the simulation cell, in order to eliminate the</p><p>vacancy–vacancy interaction in this part of the simulations. At</p><p>this low Y2O3 doping, the YSZ oxygen sublattice cannot keep its</p><p>cubic structure.34 Therefore, only in this part of the simulations,</p><p>the atoms beyond the second neighbor cation shell to Y cations</p><p>and beyond the third neighbor oxygen shell to the vacancy were</p><p>fixed to retain the cubic phase.34 The comparison of the energies</p><p>for the different YSZ models and defect configurations is presented</p><p>in Fig. 2. The energies are shifted such that the lowest</p><p>energy found is zero. The results show that the total energy can</p><p>vary by up to 0.5 eV when the Y cation is located at the third or</p><p>first neighbor site to the vacancy, instead of the second. This</p><p>result is consistent with the prior reports on cation–vacancy</p><p>interactions in YSZ, as discussed in the Introduction. Since this</p><p>large energy difference arises from local configuration differences,</p><p>the result suggests that the cation configuration</p><p>surrounding the diffusion path of the oxygen–vacancy pair</p><p>should be taken into account for accurate evaluation of the</p><p>migration barriers.</p><p>II. Distortion from the cubic lattice sites in YSZ</p><p>The atoms are not stabilized at the exact fluorite structure lattice</p><p>centers in YSZ, and large displacements of the cations and the</p><p>anions from the fluorite lattice sites are present.55–57 This is in</p><p>part due to the relaxation of the nearest neighbor oxygen anions</p><p>and cations towards and away from the vacancy, respectively.32</p><p>For an acceptable YSZ structure in this study, it is important to</p><p>capture this deviation quantitatively consistently with past</p><p>experimental and simulation reports. For this purpose, we</p><p>evaluated the dependence of the atomic displacements from the</p><p>ideal fluorite lattice sites on the local defect structures in a set of</p><p>candidate YSZ models. The fully relaxed structure of the models</p><p>constituted 9% Y2O3 doped YSZ, made of 26 Zr, 6 Y and 61 O</p><p>atoms in the simulation cell. The relaxation of the atoms were not</p><p>constrained in this and the following simulations because the</p><p>cubic structure is stable in the 9% Y2O3 doped YSZ, as long as</p><p>Fig. 2 Energy as a function of the neighboring location of the oxygen</p><p>vacancy with respect to the Y cations in YSZ. The energies were calculated</p><p>for different cation/vacancy distributions, and plotted such that the</p><p>lowest energy found is zero.</p><p>Fig. 3 Fully relaxed atomic configurations of the 9% yttria doped</p><p>zirconia with (a) nearest neighbor Y–Y pairs, and (b) uniformly</p><p>distributed Y. Yellow, gray and red spheres indicate Y, Zr and O ions.</p><p>Large distortion of the oxygen anion sublattice is evident when Y cations</p><p>are uniformly spread in the YSZ without any short-range ordering of the</p><p>Y cations.</p><p>4812 | J. Mater. Chem., 2010, 20, 4809–4819 This journal is ª The Royal Society of Chemistry 2010</p><p>the vacancies and dopant cations are distributed to enable the</p><p>stability. Fig. 3 shows the relaxed atomic configurations for two</p><p>models of YSZ. The first one (Fig. 3(a)) has first nearest neighbor</p><p>(1NN) Y–Y pairs and the second one (Fig. 3(b)) has a uniform</p><p>distribution of Y cations in the cell. In both models, the vacancies</p><p>were placed as far away as possible from each other and located</p><p>second nearest neighbor (2NN) to Y cations if allowed by the</p><p>available space. The first model retains an overall cubic structure</p><p>with oxygen anions displaced by 0.1–0.27 A from the ideal</p><p>fluorite lattice sites. The extent of the displacements calculated</p><p>here is well within the experimentally measured value of</p><p>0.36 A28 and prior simulation results of 0.20–0.31 A.34 On the</p><p>other hand, the second model shows large distortions, exemplified</p><p>at the circled region in Fig. 3(b),and has 1.2 eV higher energy</p><p>compared to the first one. This difference arises due to local</p><p>configuration differences, and thus, can affect the vacancy</p><p>migration barrier in YSZ which ranges from 0.2 to 1.4 eV.25</p><p>The lack of stability of the cubic oxygen sublattice in the second</p><p>model is due to the limited available space to accommodate the</p><p>relative positions of the vacancy and Y cations at the favorable</p><p>distances when the Y atoms are distributed uniformly into the</p><p>model. In the circled region in Fig. 3(b),a vacancy is located at</p><p>the 1NN position to the Y cations. This is because the 2NN</p><p>position for the vacancy is not available to all the Y atoms when</p><p>the Y atoms are distributed uniformly in the model. The</p><p>restriction of some vacancies to the unfavorable 1NN positions</p><p>with respect to the Y in YSZ destabilizes the cubic phase. On the</p><p>contrary, the first model configuration with the Y–Y pairs</p><p>introduces more space for the vacancies to be placed at the 2NN</p><p>positions with respect to the Y cations. This enables more</p><p>stability on the cubic oxygen sublattice. This result is consistent</p><p>with the locally ordered nature of Y cations in YSZ at high</p><p>concentrations of Y2O3, as found by Predith et al.,43 even if the</p><p>YSZ has only 9% Y2O3, here. While the effect of Y–Y pairing on</p><p>the migration barrier in a cubic YSZ is small compared to that of</p><p>the vacancy–vacancy and vacancy–cation interactions,36 it is</p><p>important to consider it to enable a stabilized cubic fluorite</p><p>structure of YSZ in the model.</p><p>III. Vacancy migration barriers in YSZ: effect of the local</p><p>defect structure</p><p>As implied in the previous section, the local defect structure</p><p>around the oxygen vacancy migration path in YSZ affects the</p><p>stability of the structure and the migration barrier,54 and this in</p><p>turn can determine the favorable paths for vacancy to migrate</p><p>through. To illustrate this, we present the differences in the</p><p>migration barrier in a model with 1NN Y–Y pairs, which lower</p><p>the energy and enable the stabilization of cubic phase. The</p><p>consequent migration paths considered were A / B / C in</p><p>Fig. 4(a). All six nearest neighbor cations (1–6 as shown in Fig. 1)</p><p>are the same, Zr, for each migration step in A/B/C. On the</p><p>other hand, the distribution of vacancies and Y cations near A</p><p>/B and B/C differ beyond the nearest neighbor cations. This</p><p>enables to investigate the variation of the migration barrier as</p><p>a function of the defect structure near the migration path beyond</p><p>the nearest neighbor cations. Fig. 4(b) shows the relative energy</p><p>in the NEB trajectories for the migration through A/B/C.</p><p>The resulting migration barriers range from 0.55 to 0.95 eV, even</p><p>though all the cations in the nearest neighbor position to the</p><p>migration path are the same. This significant difference exemplified</p><p>in the migration barriers in A/B/C indicates that it is</p><p>important to systematically take into account the longer-range</p><p>interactions between the migrating vacancy and its surrounding</p><p>cations and vacancies in YSZ.</p><p>We evaluated a subset of defect configurations near the</p><p>migration path to form a representative database of migration</p><p>barriers in YSZ. The five different configurations, models A–E,</p><p>that yield ten forward and backward jump barriers considered</p><p>here are shown in Fig. 5. This subset aims to capture the isolated</p><p>Fig. 4 (a) The vacancy positions and local structure considered for the</p><p>vacancy migration path, A / B / C. Large black spheres indicate</p><p>positions of the vacancies. Yellow, gray and red spheres indicate Y, Zr</p><p>and O ions. Numbers of Y cations at the second nearest neighbor position</p><p>to the vacancy is two at sites A and B, and one at site C. (b) NEB barriers</p><p>calculated along the A / B and B / C paths. The variation in the</p><p>migration barrier is clear in (b),even though all the six nearest neighbor</p><p>cations in the migration path are Zr. The difference in barriers is due to</p><p>the dissimilar vacancy and Y distribution near the A / B and B / C</p><p>paths beyond the nearest neighbor cations.</p><p>Fig. 5 Models A–E with varying defect configurations surrounding the</p><p>migration path. ‘‘# of 2NNY’’ denotes the number of second nearest</p><p>neighbor Y cations to the vacancy positions. ‘‘NN vacancy position’’</p><p>denotes that the nearest neighbor vacancy to the migrating vacancy</p><p>is located at n-th nearest neighbor position on the anion sublattice</p><p>(n-th ¼ 3rd or 4th in these models). Lower panel indicates the 2NN Y</p><p>position, and the 3rd and 4th nearest neighbor vacancies on the anion</p><p>sublattice with respect to the migrating vacancy V.</p><p>This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4813</p><p>and coupled effects of vacancy–vacancy and vacancy–cation</p><p>interactions in the local environment around the migration path.</p><p>The vacancy–vacancy interaction is considered up to the fourth</p><p>nearest neighbor on the anion sublattice, which was the</p><p>maximum that could be accommodated in our 2 2 2 unit cell</p><p>model when the nearest neighbor vacancy pairs are avoided. The</p><p>cation–vacancy interaction was taken up to the second nearest</p><p>neighbor. First, for considering the effect of the nearest neighbor</p><p>cations in the migration path (sites 1–6 in Fig. 1),we only</p><p>included the variation of the bridging cation pairs (sites 3 and 4</p><p>in Fig. 1). This was shown to be the most important contribution</p><p>to the vacancy migration in YSZ50 when the long-range interactions</p><p>were not considered. In models A, B, and E, the vacancy</p><p>migrates through the Zr–Zr cation pair, and in models C and D,</p><p>through the Zr–Y cation pair in the midst of the path. All the</p><p>other cation sites closest to the migration path (sites 1, 2, 5 and</p><p>6 in Fig. 1) are Zr. The path through Y–Y pair was excluded from</p><p>the simulations since it has a much higher migration barrier</p><p>compared to the Zr–Zr and Zr–Y pairs,25,50,54 making it significantly</p><p>less likely to migrate through the Y–Y cation pair. Second,</p><p>for vacancy–vacancy interaction, models A and C consider the</p><p>position of nearest neighbor vacancies, denoted as ‘‘NN vacancy</p><p>position’’ in Fig. 5, to the initial and the final position of the</p><p>migrating vacancy. The distance between the migrating vacancy</p><p>and the nearest neighbor vacancy changes from the fourth to the</p><p>third on the anion sublattice upon migrating from the initial to</p><p>the final position (and vice versa in the backward migration). # of</p><p>2NN Y are the same in models A and C. Third, for vacancy–</p><p>cation interaction beyond the nearest neighbor cations, models B</p><p>and D consider the number of Y cations that are at the second</p><p>nearest neighbor (2NN) position with respect to the initial and</p><p>the final state of the vacancy. This is denoted as ‘‘# of 2NN Y’’ in</p><p>Fig. 5. The number of 2NN Y cations changes from two to three</p><p>upon migrating from the initial to the final state (and vice versa in</p><p>the backward migration). Lastly, in model E, both the number of</p><p>2NN Y cations and the distance to the nearest vacancies at the</p><p>initial and the final state were varied. The former changes from</p><p>three to two and the latter from the fourth to the third nearest</p><p>neighbor (and vice versa in the backward migration).</p><p>In creating the database of migration paths and barriers,</p><p>a model to represent a symmetric defect configuration around the</p><p>initial and final states of the migrating vacancy was not possible.</p><p>This is because of the computational limitation in the simulation</p><p>cell size used in our DFT/NEB calculations. In computing the</p><p>effective diffusivity in the following sections, we will show that</p><p>we represented the barrier for a symmetric migration path to be</p><p>the average of the barriers from models A–E.</p><p>A simulation unit cell of 2 2 2 dimensions, containing</p><p>26 Zr, 6 Y and 61 O atoms corresponding to a 9% yttria doped</p><p>zirconia was used in the calculation of the migration energy</p><p>barriers for each of the models A–E described in Fig. 5. The</p><p>oxygen vacancies and the Y cations were distributed in the cell</p><p>satisfying the defect interactions discussed in the previous</p><p>sections. Atomic configurations and the cell vectors were fully</p><p>relaxed such that the total energy was minimized and the net</p><p>stress acting on the cell was zero. The migration barriers for each</p><p>model path were calculated using the climbing NEB technique.</p><p>The calculated barriers for each path (A–E) are shown in Fig. 6.</p><p>The energy values are relative to the initial state configuration.</p><p>The reaction coordinate of each NEB image was normalized with</p><p>respect to the distance between the initial and the final state of the</p><p>vacancy. The oxygen vacancy migration process requires</p><p>climbing over the highest energy (the saddle point) between the</p><p>initial and the final state. E!B and E B in Fig. 6 are the migration</p><p>barriers for the forward and the backward jumps, respectively,</p><p>corresponding to the migration from the initial to the final</p><p>positions of the vacancy shown in Fig. 5, and vice versa.</p><p>As a result, the migration configurations in models C and D</p><p>have higher barrier than those in models A, B and E. This is</p><p>primarily because of the difficulty in passing through the Zr–Y</p><p>pair compared to that through Zr–Zr pair. The forward jump</p><p>barrier is higher than the backward jump barrier for models</p><p>A and C, and lower for models B and D. For models A and C,</p><p>this is because the distance from the migrating vacancy to the</p><p>nearest neighbor vacancy is shorter at the final state, making it</p><p>more difficult for the migrating vacancy to approach to the final</p><p>state, primarily due the electrostatic interaction of vacancies. For</p><p>models B and D, the number of Y cations at the second nearest</p><p>neighbor position to the migrating vacancy is larger at the final</p><p>state, making this final state energetically more favorable, and</p><p>thus the reduced migration barrier in the forward migration.</p><p>In model E, the forward jump barrier is higher than the backward</p><p>one. Considering the vacancy–cation interactions, the final</p><p>state in model E is energetically more favorable than the initial</p><p>state because of the presence of a larger number of 2NN Y</p><p>cations to the migrating vacancy at the final state. This suggests</p><p>the migration barrier for the forward jump to be lower than the</p><p>backward one. However, the result shows the opposite behavior,</p><p>because the migrating vacancy is actually getting closer to its</p><p>nearest neighbor vacancy at its final state. This indicates that the</p><p>vacancy–vacancy interaction dominates the vacancy–cation</p><p>interaction in determining the energy barriers in model E. All</p><p>implications related to the defect–defect interactions leading to</p><p>the resulting migration barriers in models A–E in the unstrained</p><p>state of YSZ are consistent with those suggested by previous</p><p>Fig. 6 NEB paths and the corresponding energies in the migration path</p><p>calculated for models A–E. The relative energy and the coordinates are</p><p>normalized with respect to the initial energy, and the distance between the</p><p>initial and the final state, respectively. The saddle point energy is the</p><p>highest energy state between the initial and the final states, and E!B</p><p>and E B stand for barriers in the forward and the backward jumps,</p><p>respectively.</p><p>4814 | J. Mater. Chem., 2010, 20, 4809–4819 This journal is ª The Royal Society of Chemistry 2010</p><p>reports by Bogicevic et al.36 Therefore, the results presented here</p><p>enable a reasonable model of YSZ for assessing the effect of</p><p>biaxial lattice strain on the ionic mobility.</p><p>IV. Migration barrier as a function of biaxial lattice strain</p><p>In assessing the role of the biaxial xy-plane lattice strain on the</p><p>oxygen vacancy transport, we investigated the corresponding</p><p>changes in the migration barriers, and the underlying reasons for</p><p>these changes on models A–E in Fig. 5. The xy-plane strains</p><p>considered here were 0.00, 0.02, 0.04, 0.06 and 0.08. The range of</p><p>strain considered in this study is far above the experimentally</p><p>measured fracture strain of a single-crystal YSZ that is less than</p><p>0.5%.58 However, large strain in brittle ceramics could be achieved</p><p>by creating coherent hetero-interfaces at the nanoscale.</p><p>A maximum 7% tensile lattice strain was suggested to exist in the</p><p>YSZ thin film interfaces coherently strained between STO layers</p><p>due to the mismatch between STO and YSZ.8 On the other hand,</p><p>ability of thin-film (10–30 nm) YSZ to sustain this large tensile</p><p>strain is unclear. We considered only tensile strain to increase the</p><p>interatomic distances in YSZ, allowing more space for the</p><p>oxygen migration. Results for model A are shown in Fig. 7.</p><p>The migration barrier monotonically reduces with the increase in</p><p>the biaxial lattice strain. The barrier in the forward jump</p><p>decreases from 0.47 eV at 3 ¼ 0.00 to 0.07 eV at 3 ¼ 0.04, suggesting</p><p>a significant enhancement of the vacancy diffusivity</p><p>through this path under strain. Beyond this strain level, the</p><p>vacancy was not stable at the initial site and favored to form at</p><p>the defined final site in model A, with no barrier in the suggested</p><p>transition path. This result implies that the defect interactions</p><p>influencing the stability and energetics can also change with the</p><p>strain state of YSZ. While not explicitly included in our DFT</p><p>calculations, we implicitly took into account the effect of the</p><p>strain on the stable defect configurations in the KMC simulations,</p><p>as explained in the following section. We considered both</p><p>the forward and the backward migration barriers, and the system</p><p>has a higher probability of taking the lower energy migration</p><p>paths favoring the more stable structures.</p><p>The picture is not as straight-forward in models B–E as in</p><p>model A that shows a monotonic reduction in the migration</p><p>barrier with strain. Depending on the extent of the lattice strain,</p><p>the consequent barrier either decreases or increases in the</p><p>configurations represented by models B–E, as shown in Fig. 8.</p><p>The lack of a straight-forward trend in the magnitude of</p><p>migration barrier as a function of lattice strain suggests that there</p><p>is more than one mechanism acting simultaneously upon strain</p><p>in YSZ. Here, we explain the competing effects of strain on the</p><p>vacancy migration barrier in YSZ from an atomistic and electronic</p><p>perspective. The oxygen vacancy migration process</p><p>involves bond-breaking and bond-making between the migrating</p><p>oxygen and the surrounding cations. Therefore, we investigated</p><p>the changes in the electronic density distribution of the valence</p><p>electrons describing the Zr, Y and O atoms near the migration</p><p>path to explain the local effects of strain in YSZ. We focused on</p><p>the two planes in the migration path as shown in Fig. 9. The</p><p>charge distribution on the ‘‘oxygen plane’’ in Fig. 9 correlates to</p><p>the ‘‘space available’’ for the oxygen, denoted as O, to migrate</p><p>along the prescribed path. The migration space is also related to</p><p>the cation-cation distance midst the migration path. On the other</p><p>hand, the ‘‘bonding plane’’ in Fig. 9 shows a measure of the</p><p>‘‘bonding strength’’ between the oxygen and the nearest cations</p><p>at the initial site, quantified in terms of the oxygen–cation ‘‘bond</p><p>thickness’’. It is this set of bonds, shown between the oxygen</p><p>denoted as O and the cation denoted as C, which are broken in</p><p>the initiation of oxygen diffusion. Both the migration space and</p><p>the bond strength change when YSZ is subjected to biaxial strain.</p><p>Fig. 10 illustrates an example in which the increase in strain</p><p>increases the open space, and also increases the O–C bonding</p><p>strength to the extent of local inelastic relaxation near the initial</p><p>state of oxygen. In this case, the lattice strain is beyond just</p><p>Fig. 7 (a) NEB paths and the corresponding energies in the migration</p><p>path calculated for model A. The relative energy and the migration</p><p>coordinates are normalized with respect to the initial energy, and the</p><p>distance between the initial and the final state, respectively. The saddle</p><p>point energy is the highest energy state between the initial and the final</p><p>states. (b) The vacancy migration barrier, EB, as a function of biaxial</p><p>lattice strain. The filled circles and the open circles indicate forward and</p><p>backward jump, respectively. The barrier decrease as the strain increases</p><p>in the migration path in model A.</p><p>Fig. 8 The vacancy migration barrier, EB, as a function of biaxial lattice</p><p>strain for models B–E (in a–d). The filled circles and the open circles</p><p>indicate the forward and the backward jump barriers, respectively.</p><p>This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4815</p><p>elastically extending the oxygen–cation bond length; the oxygen,</p><p>O, was actually attracted towards the cation C and broke its</p><p>bond from the opposite cation upon local relaxation of atoms. As</p><p>a result, a stronger binding of O to C arises, shown with the</p><p>thicker O–C bond in Fig. 10. These two consequences of strain</p><p>act oppositely on the migration barrier – the former (opening of</p><p>migration space) favors the migration and the latter (local</p><p>inelastic relaxations near the oxygen) inhibits the migration.</p><p>A specific illustration of these competing consequences of</p><p>biaxial strain in YSZ represented by model E is shown in Fig. 11.</p><p>The resulting change in the migration space was quantified as the</p><p>distance between the cations midst the migration path, D(Zr–Zr)</p><p>in this configuration. The bond O–C in Fig. 11(b) is one of the</p><p>oxygen–zirconium bonds that break upon the migration of</p><p>oxygen. The large reduction in the forward migration barrier</p><p>from 0.51 eV at 3 ¼ 0.00 to 0.10 eV at 3 ¼ 0.06 as shown in</p><p>Fig. 8(d) corresponds to the increase in the D(Zr–Zr) distance</p><p>(shown in Fig. 11(a)) and the weakening of the O–C bond (shown</p><p>in Fig. 11(b)). On the other hand, this trend is not monotonic. At</p><p>higher strain, 3 ¼ 0.08, the cation–cation distance decreases and</p><p>the oxygen–zirconium bond strengthens again, resulting in</p><p>a significant increase of the migration barrier from 0.10 eV at</p><p>3 ¼ 0.06 to 0.47 eV at 3 ¼ 0.08. At this strain, the local relaxation</p><p>of the atoms particularly re-strengthens the oxygen-zirconium</p><p>bond at the initial site of the oxygen, also affecting the migration</p><p>space. If assumed independent, the oxygen–cation bond strength</p><p>is more strongly correlated to this large increase in the barrier at</p><p>large strain seen in Fig. 8(d). The increase in the cation–oxygen</p><p>bond strength is evident at 3 ¼ 0.08 in Fig. 11(b),but a very small</p><p>change occurs in the D(Zr–Zr) in Fig. 11(a). While we present this</p><p>set of results explicitly for the model E, the same qualitative</p><p>correlation of the competing effects of the migration space and</p><p>the cation–oxygen bond strength governed the migration barrier</p><p>changes with strain consistently also in the other models A–D.</p><p>V. Effective oxygen diffusivity as a function of biaxial lattice</p><p>strain</p><p>Models A–E constituted the database of the migration paths and</p><p>barriers as a function of the biaxial strain in YSZ. Using this</p><p>database to govern the overall oxygen diffusion, the macroscopic</p><p>oxygen diffusivity in the bulk of YSZ was estimated by KMC</p><p>simulations as described in the Simulation Approach section.</p><p>A 15 15 15 simulation cell of 8% yttria doped zirconia was</p><p>used, consistent with the doping ratio in YSZ reported in the</p><p>Ref. 8 for the YSZ/STO heterolayers. Y cations and vacancies</p><p>were distributed randomly in the KMC simulation cell. At a first</p><p>sight, this configuration contradicts with the implications for the</p><p>presence of 1NN Y–Y pairs in YSZ as discussed in Introduction.</p><p>On the other hand, the rules in selecting the diffusion barriers in</p><p>the KMC simulations were taken consistently with the local</p><p>defect configurations, including the effect of Y–Y pairs as evaluated</p><p>in the DFT-NEB calculations in this work. Therefore, we</p><p>considered this model in KMC to be reasonable for assessing the</p><p>overall oxygen diffusivity in YSZ. Supporting this approach,</p><p>past DFT and KMC studies of oxygen diffusivity in YSZ with</p><p>randomly distributed Y atoms and vacancies were shown to give</p><p>Fig. 9 The electronic density distribution on the ‘‘oxygen plane’’ and the</p><p>‘‘bonding plane’’ at the initial and the saddle states. The migrating</p><p>oxygen, O, moves towards the open migration space, and the bond</p><p>between oxygen O and cation C breaks in the migration.</p><p>Fig. 10 Change in the electronic density distribution on the oxygen</p><p>plane and the bonding plane under biaxial lattice strain, 3, from 0.00 to</p><p>0.07. Gray and red spheres represent Zr and O atoms, respectively.</p><p>4816 | J. Mater. Chem., 2010, 20, 4809–4819 This journal is ª The Royal Society of Chemistry 2010</p><p>the correct yttria doping ratio for the maximum oxygen</p><p>diffusivity.25,50</p><p>The effect of the local defect environment around the migration</p><p>path on the barriers input to the KMC simulations was</p><p>taken into account as follows:</p><p>Definitions. E!a and E a are the forward and backward jump</p><p>energy barriers in each model a (a ¼ A, B, C, D, E),respectively.</p><p>ri</p><p>v and rf</p><p>v stand for the distance of the first nearest neighbor (1NN)</p><p>vacancy to the migrating vacancy at the initial and final state,</p><p>respectively. ni</p><p>Y, and nf</p><p>Y are the number of second nearest</p><p>neighbor (2NN) Y cations to the migrating vacancy at the initial</p><p>and final state, respectively.</p><p>Rules. The rules in determining the energy barrier for each hop</p><p>in the KMC simulations treat the local defect configurations in</p><p>a coarse grained fashion, based on the models A–E (Fig. 5)</p><p>assumed to govern the oxygen vacancy diffusion in YSZ.</p><p>In models A and C, the vacancy forward hopping (from initial to</p><p>final state) results in a decrease of the 1NN vacancy distance to</p><p>the migrating vacancy, that is ri</p><p>v > rf</p><p>v. Then, in the KMC simulation,</p><p>E!a and E a (a ¼ A, C) are chosen for the cases of ri</p><p>v > rf</p><p>v,</p><p>and ri</p><p>v<rf</p><p>v, respectively. In models B and D, the vacancy forward</p><p>hopping (from initial to final state) results in an increase in the</p><p>number of 2NN Y to the migrating vacancy, that is ni</p><p>Y < nf</p><p>Y.</p><p>Then, in the KMC simulation, E!a and E a (a ¼ B, D) are chosen</p><p>for the cases of ni</p><p>Y <nf</p><p>Y and ni</p><p>Y > nf</p><p>Y, respectively. The explicit</p><p>definition of the rules implemented in the KMC simulations is as</p><p>follows:</p><p>(1) For the vacancy migration through a Zr–Zr pair midst the</p><p>migration path, we select the barrier among E!A, E A E!B and E B</p><p>according to the following criteria:</p><p> If ri</p><p>v > rf</p><p>v, select E!A; and if ri</p><p>v < rf</p><p>v, select E A.</p><p> If ri</p><p>v ¼ rf</p><p>v and ni</p><p>Y > nf</p><p>Y, select E!B; or if ri</p><p>v ¼ rf</p><p>v and ni</p><p>Y < nf</p><p>Y,</p><p>select E B.</p><p> If ri</p><p>v ¼ rf</p><p>v and ni</p><p>Y ¼ nf</p><p>Y, randomly choose barriers from E!A,</p><p>E A, E!B, and E B.</p><p>(2) For the vacancy migration through a Zr–Y pair midst the</p><p>migration path, we select the barrier among E!C, E C, E!D and E D</p><p>according to the following criteria:</p><p> If ri</p><p>v > rf</p><p>v, select E!C; if ri</p><p>v < rf</p><p>v, select E C.</p><p> If ri</p><p>v ¼ rf</p><p>v and ni</p><p>Y > nf</p><p>Y, select E!D; or if ri</p><p>v ¼ rf</p><p>v and ni</p><p>Y < nf</p><p>Y,</p><p>select E D.</p><p> If ri</p><p>v ¼ rf</p><p>v and ni</p><p>Y ¼ nf</p><p>Y, randomly choose barriers from E!C,</p><p>E C, E!D and E D.</p><p>(3) Vacancy migration through Y–Y pair is prohibited.</p><p>When ri</p><p>v ¼ rf</p><p>v and ni</p><p>Y ¼ nf</p><p>Y, the migration path is surrounded by</p><p>a symmetric distribution of defects. This case was not possible to</p><p>account for in the DFT-NEB calculations due to the computational</p><p>limitations on the simulation size. We assumed such</p><p>a symmetric path to be governed by a migration barrier that is</p><p>the average of the barriers from the non-symmetric models represented</p><p>in A–E here. Therefore, the random selections among</p><p>the entire set of E!a and E a (for a ¼ A, B, C, D, E) during the 108</p><p>KMC steps should provide sufficient statistics to represent the</p><p>average migration barrier for the symmetric path, ri</p><p>v ¼ rf</p><p>v and</p><p>ni</p><p>Y ¼ nf</p><p>Y. Moreover, the probability of having symmetric path in</p><p>the KMC runs was less than 3%. Therefore, the error that may</p><p>arise due to selecting an average migration barrier for the</p><p>symmetric configuration is expected to have little significance on</p><p>the effective diffusivity.</p><p>As seen from the rules above, the energy barrier for each hop</p><p>in the KMC simulations was chosen based on a coarse grained</p><p>treatment of the local defect configurations around the migrating</p><p>vacancy. We believe this approach gives a reasonable result,</p><p>particularly based on the validity of the models A–E (Fig. 5) to</p><p>govern the oxygen vacancy diffusion in YSZ. On the other hand,</p><p>the reader can refer to a more accurate treatment of the diffusion</p><p>barriers in KMC simulations using the cluster expansion</p><p>method,43,59 which can account for the energy deviations due to</p><p>the local atomic structures around each hop. In that case, the</p><p>KMC model itself has to take into account the stable distributions</p><p>of vacancies and cations (as in subsections I and II),rather</p><p>than a random distribution, making the task exhaustive.</p><p>The KMC simulations were performed as a function of</p><p>temperature ranging as 400, 600, 800 and 1000 K, and of biaxial</p><p>strain ranging as 0.00, 0.02, 0.04, 0.06 and 0.08. Ten simulations</p><p>with different distributions of dopant cations and vacancies were</p><p>conducted for each temperature and strain state. The resulting</p><p>strain dependence of the effective oxygen tracer diffusivity DO is</p><p>shown in Fig. 12. y-axis in the figure is on a relative scale with</p><p>Fig. 11 (a) Change in the cation–cation distance, D(Zr–Zr),and (b) the</p><p>bond thickness (strength) for O–C in the migration under biaxial lattice</p><p>strain from 0.00 to 0.08 for the forward jump in model E (see also</p><p>Fig. 8(d)). Increase in the migration space D(Zr–Zr),and the weakening</p><p>of the O–C bond correspond to the decrease of the migration barrier in</p><p>Fig. 8(d),and vice versa. Large gray spheres and small red spheres are Zr</p><p>and O atoms, respectively. The isoline interval is 0.5 A3.</p><p>This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4817</p><p>respect to the oxygen diffusivity, D0</p><p>O at 3 ¼ 0.00 for the lowest</p><p>temperature simulated, 400 K. The reported results show the</p><p>average values and the error bars over the ten simulations. The</p><p>result is not a simple monotonic increase of diffusivity with</p><p>increasing strain in the entire range considered. The oxygen</p><p>diffusivity exhibits an exponential increase up to a critical value,</p><p>4% of tensile strain. This increase is more significant at the lower</p><p>temperatures. The maximum enhancement in diffusivity</p><p>compared to the unstrained state in YSZ is 6.8 103 times higher</p><p>at 4% strain and at 400 K. At the higher strain states, the</p><p>diffusivity decreases. This is attributed to the local relaxations</p><p>beyond a limit of elastic bond strain, resulting in the strengthening</p><p>of the local oxygen–cation bonds that increases the</p><p>migration barrier.</p><p>Our first-principles-based simulation results indicate that</p><p>a biaxial strain in YSZ induced by a lattice-mismatch coherent</p><p>interface can significantly increase the diffusivity up to a critical</p><p>strain state, and worsen it beyond that strain state. One could</p><p>view this critical strain state as the optimal strain to induce in the</p><p>material to attain the maximum benefit of the biaxial lattice</p><p>strain on the oxygen diffusivity. We note that a recent prediction</p><p>by Schichtel et al.22 using a model based on the elastic properties</p><p>of YSZ estimated 2.5 orders of magnitude increase in the ionic</p><p>conductivity of YSZ at 7% strain at 573 K. The model of</p><p>Schichtel et al. relies on the isotropic pressure induced by the</p><p>strain to enhance the diffusivity, and does not take into account</p><p>the competing relaxations at these large tensile strain states in</p><p>YSZ. Therefore, it cannot capture the local bond-strength</p><p>changes that may mitigate the diffusivity at large strain values.</p><p>In spite of the differences in our and Schichtel’s model, the</p><p>magnitude of the maximum increase in the oxygen diffusivity in</p><p>strained YSZ predicted by these models range from 2.5 to 3.8</p><p>orders of magnitude. Therefore, the lattice strain (elastic and</p><p>inelastic) alone cannot explain the eight orders of magnitude</p><p>increase of ionic conductivity in YSZ reported by Barriocanal</p><p>et al.8 These results do not take into account the space-charge</p><p>effects in altering the anion concentrations near the YSZ interface.</p><p>On the other hand, the contribution to the increase in the</p><p>anion diffusivity due to a possible space-charge-induced increase</p><p>in the vacancy concentration is almost linear at low concentration</p><p>of vacancies (see eqn (3)) as in the 9% Y2O3 doped YSZ.</p><p>Therefore, this contribution is expected to be small compared to</p><p>the exponential increase found here due to the applied biaxial</p><p>lattice strain.60</p><p>Conclusion</p><p>This work provides a mechanistic and systematic assessment of</p><p>the effect of biaxial lattice strain on the oxygen diffusivity in</p><p>YSZ, independently from the defect chemistry near a heterointerface</p><p>with a dissimilar oxide, e.g. YSZ/STO as reported by</p><p>Barriocanal et al.8 Using DFT/NEB simulations, we uncovered</p><p>the underlying microscopic mechanism by which the biaxial</p><p>tensile strain acts on the oxygen vacancy transport in YSZ. The</p><p>migration barrier correlates fundamentally with two competing</p><p>processes acting in parallel in the presence of lattice strain: (1) the</p><p>migration space, quantified as the cation–cation distance midst</p><p>the migration path for the oxygen, and (2) the strength of the</p><p>bond between the oxygen and the nearest cation preventing the</p><p>oxygen from migrating. Increasing of the migration space and</p><p>weakening of the oxygen–cation bond contribute to the decrease</p><p>of the vacancy migration barrier, and vice versa. In assessing the</p><p>effective diffusivity as a function of strain in YSZ, we constructed</p><p>a data base of migration paths and barriers to account for the</p><p>important effects of local distributions of vacancies and dopant</p><p>cations near the migration path. KMC simulations using this</p><p>database showed that the macroscopic oxygen diffusivity</p><p>increases exponentially in YSZ up to a critical value of biaxial</p><p>tensile strain, with a more significant impact at the lower</p><p>temperatures. Beyond this critical value, the strain state in YSZ</p><p>worsens the diffusivity due to the local relaxations that</p><p>strengthen the oxygen–cation bond at the initiation of migration.</p><p>One could view this critical strain state as the optimal strain, or</p><p>the fastest strain, to attain the highest benefit of the biaxial lattice</p><p>strain on the oxygen diffusivity in the material system. The</p><p>maximum enhancement in diffusivity in YSZ was 6.8 103 times</p><p>higher at 4% strain and temperature of 400 K. This is a significant</p><p>and desirable increase. Nevertheless, it is far below the eight</p><p>orders of magnitude increase reported by Barriocanal et al. for</p><p>the thin-film YSZ layers coherently strained between STO</p><p>layers.8 Therefore, we believe that the lattice strain alone cannot</p><p>be responsible for the colossally high ionic conductivity suggested</p><p>to occur in YSZ in ref. 8. Nevertheless, given our results,</p><p>controlling the lattice strain by direct mechanical load or by</p><p>creating a coherent hetero-interface to induce an optimal strain</p><p>state can enable desirably high ionic conductivity at reduced</p><p>temperatures in the fluorite material, the YSZ in this case.</p><p>The insights gained here for the role of lattice strain on the</p><p>local bonding structure and ionic charge transport process are of</p><p>importance for modulating the transport properties in a variety</p><p>of solid state conducting material applications, beyond the</p><p>specific SOFC electrolyte material studied in this work. Interesting</p><p>examples may extend to Li-ion batteries, oxygen gas</p><p>sensors, semi-conductors and photovoltaics. We believe that the</p><p>mechanistic results obtained in this study, particularly on the</p><p>Fig. 12 The relative enhancement of the oxygen diffusivity calculated by</p><p>the kinetic Monte Carlo simulations at different temperatures and strain</p><p>states. D0</p><p>O is the oxygen diffusivity at 3 ¼ 0.00 at 400 K. The diffusivity</p><p>increases exponentially with strain, and shows a maximum at (or near)</p><p>3 ¼ 0.04 followed by a decrease at higher strains. The increase in diffusivity</p><p>is more significant at lower temperatures.</p><p>4818 | J. Mater. Chem., 2010, 20, 4809–4819 This journal is ª The Royal Society of Chemistry 2010</p><p>non-linear and competing consequences of strain on the bonding</p><p>structure in the vicinity of the ionic carrier species, provides</p><p>insights for creating functional strained nano-structures for ionic</p><p>and electronic conducting devices, including but not limited to</p><p>SOFCs.</p><p>Acknowledgements</p><p>We thank the Nuclear Regulatory Commission for financial</p><p>support, and the National Science Foundation for computational</p><p>support through the TeraGrid Advanced Support</p><p>Program under Grant No. TG-ASC090058.</p><p>References</p><p>1 J. C. Boivin and G. Mairesse, Chem. Mater., 1998, 10, 2870.</p><p>2 T. Suzuki, I. Kosacki and H. U. Anderson, Solid State Ionics, 2002,</p><p>151, 111.</p><p>3 I. Kosacki, C. M. Rouleau, P. J. Becher, J. Bentley and</p><p>D. H. Lowndes, Solid State Ionics, 2005, 176, 1319.</p><p>4 J. H. Shim, C.-C. Chao, H. Huang and F. B. Printz, Chem. Mater.,</p><p>2007, 19, 3850.</p><p>5 J. Maier, Nat. Mater., 2005, 4, 805.</p><p>6 N. Sata, K. Eberman, K. Eberl and J. Maier, Nature, 2000, 408, 946.</p><p>7 X. X. Guo, I. Matei, J.-S. Lee and J. Maier, Appl. Phys. Lett., 2007,</p><p>91, 103102.</p><p>8 J. G. Barriocanal, A. R. Calzada, M. Varela, Z. Sefrioui, E. Iborra,</p><p>C. Leon, S. J. Pennycook and J. Santamaria, Science, 2008, 321, 676.</p><p>9 S. Ramanathan, J. Vac. Sci. Technol., A, 2009, 27, 1126.</p><p>10 A. J. Jacobson, Chem. Mater., 2010, 22, 660.</p><p>11 X. Guo, Science, 2009, 324, 465; J. G. Barriocanal, A. R. Calzada,</p><p>M. Varela, Z. Sefrioui, E. Iborra, C. Leon, S. J. Pennycook and</p><p>J. Santamaria, Science, 2009, 324, 465.</p><p>12 M. P. O’Callaghan, A. S. Powell, J. J. Titman, G. Z. Chen and</p><p>E. J. Cussen, Chem. Mater., 2008, 20, 2360.</p><p>13 Y. Hoshina, K. Iwasaki, A. Yamada and M. Konagai, Jpn. J. Appl.</p><p>Phys., 2009, 48, 04C125.</p><p>14 J. Munguıa, G. Bremond, J. M. Bluet, J. M. Hartmann and</p><p>M. Mermoux, Appl. Phys. Lett., 2008, 93, 102101.</p><p>15 V. R. D’Costa, Y.-Y. Fang, J. Tolle, J. Kouvetakis and J. Menendez,</p><p>Phys. Rev. Lett., 2009, 102, 107403.</p><p>16 M. Mavrikakis, P. Stoltze and J. K. Nørskov, Catal. Lett., 2000, 64,</p><p>101.</p><p>17 J. Kilner, Nat. Mater., 2008, 7, 838.</p><p>18 J. Maier, Prog. Solid State Chem., 1995, 23, 171.</p><p>19 A. K. Ivanov-****z, Crystallogr. Rep., 2007, 52, 129.</p><p>20 D. C. Sayle, J. A. Doig, S. C. Parker and G. W. Watson, Chem.</p><p>Commun., 2003, 1804.</p><p>21 C. Korte, A. Peters, J. Janek, D. Hesse and N. Zakharov, Phys. Chem.</p><p>Chem. Phys., 2008, 10, 4623.</p><p>22 N. Schichtel, C. Korte, D. Hesse and J. Janek, Phys. Chem. Chem.</p><p>Phys., 2009, 11, 3043.</p><p>23 K. Suzuki, M. Kubo, Y. Oumi, R. Miura, H. Takaba, A. Fahmi,</p><p>A. Chatterjee, K. Teraishi and A. Miyamoto, Appl. Phys. Lett.,</p><p>1998, 73, 1502.</p><p>24 W. Araki, Y. Imai and T. Adachi, J. Eur. Ceram. Soc., 2009, 29, 2275.</p><p>25 R. Pornprasertsuk, P. Ramanarayanan, C. B. Musgrave and</p><p>F. B. J. Prinz, J. Appl. Phys., 2005, 98, 103513.</p><p>26 W. L. Roth, R. Wong, A. I. Goldman, E. Canova, Y. H. Kao and</p><p>B. Dunn, Solid State Ionics, 1986, 18–19, 1115; H. Morikawa,</p><p>Y. Shimizugawa, F. Marumo, T. Harasawa, H. Ikawa, K. Tohji</p><p>and Y. Udagawa, J. Ceram. Soc. Jpn, 1988, 96, 253; M. H. Tuilier,</p><p>J. Dexpert-Ghys, H. Dexpert and P. Lagarde, J. Solid State Chem.,</p><p>1987, 69, 153.</p><p>27 M. Weller, Z. Metallkd., 1993, 84, 6.</p><p>28 D. Steele and B. E. F. Fender, J. Phys. C: Solid State Phys., 1974, 7, 1.</p><p>29 B. W. Veal, A. G. McKale, A. P. Paulikas, S. J. Rothman and</p><p>L. J. Nowicki, Physica B+C, 1988, 150, 234; D. Komyoji,</p><p>A. Yoshiasa, T. Moriga, S. Emura, F. Kanamaru and K. Koto,</p><p>Solid State Ionics, 1992, 50, 291.</p><p>30 C. R. A. Catlow, A. V. Chadwick, G. N. Greaves and L. M. Moroney,</p><p>J. Am. Ceram. Soc., 1986, 69, 272.</p><p>31 P. Li and J. E. Penner-Hahn, Phys. Rev. B: Condens. Matter, 1993, 48,</p><p>10063; P. Li and J. E. Penner-Hahn, Phys. Rev. B: Condens. Matter,</p><p>1993, 48, 10074; P. Li and J. E. Penner-Hahn, Phys. Rev. B:</p><p>Condens. Matter, 1993, 48, 10082.</p><p>32 J. P. Goff, W. Hayes, S. Hull, M. T. Hutchings and K. N. Clausen,</p><p>Phys. Rev. B, 1999, 59, 14202.</p><p>33 K. Kawata, H. Maekawa, T. Nemoto and T. Yamamura, Solid State</p><p>Ionics, 2006, 177, 1687.</p><p>34 G. Stapper, M. Bernasconi, N. Nicoloso and M. Parrinello, Phys.</p><p>Rev. B: Condens. Matter Mater. Phys., 1999, 59, 797.</p><p>35 A. Bogicevic, C. Wolverton, G. M. Crosbie and E. B. Stechel, Phys.</p><p>Rev. B: Condens. Matter Mater. Phys., 2001, 64, 014106.</p><p>36 A. Bogicevic and C. Wolverton, Phys. Rev. B: Condens. Matter</p><p>Mater. Phys., 2003, 67, 024106.</p><p>37 A. Bogicevic and C. Wolverton, Europhys. Lett., 2001, 56, 393.</p><p>38 X. Xia, R. Oldman and R. Catlow, Chem. Mater., 2009, 21, 3576.</p><p>39 M. Sakib Khan, M. Saiful Islam and D. R. Bates, J. Mater. Chem.,</p><p>1998, 8, 2299.</p><p>40 M. O. Zacate, L. Minervini, D. J. Bradfield, R. W. Grimes and</p><p>K. E. Sickafus, Solid State Ionics, 2000, 128, 243.</p><p>41 D. A. Andersson, S. I. Simak, N. V. Skorodumova, I. A. Abrikosov</p><p>and B. Johansson, Proc. Natl. Acad. Sci. U. S. A., 2006, 103, 3518.</p><p>42 S. Garcıa-Martın, D. P. Fagg and J. T. S. Irvine, Chem. Mater., 2008,</p><p>20, 5933.</p><p>43 A. Predith, G. Ceder, C. Wolverton, K. Persson and T. Mueller, Phys.</p><p>Rev. B: Condens. Matter Mater. Phys., 2008, 77, 144104.</p><p>44 F. Pietrucci, M. Bernasconi, A. Laio and M. Parrinello, Phys. Rev. B:</p><p>Condens. Matter Mater. Phys., 2008, 78, 094301.</p><p>45 G. Kresse and J. Hafner, Phys. Rev.B: Condens. Matter, 1993, 47, 558.</p><p>46 G. Kresse and J. Furthm€uller, Phys. Rev. B: Condens. Matter, 1996,</p><p>54, 11169.</p><p>47 P. E. Bl€ochl, Phys. Rev. B: Condens. Matter, 1994, 50, 17953.</p><p>48 G. Kresse and D. Joubert, Phys. Rev. B: Condens. Matter Mater.</p><p>Phys., 1999, 59, 1758.</p><p>49 J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett., 1996, 77,</p><p>3865.</p><p>50 R. Krishnamurthy, Y.-G. Yoon, D. J. Srolovitz and R. Car, J. Am.</p><p>Ceram. Soc., 2004, 87, 1821.</p><p>51 G. Mills and H. Jonsson, Phys. Rev. Lett., 1994, 72, 1124.</p><p>52 G. Henkelman, B. P. Uberuaga and H. Jonsson, J. Chem. Phys., 2000,</p><p>113, 9901.</p><p>53 A. Van der Ven, G. Ceder, M. Asta and P. D. Tepesch, Phys. Rev. B:</p><p>Condens. Matter Mater. Phys., 2001, 64, 184307.</p><p>54 M. Martin, J. Electroceram., 2006, 17, 765.</p><p>55 F. Frey, H. Boysen and I. K- Bischoff, Z. Kristallogr., 2005, 220, 1017.</p><p>56 D. N. Argyriou, M. M. Elcombe and A. C. Larson, J. Phys. Chem.</p><p>Solids, 1996, 57, 183.</p><p>57 N. Ishizawa, Y. Matsushima, M. Hayashi and M. Ueki, Acta</p><p>Crystallogr., Sect. B: Struct. Sci., 1999, 55, 726.</p><p>58 D. Baither, M. Bartsch, B. Baufeld, A. Tikhonovsky, M. R€uhle and</p><p>U. Messerschmidta, J. Am. Ceram. Soc., 2001, 84, 1755.</p><p>59 A. Van der Ven, J. C. Thomas, Q. C. Xu, B. Swoboda and</p><p>D. Morgan, Phys. Rev. B: Condens. Matter Mater. Phys., 2008, 78,</p><p>104306.</p><p>60 A. Kushima and B. Yildiz, ECS Trans., 2009, 25, 1599.</p><p>This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4819</p><p></p><p>heres the first one..more to come..</p></blockquote><p></p>
[QUOTE="disturbed, post: 81213, member: 4763"] This ones going to be difficult to explain. I have some more information to regurgitate(love that).. I am now in contact with chemists at the local college and talking to Ivoclar to get the specifics on the molecular exchange between Zirliner and Zirconium, as I have proof that IF there is one it is VERY weak. more likely there is not. FROM THE JOURNAL OF MATERIALS CHEMISTRY accepted march 2010 sooo... recent this is complicated info..start with Section IV (fig 4) Migration barrier as a function of biaxial lattice strain. than read conclusion. this is above a lot of peoples understanding so I will be producing studies from colleges across the globe and doing my best to explain the chemistry over the next few months... and yes.. its another attempt to show you all that a bond other than mechanical retention is IMPOSSIBLE to zirconium and that their will be more chipping and failures to come. if you only care about your restorations lasting a max of 5 years than go ahead and ignore this.. it will probably only confuse and uspet you.. If I find that materials are added to zirconium oxide that allow for a chemical bond that gives an adequate bond strength I will post.. but as the vitrification of the material happens around 2700 degrees C I doubt it is possible. [ATTACH]3849.vB[/ATTACH] Oxygen ion diffusivity in strained yttria stabilized zirconia: where is the fastest strain? Akihiro Kushima and Bilge Yildiz* Received 8th January 2010, Accepted 9th March 2010 First published as an Advance Article on the web 29th April 2010 DOI: 10.1039/c000259c We present the mechanism and the extent of increase in the oxygen anion diffusivity in Y2O3 stabilized ZrO2 (YSZ) under biaxial lattice strain. The oxygen vacancy migration paths and barriers in YSZ as a function of lattice strain was assessed computationally using density functional theory (DFT) and nudged elastic band (NEB) method. Two competing and non-linear processes acting in parallel were identified to alter the migration barrier upon applied strain: (1) the change in the space, or electronic density, along the migration path, and (2) the change in the strength of the interatomic bonds between the migrating oxygen and the nearest neighbor cations that keep the oxygen from migrating. The increase of the migration space and the weakening of the local oxygen–cation bonds correspond to a decrease of the migration barrier, and vice versa. The contribution of the bond strength to the changes in the migration barrier is more significant than that of the opening of migration space in strained YSZ. A database of migration barrier energies as a function of lattice strain for a set of representative defect distributions in the vicinity of the migration path inYSZwas constructed. This database was used in kinetic Monte Carlo (KMC) simulations to estimate the effective oxygen diffusivity in strained YSZ. The oxygen diffusivity exhibits an exponential increase up to a critical value of tensile strain, or the fastest strain. This increase is more significant at the lower temperatures. At the strain states higher than the critical strain, the diffusivity decreases. This is attributed to the local relaxations at large strain states beyond a limit of elastic bond strain, resulting in the strengthening of the local oxygen–cation bonds that increases the migration barrier. The highest enhancement of diffusivity in 9%-YSZ compared to its unstrained state is 6.8 103 times at 4% strain and at 400 K. The results indicate that inducing an optimal strain state by direct mechanical load or by creating a coherent hetero-interface with lattice mismatch can enable desirably high ionic conductivity in YSZ at reduced temperatures. The insights gained here particularly on the nonlinear and competing consequences of lattice strain on the local bonding structure and charge transport process are of importance for tuning the ionic transport properties in a variety of solid-state conducting material applications, including but not limited to fuel cells. Introduction There is much interest in decreasing the Solid Oxide Fuel Cell (SOFC) operating temperature to an intermediate-to-low temperature range to improve chemical and mechanical stability and reduce the cost. Enhancing the ionic transport in the electrolyte and cathode materials is a key for this purpose. Over the past several decades, developing new compositions with high ionic mobility at low temperatures has been the primary activity.1 In recent years, however, nanoscale structures were shown as means for potentially increasing the transport properties in oxides of relevance to SOFCs.2–9 The importance of utilizing the favorable interface properties of nanoscale structures shifted the paradigm from developing new material compositions to tailoring nano-structures using existing materials for the above purpose.5,10 Among these, nanoscale hetero-layered oxides were suggested to enhance the ionic conductivity along the interfaces by orders of magnitude.5–9 The presence of an interface between dissimilar oxides is commonly the source of charged mobile defects and local anisotropic distortions and strained interatomic bonds. These contribute to changes in the local electronic structure; formation, coordination and redistribution of charged defects; and mobility of defects in the vicinity of interfaces. Most recently, an eight orders of magnitude increase in ionic conductance was reported by Barriocanal et al.8 for 1–30 nm-thick yttria stabilized zirconia (YSZ) layers coherently ‘‘strained’’ between dielectric SrTiO3 (STO) layers. The exact nature of the ionic vs. electronic conductance induced in the vicinity of the YSZ/STO interface remains debatable.11 While these observations imply that the hetero-interfaces with controlled lattice-strain and defect-chemistry can play an important role in improving the ionic mobility, results to date remain phenomenological without systematic theoretical efforts. Beyond the SOFC materials domain, similar questions for the role of hetero-interfaces on the local bonding structure and the charge transport process are open for tuning the transport properties in a variety of solid-state conducting material applications.12–15 Systematic experimental and theoretical investigation of the role of local strain states on tuning the mechanism and magnitude of charge transfer properties in these materials is prone to pursuit. Department of Nuclear Science and Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA, 02139, USA. E-mail: [email]byildiz@mit.edu[/email]; Tel: +617-324-4009 This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4809 PAPER [url=http://www.rsc.org/materials]Journal of Materials Chemistry Home-High impact applications, properties and synthesis of exciting new materials[/url] | Journal of Materials Chemistry Two main mechanisms contribute to increase the ionic transport at the hetero-interface: (1) a favorable strain state at the interface to shift and/or change the symmetry of electron energy levels to provide for improved charge transfer16 and ion mobility,17 and (2) the alteration of the defect chemistry7,18 near the interface to enhance the density and distribution of desired charge carriers. To understand the underlying mechanism of oxygen anion transport in anion conducting ceramics, a significant number of theoretical studies using atomistic simulations based on both density functional theory (DFT) and empirical potentials have been conducted (examples related to this work in ref. 19–25). However, most of them are focused on the ‘‘bulk’’ or the ‘‘surface’’ ionic transport. Only few theoretical studies exist for the anion transport at hetero-interfaces, on specific cases, i.e. CaF2/BaF2, CeO2/YSZ,20 with no generalization of results. Of particular interest for this paper is the effect of lattice strain on the ionic conductivity in the YSZ, as may be encountered due to a lattice mismatch near the YSZ/STO interface. In probing the effect of strain on conductivity in YSZ, a number of experimental and simulation studies have been conducted in the past.21–24 Even though the past results imply a potential increase in the conductivity, the magnitude of increase was far lower than what was reported recently in ref. 8 for the YSZ/STO interface. Furthermore, the reported increases in magnitude of the enhancement are inconsistent among these past studies because of the differences in conditions, structures and methods that were considered. For example, results were based on examining either incoherent hetero-interfaces which cannot sustain large lattice strain,21 or low strain conditions,24 or molecular dynamics simulation with empirical potentials23,24 which are not appropriate for interfaces, or phenomenological modeling using macroscopic elastic properties of the material.22 In essence, the mechanism by which a biaxial lattice strain alters the anion transport in the YSZ fluorite structure was not uncovered. Given these past inconsistencies and open questions, we adopt a firstprinciples- based approach in probing the mechanistic picture of the ionic conduction in strained YSZ. Our specific objectives are to elucidate the mechanism and to identify the possible extent of increase in the oxygen ion conductivity in YSZ due to a biaxial lattice strain. In order to attain reliable results in probing the effect of strain on ionic transport in YSZ, it is important to consider a realistic representation of the defect structure, in particular the distribution of vacancies and dopant cations in YSZ. This is because the oxygen vacancy migration barrier, which can be altered by strain, depends strongly on the local defect distributions.25 On the other hand, it is difficult to characterize experimentally the vacancy and cation distribution in YSZ at atomic resolution, therefore this type of information is not available from literature straightforwardly. In spite of the difficulties, past work on this topic provided useful insights on the vacancy–vacancy, vacancy– cation and cation–cation interactions, enabling a reasonable description of the defect structure in YSZ. Here we provide a summary of defect–defect interactions that we considered for constructing the YSZ model. The final YSZ model used in our work should reasonably capture these reported defect–defect interactions. First of these interactions considered, namely the vacancy–cation interaction, determines the relative position of the vacancy with respect to the Zr and Y in YSZ. In earlier studies for this interaction, there were contradicting experimental results; whether the vacancy is first nearest neighbor26–28 (1NN) or second nearest neighbor29–31 (2NN) to the Y cations. However, later experimental30–33 and theoretical34–40 studies indicated that the vacancy favors to bind to the host Zr cations, resulting in a configuration such that the vacancy is 1NN to the Zr and 2NN to the Y cations. Using DFT calculations, Bogicevic and Wolverton37 showed that the local relaxation of atoms was responsible for the vacancy to favor the 2NN position to the Y cation; opposite to what is expected from the electrostatic interactions alone. The same trend was reported in doped CeO2 when the dopant radius is larger than the host cation,41 as in YSZ. For the vacancy–vacancy and the cation– cation interactions, they concluded that the vacancies locate themselves as far away as possible from each other, and that the 1NN Y–Y pair configurations are slightly more favorable than the 2NN Y–Y pairs.36,38 More recent experimental work using selected area diffraction also indicated a random distribution of oxygen vacancies, maximizing the distance between them, in up to 10% yttria doped YSZ.42 Consistent with the findings of Bogicevic and Wolverton, theoretical results by Predith et al.43 suggest even the existence of the ordering of the dopant cations (Y) in YSZ at high doping concentrations, and this cation ordering could exist locally in YSZ at the lower Y concentration range. Among these three types of defect–defect interactions, the vacancy–vacancy interaction is the strongest followed by the vacancy–cation interaction, and the cation–cation interaction is the weakest of the three.36 Supporting the importance of the vacancy–vacancy interactions, it was shown recently that the vacancy distribution around the migration path could significantly impact the value of the migration barriers in YSZ.44 In this study, we considered only the strain state in the bulk to provide for improved ionic mobility in YSZ and did not include an interface explicitly as in Barriocanal’s work.8 Therefore, the model in this paper excludes the effect of altered defect chemistry near the interfaces of YSZ as a source of increase in oxygen diffusivity. In the following sections, we describe the simulation approach to calculate the oxygen vacancy migration path and energy barriers. We report the energetically favorable distributions of oxygen vacancies and Y cations in the YSZ model structure to enable a realistic representation of YSZ in this study. The mechanism by which the lattice strain increases or decreases the anion transport in YSZ are discussed. Because the vacancy migration barrier depends on the local defect configurations, a database of migration barrier energies as a function of lattice strain for a set of representative defect distributions in the vicinity of the migration path are constructed and presented. Finally, the effective tracer diffusivity of oxygen in YSZ at different strain states is presented and interpreted in comparison to the recent results reported for the YSZ/STO interface.8 Simulation approach The YSZ model was established to capture the charged defect interactions and the defect distribution in the bulk YSZ, as summarized in the Introduction. Upon identifying the stable vacancy–cation distribution, the effect of the biaxial lattice strain in the xy-plane on the migration barrier and on the 4810 | J. Mater. Chem., 2010, 20, 4809–4819 This journal is ª The Royal Society of Chemistry 2010 effective vacancy diffusivity were calculated and interpreted. We used first-principles-based simulations to directly and accurately probe the atomistic-scale nature of the vacancy migration in the presence of lattice strain. The kinetic Monte Carlo (KMC) method was used for calculating the effective diffusivities. The first principles calculations for identifying the energies were performed by the Vienna ab initio simulation package (VASP),45,46 employing density functional theory (DFT) using a plane-wave basis set. Projector-augmented wave (PAW) method47,48 with plane waves up to the energy cutoff at 400 eV was used and the exchange–correlation energy was evaluated by the generalized gradient approximation (GGA) using the Perdew–Burke–Ernzerhof (PBE) function.49 Y, Zr and O atoms were described by 11 (4s24p65s24d1),12 (4s24p65s24d2),and 6 (2s22p4) valence electrons, respectively. The G point in the Brillouin zone was selected in the calculation. We confirmed the energy convergence by comparing the energy obtained from each of the k-point meshes of 3 3 3, 5 5 5 and 7 7 7. The calculated energies are within 0.6% error compared to the energy obtained using only the G point. The same condition was applied to all of the simulations in this study unless stated otherwise. YSZ has the fluorite crystal structure with the cations (Y3+, Zr4+) occupying the FCC lattice sites and the O2 ions occupying the tetragonal sites. For every two Y atoms, one oxygen vacancy is created to satisfy charge neutrality. The oxygen vacancy migrates by hopping to the adjacent tetragonal site as schematically illustrated in Fig. 1. The migration barrier depends on the neighboring atoms 1–6 in the figure.25 The surrounding cations within a larger cutoff radius to the migration path must also be accounted for in order to accurately quantify the effect of the strain on the migration barrier.50 However, this leads to an impractical number of patterns (well above 10 000 configurations just by including the next nearest neighbor cations to the migration path) to be handled by the DFT simulation. Therefore, we selected a subset of defect structure patterns near the migration path, forming a database of migration paths and barriers. This database serves as input for the KMC calculations in estimating the effective diffusivity as a function of lattice strain and temperature. In constructing this database, first, we considered only two combinations of atoms at the sites 3 and 4, as Zr–Zr and Zr–Y, in the migration path shown in Fig. 1. The migration barrier through the Y–Y pair was found to be significantly higher than these two combinations, and this makes the Y–Y pair unlikely on a vacancy migration path. This selection was also validated by the work of Car et al.,50 which implied the Zr–Zr and Zr–Y subset at sites 3 and 4 to be sufficient to capture the quantitative trends in ionic conductivity in YSZ as a function of Y2O3 doping concentration. Second, we considered the defect configurations outside of the first neighbor cations to the migration path to affect the barrier. The particularly important impact of the vacancy interactions on the migration barrier was recently reported by Pietrucci et al.44 Therefore, in creating the migration path and barrier database, we varied the relative positions of the Y cations and the vacancies with respect to the initial and final positions of the vacancy on the migration path. The details of the migration path configurations in this database are provided in the following section. For each migration path configuration, nudged elastic band (NEB) method51 was used to find the vacancy migration barrier as a function of biaxial strain on the xy-plane. Given that the vacancy migration path in YSZ is rather simple with single hops between the adjacent tetragonal sites, the saddle point is expected near the midst of the migration path.25,50 Three images between the initial and the final states were taken in the simulations. Here we used the climbing image NEB method52 to calculate not only the minimum energy path but also the transition state configuration at the saddle point. The method drives the image with the highest energy up to the saddle point, by maximizing its energy along the band connected to the image and minimizing in all other directions. This way, the method allows us to accurately calculate the saddle point energy and configuration with a small number of images needed in the calculation. The NEB simulations were conducted with fixed lattice vectors. The database created for the migration barrier values at different lattice strains and the defect configurations was provided as input for kinetic Monte Carlo (KMC) simulations. KMC calculations were performed to quantify the effective diffusivity of oxygen in YSZ as a function of strain and temperature. The diffusivity calculated in this work is the oxygen tracer diffusion constant. We expect that the lattice strain has a similar order of magnitude impact on the mobility of an oxygen ion in both the chemical and the tracer diffusion constants, especially in the range of low vacancy concentrations.53 In the simulation model, oxygen atoms were assumed to hop to the nearest neighbor site on the oxygen sublattice through the path shown in Fig. 1. The rate of hops in the migration paths were obtained by the Boltzman relationship as, vAB ¼ n0 expEAB kBT (1) where EAB, found by the DFT-NEB simulations, denotes the migration barrier across the A–B cation pair at sites 3 and 4 in Fig. 1. The attempt frequency was fixed at n0 ¼ 1013 s1, which is appropriate for most metal oxide systems.50 The periodic boundary conditions were applied to all three dimensions of the simulation cell. The vacancy diffusivity, Dv was extracted from mean square displacement hR2i calculated in the KMC simulation using the Einstein relation, hR2i ¼ 6Dvt. (2) The oxygen diffusivity DO was then obtained using Dv considering the balance between the fluxes of oxygen atoms and vacancies as, Fig. 1 Migration path of the oxygen vacancy. Positions 1–6 are occupied by the cations Zr or Y. This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4811 DO ¼ cv 1 cv Dv (3) where cv is the vacancy concentration fraction. Dv exponentially depends on the effective barrier for oxygen vacancy migration, Eeff, as, Dv ¼ dv 0 exp Eeff kBT (4) where, d0v is a constant. At small values of cv, D0 scales almost linearly with cv. On the other hand, the change in the migration barrier, Eeff, contributes exponentially to Dv and D0. Results and discussion This section consists of five sub-sections, I–V. Sub-sections I–II demonstrate that the YSZ model in this work captures the reported vacancy and cation interactions in stabilizing the cubic structure in YSZ. The results in I and II are important for obtaining reliable values of the energy barriers for vacancy migration in YSZ. Sub-section III reports the effect of defect configurations in the vicinity of the migrating vacancy on the migration energy barrier. This proves that, in order to assess the effect of strain on migration barriers, one must consider a realistic subset of defect configurations in YSZ, and we report these in sub-section IV. Sub-section IV in particular illustrates the mechanism by which strain alters the diffusion barriers for oxygen vacancies in YSZ. Finally sub-section V quantifies the effective increase in oxygen tracer diffusivity as a function of biaxial lattice strain, based on the subset of defect configurations considered in this work for YSZ. I. Defect structure in YSZ: oxygen vacancy and yttrium distribution To determine the favorable location of vacancies with respect to the Y cations in the YSZ model, we tested a number of defect configurations. These models constitute 30 Zr, 2 Y and 63 O atoms, corresponding to a low doping ratio, 3% Y2O3, allowing only one vacancy in the simulation cell, in order to eliminate the vacancy–vacancy interaction in this part of the simulations. At this low Y2O3 doping, the YSZ oxygen sublattice cannot keep its cubic structure.34 Therefore, only in this part of the simulations, the atoms beyond the second neighbor cation shell to Y cations and beyond the third neighbor oxygen shell to the vacancy were fixed to retain the cubic phase.34 The comparison of the energies for the different YSZ models and defect configurations is presented in Fig. 2. The energies are shifted such that the lowest energy found is zero. The results show that the total energy can vary by up to 0.5 eV when the Y cation is located at the third or first neighbor site to the vacancy, instead of the second. This result is consistent with the prior reports on cation–vacancy interactions in YSZ, as discussed in the Introduction. Since this large energy difference arises from local configuration differences, the result suggests that the cation configuration surrounding the diffusion path of the oxygen–vacancy pair should be taken into account for accurate evaluation of the migration barriers. II. Distortion from the cubic lattice sites in YSZ The atoms are not stabilized at the exact fluorite structure lattice centers in YSZ, and large displacements of the cations and the anions from the fluorite lattice sites are present.55–57 This is in part due to the relaxation of the nearest neighbor oxygen anions and cations towards and away from the vacancy, respectively.32 For an acceptable YSZ structure in this study, it is important to capture this deviation quantitatively consistently with past experimental and simulation reports. For this purpose, we evaluated the dependence of the atomic displacements from the ideal fluorite lattice sites on the local defect structures in a set of candidate YSZ models. The fully relaxed structure of the models constituted 9% Y2O3 doped YSZ, made of 26 Zr, 6 Y and 61 O atoms in the simulation cell. The relaxation of the atoms were not constrained in this and the following simulations because the cubic structure is stable in the 9% Y2O3 doped YSZ, as long as Fig. 2 Energy as a function of the neighboring location of the oxygen vacancy with respect to the Y cations in YSZ. The energies were calculated for different cation/vacancy distributions, and plotted such that the lowest energy found is zero. Fig. 3 Fully relaxed atomic configurations of the 9% yttria doped zirconia with (a) nearest neighbor Y–Y pairs, and (b) uniformly distributed Y. Yellow, gray and red spheres indicate Y, Zr and O ions. Large distortion of the oxygen anion sublattice is evident when Y cations are uniformly spread in the YSZ without any short-range ordering of the Y cations. 4812 | J. Mater. Chem., 2010, 20, 4809–4819 This journal is ª The Royal Society of Chemistry 2010 the vacancies and dopant cations are distributed to enable the stability. Fig. 3 shows the relaxed atomic configurations for two models of YSZ. The first one (Fig. 3(a)) has first nearest neighbor (1NN) Y–Y pairs and the second one (Fig. 3(b)) has a uniform distribution of Y cations in the cell. In both models, the vacancies were placed as far away as possible from each other and located second nearest neighbor (2NN) to Y cations if allowed by the available space. The first model retains an overall cubic structure with oxygen anions displaced by 0.1–0.27 A from the ideal fluorite lattice sites. The extent of the displacements calculated here is well within the experimentally measured value of 0.36 A28 and prior simulation results of 0.20–0.31 A.34 On the other hand, the second model shows large distortions, exemplified at the circled region in Fig. 3(b),and has 1.2 eV higher energy compared to the first one. This difference arises due to local configuration differences, and thus, can affect the vacancy migration barrier in YSZ which ranges from 0.2 to 1.4 eV.25 The lack of stability of the cubic oxygen sublattice in the second model is due to the limited available space to accommodate the relative positions of the vacancy and Y cations at the favorable distances when the Y atoms are distributed uniformly into the model. In the circled region in Fig. 3(b),a vacancy is located at the 1NN position to the Y cations. This is because the 2NN position for the vacancy is not available to all the Y atoms when the Y atoms are distributed uniformly in the model. The restriction of some vacancies to the unfavorable 1NN positions with respect to the Y in YSZ destabilizes the cubic phase. On the contrary, the first model configuration with the Y–Y pairs introduces more space for the vacancies to be placed at the 2NN positions with respect to the Y cations. This enables more stability on the cubic oxygen sublattice. This result is consistent with the locally ordered nature of Y cations in YSZ at high concentrations of Y2O3, as found by Predith et al.,43 even if the YSZ has only 9% Y2O3, here. While the effect of Y–Y pairing on the migration barrier in a cubic YSZ is small compared to that of the vacancy–vacancy and vacancy–cation interactions,36 it is important to consider it to enable a stabilized cubic fluorite structure of YSZ in the model. III. Vacancy migration barriers in YSZ: effect of the local defect structure As implied in the previous section, the local defect structure around the oxygen vacancy migration path in YSZ affects the stability of the structure and the migration barrier,54 and this in turn can determine the favorable paths for vacancy to migrate through. To illustrate this, we present the differences in the migration barrier in a model with 1NN Y–Y pairs, which lower the energy and enable the stabilization of cubic phase. The consequent migration paths considered were A / B / C in Fig. 4(a). All six nearest neighbor cations (1–6 as shown in Fig. 1) are the same, Zr, for each migration step in A/B/C. On the other hand, the distribution of vacancies and Y cations near A /B and B/C differ beyond the nearest neighbor cations. This enables to investigate the variation of the migration barrier as a function of the defect structure near the migration path beyond the nearest neighbor cations. Fig. 4(b) shows the relative energy in the NEB trajectories for the migration through A/B/C. The resulting migration barriers range from 0.55 to 0.95 eV, even though all the cations in the nearest neighbor position to the migration path are the same. This significant difference exemplified in the migration barriers in A/B/C indicates that it is important to systematically take into account the longer-range interactions between the migrating vacancy and its surrounding cations and vacancies in YSZ. We evaluated a subset of defect configurations near the migration path to form a representative database of migration barriers in YSZ. The five different configurations, models A–E, that yield ten forward and backward jump barriers considered here are shown in Fig. 5. This subset aims to capture the isolated Fig. 4 (a) The vacancy positions and local structure considered for the vacancy migration path, A / B / C. Large black spheres indicate positions of the vacancies. Yellow, gray and red spheres indicate Y, Zr and O ions. Numbers of Y cations at the second nearest neighbor position to the vacancy is two at sites A and B, and one at site C. (b) NEB barriers calculated along the A / B and B / C paths. The variation in the migration barrier is clear in (b),even though all the six nearest neighbor cations in the migration path are Zr. The difference in barriers is due to the dissimilar vacancy and Y distribution near the A / B and B / C paths beyond the nearest neighbor cations. Fig. 5 Models A–E with varying defect configurations surrounding the migration path. ‘‘# of 2NNY’’ denotes the number of second nearest neighbor Y cations to the vacancy positions. ‘‘NN vacancy position’’ denotes that the nearest neighbor vacancy to the migrating vacancy is located at n-th nearest neighbor position on the anion sublattice (n-th ¼ 3rd or 4th in these models). Lower panel indicates the 2NN Y position, and the 3rd and 4th nearest neighbor vacancies on the anion sublattice with respect to the migrating vacancy V. This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4813 and coupled effects of vacancy–vacancy and vacancy–cation interactions in the local environment around the migration path. The vacancy–vacancy interaction is considered up to the fourth nearest neighbor on the anion sublattice, which was the maximum that could be accommodated in our 2 2 2 unit cell model when the nearest neighbor vacancy pairs are avoided. The cation–vacancy interaction was taken up to the second nearest neighbor. First, for considering the effect of the nearest neighbor cations in the migration path (sites 1–6 in Fig. 1),we only included the variation of the bridging cation pairs (sites 3 and 4 in Fig. 1). This was shown to be the most important contribution to the vacancy migration in YSZ50 when the long-range interactions were not considered. In models A, B, and E, the vacancy migrates through the Zr–Zr cation pair, and in models C and D, through the Zr–Y cation pair in the midst of the path. All the other cation sites closest to the migration path (sites 1, 2, 5 and 6 in Fig. 1) are Zr. The path through Y–Y pair was excluded from the simulations since it has a much higher migration barrier compared to the Zr–Zr and Zr–Y pairs,25,50,54 making it significantly less likely to migrate through the Y–Y cation pair. Second, for vacancy–vacancy interaction, models A and C consider the position of nearest neighbor vacancies, denoted as ‘‘NN vacancy position’’ in Fig. 5, to the initial and the final position of the migrating vacancy. The distance between the migrating vacancy and the nearest neighbor vacancy changes from the fourth to the third on the anion sublattice upon migrating from the initial to the final position (and vice versa in the backward migration). # of 2NN Y are the same in models A and C. Third, for vacancy– cation interaction beyond the nearest neighbor cations, models B and D consider the number of Y cations that are at the second nearest neighbor (2NN) position with respect to the initial and the final state of the vacancy. This is denoted as ‘‘# of 2NN Y’’ in Fig. 5. The number of 2NN Y cations changes from two to three upon migrating from the initial to the final state (and vice versa in the backward migration). Lastly, in model E, both the number of 2NN Y cations and the distance to the nearest vacancies at the initial and the final state were varied. The former changes from three to two and the latter from the fourth to the third nearest neighbor (and vice versa in the backward migration). In creating the database of migration paths and barriers, a model to represent a symmetric defect configuration around the initial and final states of the migrating vacancy was not possible. This is because of the computational limitation in the simulation cell size used in our DFT/NEB calculations. In computing the effective diffusivity in the following sections, we will show that we represented the barrier for a symmetric migration path to be the average of the barriers from models A–E. A simulation unit cell of 2 2 2 dimensions, containing 26 Zr, 6 Y and 61 O atoms corresponding to a 9% yttria doped zirconia was used in the calculation of the migration energy barriers for each of the models A–E described in Fig. 5. The oxygen vacancies and the Y cations were distributed in the cell satisfying the defect interactions discussed in the previous sections. Atomic configurations and the cell vectors were fully relaxed such that the total energy was minimized and the net stress acting on the cell was zero. The migration barriers for each model path were calculated using the climbing NEB technique. The calculated barriers for each path (A–E) are shown in Fig. 6. The energy values are relative to the initial state configuration. The reaction coordinate of each NEB image was normalized with respect to the distance between the initial and the final state of the vacancy. The oxygen vacancy migration process requires climbing over the highest energy (the saddle point) between the initial and the final state. E!B and E B in Fig. 6 are the migration barriers for the forward and the backward jumps, respectively, corresponding to the migration from the initial to the final positions of the vacancy shown in Fig. 5, and vice versa. As a result, the migration configurations in models C and D have higher barrier than those in models A, B and E. This is primarily because of the difficulty in passing through the Zr–Y pair compared to that through Zr–Zr pair. The forward jump barrier is higher than the backward jump barrier for models A and C, and lower for models B and D. For models A and C, this is because the distance from the migrating vacancy to the nearest neighbor vacancy is shorter at the final state, making it more difficult for the migrating vacancy to approach to the final state, primarily due the electrostatic interaction of vacancies. For models B and D, the number of Y cations at the second nearest neighbor position to the migrating vacancy is larger at the final state, making this final state energetically more favorable, and thus the reduced migration barrier in the forward migration. In model E, the forward jump barrier is higher than the backward one. Considering the vacancy–cation interactions, the final state in model E is energetically more favorable than the initial state because of the presence of a larger number of 2NN Y cations to the migrating vacancy at the final state. This suggests the migration barrier for the forward jump to be lower than the backward one. However, the result shows the opposite behavior, because the migrating vacancy is actually getting closer to its nearest neighbor vacancy at its final state. This indicates that the vacancy–vacancy interaction dominates the vacancy–cation interaction in determining the energy barriers in model E. All implications related to the defect–defect interactions leading to the resulting migration barriers in models A–E in the unstrained state of YSZ are consistent with those suggested by previous Fig. 6 NEB paths and the corresponding energies in the migration path calculated for models A–E. The relative energy and the coordinates are normalized with respect to the initial energy, and the distance between the initial and the final state, respectively. The saddle point energy is the highest energy state between the initial and the final states, and E!B and E B stand for barriers in the forward and the backward jumps, respectively. 4814 | J. Mater. Chem., 2010, 20, 4809–4819 This journal is ª The Royal Society of Chemistry 2010 reports by Bogicevic et al.36 Therefore, the results presented here enable a reasonable model of YSZ for assessing the effect of biaxial lattice strain on the ionic mobility. IV. Migration barrier as a function of biaxial lattice strain In assessing the role of the biaxial xy-plane lattice strain on the oxygen vacancy transport, we investigated the corresponding changes in the migration barriers, and the underlying reasons for these changes on models A–E in Fig. 5. The xy-plane strains considered here were 0.00, 0.02, 0.04, 0.06 and 0.08. The range of strain considered in this study is far above the experimentally measured fracture strain of a single-crystal YSZ that is less than 0.5%.58 However, large strain in brittle ceramics could be achieved by creating coherent hetero-interfaces at the nanoscale. A maximum 7% tensile lattice strain was suggested to exist in the YSZ thin film interfaces coherently strained between STO layers due to the mismatch between STO and YSZ.8 On the other hand, ability of thin-film (10–30 nm) YSZ to sustain this large tensile strain is unclear. We considered only tensile strain to increase the interatomic distances in YSZ, allowing more space for the oxygen migration. Results for model A are shown in Fig. 7. The migration barrier monotonically reduces with the increase in the biaxial lattice strain. The barrier in the forward jump decreases from 0.47 eV at 3 ¼ 0.00 to 0.07 eV at 3 ¼ 0.04, suggesting a significant enhancement of the vacancy diffusivity through this path under strain. Beyond this strain level, the vacancy was not stable at the initial site and favored to form at the defined final site in model A, with no barrier in the suggested transition path. This result implies that the defect interactions influencing the stability and energetics can also change with the strain state of YSZ. While not explicitly included in our DFT calculations, we implicitly took into account the effect of the strain on the stable defect configurations in the KMC simulations, as explained in the following section. We considered both the forward and the backward migration barriers, and the system has a higher probability of taking the lower energy migration paths favoring the more stable structures. The picture is not as straight-forward in models B–E as in model A that shows a monotonic reduction in the migration barrier with strain. Depending on the extent of the lattice strain, the consequent barrier either decreases or increases in the configurations represented by models B–E, as shown in Fig. 8. The lack of a straight-forward trend in the magnitude of migration barrier as a function of lattice strain suggests that there is more than one mechanism acting simultaneously upon strain in YSZ. Here, we explain the competing effects of strain on the vacancy migration barrier in YSZ from an atomistic and electronic perspective. The oxygen vacancy migration process involves bond-breaking and bond-making between the migrating oxygen and the surrounding cations. Therefore, we investigated the changes in the electronic density distribution of the valence electrons describing the Zr, Y and O atoms near the migration path to explain the local effects of strain in YSZ. We focused on the two planes in the migration path as shown in Fig. 9. The charge distribution on the ‘‘oxygen plane’’ in Fig. 9 correlates to the ‘‘space available’’ for the oxygen, denoted as O, to migrate along the prescribed path. The migration space is also related to the cation-cation distance midst the migration path. On the other hand, the ‘‘bonding plane’’ in Fig. 9 shows a measure of the ‘‘bonding strength’’ between the oxygen and the nearest cations at the initial site, quantified in terms of the oxygen–cation ‘‘bond thickness’’. It is this set of bonds, shown between the oxygen denoted as O and the cation denoted as C, which are broken in the initiation of oxygen diffusion. Both the migration space and the bond strength change when YSZ is subjected to biaxial strain. Fig. 10 illustrates an example in which the increase in strain increases the open space, and also increases the O–C bonding strength to the extent of local inelastic relaxation near the initial state of oxygen. In this case, the lattice strain is beyond just Fig. 7 (a) NEB paths and the corresponding energies in the migration path calculated for model A. The relative energy and the migration coordinates are normalized with respect to the initial energy, and the distance between the initial and the final state, respectively. The saddle point energy is the highest energy state between the initial and the final states. (b) The vacancy migration barrier, EB, as a function of biaxial lattice strain. The filled circles and the open circles indicate forward and backward jump, respectively. The barrier decrease as the strain increases in the migration path in model A. Fig. 8 The vacancy migration barrier, EB, as a function of biaxial lattice strain for models B–E (in a–d). The filled circles and the open circles indicate the forward and the backward jump barriers, respectively. This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4815 elastically extending the oxygen–cation bond length; the oxygen, O, was actually attracted towards the cation C and broke its bond from the opposite cation upon local relaxation of atoms. As a result, a stronger binding of O to C arises, shown with the thicker O–C bond in Fig. 10. These two consequences of strain act oppositely on the migration barrier – the former (opening of migration space) favors the migration and the latter (local inelastic relaxations near the oxygen) inhibits the migration. A specific illustration of these competing consequences of biaxial strain in YSZ represented by model E is shown in Fig. 11. The resulting change in the migration space was quantified as the distance between the cations midst the migration path, D(Zr–Zr) in this configuration. The bond O–C in Fig. 11(b) is one of the oxygen–zirconium bonds that break upon the migration of oxygen. The large reduction in the forward migration barrier from 0.51 eV at 3 ¼ 0.00 to 0.10 eV at 3 ¼ 0.06 as shown in Fig. 8(d) corresponds to the increase in the D(Zr–Zr) distance (shown in Fig. 11(a)) and the weakening of the O–C bond (shown in Fig. 11(b)). On the other hand, this trend is not monotonic. At higher strain, 3 ¼ 0.08, the cation–cation distance decreases and the oxygen–zirconium bond strengthens again, resulting in a significant increase of the migration barrier from 0.10 eV at 3 ¼ 0.06 to 0.47 eV at 3 ¼ 0.08. At this strain, the local relaxation of the atoms particularly re-strengthens the oxygen-zirconium bond at the initial site of the oxygen, also affecting the migration space. If assumed independent, the oxygen–cation bond strength is more strongly correlated to this large increase in the barrier at large strain seen in Fig. 8(d). The increase in the cation–oxygen bond strength is evident at 3 ¼ 0.08 in Fig. 11(b),but a very small change occurs in the D(Zr–Zr) in Fig. 11(a). While we present this set of results explicitly for the model E, the same qualitative correlation of the competing effects of the migration space and the cation–oxygen bond strength governed the migration barrier changes with strain consistently also in the other models A–D. V. Effective oxygen diffusivity as a function of biaxial lattice strain Models A–E constituted the database of the migration paths and barriers as a function of the biaxial strain in YSZ. Using this database to govern the overall oxygen diffusion, the macroscopic oxygen diffusivity in the bulk of YSZ was estimated by KMC simulations as described in the Simulation Approach section. A 15 15 15 simulation cell of 8% yttria doped zirconia was used, consistent with the doping ratio in YSZ reported in the Ref. 8 for the YSZ/STO heterolayers. Y cations and vacancies were distributed randomly in the KMC simulation cell. At a first sight, this configuration contradicts with the implications for the presence of 1NN Y–Y pairs in YSZ as discussed in Introduction. On the other hand, the rules in selecting the diffusion barriers in the KMC simulations were taken consistently with the local defect configurations, including the effect of Y–Y pairs as evaluated in the DFT-NEB calculations in this work. Therefore, we considered this model in KMC to be reasonable for assessing the overall oxygen diffusivity in YSZ. Supporting this approach, past DFT and KMC studies of oxygen diffusivity in YSZ with randomly distributed Y atoms and vacancies were shown to give Fig. 9 The electronic density distribution on the ‘‘oxygen plane’’ and the ‘‘bonding plane’’ at the initial and the saddle states. The migrating oxygen, O, moves towards the open migration space, and the bond between oxygen O and cation C breaks in the migration. Fig. 10 Change in the electronic density distribution on the oxygen plane and the bonding plane under biaxial lattice strain, 3, from 0.00 to 0.07. Gray and red spheres represent Zr and O atoms, respectively. 4816 | J. Mater. Chem., 2010, 20, 4809–4819 This journal is ª The Royal Society of Chemistry 2010 the correct yttria doping ratio for the maximum oxygen diffusivity.25,50 The effect of the local defect environment around the migration path on the barriers input to the KMC simulations was taken into account as follows: Definitions. E!a and E a are the forward and backward jump energy barriers in each model a (a ¼ A, B, C, D, E),respectively. ri v and rf v stand for the distance of the first nearest neighbor (1NN) vacancy to the migrating vacancy at the initial and final state, respectively. ni Y, and nf Y are the number of second nearest neighbor (2NN) Y cations to the migrating vacancy at the initial and final state, respectively. Rules. The rules in determining the energy barrier for each hop in the KMC simulations treat the local defect configurations in a coarse grained fashion, based on the models A–E (Fig. 5) assumed to govern the oxygen vacancy diffusion in YSZ. In models A and C, the vacancy forward hopping (from initial to final state) results in a decrease of the 1NN vacancy distance to the migrating vacancy, that is ri v > rf v. Then, in the KMC simulation, E!a and E a (a ¼ A, C) are chosen for the cases of ri v > rf v, and ri v<rf v, respectively. In models B and D, the vacancy forward hopping (from initial to final state) results in an increase in the number of 2NN Y to the migrating vacancy, that is ni Y < nf Y. Then, in the KMC simulation, E!a and E a (a ¼ B, D) are chosen for the cases of ni Y <nf Y and ni Y > nf Y, respectively. The explicit definition of the rules implemented in the KMC simulations is as follows: (1) For the vacancy migration through a Zr–Zr pair midst the migration path, we select the barrier among E!A, E A E!B and E B according to the following criteria: If ri v > rf v, select E!A; and if ri v < rf v, select E A. If ri v ¼ rf v and ni Y > nf Y, select E!B; or if ri v ¼ rf v and ni Y < nf Y, select E B. If ri v ¼ rf v and ni Y ¼ nf Y, randomly choose barriers from E!A, E A, E!B, and E B. (2) For the vacancy migration through a Zr–Y pair midst the migration path, we select the barrier among E!C, E C, E!D and E D according to the following criteria: If ri v > rf v, select E!C; if ri v < rf v, select E C. If ri v ¼ rf v and ni Y > nf Y, select E!D; or if ri v ¼ rf v and ni Y < nf Y, select E D. If ri v ¼ rf v and ni Y ¼ nf Y, randomly choose barriers from E!C, E C, E!D and E D. (3) Vacancy migration through Y–Y pair is prohibited. When ri v ¼ rf v and ni Y ¼ nf Y, the migration path is surrounded by a symmetric distribution of defects. This case was not possible to account for in the DFT-NEB calculations due to the computational limitations on the simulation size. We assumed such a symmetric path to be governed by a migration barrier that is the average of the barriers from the non-symmetric models represented in A–E here. Therefore, the random selections among the entire set of E!a and E a (for a ¼ A, B, C, D, E) during the 108 KMC steps should provide sufficient statistics to represent the average migration barrier for the symmetric path, ri v ¼ rf v and ni Y ¼ nf Y. Moreover, the probability of having symmetric path in the KMC runs was less than 3%. Therefore, the error that may arise due to selecting an average migration barrier for the symmetric configuration is expected to have little significance on the effective diffusivity. As seen from the rules above, the energy barrier for each hop in the KMC simulations was chosen based on a coarse grained treatment of the local defect configurations around the migrating vacancy. We believe this approach gives a reasonable result, particularly based on the validity of the models A–E (Fig. 5) to govern the oxygen vacancy diffusion in YSZ. On the other hand, the reader can refer to a more accurate treatment of the diffusion barriers in KMC simulations using the cluster expansion method,43,59 which can account for the energy deviations due to the local atomic structures around each hop. In that case, the KMC model itself has to take into account the stable distributions of vacancies and cations (as in subsections I and II),rather than a random distribution, making the task exhaustive. The KMC simulations were performed as a function of temperature ranging as 400, 600, 800 and 1000 K, and of biaxial strain ranging as 0.00, 0.02, 0.04, 0.06 and 0.08. Ten simulations with different distributions of dopant cations and vacancies were conducted for each temperature and strain state. The resulting strain dependence of the effective oxygen tracer diffusivity DO is shown in Fig. 12. y-axis in the figure is on a relative scale with Fig. 11 (a) Change in the cation–cation distance, D(Zr–Zr),and (b) the bond thickness (strength) for O–C in the migration under biaxial lattice strain from 0.00 to 0.08 for the forward jump in model E (see also Fig. 8(d)). Increase in the migration space D(Zr–Zr),and the weakening of the O–C bond correspond to the decrease of the migration barrier in Fig. 8(d),and vice versa. Large gray spheres and small red spheres are Zr and O atoms, respectively. The isoline interval is 0.5 A3. This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4817 respect to the oxygen diffusivity, D0 O at 3 ¼ 0.00 for the lowest temperature simulated, 400 K. The reported results show the average values and the error bars over the ten simulations. The result is not a simple monotonic increase of diffusivity with increasing strain in the entire range considered. The oxygen diffusivity exhibits an exponential increase up to a critical value, 4% of tensile strain. This increase is more significant at the lower temperatures. The maximum enhancement in diffusivity compared to the unstrained state in YSZ is 6.8 103 times higher at 4% strain and at 400 K. At the higher strain states, the diffusivity decreases. This is attributed to the local relaxations beyond a limit of elastic bond strain, resulting in the strengthening of the local oxygen–cation bonds that increases the migration barrier. Our first-principles-based simulation results indicate that a biaxial strain in YSZ induced by a lattice-mismatch coherent interface can significantly increase the diffusivity up to a critical strain state, and worsen it beyond that strain state. One could view this critical strain state as the optimal strain to induce in the material to attain the maximum benefit of the biaxial lattice strain on the oxygen diffusivity. We note that a recent prediction by Schichtel et al.22 using a model based on the elastic properties of YSZ estimated 2.5 orders of magnitude increase in the ionic conductivity of YSZ at 7% strain at 573 K. The model of Schichtel et al. relies on the isotropic pressure induced by the strain to enhance the diffusivity, and does not take into account the competing relaxations at these large tensile strain states in YSZ. Therefore, it cannot capture the local bond-strength changes that may mitigate the diffusivity at large strain values. In spite of the differences in our and Schichtel’s model, the magnitude of the maximum increase in the oxygen diffusivity in strained YSZ predicted by these models range from 2.5 to 3.8 orders of magnitude. Therefore, the lattice strain (elastic and inelastic) alone cannot explain the eight orders of magnitude increase of ionic conductivity in YSZ reported by Barriocanal et al.8 These results do not take into account the space-charge effects in altering the anion concentrations near the YSZ interface. On the other hand, the contribution to the increase in the anion diffusivity due to a possible space-charge-induced increase in the vacancy concentration is almost linear at low concentration of vacancies (see eqn (3)) as in the 9% Y2O3 doped YSZ. Therefore, this contribution is expected to be small compared to the exponential increase found here due to the applied biaxial lattice strain.60 Conclusion This work provides a mechanistic and systematic assessment of the effect of biaxial lattice strain on the oxygen diffusivity in YSZ, independently from the defect chemistry near a heterointerface with a dissimilar oxide, e.g. YSZ/STO as reported by Barriocanal et al.8 Using DFT/NEB simulations, we uncovered the underlying microscopic mechanism by which the biaxial tensile strain acts on the oxygen vacancy transport in YSZ. The migration barrier correlates fundamentally with two competing processes acting in parallel in the presence of lattice strain: (1) the migration space, quantified as the cation–cation distance midst the migration path for the oxygen, and (2) the strength of the bond between the oxygen and the nearest cation preventing the oxygen from migrating. Increasing of the migration space and weakening of the oxygen–cation bond contribute to the decrease of the vacancy migration barrier, and vice versa. In assessing the effective diffusivity as a function of strain in YSZ, we constructed a data base of migration paths and barriers to account for the important effects of local distributions of vacancies and dopant cations near the migration path. KMC simulations using this database showed that the macroscopic oxygen diffusivity increases exponentially in YSZ up to a critical value of biaxial tensile strain, with a more significant impact at the lower temperatures. Beyond this critical value, the strain state in YSZ worsens the diffusivity due to the local relaxations that strengthen the oxygen–cation bond at the initiation of migration. One could view this critical strain state as the optimal strain, or the fastest strain, to attain the highest benefit of the biaxial lattice strain on the oxygen diffusivity in the material system. The maximum enhancement in diffusivity in YSZ was 6.8 103 times higher at 4% strain and temperature of 400 K. This is a significant and desirable increase. Nevertheless, it is far below the eight orders of magnitude increase reported by Barriocanal et al. for the thin-film YSZ layers coherently strained between STO layers.8 Therefore, we believe that the lattice strain alone cannot be responsible for the colossally high ionic conductivity suggested to occur in YSZ in ref. 8. Nevertheless, given our results, controlling the lattice strain by direct mechanical load or by creating a coherent hetero-interface to induce an optimal strain state can enable desirably high ionic conductivity at reduced temperatures in the fluorite material, the YSZ in this case. The insights gained here for the role of lattice strain on the local bonding structure and ionic charge transport process are of importance for modulating the transport properties in a variety of solid state conducting material applications, beyond the specific SOFC electrolyte material studied in this work. Interesting examples may extend to Li-ion batteries, oxygen gas sensors, semi-conductors and photovoltaics. We believe that the mechanistic results obtained in this study, particularly on the Fig. 12 The relative enhancement of the oxygen diffusivity calculated by the kinetic Monte Carlo simulations at different temperatures and strain states. D0 O is the oxygen diffusivity at 3 ¼ 0.00 at 400 K. The diffusivity increases exponentially with strain, and shows a maximum at (or near) 3 ¼ 0.04 followed by a decrease at higher strains. The increase in diffusivity is more significant at lower temperatures. 4818 | J. Mater. Chem., 2010, 20, 4809–4819 This journal is ª The Royal Society of Chemistry 2010 non-linear and competing consequences of strain on the bonding structure in the vicinity of the ionic carrier species, provides insights for creating functional strained nano-structures for ionic and electronic conducting devices, including but not limited to SOFCs. Acknowledgements We thank the Nuclear Regulatory Commission for financial support, and the National Science Foundation for computational support through the TeraGrid Advanced Support Program under Grant No. TG-ASC090058. References 1 J. C. Boivin and G. Mairesse, Chem. Mater., 1998, 10, 2870. 2 T. Suzuki, I. Kosacki and H. U. Anderson, Solid State Ionics, 2002, 151, 111. 3 I. Kosacki, C. M. Rouleau, P. J. Becher, J. Bentley and D. H. Lowndes, Solid State Ionics, 2005, 176, 1319. 4 J. H. Shim, C.-C. Chao, H. Huang and F. B. Printz, Chem. Mater., 2007, 19, 3850. 5 J. Maier, Nat. Mater., 2005, 4, 805. 6 N. Sata, K. Eberman, K. Eberl and J. Maier, Nature, 2000, 408, 946. 7 X. X. Guo, I. Matei, J.-S. Lee and J. Maier, Appl. Phys. Lett., 2007, 91, 103102. 8 J. G. Barriocanal, A. R. Calzada, M. Varela, Z. Sefrioui, E. Iborra, C. Leon, S. J. Pennycook and J. Santamaria, Science, 2008, 321, 676. 9 S. Ramanathan, J. Vac. Sci. Technol., A, 2009, 27, 1126. 10 A. J. Jacobson, Chem. Mater., 2010, 22, 660. 11 X. Guo, Science, 2009, 324, 465; J. G. Barriocanal, A. R. Calzada, M. Varela, Z. Sefrioui, E. Iborra, C. Leon, S. J. Pennycook and J. Santamaria, Science, 2009, 324, 465. 12 M. P. O’Callaghan, A. S. Powell, J. J. Titman, G. Z. Chen and E. J. Cussen, Chem. Mater., 2008, 20, 2360. 13 Y. Hoshina, K. Iwasaki, A. Yamada and M. Konagai, Jpn. J. Appl. Phys., 2009, 48, 04C125. 14 J. Munguıa, G. Bremond, J. M. Bluet, J. M. Hartmann and M. Mermoux, Appl. Phys. Lett., 2008, 93, 102101. 15 V. R. D’Costa, Y.-Y. Fang, J. Tolle, J. Kouvetakis and J. Menendez, Phys. Rev. Lett., 2009, 102, 107403. 16 M. Mavrikakis, P. Stoltze and J. K. Nørskov, Catal. Lett., 2000, 64, 101. 17 J. Kilner, Nat. Mater., 2008, 7, 838. 18 J. Maier, Prog. Solid State Chem., 1995, 23, 171. 19 A. K. Ivanov-****z, Crystallogr. Rep., 2007, 52, 129. 20 D. C. Sayle, J. A. Doig, S. C. Parker and G. W. Watson, Chem. Commun., 2003, 1804. 21 C. Korte, A. Peters, J. Janek, D. Hesse and N. Zakharov, Phys. Chem. Chem. Phys., 2008, 10, 4623. 22 N. Schichtel, C. Korte, D. Hesse and J. Janek, Phys. Chem. Chem. Phys., 2009, 11, 3043. 23 K. Suzuki, M. Kubo, Y. Oumi, R. Miura, H. Takaba, A. Fahmi, A. Chatterjee, K. Teraishi and A. Miyamoto, Appl. Phys. Lett., 1998, 73, 1502. 24 W. Araki, Y. Imai and T. Adachi, J. Eur. Ceram. Soc., 2009, 29, 2275. 25 R. Pornprasertsuk, P. Ramanarayanan, C. B. Musgrave and F. B. J. Prinz, J. Appl. Phys., 2005, 98, 103513. 26 W. L. Roth, R. Wong, A. I. Goldman, E. Canova, Y. H. Kao and B. Dunn, Solid State Ionics, 1986, 18–19, 1115; H. Morikawa, Y. Shimizugawa, F. Marumo, T. Harasawa, H. Ikawa, K. Tohji and Y. Udagawa, J. Ceram. Soc. Jpn, 1988, 96, 253; M. H. Tuilier, J. Dexpert-Ghys, H. Dexpert and P. Lagarde, J. Solid State Chem., 1987, 69, 153. 27 M. Weller, Z. Metallkd., 1993, 84, 6. 28 D. Steele and B. E. F. Fender, J. Phys. C: Solid State Phys., 1974, 7, 1. 29 B. W. Veal, A. G. McKale, A. P. Paulikas, S. J. Rothman and L. J. Nowicki, Physica B+C, 1988, 150, 234; D. Komyoji, A. Yoshiasa, T. Moriga, S. Emura, F. Kanamaru and K. Koto, Solid State Ionics, 1992, 50, 291. 30 C. R. A. Catlow, A. V. Chadwick, G. N. Greaves and L. M. Moroney, J. Am. Ceram. Soc., 1986, 69, 272. 31 P. Li and J. E. Penner-Hahn, Phys. Rev. B: Condens. Matter, 1993, 48, 10063; P. Li and J. E. Penner-Hahn, Phys. Rev. B: Condens. Matter, 1993, 48, 10074; P. Li and J. E. Penner-Hahn, Phys. Rev. B: Condens. Matter, 1993, 48, 10082. 32 J. P. Goff, W. Hayes, S. Hull, M. T. Hutchings and K. N. Clausen, Phys. Rev. B, 1999, 59, 14202. 33 K. Kawata, H. Maekawa, T. Nemoto and T. Yamamura, Solid State Ionics, 2006, 177, 1687. 34 G. Stapper, M. Bernasconi, N. Nicoloso and M. Parrinello, Phys. Rev. B: Condens. Matter Mater. Phys., 1999, 59, 797. 35 A. Bogicevic, C. Wolverton, G. M. Crosbie and E. B. Stechel, Phys. Rev. B: Condens. Matter Mater. Phys., 2001, 64, 014106. 36 A. Bogicevic and C. Wolverton, Phys. Rev. B: Condens. Matter Mater. Phys., 2003, 67, 024106. 37 A. Bogicevic and C. Wolverton, Europhys. Lett., 2001, 56, 393. 38 X. Xia, R. Oldman and R. Catlow, Chem. Mater., 2009, 21, 3576. 39 M. Sakib Khan, M. Saiful Islam and D. R. Bates, J. Mater. Chem., 1998, 8, 2299. 40 M. O. Zacate, L. Minervini, D. J. Bradfield, R. W. Grimes and K. E. Sickafus, Solid State Ionics, 2000, 128, 243. 41 D. A. Andersson, S. I. Simak, N. V. Skorodumova, I. A. Abrikosov and B. Johansson, Proc. Natl. Acad. Sci. U. S. A., 2006, 103, 3518. 42 S. Garcıa-Martın, D. P. Fagg and J. T. S. Irvine, Chem. Mater., 2008, 20, 5933. 43 A. Predith, G. Ceder, C. Wolverton, K. Persson and T. Mueller, Phys. Rev. B: Condens. Matter Mater. Phys., 2008, 77, 144104. 44 F. Pietrucci, M. Bernasconi, A. Laio and M. Parrinello, Phys. Rev. B: Condens. Matter Mater. Phys., 2008, 78, 094301. 45 G. Kresse and J. Hafner, Phys. Rev.B: Condens. Matter, 1993, 47, 558. 46 G. Kresse and J. Furthm€uller, Phys. Rev. B: Condens. Matter, 1996, 54, 11169. 47 P. E. Bl€ochl, Phys. Rev. B: Condens. Matter, 1994, 50, 17953. 48 G. Kresse and D. Joubert, Phys. Rev. B: Condens. Matter Mater. Phys., 1999, 59, 1758. 49 J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett., 1996, 77, 3865. 50 R. Krishnamurthy, Y.-G. Yoon, D. J. Srolovitz and R. Car, J. Am. Ceram. Soc., 2004, 87, 1821. 51 G. Mills and H. Jonsson, Phys. Rev. Lett., 1994, 72, 1124. 52 G. Henkelman, B. P. Uberuaga and H. Jonsson, J. Chem. Phys., 2000, 113, 9901. 53 A. Van der Ven, G. Ceder, M. Asta and P. D. Tepesch, Phys. Rev. B: Condens. Matter Mater. Phys., 2001, 64, 184307. 54 M. Martin, J. Electroceram., 2006, 17, 765. 55 F. Frey, H. Boysen and I. K- Bischoff, Z. Kristallogr., 2005, 220, 1017. 56 D. N. Argyriou, M. M. Elcombe and A. C. Larson, J. Phys. Chem. Solids, 1996, 57, 183. 57 N. Ishizawa, Y. Matsushima, M. Hayashi and M. Ueki, Acta Crystallogr., Sect. B: Struct. Sci., 1999, 55, 726. 58 D. Baither, M. Bartsch, B. Baufeld, A. Tikhonovsky, M. R€uhle and U. Messerschmidta, J. Am. Ceram. Soc., 2001, 84, 1755. 59 A. Van der Ven, J. C. Thomas, Q. C. Xu, B. Swoboda and D. Morgan, Phys. Rev. B: Condens. Matter Mater. Phys., 2008, 78, 104306. 60 A. Kushima and B. Yildiz, ECS Trans., 2009, 25, 1599. This journal is ª The Royal Society of Chemistry 2010 J. Mater. Chem., 2010, 20, 4809–4819 | 4819 heres the first one..more to come.. [/QUOTE]
Insert quotes…
Verification
Who makes the popular shade guide?
Post reply
Forums
Lab talk, the good, the bad, and the ugly
Zirconium
not easy to explain. or understand.. or accept...
Top
Bottom