Advances in Condensed Matter Physics
Volume 2016 (2016), Article ID 9435387, 7 pages
http://dx.doi.org/10.1155/2016/9435387

Research Article

Theoretical Investigation on Structural and Electronic Properties of InN Growth on Ce-Stabilized Zirconia (111) Substrates

Yao Guo,¹ Taixuan Jia,¹ Chengbo Li,² Yongsheng Niu,¹ Shaogang Hou,¹ and Shuanjiang Liu²

¹Department of Chemistry and Environment Engineering, Anyang Institute of Technology, Anyang 455000, China
²Department of Mathematics and Physics, Anyang Institute of Technology, Anyang 455000, China

Received 21 November 2015; Accepted 23 February 2016

Academic Editor: Da-Ren Hang

Copyright © 2016 Yao Guo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The structural and electronic properties of InN on Ce-stabilized zirconia (CeSZ) (111) substrates are investigated using first-principles calculations based on density functional theory with GGA + method. Surface energy calculations indicate that the structure of Ce-segregated surface is more energetically stable than that of Ce-segregation-free surface. Adsorption energies of indium and nitrogen atoms on both Ce-segregated and Ce-segregation-free CeSZ (111) surfaces at the initial growth stage have been studied. The results suggest that the first layer of InN films consists of a nitrogen layer, which leads to epitaxial relationships between InN (0001) // CeSZ (111) and InN // CeSZ . In addition, density of states (DOS) analysis revealed that the hybridization effect plays a crucial role in determining the interface structure for the growth of InN on CeSZ (111) surfaces. Furthermore, adsorption energies of indium atoms on the nitrogen layer have also been evaluated in order to investigate the lattice polarity determination for InN films. It was found that an indium atom preferentially adsorbs at the center of three nitrogen atoms stacked on the CeSZ substrate, which results in the formation of In-polarity InN.

1. Introduction

Indium nitride (InN) is a promising material for high-speed electronic devices owing to its highest electron mobility [1–3] and saturation drift velocity [4] among the group III nitride. Growth of high-quality ultrathin InN films is necessary for the fabrication of InN-based devices since the device performance is often restricted by the threading dislocations and stacking faults at the interface [5]. Foreign substrates such as sapphire or GaN [6, 7] are usually employed for epitaxial growth of InN because of the lack of bulk InN crystals. However, the large lattice mismatches between InN and these substrates cause the formation of a high density of threading dislocations. In order to solve the problem above, development of lattice-matched substrates for InN is required. Among the various substrates, cubic phase zirconia substrates serve as a promising candidate for InN growth because they share a threefold rotational symmetry with the InN (0001) surface and have small lattice mismatch [8]. For example, the use of yttria-stabilized zirconia (YSZ) substrates provides reduced defect density and abrupt InN/YSZ interfaces [9]. In fact, the fluorite-type zirconia could be stabilized by the addition of different dopants. The stabilized cubic zirconia has its own characteristic structures which are dopant-independent [10]. The introduction of new dopants into zirconia that replaces yttrium is a straightforward strategy because the interaction between the dopant and the InN film is quite different.

Ceria-stabilized zirconia (CeSZ) is known as heterogeneous catalysis in the area of automotive pollution control [11]. On the other side, CeSZ is an ideal substrate owing to the small lattice mismatch between the InN film and the substrate. In fact, it has been experimentally determined that Ce atoms stabilize the cubic structure of zirconia and the Ce dopant segregation on the CeSZ (111) surface could be obtained due to the cation redistribution at high temperatures [12]. However, the feasibility of epitaxial growth of InN films on CeSZ (111) has not yet been investigated to date. Hence, a clear understanding of the InN/CeSZ interface microstructure is especially meaningful. To investigate the influence of the interfacial structure on electronic properties and the possible epitaxial relationships in the interface, the mechanism of growth of InN layers on CeSZ (111) substrates is discussed on the basis of first-principles calculations.

2. Computational Method

