Advances in Condensed Matter Physics
Volume 2015 (2015), Article ID 120294, 9 pages
http://dx.doi.org/10.1155/2015/120294

Research Article

Assessment of Functionals for First-Principle Studies of the Structural and Electronic Properties of -Bi₂O₃

D. H. Galván,¹ R. Núñez-González,² R. Rangel,³ P. Alemany,4,5 and A. Posada-Amarillas⁶

¹Centro de Nanociencias y Nanotecnología, Universidad Nacional Autónoma de México, Apartado Postal 2681, 22800 Ensenada, BCN, Mexico
²Departamento de Matemáticas, Universidad de Sonora, 83000 Hermosillo, SON, Mexico
³División de Estudios de Posgrado, Facultad de Ingeniería Química, Edificio V1, Ciudad Universitaria, 58060 Morelia, MICH, Mexico
⁴Instituto de Química Teórica i Computacional (IQTCUB), Universitat de Barcelona, Marti i Franqués 1-11, 08028 Barcelona, Spain
⁵Departamento de Química Física, Universitat de Barcelona, Marti i Franqués 1-11, 08028 Barcelona, Spain
⁶Departamento de Investigación en Física, Universidad de Sonora, Apartado Postal 5-088, 83190 Hermosillo, SON, Mexico

Received 19 April 2015; Revised 29 June 2015; Accepted 1 July 2015

Academic Editor: Markus R. Wagner

Copyright © 2015 D. H. Galván 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.

Fully relativistic full-potential density functional calculations with an all-electron linearized augmented plane waves plus local orbitals method were carried out to perform a comparative study on the structural and electronic properties of the cubic oxide -Bi₂O₃ phase, which is considered as one of the most promising materials in a variety of applications including fuel cells, sensors, and catalysts. Three different density functionals were used in our calculations, LDA, the GGA scheme in the parametrization of Perdew, Burke, and Ernzerhof (PBE96), and the hybrid scheme of Perdew-Wang B3PW91. The examined properties include lattice parameter, band structure and density of states, and charge density profiles. For this modification the three functionals reveal the characteristics of a metal and the existence of minigaps at high symmetry points of the band structure when spin-orbit coupling is taken into account. Density of states exhibits hybridization of Bi 6s and O 2p orbitals and the calculated charge density profiles exhibit the ionic character in the chemical bonding of this compound. The B3PW91 hybrid functional provided a better agreement with the experimental result for the lattice parameter, revealing the importance of Hartree-Fock exchange in this compound.

1. Introduction

The VI-bismuth compounds Bi₂O₃, Bi₂S₃, and Bi₂Te₃ comprise a family of technological semiconductors [1] employed in diverse areas of industry. Bi₂O₃ has been used as an effective cocatalyst for oxidation of ammonia since the beginning of the last century; also it has been used as the precursor for the production of Bi₂MoO₆, a catalyst popularly used in the selective oxidation of propylene into acrolein [2]. Galván et al. [3] showed that, using the same compound in CO to CO₂ conversion, a significant improvement in the oxidation process could be obtained, while it is likely to be used as a catalyst in catalytic converters to clean the exhaust discharge. Moreover, metal oxides, bismuth oxides, and rare-earth oxides form the basis for some of the ceramic high -superconductors [4].

Bismuth oxide is a polymorphic material that crystallizes in five modifications [5]. δ-Bi₂O₃ is formed from heating monoclinic α-Bi₂O₃ at 730°C and it is stable up to the melting point of 824°C. Tetragonal (β-Bi₂O₃) and body-centered cubic (γ-Bi₂O₃) modifications have been obtained below 650°C as metastable phases [6]. The hydrothermal method was used to obtain ε-Bi₂O₃ phase [7]. It is currently accepted that the cubic oxide δ-Bi₂O₃ is one of the most promising materials which can be used as solid electrolyte [8] in high technology equipment such as fuel cells, sensors, and membrane devices which utilize the unusual high oxide ion conductivity of this compound, a conductivity which exceeds that of stabilized zirconia [9, 10]. Several structural models for δ-Bi₂O₃ have been proposed [8, 11–13]. According to a model by Sillen [12], six possible oxygen sites of the fluorite unit cell are occupied and two others are empty and ordered in the direction. However, the large magnitude of the conductivity suggests that a high disorder on the oxygen sublattice might exist which led other researchers [13] to propose a model of statistically average occupation of all eight oxygen positions.

