## 1. Introduction

Low‐dimensional semiconductor heterostructures, otherwise known as semiconductor nanostructures, have tremendously revolutionized both the technical and the fundamental aspects of semiconductor industry in terms of device applications. With the ability to grow clean and high‐quality samples, device implementations have become a huge success [1–3]. When the dimensions of a region with free carriers (electrons) are reduced as compared to the bulk and approach the deBroglie wavelength, the electronic motion is quantized, thus resulting in carrier confinement that is quantum mechanical in origin. The phenomenon has been widely used for carrier confinement in one, two and three dimensions that consequently gives rise to nanostructures such as quantum wells, quantum wires and quantum dots, respectively. Due to the quantum confinement, the energy bands (i.e. the conduction and valence bands) are quantized into discrete energy levels/bands and are no longer continuous as in the bulk systems. These quantized energy states are known as subbands for 2D or 1D systems and sublevels for 0D systems. The energetic spacings between these quantized subbands and the sublevels are very important parameters that define the device applications both from an optical and from an electrical point of view.

The intersubband spacings in GaAs‐based 2D systems are typically in the order of 10–30 meV [4, 5], as seen in the case of two‐dimensional electron gas (2DEGs) with a triangular confinement potential formed across a GaAs/Al_{x}Ga_{1−x}As heterojunction (*x* being typically 0.3) of a high electron mobility transistor (HEMT) structure. The intersubband transitions (ISTs) typically cover the terahertz (THz) or far‐infrared region of the electromagnetic spectrum. However, in the case of a square potential well or in a different material system such as GaN/AlGaN heterojunction, these spacings can be designed to be even in the mid‐infrared or near‐infrared region. Stacking of quantum wells can further enhance the response of intersubband resonance (ISR), and such designs are the key for various applications like photodetectors or intersubband lasers [6]. One of the very common and sophisticated examples in this regard is the quantum cascade laser [7–9], which is based on the cascade phenomena and intersubband transitions across many layers of quantum wells. Such compact and powerful lasers are used for practical applications in THz spectroscopy [10–13], sensing technology [14, 15], biomedical applications [16, 17] and also in security applications [11, 18]. Structures based on quantum wells have also made significant advancement in the detector technology, for example, quantum well infrared photodetectors [19, 20]. In this chapter, we present a broad overview of the ISTs in a 2DEG formed at the GaAs‐Al_{x}Ga_{1−x}As interface of a HEMT structure. We also discuss possible methods to probe the spacing between the subbands and also to tune them significantly by applying an external bias across the sample. Furthermore, we present a fundamental study on the coupling of the ISRs with the 2DEG cyclotron resonance in the presence of tilted magnetic fields. The knowledge of ISTs and the ability of wide electrical tuning of these resonances are then exploited to study the light‐matter interaction at THz frequencies in these HEMT structures. The integrated device with 2DEG in a HEMT structure and metamaterials (frequency‐selective artificially designed structures) is electrically driven from an uncoupled to a coupled regime of light‐matter interaction and then again back to the uncoupled regime. A strong coupling is thus observed when the frequencies of both systems are brought in resonance with each other, manifested as an avoided crossing at that point.

## 2. High electron mobility transistor design

