Fields in Dispersive Media

It is not just the case that matter affects the propagation of light — or more specifically electromagnetic (EM) radiation — it is also the case that light affects the matter through which it propagates. Conversely, this affects the propagation of light through the medium, but in a much more specific way; this effect is a function of the properties of both the material and the incident EM radiation. We will additionally discuss the effects of dispersion in confined (bounded) media, i.e., where the dispersion is a function of the arrangement of certain materials and unbounded media where EM radiation is free to propagate undisturbed. This will be important when we discuss the propagation electric field signals of such media as well.


Introduction
Incident electromagnetic (EM) radiation excites the molecules in a material, and these molecules become polarized; they respond according to the direction of the electric field. Figure 1 shows the initial step of the process of material polarization.
It should be adequately noted that there are several ways that material response to light-matter interaction can be analyzed. For the purposes of relevance, this chapter will analyze the mechanical response of the atom-as illustrated in Figure 2-in response to incident EM radiation waves.
In this chapter we will discuss polarization at the level of the interacting material, the effects of dispersion on a train of pulses, various types of polarization, and methods for reducing polarization to maintain the integrity of optical signals. But first, some words on non-dispersion in unbounded media, a more general term for what would normally be called "free space." Unbounded media are the baseline of understanding EM wave propagation. In unbounded media, waves are free to propagate unperturbed. Examples of unbounded media include the ocean, the air, and outer space. An unbounded medium would be the ideal location for an isotropic antenna as the radiation would be free to propagate in all directions, only weakening in accordance with the inverse square law. Such is ideal for radio towers which produce low-frequency EM waves which can propagate for kilometers and reach many people.
Additionally, unbounded media are generally non-dispersive. This is to say that the speed of energy propagation is orthogonal to the frequency of the said propagating energy. In the previous paragraph, it was said the waves in unbounded media  The mechanical response of the material to incoming monochromatic EM radiation. The vibrational response of the molecule is a function of its properties. This same principle is also related to the color of the flame that a molecule produces [1]. Since Belt B is traveling faster than Belt A, the distance between BB À1 and RB À1 will decrease. From the perspective of Belt A, Belt B is moving with a velocity of 0.5 cm/s. At time t ¼ 5, the center of BB À1 is aligned with that of RB À1 ; however, in that time the balls on Belt B have moved a net 15 cm, while those on Belt A have moved a net 12.5 cm. However, if one were to see the conveyor belts as a group, the belts would appear to have moved 2.5 cm (measuring from location of alignment) within that time, giving the conveyor belts a "group velocity" of 0.5 m/s even though the individual components are moving faster.
The situation described is analogous to the behavior of a wave propagating in free space. As the individual waves travel with their respective phase velocities along the guide, the two peaks become disaligned for a period [2]. Because the center of RB À1 will have traveled a distance λ RB þ Δz with a phase velocity v zRB and the center of BB À1 will have traveled a distance λ BB þ Δz with phase velocity v zBB , one can then write a set of equations as We can then write Δt as In the same way, one can express Δz as Finally, one can express the group velocity as We know that the wavevector k ¼ 2π=λ, so one express the group velocity as This is to say the group velocity is the degree change in the temporal frequency for every change in the spatial frequency. What this means is that the value of the frequency is on some level a function of the k-vector and ultimately a function of the wavelength. As illustrated in Figure 4, this has implications for the propagation of light through various media as different wavelengths will travel at different velocities and take various pathways therein [3].