A number of theoretical studies have been carried out in order to understand the basic properties of metal oxide systems. For example, using the LMTO method, Medvedeva et al. [14] calculated electronic properties of δ-Bi₂O₃ with several oxygen vacancy configurations to find out the role of Bi-O and Bi-Bi bonding in the stabilization of such derived structures. Several reports have been published concerning (width of the forbidden gap between VB and CB) values for similar oxides. Gubanov and Medvedeva [15] performed an LMTO-ASA calculation, obtaining values ranging from 1.75 to 5.24 eV, respectively, for SnO₂ and ZrO₂. Different phases of Bi₂O₃ like α-Bi₂O₃, β-Bi₂O₃, and δ-Bi₂O₃ were calculated by Evarestov et al. [16] performing a SC-CNDO (Complete Neglect of Differential Overlap) calculation obtaining values ranging from 4.8 to 6.2 eV for the different phases. Reyes [17] using extended Hückel type calculations obtained 3.48 eV in the theoretical study of α-Bi₂O₃. Walsh et al. [18] reported a theoretical and experimental investigation of the electronic properties of δ-Bi₂O₃ using a combination of gradient corrected density functional theory along with X-ray photoemission and O-K shell X-ray absorption and emission spectroscopy. They reported that vacancies along result in the most energetically stable configuration of δ phase.

Despite the existence of an important number of theoretical studies using DFT, a systematic study to determine what is the most suitable exchange-correlation (XC) potential for this type of metal oxides has not yet been reported. The best performance of GGA schemes over the LDA one is usually claimed, attributing discrepancies to the analytic form of functionals. The aim of this work is to make a detailed comparative study of the structural and electronic properties of δ-Bi₂O₃ in order to assess to what extent the density functional differences are important in the qualitative analysis of theoretical results under the DFT scheme. Besides the XC potential analysis for this bismuth oxide phase, the spin-orbit (SO) interaction effect on electronic properties is also analyzed. We explore distinctive minigaps obtained in our calculations of electronic band structure which are essentially attributed to the SO effect. Conventional density functional calculations were implemented using the APW + lo method based WIEN2k program package, paying careful attention to the lattice parameter, band structure, and density of states as well as charge density profiles. In our calculations it is apparent that the SO effect is not large enough to modify the metallic nature of δ-Bi₂O₃; however, the opening of minigaps distributed along the conduction band is due to this effect. Specific calculation details are given in Section 2; results and discussion are shown in Section 3. In Section 4 we give our conclusions.

2. Structural and Computational Details

The structure of δ-Bi₂O₃ used in the present work is that reported by Sillen [12], which is similar to the cubic fluorite structure space group Pn-3m (group 224) with four Bi atoms in the following positions: 4(c): 3/4 3/4 3/4; 3/4 1/4 1/4; 1/4 3/4 1/4; 1/4 1/4 3/4 and six O atoms located at the following: 6(d): 0 1/2 1/2; 1/2 0 1/2; 1/2 1/2 0; 1/2 0 0; 0 1/2 0; 0 1/2 1/2, with the lattice parameter  Å (see Figure 1).

Figure 1: Unit cell for δ-Bi₂O₃, with distinct atomic species explicitly indicated. In the inset we show the associated Brillouin zone and the followed path in the energy bands calculation.

We performed self-consistent calculation of the electronic and atomic structure for δ-Bi₂O₃ using a relativistic full-potential method based on the augmented plane waves plus local orbitals (APW + lo) [19] as implemented in the WIEN2k code [20], within the framework of DFT [21–23]. The LDA [21, 22], PBE96 [23], and B3PW91 [24–26] energy functionals were used. The nonequivalency of PW91 and PBE functionals [27] was recently reported, which makes us cautious in the calculation of electronic properties of solid materials even though PBE is yet considered the standard approach to determine structural properties in solid state physics [28–30]. A comparison of results using LDA, PBE96, and B3PW91 functionals provides a broader insight into the use of density functionals for the prediction of physical and chemical properties in different bismuth oxide phases. It is worth mentioning that the hybrid B3PW91 functional includes 20% of Hartree-Fock exchange, which is not considered in LDA or PBE96.

