Skip to main content

REVIEW article

Front. Phys., 14 April 2022
Sec. Optics and Photonics
This article is part of the Research Topic Novel Phases in Exciton-Polariton Bose-Einstein Condensates View all 6 articles

Simulation of the Band Structure of InAs/GaSb Type II Superlattices Utilizing Multiple Energy Band Theories

Shuiliu FangShuiliu Fang1Ruiting Hao
Ruiting Hao2*Longgang ZhangLonggang Zhang1Jie GuoJie Guo1Wuming Liu,
Wuming Liu3,4*
  • 1School of Physics and Electronic Information, Yunnan Normal University, Kunming, China
  • 2School of Energy and Environment Science, Yunnan Normal University, Kunming, China
  • 3Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, China
  • 4Songshan Lake Materials Laboratory, Dongguan, China

Antimonide type II superlattices is expected to overtake HgCdTe as the preferred materials for infrared detection due to their excellent photoelectric properties and flexible and adjustable band structures. Among these compounds, InAs/GaSb type II superlattices represent the most commonly studied materials. However, the sophisticated physics associated with the antimonide-based bandgap engineering concept started at the beginning of the 1990s gave a new impact and interest in the development of infrared detector structures within academic and national laboratories. InAs/GaSb superlattices are a type II disconnected band structure with electrons and holes confined in the InAs and GaSb layers, respectively. The electron miniband and hole miniband can be regulated separately by adjusting the thickness of InAs and GaSb layers, which facilitates the design of superlattice structures and optimizes the value of band offset. In recent years, both domestic and foreign researchers have made many attempts to quickly and accurately predict the bandgaps of superlattice materials before superlattice materials grow. These works constituted a theoretical basis for the effective utilization of the InAs/GaSb system in material optimization and designing new SL structures; they also provided an opportunity for the preparation and rapid development of InAs/GaSb T2SLs. In this paper, we systematically review several widely used methods for simulating superlattice band structures, including the k·p perturbation method, envelope function approximation, empirical pseudopotential method, empirical tight-binding method, and first-principles calculations. With the limitations of different theoretical methods proposed, the simulation methods have been modified and developed to obtain reliable InAs/GaSb SL energy band calculation results. The objective of this work is to provide a reference for designing InAs/GaSb type II superlattice band structures.

Introduction

Infrared detectors and lasers are employed in both military and civil applications. As more advantages of type II superlattice (T2SL) materials are discovered, T2SL infrared lasers and detectors containing InAs, GaSb, and AlSb compounds are expected to be more widely used after HgCdTe and multi-quantum wells [1]. Unlike HgCdTe, the antimonide type II SLs exhibit a high level of reproducibility and maneuverability, large area uniformity, and low Auger recombination rates, which means that T2SL infrared lasers and detectors have lower dark currents, high-temperature operational characteristics [2], and distinct advantages in some application scenarios. For a long period, many laboratories worldwide have invested manpower and resources to perform theoretical simulations in the field of energy bands and achieved significant progress.

T2SLs were originally proposed by Esaki and Tsu in 1970 [3]. They represent periodic structures composed of two or more semiconductor layers of III–V materials with a type II band alignment and lattice constant of approximately 6.1 Å. Theoretical simulations of T2SLs related to energy band engineering applications have been initiated in the early 1990s [4]. Among these materials, InAs/GaSb SLs exhibit high flexibility in terms of bandgap adjustment and heterostructural design in the mid-wave infrared (MWIR) and long-wave infrared (LWIR) regions due to their unique properties including the lower position of the InAs conduction band than that of the GaSb valence band [5], low Auger recombination rate [6], and large effective mass [7], which attracted considerable attention from researchers [8]. However, the performance of T2SL devices is significantly lower than theoretical predictions. Moreover, InAs/GaSb SLs have several disadvantages, including small carrier lifetimes, low quantum efficiencies, and high dark currents due to the generation-recombination, which are attributed to the lack of a clear understanding of their band structure and topology [9]. Therefore, simulating SL band structures is the most important step in designing SL infrared lasers and detectors [10]. To achieve this goal, it is necessary to select appropriate theoretical methods, establish accurate device models, and propose new design improvement strategies by analyzing the physical properties of T2SL materials and/or devices. For this purpose, multiple studies on the material growth, electronic properties, and structural design of InAs/GaSb SLs have been conducted.

This article outlines the energy band structure of InAs/GaSb SLs, describes in detail several theoretical simulation methods of solid energy band commonly used for studying SL energy band structures, such as the k·p perturbation method, envelope function approximation, empirical pseudopotential method, empirical tight-binding method, density functional theory, and many-body perturbation theory. In the last section, it discusses the bottleneck problems and development trends of SL energy band simulation techniques.

Type II InAs/GaSb SL Electronic Band Structure

InAs/GaSb SLs were originally developed by Sai‐Halasz et al. in 1977 [11]. It is a periodic structure formed by InAs (a = 6.0584Å) and GaSb (a = 6.0959Å) grown alternately for several cycles. And it is worth mentioning that the superlattice bandgap is determined by the energy difference between the first electron miniband E1 and the first heavy hole miniband HH1 at the center of the Brillouin zone. The bottom of the InAs conduction band is much lower than the top of the GaSb valence band, which corresponds to a staggered type II band alignment (Figure 1A) with electrons and holes confined in the InAs and GaSb layers, respectively [12]. There is also a mutual coupling between the wavefunctions in the quantum well and the transition only occurs in the spatial region where the wave functions overlap, which broadens the electron and hole sublevels to form an energy band with a certain width (Figure 1B). It has been confirmed in the research model of Becer, et al [13]. The separation of electrons and holes in the real space not only effectively suppresses the Auger recombination of carriers, but also enables the independent adjustments of the electron and hole potential wells to achieve continuous light absorption in the wavelength range of 2–30 μm [14].

FIGURE 1
www.frontiersin.org

FIGURE 1. (A) Band structure alignment of InAs/GaSb at 0°K. In the InAs/GaSb superlattice, the conduction band of InAs is about 0.15 eV lower than the valence band of GaSb, and its heterojunction forms an off-type energy band. (B) InAs/GaSb superlattice energy band structure is a staggered type II band structure, causing electrons and holes to be confined in the InAs layer and GaSb layer, respectively.

The use of T2SLs for the fabrication of lasers and detectors depends [15] not only on the ability to grow a perfect periodic crystal structure but also on the material band-gap design [16]. During the selection of a device cut-off wavelength, the SL bandgap can be theoretically adjusted by varying the thicknesses of the InAs and GaSb layers and thus the degree of overlap of SL electronic wavefunctions. Delmas, et al. applied simulation tools to model and design high-performance InAs/GaSb T2SLs infrared detectors, showing that the SL design can improve overall device performances [17]. Today, many laboratories also showed that the experimental absorption spectra of the MWIR and LWIR InAs/GaSb and InAs/InAsSb T2SLs could be accurately simulated [4, 18]. The emergence of bandgap engineering has promoted the development of infrared detector structures.

In other words, to better understand the properties of T2SLs, their band structure must first be theoretically simulated. The theoretical methods currently used for this purpose include the k·p perturbation method, envelope function approximation (EFA), empirical pseudopotential method (EPM), empirical tight-binding method (ETBM), and first-principles calculations.

Theoretical Simulation Methods

k·p Perturbation Method

The k·p model is called a standard model because its calculation results are easily explained by the physical theory. In addition, the related calculation procedure is relatively simple and requires inputting a few parameters to solve the Schrödinger equation by changing the electron potential energy, while the SL band structure can be obtained from the wavefunction. The k·p perturbation method is based on the envelope function and effective mass approximations. It was introduced by Bardeen [19] and Seitz [20], developed by Bastard [21] in 1988, and applied to study the T2SL structures of infrared photoelectric detectors by Smith and Mailhiot [22].

After Read and Shockley extended the k·p perturbation method to an eight-fold degeneracy, the eight-band k·p matrix was used for the InAs/GaSb SL energy band simulation [23]. Klipstein derived an 8 × 8 envelope function Hamiltonian for the SL structure of Г15v and Г1c in 2010 [24]. Subsequently, Livneh et al. reduced the computational error by introducing an eight-band Hamiltonian to eliminate the terms with energies below the expected accuracy level [25]. Using this approach, the same researchers ultimately decreased the number of unknown fitting parameters to six. Afterward, they applied the eight-band k·p model for fitting the absorption spectra of InAs/GaSb SLs in the wavelength region of 4.3–10.5 μm at temperatures of 77 and 300°K to determine the six Luttinger and interfacial parameters. The fitted Luttinger parameters were very close to those originally proposed by Lawaetz [26], which indicated that the eight-band k·p model retained the high calculation accuracy. Finally, they used this model to predict the wavelengths of more than 30 SLs. The obtained photoluminescence (PL) spectra demonstrated that most errors did not exceed 0.3 μm, the maximum error was 0.6 μm, and the corresponding layer thickness error was less than 0.4 ML.

In 2010, Rejeb et al. simulated the structure of short-period InAs/GaSb/InSb SLs on GaSb substrates using the eight-band k·p method and plotted the fundamental bandgaps of the mutated and separate interfaces as functions of the period number N [27]. The obtained results indicated that in the case of interfacial interactions, the asymmetric interfacial segregation could lead to a bandgap reduction of approximately 30%. In 2012, Qiao et al. developed a band structure model using an eight-band k·p method with Dirichlet and periodic boundary conditions that included the actual interfacial layers, which was achieved by changing the previous way of adjusting valence band offset values or using potential gradient profiles [28]. Klipstein et al. used this k·p model to simulate InAs/GaSb SLs in 2014, which also considered interfacial effects and band bending [29]. The bandgap measured by PL at a temperature of 10°K corresponds well with the responsivity cut-off energy measured at 77°K [25, 30]. It is also well demonstrated that the calculated bandgap was consistent with the PL peak energy measured at a temperature of 10°K (Supplementary Figure S1, Supporting Information). If the interfacial matrix was ignored, the calculated 8.4 ML InAs/13.7 ML GaSb MWIR T2SL bandgap would exhibit a blue shift of 0.75 μm, and the 14.4 ML InAs/7.2 ML GaSb LWIR T2SL bandgap would produce a blue shift of 4.5 μm. Then, they tried to get a better fit of the measured and calculated InAs/GaSb T2SL MWIR and LWIR absorption spectra by fitting the parameter variables bandgap E0 and VBO [30]. They mentioned that E0 has a strong influence on the Luttinger parameter and VBO affects the optical transition energies between different mini-bands in the SLs and proved them. In 2019, Delmas, et al. applied the eight-band k·p method to calculate the cut-off wavelengths of InAs/GaSb SLs with four different period lengths: 10/4, 12/4, 14/4, and 17/4, and compared the results with or without considering the interface matrix HIF and the InSb interface [31]. Figure 2 showed that there was good agreement between simulation and experiment, and the error was within the deviation range. If both the HIF and the InSb interface are not considered, the model didn’t predict the measured cut-off wavelength and underestimated it. Therefore, the effect of the interfaces cannot be negligible when the eight-band k·p model is used to calculate the SL energy band structure.

FIGURE 2
www.frontiersin.org

FIGURE 2. Comparison between the calculated and measured cut-off wavelengths at 77°K for the different SL periods (red squares) along with the ideal prediction line (dashed line). Cut-off wavelengths calculated without the interface matrix (triangles) and considering neither the interface matrix nor the InSb layers (circles) are also plotted for comparison. The deviation in the predicted is represented by the solid lines and the grey area. The 10 ML (InAs)/4 ML (GaSb) SL was under compressive strain on GaSb with a large lattice mismatch and error.

The above-mentioned k·p method only considers the eight-fold degeneracy. Hence, without taking into account interfacial effects, the eight-band k·p method overestimates the SL bandgap by neglecting the interactions one neglects the interaction between near-gap bands of the GaSb layer and higher-lying bands of the InAs layer. However, a more accurate bandgap value can be calculated by a modified k·p method described below.

In 2002, Vinter [32] proposed an 18-band k·p method to describe the wavefunction and band structure of SL. The inclusion of these and only these 18 bands insures that all relevant symmetries are included in a k·p matrix. The 18-band k·p method can include the near-gap interactions. The theoretical and experimental energy values were very small differences (Supplementary Figure S2, Supporting Information). In 2015, Imbert et al. used the 18-band k·p model to simulate three devices with cut-off wavelengths of 5 μm but different lattice periods, which corresponded to the thickness ratios R between the InAs and GaSb layers equal to 0.5 (10 ML InAs/19 ML GaSb per period), 1 (10 ML), and 2 (7 ML InAs/4 MLs GaSb) [33, 34]. The relationship between the bandgap energy and the period thickness determined for the symmetric structure (R = 1) demonstrated a good agreement between the measured and those calculated bandgaps by the 18-band k·p model [35] (Supplementary Figure S3, Supporting Information).

