Abstract

In this study, the propagation characteristics of EM surface waves supported by the graphene-coated indium antimonide (InSb) planar waveguide have been investigated theoretically and numerically. The modeling of graphene has been performed by use of Kubo formalism whereas the relative permittivity of indium antimonide has been calculated using Drude’s model. The results for transverse electric (TE) and transverse magnetic (TM) polarized surface waves have been computed analytically and numerically. The major challenge is to model the atomically thick graphene sheet over the InSb grounded slab. To get the temperature-dependent characteristic equation for the electromagnetic surface waves, the surface current boundary conditions’ approach has been employed. The numerical results have been computed for both the TE and TM polarization states and reported that the TE does not support the propagation of surface waves. The dispersion relation, effective mode index, phase speed, propagation length, and field profile have been computed in Mathematica under TM polarization. The graphene and indium antimonide have been found active for low and high Terahertz regions, respectively. As temperature increases, the plasma frequency of the InSb increases due to this reason with the increase of temperature and the resonance frequency, leading to a shift in the dispersion curve. Moreover, with the increase of temperature, the effective wave number of transverse magnetic polarized surface waves also increases. Resultantly, the confinement of such surface waves supported by graphene-loaded InSb increases. It is shown that with the variation of temperature of indium antimonide, the surface waves propagating across the interface can be tuned in the Terahertz region and can be exploited for thermo-optical sensing, near-field communications waveguides, and graphene-based temperature sensor designing.

1. Introduction

In recent years, the study of surface waves has attracted significant attention due to their potential applications in a wide range of fields, including sensing, imaging, waveguides, filters, modulators, and surface communications [18]. Electromagnetic surface waves (ESWs) propagate along the interface between two media, and they can be confined and manipulated at the nanoscale using plasmonic structures. Depending upon the types of partnering surfaces, different types of ESWs can be induced. The ESWs which are controlled by the interface of a metal and any dielectric material having components of polaritons called surface plasmons-polaritons (SPPs) [9, 10]. The EM surface waves are guided by the interfaces of two partnering materials, one with isotropic and other anisotropic dielectric in nature, that type of wave is called Dyakonov surface waves [11]. The electromagnetic surface wave propagation in which at least one of the two interacting materials is often nonhomogenous is perpendicular to the interface guiding surface waves and both interacting materials are isotropic dielectric in nature. That particular type of surface wave is called the Tamm wave [11]. The interfaces of two dielectric materials, where one dielectric material is jointly anisotropic and periodically nonhomogenous in a direction perpendicular to the interface, guide a Dyakonov–Tamm wave, combining the features of Dyakonov wave and Tamm wave. The nonhomogenous electromagnetic surface wave is called the Zenneck or Sommerfeld–Zenneck wave which is incident onto a spherical or planar boundary interface with a complex Brewster angle between two homogeneous media with different dielectric constants [1214]. Electromagnetic surface waves offer efficient and low-loss transmission of information over short distances, along with light confinement in small spaces, making them valuable for compact optical devices such as sensors, filters, and waveguides [15]. Different types of schemes have been proposed by many researchers to control and manipulate the electromagnetic waves of different spectral ranges, i.e., metastructures, near-zero-index materials, and nanostructures [1618]. For the visible range, the different plasmonic platforms have been studied for strong interaction with biomolecules and nanostructures and holding potential for applications such as biosensing, drug delivery, waveguides, and nanophotonic [1921].

