# An Efficient Semianalytical Modal Analysis of Rectangular Waveguides Containing Metamaterials with Graded Inhomogeneity

Authors: Qianru Weng; Qian Lin; Haifeng Wu
Journal: International Journal of Antennas and Propagation
 (2021)

---

## Abstract

Rectangular waveguides containing inhomogeneous metamaterials with graded refractive-index profiles have potential applications in bending waveguides and radiation-enhanced antennas, and accurate eigenvalue solutions are prerequisite. Commonly used commercial electromagnetic solvers such as HFSS, COMSOL, and CST could not efficiently calculate the eigenvalues of waveguides containing graded refractive-index dielectrics. In this paper, an accurate and efficient semianalytical method based on the modal expansion has been proposed to solve these waveguides. The proposed method has been employed to calculate the eigenvalues, including the cutoff wavenumbers and dispersion relations, for metamaterials with various graded refractive-index profiles. Calculated results are then validated by comparison, using commercial solver HFSS, which indicates the superiority of the proposed method in accuracy and efficiency. Below-cutoff backward wave propagation is observed in waveguides filled with graded refractive-index metamaterials, which provides a new approach for waveguide miniaturization.


---

## Body

## 1. Introduction

Rectangular waveguides filled with inhomogeneous dielectrics are the most fundamental guided wave structures [1] to realize various functional microwave devices, such as filters [2]; couplers [3]; dividers [4]; and substrate integrated waveguides (SIWs) [5, 6]. The overall performances of such devices strongly depend on the guided wave propagation properties of the waveguides. At the same time, the propagation properties of the waveguides depend not only on the physical dimensions of the waveguides but also on the dielectric inhomogeneity in the waveguides.Most of the work reported in the publications before 2000 is concerned about the ordinary dielectrics (or right-handed materials, RHMs). As more and more kinds of metamaterials have been physically realized [7], especially the left-handed materials (LHMs) [8], filling rectangular waveguides with metamaterials has become a new area of research because intriguing propagation phenomena could be achieved. Propagation below the cutoff frequency in a rectangular waveguide partially filled with LHM has been theoretically analyzed and experimentally verified [9, 10]. This phenomenon has been applied in the design of miniaturization filters and antennas. Alternately arranged right-handed materials (RHMs) and left-handed materials (LHMs) constitute a novel kind of metamaterial, which is called composite right-handed and left-handed material (CRLHM) [11]. It has been found that a rectangular waveguide filled with CRLHM supports two propagation bands: a forward wave band (also called right-handed mode) above the cutoff frequency and a backward wave band (also called left-handed mode) below the cutoff frequency [11, 12]. A step-change in the refractive-index within the abovementioned CRLHM constitutes a basic dielectric inhomogeneity in rectangular waveguides. However, it has been already known that filling metamaterials with graded refractive-index profiles permits the manipulation of wave propagation in a waveguide, which could be used to achieve bending waveguides or waveguide splitting effect [13, 14]. Thus, it is desired to acquire the propagation characteristics of a rectangular waveguide filled with CRLHM with a graded refractive-index profile.The graded inhomogeneity of the dielectrics has severely limited the possibility to obtain the analytical eigenvalue solutions of waveguides, and accurate numerical methods need to be used. In recent years, commercial electromagnetic simulation software, including HFSS, COMSOL, and CST, are based on the full-wave numerical methods, such as the finite element method (FEM), the finite-difference time-domain (FDTD), the method of moments (MoM), and the boundary element method (BEM) [15, 16]. Although these electromagnetic solvers have powerful analysis capabilities, they require a very long calculation time and a large amount of computer memory for structures containing complex inhomogeneous dielectrics. Most of the solvers cannot even directly analyze the structures containing gradual spatially varying dielectrics, for example, metamaterials with graded refractive-index profiles.In this paper, an accurate semianalytical method based on a one-dimensional modal expansion theory has been proposed to calculate the eigenvalues of rectangular waveguides filled with CRLHM having graded refractive-index profiles. A similar method has been used to obtain the eigenvalue solutions for dielectric slab waveguides with arbitrary dielectric profiles [17] and optical graded-index planar waveguides with metal cladding [18]. This method is derived from the application of the fundamental Maxwell equations and expands the unknown electromagnetic field into a complete set of orthogonal functions, which leads to a matrix eigenvalue problem whose solution provides all the necessary eigenvalues. The eigenvalues can be applied to obtain the cutoff and dispersion characteristics. The general formulations of the proposed method are described in Section 2. Complete modal analyses, including the longitudinal-section electric (LSE) and longitudinal-section magnetic (LSM) modes, for rectangular waveguides filled with CRLHM having several representative graded refractive-index profiles are conducted using the proposed method in Section 3. Several interesting phenomena, including the less intuitive cutoff dependence on dielectric parameters and the below-cutoff propagation characteristics, are also presented in this section. The accuracy and efficiency of the proposed method are validated by comparison, using the commercial simulation software HFSS in Section 4. Section 5 concludes the development potential of the proposed method and phenomena.


## 2. Theoretical Formulation

Figure1 shows the considered rectangular waveguide filled by composite right-handed and left-handed metamaterial (CRLHM) with a graded refractive-index profile. Such waveguide is homogeneous along the z-axis, and the wave propagation is also along the z-axis. Therefore, the electromagnetic fields have an e-jβz dependence, and β is the propagation constant. The time-harmonic has an ejωt dependence. The Maxwell equations are rearranged in a source-free region.(1a)Δ×Δ×E⟶−ω2μω,xεω,xE⟶=Δμω,xμω,x×Δ×E,(1b)Δ×Δ×H⟶−ω2μω,xεω,xH⟶=Δεω,xεω,x×Δ×H⟶,(1c)Δ⋅εω,xE⟶=∂∂xεω,xEx+εω,x∂∂yEy−jβεω,xEz=0,(1d)Δ⋅μω,xH⟶=∂∂xμω,xHx+,μω,x∂∂yHy−jβμω,xEz=0,where the permittivity ϵ(ω, x) and the permeability μ(ω, x) are the frequency-dependent dielectric permittivity and magnetic permeability, respectively.Figure 1
Rectangular waveguide filled with a metamaterial with a graded refractive-index profile.


### 2.1. The Definition of CRLHM with a Graded Refractive-Index Profile