The articles studying the 18-band k·p method are sparse. Researchers are more committed to modifying the 8-band k·p method to simulate the band structure of InAs/GaSb SLs. In 2018, Machowska–Podsiadlo et al. used the eight-band k·p method to study the effects of temperature, band offset energy, strain, and interface on the edge of the SL absorption band and relative contributions of these parameters [36]. The authors studied m/8 ML and 8/n ML SLs more comprehensively (Supplementary Figure S4, Supporting Information). They found that regardless of the GaSb layer thickness, the effect of strain on the absorption edge of the 8/n ML SLs remained the same; however, for the m/8 ML SLs, this effect was more pronounced at larger thicknesses of the InAs layers. Moreover, the bandgap energy computed at Eoffset = 140 meV was 7–9 meV higher than that calculated at Eoffset = 150 meV, which also agreed with the experimental value [37]. The authors concluded that the leap energy level in SLs mainly depended on the interface type. The results of the last study can help design more complex InAs/GaSb SL structures.

In 2019, Jeffrey established an accurate 8-band k·p model by using the 8-band k·p program nextnano3 to solve for the electronic structure of the material and Python scripts to implement for the optical absorption calculations. And they presented an optimization process to optimize the parameters in the experimental data and calculated the band structure of the InAs/GaSb superlattice [38]. The calculated absorption spectrum shape of 4–10 μm was consistent with the experimental data. The calculated absorption cut-off values were consistent, confirming that the position of the quasi-Fermi level can be set appropriately (Supplementary Figure S5, Supporting Information). Du, et al. [39] investigated the electronic band structure of InAs/GaSb SLs by 8-band k·p method and compared it with density functional theory (DFT) [40], empirical pseudopotential method (EPM) [41], and experimental results [42]. It verified the accuracy of the 8-band k·p theoretical model (Supplementary Figure S6, Supporting Information). In the same year, Cui, et al. designed an M-structure T2SL detector with a cut-off wavelength of 10.5 μm based on the eight-band k·p model and studied its photoelectric performance [43]. In 2021, Mukherjee, et al. simulated the 8 ML InAs/8 ML GaSb T2SL band structure [44]. The calculated bandgap was 0.27 eV within the k·p model under the envelope function approximation at 77°K, and the corresponding cut-off wavelength was 4.59 μm. It was in good agreement with the experimental bandgap in the range of 0.269–0.275 eV. Subsequently, Mukherjee, Singh et al. proposed in their test of the above model that the key factor in designing detectors using quantum-limited SL structures for a specific wavelength range is the precise control of the constituent layers thickness. Due to the fact that InAs/GaSb T2SL is an II-type broken-band arrangement, the appropriate values of VBO are important to model these non-common atomic structures. The wavelength range of the optical response spectrum of a material depends on bandgap and bandwidth. Therefore, they investigated the effect of varying tInAs and tGaSb on the bandgap and relative band edge values (Supplementary Figure S7, Supporting Information). Kim’s theory that the band gap can be effectively adjusted by adjusting the thickness of InAs or GaSb is fully consistent with the above findings [45] (Figure 3). This study offers an engineering guideline for the design of T2SL detectors.

FIGURE 3
www.frontiersin.org

FIGURE 3. The calculated bandgaps using the modified K·P model. The x and y axes represent the number of InAs and GaSb monolayer, respectively. The calculated bandgap of 10 ML InAs/1 ML InSb/10 ML GaSb T2SL is 0.2 eV at 120°K. And as the thickness of InAs increases, while keeping the thickness of the GaSb layer constant, the cut-off wavelength of the InAs/GaSb T2SL also increases.

The presence of InSb and GaAs interfaces in InAs/GaSb SL systems gives rise to microscopic interfacial asymmetry (MIA). In 2010, Li, Xu et al. used the finite-difference method to solve the 8-band k·p model including the MIA effect, and obtained the band structure of the InAs/GaSb SL with the layer widths LA = LB = 25 Å [46]. They point out that the MIA effect leads to a significant reduction of the bandgap. Dong, Li et al. used the eight-band k·p model incorporating MIA effects to calculate the band gap of short-period InAs/GaSb T2SLs in 2015 [47]. They fitted different InAs and GaSb layer widths characterizing the MIA effect to the experimental bandgap. The calculated results are shown in Figure 4 and the bandgap calculated by the model including the MIA effect is in good agreement with the experimentally obtained bandgap. Ignoring the MIA effect often leads to an overestimation of the SL bandgap. In 2019, Le, Kamakura et al. did the same work and obtained conclusions consistent with the above [14, 48]. They showed that the MIA effect affects the hole band dispersions and the rate band shift and this interaction is most pronounced in the shorter period SLs. All the above studies indicate that MIA effects play an important role in short-period InAs/GaSb SLs.

FIGURE 4
www.frontiersin.org

FIGURE 4. Band-gap energy calculated without the MIA (dashed line) and with MIA (solid line) and measured (black symbols) bandgaps the of symmetrical 10 ML InAs/10 ML GaSb SL (R = 1) on a GaSb substrate at 77°K plotted as functions of the SL period thickness.

In the 8-band k·p method, the numerical solution of coupled ordinary differential equations mainly includes the finite element method(FEM), the finite difference method(FDM) and the transfer matrix method(TMM). FEM can handle high-order quantization and give a very accurate numerical result from the isolated system of a single quantum well to more complex periodic potential problems. However, due to the complexity of its own and variational functional for a given multiband Hamiltonian, the calculation speed of the finite element method is limited by the integration method [13, 49]. When the width of the quantum well becomes thicker, the numerical calculation of TMM sometimes fails due to growing exponentials and there is a serious numerical instability [46, 50]. However, FDM is widely used because of its simplicity and accuracy. The method is conceptually simple and uses a layer or grid envelope function to construct the Hamiltonian system. Therefore, FDM can describe as accurately as possible the most relevant parts of the electronic subband structure of a multilayer QW, and therefore can handle complex boundary conditions and arbitrary potentials and obtain stable numerical properties [27, 49].

In summary, the k·p standard model is intuitive and easy to explain. It is an effective and widely used method for determining the energy band structure near the bottom of the conduction band and the top of the valence band in semiconductors. However, this technique utilizes the effective mass approximation while ignoring atomic parameters and cannot accurately describe the electronic structure of a short-period SL [51], and it may provide a similar result in MWIR and LWIR. Therefore, it is mainly used to simulate quantum wells, SL quantum dots, and long-period materials, while a new k·p theoretical method suitable for short-period T2SLs must be developed separately. In addition, the k·p model contradicts the low-momentum hypothesis of the k·p theory in many applications, which is controversial on a theoretical basis.

Envelope Function Approximation

In the early studies on InAs/GaSb SLs, the numerical accuracy of the standard EFA algorithm was considerably increased [52, 53]. The early model still did not consider the influence of interface on the SL band structure. As a result, the bandgaps of the first conduction band and heavy hole microstrip were significantly overestimated [54].

In 2004, Szmulowicz et al. proposed an improved 8 × 8 EFA method, which considered the effects of anisotropy and interfacial couplings for non-coatomic SLs [37]. In their study, a 44.16 Å GaSb/55.46 Å InAs SL structure was combined with an InSb interface. The simulated bandgap was 111.4 meV, which was closer to the experimental value than the magnitude of 116.9 meV computed by the 14 × 14 k·p perturbation method. Therefore, the modified EFA method including interfacial effects can provide more accurate simulation results for SL systems. In 2005, J. et al. studied m ML InAs/m ML GaSb SLs with and without additional interfacial potentials using the same model [54]. As shown in Figure 5, the fitted data of the modified EFA model considering strong perturbations at the interface are in good agreement with the experimental values. Subsequently, Szmulowicz et al. selected a smoother InSb interface with higher carrier mobility. Using the modified 8 × 8 EFA model [55], they designed an MWIR superlattice with a bandgap of 310 meV and cut-off wavelength of 4 μm. The calculated bandgaps of the 23.9 Å GaSb/20.4 Å InAs SL with a shorter period and 9.9 Å GaSb/11.4 Å InAs SL with a longer period were 304.2 and 313.8 meV, respectively. The bandgap values computed without taking into account the interfacial effect were equal to 391.5 and 514.6 meV, respectively.

FIGURE 5
www.frontiersin.org

FIGURE 5. Evolution of the C1–VH1 fundamental interminiband transition energy of the InAs(m)/GaSb(m) SL structure with the layer thickness: theoretical calculations conducted with (solid line) and without (dashed line) an additional interface potential; experimental data extracted from PL measurements in this work (solid square) and from other works.

In the same year, Haugan et al. designed 4-μm short-period InAs/GaSb SLs using the EFA model [56]. The measured peak position consistently maintains a constant level of 3.757 ± 0.10 μm (330 ± 10 meV) despite the large thickness variations from 50.2 to 21.2 Å, which is slightly lower than the predicted value (Supplementary Figure S8, Supporting Information). However, the difference of 10–30 meV is considered a normal experimental error corresponding to relatively high accuracy. In 2011, Debbichi et al. investigated the electronic and optical properties of the short-period InAs/GaSb/InSb SLs grown on GaSb substrates by the modified 8 × 8 EFA model that took into account the effects of anisotropy and interfacial interactions [57]. The obtained results revealed that the bandgaps calculated at different temperatures were in good agreement with experimental data, which confirmed the high accuracy of the utilized model. In 2013, Yi Zhou et al. used the InAs/InAsSb/GaSb/InAsSb four-layer SL (including the interface) to replace the standard InAs/GaSb double-layer structure in the EFA model and optimize the n ML InAs/12 ML GaSb SL band structure [58]. The obtained fitting data presented in Figure 6 indicate that the cut-off wavelength error of the four-layered structure is less than 5%, which is much lower than that of the standard SL model.

FIGURE 6
www.frontiersin.org

FIGURE 6. Calculated band-structure of 13.5 ML/2 ML/8 ML InAs/AlSb/GaSb T2SL showing dispersion of energy with respect to electron wave vector in plane (kxy) and in the growth direction (kz) for LWIR T2SL. The in plane dispersion direction along [109] is shown on the left-hand, whereas the dispersion in the growth direction [001] is shown on the right hand. Conduction bands 1 and 2 zone edges at (π/P) are presented as verticle dashed line (P is periodical length).

Boutramine et al. applied the EFA method to study the 25 Å InAs/25 Å GaSb SL bandgap as a function of the InAs layer thicknesses d1 and temperature [59]. With increasing d1, the electronic energy E1 of the InAs layer decreased, and the heavy-hole energy HH1 of the GaSb layer increased. This result was consistent with the values predicted using the k·p method [60]. As the temperature increased, the bandgap decreased from 288.7 meV at 4.2°K to 230 meV at 300°K. The corresponding cut-off wavelengths were equal to 4.3 and 5.4 μm, respectively, which belonged to the MWIR range. These parameters were in good agreement with experimental values [61]. In 2020, Benchtaber et al. calculated the band structure and bandgap of the 21 Å InAs/24 Å GaSb SL using the EFA model [62]. The bandgaps obtained at 5 and 300°K were equal to 316 and 247 meV, respectively, and the former value was consistent with the PL spectrum recorded at 300 meV [63]. The authors suggested that the small difference of 16 meV (0.05%) might be due to a valence band shift or the low thickness measurement accuracy. Boutramine, et al. also used EFA to do a lot of research [64, 65]on InAs/GaSb SL band structure, sub-bands and effective masses of carriers with different periods and the valence band offset Λ [66]. These findings are consistent with the experimental results reported by Cervera, et al. [61]. Hostut, et al. analyzed bandgap energy and hh–lh splitting energy of the InAs/AlSb/GaSb structure T2SL (N-type) using the EFA method [67]. The bandgap of the structure was obtained as 144 meV, corresponding wavelength of 8.6 μm, which lay in the LWIR of the atmospheric window. And the hh–lh splitting energy was 166 meV, which was 22 meV higher than the bandgap (Supplementary Figure S9, Supporting Information). It is well known that the larger hh–lh splitting energy is very important for suppressing the Auger recombination and improving the minority carrier lifetime.

The EFA model requires a large number of numerical calculations and thus represents a complex theoretical approach. It is suitable not only for calculating the band structure of long-period SLs, especially the electronic states near the Г point of the Brillouin zone (BZ), but also for heterostructural modeling. Some researchers have suggested that certain interface envelope function boundary conditions cause other uncertainties [22, 68]. For example, the EFA method produces similar results for different thin-layer SLs, which complicates the analysis and drawing conclusions from calculation data.

Empirical Pseudopotential Method

EPM is an atomic calculation method that is advantageous for studying short-period or thin-layer material structures such as semiconductors and metals. First, the model determines an SL form factor from a large number of band parameters. Alternatively, the bandgap obtained by fitting PL and light absorption spectra can be used to optimize the theoretically calculated SL bandgap [69]. EPM is considered a more accurate theoretical simulation method than the k·p and EFA models [27].

The famous scientist Fermi introduced a pseudopotential concept as early as 1934 when he was studying high-level electronic states. The pseudopotential approach was implemented for cases, in which the strong Coulomb potential in the vicinity of a nucleus led to the near-free electron approximation failure [24]. Among various pseudopotential techniques, the empirical pseudopotential method proposed by Phillips and Kleinman in 1959 is the most representative one [70]. Jianbai and Baldereschi developed an EPM method for simulating long-period SLs and successfully calculated the electronic structure of type I GaAs/AlGaAs SLs in 1987 [71]. They found that this method could also effectively model T2SLs. Miao et al. and Liu et al. used this technique to introduce an imaginary crystal Hamiltonian and accurately calculated the band edge structure of InAs/GaSb T2SLs in the (001) direction with different thicknesses of the InAs and GaSb layers in the 1990s [72]. They concluded that the SL conduction band edge mainly depended on the InAs layer thickness and that the valence band edge was mainly affected by the thickness of the GaSb layer, which provided a basis for the band structural design of T2SLs.

