The Formation of an Interface and Its Energy Levels Inside a Band Gap in InAs/GaSb/AlSb/GaSb M-Structures

TL;DR

The paper explores how interface energy levels in InAs/GaSb/AlSb/GaSb M-structures affect optical absorption in materials for infrared applications.

Contribution

The study identifies a novel energy level within the band gap caused by InSb-like interfaces in M-structures.

Findings

An additional energy level within the band gap is found to originate from InSb-like interfaces.

Optical absorption increases significantly under optimized interface thickness and temperature.

HRXRD confirms the high crystal quality of the fabricated M-structures.

Abstract

We studied specially designed InAs/GaSb/AlSb/GaSb M-structures, a type-II superlattice (T2SL), that can serve as active materials for short-wavelength infrared (SWIR) applications. To obtain the dispersion relation of the investigated M-structures, k·p perturbation theory based on the eight-band model implemented in the nextnano++ v1.18.1 (nextnano GmbH, Munich, Germany) software was used. Numerical band-gap engineering and dispersion calculations for the investigated M-structures (composed of 6/1/5/1 monolayers, with InSb interfaces included) revealed the presence of an additional energy level within the energy gap. This energy level originates from the InSb-like interfaces and does not appear in structures with different layer or interface thicknesses. Its properties strongly depend on interface thickness, temperature, and strain. Numerical calculations of the probability density function |Ψ|2, absorption coefficients, and optical absorption spectra at varying temperatures demonstrate that, under specific conditions, such as an optimised interface thickness and temperature, optical absorption increases significantly. These theoretical results are based on structures fabricated using molecular-beam epitaxy (MBE) technology. High-resolution X-ray diffraction (HRXRD) measurements confirm the high crystallographic quality of these M-structures.

Figures

Full text

1. Introduction

Experimentally proposed by B-M. Nguyen et al. in 2007 [1], the superlattice InAs/GaSb/AlSb/GaSb has been deeply investigated because its M-structure has a wide spectrum of valence band-tuning and engineering possibilities [2,3,4], allowing for an effective mass to be developed and ensuring its continued growth [1]. This experimental concept aligns with theoretical descriptions [5]. The M-structure can be used in two different ways: firstly, as a material for absorbers in the short-wavelength infrared (SWIR) range and, secondly, as a barrier material operating in mid-, long-, and very-long wavelength infrared (IR) radiation ranges. As the absorber, the M-structure can easily replace the binary type-II superlattices (T2SL) InAs/GaSb. It is essential to emphasize that, in theory, the InAs/GaSb superlattice is limited, and it is difficult to produce good-quality structure during the growth process. The proper binary InAs/GaSb superlattice should be about 2/8 ML [6] or 4/8 ML [7] (both cases pertain to the SWIR region). It has been proven that in situations when there is contact between the InAs and GaSb layers, for example, InSb or GaAs, the choice of a thin layer of InAs or GaSb does not compensate for the strain (compressive or tensile type). The M-structures that are usable in the SWIR range were described by, e.g., [4,8,9,10]. In the second example, the M-structure is usable as a barrier material for the holes in the T2SL InAs/GaSb superlattice photodetectors in the mid-wavelength infrared (MWIR) range [10] and long-wavelength infrared (LWIR) range [11].

In addition to their the application potential, the structures with a broken gap, such as the M-structures, have very interesting physical properties. Theoretical and experimental investigations show that in III-V quantum well (QW) [12,13,14], as well as in T2SL, e.g., InSb/InAsSb Gal-free T2SL [15,16], Dirac-like carrier dispersion can be observed. In the case of InAs/InAsSb T2SL structures, Dirac-like dispersion can be observed for structures with ultra-thin layers within the superlattice period from 3 to 6.2 nm, providing an effective band gap from 70 meV to 0 meV [15]. In Ref. [13], the authors define novel three-layer InAs/GaSb/AlSb QWs as having an “InAs geometry” or a “GaSb geometry” if they contain two symmetrical layers of InAs or GaSb, respectively. The “GaSb geometry” of such three layer InAs/GaSb structures is due to the presence of one InAs QW with two GaSb quantum barriers (QBs) and two symmetrical AlSb layers on both sides of the structures. Such architectures are similar to those of HgTe/HgCdTe or CdTe/HgCdTe QW, in which a Dirac-like dispersion can also be observed [17].