The permittivity and the permeability of a CRLHM with a graded refractive-index profile can be written as(2a)∈ω,x=∈0∈rω,x=ε0εrω1−hx,(2b)μω,x=μ0μrω,x=μ0μrω1−hx,where εr(ω) and μr(ω) are the frequency-dependent relative permittivity and permeability, which can be described by Drude or Lorentz models [7]; ε0 and μ0 are the free space permittivity and permeability; and h(x) is a spatial-dependent profile function, which can be arbitrary. The sketches of profile functions of graded refractive-index CRLHMs are presented in Figure 2, showing a gradual transition from RHM to LHM [19]. Here, the permittivity ϵ(ω, x) and the permeability μ(ω, x) are assumed to have the identical profile function.Figure 2
Gradual spatially varying effective permittivityϵ(ω, x) and permeability μ(ω, x) of CRLHM with exponential, linear, parabolic, and constant profile functions.In this paper, three kinds of profile functionsh(x) are under consideration:(3)hx=dx,linear profile function,dx2,parabolic profile function,exp−dx,exponential profile function,where d is a positive parameter describing the steepness of the transition from RHM to LHM within a CRLHM. If d = 0, the graded refractive-index CRLHM degenerates into the aforesaid step-changed CRLHM. This specific case will be discussed in the following.


### 2.2. Hybrid LSE and LSM Modes Definition

Standard propagation modes in a rectangular waveguide filled with an inhomogeneous metamaterial as shown in Figure1 are not, in general, either a transverse electric (TE) or a transverse magnetic (TM) mode but a hybrid of TE and TM modes. This suggests that the hybrid modes in such waveguide may be derived from magnetic or electric types of Hertzian potential functions, having single components directed normal to the inhomogeneous direction (x-axis).(1)
From the magnetic Hertzian potential, the solution for a mode that has no component of the electric field normal to thex-axis is obtained. The electric field thus lies in the longitudinal interface (yz-plane) plane, and this mode is referred to as a longitudinal-section electric (LSE) mode [1].(4)Π⟶h=a⟶xϕhx,ye−jβz.(2)
From the electric Hertzian potential, the solution for a mode with no magnetic field component normal to thex-axis is obtained. This mode is called a longitudinal-section magnetic (LSM) mode [1].(5)Π⟶e=a⟶xϕex,ye−jβz.


### 2.3. Electromagnetic Fields and Boundary Conditions

#### 2.3.1. LSE Modes

In a rectangular waveguide filled with inhomogeneous metamaterials, as shown in Figure1, the whole fields of the LSE mode can be derived from Ey [1]. The y-component of Helmholtz’s equation is(6)∇μω,x∇×1μω,x∇×E⟶y+ω2μω,xE⟶y=0.The electric fieldEy could be written as(7)Ey=ϕhxcoslπybe−jβz.For LSE modes,Ey shall vanish on the perfect conducting boundaries (x = 0, a) of the rectangular waveguide. Thus, the potential function ϕh (x) satisfies the boundary condition:(8)ϕhx|x=0,a=0.


#### 2.3.2. LSM Modes

Similarly, the whole fields of the LSM mode can be derived fromHy [1]. The y-component of Helmholtz’s equation is(9)∇∈ω,x∇×1∈ω,x∇×H⟶y+ω2μω,x∈ω,xH⟶y=0.The magnetic fieldHy could be written as(10)Hy=ϕexsinlπybe−jβz.For LSM modes,∂Hy/∂x shall vanish on the perfecting conducting boundaries of the rectangular waveguide at x = 0, a. Thus, the potential function ϕe (x) satisfies the boundary condition:(11)dϕedx|x=0,a=0.


### 2.4. The Semianalytical Method Based on a Modal Expansion Theory

#### 2.4.1. LSE Modes

Substituting equation (7) into equation (6), the potential function ϕh (x) for the LSE mode satisfies the following equation:(12)d2ϕhxdx2−dlnμω,xdxdϕhxdx+ω2μω,x∈ω,x−β2−lπb2ϕhx=0.Equation (12) and the boundary condition (8) constitute a Sturm–Liouville system [1]. According to the modal expansion theory, ϕh (x) could be expanded as a complete set of orthogonal functions.(13)ϕhx=∑m=1NAmgmx,where Am is the amplitude. For the LSE mode, a suitable orthogonal function gm (x) is chosen [1].(14)gmx=sinmπxa;m=1,2,3,….As shown in equation (13), a finite series (N) is used in the derivation of eigenequation. The derivation process is as follows.Firstly, substitute equation (13) into equation (12), and then multiply the resulting equation by the orthogonal function sinnπx/a (n = 1, 2, 3, …);Secondly, integrating over the waveguide width a, equation (12) is transformed to the following:(15)−β2+lπb2An+∑m=1NAm∫0adlnμω,xdxdgmdxgndx+∑m=1NAm∫0aω2μω,x∈ω,xgmgndx=0.Thirdly, equations (15) can be written as a matrix form as follows:(16)ΡA=0,for LSE mode,where A presents a column vector with elements Am and P is a matrix whose elements are defined as follows:(17)pmm=−β2+lπb2+∫0adlnμω,xdxdgmdxgmdx+∫0aω2μω,x∈ω,xgm2dxm=n,pmn=pnm=∫0adlnμω,xdxdgmdxgndx+∫0aω2μω,x∈ω,xgmgndxm≠n.Finally, the eigenequation of LSE modes is obtained:(18)DetP=0.The accuracy of the calculated eigenvalues is closely related to the value of the series number N of the modal expansion in equation (13), and the higher the value of N, the more accurate eigenvalues could be obtained through this method. However, the larger the value of N, the higher the order of the matrix P, and the greater amount of calculation is required to solve the determinant of the matrix P as shown in equation (18). In this paper, the eigenvalue accuracy is characterized by using the change of the eigenvalues, ∆β = |1 ‒ βNβN ‒ 1|%. Through precalculation, it has been found that for rectangular waveguides filled with CRLHM with constant, linear, parabolic, and exponential profiles, ∆β = 0.1% of the LSE mode could be achieved by using N = 40 series number in general. In this paper, the series number N is chosen as 40 to calculate the eigenvalues of the LSE mode.


#### 2.4.2. LSM Modes