The theoretical calculations have been carried out in the framework of the ab initio density functional theory using the software MedeA VASP package [13]. The electron exchange and correlation were treated within the generalized gradient approximation (GGA-PBE) [14]. The wavefunctions were expanded in a plane-wave basis set with a cutoff energy of 800 eV. The electron-ion interactions were described by the full-potential PAW method [15]. The electron configurations of the pseudopotentials for In, N, Zr, Ce, and O are , , , , and , respectively. All calculations involving Ce were performed at the GGA + level with a Hubbard parameter  eV for the Ce orbitals [16]. The oxygen-terminated surface was employed because it was more stable than the Zr terminated or O-O terminated [17]. The slab contained twelve atomic layers (four O-Zr-O trilayers) separated from each other by a vacuum layer of 15 Å in the perpendicular direction to eliminate any interactions between the adjacent slabs. In the optimization process, the layers in the bottom O-Zr-O trilayer were fixed at the bulk positions, while the top layers were fully relaxed with convergence criteria of 0.05 eV/Å. The supercell consisted of a 2 × 2 array of CeSZ (111) surface unit cells. The Brillouin-zone integration was calculated with a 2 × 2 × 1 k-point grid, which was generated using the Monkhorst-Pack method [18]. We also performed the test with a denser 4 × 4 × 1 k-point grid and found that the difference in total energy is less than 1 × 10⁻⁴ eV/atom. Therefore, the 2 × 2 × 1 k-point grid is sufficient for the calculation. Similar computational details are described elsewhere [19].

3. Results and Discussion

For necessity of later calculations, the equilibrium lattice of CeSZ was investigated in total energy versus lattice parameter. A stoichiometric CeZrO₄ cubic structure in the 1 × 1 × 1 unit cell is employed for the calculation. Figure 1 shows that the total energy of CeSZ bulk was calculated using GGA + and GGA and plotted as a function of the lattice constant. The effective Hubbard parameter of 3 eV and 5 eV was used in the calculation as the GGA + approach. The  eV calculation corresponds to unmodified GGA. It allows a direct comparison of the two methods. Introducing into the GGA methods results in an increase of lattice constant so that eventually GGA + may predict the experimental lattice constant  Å [20]. Therefore, the value of  eV as suggested by Nolan et al. [16] was employed for the later calculations.

Figure 1: The equilibrium lattice constant of CeSZ calculated using the GGA and GGA + .

Based on the calculated lattice constant, the CeSZ (111) surface is constructed using a three-dimensional periodic surface slab mode for the investigation of adsorption phenomena. Both Ce-segregation-free and Ce-segregated surfaces are considered for the theoretical modeling. The optimized structure of the CeSZ (111) surface is shown in Figure 2. The segregation behavior of CeSZ (111) surface has been found in both theoretical and experimental studies [12]. Although some experimental studies such as X-ray photoelectron spectroscopy (XPS) and low energy ion spectroscopy (LEIS) measurements have found the dopant segregation in zirconia surface [21, 22], the atomic structure of the dopant stabilized zirconia surface is still unclear due to the limit of the experimental techniques at present stage. Therefore, the structure of Ce-segregation surface is assumed in this work. All the Zr atoms in the upper half of the zirconia are replaced with Ce atoms. Similar assumption has been performed in the case of the Y-segregated YSZ surface structure [19].

(a)

(b)

Figure 2: Schematic illustration of the (a) Ce-segregation-free and (b) Ce-segregated CeSZ (111) surface used for calculations.

The stability of the surface is dependent on the surface energy. The surface energy of the Ce-segregation-free and Ce-segregated surfaces was calculated by

where is the energy of the surface slab, is the energy of the bulk unit cell, is the number of atoms in the slab, is the number of atoms in the bulk unit cell, and is the unit area of the surface. According to (1), a positive and low value of corresponds to a stable surface. The calculated surface energies of both Ce-segregation-free and Ce-segregated CeSZ (111) surfaces are summarized in Table 1. Surface energy calculations revealed that the structure with the Ce atoms segregation is more energetically stable than that without Ce atoms segregation. Because the covalent character of Zr-O bond is stronger than that of Ce-O bond [12], the unsaturated CeSZ (111) surface with breaking of the metal-oxygen bonds should be more easily enriched by Ce cations than Zr cations in the top layer.

Table 1: Surface energy (in eV/A²) of Ce-segregation-free and Ce-segregated CeSZ (111) surfaces.

Surface energy

Ce-segregation-free
0.056

Ce-segregated
0.044

Calculations of the adsorption energies of In and N atoms on the two CeSZ (111) surfaces were then performed in order to discuss the growth mechanisms of InN on CeSZ substrates during the initial stage. The adsorption energy is defined as