However, the results presented above, especially when we are talking about the three QW layers presented in Ref. [13], are theoretical predictions. When using real structures obtained via molecular-beam epitaxy (MBE) technology, it is necessary to take into account the interface layers that appear during the growth process. Such interfaces can be the result of non-intentional and intentional work. When we are talking about the use of M-structures for IR purposes (as a barrier material in absorber layers), the interfaces play a crucial role in obtaining the proper lattice constant that should match the GaSb buffer layers (or GaSb substrates). In such cases, InSb-like interfaces and/or GaAs-like interfaces should be created to ensure that lattice constant mismatches are minimal [18]. Such intentional interfaces influence the energy spectrum and effective energy gap, which is important for their application. Our first approach clearly shows that, besides the minimal lattice constant mismatch, which is a result of the use of proper interfaces, additional energy levels appear inside the energy gap. These additional energy levels are proved through calculating $[eqn]$ , which derives from the interface (IF) layers [18]. Many authors have presented numerical calculations for InAs/GaSb T2SL, taking into account the parameters of each compound, which are a part of such structures, e.g., [19]. In this work, we focused on M-structures, for which no such analyses have been conducted to date.

In this work, we focused on the properties of M-structures with additional energy levels. The numerical calculations were based on the real structures obtained by MBE technology. As we can see, this additional energy level has an energy dispersion relation that is characteristic of the Dirac-like dispersion and appears in special architectures with strictly defined InSb interface widths and well-defined temperature ranges. This case is strongly different from that presented for, e.g., T2SL InAs/InAsSb [15,16], as well as the three QW layers [13], and clearly shows that the interface layers should be taken into account in such cases. In the investigated M-structures, a Dirac-like dispersion relation is observed inside the energy gap, which is about 600 meV, which increases the maxima of the optical absorption at different temperatures.

2. Materials and Methods

All presented diagrams were created from data simulated using nextnano++ software (commercial Schrödinger–Poisson current solver) [20]. M-barriers were simulated as a 6-layer structure with an infinite number of periods and also as 5-period M-structures with 0.1 nm spacing in a one-dimensional grid. The order of layers in one period of the modeled M-structure was as follows (from value 0 nm on a horizontal scale, denoting the end surface of the entire M-structure, to the right of the x-axis in the substrate): GaSb/AlSb/GaSb/InSb/InAs/GaAs. The InSb (and GaAs) layer serves as the interface (IF) between the InAs layer and adjacent GaSb layers. The band edges, probability density functions, and dispersion diagrams were calculated using an 8-band k·p model (evaluated from 6-band k·p parameters) [21,22]. The k·p model used to calculate the dispersion diagrams included 16 energy levels of electrons and holes along the [010] direction that were additionally perpendicular to the surface direction ( $[eqn]$ ) in k-space for an infinite number of periods. To calculate the optical absorption spectrum for interband and intraband transitions, the quantum approach was used with the 8-band k·p model (evaluated from 6-band k·p parameters), with 80 electron and hole energy levels and dispersion relation calculations. The size of the $[eqn]$ space was set to ≤0.5 ( $[eqn]$ ), and a 33 × 33 point grid was used for spatial integration in absorption spectrum calculations (for circularly polarised light $[eqn]$ , where x is the growth direction). The input files for absorption calculations were based on the “1D tutorial for interband transitions: Frankenberger” [23]. To reduce the computation time of the absorption spectra in nextnano++ for different parameter values, e.g., temperature, HTCondor was used as a functionality dedicated for distributed parallel computation in nextnanomat software [24], working on a cluster containing 15 compute nodes with a 4-core 3.2 GHz Intel i5-3470 processor and 16 GB RAM per node. The simulations included the evaluation of strain and lattice constants as functions of temperature. The material parameters used for simulations were taken from the material database included with the nextnano++ software.

The following parameters were used for the simulation:

• (a)M-structure with an infinite number of periods:
• -The material layers were designed in a zinc blende structure.
• -The GaSb substrate’s orientation was set to $[eqn]$ = [100] and $[eqn]$ = [010].
• -The simulation grid size for the x growth direction was set to 0.1 nm.
• -For dispersion, with respect to the grid size of the $[eqn]$ vector (for the direction perpendicular to the growth direction), it was set to about 0.05 nm^−1^, and for the $[eqn]$ direction (along the growth direction), it was set to about 0.07 nm^−1^.
• -The temperature dependence of the energy gap and the temperature dependence of the lattice constant were turned on.
• -Calculated bands: Gamma, HH, LH, and SO.
• -Strain was set as pseudomorphic.
• (b)Additional parameters were used for simulations for a 5-period M-structure (for dispersion relation and optical absorption):
• -20 nm thick AlSb buffer (with 0.1 nm grid size).
• -1 μm thick GaSb substrate (with 10 nm grid size).
• (c)Additional parameters were used for simulations for a 5-period M-structure (parameters appropriate for optical absorption):
• -0.5 size of the $[eqn]$ -space (relative to the Brillouin zone) for the computation of optical transitions.
• -16 $[eqn]$ points (along the y and z axes) where optical transitions intensities were evaluated.
• -8 points between $[eqn]$ points where optical transitions intensities were interpolated.
• -Photon energy steps were set to 0.0001 eV.
• -The Gaussian broadening of energy levels for optical transitions was set to 1.0 × $[eqn]$ eV.

3. Results and Discussion

3.1. HRXRD Measurements

The numerical calculations presented in this paper in Section 2 were carried out using the data (Table 1) obtained from the high-resolution X-ray diffraction (HRXRD) measurement of M-structures synthesized by MBE technology. Detailed information about the technological processes of investigated M-structures can be found in Ref. [18]. HRXRD measurements were performed using a Malvern Panalytical diffractometer with Cu K_α1_ radiation ( $[eqn]$ = 1.540598 Å). The 2 $[eqn]$ scan was recorded around the GaSb (004) reflection.

In the context of X-ray diffraction (XRD) analyses, the values of $[eqn]$ and $[eqn]$ , which correspond to the components of the scattering vector in inverse space, proved to be crucial: $[eqn]$ describes the component parallel to the sample surface, and $[eqn]$ describes the perpendicular component. Analyses of the signal intensity distribution along these axes provided valuable information about the stress, composition, and thickness of the epitaxial layers. The coordinates of the diffraction space are expressed in terms of angular positions as follows:

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ , and the $[eqn]$ —wave vector represents the incident beam, for which its length is $[eqn]$ . The $[eqn]$ —wave vector represents a reflected beam that has the same length $[eqn]$ .

Analyses of the changes along $[eqn]$ allows the assessment of stresses parallel to the surface and the identification of lattice mismatches in the epitaxial plane. However, the intensity distribution along $[eqn]$ provides information about perpendicular stresses, layer thicknesses, and the ordering of the crystal structure. When an epitaxial layer grows on a substrate with a different lattice constant, stresses occur due to lattice mismatch. The shift of the peaks towards $[eqn]$ indicates deformations in the plane of the layer (tension or compression). The shift of the peaks in the $[eqn]$ direction is related to the perpendicular stresses resulting from the compensation of the lattice mismatch in the growth direction. Inverse space intensity mapping can reveal crystal mosaicism, i.e., the degree of rotation and displacement of crystal domains. Blurring in the $[eqn]$ or $[eqn]$ direction indicates the presence of defects such as dislocations that are the result of lattice mismatch. The width of the peaks along $[eqn]$ and $[eqn]$ reflects the density of dislocations and other defects. A high-quality crystal is characterised by sharp and well-defined diffraction peaks. Disorder (e.g., surface roughness or unevenness in the layer thickness) causes the peaks to be blurred.

Each MBE process for each investigated sample was programmed to obtain different widths of InSb-IF, which were obtained using different Sb-soak times and different shuttering times of the In flux during IF growth. The results of the HRXRD measurements, together with the numerical simulations for all investigated samples in Figure 1, are presented. For all investigated samples, we were able to obtain a high and narrow peak, which were assigned to the GaSb substrate (see Figure 1), as well as the 0-th order peak, which originates from the T2SL M-structures. An additional peak from the 20 nm AlSb layers is also visible. Such thin films were grown between GaSb buffer layers and M-structures.