To control the THz waves for the surface communication, biomedical sensing, DNA detection, THz antenna, filters, and absorbers, graphene has been realized as a natural metasurface [2225]. Graphene is a monoatomic, transparent sheet of carbon with its exceptional optical, thermal, electrical, and mechanical properties [26]. Graphene, as proposed by Vakil and Engheta, holds immense promise as a material for actively tuning surface waves with unparalleled control. This potential opens up exciting opportunities for applications in surface communication, cloaking, and sensing [27]. Naserpour et al. theoretically proposed the active tunability and cloaking for the graphene-coated nanodimmers and nanowires [28]. Recently, the graphene-based THz metamaterials have been studied for the promising applications, i.e., biosensing, enhanced light emission, and broadband IR photodetector [2932]. The exceptional properties of graphene make it a potential candidate for its unprecedented control and manipulation in Terahertz region. Recently, the graphene-based temperature-assisted devices are the interest of many researchers for their potential applications in the wearable temperature sensors, thermographers, thermal current generation, thermal rectifiers, and IR detectors [3335]. The temperature-sensitive materials (TSM) have temperature-dependent optical, electronic, structural, and phase characteristics; i.e., with the change of external temperature, the characteristics of TSM can be tailored [36]. Indium antimonide (InSb) is realized as temperature-sensitive materials, which electronically has the characteristic of semiconductors formed by the combination of indium (In) and antimony (Sb). Due to its small bandgap, indium antimonide (InSb) growth has received interest as a possible material for infrared detectors and high-speed devices in recent years [37]. InSb-based biosensors can be used to detect microorganisms and are used in temperature-dependent imaging cameras [38]. The response of the graphene towards the thermal agitation is negligible to develop the highly sensitive graphene-based thermal platforms for the communication and sensors, a graphene-loaded indium antimonide (InSb) as a temperature-sensitive material (TSM) has been proposed in the present study. As the surface waves propagate along the interface of two dissimilar media and are sensitive to the electromagnetic characteristics of patterning material, therefore, it is hypothesized that by making an interface of temperature sensitivity, the characteristics of surface may be controlled by the external temperature. Moreover, to tune the surface waves supported by InSb is the biggest challenge and need for the realization of the TSM-based devices. To meet this need, the graphene-coated InSb planar structure has been studied analytically and numerically for the active tuning and manipulation of surface waves from lower to higher THz region. The temperature dependency and effect of graphene parameters on the dispersion curves, effective wave number, and field profiles of surface waves have been studied for both transverse electric (TE) and transverse magnetic (TM) polarizations. The organization of the manuscript is given as the analytical modeling and formulation of the problem being provided in the second section, while the numerical results and discussion and concluding remarks are given in the sections three and four, respectively.

2. Analytical Formulation

In the present research work, the influence of parametric properties, i.e., chemical potential (μ), scattering rate (τ), slab thickness (d), and temperature (T) using dispersion relation, has been analyzed. The propagation of electromagnetic surface waves which travel along the interface of two media of different dielectric constants (i.e., graphene and indium antimonide) is studied. These waves are also known as surface plasmon polariton (SPP) which decays exponentially as the distance from the interface increases. There are three regions which are shown in Figure 1, upper region (z > d) which is free space, the middle region is supposed to filled with the indium antimonide (InSb), i.e., 0 < z < d and the at z = d is the interface between the graphene layer and indium antimonide has been considered, while the lower region (z < 0) whose electrical conductivity is infinite.

To physically model such planar integrated structures, the different approaches, i.e., Dyadic green’s function approach [39], effective medium approach [40], transmission line approach [41], equivalent circuit model approach [42], and analytic absorber model approach [43], have been used in the literature. However, in the present work, the waveguide approach has been used to model the supporting electromagnetic field interaction and propagation of the surface waves supported by the proposed structure, and the proposed approach is quite simple and easy to handle as reported in different works [10, 44].

2.1. TM-Polarized Surface Waves

For transverse magnetic polarized surface waves, the electric field (E) is directed along the interface such as in the x-axis and the magnetic field (H) is directed along the y-axis. The associated field equations are as follows [11, 44].

For region ,

For region ,where and are attenuation constants for their respective regions, while A and B are unknown coefficients, is the exponential decay or damping, is the Terahertz radiation frequency, and is known as the propagation constant. To calculate the unknown constants for the transverse magnetic surface waves, the boundary conditions are applied at . Boundary conditions are given as [27]

By implementing boundary conditions on (1)–(6), the characteristic equation for TM-polarized surface wave is obtained as

2.2. TE-Polarized Surface Waves

For transverse electric polarized surface waves, the field equations are as the propagation of surface waves in which the electric field (E) is taken along y-axis and the magnetic field (H) is directed along x-axis, known as transverse electric TE-polarized surface waves [11, 44].

For region ,

For region ,

To calculate the unknown constants (A&B) for the transverse magnetic surface waves, the boundary conditions are applied at . Boundary conditions are given as

By implementing boundary conditions on (8)–(13), the characteristic equation for the TE-polarized surface wave is obtained as

3. Numerical Results and Discussion