The low‐dimensional charge carriers, trapped in the heterojunction of the HEMT design, form the core of such field‐effect transistors. This transistor design also goes by the name of modulation‐doped field‐effect transistors (MODFET). These designs are used in various high‐power [21] and high‐speed [22] electronics, high‐resolution imaging [23] and various gas, chemical and biomedical applications [24]. We begin with the design concept of this semiconductor heterostructure along with an overview of its band structure (see **Figure 1(a)**) that is obtained by solving the *Schrödinger‐Poisson’s* equations self‐consistently [25, 26] and adding the band discontinuity at the heterojunctions. A schematic of the layer sequence of the transistor structure is shown in **Figure 1(b)**. On a semi‐insulating GaAs substrate/wafer, we start the molecular beam epitaxy (MBE) growth by typically a 50‐nm‐thin GaAs layer. Then, approximately 10 periods of a GaAs/AlAs short period superlattice (SPS) are grown (not shown in the band diagram). The SPS layers help to smoothen the surface of the bare substrate for the later epitaxial growth and trap eventually surface‐segregating unintentional impurities, which have always a tendency to stick at the stoichiometric interfaces of GaAs/AlAs. Moreover, this SPS keeps unwanted charge carriers away, forbidding them to tunnel into the 2DEG layer grown on top. Since the substrate is typically undoped (or semi‐insulating), the conduction (or valence) band has no curvature at this point corresponding with Poisson’s equation, which states that the charge density is proportional to the second derivative of the potential with respect to the space coordinate. After the growth of the SPS layer, the first charged layer is the 2DEG that is formed at the heterojunction of the undoped GaAs and an undoped Al_{0.33}Ga_{0.67}As spacer layer. Since the 2DEG is essentially electrons and negatively charged, the conduction band curves downwards and reaches the maximum slope at the heterojunction between the GaAs and the Al_{0.33}Ga_{0.67}As layer, at which point the conduction band (*E _{c}*(

*z*)) jumps by Δ

*E*due to the band discontinuity. This is followed by the Al

_{c}_{0.33}Ga

_{0.67}As spacer layer where charge carriers are absent and the slope of the conduction band almost remains constant. In the doped Al

_{0.33}Ga

_{0.67}As layer, the positive charges of the donor ions cause the band to bend upwards, thus reversing the slope. Further moving to the GaAs layer,

*E*(

_{c}*z*) jumps downwards due to the band discontinuity and continues with a constant slope. On top of the GaAs layer, AlAs/GaAs SPS (also known as blocking barrier) is grown to prevent leakage of charge carriers in and out of the 2DEG channel and also to prevent leakage of surface charges into the channel. Finally, the band hits the gate grown on top of the sample with a barrier height equivalent to the Schottky barrier height. Ideally, metals (e.g. Cr or Au) are evaporated on the sample to serve as gates after the completion of the growth. The samples are typically grown by MBE. While a lot of work has been done previously using metallic Schottky gates, nonetheless, these gates suffer from huge drawbacks. These gates limit the forward bias voltage to the turn‐on voltage of the Schottky diode. Furthermore, they fail to grow lattice matched on the semiconductor, are poly‐crystalline and thus induce potentially a lot of strain on the semiconductor layer below. Moreover, they oxidize over time and thus may become highly ohmic. Due to high reflectivity and certain Drude absorption of their free charge carriers, such gates are also opaque to the incident light, thus limiting their application in optoelectronic devices. Recently, we have introduced epitaxial, complementary‐doped, electrostatic and transparent gates that are grown on top of the sample [27–29]. These gates are grown within the UHV conditions of the MBE and thus incorporate a minimum of the unwanted impurities, leading to unpreceded gate perfection, reliability and reproducibility.

These gates circumvent all the abovementioned disadvantages of Schottky gates and are typically composed of a 25‐nm‐thick bulk carbon‐doped GaAs layer (with an acceptor density of *N*_{A} = 3 × 10^{18} cm^{−3}) followed by approximately 40 periods of carbon‐delta‐doped and 0.5‐nm carbon‐doped GaAs layers with an average acceptor density *N*_{A} = 1 × 10^{19} cm^{−3}. In order to solve Poisson’s equation for the evaluation of the band structure, the knowledge of the charge density is necessary. However, it is not possible to calculate the density of charge carriers until the energy bands are known, thus requiring a self‐consistent mechanism that is otherwise adopted in the 1D Poisson solver [26].

## 3. Characteristics of HEMTs

After being introduced in 1980s, these transistors based on high‐mobility modulation‐doped heterostructures have revolutionized the semiconductor industry in terms of being the most high‐performance compound semiconductor FETs.