where is the energy of the substrate with an adatom adsorbed onto it, is the energy of the substrate without the adatom, and is the energy of an isolated adatom. According to (2), generally has a positive value, and a larger value corresponds to stable adsorption. Systematic probing of all symmetric sites on the surfaces is essential for the determination of the adsorption energy. For the Ce-segregation-free surface, nine highly symmetric adsorption sites on the surface were selected for calculations, as shown in Figure 3(a). For the Ce-segregated surface, seven highly symmetric adsorption sites on the surface were selected for calculations, as shown in Figure 3(b). These adsorption sites included the on-top sites (on top of the Ce, Zr, surface oxygen Ou, and subsurface oxygen Od), the bridge sites (in the middle between the Ce-Ou, Zr-Ou, Ce-Od, Zr-Od, and Ou-Od pairs), and the hollow sites (in the center of the Ce-Ou-Od, Zr-Ou-Od triangles). In the calculation, the surface atoms were kept frozen at the positions obtained from the clean surface optimization, and only the positions of the adatoms were fully optimized until the Hellmann-Feynman forces were below 0.05 eV/Å. The calculated adsorption energies of the In and N atoms on the CeSZ (111) surface are summarized in Tables 2 and 3.

Table 2: Adsorption energy (in eV) of In and N atoms on Ce-segregation-free CeSZ (111) surfaces.

Adatom
Adsorption site

Top Ou
Top Od
Top Ce
Top Zr
Bridge Zr-Ou
Bridge Zr-Od
Bridge Ce-Ou
Bridge Ce-Od
Bridge Ou-Od
Hollow Zr-Ou-Od
Hollow Ce-Ou-Od

In
1.86
1.49
1.17
1.61
1.61
1.41
1.54
1.51
1.82
1.75
1.48

N
3.32
1.49
1.12
1.56
3.78
2.69
3.40
1.79
3.61
3.76
3.37

Table 3: Adsorption energy (in eV) of In and N atoms on Ce-segregated CeSZ (111) surfaces.

Adatom
Adsorption site

Top Ou
Top Od
Top Ce
Bridge Ce-Ou
Bridge Ce-Od
Bridge Ou-Od
Hollow Ce-Ou-Od

In
1.74
1.40
1.26
1.47
1.41
1.70
1.55

N
3.48
1.19
1.07
3.58
1.98
4.10
3.37

(a)

(b)

Figure 3: Schematic illustration of the adsorption sites for (a) Ce-segregation-free and (b) Ce-segregated CeSZ (111) surfaces.

According to the results of adsorption energies, the most stable site for adsorption is a bridge position. It can be seen that the difference in the adsorption energies of In atoms on various sites is less than the difference in the adsorption energies of N atoms on various sites. We assumed that the barrier to diffusion can be obtained by comparing adsorption energies at different high symmetry sites. This result indicates that migration of indium atoms on the surface readily occurs. In contrast, the larger difference in the adsorption energies of N atoms on various sites suggests that N atoms tend to stay at the most stable site. These results indicate that the first layer of InN (0001) should be an N layer, and that the in-plane epitaxial relationship is InN // CeSZ , which gives a small lattice mismatch between InN (0001) and CeSZ (111), as expected.

Although the higher stability of a nitrogen atom compared to that of an indium atom adsorbed on both Ce-segregation-free and Ce-segregated CeSZ (111) surfaces is a common phenomenon, the adsorption energies of the nitrogen atoms are quite different for the different substrates. Thus, to shed light on the interaction mechanism between N atom and CeSZ (111) surface, the density of state (DOS) for nitrogen adsorption at the most favorable site (Bridge Zr-Ou) on the segregation-free surface has been performed, and the results are shown in Figure 4(a). In addition, similar DOS for nitrogen adsorption at the most favorable site (Bridge Ou-Od) on Ce-segregated surface are shown in Figure 4(b). The atomic DOS of the N and its nearest O and cation (Zr or Ce) atoms in the surface region are plotted, respectively. It can be seen that the peaks in the range from −5 to 0 eV belong to orbital; the narrow peak around −12 eV is assigned to the orbital. In all cases there exists strong hybridization among the , , and cation (Zr or Ce) orbitals within the energy range from −5 to 0 eV. This implies that the strong hybridization is the main mechanism of N atom interaction with CeSZ (111) surface and determines the value of the adsorption energy for N atom. The observation is in agreement with our previous research for YSZ (111) surface [19].