In our calculations, the muffin-tin sphere radii were fixed at  a.u. for oxygen atoms and  a.u. for bismuth atoms. The maximum value for partial waves used inside atomic spheres was , while the magnitude of the largest vector in charge density Fourier expansion is fixed at . The basis set includes 2s²2p⁴ of oxygen and 5d¹⁰6s²6p³ of bismuth as valence states, because the electronic charge of these states is not contained completely within the muffin-tin spheres. In this computational code, the relativistic effects for core states are taken into account using a scalar-relativistic basis, while the spin-orbit coupling is calculated using the second-variation method [31–33].

The band structure of δ-Bi₂O₃ is calculated using 10 × 10 × 10 -point mesh in the first Brillouin zone (FBZ) for the SCF cycle, plotting the energy versus -point values following high symmetry points from to to to to of the value of to span the three-dimensional space (see inset in Figure 1). The density of states was computed using 20 × 20 × 20 -point mesh while the Fermi energy () and the weights of each energy band were calculated using a temperature broadening scheme with a value of 0.005. For both structural and electronic properties calculation we used a plane wave cut-off of .

3. Results and Discussion

3.1. Structural Properties

Structural properties were calculated from the ground state as a function of volume. A uniform compression and expansion of the lattice with the relative atomic positions within the unit cell held constant were performed to make an isotropic variation of the cell volume, thus searching for the stable crystalline structure. The calculated total energy for δ-Bi₂O₃ is displayed as a function of the cell volume in Figures 2(a)–2(c). The total energy values were used to fit the Murnaghan [34] equation of state from which we obtained the optimized lattice parameter and the bulk modulus. Notice that an exceptional fit was achieved for the cell expansion in the three studied cases.

(a)

(b)

(c)

Figure 2: versus volume curve for δ-Bi₂O₃, obtained with the functionals (a) LDA, (b) B3PW91, and (c) PBE96.

To obtain the equilibrium unit cell, total energy versus volume (see Figure 2) was calculated without the spin-orbit interaction using the following parameters: , initial experimental parameter  Å, and 10 × 10 × 10 -points sampling the first Brillouin zone (FBZ). This procedure yielded the equilibrium lattice parameters showed in Table 1. Unit cell optimization taking into account the spin-orbit interaction was also performed using the PBE96 functional, through the same procedure to understand whether it is significant or instead it is irrelevant in the calculation of the lattice parameter. The difference between both values is of the order of 0.01 Å; that is, there is a very small effect due to spin-orbit (SO) interaction, and thus the calculations of the electronic properties were made using the calculated equilibrium lattice parameters without SO effect. The lattice parameter calculated using the B3PW91 hybrid functional is in better agreement with the experimental value  Å.

Table 1: Lattice parameters for δ-Bi₂O₃ calculated in this work. The experimental value is 5.525 Å.

XC functional
(Å) without SO

LDA
5.4742

B3PW91
5.5389

PBE96
5.5946

[15].

3.2. Electronic Properties

Figures 3(a)–3(c) display energy bands versus -points for δ-Bi₂O₃, considering spin-orbit interaction. By examining the electronic band structure, we notice that the band shape is similar and that the Fermi level () is crossed by multiple degenerate bands indicating a metallic behavior, regardless of the XC functional. On the other hand, the addition of SO coupling splits orbital degeneracies into some high symmetry points in the FBZ. Figure 3 exhibits three distinctive groups of energy bands showing a band gap between them. The first group is below −7.5 eV, a second group is between −5.6 and 2.9 eV, crossing the Fermi level, and the third group is above 3.6 eV.

(a)

(b)

(c)

Figure 3: Calculated energy bands for δ-Bi₂O₃ considering spin-orbit interaction and using density functionals (a) LDA, (b) B3PW91, and (c) PBE96. Fermi energy () is located at 0.0 eV.

