2.2: Spatial modes
- Page ID
- 113788
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Total internal reflection is not the only condition for guiding, however. Light is a wave that has physical extent in all dimensions. Thus, as demonstrated in Figure \(\PageIndex{1}\), light that has reflected from the top interface interferes with light that has reflected from the bottom interface. If the interference is not constructive, then the amplitude of the light wave will decrease with distance traveled, and eventually disappears entirely. More specifically, after two reflections, the phase change of the wave along the axis of the waveguide should be an integer multiple of \(\pi\). When this condition is satisfied, the light occupies a spatial mode of the waveguide.

The picture in Figure \(\PageIndex{1}\) is complicated by the fact that total internal reflection introduces an angle-dependent phase shift to the light wave on reflection. However, in what follows, we first apply a constant phase shift of \(\pi\) (thus treating the reflection as mirror-like), and then afterwards account for an angular dependent phase shift.
If you are not familiar with the wavefront picture above, it can be a bit confusing. The arrowed line represents the direction of travel for light (a light ray). The lines at right angle represent the peaks in amplitude, which are separated by one wavelength. This picture, then represents a frozen moment in time and allows us to compare the phase relationship between multiple waves, or, in this case, the wave with itself.
Mirror-like reflections
A geometric solution to the phase consistency problem is presented in Figure \(\PageIndex{2}\). The line \(\overline{AB}\) represents the path the ray would have taken if there was no reflection, while path \(\overline{AC}\) represents the path the light has actually taken. For self consistency the two paths can only differ by integer multiples of \(2\pi\). The accumulated phase over \(\overline{AB}\) is \(\phi_{AB} = \frac{2\pi\overline{AB}}{\lambda}\), where \(\lambda\) is the wavelength of the light. Likewise, the accumulated phase over \(\overline{AC}\) is \(\phi_{AC} = \frac{2\pi\overline{AC}}{\lambda}\). The difference between the two is \(\Delta\phi = \frac{2\pi\bar{AC}}{\lambda} - \frac{2\pi\overline{AB}}{\lambda} - 2\pi\), where the \(2\pi\) accounts for the additional phase shift due to two reflections. The difference must be \(\Delta\phi = 2m\pi\) where \(m\) is an integer and it doesn't really matter where we start \(m\), so we can subtract the \(2\pi\) and reset \(m\), leaving \(m = \frac{1}{\lambda}\left(\overline{AC}- \overline{AB}\right)\), where \(m = 1, 2, ...\).

\[\overline{AC} = \frac{d}{\sin\gamma} \label{eq:AC}\]
where \(\gamma\) is the complementary angle to \(\theta\). The triangle \(\overline{ABC}\) allows us to solve for the angles that meet the self-consistency condition.
\[\begin{align} \cos 2\gamma &= \frac{\overline{AB}}{\overline{AC}}\\ \cos 2\gamma &= 1-2\sin^2\gamma\\ 1-2\sin^2\gamma &= \frac{\overline{AB}}{\overline{AC}}\\ \overline{AC} - 2\overline{AC}\sin^2\gamma &= \overline{AB}\\ \overline{AC}-\overline{AB} &= 2\overline{AC}\sin^2\gamma\\ \overline{AC}-\overline{AB} &= \frac{2d\sin^2\gamma}{\sin\gamma}\\ m\lambda &= 2d\sin\gamma\\ \sin\gamma_m &= m\frac{\lambda}{2d} \label{eq:mirrormodes} \end{align}\]
Equation \ref{eq:mirrormodes} represents the set of angles that satisfy the self consistency relationship for any given wavelength. By setting \(\sin\gamma_m\) = 1 and solving for \(m\), we find that the maximum number of modes supported by the waveguide is
\[M \underset{\cdot}{=} \frac{2d}{\lambda}\label{eq:modeNumber}\]
where \(M\) is rounded down to the nearest integer. Note that this also means that, for this type of waveguide, there is a minimum dimension to enable a single mode to exist. If \(M=0\) that indicates that, for that wavelength, no propagating modes exist and no light will be transported by the waveguide. In this case, an attempt to couple light into the waveguide would result in a very strong scattering from the front face.
Later on, we will find that the number of modes is a severe limiting factor in the performance of an optical communications system. Indeed, for long distance communications, \(M=1\), otherwise known as single mode operation, is an important criteria.
Angle-dependent phase shift
Snell's law does not tell the full story for refraction and reflection. When light is reflected from an interface, it also undergoes a phase shift, as illustrated in Figure \(\PageIndex{3}\).