Due to spatial separation of the electrons from the ionized impurities, the scattering between them is highly reduced as compared to the bulk semiconductors, enhancing the electron mobility especially at low temperatures where the abovementioned scattering mechanism is dominant. The spacer layer further increases the electron‐to‐donor separation. While the larger separation reduces the scattering mechanism, as a negative contribution, the carrier concentration is also reduced which reduces the performance. Hence, the spacer thickness should be optimized. Typical values range from 1 to 30 nm. In order to explain how the high mobility of the electrons in HEMT makes them fast transistors, we use the *Shockley’s gradual channel approximation model* for an FET operation, which states that the rate of change of saturation drain current,

where

High carrier concentration of 2DEG in the channel

High mobility by optimization of the spacer‐layer thickness

Low access resistance by using buried/recessed gates

Better confinement of carriers in the channel due to high barriers

Reduced interface and alloy scattering mechanisms, thus enhancing mobility

Typical transistor characteristic curves of GaAs/AlGaAs HEMT structures under dark and after 1 s of illumination with a near‐infrared light emitting diode (NIR LED) are shown in **Figure 2(a)** and **(b)**, respectively. Let us now briefly discuss the transistor operation. Even under zero bias or for a small positive voltage applied to the gate, an inversion layer is formed at the semiconductor surface, the *two‐dimensional electron gas*. Now, if a small source‐drain voltage is applied, a current will flow from the source to the drain through the conducting 2DEG channel. The channel here is highly conducting, that is, offers very small resistance, and the source‐drain current (*quasi‐Fermi level*), thus depleting the channel completely. This results in a zero source‐drain current even when the source‐drain voltage is increased. As a positive bias is applied, the channel is filled with mobile electrons, and with the increase of the source‐drain voltage, the source‐drain current first increases linearly, then non‐linearly and finally reaches a saturation value, as described earlier. A schematic of these situations is depicted in **Figure 2(c)** and **(e)**, while the simplified schematic of the conduction band diagram is shown in **Figure 2(d)**. When the structure is illuminated, it becomes rather difficult to deplete the channel completely out of electrons with the previously applied negative bias, and, moreover, the saturation current also increases with the same gate bias applied before.

### 3.1. Quantum confinement and intersubband transitions in HEMTs

One of the most popular terms in nanoscience is the *quantum confinement* that results from changes in the atomic structure as a consequence of direct influence of ultra‐small length scale on the energy band structure [30]. The length scale corresponding to the regime of quantum confinement ranges from 1 to 25 nm for typical IV, III–V or II–VI semiconductors [31]. This leads to the fact that the spatial extent of the electronic wavefunction is comparable to the particle size, making the electrons *feel* the presence of the particle boundaries and respond to changes in particle size by adjusting their energy. This phenomenon is known as the quantum‐size effect. Quantization effects become most important when the particle dimension of a semiconductor is near to and below the bulk semiconductor Bohr exciton radius (in bipolar systems) or the deBroglie wavelength (in unipolar systems), making the properties of the material size-dependent.

In low‐dimensional semiconductor nanostructures, the restriction of the electronic motion in one, two and three dimensions leads to the modification of the density of states (DOS) as compared to the bulk states. The electronic DOS is defined as the number of electronic states per unit volume per unit energy, the finiteness of which is a result of the *Pauli’s exclusion principle*, which states that only two electrons with opposite spins can occupy one volume element in the phase space [32]. The confinement of electronic motion results in the quantization of the conduction and the valence band. With the knowledge of these quantized states, their filling can be explained. The number of occupied subbands depends on the electron density and also on the temperature [33]. In a 2D system, the density of electrons per unit area,

where

where

where

In the limit of low temperature, where electrons are degenerate, the 2D electron density is given by [33]:

**Figure 3(a)** and **(b)** show a triangular potential well and a schematic of the filling of the subbands, respectively. Based on the position of the Fermi level, the corresponding subbands are occupied. Under triangular confinement potential (as in the HEMT design), the energy spacing decreases for higher subband energies and finally forms the continuum. When the Fermi energy **Figure 3(c)** shows a schematic of the intersubband transition from the filled ground subband to an empty excited subband. In quantum mechanics, Fermi’s golden rule is used to calculate the transition rate (i.e. the probability of a transition to occur per unit time), from one state with a given eigenenergy to another state of higher eigenenergy or to the continuum of energy eigenstates, subjected to some kind of perturbation. According to Fermi’s golden rule, this rate of transition, *electric‐dipole approximation*[1] -
(

where

where **Figure 3(c)**. The total wavefunction can be written as the product of the lattice‐periodic Bloch wave (for electrons in a crystal), *Bloch ansatz*, the envelope function reduces to the plane wavefunction; hence, the total wavefunction is given by:

where

where the first term indicates the interband transition and the second term is the intersubband transition. The first term has the dipole matrix element of the Bloch functions that explains the selection rule for the interband transition and an overlap integral of the envelope functions. In case of transitions within the subbands of the conduction or the valence band, the first term vanishes and the second term becomes more relevant in the *one‐band model* that consists of an overlap integral of the Bloch function and a dipole matrix element of the envelope function. Further simplification of the matrix elements of the envelope function gives:

It can be observed that only the third term,

The above equation states that the electric field of the incident radiation must have a component perpendicular to the semiconductor layers or parallel to the growth direction (which is a necessary condition) in order to couple to the ISTs. This is known as the *polarization selection rule* for the ISTs. In simple words, it states that the electric field vector of the exciting electromagnetic wave or at least a finite component of it must be perpendicular to the 2DEG. Another important quantity in this regard is the oscillator strength [34] defined as:

The above quantity is used to understand and compare the strength of the transitions between initial and final states in different physical systems and obeys the *Thomas‐Reiche‐Kuhn sum rule* [37]. It is important to note that for a symmetric quantum well, only parity changing transitions (odd‐even or even‐odd) are allowed due to the inversion symmetry of the potential well. However, for asymmetric quantum wells, like that of the triangular potential well, the inversion symmetry with respect to the quantum well centre is broken by some means (i.e. internal electric fields or band structure engineering, etc.). This leads to the relaxation of the selection rule, thus allowing transitions between all the subbands.

## 4. Intersubband‐Landau coupling under tilted magnetic fields

When a magnetic field is applied in a plane perpendicular to the semiconductor surface, the free electrons that carry the electric charge perform an orbital motion in the plane perpendicular to the magnetic field direction. This motion is quantized, and equally spaced levels (called the Landau levels) separated in energy are formed. The Hamiltonian of the quantum mechanical system thus gets decoupled into a magnetic and an electric component, and the energy spectrum consists of a series of Landau ladders for each subband. In the presence of a magnetic field, the Drude conductivity, normalized to

where

The transmission,

where

where

where

Using the perturbation theory, one can solve the above Hamiltonian treating

(19) |

where *i*th subband to *i*′th subband. The first two terms represent the linear zeroth‐order terms, corresponding to the subband energies and the Landau energies respectively. The third term represents the diamagnetic shift, and the fourth term results from the coupling Hamiltonian. This non‐resonant regime is known as the half‐field coupling, where the splitting is proportional to *T*(*B*)), and then the magnetic field is turned off, during which the same number of transmission scans (*T*(0)) are taken. At first, the transmission experiments are performed under no tilt of the magnetic field. A contour plot of the normalized transmission spectra for different fields in the range 3.6–4.2 T is plotted in **Figure 4(a)**. Clearly, the only visible resonance observed is the cyclotron resonance under perpendicular magnetic fields that scales linearly with the field. On introduction of the tilt (approximately 30°), a clear anti‐crossing is observed at around 3.9 T (see **Figure 4(b)**). The apparent observation of the satellite peaks (shown by black arrows in **Figure 4(c)**) in the presence of the magnetic field at the anti‐crossing point is a manifestation of the subband‐Landau coupling and hence the resonance splitting (the splitting between the two satellite peaks across the dominant cyclotron resonance). The spacings of the subbands are twice the value of the splitting at the anti‐crossing frequency (according to the half‐field coupling regime). This corresponds to

## 5. Tuning and probing of intersubband transitions electrically

