The structural and electronic properties of the ternary Cd1-xBexSe alloys have been calculated using the full-potential linear muffin-tin-orbital (FP-LMTO) method based on density functional theory within local density approximation (LDA). The calculated equilibrium lattice constants and bulk moduli are compared with previous results. The concentration dependence of the electronic band structure and the direct and indirect band gaps are investigated. Moreover, the refractive index and the optical dielectric constant for Cd1-xBexSe are studied. The thermodynamic stability of the alloys of interest is investigated by means of the miscibility. This is the first quantitative theoretical prediction to investigate the effective masses, optical and thermodynamic properties for Cd1-xBexSe alloy, and still awaits experimental.
1. Introduction
In recent years, the wide-gap II–VI compounds are widely investigated because of the attractive applications in fabricating blue-green and blue optoelectronic devices, such as light-emitting diodes and laser diodes [1–4]. The applications include the use of II–VI compound based materials as light sources, in full colour displays, and for increasing the information density in optical recording [5, 6]. The beryllium containing II–VI compounds had been found to possess an enhanced ability to significantly reduce the defect propagation due to a greater prevalence of strong covalent bonding and lattice hardening in the materials [7, 8]. The strong covalent bonding in beryllium-based II–VI compounds achieves a considerable lattice hardening, which avoids multiplication of defects during the operation of II–VI semiconductor laser devices [9, 10] BexCd1-xSe alloys have attracted great attention because they are promising for the fabrication of full-colour visible optical devices due to a large difference in the energy gaps Eg of the binary constituents (CdSe, Eg=1.74 eV; BeSe, Eg=5.5 eV) [11].
In the present theoretical work, band gap of zinc-blende CdSe is varied systematically by alloying with Be. In order to investigate the optoelectronic nature of these alloys, their structural, electronic, and optical properties are calculated. All calculations are based on density functional full-potential linear muffin-tin orbital (FP-LMTO) method with perdew-wang local density approximation (LDA).
2. Method of Calculations
The calculations reported here were carried out using the ab initio full-potential linear muffin-tin orbital (FP-LMTO) method [12–15] as implemented in the Lmtart code [16]. The exchange and correlation potential was calculated using the local density approximation (LDA) [17]. The FP-LMTO is an improved method compared to previous LMTO techniques, and it treats muffin-tin spheres and interstitial regions on the same footing, leading to improvements in the precision of the eigenvalues. At the same time, the FP-LMTO method, in which the space is divided into an interstitial region (IR) and nonoverlapping muffin-tin spheres (MTS) surrounding the atomic sites, uses a more complete basis than its predecessors. In the IR regions, the basis set consists of plane waves. Inside the MT spheres, the basis set is described by radial solutions of the one particle Schrödinger equation (at fixed energy) and their energy derivatives multiplied by spherical harmonics. The charge density and the potential are represented inside the MTS by spherical harmonics up to lmax=6. The integrals over the Brillouin zone are performed up to 35 special k-points for binary compounds and 27 special k points for the alloys in the irreducible Brillouin zone (IBZ) using Blochl’s modified tetrahedron method [18]. The self-consistent calculations are considered to be converged when the total energy of the system is stable within 10−6 Ry. In order to avoid the overlap of atomic spheres, the MTS radius for each atomic position is taken to be different for each composition. We point out that the use of the full-potential calculation ensures that the calculation is not completely independent of the choice of sphere radii.
Structural properties of Cd1-xBexSe are calculated using Murnaghan’s equation of state [19] as follows:
(1)E(V)=E0+B0VB0′((V0/V)B0′B0′-1+1)-B0V0B0′-1,
where E0 is the total energy of the supercell, V0 is the unit volume, B0 is the bulk modulus at zero pressure, and B0′ is the derivative of bulk modulus with pressure.
Optical properties of Cd1-xBexSe are calculated using a fine k mesh of 1500 points for the present calculation. The dielectric function of a crystal depends on the electronic band structure and its investigation by optical spectroscopy which is a powerful tool in the determination of the overall optical behavior of a compound. It can be divided into two parts, real and imaginary as follows:
(2)ε(ω)=ε1(ω)+iε2(ω).
The imaginary part of the complex dielectric function, ε2(ω), in cubic symmetry compounds can be calculated by the following relation [20, 21]:
(3)ε2(ω)=82πω2∑nn′∫|pnn′(k)|2dSk∇ωnn′(k),
while ε1(ω) is used to calculate the real part of the complex dielectric function as follows:
(4)ε1(ω)=1+2πp∫0∞ω′ε2(ω′)ω′2-ω′2dω′.
Refractive index is calculated in terms of real and imaginary parts of dielectric function by the following relation
(5)n(ω)=12[{ε1(ω)2+ε2(ω)2}1/2+ε1(ω)]1/2.
3. Result and Discussion
In order to study the structural properties of Cd1-xBexSe(0<x<1), alloys are modeled at various compositions of Be with a step of 0.25. Structure optimization of each compound is performed by minimizing the total energy with respect to the unit cell volume and c/a ratio using Murnaghan’s equation of state [19]. The crystal structure of CdSe and BeSe zinc-blende with space group F4-3m (no. 216). Structural parameters such as lattice constant, a (Å), are calculated from the stable volume and are presented in Table 1. It is clear from the table that our calculated results for the binary compounds are in good agreement with the available experimental and calculated data.
Lattice constants a and bulk modulus B of Cd1-xBexSe compared with experimental results, Vegard’s law, and other theoretical calculations.
Figures 1 and 2 show the variation of the calculated equilibrium lattice constant and bulk modulus as a function of concentrations x for the Cd1-xBexSe alloy. The obtained results for the composition dependence of the calculated equilibrium lattice parameter almost follow Vegard’s law [22]. In going from CdSe to BeSe, when the Be-content increases, the values of the lattice parameters of the Cd1-xBexSe alloy decrease. This is due to the fact that the size of the Be atom is smaller than the Cd atom. On the opposite side, one can see from Figure 2 that the value of the bulk modulus increases as the Be concentration increases.
Variation in lattice constants a and bulk modulus B of Cd1-xBexSe as a function of composition x.
Calculated band structure of (a) Cd0.75Be0.25Se, (b) Cd0.5Be00.5Se, (c) Cd0.25Be0.75Se, and (d) band gap as a function of x.
The calculated band structure and partial density of states for Cd1-xBexSe(0<x<1) are presented in Figure 3. It is clear from the figure that Cd1-xBexSe(0<x<1) is a direct band gap material. The substitution of Be does not affect the direct band gap nature of the compound but increases the gap, which is clear from Figure 2(d). The direct band gap varies from 0.39 to 4.45 eV and the indirect band gap also increases from 1.89 to 3.71 eV with the increase in Be concentration between 0.25 and 0.75. It is obvious from the data presented in Table 2 that our calculated values for the band gaps of CdSe and BeSe are closer to the experimental results than the other calculated ones. The reason for our better results is the use of effective Perdew and Wang potential in the LDA scheme [17] and high k-points (1500). DFT always underestimates the band gaps; the origin of band structures is presented in Figure 3. In general, the band structures of these three compositions x (0.25, 0.5, 0.75) are very similar. The uppermost valence band is mainly formed by Se 4p states.
Fundamental direct and indirect band gaps of Cd1-xBexSe compared with the experimental and other calculations.
Frequency dependent imaginary part and real part of dielectric functions of Cd1-xBexSe.
The valence band maximum appears to be almost degenerate at the Γ and Xk points for three concentrations (0.25, 0.5, and 0.75), the energy at X being (1.04, 1.29, and 1.24) eV lower than that at Γ, respectively. On the other hand, the conduction band minimum occurs at Γ, so there is a direct gap of 0.84, 1.44, and 2.47 eV at Γ and an indirect Γ–X gap of 1.89, 2.74, and 3.71 eV for three concentrations (0.25, 0.5, and 0.75), respectively. To the best of our knowledge, there are no theoretical or experimental data on the energy band gaps for x=0.25, 0.5, and 0.75 available in the literature to make a meaningful comparison. It is clear from the results that the conduction band is mainly composed of Be-2p state for all ternary alloys. Figures 3(a) and 3(b) show that the lower part of the valence band is composed of Cd-4d and the upper part is mainly dominated by Se-4p state. This is due to the nature of bonding. At x=0.25, 0.5, and 0.75 partial covalent bond is stronger so the charge is shared by Se-4p and Be-2s states.
The calculated band gap versus concentration was fitted by a polynomial equation. The results are shown in Figure 2(d) and are summarised as follows:
(6)Cd1-xBexSe⟶{EΓ-Γ=0.466-0.088x+3.99x2EΓ-X=1.978+2.593x-1.737x2.
The variation in the band gap of Cd1-xBexSe provides promising results of the use of the compound in optoelectronic devices working in visible to ultraviolet regions. Depending on the need and requirement of a particular application, any desired band gap between 0.39 and 4.45 eV can be achieved.
We have calculated the frequency dependent imaginary dielectric function and real dielectric function. The effects of using k points in the BZ have already been discussed in the earlier work by Khan et al. [23]. The knowledge of both the real and the imaginary parts of the dielectric function allows the calculation of important optical functions. In this work, we also present and analyse the refractive index n(ω) given by (4).
The calculated imaginary part of the dielectric function for Cd1-xBexSe (x=0.25, 0.5 and 0.75) in the energy range 0–12 eV is shown in Figure 3. It is clear from the figure that for x=0.25, 0.50, 0.75, and 1.0 the critical points in the imaginary part of the dielectric function occur at about 0.85, 1.45, 2.48, and 4.45 eV, respectively. These points are closely related to the direct band gaps EGΓ-Γ; 0.84, 1.44, 2.47, and 4.45 eV of Cd1-xBexSe for the corresponding values of x = 0. 25, 0.50, 0.75 and 1.
The calculated real parts of the complex dielectric function ε1(ω) for Cd1-xBexSe are presented in Figure 3. It is clear from the figure that the static dielectric constant, ε1(ω), is strongly dependent on the band gap of the compound. The calculated values of ε1(ω) for Cd1-xBexSe at x=0, 0.25, 0.50, 0.75 and 1.0 are 4.35, 3.94, 3.74, 3.57, and 3.42 for corresponding direct band gaps 0.39, 0.84, 1.44, 2.47, and 4.45 eV, respectively. These data explain that the smaller energy gap yields larger ε1(0) value. This inverse relation of ε1(ω) with the band gap can be explained by the Penn model [24] as follows:
(7)ε1(0)≈1+(ℏωpEg)2.
The calculated values of the optical dielectric constant ε(ω) and refractive index n(ω) are listed in Table 3; comparison with the available data has been made where possible. As compared with other calculations, it seems that the values of n(ω) obtained from FP-LMTO method for the end-point compounds (i.e., CdSe and BeSe) are in good agreement with the theoretical results, together with the refractive index n(ω)=ε at zero pressure. Note that ε is obtained from the zero-frequency limit of ε1(ω), and it corresponds to the electronic part of the static dielectric constant of the material, a parameter of fundamental importance in many aspects of materials properties. It is clear from Figure 4 that the refractive index of the material decreases with the increase in the Be concentration.
Refractive index, optical dielectric constant of Cd1-xBexSe alloys for different compositions x.
Computed static optical dielectric constant and static refractive index as function of composition for Cd1-xBexSe.
Figure 4 shows the variation of the computed static optical dielectric constant and static refractive index versus composition for Cd1-xBexSe alloys. The computed static optical dielectric constant and static refractive index versus composition were fitted by polynomial equation. The results are summarized as follows:
(8)Cd1-xBexSe⟶{ε(0)=4.3400-1.5610x+0.6662x2n(0)=2.0745-0.3685x+0.1485x2.
4. Conclusions
Density functional calculations are carried out for the first time to investigate structural and optoelectronic properties Cd1-xBexSe. Structure as well as bonding nature of the material significantly varies with Be concentration. The lattice constant of the crystal decreases linearly with x. The calculated band structure predicts that the alloys have direct band gaps, which increase with the increase in x. On the basis of wide range of fundamental direct band gaps (0.39–4.45 eV) and indirect band gaps (1.89–3.71 eV) between 0.25 and 0.75, it can be concluded that the material can be used in optoelectronic devices working in the IR, visible, and UV regions of the spectrum.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
LuoH.FurdynaJ. K.The II-VI semiconductor blue-green laser: challenges and solution1995108104110.1088/0268-1242/10/8/001TamargoM. C.LinW.GuoS. P.GuoY.LuoY.ChenY. C.Full-color light-emitting diodes from ZnCdMgSe/ZnCdSe quantum well structures grown on InP substrates2000214-21510581063HaaseM. A.QiuJ.DePuydtJ. M.ChengH.Blue-green laser diodes19915911127212742-s2.0-1154431159310.1063/1.105472KatoE.NoguchiH.NagaiM.OkuyamaH.KijimaS.IshibashiA.Significant progress in II-VI blue-green laser diode lifetime19983432822842-s2.0-0032484826ViguéF.TourniéE.FaurieJ. P.Zn(Mg)BeSe-based p-i-n photodiodes operating in the blue-violet and near-ultraviolet spectral range20007624210.1063/1.125715BousquetV.TourniéE.LaugtM.VenneguesP.FaurieJ. P.Structural and optical properties of lattice-matched ZnBeSe layers grown by molecular-beam epitaxy onto GaAs substrates199770356410.1063/1.119234WaagA.FischerF.SchüllK.BaronT.LugauerH.-J.LiteTh.ZehnderU.OssauW.GerhardT.KeimM.ReuscherG.LandwehrG.Laser diodes based on beryllium-chalcogenides19977032802822-s2.0-0030783284ZhangJ. Y.ShenD. Z.FanX. W.YangB. J.ZhangZ. H.ZnBeSe epitaxy layers grown by photo-assisted metalorganic chemical vapor deposition2000214-21510010310.1016/S0022-0248(00)00036-1MuñozA.Rodríguez-HernándezP.MujicaA.Ground-state properties and high-pressure phase of beryllium chalcogenides BeSe, BeTe, and BeS1996541186110.1103/PhysRevB.54.11861WaagA.FischerF.LugauerH. J.LitzTh.LaubenderJ.LunzU.ZehnderU.OssauW.GerhardtT.MöllerM.LandwehrG.Molecular-beam epitaxy of beryllium-chalcogenide-based thin films and quantum-well structures19968027927962-s2.0-0000479565MaksimovO.GuoS. P.TamargoM. C.Be-Chalcogenide Alloys for Improved R-G-B LEDs: BexZnyCd1−x−ySe on InP2002229210051009SavrasovS. Y.SavrasovD. Y.Full-potential linear-muffin-tin-orbital method for calculating total energies and forces1992461912181121952-s2.0-3374633023210.1103/PhysRevB.46.12181SavrasovS. Y.Linear-response theory and lattice dynamics: a muffin-tin-orbital approach1996541647010.1103/PhysRevB.54.16470HohenbergP.KohnW.Inhomogeneous electron gas19641363B864B8712-s2.0-1064425025710.1103/PhysRev.136.B864KohnW.ShamL. J.Self-consistent equations including exchange and correlation effects19651404A1133A11382-s2.0-004211315310.1103/PhysRev.140.A1133SavrasovS. Y.Program LMTART for electronic structure calculations20052205-65555572-s2.0-2014437316810.1524/zkri.220.5.555.65067PerdewJ. P.WangY.Pair-distribution function and its coupling-constant average for the spin-polarized electron gas1992462012947129542-s2.0-000056464410.1103/PhysRevB.46.12947BlochlP.JepsenO.AndersenO. K.Improved tetrahedron method for Brillouin-zone integrations1994491622310.1103/PhysRevB.49.16223MurnaghanF. D.The compressibility of media under extreme pressures1944309244247WootenF.1972New York, NY, USAAcademic PressFoxM.2001Oxford University PressVegardL.Die konstitution der mischkristalle und die raumfüllung der atome19215117262-s2.0-2784458380010.1007/BF01349680KhanM. A.KashyapA.SolankiA. K.NautiyalT.AuluckS.Interband optical properties of Ni3Al1993482316974169782-s2.0-000049997310.1103/PhysRevB.48.16974PennD.Wave-number-dependent dielectric function of semiconductors1962128209310.1103/PhysRev.128.2093ZakharovO.RubioA.BlaseX.CohenM. L.LouieS. G.Quasiparticle band structures of six II-VI compounds: ZnS, ZnSe, ZnTe, CdS, CdSe, and CdTe1994501510780107872-s2.0-000167742110.1103/PhysRevB.50.10780HeydJ.PeraltaJ. E.ScuseriaG. E.MartinR. L.Energy band gaps and lattice parameters evaluated with the Heyd-Scuseria-Ernzerhof screened hybrid functional200512317410110.1063/1.2085170YuP. Y.CardonaM.2001Berlin, GermanySpringerChangY. H.ParkC. H.SatoK.Katayama-YoshidaH.First-principles study of the effect of the superexchange interaction in (Ga, Mn)V (V = N, P, As and Sb)20064912032082-s2.0-33747128580LuoH.GhandehariK.GreeneR. G.RuoffA. L.TrailS. S.DiSalvoF. J.Phase transformation of BeSe and BeTe to the NiAs structure at high pressure199552705810.1103/PhysRevB.52.7058NarayanaC.NesamonyV. J.RuoffA. L.Phase transformation of BeS and equation-of-state studies to 96 GPa1997561433810.1103/PhysRevB.56.14338LarefS.LarefA.Thermal properties of BeX (X = S, Se and Te) compounds from ab initio quasi-harmonic method20125111351402-s2.0-8005203726010.1016/j.commatsci.2011.07.016Gonzalez-DiazM.Rodriguez-HernandezP.MunozA.Elastic constants and electronic structure of beryllium chalcogenides BeS, BeSe, and BeTefrom first-principles calculations1997551404310.1103/PhysRevB.55.14043SamarthN.LuoH.FurdynaJ. K.QadriS. B.LeeY. R.RamdasA. K.OtsukaN.Growth of cubic (zinc blende) CdSe by molecular beam epitaxy19895426268026822-s2.0-3654909730510.1063/1.101033TsidilkovskiI. M.1982Amsterdam, The NetherlandsElsevier Science & Technology BooksCardonaM.Fundamental reflectivity spectrum of semiconductors with zinc-blende structure19613210215121552-s2.0-374314338910.1063/1.1777034Salcedo-ReyesJ. C.Electronic band structure of the ordered Zn0.5Cd0.5Se alloy calculated by the semi-empirical tight-binding method considering second-nearest neighbor132198207ZakharovO.RubioA.BlaseX.CohenM. L.LouieS. G.Quasiparticle band structures of six II-VI compounds: ZnS, ZnSe, ZnTe, CdS, CdSe, and CdTe1994501510780107872-s2.0-000167742110.1103/PhysRevB.50.10780HuangP. J.HuangY. S.FirsztF.Photoluminescence and contactless electroreflectance characterization of BexCd1−xSe alloys200719202620810.1088/0953-8984/19/2/026208IvanovS. V.Abstracts of the 9th International Conference on ll-VI Compounds1999Kyoto, Japan209(to be published in Journal of Crystal Growth, 2000)IvanovS. V.NekrutkinaO. V.KaygorodovV. A.Optical and structural properties of BeCdSe/ZnSe QW heterostructures grown by MBEProceedings of the 8th Nanostructures: Physics and Technology International SymposiumJune 2000St.Petersburg, RussiaGonzalez-DiazM.Rodriguez-HernandezP.MunozA.Elastic constants and electronic structure of beryllium chalcogenides BeS, BeSe, and BeTefrom first-principles calculations199755140431404610.1103/PhysRevB.55.14043FleszarA.HankeW.Electronic excitations in beryllium chalcogenides from the ab initio GW approach2000624246624742-s2.0-4243658375YimW. M.DismakesJ. B.StofkoE. J.PaffR. J.Synthesis and some properties of BeTe, BeSe and BeS197233250150510.1016/0022-3697(72)90032-7HamiziN. A.JohanM. R.Optical properties of CdSe quantum dots via Non-TOP based route2012784588467OuendadjiS.GhemidS.MeradjiH.HassanF. E. H.Density functional study of CdS1−xSex and CdS1−xTex
alloys20104812062112-s2.0-7644909204310.1016/j.commatsci.2009.12.029BieniewskiT. M.CzyzakS. J.Refractive indexes of single hexagonal ZnS and CdS crystals196353449649710.1364/JOSA.53.000496KootstraF.De BoeijP. L.SnijdersJ. G.Application of time-dependent density-functional theory to the dielectric function of various nonmetallic crystals200062707110.1103/PhysRevB.62.7071SrivastavaG. P.TutuncuH. M.GunhanN.First-principles studies of structural, electronic, and dynamical properties of Be chalcogenides20047008520610.1103/PhysRevB.70.085206WagnerV.LiangJ. J.KruseR.GundelS.KeimM.WaagA.GeurtsJ.Lattice dynamics and bond polarity of be-chalcogenides a new class of II-VI materials1999215187912-s2.0-0033242081