The wave properties of surface waves on graphene-loaded indium antimonide (InSb), i.e., dispersion relation, effective mode index, penetration depth, propagation length, phase velocity, and field profile, have been computed numerically. The characteristics of surface waves under the temperature, thickness of the InSb slab, and chemical potential variations have been studied. The numerical solution for each characteristic has been computed using the contour plot technique. The scattering rate of electron-photon has been taken as . The modeling of graphene monolayer has been carried out [44, 45], using wolfram language-based Mathematica software. The electromagnetic modeling of graphene has been derived using Kubo’s formalism. The optical surface conductivity of graphene is the result of the contribution of two factors, intraband and interband conduction [46]. The intraband conductivity supports graphene as a strong candidature for the tunability of surface waves in the Terahertz regime. It is caused by electron-photon processes. The intraband conductivity [46] is given aswhereas the interband conductivity facilitates graphene as a potential candidate for the tunability of surface waves in the visible and infrared regime. It is caused by direct electron interband transitions. The interband conductivity term [46] is given as

The interband conductivity dominates at THz and infrared frequencies, and graphene gives constant conductivity which is given as . It means that a monolayer graphene sheet can absorb 2.3% of the incident light over a large frequency gap [47]. The total conductivity of graphene can be taken as the sum of intraband and interband conduction, i.e., .where the terms, , , , , and stand for the incident frequency, chemical potential, temperature, Boltzmann constant, reduced Plank’s constant, and electron-photon scattering rate, respectively.

3.1. Surface Waves on Graphene-Covered InSb

The effect of tuning the parameters, i.e., thickness, temperature, and chemical potentials on propagation of surface waves for graphene-covered indium antimonide slab, has been studied and calculated using Mathematica. The amplitude of surface waves decreases as we move away from the indium antimonide slab. With respect to polarization, two types of modes can be discussed for such types of surface waves which are given as(i)TE-polarized surface wave(ii)TM-polarized surface wave

Figure 2 shows the dispersion curves of TE-polarized surface waves with the variation of temperature, i.e., T = 200 K, T = 220 K, T = 240 K, T = 260 K, T = 280 K, and T = 300 K. It is predicted that the TE mode for such type of surface waves facilitated by graphene and indium antimonide structure does not exist since no variation in the straight line has been found on the tuning of temperature. All the lines obtained at various temperature variations follow a light line. So, there exists no solution for such temperature variations. Therefore, it is concluded that TE mode for electromagnetic surface waves does not exist unlike TM mode, which is consistent with the literature as given in [10, 44], while in the subsequent section, the numerical results for the surface wave under TM mode have been given.

3.2. Dispersion Curve Analysis for TM Polarization

In this section, propagation of the transverse magnetic polarized surface for indium antimonide and graphene-coated indium antimonide has been studied graphically and calculated using software pack. The electromagnetic modeling of temperature-sensitive materials such as indium antimonide can be tunable by controlling the temperature values. For Terahertz regime, the relative permittivity of InSb is given by the Drude model [44, 48].

Herein, is the high frequency relative permittivity, is the plasma frequency, and is the damping constant. Plasma frequency can be calculated as ; here, the plasma frequency is affected by the intrinsic carrier density “N,” electronic charge “,” and mass of electron “”, i.e., , , and . The intrinsic carrier density “N” depends upon the temperature “T,” bandgap “,” and Boltzmann constant “” and can be calculated as where the values of band gap “” and Boltzmann constant “” are and As intrinsic carrier density “N” depends upon the temperature “T,” we can say that the plasma frequency is also dependent upon the temperature. Therefore, by tuning the temperature of indium antimonide (InSb), the plasma frequency can be adjusted. Thus, the relative permittivity of InSb can be controlled in Figure 3 by varying the temperature parameters.

The dispersion curves of TM-polarized surface waves supported by graphene-loaded InSb under the variation of temperature, chemical potential, and thickness of InSb slab have been presented graphically in Figure 4. Figure 4(a) depicts that with the increase in temperature, high resonance frequency-polarized surface waves on an InSb slab coated with graphene can be achieved. The impact of chemical potential on the dispersion curve by keeping the temperature constant has been presented in Figure 4(b). It has been found that with the increase of chemical potential, i.e., 0.2 eV, 0.4 eV, 0.6 eV, and 0.8 eV, the resonance frequency decreases. It is obvious from the figure that the chemical potential () affects the conductivity of graphene and the propagation constant starts decreasing with the increase of chemical potential, which shows that for the highly doped graphene layer, the confinement of the surface waves is low. Further, it can be deduced from this result that the skin depth of THz for graphene-coated InSb slab starts decreasing. Due to this reason, the confinement of the surface wave also decreases with the increase of chemical potential. While the effects of the thickness of InSb slab on the propagation of surface waves supported by graphene-coated InSb slab have been depicted in Figure 4(c), it is clear from the figure that with the increase in the thickness d of the InSb slab, the resonance frequency decreases rapidly even at low THz frequency.