Substituting equation (10) into equation (9), the potential function ϕe (x) for the LSM mode satisfies the following equation:(19)d2ϕexdx2−dln∈ω,xdxdϕexdx+ω2μω,x∈ω,x−β2−lπb2ϕex=0.The eigenequation derivation of the LSM mode follows the same line as that of the LSE mode described above. Considering the boundary condition for the LSM mode,ϕex takes the following form:(20)ϕex=∑m=0NBmfmx.For the LSM mode, a suitable functionfm (x) is chosen as follows [1]:(21)fmx=cosmπxa;m=0,1,2,3,….For the LSM mode, substituting equation (20) into equation (19), then multiplying the resulting equation by the orthogonal function cosnπxa (n = 0, 1, 2, 3, …) for the LSM mode, and then integrating over the waveguide width a, equation (19) is transformed to(22)−β2+lπb2Bn+∑m=1NBm∫0adlnεω,xdxdfmdxfndx+∑m=1NBm∫0aω2μω,x∈ω,xfmfndx=0.Equation (22) can be written as a matrix form:(23)QB=0,for LSM mode,where B presents a column vector with elements Bm and Q is a matrix whose elements are defined as(24)qmm=−β2+lπb2+∫0adInεω,xdxdfmdxfmdx+∫0aω2μω,x∈ω,xfm2dxm=n,qmn=qnm=∫0adlnεω,xdxdfmdxfndx+∫0aω2μω,x∈ω,xfmfndxm≠n.Thus, the eigenfunction for LSM modes of a rectangular waveguide filled inhomogeneous metamaterials can be written as(25)DetQ=0.Therefore, the cutoff frequencies and the dispersion relations of the LSM modes could be obtained by calculating eigenvalues of the matrixQ. The accuracy of the calculated eigenvalues of the LSM modes is also closely related to the series number N of the modal expansion in equation (20). From the precalculation, to obtain an error within Δβ = 0.1% for the LSM modes, N=50 series number needs to be chosen to satisfy this accuracy requirement. Thus, in this paper, the series number N is chosen as 50 to calculate the eigenvalues of the LSM mode.


### 2.5. A Specific Case: A Rectangular Waveguide Filled with Homogeneous Dielectric

Forhx=0, a rectangular waveguide filled with graded-index dielectric degenerates into a rectangular waveguide filled with homogeneous dielectric. Thus, the dielectric parameters are changed into(26)εω=ε0εrω,μω=μ0μrω.Substitute equation (26) into the dispersion equations (18) and (25) for LSE and LSM modes. For these two kinds of modes, the dispersion equation becomes(27)β=±k0εrωμrω1−kck02,kc=mπa2+nπb2,which corresponds precisely to the definition of the dispersion constants for TE and TM modes in a homogeneous waveguide.


## 2.1. The Definition of CRLHM with a Graded Refractive-Index Profile

The permittivity and the permeability of a CRLHM with a graded refractive-index profile can be written as(2a)∈ω,x=∈0∈rω,x=ε0εrω1−hx,(2b)μω,x=μ0μrω,x=μ0μrω1−hx,where εr(ω) and μr(ω) are the frequency-dependent relative permittivity and permeability, which can be described by Drude or Lorentz models [7]; ε0 and μ0 are the free space permittivity and permeability; and h(x) is a spatial-dependent profile function, which can be arbitrary. The sketches of profile functions of graded refractive-index CRLHMs are presented in Figure 2, showing a gradual transition from RHM to LHM [19]. Here, the permittivity ϵ(ω, x) and the permeability μ(ω, x) are assumed to have the identical profile function.Figure 2
Gradual spatially varying effective permittivityϵ(ω, x) and permeability μ(ω, x) of CRLHM with exponential, linear, parabolic, and constant profile functions.In this paper, three kinds of profile functionsh(x) are under consideration:(3)hx=dx,linear profile function,dx2,parabolic profile function,exp−dx,exponential profile function,where d is a positive parameter describing the steepness of the transition from RHM to LHM within a CRLHM. If d = 0, the graded refractive-index CRLHM degenerates into the aforesaid step-changed CRLHM. This specific case will be discussed in the following.


## 2.2. Hybrid LSE and LSM Modes Definition

Standard propagation modes in a rectangular waveguide filled with an inhomogeneous metamaterial as shown in Figure1 are not, in general, either a transverse electric (TE) or a transverse magnetic (TM) mode but a hybrid of TE and TM modes. This suggests that the hybrid modes in such waveguide may be derived from magnetic or electric types of Hertzian potential functions, having single components directed normal to the inhomogeneous direction (x-axis).(1)
From the magnetic Hertzian potential, the solution for a mode that has no component of the electric field normal to thex-axis is obtained. The electric field thus lies in the longitudinal interface (yz-plane) plane, and this mode is referred to as a longitudinal-section electric (LSE) mode [1].(4)Π⟶h=a⟶xϕhx,ye−jβz.(2)
From the electric Hertzian potential, the solution for a mode with no magnetic field component normal to thex-axis is obtained. This mode is called a longitudinal-section magnetic (LSM) mode [1].(5)Π⟶e=a⟶xϕex,ye−jβz.


## 2.3. Electromagnetic Fields and Boundary Conditions

### 2.3.1. LSE Modes

In a rectangular waveguide filled with inhomogeneous metamaterials, as shown in Figure1, the whole fields of the LSE mode can be derived from Ey [1]. The y-component of Helmholtz’s equation is(6)∇μω,x∇×1μω,x∇×E⟶y+ω2μω,xE⟶y=0.The electric fieldEy could be written as(7)Ey=ϕhxcoslπybe−jβz.For LSE modes,Ey shall vanish on the perfect conducting boundaries (x = 0, a) of the rectangular waveguide. Thus, the potential function ϕh (x) satisfies the boundary condition:(8)ϕhx|x=0,a=0.


### 2.3.2. LSM Modes

Similarly, the whole fields of the LSM mode can be derived fromHy [1]. The y-component of Helmholtz’s equation is(9)∇∈ω,x∇×1∈ω,x∇×H⟶y+ω2μω,x∈ω,xH⟶y=0.The magnetic fieldHy could be written as(10)Hy=ϕexsinlπybe−jβz.For LSM modes,∂Hy/∂x shall vanish on the perfecting conducting boundaries of the rectangular waveguide at x = 0, a. Thus, the potential function ϕe (x) satisfies the boundary condition:(11)dϕedx|x=0,a=0.


## 2.3.1. LSE Modes

In a rectangular waveguide filled with inhomogeneous metamaterials, as shown in Figure1, the whole fields of the LSE mode can be derived from Ey [1]. The y-component of Helmholtz’s equation is(6)∇μω,x∇×1μω,x∇×E⟶y+ω2μω,xE⟶y=0.The electric fieldEy could be written as(7)Ey=ϕhxcoslπybe−jβz.For LSE modes,Ey shall vanish on the perfect conducting boundaries (x = 0, a) of the rectangular waveguide. Thus, the potential function ϕh (x) satisfies the boundary condition:(8)ϕhx|x=0,a=0.


## 2.3.2. LSM Modes