(a)

(b)

Figure 4: The DOS for the structure of a nitrogen atom on (a) Ce-segregation-free and (b) Ce-segregated CeSZ (111) surfaces.

Furthermore, for a better characterization of the InN/CeSZ interfacial structure, we performed the calculation of the adsorption energies of an indium atom on the first nitrogen layer covered CeSZ (111) surfaces to compare the stability of In-polar InN with that of N-polar InN on Ce-segregation-free and Ce-segregated surfaces, respectively. In fact, the lattice polarity of the InN is determined by the adsorption site for In atom. The In-polarity InN is modeled by the adsorption of an indium atom at the center of three nitrogen atoms (Hollow N site). By comparison, In atom adsorption on the top of a nitrogen atom (Top N site) corresponds to N-polarity InN [19].

The adsorption energies of an In atom at the “Hollow N” and “Top N” sites on both Ce-segregation-free and Ce-segregated CeSZ (111) surfaces are summarized in Table 4. It is found that the adsorption energy for In-polar InN is larger than that for N-polar InN on both two CeSZ surfaces. The growth of In-polar InN is energetically preferred. Moreover, the difference between In- and N-polar InN for Ce-segregated surface is smaller than that for Ce-segregation-free surface. It indicates that the polarity conversion between In- and N-polar InN on Ce-segregated substrate during growth easily occurred due to the low energy barrier. The difference in polarity determination on Ce-segregation-free and Ce-segregated surfaces is generally attributed to the Ce-segregation on the surface. The different -orbital electrons of Zr and Ce will affect the adsorption of a nitrogen atom on CeSZ surface. Moreover, the adsorption of an indium atom on the first nitrogen layer covered CeSZ surface. The lattice polarity determination is dependent on the surface segregated Ce atoms.

Table 4: Adsorption energies (in eV) of an In atom on the nitrogen layer covering CeSZ (111) surfaces.

Surface
Adsorption energy of In atom

In-polarity
N-polarity

Ce-segregated
3.32
3.26

Ce-segregation-free
3.95
3.66

4. Conclusion

In this study, we have investigated the growth mechanism of InN on a CeSZ (111) surface by first-principles calculations. N atoms have larger adsorption energies than In atoms on both Ce-segregation-free and Ce-segregated CeSZ (111) surface. The results imply that the first layer at the initial stage of growth should be a nitrogen layer. The theoretical studies have revealed that the epitaxial relationships of InN (0001) // CeSZ (111) and InN // CeSZ are determined by the adsorption process of nitrogen atoms on the surfaces at the initial stage of the growth. The DOS analysis reveals that there is a clear strong hybridization among the , , and cation (Zr or Ce) orbitals. It is responsible for the main mechanism of N atom interaction with CeSZ (111) surface and determines the interface structure. Moreover, the mechanism of the determination of InN polarity on CeSZ (111) substrates was studied by theoretical calculations. It was found that the growth of In-polar InN is energetically more stable than that of N-polar InN. The difference in polarity determination is determined by different -orbital electrons of Zr and Ce atoms.

Competing Interests

The authors declare that there is no competing interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. U1404217), the Foundation of He’nan Educational Committee (Grant no. 16B430001), and the Program of Department of Science and Technology of He’nan Province (Grant no. 152102210195).

References

S. L. Wang, H. X. Liu, X. Song, Y. L. Guo, and Z. N. Yang, “An analytical model of anisotropic low-field electron mobility in wurtzite indium nitride,” Applied Physics A: Materials Science and Processing, vol. 114, no. 4, pp. 1113–1117, 2014.

K. Okubo, A. Kobayashi, J. Ohta, M. Oshima, and H. Fujioka, “Electron mobility of ultrathin InN on yttria-stabilized zirconia with two-dimensionally grown initial layers,” Applied Physics Letters, vol. 102, no. 2, Article ID 022103, 2013.

X. Wang, S. Liu, N. Ma et al., “High-electron-mobility InN layers grown by boundary-temperature-controlled epitaxy,” Applied Physics Express, vol. 5, no. 1, Article ID 015502, 2012.