The atomic empirical pseudopotential method (AEPM) was originally proposed by Dente and Tilton in 1999 [73]. It uses the exact superposition of atomic pseudopotentials to calculate the energy and wavefunctions by solving a block equation for each atom and thus can simulate the interfacial structure very accurately. The authors successfully applied this technique to match the bandgap of the InAs/GaSb SLs reported in the study of Ram-Mohan [74] with an experimental PL spectrum. Finally, the AEPM was applied to the LWIR W-type laser AlSb/AlAsSb/InAs/GaSb/InAs SL structures with different InAs layer thicknesses. The cut-off wavelengths obtained from the corresponding PL spectra were 3.40, 3.85, and 4.40 μm, while the calculated values were 3.53, 4.01, and 4.52 μm, respectively. Hence, both datasets were almost identical within an error bar.

Subsequently, Dente and Tilton replaced the SL block function with an SL pseudopotential using pseudopotential form factors [75]. This led to a simpler version of the AEPM named the superlattice empirical pseudopotential method (SEPM), which assumed the redistribution of charges at the heterogeneous interface, making the SL components as bulky as possible and the energy as low as possible. The authors found that the bandgaps of InAs and GaSb bulk materials were equal to 0.368 and 0.8 eV at 77°K, which slightly differed from the actual values of 0.371 and 0.8 eV, respectively. This confirmed the SEPM ability to accurately model SL electronic structures. In 2013, Masur et al. compared the bandgap of a set of different InAs/GaSb SL components and found that the two-component standard SEPM method disregards the interface layers and the actual material composition of the individual layers. So they extended the two-component SEPM into a four-component model. That is, the structure uses a four-layer structure of InAs/InAsSb/GaAsSb/InAsSb including the interfaces as a full cycle [76]. As shown in Figures 7A,B, the two-component model is only suitable for SLs with layer thicknesses greater than 7 ML, whereas the four-component SEPM is significantly more accurate, and its calculated bandgaps closely match experimental values.

FIGURE 7
www.frontiersin.org

FIGURE 7. Bandgaps calculated by the (A) two-component (InAs–GaSb) and (B) four-component (InAs–InAs0.17Sb0.83–GaAs0.17Sb0.83–InAs0.17Sb0.83) SEPMs and their experimental values obtained from the results of PL measurements conducted at 10°K for T2SLs with variable layer thicknesses. The circles and squares represent thick individual layers (SL samples ≧ 7 MLs) and thin individual layers (SL samples ≦ 7 MLs), respectively. The dashed line represent that the measured and calculated energies match the ideal behavior to within an accuracy of ± kB at 77 K. The black symbols represents bisecting.

In 2001, Ongstad et al. calculated the band structure of InAs/GaSb SLs by EPM [42]. The obtained results showed that the peak wavelength of the obtained PL spectra varied from 4.2 μm for an 8 ML/8 ML sample to 3.35 μm for an 8 ML/40 ML sample, leading to a blue shift. The EPM data were in good agreement with experimental absorption and PL spectra, and the bandgap ultimately converged to a constant value. In 2003, Magri et al. adopted the EPM to change the composition of interfacial bonds by exchanging only one interfacial anion plane between Sb and As atoms and determine the relationship between the bandgap and interfacial composition [77]. They theoretically predicted the separation interface blueshift of 64 meV and the mutation interface blueshift of 95 meV (the corresponding experimental value was 70 meV). This result indicates that EPM is suitable only for the SL separation interface. The same research team published another study that compared several EPM methods (Figure 8). Here, the bandgap calculated by AEPM was much smaller than that determined by the Dente and Tilton EPM method [75], which was close to the experimental value [78]. Piquini et al. used EPM to simulate the band edges and bandgaps of n ML InAs/m ML GaSb SLs for GaSb and InAs substrates, C2v and D2d dot groups, and (001) and (110) growth directions, respectively, and compared them with experimental data [79]. The obtained results revealed that EPM could accurately predict the energy band structures of T2SLs containing thin layers.

FIGURE 8
www.frontiersin.org

FIGURE 8. Comparison of the experimental data (solid diamonds) for the (InAs)8/(GaSb)n Cut-off wavelengths with the values given by (a) the EPM of Dente et al. (solid triangles) and by the EPM for superlattices with abrupt (solid circles) and segregated (solid squares) interfaces.

In 2015, Çakan et al. used EPM to calculate the band structures of InAs, GaSb, and GaAs interfaces as well as InSb interface under strain [80]. To verify the obtained EPM data, the latter were compared with the results obtained by a hybrid density functional theory (DFT) HSEsol method [81] and experimental bandgaps. The EPM and measured bandgaps were exactly the same and equal to 0.41, 0.81, 1.51, and 0.23 eV [82]. The bandgap values calculated by the HSE method were 0.34, 0.81, 1.36, and 0.27 eV, respectively, which significantly differed from the experimental ones. In 2018, Akel et al. simulated the energy band structure of InAs/GaSb and InAs/AlSb/GaSb N-type T2SLs by EPM to determine the dependences of the SL bandgap and hh–lh splitting energy on the AlSb/GaSb and InAs layer thicknesses [69]. Note that the hh–lh splitting energy, which is larger than the bandgap, can effectively suppress the Auger recombination process. The best results with a minimum bandgap of 128 meV and hh–lh splitting energy of 194 meV were achieved for the 17 ML InAs/3 ML AlSb/6 ML GaSb SLs. Figure 9 shows the bandgaps of two T2SL structures, x ML InAs/3 ML AlSb/6 ML GaSb and x ML InAs/9 ML GaSb, which decrease with an increase in the InAs layer thickness. At x = 6, the bandgaps and the corresponding wavelengths were equal to 354 meV and 3.5 μm, and 248 meV and 5 μm, respectively, which were within the MWIR range. Hence, this study accurately calculated the bandgaps and hh–lh splitting energies of the MWIR and LWIR bands, which could be potentially used for designing photodetectors operating in these ranges. In 2021, Akel, et al. calculated the interband optical absorption of the InAs/GaSb T2SL structures [83]. The EPM method prediction showed that the SLs bandgap was underestimated about 0.4 μm, which corresponded to an uncertainty of less than 0.3 ML in the layer width, and it was in full agreement with the findings of Livneh et al. [25], Hostut and Ergun applied the method to calculate the band structure of the InAs/GaSb based T2SL [84]. The energy gap was measured as 246 meV with only ±10 meV variation compared to the experimental result at Γ point (k = 0).

FIGURE 9
www.frontiersin.org

FIGURE 9. Eg and hh–lh splitting energies plotted as functions of x (InAs ML) for the (InAs)x/(GaSb)9 and (InAs)x/(AlSb)3/(GaSb)6 T2SL structures with x varying from 6 to 19 ML. when InAs thickness is 17 ML, Eg values intersect at 124 meV with a corresponding wavelength of 10 μm within the long-wavelength range. Bandgap energy decreases below the hh-lh splitting energy as InAs ML increases.

It is noteworthy that EPM is a non-self-consistent atomic method that takes into account interfacial effects. The greatest advantage of this technique is its simple form, which is easy to implement and allows calculating the energy bands of any crystal with minimum manpower and material resources. SEPM exhibits high calculation efficiency, and the obtained wavelength tuning data for type II antimonide lasers are in good agreement with experimental values [85]. However, the EPM method cannot solve application problems in different chemical environments, especially when dealing with extremely delicate situations such as the charge transfer near the interface of a thin layer.

Empirical Tight-Bound Method

When studying SL properties, ETBM not only considers the effects of strain, interface, and antimony segregation but also uses first principles. This model can estimate specific bandgaps according to the requirements of heterostructural calculations. Thus, it may reproduce important band structural features better than the standard model and is suitable for modeling quantum well heterostructures [86].

ETBM was initially proposed by Bloch in 1929 and subsequently used to determine the periodic potentials of solids by Slater and Koster in 1954 [87, 88]. It was originally called a linear combination of atomic orbitals. In 1983, P. et al. proposed a nearest-neighbor semi-empirical tight-binding theory of sphalerite materials [89]. They assumed that the ETBM model could solve the problem of material variation at the atomic scale and retain the complete crystal and electronic symmetries of semiconductor materials. This theory was eventually applied to studying bulk materials such as GaSb and InAs.

In 2000, Klimeck et al. fitted the orbital interaction energies of nine binary compounds, including InAs and GaSb, with the sp3s* ETBM model at room temperature using literature data [89, 90] as target values [91]. The calculated energies of the lowest conduction bands of InAs and GaSb were 0.368 and 0.751 eV, while the corresponding target values were 0.370 and 0.750 eV, respectively. The calculated energies of the three highest valence bands of InAs and GaSb were −12.159, 4.126, and 4.543 eV and −12.683, 3.123, and 4.033 eV with the corresponding target values of −12.300, 4.390, and 4.630 eV and −12.000, 3.400, and 4.700 eV, respectively. These results indicate that the utilized method is very effective in predicting bandgaps. Subsequently, Wei and Razeghi modeled InAs/GaSb SLs with InSb interfaces using the more accurate sp3s* ETBM by considering the antimony bias in the InAs layer [92]. The authors computationally fixed the thickness of the GaSb layer at 40 Å (13 ML) and varied the InAs layer thickness from 40 Å (13 ML) to 66 Å (22 ML). The calculated bandgaps were compared with experimental values. The experimental data points are closely scattered around the calculated curve within an uncertainty range (Supplementary Figure S10, Supporting Information). This confirms that ETBM, which considers the interfacial and antimony segregation effects, is a reliable method for the SL design process. They subsequently studied the [(AsIn)6/AsGax1In1-x1-(SbGa)10-SbGax2In1-x2]N SLs with GaxIn1-x type mixed interfaces [93]. The maximum cutoff wavelength can be obtained when x1 = 1 and x2 = 0. The minimum cutoff wavelength can be obtained when x1 ≈ 0.5 and x2 ≈ 0.5 by the 3-dimensional view. In 2005, Wei, Hood et al. designed two types of SLs were with a cutoff wavelength around 5 μm by controls of superlattice design using ETBM: 9 ML InAs/10 ML GaSb with mixed interfaces, and 8 ML InAs/11 ML GaSb with mixed interfaces [94]. The calculation showed that the lattice constant changes do not lead to a significant change in the bandgap between 77 and 300°K.

The concept of M-type structure was initially proposed by Nguyen et al. for the InAs/GaSb and InAs/GaSb/AlSb/GaSb SLs in 2007 [95]. It maintains the type II band arrangement, which can significantly reduce the Auger transition of the electron–hole splitting band, while the AlSb barrier decreases the dark current serving as a barrier. Using this approach, the energies of the valence and conduction bands may be significantly adjusted. The authors calculated the electronic energy bands of such structures by ETBM and grew six M-type SL structures with different layer thicknesses by molecular beam epitaxy under the same conditions. The obtained PL spectra were consistent with the theoretical energy gaps, indicating that M-type structures could exhibit cut-off wavelengths greater than 11 μm (Supplementary Figure S11, Supporting Information). Nguyen et al. applied the above-mentioned sp3s* ETBM and determined the minimum and maximum values of the conduction and valence bands of the InAs/GaSb and M-type SLs by performing energy band structure calculations (Figure 10), confirming the ability of this method to adjust the band edges of the M-type structure [96]. Razeghi and Nguyen also simulated the M-type structure of InAs/GaSb SLs using ETBM and concluded that the utilized method could effectively simulate the M-type structure of type II SLs [97]. In 2012, Hoang, Chen et al. considered the GaAs and InSb interfaces to balance the compressive stresses due to mismatch induced by thick AlSb and thin InAs layers. They calculated the 50% cutoff wavelengths of the M-structure (InAs)n/GaSb/(AlSb)m/GaSb SLs using the ETBM [98]. The results show that the cutoff wavelength varies with the thickness of AlSb and InAs layers. The cutoff wavelengths of M-structure SLs are as short as 1.5 μm in the acceptable lattice mismatch range. Then, they designed one that gives a cut-off wavelength of around 2.0 μm by (InAs)6/(GaSb)1/(AlSb)5/(GaSb)1 SLs. Pour, Huang et al. designed (InAs)10/(GaSb)1/(AlSb)5/(GaSb)1 SLs with one GaAs and one InSb and (InAs)14/(GaSb)1/(AlSb)5/(GaSb)1 SLs with two InSb interfaces [99]. The cut-off wavelengths are 2.26 and 2.97 μm, respectively.

FIGURE 10
www.frontiersin.org

FIGURE 10. Energy distributions of various SL designs: (A) type II InAs/GaSb conduction band (•), (B) M-type structure conduction band (•), (C) type II InAs/GaSb valence band (light blue squares), and (D) M-type structure valence band (dark blue squares). The M-type structure exhibits a larger bandwidth.