The total and projected density of states (DOS) for δ-Bi₂O₃, Bi, and O atoms, considering spin-orbit interaction, are provided in Figures 4(a)–4(c) where energy versus DOS is plotted. Notice that the d and f contributions for Bi are not present in the DOS graphs but are taken into account on its shape because they are included in the relativistic approach. Analysis of the band structure and partial DOS permits us to gain insight into their electronic (orbital) origin. The main contribution to the band from −8.0 to −10.0 eV comes from Bi 6s, identified as the bulk response given by the electrons closer to the nucleus. The following band −6.0 to +3.0 eV is formed by Bi 6s (0.4 states/eV) and Bi 6p (0.2–0.4 states/eV) orbitals, while O 2p orbital contributes with 1.4 states/eV. The Fermi level was shifted to 0.0 eV to perform an analysis of the contributions from each atom in the vicinity of . According to the orbital analysis, it is mainly formed by hybridized Bi 6s (0.4 states/eV) and O 2p (1.4 states/eV) orbitals. Finally, the conduction bands from +3.5 to 6.0 eV are mainly formed by Bi 6p (0.2–1.0 states/eV) contribution.

(a)

(b)

(c)

Figure 4: Total and projected density of states (DOS) for δ-Bi₂O₃ for the case when spin-orbit interaction is considered and using functionals (a) LDA, (b) B3PW91, and (c) PBE96.

To test the degree of covalency or ionicity in Bi₂O₃, electronegativity values (in the Pauling scale) of Bi (2.02) and O (3.44) give us a clue on the bonding character. Thus, based on the Pauling electronegativity scale [35], the electronegativity difference is about 1.42 which indicates to us that the amount of ionic character is ~40%. This suggests that the ionic character is significant in the bond between Bi and O but the covalent character is greater, and therefore we expect to find a fraction of electron density shared in-between. In order to verify the predicted ionic contribution on bonding type we carried out valence charge density calculations whose contours for δ-Bi₂O₃ are plotted in (100) and (110) planes, depicted in Figures 5(a)–5(f), where each atom has been labeled by its chemical symbol. These figures exhibit a similar charge distribution for the different XC potentials used and the existence of almost uniform charge distribution between Bi and O, which provides additional evidence of the ionic bonding feature of this compound. A quantitative analysis of the charge density between Bi and O atoms shows that the charge density is, at least, of 0.2 e/Å³. All calculations were carried out considering spin-orbit interaction effects.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5: Charge density profile for δ-Bi₂O₃ obtained by using LDA (a, b), B3PW91 (c, d), and PBE96 (e, f) exchange-correlation functionals. The charge profiles are plotted in (100) and (110) planes. Each atom has been identified with the corresponding chemical symbol.

4. Concluding Remarks

The structural and electronic properties of δ-Bi₂O₃ were calculated using WIEN2K employing relativistic considerations for Bi and O atoms and the exchange-correlation functionals LDA, PBE96, and the hybrid B3PW91. Structural optimization yielded an equilibrium lattice parameter of  Å,  Å, and  Å. The overestimation using PBE96 and underestimation using LDA with regard to the experimental value are within the expected on a GGA and LDA approach, respectively. A meticulous analysis of electronic properties was also performed. The energy bands were identified to have s and p character. Our results also provide evidence that the system behaves like a metal, considering the spin-orbit interaction. The major contributions to each band were identified from the projected DOS. The bands closer to (0 to −1.0 eV) are formed mainly from O 2p (1.4 states/eV) and Bi 6s (1.5 states/eV) hybridized orbitals. Finally, the charge density profiles exhibit the partial ionic chemical bonding features of this compound, with a small amount of charge between Bi and O atoms as a consequence of the partial covalent bonding. Even though most of the electronic properties describe similar behavior for each of the XC functionals used, the lattice parameter obtained with the B3PW91 functional is in better agreement with the experimental value, probably due to the existence of Hartree-Fock exchange in this functional.

Conflict of Interests

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

Acknowledgments

