Recently, efforts have been made to distinguish semiconductor materials that are useful in scientific applications, as well as their relevant intrinsic properties. Semiconductors are used in a variety of remarkable device applications. Superlattices have been made available for use in many applications owing to developments in their production techniques, such as the strain-induced lateral ordering process [1, 2], molecular beam epitaxy [2, 3] and low-pressure chemical-vapor deposition [2, 4]. The development of superlattices is actively being researched by studying the properties of superlattices with reduced size and dimensionality and by understanding their properties within the context of designing new electronic and optoelectronics devices.
Superlattices (SLs) are structures that are made up of two types of semiconducting (s.c) materials, with one type (s.c1) acting as a quantum well and the other (s.c2) acting as a quantum barrier [2]. In the present work, we explore the structural and electronic properties of the II-VI semiconductor compounds BeTe and ZnSe. The most widely studied telluride/selenide superlattice (SL) structures are ZnTe/ZnSe and ZnSe/BeTe [5].
The objective of the present work is to extract physical parameters from the structural and electronic properties of (BeTe)n/(ZnSe)m superlattices (where n and m are numbers of monolayers; n = 1, 2, 3, 4, or 5; m = 1, 2, 3, 4, or 5) and to compare them with previous theoretical and experimental results. We seek to more carefully and accurately assess the effects of different superlattice configurations on electronic properties and, in particular, to observe the dependence of band gap behavior on the layers used.
We thus investigate the structural and electronic properties of nine (BeTe)n/(ZnSe)m superlattices with a tetragonal structure by using the full-potential linear muffin-tin orbital (FP-LMTO) method in the framework of density functional theory (DFT) within the local density approximation (LDA) for the exchange correlation functional. The organization of the present study is as follows: the adopted computational method is discussed in Section 2, the results are presented and discussed in Section 3, and conclusions and brief remarks are given in Section 4.
Density functional theory (DFT) is a powerful tool that is widely employed for the calculation of electronic and structural properties of solids and has been shown to yield relevant information about condensed matter phases and materials worth computing. In this work, the FP-LMTO method as implemented in the LmtART computer code [6, 7] was applied to perform first-principles total-energy calculations. This method is based on DFT, which is a universal quantum mechanical approach for many-body problems. In this approach, the quantum many-body problem of an interacting electron gas is mapped exactly onto a set of single particles moving in an effective local potential with the same density as the real system; the obtained one-electron equations are called the Kohn-Sham equations [8, 9].
In the LMTO method, space is divided into an interstitial region (IR) and non-overlapping muffintin (MT) spheres centered at the atomic sites. In the IR regions, the basis set consists of plane waves. Inside the MT spheres, the basis sets are described by radial solutions of the one-particle Schrödinger equation (at fixed energy), and their energy derivatives are multiplied by spherical harmonics. We have used the recently developed LmtART package (LmtART 7) with the electrons exchange-correlation energy described using the Perdew-Wang parameterization of the local density approximation (LDA) [10].
The details of the calculations are as follows: the charge density and the potential are represented inside the muffin-tin sphere (MTS) spherical harmonics up to lmax = 6. The k integration over the Brillouin zone is performed using the tetrahedron method [11] and is set up differently following the case. For SL(m, n), meshes of (6, 6, 6), (8, 8, 8) and (10, 10, 10) are utilized for m + n = 2, m + 2 = 4 and m + n = 6, respectively. The self-consistent calculations are considered to be converged within 10−6 for the total energy. The values of the sphere radii (MTS) and the number of plane waves (NPLW) used in the present calculations are listed in Table 1. We observed that NPLW does not vary with n + m; by contrast, the RMTS and Ecut vary with different configurations of the layers.
Input parameters: number of plane waves, energy cut-off and muffin-tin radii.
Compounds | NPLW (Total) | Ecut [Ryd] | MTS [a.u.] | ||||
---|---|---|---|---|---|---|---|
Be | Te | Zn | Se | ||||
(BeTe)1/(ZnSe)1 | 16242 | 137.1748 | 2.122 | 2.476 | 2.122 | 2.476 | |
(BeTe)2/(ZnSe)2 | 32458 | 137.0922 | 2.024 | 2.526 | 2.165 | 2.392 | |
(BeTe)3/(ZnSe)3 | 48690 | 137.6625 | 2.020 | 2.540 | 2.158 | 2.401 | |
(BeTe)1/(ZnSe)3 | 32458 | 137.0651 | 2.123 | 2.477 | 2.180 | 2.392 | |
(BeTe)3/(ZnSe)1 | 32458 | 137.3635 | 2.022 | 2.540 | 2.121 | 2.474 | |
(BeTe)2/(ZnSe)4 | 48690 | 137.3247 | 2.022 | 2.528 | 2.172 | 2.401 | |
(BeTe)4/(ZnSe)2 | 48690 | 137.6141 | 2.020 | 2.546 | 2.140 | 2.431 | |
(BeTe)1/(ZnSe)5 | 48690 | 137.0638 | 2.123 | 2.477 | 2.182 | 2.401 | |
(BeTe)5/(ZnSe)1 | 48690 | 138.0711 | 2.017 | 2.548 | 2.108 | 2.475 |
Electronic configuration of superlattices (BeTe)n/(ZnSe)m in the ground state is Be:[He] 2s2, Te:[Kr] 4d10 5s2 5p4, Zn:[Ar] 3d10 4s2 and Te:[Se] 4d10 4p4 4s2. As a first step, we determined the structural properties of the binary compounds in the zinc-blende (ZB) structure. The lattice constant (a0) was obtained by fitting the total energy as a function of volume to Birch’s [12] equation of state. We obtained the lattice constant, the bulk modulus and its pressure derivative from this numerical fitting procedure. The calculated structural parameters of binary compounds for LDA are presented in Table 2, along with the previous theoretical calculations and experimental data. Inspection of Table 2 shows that our LDA results are in reasonable agreement with experimental and available theoretical values. In the second step, we were interested in the quantum well superlattice consisting of binary compounds. In our case, (ZnSe) played the role of the barrier while (BeTe) acted as the well.
Calculated lattice parameter a ( Å), bulk modulus B0 (GPa), derived modulus B
Lattice constant (Å) | Bulk modulus | Eg(eV) | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
This work | Exp | Other works | This work | Exp | Other works | This work | Exp | Other works | This work | Exp (Γ-X) | Other works |
BeTe 5.588 | 5.617 a[13], | 5.581 b[14], c[15], | 56.448 | 66.8 a[13], | 62.36 b[14], c[15], d[16], e[17], f[18], g[19], | 3.62 | 4 a[13], | 3.69 b[14], c[15], | 1.907 | 2.7 h[20], | 1.76 g[19], d[16], |
5.53 d[16], e[17], | 71 d[16], e[17], | 3.38 d[16], e[17], | |||||||||
5.531 f[18], g[19], | 71 f[18], g[19], | 3.377 f[18], g[19], | |||||||||
Exp (Γ-Γ) | |||||||||||
ZnSe 5.65 | 5.667 i[21], j[22], | 5.624 k[23], l[24], m[25], | 59.68 | 64.7 i[21], | 71.82 k[23], l[24], | 3.96 | 4.77 i[21], | 4.88 k[23], l[24], | 1.0963 | 2.82 q[29], | 1.31 k[23], |
5.666n | 67.32 m[25], | 4.05n, 599 o[27], | 1.863 r[30], | ||||||||
5.578 o[27], p[28], | 71.84 o[27], p[28], | 4.57 p[28], | 1.83 s[31]. |
In the present work, the superlattices consist of a sequence of alternating n and m layers of ZnSe and BeTe along a specified growth direction. X-axis as the growth axis and supperlattice (m, n) with tetragonal symmetry have been chosen (n is the number of ZnSe monolayers, and m is the number of BeTe monolayers). The superlattice SL(1,1) is made up of alternate monolayers of BeTe and ZnSe, each monolayer containing two atoms; this alternation increases with m and n. The elementary cell volume is proportional to the number of monolayers. For the considered structures, we have performed structural optimization of the parameter of the host compound by minimizing the total energy with respect to the cell parameters. The total energies calculated as a function of the unit cell volume are fitted to the Birch-Murnaghan equation of state [12].
The deduced results are illustrated in Table 3. The data show that the SL11 lattice parameter interpolates the bulk material values, giving a good epitaxial interface and minimizing the creation of defects; for the other SL with m = n, the lattice parameter increases with the number of layers used. Finally, for n ≠ m, we observe a slight difference of the lattice parameter, especially when the number of BeTe layers increases. This difference is due to the atomic radius of the compound, which causes compression of the superlattice.
The calculated equilibrium constant a ( Å), bulk modulus B0 (GPa) and B
Compounds | a(Å ) | B0 (GPa) | |
---|---|---|---|
(BeTe)1/(ZnSe)1 | 5.619 | 58.947 | 4.99803 |
(BeTe)2/(ZnSe)2 | 11.241 | 59.829 | 4.14331 |
(BeTe)3/(ZnSe)3 | 16.826 | 62.328 | 4.16186 |
(BeTe)1/(ZnSe)3 | 11.242 | 64.827 | 4.93767 |
(BeTe)3/(ZnSe)1 | 11.230 | 60.270 | 3.47019 |
(BeTe)2/(ZnSe)4 | 16.847 | 64.97 | 4.38255 |
(BeTe)4/(ZnSe)2 | 16.830 | 61.299 | 3.17126 |
(BeTe)1/(ZnSe)5 | 16.863 | 65.856 | 4.64635 |
(BeTe)5/(ZnSe)1 | 16.802 | 58.065 | 3.36946 |
For the binary compounds and their super-lattices, the electronic band structures have been calculated at the equilibrium lattice. The most prominent features of the calculated band gaps are shown in Table 2 and Table 4. It is interesting to note that BeTe has an indirect band gap with the valence band maximum at G and conduction band minimum at X, whereas ZnSe has a direct band gap with the valence band maximum (VBM) and conduction band minimum (CBM) at G. The calculated band gaps for the binary compounds are in good agreement with the available theoretical results. We can see also that, due to the approximations used in our work, the band gap values are underestimated compared with the experimental data.
Calculated energy band gaps.
Compounds | Eg (eV) | Nature |
---|---|---|
(BeTe)1/(ZnSe)1 | 1.829038 | Gap direct |
(BeTe)2/(ZnSe)2 | 2.102459 | Gap direct |
(BeTe)3/(ZnSe)3 | 1.933643 | Gap direct |
(BeTe)1/(ZnSe)3 | 1.582407 | Gap direct |
(BeTe)3/(ZnSe)1 | 1.858339 | Gap indirect |
(BeTe)2/(ZnSe)4 | 1.706723 | Gap direct |
(BeTe)4/(ZnSe)2 | 1.971915 | Gap direct |
(BeTe)1/(ZnSe)5 | 1.427891 | Gap direct |
(BeTe)5/(ZnSe)1 | 1.852561 | Gap indirect |
In Fig. 1 and Fig. 2, we show the calculations of band structures for superlattices at m = n and m ≠ n. Examination of these figures shows that they exhibit both types of band gaps: for m = n superlattices there is a direct band gap at Γ whereas for m > n, the superlattice band gap is indirect, with VBM located at G and CBM at R (except for m = 4 and n = 2, where the gap is direct but is close to G-R). The values of these gaps are 1.427 and 2.10 eV, respectively. No previous superlattice results are available for comparison; from the binary compound data, we conclude that the underestimation of band gaps by LDA also applies to superlattices. This underestimation is due to strongly correlated 3d electrons of Zn.
To describe the number of states that are available to be occupied by electrons per interval of energy at each energy level, we have also calculated the total and the partial densities of states (DOS) of these compounds, as displayed in Fig. 3 and Fig. 4.
Inspection of Fig. 3 shows that there are differences in the conduction band densities of states of the two binary compounds, with 3 regions present in the ZnSe DOS and two regions present in the BeTe DOS. This difference is due to the existence of the middle region dominated by the 3d-Zn states of the ZnSe compound; by contrast, the Be d sates are empty for the BeTe compound.
Examination of Fig. 4 shows that the topology of the superlattices densities of states is the same as that of the ZnSe binary compound. Two examples are shown in Fig. 4: m = n =1, 2, 3 and m ≠ n with m = 3 and n = 1. The lower occupied bands, located between ∼ 13 and
In this section, we discuss the optical properties of a material that must be investigated to determine its potential usefulness in optoelectronic applications. For this reason, we only chose the materials that showed a direct band gap character in our LDA study. When examining the optical response of the compounds under investigation, it is convenient to take into account the transitions of electrons from the occupied energy bands to the unoccupied energy bands, particularly at the high symmetry points in the Brillouin zone. The real part
To calculate the optical spectra of the dielectric function,
To illustrate a minor difference in the description of these properties, the dielectric function ε2(
In Fig. 6 and Fig. 7, the refractive index n(
We have performed detailed investigation on the structural, elastic, electronic, and optical properties of (BeTe)n/(ZnSe)m superlattices using the full-potential linear muffin-tin orbital method within the Perdew-Wang LDA. The calculated ground-state properties of binary compounds are in good agreement with the available experimental data. Our results for band structure and DOS show that our superlattices are semiconductors, but their band gaps change from direct to indirect (Γ-Γ to Γ-R) depending on the configuration of the layers. The results obtained for the energy band gaps using LDA show a strong dependence on the number of layers used. The imaginary and real parts of the dielectric function were investigated and analyzed to identify the optical transitions. The static dielectric constants