Figure 2 shows the reciprocal space maps (RSMs) for the two superlattices, which showed a smaller mismatch in the crystal lattice parameters. To make the results comparable, all structures were measured with the same parameters. The maps show that the reciprocal lattice points of the GaSb substrate and the ’0’ point are very close to each other, indicating minimal lattice mismatch. The lattice points of the zero-order satellite peaks are located in narrow ranges for the $[eqn]$ and $[eqn]$ axes, with values of −0.00040 Å^−1^ $[eqn]$ 0.00071 Å^−1^, and −0.00028 Å^−1^ ≤ $[eqn]$ ≤ 0.00059 Å^−1^. Both maps, (a) and (b), indicate an excellent fit of the structure to the substrate, with no evidence of relaxation in the overlying layer. The structures observed are of high crystal quality, as confirmed by the presence of a surface streak, which is known as a crystal truncation rod (CTR), on which Pendellösung peaks are visible along the $[eqn]$ axis. Also noticeable on the maps are the streaks associated with the analyser, which run parallel to the diffraction beam (A). In addition, subtle diffuse scattering (DS) can be seen, and it is visible as red ellipses on the maps.

3.2. Nextnano++ Computer Simulations

The dispersion relation obtained from nextnano++ simulations for all investigated samples (with an infinite number of periods) using the data from Table 1 is shown in Figure 3. For values of k greater than zero, this concerns the dispersion relation for the k-vector parallel to the layers in the [010] direction. For negative values of k, the dispersion relation is for the k-vector perpendicular to the layer plane, and it is denoted by $[eqn]$ (the 0 eV on the vertical energy scale is set to the energy value in the highest hole level at k = 0). The values shown for the energy gap ( $[eqn]$ ) were obtained at the $[eqn]$ point, calculated as the difference between the energy of the lowest electron level and the energy of the highest hole level (excluding the electron energy level coming from the InSb interface). The red/blue lines show hole/electron energy levels, and the blue dashed line shows the electron levels located inside the energy gap and those coming from the InSb interface (EGIF)—which occur at 300 K only for samples 35E183B1M and 35E183B1M, as shown in Figure 3b and Figure 3d, respectively.

For the calculation of the optical absorption spectrum, the samples were modelled as five-period M-structures with a 20 nm thick AlSb buffer and a 1 μm thick GaSb substrate (for quantum and optical calculations, the buffer and substrate are not included). For this reason, dispersion calculations were also performed for a five-period M-structures to compare with the results obtained for M-structures with an infinite number of periods, and the results are shown in Figure 4. On the left-hand side of each subfigure in Figure 4, there are results of calculations for the five-period model (with 80 electron and hole energy levels); on the right-hand side, there are results for the infinite number of period models (0 eV on the vertical energy scale is set to the energy value in the highest hole level at k = 0). It is worth noting that in the five-period model, the number of electron levels from the IF located inside the energy gap is equal to 10 compared to two levels for an infinite number of period models (this is visible as wider lines on the left side of the subfigures in Figure 4).

Figure 5 (corresponding to the results in Figure 3) shows the calculated band edge diagram and probability density functions ( $[eqn]$ ) for electrons and holes obtained for a temperature of 300 K. The dotted blue line represents the $[eqn]$ for EGIF and is only visible (at a temperature of 300 K) for the 35E183B1M and 35E185B1M samples. The $[eqn]$ simulation results for the investigated samples and for a wide temperature range (from 10 K to 400 K) are shown in the form of a colour map in Figure 6, where the blue component of the colour represents $[eqn]$ , corresponding with EGIF. In this figure, it can be seen that the maximum value of $[eqn]$ for the IF InSb region occurs within a limited specific temperature range, and it is visible in purple for lower temperatures (as a sum of red and blue colours) and in cyan (as a sum of green and blue colours) on the map.