D. H. Galván acknowledges CONACYT under Grant 3027P-E and Proyecto de Supercomputo SC-004696 for the support provided. R. Núñez-González acknowledges Área de Cómputo de Alto Rendimiento (ACARUS) at Universidad de Sonora for giving access to their supercomputer system. A. Posada-Amarillas is grateful to Dr. L. Oliver Paz-Borbón from Chalmers University of Technology for useful discussions and valuable suggestions throughout this investigation.

References

A. F. Ioffe, Poluprovodnikovie Termoelementary, Izvestiya Akademii Nauk SSSR, Moskova, Russia, 1960 (Russian).

J. M. Thomas, D. A. Jefferson, and G. R. Willward, “The nanostructure of heterogeneous catalysts,” JEOL News E, vol. 23, no. 7, p. 7, 1985.

D. H. Galván, S. Fuentes, M. Avalos-Borja et al., “Structure and catalytic activity characterization of bismuth molybdate catalysts,” Catalysis Letters, vol. 18, no. 3, pp. 273–281, 1993.

K. Fossheim, E. D. Tuset, T. W. Ebbesen, M. M. J. Treacy, and J. Schwartz, “Enhanced flux pinning in Bi₂Sr₂CaCu₂O8+x superconductor with embedded carbon nanotubes,” Physica C: Superconductivity, vol. 248, no. 3-4, pp. 195–202, 1995.

A. J. Salazar-Pérez, M. A. Camacho-López, R. A. Morales-Luckie, V. Sánchez-Mendieta, F. Ureña-Ñúñez, and J. Arenas-Alatorre, “Structural evolution of Bi₂O₃ prepared by thermal oxidation of bismuth nano-particles,” Superficies y vacío, vol. 18, pp. 4–8, 2005.

B. Y. Liaw and W. Weppner, “Low temperature limiting-current oxygen sensors based on tetragonal zirconia polycrystals,” Journal of the Electrochemical Society, vol. 138, no. 8, pp. 2478–2483, 1991.

N. Cornei, N. Tancret, F. Abraham, and O. Mentré, “New ε-Bi₂O₃ metastable polymorph,” Inorganic Chemistry, vol. 45, no. 13, pp. 4886–4888, 2006.

E. C. Subbarao and H. S. Maiti, “Solid electrolytes with oxygen ion conduction,” Solid State Ionics, vol. 11, no. 4, pp. 317–338, 1984.

V. P. Zhuk, A. A. Vecher, and V. V. Samokhval, Vestn. Beloruss. Gos. Univ., vol. 2, no. 1, pp. 8–15, 1984.

T. Takahashi, H. Iwahara, and T. Arao, “High oxide ion conduction in sintered oxides of the system Bi₂O₃-Y₂O₃,” Journal of Applied Electrochemistry, vol. 5, no. 3, pp. 187–195, 1975.

H. A. Harwig and Z. Anorg, “On the structure of bismuthsesquioxide: the α, β, γ, and δ-phase,” Zeitschrift für Anorganische und Allgemeine Chemie, vol. 444, no. 1, pp. 151–166, 1978.

L. G. Sillén, “X-ray studies on Bismuth-Trioxili,” Arkiv för Kemi, Mineralogi och Geologi A, vol. 12, no. 18, pp. 1–15, 1937.

G. Gattow and H. Schröder, “Über wismutoxide. III. Die kristallstruktur der hochtemperaturmodifikation von Wismut(III)-oxid (δ-Bi₂O₃),” Zeitschrift für Anorganische und Allgemeine Chemie, vol. 318, no. 3-4, pp. 176–189, 1962.

N. I. Medvedeva, V. P. Zhukov, V. A. Gubanov, D. L. Novikov, and B. M. Klein, “Electronic structure and chemical bonding in δ-Bi₂O₃,” Journal of Physics and Chemistry of Solids, vol. 57, no. 9, pp. 1243–1250, 1996.

V. A. Gubanov and N. I. Medvedeva, “Electronic band structure and chemical bonding in the transition metal dioxides,” Physica B: Condensed Matter, vol. 172, no. 1-2, pp. 285–288, 1991.