Note that the phase shift also depends on the orientation of the electric field with respect to the interface. In what follows, we only discuss the case for the electric field parallel to the interface, called transverse electric (TE) modes (when the magnetic field is parallel to the interface, this is called the transverse magnetic (TM) mode). The phase shift due to reflection is given by
\[\tan\frac{\phi_r}{2} = \sqrt{\frac{\sin^2\gamma_c}{\sin^2\gamma}-1} \label{eq:phaseShift}\]
where \(\gamma_c\) is the complementary angle to the critical angle, \(\theta_c\), and \(\phi_r\) is the phase shift due to reflection. Note that equation \ref{eq:phaseShift} is not valid for \(\gamma>\gamma_c\), where, instead, one must use the Fresnel amplitude reflection coefficients, that were used to in Figure \(\PageIndex{3}\).
The self-consistency expression given in equation \ref{eq:mirrormodes} becomes
\[\tan\left(\frac{\pi d}{\lambda}\sin\gamma - \frac{m\pi}{2}\right) = \sqrt{\frac{\sin^2\gamma_c}{\sin^2\gamma}-1} \label{eq:dielectricModes}\]
where \(m = 0, 1, 2, ...\). Equation \ref{eq:dielectricModes} is best solved graphically by looking for intercepts between the left hand side and right hand side, as shown in Figure \(\PageIndex{4}\). The total number of TE modes supported by a waveguide is given by
\[M \dot{=} \frac{2d\sin\gamma_c}{\lambda}\label{eq:dielectricModeNumber}\]
which is related to the NA by \(M \dot{=} \frac{2d}{\lambda_0}\textrm{NA}\), where \(M\) is now always rounded up, and \(\lambda_0\) denotes the wavelength in vacuum. In contrast to the mirror-like case, there is no cutoff for any wavelength. However, to achieve single mode operation, the dimensions have to be quite close to the operational wavelength.
In order to know precisely which mode is being discussed, a labeling scheme is used that tells us the polarization of the mode and its spatial character. For a slab waveguide, as discussed here, the modes are denoted as TE\(_0\), TE\(_1\), TE\(_2\), ... and TM\(_0\), TM\(_1\), TM\(_2\), ... for the transverse electric and transverse magnetic modes.

Going beyond slab waveguides
So far, we have looked at a waveguide that is an infinite wide slab so that we only have to deal with a single dimension. For square waveguides, the mathematics remains the same, but now there are modes associated with both spatial dimensions, thus, if the axis of propagation of a waveguide is along the \(z\) axis, then the two dimensions \(d_x\), and \(d_y\) each support \(M_x\) and \(M_y\) modes (calculated from equation \ref{eq:dielectricModeNumber}). The modes along each axis are independent of each other, so the total number of modes, \(M\) is all possible combinations \(M = M_xM_y\).
Optical fibers are cylinders, meaning that the geometric approach used above cannot be used. The electric and magnetic fields are never quite parallel to the interface of an optical fiber, resulting in so-called linearly polarized (LP) modes that are denoted by two integers \(l\) and \(m\). As with the dielectric slab waveguide, it is simplest to find the modes by graphically finding the intercepts between two equations (both Bessel's functions) that represent the electric field in the core and the cladding of the fiber.
Even the total number of modes supported by an optical fiber can be difficult to calculate. However, we are mostly interested in whether the fiber is single mode, multi-mode, and occasionally few-mode. To estimate this, the \(V\) parameter, which (sort of) represents the mismatch between light propagation speed in the core and cladding materials.
\[V=\frac{2\pi r_{co}}{\lambda_0}\textrm{NA} \label{eq:vParam}\]
where \(r_{co}\) is the radius of the core. Single mode operation is guaranteed for \(V<\)2.405. For multimode operation, the total number of modes can be estimated from \(M\approx\frac{4V^2}{\pi^2}\).