Czuba, et al. also used the sp3s* tight-binding model to calculate the electronic structure of the inter-band cascade infrared detector with a cut-off wavelength of 10.7 μm, and obtained its effective bandgap [100]. Kato and Souma investigated the bandgap variations of (InAs)n(GaSb)n, (InAs)n(InSb)(GaSb)n and (InAs)n(GaSb)n(GaAs) T2SLs with different interface structures using the sp3s* ETBM method in 2018 [101]. They found that the shape of the band edge energy curve of (InAs)n(Gasb)n SLs is between that of (InAs)n(Insb)(Gasb)n and (InAs)n(Gasb)n(GaAs) SLs. And the intermediate property of (InAs)n(Gasb)n is thought to be due to the presence of both Ga-As and In-Sb interfacial bonds. In 2021, Zhu, et al. considered the thermal strain in the sp3s* ETBM method and calculated the energy band structure of [(AsIn)11–AsGa0.5In0.5–(SbGa)11–SbGa0.2In0.8]20 T2SLs devices. The calculated value is the bandgap of 0.2596 eV. The cutoff wavelength is 4.773 μm and the lattice mismatch is 0.46%. The measured PL spectrum and spectral responsivity verified the simulation results [102]. In the device design stage, the performance of the device can be simulated to guide the design and achieve the required performance.

In 2018, Akitaka et al. determined the parameters of the sp3d5s* ETBM model established for InAs and GaSb semiconductors and used the sp3s* and sp3d5s* models to fit their bandgaps calculated by the hybrid QSGW and ETBM approaches [103]. As seen in Figure 11, the lowest valence band energy determined the sp3s* ETBM is almost flat between points X, W, U, and K, while the sp3d5s* method successfully solves this problem and matches the results of hybrid QSGW calculations. This shows that the sp3d5s* EBTM model significantly enhances the sp3s* model and is suitable for guiding the SL design process.

FIGURE 11
www.frontiersin.org

FIGURE 11. Band structures of (A) InAs and (B) GaSb obtained by the hybrid QSGW (black) and ETBM calculations (the blue color denotes the sp3s model, and the red color represents the sp3d5s model). Whereas the lowest conduction band is fitted moderately well by the sp3s model along with the left two segments (L→Γ→X), it is almost flat between the X, W, and U, K points, reflecting the poorly described transverse mass. The sp3d5s method resolves this problem.

In 2017, Jiang, et al. reported the flexibility of adjusting the valence band level by inserting a thin (0.6 ML) InSb layer in the middle of the GaSb layer of 15 ML InAs/7 ML GaSb T2SL [104]. They used the ETBM method to calculate the tunability of the bandgap in the VLWIR SLs, and the results showed that this method extended the cut-off wavelength from 14.5 to 18.2 μm. The consistency of theoretical prediction and experimental measurement shows that this advantage can be used to achieve very long-wave infrared detection without increasing the thickness of the InAs layer. Nejad and Sheshkelani modeled the 15 ML InAs/4 ML AlSb T2SL short-wave infrared detector using the ETBM model and proposed the corresponding band structure extraction algorithm in 2020 [105]. The result showed that the ETBM simulated cut-off wavelength differs from the experiment about 30 nm, which was a good match. It is an accurate method to model the SWIR T2SL bandgap.

Compared with other methods, ETBM is advantageous for calculating the structure of the entire BZ and accurately describes the parameters of various structures from atoms to SLs. It also considers material imperfections during growth, does not require a large number of complex numerical calculations, and is relatively fast. ETBM is quite suitable for the band structure calculations of short-period SLs. However, proper parameter selection is very important for its practical implementation. The overlapping parameters have a clear and simple physical meaning, and their number increases rapidly with an increase in the number of neighbors, which negatively affects computational accuracy.

First-Principles Calculations

Density Functional Theory

The DFT was developed by Hohenberg, Kohn, and Sham in the early 1960s [106, 107]. It has been used as a standard approach for calculating the electronic structures of solids. Typically, the bandgap of an SL has a small positive value. However, DFT and other DFT-based first-principles calculation methods such as local density approximation (LDA) and generalized gradient approximation (GGA) often produce zero or even negative bandgap values due to the presence of a non-physical Coulombic self-repulsion term, which leads to a systematic bandgap underestimation. Because DFT methods are fast and widely available, many researchers have attempted to improve them to solve this problem.

In 2011, Sun, et al. calculated the energy band structure of GaAs, InSb type interfaces and their alternating superlattices based on the relaxations of both the lattice constants and atomic positions using the GGA method of DFT all-electron relativity [108, 109]. The bandgap of InSb interface superlattices calculated by virtual–crystal approximation (VCA), lattice relaxation and full relaxation are 0.2386, 0.2077, 0.1707 eV, and GaAs interface superlattices bandgap are 0.2636, 0.2511, 0.2396 eV, respectively. VCA bandgap and the fundamental optical bandgap for the alternative InSb and GaAs interfaces superlattices are 0.2243, 0.2039 eV, respectively. The results showed that VCA does not correctly describe the electronic structure of the quaternary superlattice containing the Sb chemical bond interface. It greatly overestimates the electronic bandgap. The GGA with full relaxation in lattice constants and atomic positions GGA method calculated predicts a superlattice with a much smaller bandgap. In 2012, Sun and Zheng used the method to calculate the Г-point bandgaps of InAs, GaSb, GaAs, and InSb semiconductors related to InAs/GaSb SLs [109]. The obtained results showed that for InAs and InSb with very small bandgaps, the computed values amounted to 0 eV, while the calculated bandgaps of GaAs and GaSb equal to 1.258 and 0.603 eV, respectively, were much smaller than experimental values. Thus, the electronic bandgaps calculated by the GGA method are generally lower than experimental values, and the calculation accuracy depends on a particular compound.

Caid et al. studied InAs and GaSb compounds using a full-potential linear muffin-tin orbital approach based on LDA in 2019 [110]. The obtained electronic band structure was in good agreement with the results of previous calculations, and the corresponding bandgap was almost 0 eV. Because the corresponding experimental values were 0.42 and 0.72 eV, respectively, the LDA method significantly underestimated the bandgaps of these compounds.

In 2006, Becke and Johnson proposed a very simple and effective local potential that considerably increased the bandgap computational accuracy [111]. Fabien and Peter developed a modified Becke–Johnson potential combined with a local density approximation (MBJLDA) method for bandgap calculations in 2009 [112]. As shown in Figure 12, the bandgaps calculated by the MBJLDA technique much better matched experimental values than the magnitudes computed by the LDA, hybrid functional, and GW methods [113, 114]. Subsequently, Kim et al. obtained the conduction and valence bands of five semiconductor bulk materials, including InAs and GaSb, at their Г, X, and L points by the MBJLDA, GW, and hybrid generalized methods [115]. They found that this MBJLDA method with a modified potential produced more accurate band topologies and bandgap values than the other methods. In 2016, Gmitra and Fabian used the first-principles full-potential linear augmented plane wave method to reproduce the experimental bandgaps of GaSb and InAs [116]. In that study, the TB–MBJ exchange potential was used as the exchange energy term. The calculated bandgap was 0.822 and 0.417 eV, and the corresponding experimental value was equal to 0.812 and 0.417 eV, respectively [82]. Hence, the calculated values were consistent with the experimental ones within a certain error. The authors also compared their data with the bandgap computed by Chantis et al. [117] using the GW method, which were equal to 1.16 and 0.68 eV, respectively. Thus, it can be concluded that the TB–MBJ exchange potential produces more accurate bandgap values than other DFT techniques. In recent years, Patra, et al. used LDA and MBJLDA to perform theoretical analysis and prediction of superlattice heterostructures, verifying that change in the thickness of the InSb layer could lead to the change between the direct and indirect bandgap [118]. Furthermore, the MBJLDA method is reliable and may provide a benchmark for the empirical fitting of energy band structures.

FIGURE 12
www.frontiersin.org

FIGURE 12. Theoretical and experimental bandgaps of 23 different solid compounds calculated by various DFT methods, including the LDA (plus sign), MBJLDA potentials (hollow circle) and other methods, include HSE03 (blue asterisk) and HSE06 (purple asterisk), all-electron non-self-consistent G0W0 (triangles) and self-consistent GW (crosses). For most cases, the MBJLDA potential yields bandgaps which are in good agreement with the experiment leading to typical errors of less than 10%.

In 2014, Wang and Zhang attempted to predict the electronic structure of GaSb/InAs SLs by first-principles methods [40]. They empirically corrected the s, p, and d components of the atomic pseudopotential constants by performing self-consistent calculations and adjusted the DFT electronic structures of GaSb and InAs to match experimental values or quasi-particle calculation data and thus correct the bandgap error caused by the LDA method. Subsequently, the authors studied n/8 ML SLs and compared the calculated values with experimental and EPM data to verify their accuracy. Although the calculated values are slightly higher than the experimental ones, they are still closer to the measured values than the magnitudes obtained by the EPM model (Supplementary Figure S12, Supporting Information).

Castaño-González et al. studied the electronic structure properties of GaSb using traditional DFT methods, including LDA, GGA (PBE, PBEsol, PW91, rPBE, AM05), and meta-GGA (TPSS, RTPSS, MBJ) techniques [119]. Their work provided a reference for selecting the most suitable DFT method. It was found that the bandgaps obtained by the LDA and GGA approaches were close to zero. Meanwhile, the values calculated by the MBJ, and RTPSS methods based on the meta-GGA functional and AM05 method using the GGA functional were equal to 1.015, 0.513, and 0.482 eV, respectively. These magnitudes were very close to the experimental values reported in the literature, and the MBJ and RTPSS techniques produced the lowest relative errors.

The hybrid functional incorporates a part of the Hartree–Fock exchange energy term to solve the bandgap problem. In 2010, Tomić et al. studied the hybrid exchange generalized functions B3LYP and PBE0, which contained three widely used generalized functions (LDA, GGA, and PBE) [120]. The PBE0 method overestimated the bandgaps, while the B3LYP approach produced slightly more accurate data (Supplementary Figure S13, Supporting Information). In 2016, Garza and Scuseria predicted the bandgaps of InAs and GaSb bulk materials and compared the performances of four commonly used hybrid density functionals (HSE06, B3PW91, B3LYP, and PBE0) [121]. As shown in Figure 13, these four methods somewhat overestimate the bandgap values; however, the obtained results are relatively concentrated, and the HSE06, B3PW91, and B3LYP data are closer to the experimental values within an error. Moreover, although PBE0 overestimates the bandgap, the error is systematic and can be corrected by linear fitting. Considering the cost and accuracy of this method, it may potentially have a wide application range. Yao et al. calculated the electron band structures of InAs/GaSb SLs in the (111) direction and compared the results of hybrid functional calculations with those obtained by the ordinary DFT methods [122]. The bandgaps of InAs and GaSb were also determined by the traditional PBE method and hybrid functional HSE06-PBE, HSE06-PBEsol, B3LYP as well as experimentally [123]. The obtained magnitudes were 0.00, 0.36, 0.45, 0.43, and 0.42 eV (for InAs) and 0.00, 1.02, 0.99, 0.34, and 0.81 eV (for GaSb). Hence, the traditional PBE method cannot accurately predict the bandgaps of InAs and GaSb, while the three hybrid functional methods produced relatively narrow bandgaps of these two compounds (the best results were obtained by the HSE06-PBEsol). Garwood, et al. calculated the bandgaps of InAs/GaSb T2SLs using the PBE0 hybrid functional, which was consistent with the calculated value using the 18% exact exchange function, and the deviation range from the experimental bandgap was 3–11% [124]. In 2021, Yang, et al. applied Bayesian optimization (BO) machine learning to introduce a new method of DFT with Hubbard U correction in the calculational methods [125]. According to report [126], InAs and GaSb the optimal values were shown as follows: UeffIn,p=0.5eV, UeffAs,p=7.5eV, UeffGa,p=0.8eV, UeffSb,p=6.9eV. They compared the band structures calculated by PBE+(BO) and HSE. It could be seen that PBE+U(BO) and HSE03 [127, 128] are generally in good agreement, but PBE+U(BO) underestimated the bandgap (Supplementary Figure S14, Supporting Information).

FIGURE 13
www.frontiersin.org

FIGURE 13. Calculated vs experimental bandgaps for various compounds using DFT functionals. B3LYP and B3PW91 give extremely similar bandgaps. This is most likely due to the incorporation of similar amounts of HF exchange, B3LYP and B3PW91 both have 20% full-range nonlocal exchange, while HSE06 includes 25% but only in the short-range, PBE0 includes 25% full-range HF exact exchange, which leads to an overestimation of the bandgaps.

Asadi and Nourbakhsh calculated the bandgap of InAs using the Trouiller–Martins pseudopotential combined with LDA (FHI.LDA), Hartwigsen–Goedecker–Hutter pseudopotential with parametric conservation (HGH.G0W0), and the FHI.G0W0 method including the In4d state [129]. The calculated and experimental bandgap values were 0.181, 0.845, 0.483, and 0.417 eV, respectively, indicating that the FHI.G0W0 data obtained using the multi-body perturbation theory were in good agreement with the experimental results. In 2020, the same studies the electronic structure of InAs in the DFT framework for the first time, using LDA and the norm-conserving pseudopotential that treated the In 4d electrons as valence electrons [130]. The computed bandgap without spin coupling was 0.368 eV, which underestimated the experimentally determined value of 0.420 eV by only 12.38%.

In general, the DFT approach does not rely on empirical or experimental parameters. However, DFT methods typically require performing a large number of calculations when simulating nanostructures composed of 1000s of atoms, leading to very small bandgap values. Several methods mentioned above were able to significantly improve the prediction accuracy of the electronic structures and bandgaps of various semiconductors. They included the TB–MBJ exchange potential with a localized potential, pseudopotential correction, treatment of the In 4d electrons in pseudopotentials as valence electrons, and hybridized generalization techniques. However, the problem of bandgap underestimation by DFT methods has not been completely resolved yet.