Similarly, the whole fields of the LSM mode can be derived fromHy [1]. The y-component of Helmholtz’s equation is(9)∇∈ω,x∇×1∈ω,x∇×H⟶y+ω2μω,x∈ω,xH⟶y=0.The magnetic fieldHy could be written as(10)Hy=ϕexsinlπybe−jβz.For LSM modes,∂Hy/∂x shall vanish on the perfecting conducting boundaries of the rectangular waveguide at x = 0, a. Thus, the potential function ϕe (x) satisfies the boundary condition:(11)dϕedx|x=0,a=0.


## 2.4. The Semianalytical Method Based on a Modal Expansion Theory

### 2.4.1. LSE Modes

Substituting equation (7) into equation (6), the potential function ϕh (x) for the LSE mode satisfies the following equation:(12)d2ϕhxdx2−dlnμω,xdxdϕhxdx+ω2μω,x∈ω,x−β2−lπb2ϕhx=0.Equation (12) and the boundary condition (8) constitute a Sturm–Liouville system [1]. According to the modal expansion theory, ϕh (x) could be expanded as a complete set of orthogonal functions.(13)ϕhx=∑m=1NAmgmx,where Am is the amplitude. For the LSE mode, a suitable orthogonal function gm (x) is chosen [1].(14)gmx=sinmπxa;m=1,2,3,….As shown in equation (13), a finite series (N) is used in the derivation of eigenequation. The derivation process is as follows.Firstly, substitute equation (13) into equation (12), and then multiply the resulting equation by the orthogonal function sinnπx/a (n = 1, 2, 3, …);Secondly, integrating over the waveguide width a, equation (12) is transformed to the following:(15)−β2+lπb2An+∑m=1NAm∫0adlnμω,xdxdgmdxgndx+∑m=1NAm∫0aω2μω,x∈ω,xgmgndx=0.Thirdly, equations (15) can be written as a matrix form as follows:(16)ΡA=0,for LSE mode,where A presents a column vector with elements Am and P is a matrix whose elements are defined as follows:(17)pmm=−β2+lπb2+∫0adlnμω,xdxdgmdxgmdx+∫0aω2μω,x∈ω,xgm2dxm=n,pmn=pnm=∫0adlnμω,xdxdgmdxgndx+∫0aω2μω,x∈ω,xgmgndxm≠n.Finally, the eigenequation of LSE modes is obtained:(18)DetP=0.The accuracy of the calculated eigenvalues is closely related to the value of the series number N of the modal expansion in equation (13), and the higher the value of N, the more accurate eigenvalues could be obtained through this method. However, the larger the value of N, the higher the order of the matrix P, and the greater amount of calculation is required to solve the determinant of the matrix P as shown in equation (18). In this paper, the eigenvalue accuracy is characterized by using the change of the eigenvalues, ∆β = |1 ‒ βNβN ‒ 1|%. Through precalculation, it has been found that for rectangular waveguides filled with CRLHM with constant, linear, parabolic, and exponential profiles, ∆β = 0.1% of the LSE mode could be achieved by using N = 40 series number in general. In this paper, the series number N is chosen as 40 to calculate the eigenvalues of the LSE mode.


### 2.4.2. LSM Modes

Substituting equation (10) into equation (9), the potential function ϕe (x) for the LSM mode satisfies the following equation:(19)d2ϕexdx2−dln∈ω,xdxdϕexdx+ω2μω,x∈ω,x−β2−lπb2ϕex=0.The eigenequation derivation of the LSM mode follows the same line as that of the LSE mode described above. Considering the boundary condition for the LSM mode,ϕex takes the following form:(20)ϕex=∑m=0NBmfmx.For the LSM mode, a suitable functionfm (x) is chosen as follows [1]:(21)fmx=cosmπxa;m=0,1,2,3,….For the LSM mode, substituting equation (20) into equation (19), then multiplying the resulting equation by the orthogonal function cosnπxa (n = 0, 1, 2, 3, …) for the LSM mode, and then integrating over the waveguide width a, equation (19) is transformed to(22)−β2+lπb2Bn+∑m=1NBm∫0adlnεω,xdxdfmdxfndx+∑m=1NBm∫0aω2μω,x∈ω,xfmfndx=0.Equation (22) can be written as a matrix form:(23)QB=0,for LSM mode,where B presents a column vector with elements Bm and Q is a matrix whose elements are defined as(24)qmm=−β2+lπb2+∫0adInεω,xdxdfmdxfmdx+∫0aω2μω,x∈ω,xfm2dxm=n,qmn=qnm=∫0adlnεω,xdxdfmdxfndx+∫0aω2μω,x∈ω,xfmfndxm≠n.Thus, the eigenfunction for LSM modes of a rectangular waveguide filled inhomogeneous metamaterials can be written as(25)DetQ=0.Therefore, the cutoff frequencies and the dispersion relations of the LSM modes could be obtained by calculating eigenvalues of the matrixQ. The accuracy of the calculated eigenvalues of the LSM modes is also closely related to the series number N of the modal expansion in equation (20). From the precalculation, to obtain an error within Δβ = 0.1% for the LSM modes, N=50 series number needs to be chosen to satisfy this accuracy requirement. Thus, in this paper, the series number N is chosen as 50 to calculate the eigenvalues of the LSM mode.


## 2.4.1. LSE Modes

Substituting equation (7) into equation (6), the potential function ϕh (x) for the LSE mode satisfies the following equation:(12)d2ϕhxdx2−dlnμω,xdxdϕhxdx+ω2μω,x∈ω,x−β2−lπb2ϕhx=0.Equation (12) and the boundary condition (8) constitute a Sturm–Liouville system [1]. According to the modal expansion theory, ϕh (x) could be expanded as a complete set of orthogonal functions.(13)ϕhx=∑m=1NAmgmx,where Am is the amplitude. For the LSE mode, a suitable orthogonal function gm (x) is chosen [1].(14)gmx=sinmπxa;m=1,2,3,….As shown in equation (13), a finite series (N) is used in the derivation of eigenequation. The derivation process is as follows.Firstly, substitute equation (13) into equation (12), and then multiply the resulting equation by the orthogonal function sinnπx/a (n = 1, 2, 3, …);Secondly, integrating over the waveguide width a, equation (12) is transformed to the following:(15)−β2+lπb2An+∑m=1NAm∫0adlnμω,xdxdgmdxgndx+∑m=1NAm∫0aω2μω,x∈ω,xgmgndx=0.Thirdly, equations (15) can be written as a matrix form as follows:(16)ΡA=0,for LSE mode,where A presents a column vector with elements Am and P is a matrix whose elements are defined as follows:(17)pmm=−β2+lπb2+∫0adlnμω,xdxdgmdxgmdx+∫0aω2μω,x∈ω,xgm2dxm=n,pmn=pnm=∫0adlnμω,xdxdgmdxgndx+∫0aω2μω,x∈ω,xgmgndxm≠n.Finally, the eigenequation of LSE modes is obtained:(18)DetP=0.The accuracy of the calculated eigenvalues is closely related to the value of the series number N of the modal expansion in equation (13), and the higher the value of N, the more accurate eigenvalues could be obtained through this method. However, the larger the value of N, the higher the order of the matrix P, and the greater amount of calculation is required to solve the determinant of the matrix P as shown in equation (18). In this paper, the eigenvalue accuracy is characterized by using the change of the eigenvalues, ∆β = |1 ‒ βNβN ‒ 1|%. Through precalculation, it has been found that for rectangular waveguides filled with CRLHM with constant, linear, parabolic, and exponential profiles, ∆β = 0.1% of the LSE mode could be achieved by using N = 40 series number in general. In this paper, the series number N is chosen as 40 to calculate the eigenvalues of the LSE mode.