Fermi velocity $[eqn]$ can be obtained from the $[eqn]$ relation, $[eqn]$ (in the linear region of the dispersion relation), as shown in Figure 7, and this type of dispersion is represented by one and the two dotted blue lines (the one with higher energy) for the $[eqn]$ direction and at the lowest temperature when EGIF occurs. Similarly to Figure 7, Figure 8 shows that the dispersion relation for the highest temperature of EGIF still exists inside the energy gap; however, its slope is smaller, and the linear dispersion is also not clearly visible. The values of the Fermi velocity $[eqn]$ for all structures were calculated for $[eqn]$ = 0.05 nm^−1^ and plotted in Figure 9 as a function of $[eqn]$ for three temperatures: lowest, middle, and highest temperatures at which EGIF occurs. For $[eqn]$ , the calculated value of $[eqn]$ is around 200 km/s (or reaches negative values slightly below 0 km/s). For k in the range of 0–0.2…0.5 nm^−1^, the value of $[eqn]$ increases almost linearly to the maximum value of around 450–800 km/s, and only for higher values of k (within the range of 0.5–1.0 nm^−1^ does it reach a fairly constant level that is appropriate for the linear dependence of $[eqn]$ —(with a slight decrease in value for $[eqn]$ with respect to increased k values). However, when the temperature is lower (the energy values of EGIF are closer to the hole levels), the Fermi velocity $[eqn]$ reaches a relatively constant level for lower values of $[eqn]$ and has a slightly higher value than for the higher temperatures (where the values of energy of EGIF are closer to the electron levels).

The computed energy spectrum at the $[eqn]$ point as a function of temperature for the same investigated samples is presented in Figure 10. The EGIF and their ‘formation’ are visible in this figure. In the case of increasing temperatures, at lower temperatures, some of the hole levels (the red dots) rapidly change their energies; if the temperature is high enough, they appear inside the energy gap as a blue square, and finally, at the highest temperatures, they continue to move outside the energy gap at a ’normal’ electron level (the blue dots). The temperature range for the levels occurring inside the energy gap is the same as for the temperature range for $[eqn]$ corresponding to EGIF in Figure 6. Figure 6 also shows the temperature dependence of EGIF in the IF InSb layer of each of six samples with different IF InSb widths—except samples 35E184B1M and 35E187B1M with similar IF InSb widths.

For this purpose, a series of nextnano++ simulations were performed for a sample with parameters similar to 35E180B1M, but only the IF InSb width was varied. Figure 11 shows the energy spectrum (similarly to Figure 10) versus the width of IF InSb in the range of 0.01–1.00 nm, with 0.01 nm intervals at three temperatures of (a) 200 K, (b) 300 K, and (c) 400 K for parameters that are similar to the 35E180B1M sample, modelled as an infinite number of period structures. For the temperatures of 200 K/300 K/400 K/in Figure 11, the energy level of the electrons in the energy gap occurs in 4/3/2 ‘groups’, respectively, and it is located at the larger InSb IF width at higher temperatures. For the temperature of 200 K in Figure 11a, the first group of these levels occurs at an IF InSb width within the range of 0.12–0.15 nm, the second group at 0.37–0.40 nm, and the third group at 0.66–0.67 nm (for temperatures of 300 K and 400 K, the first group is 0.16–0.18 nm and 0.21–0.27 nm, respectively)—which is consistent with the results for all investigated samples presented in Figure 6, Figure 7, and Figure 8. The groups of levels are placed (more or less) at odd multiples of the IF InSb’s width corresponding to the first group of levels (e.g., at a temperature if 200 K: 0.12 nm, 3 × 0.12 nm = 0.36 nm, and 5 × 0.12 nm = 0.60 nm).

Figure 12 shows the band edge’s diagram and the probability density function $[eqn]$ of electrons and holes obtained from simulations at a temperature of 200 K for a modelled sample similar to 35E180B1M with an infinite number of period structures with different $[eqn]$ widths: (a) 0.13 nm, (b) 0.39 nm, and (c) 0.66 nm. The red line represents the $[eqn]$ holes, and the blue line represents the $[eqn]$ electrons. The dotted blue line represents electron $[eqn]$ coming from EGIF, with one/two/three maxima with respect to $[eqn]$ for a width $[eqn]$ equal to 0.13/0.39/0.66 nm respectively.

It is worth noting that the simulation results indicate that the energy gap decreases with an increase in temperature (Figure 13) or IF InSb width (Figure 14). However, when the temperature increases, a decrease in the electron level contributes to gap reductions (Figure 10), and when the IF InSb’s width increases, it contributes most to the increase in the hole level (Figure 11).

In order to check the optical properties and the influence of EGIF levels on them, the optical absorption coefficient was simulated using the quantum approach for all investigated samples modelled as five-period M-structures with a 20 nm AlSb buffer and a 1 μm GaSb substrate (however, for the quantum and the optical calculation, the buffer and substrate are not included). The simulation grid was 0.1 nm for the periodic M-structures and buffer and 10 nm for the substrate. The results for the optical absorption coefficient as a function of temperature within the range of 10–400 K, with intervals of 10 K in photon energy within the range 0–2 eV for all investigated samples, are shown as a colour maps in Figure 15. Two important phenomena are visible. Firstly, the optical absorption edges are around 0.6–0.7 eV (as a vertical, slightly diagonal line separating the black and blue areas), and they are comparable to the $[eqn]$ previously calculated with dispersion relation data. Secondly, a significant increase in the absorption coefficient is visible for temperature ranges close to the temperature range where EGIF levels occur for each sample with respect to the 0–2 eV photon energy range (for sample 35E186B1M, the temperature at 400 K is too low for exhibiting a decrease in absorption at high temperatures). The increase in the absorption coefficient (especially for photon energies higher than the absorption edge) occurs at a slightly lower temperature than the temperatures assumed to be the lowest for IF, as shown in Figure 7. On the other hand, the maximum absorption for the lowest photon energy occurs at a temperature close to the middle between the lowest and highest temperatures for the EGIF shown in Figure 7 and Figure 8, respectively. Above this temperature, there is a sharp decrease in absorption, especially for photon energies that are lower than the absorption edge (and with absorption decaying relatively slowly for photon energies of >1 eV). To better visualise the absorption spectrum near the absorption edge as a function of temperature, the same data are plotted on a logarithmic scale within the photon energy range of 0–1 eV in Figure 16.

Figure 17 shows, on a logarithmic scale, the optical absorption spectrum obtained for several temperatures within the range of 10–450 K for all samples investigated in order to better visualise the shape of the absorption spectra, the absorption edge for the temperatures corresponding to EGIF (the three solid lines), and the slightly lower and higher (the dashed lines) temperatures. The same data as shown in Figure 17 are presented in Figure 18 on a linear scale where a significant increase in the optical absorption coefficient is shown for the lowest temperatures when EGIF exists (presented in Figure 7) and when temperatures are close to the middle values, between the lowest and highest temperatures (presented in Figure 8). The absorption spectra for the highest temperatures when EGIF still exists are comparable to the spectra obtained at higher temperatures.

The contribution of optical transitions involving EGIF levels in the absorption spectrum was analysed. Figure 19 shows, relative to a logarithmic scale (and Figure 20 relative to linear scales), the optical absorption spectrum obtained from nextnano++ simulations for temperatures of (a) 200 K, (b) 220 K, and (c) 250 K for the investigated sample 35E180B1M. The solid black line (denoted as ‘All’) shows the total optical absorption coefficient as a sum of the individual contributions shown in dashed lines. For the lower and middle temperatures of (a) 200 K and (b) 220 K, the main contributions come from transitions from EGIF to other levels (the red dashed line denoted as ‘IF-other’). However, for the highest temperature at which EGIF is still present ((c) 250 K), the main contribution comes from transitions from other levels to EGIF (the green dashed line denoted as ‘other-IF’), and contributions from EGIF to other levels (the red dashed line) are completely damped for photon energies greater than about 0.15 eV (for photon energies < 0.15 eV, relatively low absorption occurs; see Figure 19). The contributions of intraband transitions (the blue dashed line denoted as ‘IF-IF’) for photon energy within the range of 0.1–2.0 eV are absent at these three temperatures (only visible as sharp absorption peaks for photon energy < 0.1 eV, with the highest level at a temperature of 200 K and the lowest level at a temperature of 250 K; see Figure 19). The contribution of transitions between non-EGIF levels (the grey dashed line denoted as ‘other-other’) is approximately at the same level for all these three temperatures.

In order to check the influence of the IF InSb’s width on the absorption spectrum, a series of simulations for five-period M-structures with 80 electron and hole levels (and a 33 × 33 point grid in k-space) was performed for a sample with parameters similar to 35E180B1M, but only the IF InSb’s width was varied within the range 0.01–0.50 nm, with 0.01 nm intervals at three temperatures. The selected range of the InSb interface’s thickness is based on actual measurements of this thickness using the HRXRD method. In reality, 0.5 nm corresponds to an interfacial thickness of approximately 1.5 monolayer (ML). Such thick interfaces are not actually used, which is why this is the upper limit considered in the calculations presented. As can be observed in Table 1, the actual interface’s thicknesses are at a level of less than one atomic monolayer, i.e., a thickness of about 0.3 nm. The determined thicknesses are accurate relative to 0.01 nm; hence, this value was also adopted in calculations as the smallest step. Figure 21 shows the results obtained for the following temperatures: (a) 200 K, (b) 300 K, and (c) 400 K. Similarly to Figure 15 and Figure 16, the optical absorption edges are clearly visible. The first maximum of the optical absorption at temperatures of 200/300/400 K exists at 0.14/0.16/0.24 nm of $[eqn]$ —which is consistent with the data presented in Figure 11. The second maximum of the optical absorption at a temperature of 200 K is also visible (but weaker) for IF InSb widths of approximately 0.38 nm, and this value is also comparable with the data presented in Figure 11 (at a temperature of 300 K, higher absorption is visible for IF InSb widths within the range of 0.45–0.48 nm).

4. Conclusions

In this paper, we investigated specially designed InAs/GaSb/AlSb/GaSb M-structures modelled as T2SL structures, and they were proposed as active materials in IR detectors for the SWIR region. We calculated dispersion relations, energy gaps, Fermi velocities, probability density functions $[eqn]$ , and absorption coefficients for all six investigated M-structures using nextnano++ software over a wide temperature range from 10 K to 400 K with 10 K steps (or for several temperatures). In addition, we performed calculations of the absorption coefficient and energy spectrum for one M-structure (35E180B1M) with varying InSb-like interface widths over a wide range (up to 0.5 nm with 0.01 nm steps) for the three temperatures. Real material samples of the investigated M-structures were grown using MBE technology, and HRXRD measurements were performed: 004 reflection $[eqn]$ scans for all samples and RSM maps for the two investigated samples (35E183B1M and 35E185B1M). The results of the HRXRD measurements indicated the thickness values of each layer and confirmed the very good quality of the crystallographic properties of the grown M-structure material samples.

Detailed conclusions from the results obtained from the nextnano++ simulations of the investigated samples indicate the following:

1. EGIF levels form within a specific temperature range and a specific range of IF InSb widths. The specific width of IF InSb for EGIF is higher at higher temperatures. The specific temperature for EGIF is higher for higher IF InSb widths.
2. In simulations with an infinite number of M-structure periods, the number of EGIF levels is equal to two. In simulations with a five-period M-structure, the number of EGIF levels is equal to 10. The results of the dispersion relation and the EGIF are comparable between these two types of simulations.
3. The values of the calculated Fermi velocity $[eqn]$ for the EGIF levels for the $[eqn]$ direction are as follows: for $[eqn]$ , they range from slightly below 0 to about 500 km/s (average about 200 km/s) depending on the temperature and M-structure parameters. For $[eqn]$ in the range of 0–0.2…0.5 nm^−1^, it increases almost linearly for a maximum value of approximately 450…800 km/s (for lower temperatures, the maximum value of $[eqn]$ is higher and is achieved for smaller values of $[eqn]$ ). For the value of $[eqn]$ in the range of 0.2–1.0 nm^−1^, velocity $[eqn]$ reaches a relatively constant level corresponding to linear dependence $[eqn]$ —with a slight decrease in the value of $[eqn]$ with an increase in the value of $[eqn]$ .
4. The EGIF manifests as the maxima of $[eqn]$ for electrons located in the IF InSb layer. At specific temperatures, EGIF exists for an odd multiplication of the smallest IF InSb base’s width when the first EGIF occurs at this temperature. For the ‘first’ EGIF, there is only one maximum for the EGIF $[eqn]$ function. For the ‘second’ EGIF, there are two maxima, and for the ‘third’ EGIF, there are three maxima. The amplitudes of the ‘second’/’third’ maxima are approximately two/three times lower than the ‘first’ maxima.
5. The computed optical absorption coefficient indicates the optical absorption edge corresponding to the computed $[eqn]$ . The maxima of absorption versus temperature at specific widths of the IF InSb correspond to the computed EGIF—a high increase in optical absorption at the lowest temperature (and also at lower temperatures) and the middle temperature where EGIF is present. The increase in optical absorption within this temperature range is mainly due to the transitions from EGIF to other electron levels (at the highest temperature where EGIF still exists, the main contribution to optical absorption is from transitions from hole levels to the EGIF).
6. The absorption maxima as a function of the IF InSb’s width at specific temperatures correspond to the calculated EGIF. At temperatures of 200 K and 300 K, two optical absorption maxima occur for the IF InSb width within the range of 0.01–0.50 nm.