Many-Body Perturbation Theory (GW Method)

The GW method based on the many-body perturbation theory can be considered a DFT method with some modifications and perturbations, which make it suitable for calculating the excited states of a multi-body system. The widely used GW techniques include hybrid quasi-particle self-consistent GW (QSGW) and fast non-self-consistent GW (G0W0) methods. They are able to reduce the self-interaction error and solve the problem of very low theoretical bandgaps obtained by DFT methods due to the utilization of the Green function that accurately calculates the excitation spectra of quasiparticles, while the Kohn–Sham (KS) equation employed by DFT is not strictly physical.

The GW approximation was proposed in 1965 by Hedin to describe the kinetic response of a system to external perturbations [131]. The authors also derived a set of more accurate self-consistent equations containing the Green function. This step was critical for describing the quasiparticles of many-body systems and allowed a more accurate prediction of the energy band structures of solids. In 2014, Hinuma et al. performed first-principles calculations using the PBE semi-local approximation, HSE hybrid approximation, and GW approximation to study the band arrangements of various semiconductors in the sphalerite structure [123]. They compared the bandgaps computed by the GW0@PBE, GWTC–TC@HSE03 and GWГ1@HSE03 methods with the PBE, HSE03 [132], and experimental values. Figure 14 shows that 1) GWГ1@HSE most accurately reproduces the experimental data; 2) the HSE approximation can reasonably predict bandgap energies; and 3) all three GW approximations produce small deviations from experimental values. In 2014, Kotani developed the PMT–QSGW method based on the hybrid all-electron total potential method (PMT), which used both augmented plane waves and muffin-tin orbitals [133, 134]. They studied the dependence of the GaAs bandgap on the number of k points in the first BZ for self-consistent energy calculations. The GaAs bandgap at the Γ point smoothly converges with an increasing number of k points, and the value obtained for the 4 × 4 × 4 k-point mesh is only 0.1 eV higher than that determined using the 10 × 10 × 10 mesh. These results can help select an optimal number of k points to maximize the computational accuracy at limited resources.

FIGURE 14
www.frontiersin.org

FIGURE 14. Bandgaps obtained by various DFT approximations (PBE, HSE03, GW0@PBE, GWTC-TC@HSE03, and GWГ1@HSE03) and the corresponding experimental values. PBE semilocal functional has the largest deviation. The three fitting results related to GW method are consistent with experimental bandgap values.

The computational cost of the complete GW calculation is too high, which led to the development of the fast non-self-consistent GW approximation (G0W0). In 2013, Malone and Cohen used the first-principles plane-wave pseudopotential and G0W0 methods to calculate the quasiparticle band structures of InAs and GaSb semiconductors [135]. The bandgap values Г15v and Г1c of the InAs bulk material computed at the Г point by the LDA plane-wave pseudopotential, G0W0, and G0W0+SO methods, and their corresponding experimental magnitudes were equal to 0 and −0.46 eV, 0 and 0.55 eV, −0.38 and 0.42 eV, and −0.37 and 0.36 eV, respectively. The bandgap values of the GaSb bulk material were 0 and −0.13 eV, 0 and 0.94 eV, −0.73 and 0.70 eV, and −0.76 and 0.81 eV, respectively. Hence, the G0W0+SO data most closely matched the experimental values.

In 2016, Deguchi et al. proposed a hybrid QSGW method that consisted of 80% QSGW and 20% LDA and used it to calculate the band structures of various compounds, including InAs and GaSb [136]. The authors also compared the bandgaps estimated by the LDA and LDA+SO (0.00 eV), QSGW (0.8 and 1.2 eV), QSGW+SO (0.68 and 0.99 eV), QSGW80 (0.48 and 0.99 eV), and QSGW80+SO (0.36 and 0.77 eV) methods with the experimental values of 0.42 and 0.82 eV, respectively. The hybrid QSGW80+SO technique produced the closest match to the experimental data, indicating that it was one of the most accurate first-principles methods. Subsequently, Otsuka et al. applied this method for simulating a short-period InAs/GaSb infrared sensor [137]. They performed self-consistent calculations using the 4 × 4 × 3, 4 × 4 × 2, and 4 × 4 × 1 k-point mesh and calibrated the obtained bandgaps. They compared the bandgaps calculated by the QSGW, EPP [79] and ETB [92] methods at the Γ points and the corresponding experimental values derived from PL spectra [42] as functions of n (Supplementary Figure S15, Supporting Information). The bandgap energies calculated by the QSGW method are in good agreement with the ETB results. This result was in full agreement with the comparison result of the empirical sp3s* tight-binding method (TB) and eight-band k·p method, the empirical pseudopotential method (EP), the QSGW method by Kato and Souma in 2018 [101]. Although the cited research study focused on small SLs, its findings can be used as a basis for applying ETB methods to larger SLs.

Considering the importance of the InSb interfacial layer, Taghipour et al. first applied the GW approximation to study the electronic structure of long-period InAs/GaSb T2SLs in 2018 [138]. Figure 15 shows that DFT–LDA underestimates the bandgaps of all structures, while the bandgaps determined by the G0W0 method are close to the experimental values [139] and slightly differ from the GW0 data. Although the bandgap decreases with the increasing thickness of the InAs layer, the simulation results do not follow this trend due to calculation errors. The electronic bandgaps of T2SLs predicted by the GW methods are relatively consistent with the experimental data.

FIGURE 15
www.frontiersin.org

FIGURE 15. Simulated bandgaps of (4,7), (6,7), (8,8), (10,8), and (10,10) InAs/GaSb SLs in the first Brillouin zone using the DFT–LDA, G0W0+LDA, and GW0+LDA approximations with the corresponding experimental values. GW+LDA lead to results in good agreement with the experiment.

Multi-body interactions are included in the GW methods to modify the calculated LDA energy band by considering the dynamic shielding exchange, self-energy caused by Coulombic holes, local field, and dynamic shielding effect. As a result, these techniques produce more accurate bandgaps in the calculations of semiconductor electronic band structures. The number of computational steps performed during these calculations is very large; therefore, all GW methods are relatively expensive and even more expensive than the hybrid generalized function approximation. Furthermore, their SL periods are limited by a value much lower than that of the real sensor, which is a major disadvantage of these techniques.

Method Limitations and General Outlook

As the most promising materials for infrared lasers and detectors, antimonide type II SLs have very impressive application prospects. The key to the rapid advancement of type II SL technology is the realization of the importance of the design and prediction steps before material growth. In particular, designing and predicting the electronic structures and bandgap of material quickly and accurately is a time-saving and economical approach for both theoretical researchers and experimentalists. In recent decades, significant progress has been made in the prediction of class II SL energy bands. However, type II SL technology has failed to reach a theoretical prediction level comparable to that of the mature HgCdTe technology, and the calculated bandgaps of InAs/GaSb SLs are not very accurate, which significantly limits their practical use.

From the data discussed above, the following conclusions can be drawn.

• First, the quantum confinement of electrons and holes in the InAs/GaSb system is only possible for thin semiconductor layers. Second, only thin InAs and GaSb layers may produce a sufficient wavefunction overlap; however, the k·p method cannot adequately describe the electronic structure of a thin SL. Third, the k·p method does not take into account the coupling of different sub-bands, and its application to InAs/GaSb T2SLs produces a significant error in the case of strong energy band coupling. In addition, the k·p method, EPM, and ETBM are still unable to completely resolve the point defect and interfacial problems.

• EFA calculations are very complex. Although EFA considers interfacial effects, it ignores the difference between the block functions of two constituent materials in the center of the BZ at the interface, and the wavefunction is directly matched at the interface. In addition, the boundary conditions of the model may cause other uncertainties.

• The first principles use the basic laws of physics to solve Schrödinger’s equation. However, there are strong interactions between electrons, and the exact solution of the Schrodinger equation cannot be obtained for complex InAs/GaSb T2SLs. More advanced functionals, such as the hybrid density and meta-GGA ones, depend on KS orbitals. The excited states of these models are simply treated as the differences between the KS energy levels of quasi-particles, and the calculations of the excited states and optical properties produce large errors.

• The first principle calculations give very accurate results, but they require considerable execution time. They are limited to systems with a small number of atoms. In contrast, empirical methods are usually used for band structures calculation of large semiconductor heterostructures. These experimentally fitted parameters give a tremendous reduction to the computational cost. Of the various empirical methods, only EPM and ETBM are quite suitable for the band structure calculations of short-period SLs.

Interfacial effects that are generally ignored in the calculations of the InAs/GaSb SL energy bands will be the focus of future studies on bandgap prediction methods. The problems caused by electronic interactions and the first-principles underestimation of bandgap values remain important challenges in the field of energy band calculations. In the field of semiconductor materials science and engineering, it is a promising research direction that should adopt a suitable theoretical method for calculating the electronic structure bandgap and adjusting the InAs/GaSb SL components to avoid the formation of defective energy levels in the forbidden band and optimize the energy band structure. In addition, the use of experimental data for guiding the theoretical design process is also an important direction of T2SL energy band research studies. At present, experimental results are mostly used to verify the accuracy of theoretical calculations, and a large gap exists between the synergetic realization of these two aspects.

Conclusion

Type II superlattices have received a lot of attention due to their special energy band structure and excellent device performance. InAs/GaSb type II superlattices is a broken-gap band structure with the ability to adjust the positions of the conduction and valence band edges independently, providing abundant operating space for antimonide superlattices to carry out band engineering design. The energy gap falls well within the infrared regime by simply changing the components of the material and the thickness of each thin layer to be tuned, which is a key in infrared detection application. Simulating the energy band structure of InAs/GaSb T2SLs by various calculation methods, changing the bandgap, can make the theoretical calculations and experimental results complement and confirm each other. As a result, the performance of antimonide infrared detectors has reached that of HgCdTe-based systems in a relatively short period and even surpassed it in some aspects. The review covers the energy band structure of InAs/GaSb T2SLs, several commonly simulation calculation methods. And on this basis, the development direction of more suitable calculation methods oriented to interface effects and electronic interactions is proposed, which provides a reliable reference for the simulations of different superlattices structures.

Author Contributions

SF: Conceptualization, Data curation, Writing - original draft, Writing - review & editing. RH: Funding acquisition, Investigation, Methodology, Supervision, Project administration. LZ: Writing - review & editing. JG: Writing - review & editing. WL: Funding acquisition, Investigation, Methodology, Supervision, Project administration.

Funding

This work was supported by the National Key R&D Program of China under grants No. 2021YFA1400900, 2021YFA0718300, 2021YFA1400243, NSFC under grants Nos. 61835013, 61774130.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2022.822800/full#supplementary-material

Abbreviations

AEPM, atomic empirical pseudopotential method; BZ, Brillouin zone; DFT, density-functional theory; EFA, envelope function approximation; EPM, empirical pseudopotential method; ETBM, empirical tight-binding method; GGA, generalized gradient approximation; LDA, local density approximation; LWIR, long-wave infrared; MBJLDA, modified Becke–Johnson local density approximation; MWIR, mid-wave infrared; SEPM, superlattice empirical pseudopotential method; T2SL, type II superlattice.

References

1. Lv Y, Lu X, Lu Z, Li M. Review of Antimonide Infrared Detectors Development at Home and Abroad. Aviation weapon (2020) 27:1–12.

Google Scholar

2. Chang FR, Xu YQ, Niu ZC, Hao RT, Qi TT, Zhao QC, et al. High Material Quality Growth of AlInAsSb Thin Films on GaSb Substrate by Molecular Beam Epitaxy. Chin Phys B (2019) 028:445–9. doi:10.1088/1674-1056/28/2/028503

CrossRef Full Text | Google Scholar

3. Esaki L, Tsu R. Superlattice and Negative Differential Conductivity in Semiconductors. IBM J Res Dev (1970) 14:61–5. doi:10.1147/rd.141.0061

CrossRef Full Text | Google Scholar

4. Rogalski A, Martyniuk P, Kopytko M. InAs/GaSb Type-II Superlattice Infrared Detectors: Future prospect. Appl Phys Rev (2017) 4:031304. doi:10.1063/1.4999077

CrossRef Full Text | Google Scholar

5. Wei Y, Hood A, Yau H, Yazdanpanah V, Razeghi M, Tidrow MZ, et al. High-performance Type-II InAs/GaSb Superlattice Photodiodes with Cutoff Wavelength Around 7 μm. Appl Phys Lett (2005) 86:091109–0911093. doi:10.1063/1.1879113

CrossRef Full Text | Google Scholar

6. Youngdale ER, Meyer JR, Hoffman CA, Bartoli FJ, Grein CH, Young PM, et al. Auger Lifetime Enhancement in InAs-Ga1−xInxSb Superlattices. Appl Phys Lett (1994) 64:3160–2. doi:10.1063/1.111325

CrossRef Full Text | Google Scholar

7. Smith DL, Mailhiot C. Proposal for Strained Type II Superlattice Infrared Detectors. J Appl Phys (1987) 62:2545–8. doi:10.1063/1.339468

CrossRef Full Text | Google Scholar

