7.6: Waveguides- H and E Waves
( \newcommand{\kernel}{\mathrm{null}\,}\)
Let us now return to Eqs. (100) and explore the H- and E-waves – with, respectively, either Hz or Ez different from zero. At the first sight, they may seem more complex. However, Eqs. (101), which determine the distribution of these longitudinal components over the cross-section, are just the 2D Helmholtz equations for scalar functions. For simple cross-section geometries, they may be readily solved using the methods discussed for the Laplace equation in Chapter 2, in particular the variable separation. After the solution of such an equation has been found, the transverse components of the fields may be calculated by differentiation, using the simple formulas,
Et=ik2t[kz∇tEz−kZ(nz×∇tHz)],Ht=ik2t[kz∇tHz+kZ(nz×∇tEz)],
which follow from the two equations in the first line of Eqs. (100).52
In comparison with the boundary problems of electro- and magnetostatics, the only conceptually new feature of Eqs. (101) is that they form the so-called eigenproblems, with typically many solutions (eigenfunctions), each describing a specific wave mode, and corresponding to a specific eigenvalue of the parameter kt. The good news here is that these values of kt are determined by this 2D boundary problem and hence do not depend on kz. As a result, the dispersion law ω(kz) of any mode, which follows from the last form of Eq. (102),
Universal dispersion relationω=(k2z+k2tεμ)1/2≡(ν2k2z+ω2c)1/2,
is functionally similar for all modes. It is also similar to that of plane waves in plasma (see Eq. (38), Fig. 6, and their discussion in Sec. 2), with the only differences that the speed in light c is generally replaced with ν=1/(εμ)1/2, i.e. the speed of the plane (or any TEM) waves in the medium filling the waveguide, and that ωp is replaced with the so-called cutoff frequency
ωc≡νkt,
specific for each mode. (As Eq. (101) implies, and as we will see from several examples below, kt has the order of 1/a, where a is the characteristic dimension of waveguide’s cross-section, so that the critical value of the free-space wavelength λ≡2πc/ω is of the order of a.) Below the cutoff frequency of each particular mode, such wave cannot propagate in the waveguide.53 As a result, the modes with the lowest values of ωc present special practical interest, because the choice of the signal frequency ω between the two lowest values of the cutoff frequency (123) guarantees that the waves propagate in the form of only one mode, with the lowest kt. Such a choice enables engineers to simplify the excitation of the desired mode by wave generators, and to avoid the parasitic transfer of electromagnetic wave energy to undesirable modes by (virtually unavoidable) small inhomogeneities of the system.
The boundary conditions for the Helmholtz equations (101) depend on the propagating wave type. For the E-modes, with Hz=0 but Ez≠0, the condition Eτ=0 immediately gives
Ez|C=0,
where C is the inner contour limiting the conducting wall’s cross-section. For the H-modes, with Ez=0 but Hz≠0, the boundary condition is slightly less obvious and may be obtained using, for example, the second equation of the system (100), vector-multiplied by nz. Indeed, for the component normal to the conductor surface, the result of such multiplication is
ikz(Ht)n−ikZ(nz×Et)n=∂Hz∂n.
But the first term on the left-hand side of this relation must be zero on the wall surface, because of the second of Eqs. (104), while according to the first of Eqs. (104), the vector Et in the second term cannot have a component tangential to the wall. As a result, the vector product in that term cannot have a normal component, so that the term should equal zero as well, and Eq. (125) is reduced to
∂Hz∂n|C=0.
Let us see how does all this machinery work for a simple but practically important case of a metallic-wall waveguide with a rectangular cross-section – see Fig. 22