The ratio of the propagation constant of medium to free space is known as the effective mode index . The normalized effective wave number determines how surface waves are being confined at the interface and calculated as . On the graphene-coated THz-InSb slab, the confinement of the TM-polarized surface waves has been reported under temperature range, i.e., T = 200 K, 250 K, 300 K, 350 K, and 400 K in Figure 5. It is clear from the figure that with increasing temperature, the effective wave number of transverse magnetic-polarized surface waves also increases. Resultantly, the confinement of such surface waves supported by graphene-loaded indium antimonide THz slab also increases.

To further study the confinement of the surface waves on the graphene supported THz-InSb slab, the influence of thickness of InSb slab on the effective mode index () under the function of frequency has been explained graphically in Figure 6. The possible solution has been calculated using the contour plot technique. Higher values of effective mode index mean higher confinement. It is clear from the figure that thickness of the slab plays an important role in tuning the confinement of surface waves. The effective mode index under the different values of thickness, i.e., d = 20 nm, d = 40 nm, d = 60 nm, d = 80 nm, and d = 100 nm, has been analyzed. It is obvious from the figure that with the increase of thickness of the InSb slab, the effective mode index decreases, and consequently, the confinement of such waves decreases from higher THz frequencies. Figure 7.

Propagation length () is the measurement of the distance travelled by a surface wave while travelling along the interface of two media and is given by the imaginary part of the propagation constant. The effect of temperature on the propagation length has been presented in Figure 8. The propagation length is normalized by wave number . It depicts that propagation length increases with the increase in frequency. It is also concluded that with the increase in temperature, it shifts towards a high-frequency range.

In Figure 9, the influence of temperature on the phase velocity has been shown graphically. Phase velocity is given as . To estimate the phase velocity within the temperature range of T = 200 K to 400 K, it is evident from Figure 9 that the phase velocity decreases as the frequency increases and can be adjusted by varying the temperature.

The comparison among the field profile of electromagnetic surface waves under different temperatures, i.e., T = 200 K and T = 400 K, has been presented in Figure 10. The normalized field profiles of have been taken as a function of transverse distance (z) from the interface. It is obvious from the figure that the field profiles are temperature-sensitive; i.e., high temperature has more intensity as compared to lower temperature. Further, it can be verified from the figure that as the transverse distance from the interface increases, the field profiles follow the exponential decay, which confirms the basic characteristic of surface waves.

In view of the practical implementation, the proposed structure can be easily realizable. The synthesis of graphene can be carried out by the chemical vapor deposition on the indium antimonide substrate [44]; however, the major limitation is the excitation scheme for the surface waves in the THz regime, and to bridge this limitation, many approaches have been studied for the excitation of such evanescent modes in THz, which may be used for the realization of the work [49]. The proposed results have potential applications in the designing of the temperature-assisted THz optical devices for the communication, thermal imaging, and sensing applications.

4. Concluding Remarks

The electromagnetic surface wave propagation supported by graphene-coated indium antimonide (InSb) has been studied in detail. Maxwell’s equations have been solved using surface waves propagating on graphene supported by temperature-sensitive materials, such as indium antimonide, for transverse electric (TE) and transverse magnetic (TM) modes of polarization. Analytical and numerical studies have been conducted on the electromagnetic surface wave propagation on graphene coated with indium antimonide. For mathematical and graphical analysis, the Wolfram Mathematica pack has been utilized. The subsequent conclusions from the numerical results can be derived as the electromagnetic surface waves exit for only TM polarization, while in the case of TE polarization, there is no existence of surface waves. The propagation characteristics of surface waves are found to be temperature-sensitive. The propagation length, penetration depth, effective mode index, and phase velocity of surface waves can be varied by tuning the temperature. It is observed that the resonance frequency of surface waves is found to increase with temperature. The confinement of the surface waves can be tuned by the operating frequency and temperature of InSb. The field profile of the surface wave decreases progressively and decreases as the distance from the interface increases. The proposed results can be used for the realization and future application of the graphene-coated TSM structure for Terahertz applications in designing thermo-optical sensing, THz communication, and near-field thermal imaging platforms. Moreover, the work can be extended for the graphene-based curvilinear temperature-sensitive metasurfaces and bilayer asymmetric temperature-sensitive surfaces for the flexible thermo-optical sensors in THz regime as discussed in [50, 51].

Data Availability

The data used to support the findings of this study are included in the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the researchers supporting project number RSP2024R416, King Saud University, Riyadh, Saudi Arabia.