8. Xie LL, Wang GY, Sun QY. Advances and Trends of Type-II Superlattice Infrared Detectors. In: Applied Optics and Photonics China (AOPC) Conference - Infrared Device and Infrared Technology; Beijing, CHINA. Spie-Int Soc Optical Engineering (2020).

CrossRef Full Text | Google Scholar

9. Chen YH, Tao HS, Yao DX, Liu WM. Kondo Metal and Ferrimagnetic Insulator on the Triangular Kagome Lattice. Phys Rev Lett (2012) 108:246402.1–246402.5. doi:10.1103/physrevlett.108.246402

CrossRef Full Text | Google Scholar

10. Zhang XL, Liu LF, Liu WM. Quantum Anomalous Hall Effect and Tunable Topological States in 3d Transition Metals Doped Silicene. Sci Rep (2013) 3:2908–8. doi:10.1038/srep02908

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Sai‐Halasz GA, Tsu R, Esaki L. A New Semiconductor Superlattice. Appl Phys Lett (1977) 30:651–3.

Google Scholar

12. Özcelik VO, Azadani JG, Yang C, Koester SJ, Low T. Band Alignment of Two-Dimensional Semiconductors for Designing Heterostructures with Momentum Space Matching. Phys Rev B (2016) 94:035125.

Google Scholar

13. Becer Z, Bennecer A, Sengouga N. Modeling Energy Bands in Type II Superlattices. Crystals (2019) 9:629. doi:10.3390/cryst9120629

CrossRef Full Text | Google Scholar

14. Thi YL, Kamakura Y, Mori N. Simulation of Dark Current Characteristics of Type-II InAs/GaSb Superlattice Mid-wavelength Infrared P-I-N Photodetector. Jpn J Appl Phys (2019) 58:044002. doi:10.7567/1347-4065/ab03ca

CrossRef Full Text | Google Scholar

15. Jiang M, Xiao H, Peng S, Qiao L, Yang G, Liu Z, et al. First-Principles Study of Point Defects in GaAs/AlAs Superlattice: the Phase Stability and the Effects on the Band Structure and Carrier Mobility. Nanoscale Res Lett (2020) 13:301. doi:10.1186/s11671-018-2719-7

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Klipstein PC, Benny Y, Cohen Y, Fraenkel N, Fraenkel R, Gliksman S. Type II Superlattice Detectors at SCD. In: Conference on Infrared Technology and Applications XLVII. Spie-Int Soc Optical EngineeringElectr Network (2021).

CrossRef Full Text | Google Scholar

17. Delmas M, Liang BL, Huffaker DL. A Comprehensive Set of Simulation Tools to Model and Design High Performance Type-II InAs/GaSb Superlattice Infrared Detectors. In: Conference on Quantum Sensing and Nano Electronics and Photonics XVI; San Francisco, CA. Spie-Int Soc Optical Engineering (2019).

CrossRef Full Text | Google Scholar

18. Rogalski A, Kopytko M, Martyniuk P. InAs/GaSb Type-II Superlattice Infrared Detectors: Three Decades of Development. In: Conference on Infrared Technology and Applications XLIII; Anaheim, CA (2017).

CrossRef Full Text | Google Scholar

19. Bardeen J. An Improved Calculation of the Energies of Metallic Li and Na. J Chem Phys (1938) 6:367–71. doi:10.1063/1.1750270

CrossRef Full Text | Google Scholar

20. Seitz F, Johnson RP. Modern Theory of Solids. I. J Appl Phys (1937) 8:84–97. doi:10.1063/1.1710273

CrossRef Full Text | Google Scholar

21. Jaros M. Wave Mechanics Applied to Semiconductor Heterostructures. Optica Acta Int J Opt (1992) 38:1211–2.

Google Scholar

22. Smith DL, Mailhiot C. Theory of Semiconductor Superlattice Electronic Structure. Rev Mod Phys (1990) 62:173–234. doi:10.1103/revmodphys.62.173

CrossRef Full Text | Google Scholar

23. Read WT, Shockley W. Dislocation Models of Crystal Grain Boundaries. Phys Rev (1950) 78:275–89. doi:10.1103/physrev.78.275

CrossRef Full Text | Google Scholar

24. Klipstein PC. Operator Ordering and Interface-Band Mixing in the Kane-like Hamiltonian of Lattice-Matched Semiconductor Superlattices with Abrupt Interfaces. Phys Rev (2010) 81:235314.1–235314.21. doi:10.1103/physrevb.81.235314

CrossRef Full Text | Google Scholar

25. Livneh Y, Klipstein PC, Klin O, Snapi N, Grossman S, Glozman A, et al. k·pmodel for the Energy Dispersions and Absorption Spectra of InAs/GaSb Type-II Superlattices. Phys Rev B (2012) 86:235311. doi:10.1103/physrevb.86.235311

CrossRef Full Text | Google Scholar

26. Lawaetz P. Valence-band Parameters in Cubic Semiconductors. Phys Rev B (1971) 4:3460–7. doi:10.1103/physrevb.4.3460

CrossRef Full Text | Google Scholar

27. Rejeb SB, Debbichi M, Said M, Gassenq A, Tournié E, Christol P. Modelling of an InAs/GaSb/InSb Short-Period Superlattice Laser Diode for Mid-infrared Emission by the k.P Method. J Phys D: Appl Phys (2010) 43:325102. doi:10.1088/0022-3727/43/32/325102

CrossRef Full Text | Google Scholar

28. Qiao P-F, Mou S, Chuang SL. Electronic Band Structures and Optical Properties of Type-II Superlattice Photodetectors with Interfacial Effect. Opt Express (2012) 20:2319–34. doi:10.1364/oe.20.002319

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Klipstein PC, Livneh Y, Glozman A, Grossman S, Klin O, Snapi N, et al. Modeling InAs/GaSb and InAs/InAsSb Superlattice Infrared Detectors. J Elec Materi (2014) 43:2984–90. doi:10.1007/s11664-014-3169-3

CrossRef Full Text | Google Scholar

30. Klipstein PC, Livneh Y, Klin O, Grossman S, Snapi N, Glozman A, et al. A K·p Model of InAs/GaSb Type II Superlattice Infrared Detectors. Infrared Phys Technology (2013) 59:53–9. doi:10.1016/j.infrared.2012.12.009

CrossRef Full Text | Google Scholar

31. Delmas M, Kwan DCM, Debnath MC, Liang BL, Huffaker DL. Flexibility of Ga-containing Type-II Superlattice for Long-Wavelength Infrared Detection. J Phys D: Appl Phys (2019) 52:475102. doi:10.1088/1361-6463/ab3b6a

CrossRef Full Text | Google Scholar

32. Vinter B. Auger Recombination in Narrow-gap Semiconductor Superlattices. Phys Rev B (2002) 66:045324 1–9. doi:10.1103/physrevb.66.045324

CrossRef Full Text | Google Scholar

33. Taalat R, Rodriguez JB, Delmas M, Christol P. Influence of the Period Thickness and Composition on the Electro-Optical Properties of Type-II InAs/GaSb Midwave Infrared Superlattice Photodetectors. J Phys D Appl Phys (2014) 47:5101. doi:10.1088/0022-3727/47/1/015101

CrossRef Full Text | Google Scholar

34. Taalat R, Rodriguez J-B, Delmas M, Christol P. Influence of the Period Thickness and Composition on the Electro-Optical Properties of Type-II InAs/GaSb Midwave Infrared Superlattice Photodetectors. J Phys D: Appl Phys (2013) 47:015101. doi:10.1088/0022-3727/47/1/015101

CrossRef Full Text | Google Scholar

35. Imbert J, Trinite V, Derelle S, Jaeck J, Giard E, Delmas M, et al. Electronic Structure of InAs/GaSb Superlattice for the Modelling of MWIR Pin Photodiode. Infrared Phys Technology (2015) 70:81–6. doi:10.1016/j.infrared.2014.09.035

CrossRef Full Text | Google Scholar

36. Machowska-Podsiadlo E, Bugajski M. Influence of Various Parameters and Phenomena on the Absorption Edge of InAs/GaSb Superlattices. Superlattices Microstructures (2018) 125:214–9. doi:10.1016/j.spmi.2018.06.004

CrossRef Full Text | Google Scholar

37. Szmulowicz F, Haugan H, Brown GJ. Effect of Interfaces and the Spin-Orbit Band on the Band Gaps of InAs/GaSb Superlattices beyond the Standard Envelope-Function Approximation. Phys Rev B (2004) 69:155321. doi:10.1103/PhysRevB.69.155321

CrossRef Full Text | Google Scholar

38. Jeffrey S. Modeling of 6.1Å Family Strained Layer Superlattice Structures (2019).

Google Scholar

39. Du Y-n., Xu Y, Song G-f. Theoretical Analysis on the Energy Band Properties of N- and M-Structure Type-II Superlattices. Superlattices and Microstructures (2020) 145:106590. doi:10.1016/j.spmi.2020.106590

CrossRef Full Text | Google Scholar

40. Wang J, Zhang Y. Band-gap Corrected Density Functional Theory Calculations for InAs/GaSb Type II Superlattices. J Appl Phys (2014) 116:214301. doi:10.1063/1.4903063

CrossRef Full Text | Google Scholar

41. Magri R, Zunger A. Effects of Interfacial Atomic Segregation and Intermixing on the Electronic Properties of InAs/GaSb Superlattices. Phys Rev B (2002) 65:165302. doi:10.1103/physrevb.65.165302

CrossRef Full Text | Google Scholar

42. Ongstad AP, Kaspi R, Moeller CE, Tilton ML, Gianardi DM, Chavez JR, et al. Spectral Blueshift and Improved Luminescent Properties with Increasing GaSb Layer Thickness in InAs-GaSb Type-II Superlattices. J Appl Phys (2001) 89:2185–8. doi:10.1063/1.1337918

CrossRef Full Text | Google Scholar

43. Cui S-N, Jiang D-W, Sun J, Jia Q-X, Li N, Zhang X, et al. Investigation of Active-Region Doping on InAs/GaSb Long Wave Infrared Detectors*. Chin Phys. B (2020) 29:048502. doi:10.1088/1674-1056/ab773c

CrossRef Full Text | Google Scholar

44. Mukherjee S, Singh A, Bodhankar A, Muralidharan B. Carrier Localization and Miniband Modeling of InAs/GaSb Based Type-II Superlattice Infrared Detectors. J Phys D: Appl Phys (2021) 54:345104. doi:10.1088/1361-6463/ac0702

CrossRef Full Text | Google Scholar

45. Kim HS. Optical and Electrical Study of InAs/GaSb Type II Strained Layer Superlattice for Mid-wave Infrared Detector. J Korean Phys Soc (2019) 74:358–62. doi:10.3938/jkps.74.358

CrossRef Full Text | Google Scholar

46. Li LL, Xu W, Peeters FM. Intrinsic Optical Anisotropy of [001]-grown Short-Period InAs/GaSb Superlattices. Phys Rev B (2010) 82:235422. doi:10.1103/physrevb.82.235422

CrossRef Full Text | Google Scholar

47. Dong HM, Li LL, Xu W, Han K. Effect of Microscopic Interface Asymmetry on Optical Properties of Short-Period InAs/GaSb Type-II Superlattices. Thin Solid Films (2015) 589:388–95. doi:10.1016/j.tsf.2015.05.066

CrossRef Full Text | Google Scholar

48. Le YT, Kamakura Y, Mori N. Evaluation of the Optical Characteristics of Type-II InAs/GaSb Superlattice Infrared P-I-N Photodetector. Jpn J Appl Phys (2019) 58.

Google Scholar

49. Jiang Y, Ma X, Xu Y, Song G. Finite Difference Method for Analyzing Band Structure in Semiconductor Heterostructures without Spurious Solutions. J Appl Phys (2014) 116:173702.

CrossRef Full Text | Google Scholar

50. Shams MIB, Xie Y, Lu Y, Fay P. An Accurate Interband Tunneling Model for InAs/GaSb Heterostructure Devices. Phys Status Solidi C (2013) 10:740–3. doi:10.1002/pssc.201200624

CrossRef Full Text | Google Scholar

51. Zunger A. On the Farsightedness (Hyperopia) of the Standard K . P Model. Phys Stat Sol (A) (2002) 190:467–75. doi:10.1002/1521-396x(200204)190:2<467::aid-pssa467>3.0.co;2-4

CrossRef Full Text | Google Scholar

52. Wang L-W, Wei S-H, Mattila T, Zunger A, Vurgaftman I, Meyer JR. Multiband Coupling and Electronic Structure of(InAs)n/(GaSb)nsuperlattices. Phys Rev B (1999) 60:5590–6. doi:10.1103/physrevb.60.5590

CrossRef Full Text | Google Scholar

53. Kitchin MR, Hagon JP, Jaros M. Models of GaSb/InAs Type-II Infrared Detectors at Very Long Wavelengths: Band Offsets and Interface Bonds. Semicond Sci Technol (2003) 18:225–33. doi:10.1088/0268-1242/18/4/306

CrossRef Full Text | Google Scholar

54. Rodriguez JB, Christol P, Chevrier F, Joullié A. Optical Characterization of Symmetric InAs/GaSb Superlattices for Detection in the 3-5μm Spectral Region. Physica E: Low-dimensional Syst Nanostructures (2005) 28:128–33. doi:10.1016/j.physe.2005.02.007

CrossRef Full Text | Google Scholar

