Joohwi Leeab,
Seung-Cheol Leea,
Cheol Seong Hwangb and
Jung-Hae Choi*a
aElectronic Materials Research Center, Korea Institute of Science and Technology, Seoul 136-791, Korea. E-mail: choijh@kist.re.kr; Fax: +82 2 958 6658; Tel: +82 2 958 5488
bWCU Hybrid Materials Program, Department of Materials Science and Engineering, Inter-university Semiconductor Research Center, Seoul National University, Seoul 151-744, Korea
First published on 5th August 2013
Thermodynamic stabilities of various phases in ZnO–SnO2 systems were investigated based on the Gibbs energy obtained from density functional theory (DFT) calculations. The pressure–temperature (p–T) phase diagram was determined; the coexistence of ZnO and SnO2 was the most stable phase in the low temperature region at zero external pressure, while Zn2SnO4 with the inverse spinel structure and ZnSnO3 with the lithium niobate structure were stable at the high temperature and high pressure region. Various octahedral configurations of the inverse spinel structures of Zn2SnO4 were considered. The calculated results showed feasible agreement with experimental data on the phase stability and explained well the experimental observation of the mixed state of Zn2SnO4, ZnO and SnO2 at mid-range temperatures and pressures. Considering the atomic structures, bulk moduli and thermodynamic stabilities, the local density approximation calculations were found to describe experimental observations more accurately than the generalized gradient approximation calculations. The phase transitions in the ZnO–SnO2 system were found to be dominated by the changes in both the Zn–O bond length and the coordination number of Zn, rather than changes in the bond length of Sn–O and the coordination number of Sn.
Making AOS with high mobility basically depends on several principles; firstly, it is necessary to mix component oxides whose crystal structures are quite distinctive. Secondly, the use of as many component oxides as possible, such as in InGaZnO4, is desirable. Thirdly, s orbitals originating from metal ions must be the main components of the lower edge (or mobility edge) of the conduction band, considering their non-directionality and spatial delocalization. Therefore, zinc tin oxide is a feasible candidate for high performance TFTs. However, there could be several complications in the thermodynamic stability of the various stable (metastable) phases in the ZnO–SnO2 system.
An oxide containing Zn and/or Sn can exist as a binary oxide and as a ternary oxide, such as Zn2SnO4 (ZTO214) and ZnSnO3 (ZTO113). The inverse spinel phase of ZTO214 (ref. 12–15) was reported to exist as a single phase at 1275 °C.16,17 Polycrystalline ZTO214 thin films were also reported, which were deposited by sputtering onto glass substrates.12,13 Additionally, two rhombohedral phases of ZTO113 were observed: a lithium niobate phase (LN) and an ilmenite phase (IL). LN-ZTO113 was synthesized from a 1:1 mixture of ZnO and SnO2 by quenching from a pressure greater than 7 GPa and from a high temperature of 1000 °C.18 Meanwhile, IL-ZTO113 was synthesized by the ion-exchange reaction of Li2SnO3, which has the same structure as IL-ZTO113.19 Zinc tin oxide has various phases and shows a tendency to form an amorphous phase.3,4,7,8 Therefore, understanding the crystalline phases and their thermodynamic stability is a prerequisite to figure out the amorphous structure and nanostructure.
Previous ab initio calculations reported that pressurization changes the thermodynamic stability from the coexistence of ZnO and SnO2 to the coexistence of ZTO214 and SnO2.20–22 These calculations only investigated the pressure effect at 0 K based on the enthalpy, and the reported transition pressure varies from 0 to 12 GPa. However, temperature is in fact a much more important thermodynamic variable for determining the stability of materials and for understanding their properties under practical synthesis and working conditions. Furthermore, the space group of the inverse spinel structure varies according to the occupancy of the octahedral sites by Zn and Sn at the unit cell scale when an inverse spinel is assumed to be ordered.15 Seko et al. reported that the most stable structure of ZTO214 belongs to the space group P4122,15 while Gracia et al. calculated the stability of the ZTO214 phase that belongs to the space group Imma.22 Nakayama et al. reported the most stable structure of ZTO214 from their calculations but did not specify the space group.20 Additionally, the different exchange-correlation functionals (xc-functionals) used in the previous calculations led to different optimized structures and energetics, which could affect the thermodynamic stability.
In this paper, ab initio calculations were performed for various zinc tin oxide phases using two xc-functionals. The ab initio phonon calculations were carried out to obtain the enthalpy and Gibbs energy, which were used to investigate the effects of both pressure and temperature on the thermodynamic stabilities. The pressure–temperature (p–T) phase diagram was obtained for the composition with an atomic ratio of Zn to Sn to O of 1:1:3 and compared with previous synthesis or existence conditions (stability range) of various zinc tin oxide phases. The composition of Zn to Sn to O of 2:1:4 was also considered. For these two compositions, several space groups of the ordered and disordered inverse spinel structure were examined. Finally, the changes in the bond length of cation–oxygen and coordination numbers (CN) of the cations were discussed.
Fig. 1 Schematic illustrations of the Gibbs energy as a function of the cell volume. The energy can be extracted using (a) E + pV to consider the pressure effect (T = 0) and (b) E + Fvib to consider the temperature effect (p = 0). |
The thermodynamic stability under external pressure can be determined by the enthalpy, H(p), obtained from the energy–volume curve (E–V curve) calculated from ab initio calculations at 0 K. Fig. 1a illustrates how to add the pV term to the total energy of the electronic structure for a given volume at 0 K, E(V),23 as shown in eqn (1), to consider the external hydrostatic pressure.
H(p) = E(V) + pV | (1) |
The solid circles indicate the E–V curve at p = 0 at 0 K, which are the total energies of the electronic structure for a given volume. The open circles indicate the curves under external p. Note that a positive pressure denotes hydrostatic compression. The minimum thermodynamic values for the given pressures were found and linked by a dashed line, which indicates that the higher the external compressive pressure, the smaller the equilibrium volume.
The temperature effects on the thermodynamic stability can be included by including phonon calculations in the ab initio study. For any point on the E–V curve, the phonon frequency can be obtained using the linear response method based on density functional perturbation theory (DFPT).24 Then, the Helmholtz phonon vibrational energy, Fvib, can be calculated for a given volume and temperature using the quasi-harmonic approximation, as shown in eqn (2),
(2) |
Fig. 1b shows how to calculate the temperature effect shown in eqn (2). The curve composed of solid circles is identical to that in Fig. 1a, whereas the curve composed of solid rhombi was obtained by adding the first term of eqn (2) to account for the temperature-independent vibrational property. The other E–V curves, represented by open rhombi, were obtained for various temperatures, and the minimum energy values for the given temperatures were found and linked by a dashed line. It is noted that the magnitude of the Fvib(T) term increases as the temperature increases, and that this increasing tendency is more significant at larger volumes. These behaviors result in an increase in the equilibrium volume, which describes the thermal expansion.
In addition, for the disordered inverse spinel, the configurational energy is given by eqn (3),15
(3) |
By combining eqn (1)–(3), the Gibbs energy is defined as the minimum value with respect to the volume for a given p and T, as in eqn (4).25,26
G(p, T) = minV[E(V) + pV + Fvib(V, T) + Fconf(T)] | (4) |
The points on the E–V curve were fitted to the Birch–Murnaghan equation of state (BM-EOS)29 before determining the minimum value. The Gibbs energy can be used to determine the phase stability and phase transition. In this paper, the term “transition” is used to describe the onset of changing thermodynamic stability, which is the prerequisite for the real occurrence of a phase transformation.
To obtain the E–V curves, the energies at 15 different volumes under external pressure were calculated for all phases considered. The calculations were performed in two steps, as shown schematically in Fig. 2a. In the first step, full relaxation of the cell parameters and internal atomic coordinates was performed with an applied external p, ranging from −14 GPa to 14 GPa with an interval of 2 GPa. In this step, the reduced or expanded volumes of the unit cells were obtained, which are indicated as the open circles with x-marks in Fig. 2a. Electronic iterations were performed until the energy difference was less than 10−6 eV, and atomic relaxations were performed until the maximum force between atoms was less than 10−3 eV Å−1. All the phases considered in this study show an anisotropic change in their cell parameters under hydrostatic pressure due to their non-cubic symmetry. In the second step, static calculations were performed without changing the volume and shape of the unit cell and the atomic coordinates. The external pressure was eliminated in this step. Then, the total energy of the electronic structure for a given V, E(V), was obtained, which is shown by the solid circles in Fig. 2a. Finally, the E–V curve with 15 points was fitted to the BM-EOS to calculate the optimized unit cell volume and the bulk modulus of each phase at 0 K.
Fig. 2 (a) Schematic illustration of the E–V curve. The volume at step 1 is maintained at step 2. (b) The summation of energy and volume were performed at the same hydrostatic pressure. The positive and negative values correspond to the hydrostatic compressive and tensile stresses, respectively. |
The coexistence of two phases under the same external pressure was obtained by the summation of the E–V curves of the two phases, as shown schematically in Fig. 2b. The x- and y-axes represent the sum of the unit cell volumes and the sum of the electronic structure total energies of the two phases, respectively. The summation of the unit cell volume and the total energy was performed at each external hydrostatic pressure. For example, the energy of the coexistence of phases A and B at pressure p2 is the sum of the energies of phases A and B at pressure p2.
The linear response method based on DFPT was used to calculate the force constants of the supercells.24 Iterations of the electronic steps were performed until the energy differences were less than 10−8 eV. The Helmholtz vibrational energies were used to confirm the convergence of supercells and k-points, which are less than 1 meV f.u.−1. The phonon vibrational frequency and phonon density of states were calculated from the force constants using Phonopy.26,27 Table 1 summarizes the space groups, number of atoms in the unit cells, supercell sizes for the phonon calculations and k-points mesh for the unit cells and supercells.
Composition | Space group (abbreviation) | Unit cell | Number of atoms in unit cell (Zn/Sn/O) | k-Points for unit cell (for supercell) | Supercell size for phonon calculation |
---|---|---|---|---|---|
a Γ-centered.b Monkhorst–Pack.28c This cell is slightly deviated from the primitive cell of the normal spinel, which belongs to the space group Fdm. | |||||
ZnO | P63mc | Hexagonal (primitive) | 2/0/2 | 7 × 7 × 5 − Γa (3 × 3 × 3 − Γ) | 3 × 3 × 2 |
SnO2 | P42/mnm | Tetragonal (primitive) | 0/2/4 | 4 × 4 × 6 − MPb (2 × 2 × 2 − MP) | 2 × 2 × 3 |
ZnSnO3 | R3c (LN) | Hexagonal | 6/6/18 | 5 × 5 × 2 − Γ (3 × 3 × 2 − Γ) | 2 × 2 × 1 |
R (IL) | Hexagonal | 6/6/18 | 5 × 5 × 2 − Γ (3 × 3 × 2 − Γ) | 2 × 2 × 1 | |
Zn2SnO4 | P4122 (tetra) | Tetragonal (primitive) | 8/4/16 | 4 × 2 × 4 − MP (2 × 2 × 2 − MP) | 2 × 1 × 2 |
Imma (ortho) | FCC-likec (primitive) | 4/2/8 | 5 × 5 × 5 − Γ (3 × 3 × 3 − Γ) | 2 × 2 × 2 | |
P1 (mono) | Monoclinic | 16/8/32 | 2 × 2 × 2 − MP (2 × 2 × 2 − MP) | 1 × 1 × 1 |
The Gibbs energy was determined in the range of 0–14 GPa with an interval of 0.1 GPa and in the range of 0–1600 K with an interval of 5 K. Then, the p–T phase diagram was constructed.
For ZTO113, two phases were examined: the lithium niobate phase (R3c, space group 161; LN-ZTO113) and the ilmenite phase (R, space group 148; IL-ZTO113). They have rhombohedral primitive cells;38 however, hexagonal unit cells were used in the ab initio calculations for convenience. In the hexagonal cell, Zn, Sn and a vacancy occupy the octahedral sites in distinctive ways. In IL-ZTO113, a plane normal to the 〈0001〉 direction is occupied by one cation species or a vacancy, whereas it is occupied by the mixture of Zn, Sn and vacancies in LN-ZTO113.
The stable phase of ZTO214 was reported to be an inverse spinel.12–15 If all of the octahedral sites are treated as symmetrically identical due to cation disordering, it is called a disordered inverse spinel ZTO214 (DIS-ZTO214), which belongs to the space group Fdm as the normal spinel.15 However, on the atomic scale, the occupancy of octahedral sites by two different metal elements might lower the symmetry of the structure, which is called an ordered inverse spinel ZTO214 (OIS-ZTO214).
For OIS-ZTO214, three configurations were investigated: a tetragonal phase (P4122, space group 91; tetra-ZTO214),15 an orthorhombic phase (Imma, space group 74; ortho-ZTO214)22 and a monoclinic phase (P1, space group 1; mono-ZTO214). Other unit cells, belonging to the space groups Pm2 and Cm, were found to change into ortho-ZTO214 via relaxation of the shape and volume of the cells. Meanwhile, the volume, energy, and the vibrational properties of DIS-ZTO214 were obtained by averaging over eight different OIS-ZTO214 cells. Six of them were constructed by generating random numbers to have different fillings of the octahedral sites by Zn and Sn, where each individual structure was a mono-ZTO214. Tetra-ZTO214 and ortho-ZTO214 were also included as samples to obtain the average properties of the DIS-ZTO214. Note that tetra-ZTO214 was calculated to have the lowest energy,15 and this was also confirmed by this calculation. The primitive cell of each phase was defined in comparison with the conventional unit cell of the cubic spinel, as shown in Fig. 3.
Fig. 3 Unit cells of various inverse spinel structures of ZTO214. u, v and w are the lattice vectors of the conventional unit cell of the normal spinel structure. Dashed vectors of a1, a2, and a3 indicate the unit vectors of the primitive cell for (a) DIS-ZTO214 and several OIS-ZTO214 structures, which are (b) tetra-ZTO214 (P4122), (c) ortho-ZTO214 (Imma) and (d) mono-ZTO214 (P1), respectively. |
Tables 2 and 3 show the volume of the unit cell and the bulk modulus for each phase from the LDA and GGA calculations, respectively. They also show the comparison with experimental values and some previous calculation results. ZnO in the wurtzite structure and SnO2 in the rutile structure were also calculated as references. The volumes and bulk moduli were obtained from the BM-EOSs.29 The difference between the energy obtained using the relaxed volume and that obtained using the volume fitted to the EOS was negligible, below 1 meV f.u.−1. The calculated volumes of the optimized cells agree well with the experimental data, to within an error of 5%. The volumes obtained from the LDA calculations are slightly smaller and those obtained from the GGA calculations are slightly larger than the experimental values, as observed for many other materials.25 The volumes from the LDA calculations were closer to the experimental data than those from the GGA calculations.
Composition | Space group (abbreviation) | Volume (Å3 f.u.−1) | |||||
---|---|---|---|---|---|---|---|
This work | References | ||||||
LDAa | GGAb | Experimental | LDA | GGA | B3LYPd | ||
a LDA parameterized by Ceperley and Alder.33b GGA parameterized by Perdew et al.34c LDA parameterized by Vosko et al.36d B3LYP parameterized by Lee et al.37e Disordered inverse spinel.f Ordered inverse spinel. | |||||||
ZnO | P63mc | 22.86 | 24.88 | 23.80 [ref. 45] | 22.87a [ref. 23] | 24.83 [ref. 23] | 24.38 [ref. 22] |
SnO2 | P42/mnm | 35.78 | 37.80 | 35.75 [ref. 40] | 34.94c [ref. 39] | 35.50 [ref. 44] | |
ZnSnO3 | R3c (LN) | 54.99 | 58.72 | 55.95 [ref. 18] | 58.60a [ref. 42] | 60.08 [ref. 41] | |
54.15c [ref. 43] | |||||||
R (IL) | 55.50 | 59.68 | 56.78 [ref. 19] | 56.29c [ref. 43] | 60.82 [ref. 41] | ||
Zn2SnO4 | Fdm (DIS)e | 79.15 | 85.34 | 81.11 [ref. 14 and 46] | |||
P4122 (tetra)f | 78.89 | 84.96 | |||||
Imma (ortho)f | 78.91 | 85.02 | 83.74 [ref. 22] | ||||
P1 (mono)f | 79.12 | 85.39 |
Composition | Space group (abbreviation) | Bulk modulus (GPa) | |||||
---|---|---|---|---|---|---|---|
This work | References | ||||||
LDAa | GGAb | Experimental | LDA | GGA | B3LYPd | ||
a LDA parameterized by Ceperley and Alder.33b GGA parameterized by Perdew et al.34c LDA parameterized by Vosko et al.36d B3LYP parameterized by Lee et al.37e Disordered inverse spinel.f Ordered inverse spinel.g No method was reported. | |||||||
ZnO | P63mc | 161 | 128 | 183 [ref. 45] | 162a [ref. 23] | 134 [ref. 23] | 114 [ref. 22] |
SnO2 | P42/mnm | 202 | 173 | 205 [ref. 40] | 245c [ref. 39] | 221 [ref. 44] | |
ZnSnO3 | R3c (LN) | 195 | 167 | 161 [ref. 41] | |||
R (IL) | 173 | 137 | 141 [ref. 41] | ||||
Zn2SnO4 | Fdm (DIS)e | 169 | 135 | 169g [ref. 51] | |||
P4122 (tetra)f | 173 | 142 | |||||
Imma (ortho)f | 170 | 139 | 186 [ref. 22] | ||||
P1 (mono)f | 167 | 133 |
To the best of the authors' knowledge, the bulk moduli of the ternary oxides of ZTO have scarcely been reported experimentally, making the verification of the calculation results by comparing them with experimental results improbable. In order to verify the validity of the calculation of the bulk moduli, therefore, the calculated bulk moduli of the binary oxides, ZnO and SnO2, were compared with experimental data. The bulk moduli from the LDA calculations were close to the experimental values, however, those from the GGA calculations were quite a bit lower than the experimental values. The bulk moduli of the three OIS-ZTO214 phases and DIS-ZTO214 are similar regardless of the space group, i.e., regardless of the configuration of the octahedral sites within the same exchange-correlation functionals. In contrast, the LN-ZTO113 phase shows a slightly higher bulk modulus than the IL-ZTO113 phase.
Fig. 4 E–V curves obtained from the (a) LDA and (b) GGA calculations. The negative common tangent corresponds to the transition pressure. The arrow in (b) indicates that the curve becomes meaningless for IL-ZTO113 under an external pressure over −10 GPa. |
The general tendencies of the E–V curve and the thermodynamic stability are the same for both exchange-correlation functionals: a transition from the coexistence of ZnO and SnO2 at zero external pressure to the coexistence of tetra-ZTO214 and SnO2, and then to LN-ZTO113 as the external hydrostatic pressure increases. The transition reactions for the ZTO113 composition were defined as follows:
Transition A: ZnO + SnO2 → ½(ZTO214 + SnO2) | (5) |
Transition B: ½(ZTO214 + SnO2) → LN-ZTO113 | (6) |
Transition C: ZnO + SnO2 → LN-ZTO113 | (7) |
The transition reactions for the ZTO214 composition were also defined:
Transition A′: 2ZnO + SnO2 → ZTO214 | (8) |
Transition D: ZTO214 → LN-ZTO113 + ZnO | (9) |
Transition C′: 2ZnO + SnO2 → LN-ZTO113 + ZnO | (10) |
Fig. 5 PhDOSs of (a) ZnO, (b) SnO2, (c) tetra-ZTO214, (d) ortho-ZTO214, (e) LN-ZTO113 and (f) IL-ZTO113 at −14, 0 and 14 GPa. The arrow in (f) indicates the negative phonon frequency in IL-ZTO113 at −14 GPa. |
Fig. 6 Gibbs energies of various phases as a function of temperature for the ZTO113 composition using the LDA calculations when the external pressure is (a) p = 0, (b) p = 4 GPa and (c) p = 8 GPa. |
Fig. 6b and c show the pressure effects on the variation of Gibbs energy as a function of temperature. Fig. 6b at 4 GPa shows that the coexistence of tetra-ZTO214 and SnO2 is the most stable below 1035 K, while DIS-ZTO214 replaces tetra-ZTO214 above that temperature. Note that when the external pressure is 4 GPa, the Gibbs energy of LN-ZTO113 drastically decreases, and becomes comparable to the energies of the other phases. In contrast, at 8 GPa in Fig. 6c, LN-ZTO113 is the most stable phase for the whole temperature range, which agrees well with the experimental observation that LN-ZTO113 is the high-pressure phase.18
Fig. 7 p–T phase diagrams for (a) the ZTO113 composition treating ZTO214 as tetra-ZTO214, (b) the ZTO113 composition treating ZTO214 as ortho-ZTO214, (c) the ZTO113 composition considering DIS-ZTO214 instead of OIS-ZTO214s, (d) the ZTO214 composition treating ZTO214 as tetra-ZTO214, and (e) the ZTO214 composition considering DIS-ZTO214 instead of OIS-ZTO214s. |
Fig. 7b shows the p–T phase diagram of the ZTO113 composition when ZTO214 was assumed to be ortho-ZTO214, which is another OIS-ZTO214. As Fig. 6 shows, the Gibbs energy of ortho-ZTO214 was higher than that of tetra-ZTO214 by approximately 0.04 eV f.u.−1. As a result, the coexistence region of ZTO214 and SnO2 does not appear and a phase boundary between the simultaneous presence of ZnO and SnO2 phases and the LN-ZTO113 phase appears (transition C).
Fig. 7c shows the p–T phase diagram of the ZTO113 composition assuming that ZTO214 is DIS-ZTO214. The coexistence region of ZTO214 and SnO2 does not appear at low temperatures due to the high energy state of DIS-ZTO214. However, the additional configurational energy enhances transition A as the temperature increases, therefore, the temperature of transition A at zero external pressure is lower than that of tetra-ZTO214 by more than 300 K. This agrees with the previous simulation result which reported that DIS-ZTO214 is more stable than OIS-ZTO214 at high temperatures over ∼1000 K.15 The phase boundary between the coexistence of ZnO and SnO2 phase and the LN-ZTO113 phase appears as the pressure increases to 5.4 GPa at 0 K.
Fig. 7d shows the p–T phase diagram for the ZTO214 composition when ZTO214 was assumed to be tetra-ZTO214. The thermodynamically stable regions for the coexistence of ZnO and SnO2, tetra-ZTO214, and the coexistence of LN-ZTO113 and ZnO phases are presented. The general tendency of Fig. 7d is similar to that of Fig. 7a, except for the much higher transition pressure required for LN-ZTO113 phase to appear, suggesting that the synthesis pressure should be increased at this ZTO214 composition (transition D compared to transition B). The single phase region of ZTO214 reported in previous phase diagrams16,17 is shown as an inverted triangle, which agrees with this p–T phase diagram.
Fig. 7e shows the p–T phase diagram for the ZTO214 composition when ZTO214 was assumed to be DIS-ZTO214. In the low temperature region, the minimum pressure for LN-ZTO113 to be stable decreased obviously compared to Fig. 7d (from transition D to transition C′). However, due to the stability of DIS-ZTO214 at high temperature, the transition pressure for LN-ZTO113 drastically increases as the temperature increases, which is a remarkable contrast to transition B in Fig. 7c. The transition pressures determined from our calculation results are summarized in comparison with the previous reports in Table 4.
Thermodynamic variable | Ref. | Space group of ZTO214 | xc-functional | Minimum pressure for ZTO214a (GPa) | Minimum pressure for LN-ZTO113b (GPa) |
---|---|---|---|---|---|
a Transition A.b Transtition B or C.c LDA parameterized by Ceperley and Alder.33d GGA parameterized by Perdew et al.34e B3LYP parameterized by Lee et al.37f NA: not available.g The authors of this reference reported that they calculated the most stable structure of ZTO214 but did not specify the space group. | |||||
G | This work | P4122 | LDAc | 2.3 (0 K), 0.6 (1273 K) | 7.3 (0 K), 5.1 (1273 K) |
GGAd | 2.4 (0 K), 0 (1273 K) | 11.7 (0 K), 9.9 (1273 K) | |||
Fdm | LDAc | No transition (0 K), 0 (1273K) | 5.4 (0 K), 6.2 (1273K) | ||
H | This work | P4122 | LDAc | 3.1 (0 K) | 8.5 (0 K) |
GGAd | 3.2 (0 K) | 12.3 (0 K) | |||
H | 20 | NAf,g | GGA | ∼0 (0 K) | 0.8 (0 K) |
H | 21 | NAf | GGA | 5.9 (0 K) | 7.1 (0 K) |
H | 22 | Imma | B3LYPe | 12c (0 K) |
Our calculations describe well the experimental observations; at zero-pressure and low temperature, the simultaneous presence of ZnO and SnO2 is the most stable irrespective of the overall composition and atomic configuration in ZTO214. Experimentally, the mixture films of ZnO and SnO2 were found in an amorphous state below the crystallization temperature of ZTO214, which was reported to be approximately 900 K.7,12,19,47 In other words, they do not form any ternary crystalline compounds in the usual film growth temperature range (<900 K) according to thermodynamic considerations. When the amorphous mixture film was heated up to ∼900 K, they seemed to be separated into the ZnO and SnO2 binary oxides by atomic force microscopy and energy dispersive X-ray spectroscopy.7 In contrast, in the high temperature region, ZTO214 is stable regardless of the disordering of ZTO214. For the ZTO113 composition, Inaguma et al. reported that LN-ZTO113 decomposed into mixed phases comprised of ZTO214 and SnO2 above 1073 K.18 The decomposition of IL-ZTO113 into the coexistence of ZTO214 and SnO2 was also observed for a film above 873 K (ref. 47) and in the bulk above 1273 K.18 Additionally, for the ZTO214 composition, Wang et al. reported that the simultaneously present ZnO and SnO2 phases gradually changed into ZTO214 at approximately 1000 K and finally disappeared at 1200 K.49,50 At the high temperature and high pressure region, the synthesis condition of LN-ZTO113 (ref. 18) was represented as a triangle in Fig. 7a–c. The results also agree with the calculation results, regardless of the local atomic configuration of ZTO214.
There are certain discrepancies between the theoretical and experimental results for the phase transition temperatures and pressures, which could be attributed to several factors, although the general trends are well matched. One of the reasons is that the thermodynamic calculations using the ab initio method cannot account for the occurrence of a transformation route due to the kinetic activation barrier. In addition, our assumption to describe the cation disordering for DIS-ZTO214 includes the averaging effect with a limited number of structures, which may bring about some deviation of transition temperature and pressure. As shown in Fig. 7, the thermodynamically stable region can change with the atomic distribution such that ZTO214 can be the stable phase for mid-range pressures depending on the distribution of Zn and Sn. In fact, it seems plausible that actually ZTO214 has a mixture of a disordered inverse spinel structure and several ordered inverse spinel structures.
Another factor is that the experimental data were almost always obtained from thin films, whose structures were amorphous in most cases. During crystallization, they are under biaxial compressive or tensile stress depending on the mismatch of both the lattice parameters and the thermal expansion with the substrates. The compressive stress on the order of GPa (about 104 atm) is practically quite high. However, such a level of stress is often observed in thin film systems. For example, 6 GPa corresponds to approximately 2% biaxial strain in LN-ZTO113 thin films, which is not an uncommon level of strain in several thin film systems. The diverse experimental observations and the coexistence of multi-phases can be explained by the various configurations of ZTO214 and kinetics. In fact, Hoel reported the simultaneous presence of ZTO214, ZnO and SnO2 phases in a polycrystalline film when the film was formed at a high temperature of over 927 K.47,48
Table 4 shows the transition pressures from our calculation results in comparison with the previous reports. For the ZTO113 composition, the sequence of transition A and B along the pressure increase was the same as the previous reports which only considered the enthalpy.20,21 Nakayama et al. calculated that the transition pressures of transitions A and B are 0 and 0.8 GPa, respectively.20 However, they mentioned three plausible reasons for their underestimation of the transition pressures: the softness of the GGA functional, the exclusion of the temperature effect, and the method of selecting the occupancy of Zn and Sn in the unit cell. It should be pointed out that this study took these three drawbacks into account by testing the LDA functional as well as the GGA functional, by using the Gibbs energy to estimate the thermodynamic stability, and by investigating several atomic configurations in both ZTO113 and ZTO214.
The GGA results in Table 4 overestimated the transition pressure and temperature and disagreed with this experimental observation. The use of the exchange-correlation functional of the DFT+U method35 was also attempted. However, it was not helpful for the thermodynamic calculation of zinc tin oxide systems; when the LDA+U scheme was used, all the volumes of the zinc tin oxide phases decreased slightly. This means that the differences between the calculated and experimental values of the bond lengths and the cell volumes increased. Because the bond lengths are related to the phonon vibrations, the LDA+U scheme was not appropriate to use. Conversely, when the GGA+U scheme was adopted, the calculated volume became closer to the experimental value. However, the failure to describe the tensile region of IL-ZTO113 and the stable region of LN-ZTO113 also remained in the GGA+U calculation.
Fig. 8 Bond lengths of (a) Zn–O and (b) Sn–O as a function of the external pressure at 0 K. Black, blue and red symbols indicate the bonds in the binary phases, tetra-ZTO214 and LN-ZTO113 phases, respectively. The squares and stars in (a) represent the Zn–O bonds where the CNs of Zn are 4 and 6, respectively. The hexagons in (b) represent the Sn–O bonds where the CN of Sn is 6. |
Composition | Space group (abbreviation) | xc-functional | Bond length of Zn–O (Å) | CN of Zn | Bond length of Sn–O (Å) | CN of Sn |
---|---|---|---|---|---|---|
a LDA parameterized by Ceperley and Alder.33b GGA parameterized by Perdew et al.34 | ||||||
ZnO | P63mc | LDAa | 4 × 1.951 (±0.003) | 4 | — | — |
GGAb | 4 × 2.008 (±0.004) | — | ||||
SnO2 | P42/mnm | LDA | — | — | 6 × 2.052 (±0.006) | 6 |
GGA | — | 6 × 2.093 (±0.001) | ||||
ZnSnO3 | R3c (LN) | LDA | 3 × 1.97, 3 × 2.29 | 6 | 3 × 2.03, 3 × 2.09 | 6 |
GGA | 3 × 2.03, 3 × 2.35 | 3 × 2.07, 3 × 2.13 | ||||
Zn2SnO4 | P4122 (tetra) | LDA | 2 × 1.95, 2 × 1.99 | 4 | 2 × 2.03, 2 × 2.07, 2 × 2.10 | 6 |
2 × 2.03, 2 × 2.08, 2 × 2.19 | 6 | |||||
GGA | 2 × 2.01, 2 × 2.05 | 4 | 2 × 2.06 2 × 2.11, 2 × 2.14 | 6 | ||
2 × 2.08, 2 × 2.15, 2 × 2.25 | 6 |
In Fig. 8a, all Zn–O bond lengths decrease almost linearly as the pressure increases, and the longer bonds decrease slightly more in length than the shorter bonds. In the low pressure region, all Zn atoms in the stable phase, ZnO, are located at tetrahedral sites, while in the high pressure region, all Zn atoms in the stable LN-ZTO113 are located at octahedral sites. At mid-range pressure, in contrast, half of the Zn atoms occupy the tetrahedral sites and the other half occupy the octahedral sites in tetra-ZTO214. This increasing tendency of the CN of Zn with increasing pressure is the same as for ZnO. In fact, the CN of Zn in ZnO increases from four to six with the change in the stable phase from the wurtzite to the rock salt structure at 9 GPa.23
The Sn–O bond lengths in Fig. 8b also decrease as the external pressure increases, however, they show a very different tendency compared to the Zn–O bonds. All Sn atoms are located at the octahedral sites in all the phases considered, hence their CN is kept at six. In addition, the discrepancy in the bond length of Sn–O in each phase is remarkably small compared to the case of Zn. It is also noted that the rutile structure SnO2 with the CN of 6 is stable up to 21 GPa,40 which is much higher than the transition pressure of ZnO. This means that Zn–O bonds are more sensitive to the pressure change compared to Sn–O bonds, and this can be inferred from the lower bulk modulus of ZnO than of SnO2 in Table 3. A similar tendency was also observed in ZTO as shown in Fig. 8a and b. Therefore, the phase stabilities and phase transitions in the ZTO system are dominated by changes in both the bond length of Zn–O and the CN of Zn, rather than changes in the bond length of Sn–O and the CN of Sn.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c3tc30960f |
This journal is © The Royal Society of Chemistry 2013 |