Plasmonic dispersion
The Lorentz force is the force on a point charge due to EM fields. As discussed in the Introduction, EM radiation causes motion among the particles in a material as shown in Figure 2. The direction (vector) of this force is called the polarization. Assuming a sea of free particles in vacuum [4], the strength of the Lorentz force on a point charge is expressed as where , and B ! are the force, the charge, the electric field strength, the particle velocity, and the magnetic flux density, respectively. For our purposes, we consider an isotropic material where the electric permittivity is a simple scalar ϵ. This means that the Lorentz force vector will be unidirectional. In birefringent and anisotropic, the permittivity would be expressed as a tensor, and the Lorentz force vector would be multidirectional.
We additionally assume a dielectric unmagnetized material such that the value of B ! is assumed to be zero. This makes the response of the atoms into a simple harmonic oscillator. In this case we can express the Lorentz force as simply where m, ω, and d ! are the mass, the angular frequency, and dipole vector, respectively. We can simplify to get the dipole vector as where P ! is the polarization vector and N is the dipole density. The electric displacement field D ! is the permittivity ϵ multiplied by the electric field E ! , but in dielectric materials, the electric field induces a response in the material known as the polarization. This effect is added to the original displacement field to make it so that As a result, the effective permittivity is expressed in (11) as where the plasma frequency ω p ¼ q ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi N=ϵ 0 m p . The plasma frequency is the natural oscillation frequency of the displaced electrons in a neutral plasma [4,5] of free electrons where it is assumed that collisions are rare. Altogether we can express the dispersion relation thusly This means that the group velocity is expressed as In a material like germanium (Ge) the intrinsic carrier density N Ge is 2:5 Â 10 19 m À3 [6]; the rest of the variables are known constants. This gives germanium a plasma frequency of f p,Ge ¼ 44:8 GHz. We apply the relation in (13) to germanium to measure the frequency-dependent group velocity of light therein as demonstrated in Figure 5.
Below the plasma frequency, the value of the k-vector becomes purely imaginary, meaning that the electric field is purely evanescent-and non-propagatingat that point. The group velocity is the slope of the dispersion curve in Figure 6. For both curves, the lowest allowable frequency is the plasma frequency. Below this frequency, the permittivity is negative, the k-vector is imaginary, and, consequently, the field simply evanesces into the material and propagation completely ceases.

Dispersion in conductive media
In the previous section, we derived plasmonic dispersion from understanding the Lorentz force law as it pertains to electrons in a dielectric material. Dispersion can also occur in a conducting material where the charges therein are unbound. When we apply our knowledge of electromagnetics and electrostatics to a conductor, we discover a few things: 1. There is zero net charge within a conductor, although charges can still randomly move around therein as seen in Figure 7.
2. All the charge rests on the surface of the conductor, perpendicular to the surface.
3. The presence of an electric field induces the motion of electrons, creating a conduction current [7].
Additionally, in a conductor the free charges are always pushed to the surface so that electrostatic equilibrium is maintained. This means that the gradient of the electric field ∇ ! Á E ! equals zero. Things are different in the presence of an electric field, e.g., Figure 8, where the electrons are pulled in the opposite direction of the electric field.
With regard to conductive media, we can do the same thing through an understanding and application of Faraday's law which states that a time-varying magnetic flux density always accompanies a nonconservative electric field, i.e., a field E This is demonstrated by the equation We can combine this with Ampère's law where the curl of the magnetic field generates a time-varying electric field such that where J ! f is the free current density vector produced in the presence of an electric field as exemplified in Figure 8 and σ is the conductivity. We can calculate the effective permittivity this way such that . In this way we can write the dispersion relation.
This wavenumber k is also complex where k ¼ k 0 À jk 00 , where The real part of the wavenumber k 0 represents wavenumber inside the material and the imaginary part k 00 represents the attenuation coefficient inside the material.
Of course, there are various conditions that could simplify Eq. (16). The σ=ωϵ component in (15) is known as the dielectric loss tangent tan δ [8]. This figure δ describes a materials inherent ability to dissipate electromagnetic energy within it. For example, a material with a high conductivity will absorb-dissipate-EM radiation more quickly. Such a material has would have a high loss tangent compared to a good dielectric.
We can express the dielectric relaxation time τ d as ϵ=σ. In cases that δ ≪ 1 or τ À1 d ≪ ω-as would be the case in a lossy dielectric-we can approximate the permittivity as ϵ eff ≈ 1 À jσ 2ωϵ using the binomial expansion theorem. The material behaves like a standard dielectric as v g ¼ c=n, where n is the index of refraction, and k 00 ¼ ση=2, where η is the impedance of the material ffiffiffiffiffiffiffiffiffi ffi μ 0 =ϵ p . In cases where 1 ≪ δ or ω ≪ τ À1 d , the case of a good conductor, k ≈ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi σμ 0 ω=2 . In the latter case, the field penetrates only slightly into the material as O k 0 Additionally, the attenuation is a new phenomenon resulting from the fact that the k-vector is complex. When we take the reciprocal of this coefficient get another quantity known as the skin depth or penetration depth δ p , the extent to which the field penetrates the material. In a good conductor, δ p takes a value ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi 2=σμ 0 ω p . The values of ϵ and σ are listed in Table 1 for some good conductors [9]. Good conductors are characterized by an abundance of free electrons. Good conductors (good dielectrics), on the other hand, are characterized by a paucity of such electrons. Examples of good dielectrics are enumerated in Table 2 [10]. Figure 9 shows the interaction between an electric field and an good conductor (e.g., gold) at a wavelength of 1.55μm. The field decays very rapidly as gold has a skin depth of 5.66 nm and a dielectric relaxation time of 1.49 attoseconds (10 À18 s) at this wavelength. This contrasts with a lossy dielectric where the skin depth is much larger, and there is thus more field penetration as displayed in Figure 10 with germanium where the skin depth is roughly a centimeter and the dielectric relaxation time is 66.1 picoseconds.
Conductive media have more free electrons that move in the presence of an electric field than those that simply vibrate. Conductive media produce a loss (attenuation) component in propagation of incoming EM radiation. Additionally, the drift movement of electrons in conductive media causes heat to which the electrical energy from the incoming field is transferred. This energy conversion causes the dissipation of the field. Nevertheless, in the context of the phenomenon of dispersion, the degree of this attenuation is a function of the frequency of the incoming EM radiation.