55. Szmulowicz F, Haugan H, Brown G, Mahalingam K, Ullrich B, Munshi S, et al. Interfaces as Design Tools for Short-Period InAs/GaSb Type-II Superlattices for Mid-infrared Detectors. Opto-Electronics Rev (2006) 14:69–75. doi:10.2478/s11772-006-0010-4

CrossRef Full Text | Google Scholar

56. Haugan HJ, Szmulowicz F, Brown GJ, Ullrich B, R Munshi S, Grazulis L, et al. Pushing the Envelope to the Maximum: Short-Period InAs/GaSb Type-II Superlattices for Mid-infrared Detectors. Physica E: Low-dimensional Syst Nanostructures (2006) 32:289–92. doi:10.1016/j.physe.2005.12.072

CrossRef Full Text | Google Scholar

57. Debbichi M, Ben Rejeb S, Debbichi L, Said M. Interfaces as Design Tools for the InAs/GaSb/InSb Short-Period Superlattice for Mid-infrared Emission. Semicond Sci Technol (2011) 26:095010. doi:10.1088/0268-1242/26/9/095010

CrossRef Full Text | Google Scholar

58. Zhou Y, Chen J-X, Xu Q-Q, Xu Z-C, Jin C, Xu J-J, et al. Long Wavelength Infrared Detector Based on Type-II InAs/GaSb Superlattice. J Infrared Millimeter Waves (2013) 32:210–3+224. doi:10.3724/sp.j.1010.2013.00210

CrossRef Full Text | Google Scholar

59. Boutramine A, Nafidi A, Barkissy D, Hannour A, Massaq M, Chaib H. Electronic Bands Structure and gap in Mid-infrared Detector InAs/GaSb Type II Nanostructure Superlattice. J Eng Res Appl (2014) 2014:132–5.

Google Scholar

60. Li LL, Xu W, Zhang J, Shi YL. Midinfrared Absorption by InAs/GaSb Type-II Superlattices. J Appl Phys (2009) 105:013115. doi:10.1063/1.3058692

CrossRef Full Text | Google Scholar

61. Cervera C, Rodriguez JB, Perez JP, Aït-Kaci H, Chaghi R, Konczewicz L, et al. Unambiguous Determination of Carrier Concentration and Mobility for InAs/GaSb Superlattice Photodiode Optimization. J Appl Phys (2009) 106:033709. doi:10.1063/1.3191175

CrossRef Full Text | Google Scholar

62. Benchtaber N, Nafidi A, Barkissy D, Boutramine A, Benaadad M, Melkoud S, et al. Theoretical Electronic Band Structures and Transport in InAs/GaSb Type II Nanostructure Superlattice for Medium Infrared Detection. Mater Today Proc (2020) 22:41–4. doi:10.1016/j.matpr.2019.08.069

CrossRef Full Text | Google Scholar

63. Haugan HJ, Elhamri S, Ullrich B, Szmulowicz F, Brown GJ, Mitchel WC. Optimizing Residual Carriers in Undoped InAs/GaSb Superlattices for High Operating Temperature Mid-infrared Detectors. J Cryst Growth (2009) 311:1897–900. doi:10.1016/j.jcrysgro.2008.09.141

CrossRef Full Text | Google Scholar

64. Boutramine A, Nafidi A, arkissy DB, Hannour A, Elanique A, Gouti TE. Application of the Transition Semiconductor to Semimetal in Type II Nanostructure Superlattice for Mid-infrared Optoelectronic Devices. Appl Phys A (2016) 122:1–6. doi:10.1007/s00339-016-9911-3

CrossRef Full Text | Google Scholar

65. Boutramine A, Nafidi A, Barkissy D, El-Frikhe E-S, Charifi H, Elanique A, et al. Electronic Band Structure and Shubnikov-De Haas Effect in Two-Dimensional Semimetallic InAs/GaSb Nanostructure Superlattice. Appl Phys A (2016) 122:70. doi:10.1007/s00339-015-9561-x

CrossRef Full Text | Google Scholar

66. Boutramine A, Nafidi A, Barkissy D, Bellioua M, Khalal A. Correlation between Electronic Bands Structure and Magneto-Transport Properties of Nanostructure Type II Superlattice for Terahertz Detection. Superlattices and Microstructures (2019) 127:151–6. doi:10.1016/j.spmi.2017.12.036

CrossRef Full Text | Google Scholar

67. Hostut M, Tansel T, Kilic A, Akin T, Ergun Y. The Detailed Analysis of Wavefunction Overlaps for InAs/AlSb/GaSb Based N-Structure Type-II SL Pin Photodetectors. Phys Scr (2019) 94:075007. doi:10.1088/1402-4896/ab13f9

CrossRef Full Text | Google Scholar

68. Wood DM, Zunger A. Successes and Failures of Thek⋅pmethod: A Direct Assessment for GaAs/AlAs Quantum Structures. Phys Rev B (1996) 53:7949–63. doi:10.1103/physrevb.53.7949

CrossRef Full Text | Google Scholar

69. Akel K, Hostut M, Tansel T, Ergun Y. Large Hh-Lh Splitting Energy for InAs/AlSb/GaSb Based N-Structure Photodetectors. J Appl Phys (2018) 123:025703. doi:10.1063/1.4999632

CrossRef Full Text | Google Scholar

70. Phillips JC, Kleinman L. New Method for Calculating Wave Functions in Crystals and Molecules. Phys Rev (1959) 116:287–94. doi:10.1103/physrev.116.287

CrossRef Full Text | Google Scholar

71. Xia J, Baldereschi A. Pseudopotential Calculation for the Electronic Structures of Superlattices. Chin J Semiconductors (1987) 574–84.

Google Scholar

72. Liu S, Long F, Mei F. Pseudopotential Calculation of the Band Edge Structures of the Type-Ⅱ InAs/GaSb(001) Superlattices. J JLANGXI NORMAL UNIVERSITY (1996) 83–6.

Google Scholar

73. Dente GC, Tilton ML. Pseudopotential Methods for Superlattices: Applications to Mid-infrared Semiconductor Lasers. J Appl Phys (1999) 86:1420–9. doi:10.1063/1.370905

CrossRef Full Text | Google Scholar

74. Ram-Mohan . Private Communication (2022).

Google Scholar

75. Dente GC, Tilton ML. Comparing Pseudopotential Predictions for InAs/GaSb Superlattices. Phys Rev B (2002) 66:165307. doi:10.1103/physrevb.66.165307

CrossRef Full Text | Google Scholar

76. Masur J-M, Rehm R, Schmitz J, Kirste L, Walther M. Four-component Superlattice Empirical Pseudopotential Method for InAs/GaSb Superlattices. Infrared Phys Technology (2013) 61:129–33. doi:10.1016/j.infrared.2013.07.014

CrossRef Full Text | Google Scholar

77. Magri R, Zunger A. Theory of Optical Properties of 6.1 Å III-V Superlattices: The Role of the Interfaces. J Vac Sci Technol B (2003) 21:1896–902. doi:10.1116/1.1589519

CrossRef Full Text | Google Scholar

78. Magri R, Zunger A. Theory of Optical Properties of Segregated InAs/GaSb Superlattices. IEE Proc Optoelectron (2003) 150:409–14. doi:10.1049/ip-opt:20030843

CrossRef Full Text | Google Scholar

79. Piquini P, Zunger A, Magri R. Pseudopotential Calculations of Band Gaps and Band Edges of Short-period(InAs)n∕(GaSb)msuperlattices with Different Substrates, Layer Orientations, and Interfacial Bonds. Phys Rev B (2008) 77:115314. doi:10.1103/physrevb.77.115314

CrossRef Full Text | Google Scholar

80. Çakan A, Sevik C, Bulutay C. Strained Band Edge Characteristics from Hybrid Density Functional Theory and Empirical Pseudopotentials: GaAs, GaSb, InAs and InSb. J Phys D: Appl Phys (2016) 49:085104. doi:10.1088/0022-3727/49/8/085104

CrossRef Full Text | Google Scholar

81. Laurids S, Judith H, Georg K. Improved Hybrid Functional for Solids: the HSEsol Functional. J Chem Phys (2011) 134:024116. doi:10.1063/1.3524336

PubMed Abstract | CrossRef Full Text | Google Scholar

82. Vurgaftman I, Meyer JR, Ram-Mohan LR. Band Parameters for III-V Compound Semiconductors and Their Alloys. J Appl Phys (2001) 89:5815–75. doi:10.1063/1.1368156

CrossRef Full Text | Google Scholar

83. Akel K, Hoştut M, Akın T, Ergün Y. Interband Optical Absorption Obtained by Pseudopotential Method for Type-II InAs/GaSb SL Photodetectors. J Phys D: Appl Phys (2021) 54:195103. doi:10.1088/1361-6463/abe3ae

CrossRef Full Text | Google Scholar

84. Hostut M, Ergun Y. Quantum Efficiency Contributions for Type-II InAs/GaSb SL Photodetectors. Physica E: Low-dimensional Syst Nanostructures (2021) 130:114721. doi:10.1016/j.physe.2021.114721

CrossRef Full Text | Google Scholar

85. Dente GC, Tilton ML, Ongstad AP, Kaspi R. Wavelength Tuning Predictions and Experiments for Type II Antimonide Lasers. J Appl Phys (2008) 103:041122. doi:10.1063/1.2836785

CrossRef Full Text | Google Scholar

86. Boykin TB, Klimeck G, Bowen RC, Lake R. Effective-mass Reproducibility of the Nearest-Neighborsp3s*models: Analytic Results. Phys Rev B (1997) 56:4102–7. doi:10.1103/physrevb.56.4102

CrossRef Full Text | Google Scholar

87. Bloch F. Über die Quantenmechanik der Elektronen in Kristallgittern. Z Physik (1929) 52:555–600. doi:10.1007/bf01339455

CrossRef Full Text | Google Scholar

88. Slater JC, Koster GF. Simplified LCAO Method for the Periodic Potential Problem. Phys Rev (1954) 94:1498–524. doi:10.1103/physrev.94.1498

CrossRef Full Text | Google Scholar

89. Vogl P, Hjalmarson HP, Dow JD. A Semi-empirical Tight-Binding Theory of the Electronic Structure of Semiconductors†. J Phys Chem Sol (1983) 44:365–78. doi:10.1016/0022-3697(83)90064-1

CrossRef Full Text | Google Scholar

90. Jancu J-M, Scholz R, Beltram F, Bassani F. Empiricalspds*tight-binding Calculation for Cubic Semiconductors: General Method and Material Parameters. Phys Rev B (1998) 57:6493–507. doi:10.1103/physrevb.57.6493

CrossRef Full Text | Google Scholar

91. Klimeck G, Chris Bowen R, Boykin TB, Cwik TA. sp3s*Tight-binding Parameters for Transport Simulations in Compound Semiconductors. Superlattices and Microstructures (2000) 27:519–24. doi:10.1006/spmi.2000.0862

CrossRef Full Text | Google Scholar

92. Wei Y, Razeghi M. Modeling of Type-II InAs/GaSb Superlattices Using an Empirical Tight-Binding Method and Interface Engineering. Phys Rev B (2004) 69:085316. doi:10.1103/physrevb.69.085316

CrossRef Full Text | Google Scholar

93. Wei Y, Razeghi M, Brown GJ, Tidrow MZ. Modeling Type-II InAs/GaSb Superlattices Using Empirical Tight-Binding Method: New Aspects. In: Quantum Sensing and Nanophotonic Devices; 26-29 JANUARY 2004; San Jose, CA, US. International Society for Optics and Photonics (2004). p. 301–8. doi:10.1117/12.528297

CrossRef Full Text | Google Scholar

94. Wei Y, Hood A, Yau H, Gin A, Razeghi M, Tidrow MZ, et al. Uncooled Operation of Type-II InAs∕GaSb Superlattice Photodiodes in the Midwavelength Infrared Range. Appl Phys Lett (2005) 86:233106. doi:10.1063/1.1947908

CrossRef Full Text | Google Scholar

95. Nguyen B, Razeghi M, Nathan V, Brown GJ. Type-II M Structure Photodiodes: an Alternative Material Design for Mid-wave to Long Wavelength Infrared Regimes. In: Quantum Sensing and Nanophotonic Devices IV; 22-25. Jan 2007; San Jose, CA, US. International Society for Optics and Photonics (2007). p. 64790S.

CrossRef Full Text | Google Scholar

96. Nguyen B-M, Hoffman D, Delaunay P-Y, Huang EK-W, Razeghi M, Pellegrino J. Band Edge Tunability of M-Structure for Heterojunction Design in Sb Based Type II Superlattice Photodiodes. Appl Phys Lett (2008) 93:163502. doi:10.1063/1.3005196

CrossRef Full Text | Google Scholar

97. Razeghi M, Nguyen B-M. Band gap Tunability of Type II Antimonide-Based Superlattices. Phys Proced (2010) 3:1207–12. doi:10.1016/j.phpro.2010.01.164

CrossRef Full Text | Google Scholar

98. Hoang AM, Chen G, Haddadi A, Abdollahi Pour S, Razeghi M. Demonstration of Shortwavelength Infrared Photodiodes Based on Type-II InAs/GaSb/AlSb Superlattices. Appl Phys Lett (2012) 100:211101. doi:10.1063/1.4720094

CrossRef Full Text | Google Scholar