## 2.4.2. LSM Modes

Substituting equation (10) into equation (9), the potential function ϕe (x) for the LSM mode satisfies the following equation:(19)d2ϕexdx2−dln∈ω,xdxdϕexdx+ω2μω,x∈ω,x−β2−lπb2ϕex=0.The eigenequation derivation of the LSM mode follows the same line as that of the LSE mode described above. Considering the boundary condition for the LSM mode,ϕex takes the following form:(20)ϕex=∑m=0NBmfmx.For the LSM mode, a suitable functionfm (x) is chosen as follows [1]:(21)fmx=cosmπxa;m=0,1,2,3,….For the LSM mode, substituting equation (20) into equation (19), then multiplying the resulting equation by the orthogonal function cosnπxa (n = 0, 1, 2, 3, …) for the LSM mode, and then integrating over the waveguide width a, equation (19) is transformed to(22)−β2+lπb2Bn+∑m=1NBm∫0adlnεω,xdxdfmdxfndx+∑m=1NBm∫0aω2μω,x∈ω,xfmfndx=0.Equation (22) can be written as a matrix form:(23)QB=0,for LSM mode,where B presents a column vector with elements Bm and Q is a matrix whose elements are defined as(24)qmm=−β2+lπb2+∫0adInεω,xdxdfmdxfmdx+∫0aω2μω,x∈ω,xfm2dxm=n,qmn=qnm=∫0adlnεω,xdxdfmdxfndx+∫0aω2μω,x∈ω,xfmfndxm≠n.Thus, the eigenfunction for LSM modes of a rectangular waveguide filled inhomogeneous metamaterials can be written as(25)DetQ=0.Therefore, the cutoff frequencies and the dispersion relations of the LSM modes could be obtained by calculating eigenvalues of the matrixQ. The accuracy of the calculated eigenvalues of the LSM modes is also closely related to the series number N of the modal expansion in equation (20). From the precalculation, to obtain an error within Δβ = 0.1% for the LSM modes, N=50 series number needs to be chosen to satisfy this accuracy requirement. Thus, in this paper, the series number N is chosen as 50 to calculate the eigenvalues of the LSM mode.


## 2.5. A Specific Case: A Rectangular Waveguide Filled with Homogeneous Dielectric

Forhx=0, a rectangular waveguide filled with graded-index dielectric degenerates into a rectangular waveguide filled with homogeneous dielectric. Thus, the dielectric parameters are changed into(26)εω=ε0εrω,μω=μ0μrω.Substitute equation (26) into the dispersion equations (18) and (25) for LSE and LSM modes. For these two kinds of modes, the dispersion equation becomes(27)β=±k0εrωμrω1−kck02,kc=mπa2+nπb2,which corresponds precisely to the definition of the dispersion constants for TE and TM modes in a homogeneous waveguide.


## 3. Results and Discussion

From the waveguide modal theory presented in the previous section, a specific code that can compute the eigenvalues of a rectangular waveguide filled with a graded-index CRLHM has been developed. In this section, the cutoff characteristics and dispersion diagrams are obtained from the calculated eigenvalues. In order to extend the results for a general use, the solution data in this paper are normalized for all length scales.


### 3.1. Cutoff Characteristics

The fundamental LSE mode in an inhomogeneous rectangular waveguide is the LSE10 mode, while the fundamental LSM mode is the LSM01 mode. The relationships between the normalized cutoff wavenumber kca and the profile parameter d for the LSE10 mode and LSM01 mode are shown in Figures 3(a) and 3(b), respectively. Notably, d represents the slope of the profile function. The higher the value of d is, the larger the area is occupied by LHM. The representative waveguide has physical dimensions a = 1 and b = 0.5.Figure 3
The normalized cutoff wavenumberskca versus the profile parameter d for various graded-index profile functions, a = 1, b = 0.5: (a) LSE10 mode and (b) LSM01 mode.
(a)(b)As shown in Figure3(a), when d=0, the rectangular waveguide is filled with a homogeneous RHM. Thus, the normalized cutoff wave number of the LSE10 mode is kca=3.14. In the region of d<2, RHM occupies a larger area, and kca increases as the value of d increases (i.e., the occupied area of RHM increases). This behavior follows the well-known result that the cutoff wavenumber of a rectangular waveguide increases as the filling area of the RHM decreases [1].In the region ofd>2, the graded-index LHM occupies the larger area, and a larger profile parameter d in this range produces an increased kca; as d⟶2, kca of LSE10 mode is dramatically lowered which is very different from the cutoff dependence of RHM, suggesting that this type of rectangular waveguide may support propagation at very low wavenumbers.In Figure3(a), the results for constant, linear, parabolic, and exponential profile functions are presented. The shapes of these curves are similar in general. A constant profile, which means filling a rectangular waveguide with homogeneous RHM and homogeneous LHM (abrupt transitions between LHM and RHM), has exact analytical solutions [11]. In Figure 3(a), for the constant profile, the approximate solution using the proposed technique (green dot-dash line) is compared with the exact analytical solution (empty black circles). These two solutions agree well with each other. In Figure 3(b), the relationship between the normalized cutoff wavenumbers kca versus d for the LSM01 modes is presented.In the region ofd⟶2, the cutoff wavenumbers can be dramatically lowered, and the transversal dimension of the designed rectangular waveguides working in this region can be extremely small. Thus, operation in these frequency-reduced regions of the LSE10 and LSM01 modes offers the potential for the miniaturization of rectangular waveguide components.


### 3.2. Dispersion Relations

#### 3.2.1. LSE Modes