Modal dispersion
So far, this chapter has covered the form that dispersion takes in media (i.e., material dispersion). However, there are other types of dispersion that exist which ultimately stem from the properties of the material. One of these other ways is  modal dispersion [11]. This is when the velocity of various modes changes as a function of the input wavelength. As displayed in Figure 11, the different modes that propagate at different group velocities and angles, which themselves are functions of an input wavelength. These modes that arise in a waveguide are the different forms that the electric field can take in a geometrically confined structure. These forms are functions of satisfying the Helmholtz equation for the geometry of the confining structure. Only a finite number of fields do this and produce the appropriate wavevector in the z-direction, β. Consequently, only a finite number of modes will be allowed in the structure.
The angle at which certain modes propagate within makes a difference as to the speed of propagation because shallower angles have less deviation from the center equilibrium as shown in Figure 12.
In a rectangular waveguide, the total distance a mode must travel is expressed as where θ x and θ y are the angles of deviation in the x-and y-directions, respectively. The higher the order of the mode, the steeper the angle of the travel, and the more distance the mode must traverse through the waveguide. In Figure 13, the same principle holds true in an optical fiber. The rays of the fundamental modes (TM and TE) take a meridional form (a) passing through origin as they propagate and the hybrid modes-HE and EH-spiral in a helical form taking a twisted path through the fiber.
Additionally, modes that have a shallower angle of interaction with the corecladding boundary will have higher group velocities. A higher-order mode has a  lower β value as its k-vectors in other directions are higher, leading to steeper interaction angles with the core-cladding boundary.
The way that the modes propagate through the fiber at different velocities makes the waveguide an effective multipath propagation channel where the received signal is a sum of the various scaled "echoes" of the original signal carried by the various modes that propagate therein.
As a result, this dispersion distorts the pulse and corrupts the information. As shown in Figures 14 and 15, the pulses deform from their respective Gaussian and rectangular forms as a result of the different path delays, a phenomenon known as delay distortion [13]. Calculation of modal delays is used in order to set various bit rates for optical pulses.
One way to ameliorate modal (intermodal) dispersion is to use a single-mode waveguide where the only one mode is carried through. Another way is to make use of waveguides with a parabolic index of refraction, producing a quasi-sinusoidal pathway through the waveguide as demonstrated in Figure 16. Rays taking meridional (a) and helical (b) forms in an optical fiber [12]. The ray switches paths as it passes the region of highest index, meaning that it spends less of its time in the region of lowest index. This is to say that the phase velocity of the wave changes within the waveguide. From a ray optics perspective, the crest of the wave has a position-varying velocity, effectively making continuous switches through "different materials." The dispersion characteristic for graded-index (GRIN) waveguides is quite different as well. Figure 17 shows that higher wavelengths are correlated with less dispersion for GRIN waveguides. The step-index, on the other hand, maintains a positive wavelength dispersion characteristic.
It should be noted that modal dispersion is more of a function of the arrangement of media, not the sole medium itself. Changes in size, cladding material, and wavelength can affect the propagation of light therein.

Chromatic dispersion
Changes in wavelength are important in understanding how light will propagate. As shown in Figure 18, the way these β values vary with the wavelength produces various group velocities for different modes. The higher-order modes have a lower group velocity as they travel. The index of the fused silica core in the waveguide  The pathway through the waveguide is sinusoidal [14] as it is bound to its axis. described in Figure 18 is 1.5, meaning that light will travel roughly 2 Â 10 8 m/s therein. However, the contribution from dispersion increases the effective index n eff "experienced" by each mode at a particular wavelength.  The plot in Figure 18 can be reconfigured to show that each mode has an effective index of refraction in the core that varies with its wavelength. This means that different "colors" of light will have different refractive indices in the waveguide. Figure 19 shows that as λ increases, the group velocity decreases. Figure 19 shows that there is an added index of refraction that comes from the dispersion that contributes the delay of the modes; this is also known as the group index N g .