R. A. Evarestov, V. O. Shapovalov, and V. A. Veryazov, “Electronic structure and chemical bonding in Bi₂O₃,” Physica Status Solidi (b), vol. 183, no. 1, pp. K15–K17, 1994.

J. A. Samaniego-Reyna, Estructura electronic del α-Bi₂O₃ con el metodo de enlace-fuerte Hückel Extendido [M.S. thesis], CICESE, Ensenada, Mexico, 1997.

A. Walsh, G. W. Watson, D. J. Payne et al., “Electronic structure of the α and δ phases of Bi₂O₃: a combined ab initio and x-ray spectroscopy study,” Physical Review B, vol. 73, no. 23, Article ID 235104, 2006.

G. K. H. Madsen, P. Blaha, K. Schwarz, E. Sjöstedt, and L. Nordström, “Efficient linearization of the augmented plane-wave method,” Physical Review B, vol. 64, Article ID 195134, 2001.

P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, and J. Luiz, WIEN2K, An Augumented Plane Wave + Local Orbital Program for Calculating Crystal Properties, WIEN2K 09.2, Vienna University of Technology, Vienna, Austria, 2009.

P. Hohenberg and W. Kohn, “Inhomogeneous electron gas,” Physical Review, vol. 136, article B864, 1964.

W. Kohn and L. J. Sham, “Self-consistent equations including exchange and correlation effects,” Physical Review, vol. 140, no. 4, pp. A1133–A1138, 1965.

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

A. P. Scott and L. Radom, “Harmonic vibrational frequencies: an evaluation of Hartree-Fock, Møller-Plesset, quadratic configuration interaction, density functional theory, and semiempirical scale factors,” Journal of Physical Chemistry, vol. 100, no. 41, pp. 16502–16513, 1996.

P. Geerlings, F. De Proft, and J. M. L. Martin, Theoretical and Computational Chemistry, vol. 4 of Recent Developments and Applications of Modern Density Functional Theory, Elsevier, New York, NY, USA, 1996, edited by: J. Seminario.

D. Torumba, P. Novák, and S. Cottenier, “Hybrid exchange-correlation functionals applied to hyperfine interactions at lanthanide and actinide impurities in Fe,” Physical Review B, vol. 77, no. 15, Article ID 155101, 2008.

A. E. Mattsson, R. Armiento, P. A. Schultz, and T. R. Mattsson, “Nonequivalence of the generalized gradient approximations PBE and PW91,” Physical Review B—Condensed Matter and Materials Physics, vol. 73, no. 19, Article ID 195123, 2006.

L. S. Pedroza, A. J. R. da Silva, and K. Capelle, “Gradient-dependent density functional of the Perdew-Burke-Ernzenhof type for atoms, molecules and solids,” Physical Review B, vol. 79, no. 20, Article ID 201106, 4 pages, 2009.

P. Hass, F. Tran, P. Blaha, K. Schwarz, and R. Laskowski, “Insight into the performance of GGA functionals for solid-state calculations,” Physical Review B, vol. 80, no. 19, Article ID 195109, 13 pages, 2009.

P. Haas, F. Tran, P. Blaha et al., “Systematic investigation of a family of gradient-dependent functionals for solids,” Physical Review B—Condensed Matter and Materials Physics, vol. 81, no. 12, Article ID 125136, 2010.

D. D. Koelling and B. N. Harmon, “A technique for relativistic spin-polarised calculations,” Journal of Physics C: Solid State Physics, vol. 10, no. 16, article 3107, 1977.

A. H. MacDonald, W. E. Picket, and D. D. Koelling, “A linearised relativistic augmented-plane-wave method utilising approximate pure spin basis functions,” Journal of Physics C: Solid State Physics, vol. 13, no. 14, article 2675, 1980.

D. Sing, Plane Waves, Pseudopotentials and the LAPW Method, Kluwer Academic Publishers, 1994.

F. D. Murnaghan, “The compressibility of media under extreme pressures,” Proceedings of the National Academy of Sciences of the United States of America, vol. 30, pp. 244–247, 1944.

L. Pauling, General Chemistry, chapter 6, Dover Publications, New York, NY, USA, 3rd edition, 1988.