For the chosen waveguide dimensions,a = 1 and b = 0.5, it can be observed from Figure 3(a) that a profile parameter of d=2.1 will result in reduced normalized cutoff wavenumbers kca=0.45,0.51,0.75,and0.92 for the exponential, linear, parabolic, and constant profiles of the LSE10 mode, respectively. In Figure 4, the dispersion relations for various profile functions of the fundamental LSE10 mode and the higher-order LSE20 mode are presented. Here, the curves are represented by different symbols, indicating the normalized propagation characteristics βa for different profiles. The TE10 mode dispersion relation of a hollow rectangular waveguide with the same dimensions is also shown for comparison, where the dispersion curve is represented using solid gray circles. Interestingly, the rectangular waveguide filled with graded-index CRLHM supports two propagation bands: a forward wave band (∂ω/∂β>0) and a back wave band (∂ω/∂β>0) below the cutoff frequency. There is a stopband between these two propagation bands, which is influenced by the profile functions. The backward wave propagation band is significantly below the fundamental TE10 mode of a hollow rectangular waveguide. The analytical dispersion solutions have been calculated for a rectangular waveguide filled by the CRLHM with a constant profile, as indicated in Figure 4, and good agreement is observed between the analytical solution and the approximate solution presented in this paper.Figure 4
The dispersion relations of the CRLHM-filled rectangular waveguide’s LSE modes for exponential, linear, parabolic, and constant profiles and the dispersion relation of the TE10 modes of a hollow rectangular waveguide with the same dimension, d=2.1.It is known that the dispersion relations of the guided waves are dependent on the profile parameterd. In Figure 5, the dispersion relations of the fundamental LSE10 mode and the higher-order LSE20 mode are shown for the exponential, linear, parabolic, and constant profile functions in the region where the value of d is near the frequency-reduced region. In Figure 5, backward wave propagation of the LSE10 mode below the cutoff frequency could be observed for the four profile functions, while the LSE20 modes support forward wave propagation. Some dispersion curves of the LSE20 modes in Figure 5, such as the dispersion curve for d=3.3 of the exponential profile and the dispersion curves for d=2.9and3.3 of the parabolic profile, exhibit negative group velocities (∂ω/∂β<0), which could be observed in metallic waveguides containing metamaterials [11].Figure 5
The dispersion relations of the LSE modes with various values of the profile parameterd and the step width of Δd=0.4: (a) exponential profile; (b) linear profile; (c) parabolic profile; (d) constant profile.
(a)(b)(c)(d)From Figure5, it can be found that the bandwidths of the stopband are affected by the profile functions and the profile parameter d. To clearly illustrate the parameter dependence of the stopband bandwidths, Figure 6 shows the stopband bandwidths versus the profile parameter d for different profile functions: constant, linear, parabolic, and exponential. The results show that the constant profile has the widest stopband bandwidth and increases as the parameter d increases. For the exponential, linear, and parabolic profiles, the bandwidths decrease as d increases at first and then increase as d increases.Figure 6
The stopband bandwidths between the backward LSE10 mode and the forward LSE20 mode versus the profile parameter d for the constant, linear, parabolic, and exponential profiles.


#### 3.2.2. LSM Modes

Figure7 presents the dispersion relations of the fundamental LSM01 mode and the higher-order LSM11 mode for various profile functions when d=2.1. The TM01 mode dispersion of a hollow rectangular waveguide of the same dimension is also plotted for comparison. Both of LSM01 and LSM11 modes propagate forward (∂ω/∂β>0), and the curve slope of the LSM01 mode is larger than that of the LSM11 mode. The cutoff wavenumbers of the higher-order LSM11 mode in the CRLHM-filled rectangular waveguide are near the cutoff wavenumber of the TM01 mode of a hollow rectangular waveguide. Figure 8 shows the dispersion relations for different profile functions versus different values of d. The negative group velocities (∂ω/∂β>0) could be observed in a constant profile.Figure 7
The dispersion relations of LSM modes for constant, linear, parabolic, and exponential profiles and the dispersion of the TM01 modes of a hollow rectangular waveguide with the same dimension.Figure 8
The dispersion relations of the LSM modes with various values of the profile parameterd and step width of Δd=0.4: (a) exponential profile; (b) linear profile; (c) parabolic profile; (d) constant profile.
(a)(b)(c)(d)


## 3.1. Cutoff Characteristics

The fundamental LSE mode in an inhomogeneous rectangular waveguide is the LSE10 mode, while the fundamental LSM mode is the LSM01 mode. The relationships between the normalized cutoff wavenumber kca and the profile parameter d for the LSE10 mode and LSM01 mode are shown in Figures 3(a) and 3(b), respectively. Notably, d represents the slope of the profile function. The higher the value of d is, the larger the area is occupied by LHM. The representative waveguide has physical dimensions a = 1 and b = 0.5.Figure 3
The normalized cutoff wavenumberskca versus the profile parameter d for various graded-index profile functions, a = 1, b = 0.5: (a) LSE10 mode and (b) LSM01 mode.
(a)(b)As shown in Figure3(a), when d=0, the rectangular waveguide is filled with a homogeneous RHM. Thus, the normalized cutoff wave number of the LSE10 mode is kca=3.14. In the region of d<2, RHM occupies a larger area, and kca increases as the value of d increases (i.e., the occupied area of RHM increases). This behavior follows the well-known result that the cutoff wavenumber of a rectangular waveguide increases as the filling area of the RHM decreases [1].In the region ofd>2, the graded-index LHM occupies the larger area, and a larger profile parameter d in this range produces an increased kca; as d⟶2, kca of LSE10 mode is dramatically lowered which is very different from the cutoff dependence of RHM, suggesting that this type of rectangular waveguide may support propagation at very low wavenumbers.In Figure3(a), the results for constant, linear, parabolic, and exponential profile functions are presented. The shapes of these curves are similar in general. A constant profile, which means filling a rectangular waveguide with homogeneous RHM and homogeneous LHM (abrupt transitions between LHM and RHM), has exact analytical solutions [11]. In Figure 3(a), for the constant profile, the approximate solution using the proposed technique (green dot-dash line) is compared with the exact analytical solution (empty black circles). These two solutions agree well with each other. In Figure 3(b), the relationship between the normalized cutoff wavenumbers kca versus d for the LSM01 modes is presented.In the region ofd⟶2, the cutoff wavenumbers can be dramatically lowered, and the transversal dimension of the designed rectangular waveguides working in this region can be extremely small. Thus, operation in these frequency-reduced regions of the LSE10 and LSM01 modes offers the potential for the miniaturization of rectangular waveguide components.


## 3.2. Dispersion Relations

### 3.2.1. LSE Modes