Intramodal dispersion
Dispersion can occur even in the context of a single-mode waveguide (of whatever geometry). In the previous section, we demonstrated a difference in n eff as it pertains to various wavelengths stimulating a particular mode. However, we can observe a similar phenomenon in a single-mode fiber.
In the Introduction, we discussed how one can model the plasmonic impulse response of a material in the presence of an electric field as a mechanical system for which where E 0 , ω 0 , e, m, and γ are the electric field amplitude, the resonant frequency, the electric charge, the mass of the electron, and the damping constant, respectively [11]. We can multiply the quantity described in (21) be e to get the dipole moment vector. We additionally multiply by the dipole density N to get the bulk polarization. We see in (9) that the displacement field is equal to the added where Δν N is the FWHM frequency point or the natural linewidth [17]. This is the spectral profile of a typical laser centered around a frequency ν ¼ ν 0 . We can convert this line-shape function to one of the k-vectork where where k 0 is the center spatial frequency. The spectrum described in (28) shows that there is an uncountable infinity-i.e., a continuum-of k-vectors (and thus wavelengths) that is allowed.
In the case of a gas, the same can be done with the Boltzmann relation that is used to compute the statistical distribution of velocities (and thus k-vectors) under Doppler [18] broadening where the velocity component v is normally distributed N 0, k B T=m ð Þwhile the energy E undergoes a chi-squared distribution χ 2 1 E ð Þ as shown in (28): where m, k B , T, and v are the particle's mass, Boltzmann's constant, the temperature, and the velocity, respectively.
Because of the Doppler effect, the radiation will be shifted to a higher frequency ν which can be expressed as where ν 0 is the center frequency. With (30) in mind, we can rewrite (28) in terms of ν and convert to the k-domain like so What this means is that in the case of a realistic laser, dispersion within a confined medium will result even if the medium allows for only one mode to be activated. This is to say that Lorentzian homogenous broadening and Gaussian Doppler broadening both result in intramodal pulse broadening.
What is interesting is the counterintuitive connection between the frequency broadening and the temporal pulse broadening as their behaviors mirror each other. However, this is because we are dealing with two separate systems that happen to be connected. The first system is that of the laser linewidth and the decaying atoms and spontaneous thermal processes that diminish monochromaticity. These processes spread the spectrum out and a broader range of frequency (and wavelength) contribution to the spectrum. The second system is that of the delays accrued within the fiber-waveguide channel according to the wavelength contributions of the input spectrum. The input laser spectrum allows for a larger contribution to delays that will be "experienced" by the larger range of frequency and wavelength components, which broaden the temporal profile of the signal.
As noise is intrinsic to a lasing system so is dispersion when that system is directed into a medium; this is especially true in the context of confined or bounded media and their applications in modern communication networks.

Conclusion
Dispersion has been discussed in this chapter as it relates to electric fields in various media. We first begin with an understanding of electron responses to EM radiation; i.e., the light-matter interaction as understood in classical electromagnetics. When then used this concept to understand how the velocity of a signal through a medium is a function of the frequency of that signal.
Additionally, we explained and developed the concept of a group velocity with a simple analogy and then moved on to discuss that concept at length in various media. We contrasted group velocity with phase velocity and developed dispersion relations which contrasted the two.
We took the discussion of dispersive materials to light-confining structures. Optical waveguides and fibers function as dispersive materials not only because of their material properties but also because these materials are arranged. The arrangement of these materials-i.e., core-cladding indices, size, index distribution, etc.-all plays a role in the phenomenon of dispersion. This also influences the allowable forms of light that propagate in one of these confined media; the properties of the medium affect these modes of propagation. We additionally learned about ways in which this effect can be eliminated for engineering applications.
We finally ended this chapter with a discussion of effect of intramodally dispersive media on the frequency spectrum and time evolution of electric field signals. This was demonstrated by first dismantling the "myth of monochromaticity" and understanding light itself as a spectrum. This helped us better understand how a single mode can undergo dispersion with itself. It most importantly demonstrated the degree to which dispersion in inherent in EM propagation through all media.