W. A. Hadi, P. K. Guram, M. S. Shur, and S. K. O'Leary, “Steady-state and transient electron transport within wurtzite and zinc-blende indium nitride,” Journal of Applied Physics, vol. 113, no. 11, Article ID 113709, 2013.

J. Kioseoglou, E. Kalesaki, G. P. Dimitrakopulos, T. Kehagias, P. Komninou, and T. Karakostas, “Atomistic modeling and HRTEM analysis of misfit dislocations in InN/GaN heterostructures,” Applied Surface Science, vol. 260, pp. 23–28, 2012.

N. Nepal, N. A. Mahadik, L. O. Nyakiti et al., “Epitaxial growth of cubic and hexagonal InN thin films via plasma-assisted atomic layer epitaxy,” Crystal Growth and Design, vol. 13, no. 4, pp. 1485–1490, 2013.

P. Wang, Y. Yuan, C. Zhao et al., “Lattice-polarity-driven epitaxy of hexagonal semiconductor nanowires,” Nano Letters, vol. 16, no. 2, pp. 1328–1334, 2016.

T. Honke, H. Fujioka, J. Ohta, and M. Oshima, “InN epitaxial growths on Yttria stabilized zirconia (111) step substrates,” Journal of Vacuum Science & Technology A, vol. 22, no. 6, pp. 2487–2489, 2004.

M. Oseki, K. Okubo, A. Kobayashi, J. Ohta, and H. Fujioka, “Field-effect transistors based on cubic indium nitride,” Scientific Reports, vol. 4, article 3951, 2014.

M. O. Zacate, L. Minervini, D. J. Bradfield, R. W. Grimes, and K. E. Sickafus, “Defect cluster formation in M₂O₃-doped cubic ZrO₂,” Solid State Ionics, vol. 128, no. 1–4, pp. 243–254, 2000.

D. R. Sellick, A. Aranda, T. García et al., “Influence of the preparation method on the activity of ceria zirconia mixed oxides for naphthalene total oxidation,” Applied Catalysis B: Environmental, vol. 132-133, pp. 98–106, 2013.

R. G. Crespo, N. H. De Leeuw, S. Hamad, and U. V. Waghmare, “Phase separation and surface segregation in ceria-zirconia solid solutions,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 467, no. 2131, pp. 1925–1938, 2011.

X. Wei, C. Dong, Z. Chen, K. Xiao, and X. Li, “Density functional theory study of SO₂-adsorbed Ni(1 1 1) and hydroxylated NiO(1 1 1) surface,” Applied Surface Science, vol. 355, pp. 429–435, 2015.

J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, no. 18, pp. 3865–3868, 1997.

P. E. Blöchl, “Projector augmented-wave method,” Physical Review B, vol. 50, no. 24, pp. 17953–17979, 1994.

M. Nolan, “Modifying ceria (111) with a TiO₂ nanocluster for enhanced reactivity,” The Journal of Chemical Physics, vol. 139, no. 18, Article ID 184710, 2013.

A. Christensen and E. A. Carter, “First-principles study of the surfaces of zirconia,” Physical Review B, vol. 58, no. 12, pp. 8050–8064, 1998.

H. J. Monkhorst and J. D. Pack, “Special points for Brillouin-zone integrations,” Physical Review B, vol. 13, no. 12, pp. 5188–5192, 1976.

Y. Guo, S. Inoue, A. Kobayashi, J. Ohta, and H. Fujioka, “Theoretical investigation of the polarity determination for c-plane InN grown on yttria-stabilized zirconia (111) substrates with yttrium surface segregation,” Applied Physics Express, vol. 6, no. 2, Article ID 021002, 2013.

M. Kuhn, S. R. Bishop, J. L. M. Rupp, and H. L. Tuller, “Structural characterization and oxygen nonstoichiometry of ceria-zirconia solid solutions,” Acta Materialia, vol. 61, no. 11, pp. 4277–4288, 2013.

A. Bernasik, K. Kowalski, and A. Sadowski, “Surface segregation in yttria-stabilized zirconia by means of angle resolved X-ray photoelectron spectroscopy,” Journal of Physics and Chemistry of Solids, vol. 63, no. 2, pp. 233–239, 2002.

M. de Ridder, R. G. van Welzenis, A. W. Denier van der Gon et al., “Subsurface segregation of yttria in yttria stabilized zirconia,” Journal of Applied Physics, vol. 92, no. 6, pp. 3056–3064, 2002.