For the chosen waveguide dimensions,a = 1 and b = 0.5, it can be observed from Figure 3(a) that a profile parameter of d=2.1 will result in reduced normalized cutoff wavenumbers kca=0.45,0.51,0.75,and0.92 for the exponential, linear, parabolic, and constant profiles of the LSE10 mode, respectively. In Figure 4, the dispersion relations for various profile functions of the fundamental LSE10 mode and the higher-order LSE20 mode are presented. Here, the curves are represented by different symbols, indicating the normalized propagation characteristics βa for different profiles. The TE10 mode dispersion relation of a hollow rectangular waveguide with the same dimensions is also shown for comparison, where the dispersion curve is represented using solid gray circles. Interestingly, the rectangular waveguide filled with graded-index CRLHM supports two propagation bands: a forward wave band (∂ω/∂β>0) and a back wave band (∂ω/∂β>0) below the cutoff frequency. There is a stopband between these two propagation bands, which is influenced by the profile functions. The backward wave propagation band is significantly below the fundamental TE10 mode of a hollow rectangular waveguide. The analytical dispersion solutions have been calculated for a rectangular waveguide filled by the CRLHM with a constant profile, as indicated in Figure 4, and good agreement is observed between the analytical solution and the approximate solution presented in this paper.Figure 4
The dispersion relations of the CRLHM-filled rectangular waveguide’s LSE modes for exponential, linear, parabolic, and constant profiles and the dispersion relation of the TE10 modes of a hollow rectangular waveguide with the same dimension, d=2.1.It is known that the dispersion relations of the guided waves are dependent on the profile parameterd. In Figure 5, the dispersion relations of the fundamental LSE10 mode and the higher-order LSE20 mode are shown for the exponential, linear, parabolic, and constant profile functions in the region where the value of d is near the frequency-reduced region. In Figure 5, backward wave propagation of the LSE10 mode below the cutoff frequency could be observed for the four profile functions, while the LSE20 modes support forward wave propagation. Some dispersion curves of the LSE20 modes in Figure 5, such as the dispersion curve for d=3.3 of the exponential profile and the dispersion curves for d=2.9and3.3 of the parabolic profile, exhibit negative group velocities (∂ω/∂β<0), which could be observed in metallic waveguides containing metamaterials [11].Figure 5
The dispersion relations of the LSE modes with various values of the profile parameterd and the step width of Δd=0.4: (a) exponential profile; (b) linear profile; (c) parabolic profile; (d) constant profile.
(a)(b)(c)(d)From Figure5, it can be found that the bandwidths of the stopband are affected by the profile functions and the profile parameter d. To clearly illustrate the parameter dependence of the stopband bandwidths, Figure 6 shows the stopband bandwidths versus the profile parameter d for different profile functions: constant, linear, parabolic, and exponential. The results show that the constant profile has the widest stopband bandwidth and increases as the parameter d increases. For the exponential, linear, and parabolic profiles, the bandwidths decrease as d increases at first and then increase as d increases.Figure 6
The stopband bandwidths between the backward LSE10 mode and the forward LSE20 mode versus the profile parameter d for the constant, linear, parabolic, and exponential profiles.


### 3.2.2. LSM Modes

Figure7 presents the dispersion relations of the fundamental LSM01 mode and the higher-order LSM11 mode for various profile functions when d=2.1. The TM01 mode dispersion of a hollow rectangular waveguide of the same dimension is also plotted for comparison. Both of LSM01 and LSM11 modes propagate forward (∂ω/∂β>0), and the curve slope of the LSM01 mode is larger than that of the LSM11 mode. The cutoff wavenumbers of the higher-order LSM11 mode in the CRLHM-filled rectangular waveguide are near the cutoff wavenumber of the TM01 mode of a hollow rectangular waveguide. Figure 8 shows the dispersion relations for different profile functions versus different values of d. The negative group velocities (∂ω/∂β>0) could be observed in a constant profile.Figure 7
The dispersion relations of LSM modes for constant, linear, parabolic, and exponential profiles and the dispersion of the TM01 modes of a hollow rectangular waveguide with the same dimension.Figure 8
The dispersion relations of the LSM modes with various values of the profile parameterd and step width of Δd=0.4: (a) exponential profile; (b) linear profile; (c) parabolic profile; (d) constant profile.
(a)(b)(c)(d)


## 3.2.1. LSE Modes

For the chosen waveguide dimensions,a = 1 and b = 0.5, it can be observed from Figure 3(a) that a profile parameter of d=2.1 will result in reduced normalized cutoff wavenumbers kca=0.45,0.51,0.75,and0.92 for the exponential, linear, parabolic, and constant profiles of the LSE10 mode, respectively. In Figure 4, the dispersion relations for various profile functions of the fundamental LSE10 mode and the higher-order LSE20 mode are presented. Here, the curves are represented by different symbols, indicating the normalized propagation characteristics βa for different profiles. The TE10 mode dispersion relation of a hollow rectangular waveguide with the same dimensions is also shown for comparison, where the dispersion curve is represented using solid gray circles. Interestingly, the rectangular waveguide filled with graded-index CRLHM supports two propagation bands: a forward wave band (∂ω/∂β>0) and a back wave band (∂ω/∂β>0) below the cutoff frequency. There is a stopband between these two propagation bands, which is influenced by the profile functions. The backward wave propagation band is significantly below the fundamental TE10 mode of a hollow rectangular waveguide. The analytical dispersion solutions have been calculated for a rectangular waveguide filled by the CRLHM with a constant profile, as indicated in Figure 4, and good agreement is observed between the analytical solution and the approximate solution presented in this paper.Figure 4
The dispersion relations of the CRLHM-filled rectangular waveguide’s LSE modes for exponential, linear, parabolic, and constant profiles and the dispersion relation of the TE10 modes of a hollow rectangular waveguide with the same dimension, d=2.1.It is known that the dispersion relations of the guided waves are dependent on the profile parameterd. In Figure 5, the dispersion relations of the fundamental LSE10 mode and the higher-order LSE20 mode are shown for the exponential, linear, parabolic, and constant profile functions in the region where the value of d is near the frequency-reduced region. In Figure 5, backward wave propagation of the LSE10 mode below the cutoff frequency could be observed for the four profile functions, while the LSE20 modes support forward wave propagation. Some dispersion curves of the LSE20 modes in Figure 5, such as the dispersion curve for d=3.3 of the exponential profile and the dispersion curves for d=2.9and3.3 of the parabolic profile, exhibit negative group velocities (∂ω/∂β<0), which could be observed in metallic waveguides containing metamaterials [11].Figure 5
The dispersion relations of the LSE modes with various values of the profile parameterd and the step width of Δd=0.4: (a) exponential profile; (b) linear profile; (c) parabolic profile; (d) constant profile.
(a)(b)(c)(d)From Figure5, it can be found that the bandwidths of the stopband are affected by the profile functions and the profile parameter d. To clearly illustrate the parameter dependence of the stopband bandwidths, Figure 6 shows the stopband bandwidths versus the profile parameter d for different profile functions: constant, linear, parabolic, and exponential. The results show that the constant profile has the widest stopband bandwidth and increases as the parameter d increases. For the exponential, linear, and parabolic profiles, the bandwidths decrease as d increases at first and then increase as d increases.Figure 6
The stopband bandwidths between the backward LSE10 mode and the forward LSE20 mode versus the profile parameter d for the constant, linear, parabolic, and exponential profiles.


