Adf 8
Adf 8
Adf 8
Density functional theory of cubic zirconia and 6–15 mol% doped yttria-stabilized
zirconia: structural and mechanical properties
Abstract: Results of ab-initio density-functional theory calculations within the generalized gradient approximation
(GGA-PBE) of structural and mechanical properties of cubic zirconia and yttria-stabilized zirconia (YSZ) with yttria
(Y 2 O 3 ) concentrations of 6.67, 10.34, and 14.28 mol% are reported. It is found that the calculated structural and
mechanical parameters of all considered structures are highly consistent with the existing experimental data and the other
theoretical values. The doping concentration of yttria-stabilized zirconia has critical importance in ionic conductivity and
stabilization of high temperature down to room temperature. The lattice parameter and cell volume linearly decrease
with increases in doping concentration. Moreover, the effects of doping of yttria-stabilized zirconia on elastic constants
are studied. Elastic constants of cubic zirconia and 6–15 mol% doped yttria-stabilized zirconia are calculated using the
strain-stress approach and linear response theory.
Key words: First-principle calculations, DFT, structural properties, mechanical properties, cubic zirconia, yttria-
stabilized zirconia
1. Introduction
Renewable energy plays an important role in providing sustainable energy to meet demand. Fuel cell technology
has been gaining popularity amongst renewable energy technology in recent years [1–3]. A fuel cell is an
electrochemical device that converts electricity from a chemical reaction. Fuel cells have some advantages; to
name a few: their significant environmental benefits, high efficiency converting chemical energy into usable form,
and hence expanding the alternatives of providing sustainable energy to meet demand [4]. Solid oxide fuel cells
(SOFCs) are quite likely to become commercially viable among fuel cell technologies because of operating at
intermediate temperature (400–700 ◦ C). SOFCs use a solid ceramic electrolyte, such as stabilized yttria zirconia,
instead of a liquid or membrane. Zirconium oxide (ZrO 2 ) is known, as zirconia has been of great scientific and
technological material interest in gas sensors, high-temperature electrolysis, thermal barrier coatings, solid oxide
fuel cells, etc. ZrO 2 exists in three phases: monoclinic zirconia at temperatures less than 1170 ◦ C, tetragonal
from 1170 ◦ C to 2370 ◦ C, and cubic from 2370 ◦ C to 2680 ◦ C polymorphs at ambient pressure [5]. The
highest ionic conductivity has been shown in the cubic phase [6]. Cubic zirconia can be stabilized by addition of
aliovalent oxides such as Y 2 O 3 at ambient temperature. The addition of yttria to pure zirconia replaces some
Zr +4 ions with Y +3 ions that create charge compensating oxygen vacancies for balancing the valence charge.
∗ Correspondence: berna.akgenc@klu.edu.tr
223
AKGENÇ and ÇAĞIN/Turk J Phys
In Kröger–Vink notation:
′
R2 O 3 −
−−→ 2RZr + V Ö +3 O x
O (1)
ZrO2
2. Calculation methods
The first-principles calculations are performed by the Vienna ab-initio simulation package program (VASP),
based on density functional theory (DFT) with the generalized gradient approximation (GGA) using a plane-
wave basis set [14]. We have started with calculating the structural properties of the cubic form of ZrO 2 .
The Zr +4 atoms have been described by twelve (4s 2 4p 6 4d 2 5s 2 ) , Y +3 atoms by eleven (4s 2 4p 6 4d 1 5s 2 ), and
O −2 atoms by six (2s 2 2p 4 ) valence electrons, respectively. Kinetic energy cut-off for plane waves, 600 eV,
and a 4x4x4 k-point mesh for the total energy calculations have been determined as optimized values. The
Monkhorst–Pack method is used for k-mesh. The geometrical relaxation calculations are also performed with
the following common parameters: convergence criterion for energy is 10−5 eV ; the convergence criterion for
the maximal force between the atoms is 0.005 eV/Å; the maximum ionic steps are 200. The crystal structure
of the cubic form of ZrO 2 is the calcium fluoride structure with space group Fm3m (space number 225). It
is based on the fcc lattice and basis with Zr atom (0,0,0) and O’s at (0.25,0.25,0.25) and (0,75,0.75,0,75) to
generate 4 formula units of ZrO 2 in the primitive unit cell. The conventional bulk cell of the fluorite structure
consists of four zirconium and eight oxygen ions. Yttria stabilized zirconia (YSZ) is generated by starting with
the cubic form of ZrO 2 in Figure 1.
The doped yttria on cubic zirconia not only results in high ionic conductivity but also in the stabilization
of the high-temperature fluorite structure down to ambient temperature. To get doped 6–15 mol% yttria of
zirconia, we have studied several sizes of super cells: 2 × 1 × 1 (Model 1- 24 atoms), 2 × 2 × 1 (Model
2- 48 atoms) and 2 × 2 × 2 (Model 3- 96 atoms). The addition of yttria to the pure zirconia for each
224
AKGENÇ and ÇAĞIN/Turk J Phys
Figure 2. The crystal structure of 14.28 mol% Y 2 O 3 doped YSZ for 2 × 1 × 1 (a), 6.67 mol% Y 2 O 3 doped YSZ for
2 × 2 × 1 (b), 10.34 mol% Y 2 O 3 doped YSZ for 2 × 2 × 2 (c).
The ionic relaxation is performed and the equilibrium volume of the primitive cell is calculated at zero pressure
for all structures and configurations.
The structure of 2 × 1 × 1 supercell is created by two Y +3 ions interstitially replaced by two Zr +4
ions with one oxygen vacancy. According to 14.28 mol% yttria concentration of 24 atoms for 2 × 1 × 1
8!
unit cell, all the possible configurations could easily be studied ( 6!2! 28 different configurations, 8 different
interstitial sites for 2 yttrium atoms). We have found the configuration that gives the minimum energy among
225
AKGENÇ and ÇAĞIN/Turk J Phys
all possible configurations. We have calculated structural and mechanical properties for the configuration that
gives the minimum energy. With increasing number of supercells as 2 × 2 × 1 and 2 × 2 × 2 supercells, the
possible configurations of yttrium atoms in the crystal is hard to calculate. We have randomly created oxygen
vacancy for the rest of the crystal structures. Because of the limitation of atom numbers of quantum chemistry
calculations, increasing number of atoms in supercells is balanced by reducing the number of k-point sampling.
For optimal computational accuracy and efficient use of computational resources, we have performed 2 × 2 ×
2 k-point sampling for 2 × 2 × 2 supercell.
3. Results
3.1. Structural properties
The calculated total energies as a function of a primitive cell volume (12 atoms) for the cubic fluorite phase of
zirconia and 6–15 mol% Y 2 O 3 on cubic zirconia are used to determine structural and mechanical properties
by a least-squares fit of third-order Birch–Murnaghan equation, which is given as [12]
[ ]3 [( ) ]2 ⌈
( )2/3 2/3 ( )2/3 ⌉
9V0 B0 V0 ′ V0 V0
E (V ) = E0 + − 1 B0 + −1 6−4 , (3)
16 V V V
′
where E0 is the total energy, V0 is the equilibrium volume, B0 is the bulk modulus at 0 GPa pressure, and B0
is the first derivative of bulk modulus with respect to pressure in Figure 3.
Figure 3. Total energy as a function of conventional cell volume (12 atoms) for cubic fluorite phase of zirconia from
ab-initio calculations. The calculated values are shown by circles. The dashed line represents the fit of the data to the
Birch–Munganhan EOS.
Table 1 shows the structural parameters obtained after the fitting procedure with respect to pressure at
0 GPa. It also shows the calculated equilibrium energy (E 0 in eV), the lattice parameters ( a in Å), the bulk
modulus B in GPa, the pressure derivative of bulk modulus B‘, the Zr–O distance (dZr−O in Å) , and the Y–O
distance ( dY −O in Å) after relaxation. We find the lattice constant of pure zirconia to be 5.15 Å , which is
very close to the experimental value of 5.09 Å within the GGA-PBE approximation.
Vacancies or dopant/substitution atoms have critical importance in material properties for YSZ. While
substituent Y atoms decrease over all positive charge in the cubic form of zirconia, through the introduction of
oxygen vacancy models maintain charge neutrality. As the oxygen vacancies prefer to be close to the smallest
226
AKGENÇ and ÇAĞIN/Turk J Phys
Table 1. The calculated structural parameters of cubic zirconia and 6–15 mol% doped yttria-stabilized zirconia.
3
Structure E0 (eV) V0 (Å ) a (Å) B (GPa) B‘ dZr−O (Å) dY −O (Å)
–113.23 136.72 5.15 233.44 4.26 2.21 -
ZrO2
–91.80c 129.97d 5.07a 5.14c 209b 228c 4.43e 2.19d -
Y2 Zr14 O31
–111.52 140.39 5.19 151.08 5.96 2.97 2.41
(6.67 mol%)
Y6 Zr26 O61
–110.11 141.08 5.20 155.68 5.67 2.88 2.56
(10.24 mol%)
Y2 Zr6 O15
–108.42 141.36 5.21 14 1.44 6.15 2.42 2.48
(14.28 mol%)
a
Ref [15], b Ref [16], c Ref [10], d Ref [17], e Ref [18]
cation because of smaller Pauling radius (since Zr +4 has a Pauling radius of 0.80 Å and Y +3 of 0.93 Å,
respectively), the distance of cation atom (Zr, Y) and oxygen atom (O) decreases. Replacing Zr atom with
Y atom leads to an increase in volume of the cell at ground state. Since the atomic radius of yttrium and
zirconium is larger than the atomic radius of oxygen, distances of dZr−O and dY −O can change with creating
oxygen vacancy.
Table 2. The calculated elastic constants ( Cij in GPa) for cubic zirconia with available theoretical calculations and
experimental measurements.
In accordance with finite strain continuum elastic theory, the deformation corresponding to small strain
is imposed on a model crystal in a linear elastic manner. The calculated energy of the strained system may
then be expressed by a Taylor series expansion in terms of strain as follows:
∑
6
1 ∑6 1 ∑6
E (V, ε) = E (V0 , 0) + V0 σi εi + V0 Cij εi εj + V0 Cijk εi εj εk , (4)
i=1
2! i,j=1 3! i,j,k=1
227
AKGENÇ and ÇAĞIN/Turk J Phys
where ε and σ are the strain and stress tensors, respectively. E(V0 , 0), is the zero strain total energy and Cij
are the second order elastic constants.
As noted above, the elastic stiffness coefficients Cij can be obtained by straining the lattice. According
to the symmetry of the cubic crystal, it is completely described by three independent constants, C 11 , C 12 , and
C 44 , obtained by the three strains C1, C2, and C3 given below, respectively.
δ2
δ 0 0 δ 0 0 1−δ 2 0 0
C1 = 0 δ 0 , C2 = 0 −δ 0 , C3 = 0 0 δ (5)
2
δ
0 0 δ 0 0 1−δ 2 0 δ 0
The calculations of the elastic properties of the tetragonal structures are continued according to the symmetry of
the tetragonal crystal using linear response theory performing six finite distortions of the lattice; it is completely
described by six independent constants given in Table 3.
Table 3. The calculated elastic constants ( Cij in GPa) for 6–15 mol% doped yttria-stabilized zirconia with linear
response theory.
The mechanical stability criteria for cubic and tetragonal structures are given as
C11 > 0, C33 > 0, C44 > 0, C66 > 0, C11 − C12 > 0, C11 + C33 − 2C13 > 0,
As described below, we have also calculated bulk modulus (B) and shear modulus (G) for cubic and tetragonal
structures from anisotropic elastic constants given as follows:
Cubic structure
BV = BR = (C11 + 2C12 )/3,
Tetragonal structure
1
BV = [2C 11 + C12 ) + C33 + 4C13 ] (7)
9
228
AKGENÇ and ÇAĞIN/Turk J Phys
1
GV = [4C 11 − 2C12 − 4C13 + 2C33 + 12C44 + 6C66 ]
30
[ ]−1
18BV 6 6 3
GR = 15 2 + C + +
(C 11 + C12 )C33 − 2C13 11 − C12 C44 C66
B = (1/2)(BV + BR ),
G = (1/2)(GV + GR ),
where the subscripts V and R stand for Voigt and Reuss forms, respectively.
Young’s modulus E and Poisson’s ratio υ are also obtained using
δ is strain parameter and δ = 0 is referred to the equilibrium volume V 0 at 0 GPa. The strain parameter
is varied from 0.07 to –0.07 in steps of 0.015 and the total energy is calculated at every single increment. A
quadratic behavior in strain energy as a function of strain is used to determine the elastic modulus in Figure 4.
229
AKGENÇ and ÇAĞIN/Turk J Phys
In addition to the structural and energetic stability, mechanical properties are determined in this study.
Extensive systematic studies on the oxygen point defects in cubic ZrO 2 [9] exist, but similar studies are still
lacking for yttria doped zirconia. In this study, we have shown the structural and mechanical properties of 6–15
mol% doped yttria-stabilized zirconia and they are mechanically stable. An understanding of the diffusion of
oxygen vacancies with migration energy barriers is also essential for technological applications [20–23].
Acknowledgments
This work is also supported by Kırklareli University Research Project Unit (BAP) under Project No 125. Partial
support for TÇ is provided by NSF EAGER and NSF IMI project.
References
[1] Fabbri, E.; Pergolesi, D.; Traversa, E. Chem. Soc. Rev. 2010, 39, 4355-4369.
[2] Fabbri, E.; Pergolesi, D.; Traversa, E. Sci. Technol. Adv. Mater. 2010, 11, 054503.
[3] Pergolesi, D.; Fabbri, E.; D’Epifanio, A. Nature Materials 2010, 9, 846-852.
[4] Dyer, C. K. Sci. Am. 1999, 281, 70.
[5] Andersen, N. H.; Clausen, K.; Hackett, M. A.; Hayes, W.; Hutchings M. T.; Macdonald, J. E.; Osborn, R. Physica
1986, 136, 315-317.
[6] Aldert, P.; Traverse, L. P. J. Am. Ceram. Soc. 1985, 68, 34.
[7] Stefanivich, E. V., Shluger, A. L., Catlow, C. R. A. Phys. Rev. B 1994, 49, 11560-11571.
[8] Jomard, G.; Petit, T.; Pasturel, A.; Magaud, L.; Kresse, G.; Hafner, J. Phys. Rev. B 1999, 59, 4044-4052.
[9] Malyi, O. I.; Wu, P.; Kulish, V.V.; Bai, K.; Chen, Z. Solid State Ionics 2012, 212, 117-122.
[10] Cousland, G. P.; Cui, X. Y.; Ringer, S.; Smith, A. E.; Stampfl, A. P. J.; Stampfl, C. M. J. Phys. Chem. Solids
2014, 75, 1252-1264.
[11] Pomfret, M. B.; Chad, S.; Varughese, B.; Walker, R. Anal. Chem. 2005, 77, 1791-1795.
[12] Liang, Z.; Wang, W.; Zhang, M.; Wu, F.; Chen, J. F.; Xue, C.; Zhao, H. Physica B 2017, 511, 10-19.
[13] Xu, Y. N.; Gu, Z. Q.; Ching, W. Y. Phys. Rev. B 1997, 56, 14993-15000.
230
AKGENÇ and ÇAĞIN/Turk J Phys
231