99. Pour SA, Huang EK, Chen G, Haddadi A, Nguyen B-M, Razeghi M. High Operating Temperature Midwave Infrared Photodiodes and Focal Plane Arrays Based on Type-II InAs/GaSb Superlattices. Appl Phys Lett (2011) 98:143501. doi:10.1063/1.3573867

CrossRef Full Text | Google Scholar

100. Czuba K, Sankowska I, Jasik A, Papis-Polakowska E, Kaniewski J, The Analysis of Functional Regions in the Long-Wavelength-Infrared Interband cascade Photodetector, Electro-Optical and Infrared Systems: Technology and Applications XIV, Int Soc Opt Photon (2017) 10433:104330K.

Google Scholar

101. Kato T, Souma S. sp3s∗ Tight-Binding Calculation of Band Edges and Effective Masses of InAs/GaSb Superlattices with Different Interface Structures. Superlattices and Microstructures (2018) 122:492–500. doi:10.1016/j.spmi.2018.06.060

CrossRef Full Text | Google Scholar

102. Zhu X, Jie W, Lyu Y, Peng Z, Cao X, Li M, et al. Band Structure Calculation in Type-II InAs/GaSb Superlattices Detector with thermal Strain Using an Empirical Tight-Binding Method. Appl Phys A (2021) 127:1–7. doi:10.1007/s00339-020-04158-1

CrossRef Full Text | Google Scholar

103. Sawamura A, Otsuka J, Kato T, Kotani T, Souma S. Nearest-neighbor Sp3d5s* Tight-Binding Parameters Based on the Hybrid Quasi-Particle Self-Consistent GW Method Verified by Modeling of Type-II Superlattices. Opt Mater Express (2018) 8:1569–84. doi:10.1364/ome.8.001569

CrossRef Full Text | Google Scholar

104. Jiang D, Han X, Hao H, Sun Y, Jiang Z, Lv Y, et al. Significantly Extended Cutoff Wavelength of Very Long-Wave Infrared Detectors Based on InAs/GaSb/InSb/GaSb Superlattices. Appl Phys Lett (2017) 111:161101. doi:10.1063/1.4998502

CrossRef Full Text | Google Scholar

105. Nejad SM, Sheshkelani NR. The Modeling of a SWIR Type-II InAs/AlSb Superlattice Using an ETBM and Interface Engineering. Superlattices and Microstructures (2020) 143:106523. doi:10.1016/j.spmi.2020.106523

CrossRef Full Text | Google Scholar

106. Kohn W, Sham LJ. Self-consistent Equations Including Exchange and Correlation Effects. Phys Rev (1965) 140:A1133–A1138. doi:10.1103/physrev.140.a1133

CrossRef Full Text | Google Scholar

107. Hohenberg P, Kohn W. Inhomogeneous Electron Gas. Phys Rev (1964) 136:B864–B871. doi:10.1103/physrev.136.b864

CrossRef Full Text | Google Scholar

108. Sun W-F, Li M-C, Zhao L-C. First-principles Study of Interface Relaxation Effect on Interface and Electronic Structures of InAs/GaSb Superlattices with Different Interface Types. Superlattices and Microstructures (2011) 49:81–90. doi:10.1016/j.spmi.2010.11.006

CrossRef Full Text | Google Scholar

109. Sun W-F, Zheng X-X. First-principles Study of Interface Relaxation Effect on Influence Structure, Band Structure and Optical Property of InAs/GaSb Superlattices. Acta Phys (2012) 61:461–70. doi:10.7498/aps.61.117301

CrossRef Full Text | Google Scholar

110 Caid M, Rached D, Cheref O, Righi H, Rached H, Benalia S, et al. Full Potential Study of the Structural, Electronic and Optical Properties of (InAs)m/(GaSb)n Superlattices. Comput Condensed Matter (2019) 21:e00394. doi:10.1016/j.cocom.2019.e00394

CrossRef Full Text | Google Scholar

111 Becke AD, Johnson ER. A Simple Effective Potential for Exchange. J Chem Phys (2006) 124:221101. doi:10.1063/1.2213970

PubMed Abstract | CrossRef Full Text | Google Scholar

112 Tran F, Blaha P. Accurate Band Gaps of Semiconductors and Insulators with a Semilocal Exchange-Correlation Potential. Phys Rev Lett (2009) 102:226401. doi:10.1103/physrevlett.102.226401

PubMed Abstract | CrossRef Full Text | Google Scholar

113 Marsman M, Paier J, Stroppa A, Kresse G. Hybrid Functionals Applied to Extended Systems. J Phys Condens Matter (2008) 20:064201. doi:10.1088/0953-8984/20/6/064201

PubMed Abstract | CrossRef Full Text | Google Scholar

114 Shishkin M, Marsman M, Kresse G. Accurate Quasiparticle Spectra from Self-ConsistentGWCalculations with Vertex Corrections. Phys Rev Lett (2007) 99:246403. doi:10.1103/physrevlett.99.246403

PubMed Abstract | CrossRef Full Text | Google Scholar

115 Kim Y-S, Marsman M, Kresse G, Tran F, Blaha P. Towards Efficient Band Structure and Effective Mass Calculations for III-V Direct Band-gap Semiconductors. Phys Rev B (2010) 82:205212. doi:10.1103/physrevb.82.205212

CrossRef Full Text | Google Scholar

116 Gmitra M, Fabian J. First-principles Studies of Orbital and Spin-Orbit Properties of GaAs, GaSb, InAs, and InSb Zinc-Blende and Wurtzite Semiconductors. Phys Rev B (2016) 94:165202. doi:10.1103/physrevb.94.165202

CrossRef Full Text | Google Scholar

117 Chantis AN, van Schilfgaarde M, Kotani T. Ab Initio prediction of Conduction Band Spin Splitting in Zinc Blende Semiconductors. Phys Rev Lett (2006) 96:086405. doi:10.1103/PhysRevLett.96.086405

PubMed Abstract | CrossRef Full Text | Google Scholar

118 Patra A, Chakraborty M, Roy A. Electronic Band Structure Engineering in InAs/InSbAs and InSb/InSbAs Superlattice Heterostructures. J Appl Phys (2019) 125:025704. doi:10.1063/1.5056196

CrossRef Full Text | Google Scholar

119 Castaño-González EE, Seña N, Mendoza-Estrada V, González-Hernández R, Dussan A, Mesa F. First-principles Calculations of the Electronic and Structural Properties of GaSb. Semiconductors (2016) 50:1280–6.

Google Scholar

120 Tomić S, Harrison NM, Caldas Ml., Studart N. Electronic Structure of III–V’s Semiconductors from B3LYP and PBE0 Functionals. In: 29th International Conference on the Physics of Semiconductors. College Park, Maryland, US: American Institute of Physics (2010). p. 65–6.

Google Scholar

121 Garza AJ, Scuseria GE. Predicting Band Gaps with Hybrid Density Functionals. J Phys Chem Lett (2016) 7:4165–70. doi:10.1021/acs.jpclett.6b01807

PubMed Abstract | CrossRef Full Text | Google Scholar

122 Yao LC, Zhou XH, Chen XS. Hybrid Functional Calculation of Electronic Structure of InAs/GaSb Superlattice in (111) Orientation. J Infrared Millim W (2016) 35:646–51.

Google Scholar

123 Hinuma Y, Grüneis A, Kresse G, Oba F. Band Alignment of Semiconductors from Density-Functional Theory and many-body Perturbation Theory. Phys Rev B (2014) 90:155405. doi:10.1103/physrevb.90.155405

CrossRef Full Text | Google Scholar

124 Garwood T, Modine NA, Krishna S. Electronic Structure Modeling of InAs/GaSb Superlattices with Hybrid Density Functional Theory. Infrared Phys Technology (2017) 81:27–31. doi:10.1016/j.infrared.2016.12.007

CrossRef Full Text | Google Scholar

125 Yang S, Dardzinski D, Hwang A, Pikulin DI, Winkler GW, Marom N. First Principles Feasibility Assessment of a Topological Insulator at the InAs/GaSb Interface. Phys Rev Mater (2021) 5:10. doi:10.1103/physrevmaterials.5.084204

CrossRef Full Text | Google Scholar

126 Yu M, Yang S, Wu C, Marom N. Machine Learning the Hubbard U Parameter in DFT+U Using Bayesian Optimization. NPJ COMPUTATIONAL MATERIALS (2020) 6. doi:10.1038/s41524-020-00446-9

CrossRef Full Text | Google Scholar

127 Heyd J, Scuseria GE, Ernzerhof M. Hybrid Functionals Based on a Screened Coulomb Potential. J Chem PhysicsJournal Chem Phys (20032006) 118118:82078207–8215. doi:10.1063/1.1564060

CrossRef Full Text | Google Scholar

128 Krukau AV, Vydrov OA, Izmaylov AF, Scuseria GE. Influence of the Exchange Screening Parameter on the Performance of Screened Hybrid Functionals. J Chem Phys (2006) 125:224106. doi:10.1063/1.2404663

PubMed Abstract | CrossRef Full Text | Google Scholar

129 Asadi Y, Nourbakhsh Z. First Principle Characterization of Structural, Electronic, Mechanical, Thermodynamic, Linear and Nonlinear Optical Properties of Zinc Blende InAs, InSb and Their InAsxSb1-X Ternary Alloys. J Phys Chem Sol (2019) 132:213–21. doi:10.1016/j.jpcs.2019.04.031

CrossRef Full Text | Google Scholar

130 Mahmood W, Bukhtiar A, Haroon M, Dong B. First Principles Calculations on Theoretical Band gap Improvement of IIIA-VA Zinc-Blende Semiconductor InAs. Int J Mod Phys C (2020) 31:2050178. doi:10.1142/s0129183120501788

CrossRef Full Text | Google Scholar

131 Hedin L. New Method for Calculating the One-Particle Green's Function with Application to the Electron-Gas Problem. Phys Rev (1965) 139:A796–A823. doi:10.1103/physrev.139.a796

CrossRef Full Text | Google Scholar

132 Heyd J, Scuseria GE. Assessment and Validation of a Screened Coulomb Hybrid Density Functional. J Chem Phys (2004) 120:7274–80. doi:10.1063/1.1668634

PubMed Abstract | CrossRef Full Text | Google Scholar

133 Kotani T. Quasiparticle Self-Consistent GW Method Based on the Augmented Plane-Wave and Muffin-Tin Orbital Method. J Phys Soc Jpn (2014) 83:094711. doi:10.7566/jpsj.83.094711

CrossRef Full Text | Google Scholar

134 Kotani T, Kino H, Akai H. Formulation of the Augmented Plane-Wave and Muffin-Tin Orbital Method. J Phys Soc Jpn (2015) 84:034702. doi:10.7566/jpsj.84.034702

CrossRef Full Text | Google Scholar

135 Malone BD, Cohen ML. Quasiparticle Semiconductor Band Structures Including Spin-Orbit Interactions. J Phys Condens Matter (2013) 25:105503. doi:10.1088/0953-8984/25/10/105503

PubMed Abstract | CrossRef Full Text | Google Scholar

136 Deguchi D, Sato K, Kino H, Kotani T. Accurate Energy Bands Calculated by the Hybrid Quasiparticle Self-consistentGWmethod Implemented in the Ecalj Package. Jpn J Appl Phys (2016) 55:051201. doi:10.7567/jjap.55.051201

CrossRef Full Text | Google Scholar

137 Otsuka J, Kato T, Sakakibara H, Kotani T. Band Structures for Short-Period (InAs)n(GaSb)nsuperlattices Calculated by the Quasiparticle Self-Consistent GW Method. Jpn J Appl Phys (2017) 56:021201. doi:10.7567/jjap.56.021201

CrossRef Full Text | Google Scholar

138 Taghipour Z, Shojaee E, Krishna S. Many-body Perturbation Theory Study of Type-II InAs/GaSb Superlattices within the GW Approximation. J Phys Condens Matter (2018) 30:325701. doi:10.1088/1361-648x/aacdce

PubMed Abstract | CrossRef Full Text | Google Scholar

139 Huang J, Ma W, Wei Y, Zhang Y, Cui K, Shao J. Interface Effect on Structural and Optical Properties of Type II InAs/GaSb Superlattices. J Cryst Growth (2014) 407:37–41. doi:10.1016/j.jcrysgro.2014.08.020

CrossRef Full Text | Google Scholar

Keywords: InAs/GaSb, type II superlattice, band structure, bandgap, theoretical calculation and simulation

Citation: Fang S, Hao R, Zhang L, Guo J and Liu W (2022) Simulation of the Band Structure of InAs/GaSb Type II Superlattices Utilizing Multiple Energy Band Theories. Front. Phys. 10:822800. doi: 10.3389/fphy.2022.822800

Received: 26 November 2021; Accepted: 10 March 2022;
Published: 14 April 2022.

Edited by:

Zhaoxin Liang, Zhejiang Normal University, China

Reviewed by:

Souraya Goumri-Said, Alfaisal University, Saudi Arabia
Teng Yang, Chinese Academy of Sciences (CAS), China
Bhaskaran Muralidharan, Indian Institute of Technology Bombay, India

Copyright © 2022 Fang, Hao, Zhang, Guo and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ruiting Hao, cnVpdGluZ2hhb0BzZW1pLmFjLmNu; Wuming Liu, d2xpdUBpcGh5LmFjLmNu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.