In the natural Cartesian coordinates, shown in this figure, both Eqs. (101) take the simple form
(∂2∂x2+∂2∂y2+k2t)f=0, where f={Ez, for E− modes, Hz, for H− modes.
From Chapter 2 we know that the most effective way of solution of such equations in a rectangular region is the variable separation, in which the general solution is represented as a sum of partial solutions of the type
f=X(x)Y(y).
Plugging this expression into Eq. (127), and dividing each term by XY, we get the equation,
1Xd2Xdx2+1Yd2Ydy2+k2t=0,
which should be satisfied for all values of x and y within the waveguide’s interior. This is only possible if each term of the sum equals a constant. Taking the X-term and Y-term constants in the form (−k2x) and (−k2y), respectfully, and solving the corresponding ordinary differential equations,54 for the eigenfunction (128) we get
f=(cxcoskxx+sxsinkxx)(cycoskyy+sysinkyy), with k2x+k2y=k2t,
where the constants c and s should be found from the boundary conditions. Here the difference between the H-modes and E-modes comes in.
For the H-modes, Eq. (130) is valid for Hz, and we should use the boundary condition (126) on all metallic walls of the waveguide, i.e. at x=0 and a; and y=0 and b – see Fig. 22. As a result, we get very simple expressions for eigenfunctions and eigenvalues:
(Hz)nm=Hlcosπnxacosπmyb,
kx=πna,ky=πmb, so that (kt)nm=(k2x+k2y)1/2=π[(na)2+(mb)2]1/2,
where Hl is the longitudinal field’s amplitude, and n and m are two integer numbers – arbitrary besides that they cannot be equal to zero simultaneously.55 Assuming, just for certainty, that a≥b (as shown in Fig. 22), we see that the lowest eigenvalue of kt, and hence the lowest cutoff frequency (123), is achieved for the so-called H10 mode with n=1 and m=0, and hence with
Fundamental mode’s cutoff(kt)10=πa,
thus confirming our prior estimate of kt.
Depending on the a/b ratio, the second-lowest kt (and hence ωc) belongs to either the H11 mode with n=1 and m=1:
(kt)11=π(1a2+1b2)1/2≡[1+(ab)2]1/2(kt)10,
or to the H20 mode with n=2 and m=0:
(kt)20=2πa≡2(kt)10.
These values become equal at a/b=√3≈1.7; in practical waveguides, the a/b ratio is not too far from this value. For example, in the standard X-band (∼10−GHz) waveguide WR90, a≈2.3 cm (fc≡ωc/2π≈6.5 GHz), and b≈1.0 cm.
Now let us have a look at the alternative E-modes. For them, we still should use the general solution (130) with f=Ez, but now with the boundary condition (124). This gives us the eigenfunctions
(Ez)nm=Elsinπnxasinπmyb,
and the same eigenvalue spectrum (132) as for the H modes. However, now neither n nor m can be equal to zero; otherwise Eq. (136) would give the trivial solution Ez(x,y)=0. Hence the lowest cutoff frequency of TM waves is achieved at the so-called E11 mode with n=1, m=1, and with the eigenvalue given by Eq. (134), always higher than (kt)10.
Thus the fundamental H10 mode is certainly the most important wave in rectangular waveguides; let us have a better look at its field distribution. Plugging the corresponding solution (131) with n=1 and m=0 into the general relation (121), we easily get
(Hx)10=−ikzaπHlsinπxa,(Hy)10=0,
(Ex)10=0,(Ey)10=ikaπZHlsinπxa.
This field distribution is (schematically) shown in Fig. 22. Neither of the fields depends on the coordinate y – the feature very convenient, in particular, for microwave experiments with small samples. The electric field has only one (in Fig. 22, vertical) component that vanishes at the side walls and reaches its maximum at the waveguide’s center; its field lines are straight, starting and ending on wall surface charges (whose distribution propagates along the waveguide together with the wave). In contrast, the magnetic field has two non-zero components ( Hx and Hz), and its field lines are shaped as horizontal loops wrapped around the electric field maxima.
An important question is whether the H10 wave may be usefully characterized by a unique impedance introduced similarly to ZW of the TEM modes – see Eq. (115). The answer is not, because the main value of ZW is a convenient description of the impedance matching of a transmission line with a lumped load – see Fig. 19 and Eq. (118). As was discussed above, such a simple description is possible (i.e., does not depend on the exact geometry of the connection) only if both dimensions of the line’s cross-section are much less than λ. But for the H10 wave (and more generally, any non-TEM mode) this is impossible – see, e.g., Eq. (129): its lowest frequency corresponds to the TEM wavelength λmax=2π/(kt)min=2π/(kt)10=2a. (The reader is challenged to find a simple interpretation of this equality.)
Now let us consider metallic-wall waveguides with a round cross-section (Fig. 23a). In this single-connected geometry, the TEM waves are impossible again, while for the analysis of H-modes and E-modes, the polar coordinates {ρ,φ} are most natural. In these coordinates, the 2D Helmholtz equation (101) takes the following form:
[1ρ∂∂ρ(ρ∂∂ρ)+1ρ2∂2∂φ2+k2t]f=0, where f={Ez, for E - modes, Hz, for H - modes.
Separating the variables as f=R(ρ)F(φ), we get
1ρRddρ(ρdRdρ)+1ρ2Fd2Fdφ2+k2t=0.
But this is exactly the Eq. (2.127) that was studied in Sec. 2.7 in the context of electrostatics, just with a replacement of notation: γ→kt. So we already know that to have 2π-periodic functions F(φ) and finite values R(0) (which are evidently necessary for our current case – see Fig. 23a), the general solution must have the form given by Eq. (2.136), i.e. the eigenfunctions are expressed via integer-order Bessel functions of the first kind:
fnm=Jn(knmρ)(cncosnφ+snsinnφ)≡const×Jn(knmρ)cosn(φ−φ0),
with the eigenvalues knm of the transverse wave number kt to be determined from appropriate boundary conditions, and an arbitrary constant φ0.

As for the rectangular waveguide, let us start from the H-modes (f=Hz). Then the boundary condition on the wall surface (ρ=R) is given by Eq. (126), which, for the solution (141), takes the form
ddξJn(ξ)=0, where ξ≡kR.
This means that eigenvalues of Eq. (139) are
kt=knm=ξ′nmR,
where ξ′nm is the mth zero of the function dJn(ξ)/dξ. Approximate values of these zeros for several lowest n and m may be read out from Fig. 2.18; their more accurate values are given in Table 1 below.
m=1 | 2 | 3 | |
n=0 | 3.83171 | 7.015587 | 10.1735 |
1 | 1.84118 | 5.33144 | 8.53632 |
2 | 3.05424 | 6.70613 | 9.96947 |
3 | 4.20119 | 8.01524 | 11.34592 |
The table shows, in particular, that the lowest of the zeros is ξ′11≈1.84.56 Thus, perhaps a bit counter-intuitively, the fundamental mode, providing the lowest cutoff frequency ωc=νknm, is H11, corresponding to n=1 rather than n=0:
Hz=HlJ1(ξ′11ρR)cos(φ−φ0).
It has the transverse wave number is kt=k11=ξ′11/R≈1.84/R, and hence the cutoff frequency corresponding to the TEM wavelength λmax=2π/k11≈3.41R. Thus the ratio of λmax to the waveguide’s diameter 2R is about 1.7, i.e. is close to the ratio λmax/a=2 for the rectangular waveguide. The origin of this proximity is clear from Fig. 24, which shows the transverse field distribution in the H11 mode. (It may be readily calculated from Eqs. (121) with Ez=0, and Hz given by Eq. (144).)

One can see that the field structure is actually very similar to that of the fundamental mode in the rectangular waveguide, shown in Fig. 22, despite the different nomenclature (which is due to the different coordinates used for the solution). However, note the arbitrary constant angle φ0, indicating that in circular waveguides the transverse field’s polarization is arbitrary. For some practical applications, such degeneracy of these “quasi-linearly-polarized” waves creates problems; some of them may be avoided by using waves with circular polarization.
As Table 1 shows, the next lowest H-mode is H21, for which kt=k21=ξ′21/R≈3.05/R, almost twice larger than that of the fundamental mode, and only then comes the first mode with no angular dependence of any field, H01, with kt=k01=ξ′01/R≈3.83/R,57 followed by several angle-dependent modes: H31, H12, etc.
For the E modes, we may still use Eq. (141) (with f=Ez), but with the boundary condition (124) at ρ=R. This gives the following equation for the problem eigenvalues:
Jn(knmR)=0, i.e. knm=ξnmR,
where ξnm is the mth zero of function Jn(ξ) – see Table 2.1. The table shows that the lowest kt is equal to ξ01/R≈2.405/R. Hence the corresponding mode (E01), with no angular dependence of its fields, e.g.
Ez=ElJ0(ξ01ρR),
has the second-lowest cutoff frequency, ~30% higher than that of the fundamental mode H11.
Finally, let us discuss one more topic of general importance – the number N of electromagnetic modes that may propagate in a waveguide within a certain range of relatively large frequencies ω>>ωc. It is easy to calculate for a rectangular waveguide, with its simple expressions (132) for the eigenvalues of {kx,ky}. Indeed, these expressions describe a rectangular mesh on the [kx,ky] plane, so that each point corresponds to the plane area ΔAk=(π/a)(π/b), and the number of modes in a large k-plane area Ak>>ΔAk is N=Ak/ΔAk=abAk/π2=AAk/π2, where A is the waveguide’s cross-section area.58 However, it is frequently more convenient to discuss transverse wave vectors kt of arbitrary direction, i.e. with an arbitrary sign of their components kx and ky. Taking into account that the opposite values of each component actually give the same wave, the actual number of different modes of each type ( E- or H-) is a factor of 22≡4 lower than was calculated above. This means that the number of modes of both types is
N=2AkA(2π)2.
Let me leave it for the reader to find hand-waving (but convincing :-) arguments that this mode counting rule is valid for waveguides with cross-sections of any shape, and any boundary conditions on the walls, provided that N>>1.
Reference
52 For the derivation of Eqs. (121), one of two linear equations (100) should be first vector-multiplied by nz. Note that this approach could not be used to analyze the TEM waves, because for them kt=0,Ez=0,Hz=0, and Eqs. (121) yield uncertainty.
53 An interesting recent twist in the ideas of electromagnetic metamaterials (mentioned in Sec. 5 above) is the so-called ε-near-zero materials, designed to have the effective product εμ much lower than ε0μ0 within certain frequency ranges. Since at these frequencies the speed ν (4) becomes much lower than c, the cutoff frequency (123) virtually vanishes. As a result, the waves may “tunnel” through very narrow sections of metallic waveguides filled with such materials – see, e.g., M. Silveirinha and N. Engheta, Phys. Rev. Lett. 97, 157403 (2006).
54 Let me hope that the solution of equations of the type d2X/dx2+k2xX=0 does not present any problem for the reader, at least due to their prior experience with problems such as standing waves on a guitar string, wavefunctions in a flat 1D quantum well, or (with the replacement x→t) a classical harmonic oscillator.
55 Otherwise, the function Hz(x,y) would be constant, so that, according to Eq. (121), the transverse components of the electric and magnetic field would equal zero. As a result, as the last two lines of Eqs. (100) show, the whole field would be zero for any kz≠0.
56 Mathematically, the lowest root of Eq. (142) with n=0 equals 0. However, it would yield k=0 and hence a constant field Hz, which, according to the first of Eqs. (121), would give a vanishing electric field.
57 The electric field lines in the H01 mode (as well as all higher H0m modes) are directed straight from the symmetry axis to the walls, reminding those of the TEM waves in the coaxial cable. Due to this property, these modes provide, at ω>>ωc, much lower energy losses (see Sec. 9 below) than the fundamental H11 mode, and are sometimes used in practice, despite the inconvenience of working in the multimode frequency range.
58 This formula ignores the fact that, according to the above analysis, some modes (with n=0 and m=0 for the H modes, and n=0 or m=0 for the E modes) are forbidden. However, for N>>1, the associated corrections of Eq. (147) are negligible.