Intersubband transitions are the most fundamental optical transitions that can be excited in low‐dimensional semiconductor nanostructures. The observation of ISRs is a result of the fact that the component of the incident infrared electric field perpendicular to the semiconductor layers or parallel to the growth direction selectively couples, thus exciting the electrons from the lower occupied subband to the higher empty subband. By applying a voltage across the structure, it is possible to deplete and selectively populate the subbands. Thus, a more direct scheme of transmission measurement is proposed to study the intersubband spacing in such semiconductor nanostructures (viz. HEMT) even in the absence of an external magnetic field.

The change in the transmission,

where

where

where *i*th state to the *i*′th state. The ISR frequency can however differ from the observed one due to resonance screening (depolarization shift) arising from the many body effects [46]. In the density‐chopping scheme, a certain number of scans are taken at the reference voltage (a voltage much below the threshold voltage), when the 2DEG is completely depleted of charge carriers. Then, the gate voltage is slowly increased to a value when the subbands start populating (this can be well seen from the capacitance‐voltage spectrum in **Figure 5(a)** where the change in capacitance is measured during the broadband absorption onset upon changing the gate voltage and modulating it with the LockIn technique). The same number of scans is taken at this higher gate voltage. The voltages are then changed alternatively, and the respective scans are co‐added and averaged over long measurement times. The long measurement time ensures that any drift arising from the complicated experimental setup can be averaged out to zero and a good signal‐to‐noise ratio is obtained. **Figure 5(b)** shows the density‐chopped transmission spectra of a HEMT sample. With the increase of the gate voltage, the conduction band is pulled below the quasi‐Fermi level. This results in the steepening of the triangular potential well, resulting in an increase of the intersubband spacing. Thus, by applying a bias on the gate, the intersubband resonance can be significantly tuned over a wide frequency (1–3 THz) or energy (4–12 meV) range, as can be seen in **Figure 5(c)**.

## 6. An access to the interior of HEMT via artificial structures

Artificial structures such as metamaterials are engineered in the sub‐wavelength sizes for certain desired properties. They are designed in assemblies of multiple individual elements called unit cells. These structures possess unique properties such as negative [47] or zero [48] refractive indices, magnetism at optical frequencies [49], etc. The special properties are not inherent to the materials but the design of the structures and the way electromagnetic field interacts with them. They can also be treated as planar cavities with certain resonance frequencies. When electromagnetic radiation with a certain polarization is incident on these structures, the electric or the magnetic field couples to the cavity and exhibits a resonance that is known as the cavity resonance or resonance frequency of the metamaterials. In transmission measurements, this appears as a dip at that particular resonance frequency. An array of interconnected double split‐ring resonators (see **Figure 6(a–d)**) is adopted for the metamaterial design, whose dimensions and the characteristic frequency response are first simulated by the standard finite difference time domain solver (like CST microwave studio). For simplicity and small computation time, only one unit cell, as shown in **Figure 6(e)**, is used for the simulation with a periodic boundary condition in the planar directions. Moreover, these meta‐atoms are placed far apart from each other to avoid any influence of inter‐meta‐atom interactions. For the right coupling of the electromagnetic radiation, the electric fields are confined in the two narrow capacitor arms of the double split‐ring resonator (see **Figure 6(e)**). Moreover, the fringing field effect ensures that there is a strong electric field component along the growth direction that extends over a few 100 nm [29]. This component of the electric field couples with the HEMT to excite the ISRs in accordance with the polarization selection rule as discussed before.

Two transmission minima (or dips) are observed—one at 1.2 THz and the other at 2.4 THz (see **Figure 6(f)**). The experimental characterization of the metamaterial array is performed by a standard THz time‐domain spectroscopy at room temperature, where the thermal energy, **Figure 6(f)**), the incident infrared radiation couples into the structures. When the metamaterials are oriented at 90° (black-dashed arrow in **Figure 6(f)**), the field does not couple. By normalizing the transmitted spectrum of the metamaterial at 0° with respect to the one at 90°, two transmission dips are obtained—one at 1.2 THz and the other at 2.4 THz as shown in **Figure 6(f)**.