## 3.2.2. LSM Modes

Figure7 presents the dispersion relations of the fundamental LSM01 mode and the higher-order LSM11 mode for various profile functions when d=2.1. The TM01 mode dispersion of a hollow rectangular waveguide of the same dimension is also plotted for comparison. Both of LSM01 and LSM11 modes propagate forward (∂ω/∂β>0), and the curve slope of the LSM01 mode is larger than that of the LSM11 mode. The cutoff wavenumbers of the higher-order LSM11 mode in the CRLHM-filled rectangular waveguide are near the cutoff wavenumber of the TM01 mode of a hollow rectangular waveguide. Figure 8 shows the dispersion relations for different profile functions versus different values of d. The negative group velocities (∂ω/∂β>0) could be observed in a constant profile.Figure 7
The dispersion relations of LSM modes for constant, linear, parabolic, and exponential profiles and the dispersion of the TM01 modes of a hollow rectangular waveguide with the same dimension.Figure 8
The dispersion relations of the LSM modes with various values of the profile parameterd and step width of Δd=0.4: (a) exponential profile; (b) linear profile; (c) parabolic profile; (d) constant profile.
(a)(b)(c)(d)


## 4. Full-Wave Simulation Validation

Validation of the proposed numerical method has been conducted by the commercial finite element method (FEM) solver HFSS. HFSS is based on the discretization of the electromagnetic structures which has very high analysis capability, but the discretization is very time- and memory-consuming. Since the spatial-dependent permittivity and permeability for the graded-index metamaterials cannot be set in HFSS, discretization on the graded-index metamaterial is conducted. As shown in the inset of Figure9, the studied region 0,a is divided into n subregions, and the subregion length is xi−xi−1. The division of subintervals can be equal or unequal. If the value of n is large enough, the permittivity and the permeability of each subinterval can be viewed as constants. The constants could be replaced by the left endpoint’s values, the right endpoint, or the midpoint. Here, the values of the right endpoints are employed.Figure 9
Convergence analysis conducted by HFSS and the inset presents the treatment of the discretization of the graded-index metamaterials.The accuracy of the HFSS simulation results is directly related to the number of the discretized subareas for the graded-index metamaterials. Convergence analysis for the subregion number is shown in Figure9 for a metamaterial with a parabolic profile. The horizontal axis is the iteration step number, and the vertical axis is the change in the propagation constant β, which is defined as Δβ=1−βs/βs−1%. For each iteration, the length of the subinterval is half of the previous iteration. The starting discretization of the graded-index metamaterials is ten subregions, and convergence to 0.1% on Δβ is achieved in 7 steps (using the eigenmode simulation).Compared with the results obtained by the proposed numerical, the convergent results of the propagation constants for the graded-index metamaterials are compared. The dispersion diagrams for constant and linear dielectric profile functions are proposed in Figure10, and the results agree well with each other.Figure 10
Dispersion diagrams of a rectangular waveguide for constant and parabolic profile functions; the curves illustrated with green or red lines are obtained with the proposed technique; curves with green or red circles are obtained with HFSS simulation.Although the comparisons between the proposed numerical and the HFSS agree well with each other, the consumed calculation time is very different. In this paper, the code of the proposed method and the HFSS simulation are both performed by using a regular PC with an Intel Corei7-8500 CPU at 1.80 GHz (eight CPUs) and 16 GB of RAM for rectangular waveguide filled with graded-index CRLHM. The computation time comparison of the proposed method and the HFSS simulation is presented in Table1. 30 points of β have been computed in both methods. From Table 1, for constant profile function, the calculation time of the proposed technique needs overall tens of seconds to compute the eigenvalues, while the HFSS simulation time is about five minutes. For parabolic profile functions, not considering the time-consuming of the convergence analysis, the HFSS simulation time is up to twenty minutes, while the calculation time of the proposed techniques only takes up to nearly one minute. Therefore, the numerical method proposed in this paper has both accuracy and efficiency.Table 1
Computation time comparison of the semianalytical method and the HFSS simulation.
Semianalytical methodHFSS simulationLSE10LSE20LSE10 (s)LSE20 (s)Constant profile15375.66.2Parabolic profile421.223.425.9


## 5. Conclusion

In this paper, an accurate semianalytical modal theory for propagation in inhomogeneous waveguides has been developed for rectangular waveguides filled with graded-index metamaterials. Based on this theory, the cutoff frequencies and the dispersion relations have been computed for waveguides with various graded-index profile functions. These results have been validated by comparison with the results obtained by the commercial simulation software HFSS. It should be noted that HFSS cannot set the spatially dependent dielectric parameters, and the discretization and convergence analysis must be conducted in the simulation of graded-index metamaterials. The technique proposed in this paper can be employed for arbitrary inhomogeneity analysis. Moreover, the proposed technique is far more efficient than the HFSS simulation. Finally, exciting phenomena of such waveguides containing composite LHM and RHM with graded-index profiles have been highlighted. By choosing appropriate profile functions, the lower cutoff frequency might be substantially reduced compared with classical waveguides, and backward wave propagation below the cutoff frequency is obtained, which allows waveguide miniaturization. The losses and conductivity of the metamaterial are not considered in this paper. The modal theory presented in this paper is still verified if the losses and conductivity are known. Further work is related to the specific implementation of the graded-index metamaterials with consideration of dispersion and losses properties. It has been experimentally verified that a linear graded-index metamaterial could be realized by an array of split-ring resonators (SRRs) with a linearly varying substrate depth, and the graded permeability might be realized by arranging a set of wires that gradually changes their density. The technique may be easily applied to the circular waveguides, the elliptical waveguides, and the polygonal waveguides, by replacing the trigonometry functions in potential functions into the Bessel functions and Mathieu functions.



---

*Source: 6107378-2021-02-02.xml*