Once the sample is cooled down to liquid helium temperatures, at first, the characterization of the voltage range is performed over which the density‐chopping measurements are to be taken. The change in capacitance with the gate voltage is measured by capacitance‐voltage spectroscopy, mentioned before, as shown in **Figure 7(a)**. A typical charging spectrum of 2DEG has a capacitance close to zero in the beginning and then as the gate voltage is increased, the conduction band is pulled below the quasi‐Fermi level and subsequently the 2DEG subbands are filled with electrons. The filling is observed as a steep increase in the capacitance. The region of interest is the steep slope, where increasing the gate voltage increases the 2DEG ISRs. This is due to the fact that with more positive gate voltage, the slope of the triangular potential confinement steepens and hence increases the subband spacings, thus shifting the ISRs to higher energies. This phenomenon is well known as the *quantum‐confined Stark effect*. It is also necessary to completely ionize the donor‐exchange centres (*DX centres*[2] -
). As more DX centres are ionized, less forward bias is required to charge the 2DEG subbands with electrons. With longer illumination, all the DX centres are successively ionized, leading to the shift of the charging slope in the capacitance‐voltage spectra towards more negative biases. The spectrum shown in **Figure 7(a)** is obtained after 3 h of continuous illumination. The shaded region in the charging spectra, shown in **Figure 7(a)**, indicates the region where the density‐chopped infrared transmission measurements are performed. The density‐chopping scheme is similar to that explained before, where the change in transmission, **Figure 7(b)**. It is observed that at low temperatures, the cavity resonance slightly shifted to a higher frequency as a result of the lower losses in the cavity in comparison to the room temperature measurements. At

where

where

where

## 7. Conclusion

In conclusion, we have reviewed the quantum mechanical phenomenon that governs various electrical and optical properties in the low‐dimensional semiconductor nanostructures such as a HEMT. We have demonstrated how one could electrically, or in combination with magnetic fields, probe and tune the intersubband transitions in the heterojunction of a HEMT structure. Such structures primarily have a triangular confinement potential. In the presence of a magnetic field, each subband is further split into a series of Landau levels or cyclotron orbits. Upon optical excitation with an infrared source, the intersubband resonances couple to the cyclotron resonance under tilted magnetic fields. This leads to the appearance of satellite peaks at the anti‐crossing point. From the values of splitting at the anti‐crossing points, the spacing between the corresponding subbands can be evaluated. Experiments performed in the absence of magnetic fields demonstrated that it is also possible to directly measure and tune these spacings via density‐chopped infrared transmission spectroscopy. The subband spacings are measured directly and found to be in the far‐infrared region (wide electrically tunable from 6 to 12 meV) of the electromagnetic spectrum. New epitaxial, complementary‐doped, semi‐transparent electrostatic gates that have better optical transmission are introduced [29]. The integrated device with a 2DEG in a high electron mobility transistor structure and artificial metamaterials forms a strongly coupled system that can be electrically driven from an uncoupled to a coupled and again back to the uncoupled regime. In the strongly coupled regime, a periodic exchange of energy between the two systems is observed as a splitting of 0.47 THz at the point of avoided crossing. This is a very high‐energy separation, considering the fact that only one quantum well is employed and thus the achievement of a strong coupling regime can be safely claimed. The tuning mechanism is attributed to the quantum‐confined Stark effect. This device architecture is particularly interesting in designing devices like modulators and detectors specifically in the THz regime. The integrated device has the high‐speed dynamic characteristics of the HEMT design and the appropriate frequency‐controlling ability of the metamaterials. From the design perspective of the metamaterials, they can be made particularly for the THz regime with appropriate dimensions (like the one used in this chapter). Upon excitation with a broadband source, this layer selects the desired frequency for which it is designed, and under the application of an external electrical field across the structure, the transmission of this frequency can be controlled and also modulated. This control dynamics can be very fast, simply owing to the fast dynamics of the HEMT design [52, 53]. Furthermore, this design can also be used to detect THz frequencies. Various other 2D materials (like graphene [54–56] or black phosphorous [57, 58]) are also used these days in the transistor configuration for developing THz detectors, simply utilizing the fast dynamics of the transistor design. These novel devices have thus helped to reduce the long‐debated THz gap in the electromagnetic spectrum, where there is a severe lack of fast electronic devices.