## 1. Introduction

Up to now, the photon is understood as the quantum of electromagnetic radiation. In 1905, Einstein proposed that energy quantization was a property of electromagnetic radiation itself. Accepting the validity of Maxwell's theory, he pointed out that several experiments with results unpredictable by the classical electrodynamics theory could be explained if the energy of a light wave was localized into point‐like quanta moving independently of one another. A very simple and intuitive interpretation at the level of undergraduate teaching of quantum physics is that photons are the fundamental particles of light having the property that they behave both as a particle and a wave (wave‐particle duality). They also have characteristics, which make them different from other particles. One of these characteristics is that, as theorized up to now, when freely propagating, they behave as massless particles not interacting between them and carrying linear and intrinsic angular momentum.

In modern terms, a photon is considered as an elementary excitation of the quantized electromagnetic field, and it can be treated as a (quasi‐) particle, roughly analogous to an electron. It has unique properties, arising from its zero rest mass and its spin‐one nature. In particular, since the early days of quantum mechanics, it has been argued that there is no position operator for a photon, leading someone to conclude that there can be no properly defined wave function, in the Schrödinger sense, which allows to know the probability of finding the particle in a given spatial region. Nevertheless, photon position operators have been postulated whose eigenvectors form bases of localized states, as in Ref. [1].

The aim of this chapter is study the propagation of photons through dispersive media. This chapter is organized as follows. In Section 2, a semiclassical description of the dynamics of free photons is presented using a Dirac‐like equation. In Section 3, the positive and negative energy solutions arising from these equations are analysed. A comparison with the case of a free electron is made. Section 4 presents an analysis of the interaction of photons with the medium considering both a classical as a quantum treatment of light. Section 5 includes an analysis of the propagation of photons along a waveguide where they behave as if they did have a finite mass. In Section 6, some technological applications of the theoretical frame here presented are shown, such as the use of the properties of metamaterials to control the propagation of waves through waveguides filled with this kind of materials.

## 2. A Dirac‐like equation for the photon

Maxwell's equations can be considered as a classical field theory for a single photon that can be field (or “second”) quantized to obtain a quantum field theory of many photons.

In Ref. [2], it has been shown that in a region without sources Maxwell's equations can be written in the form of a Schrödinger‐like equation for a single photon adding a transversality condition. Although in quantum mechanics Schrödinger's equation is valid for describing the dynamics of a nonrelativistic particle, its application for the case of a photon must be considered only within the context of classical electrodynamics and taking into account that we are dealing with an equation which has the form of Schrödinger's equation and that it is equivalent to Maxwell's equations.

Considering that one important requirement of a quantum theory for describing the dynamics of photons is Lorentz invariance, in this section, we study the application of a Dirac‐like equation. In Refs. [3, 4] it has been shown that Maxwell's equations without sources can be written in a form analogous to that of Dirac's equation for a free electron. These last works also show that optical spin and light orbital angular momentum can be obtained from this Dirac‐like equation. As an extension of these works we give arguments for obtaining this equation in a similar form to those used for the deduction of Dirac's equation for an electron, starting from the relativistic expression for the energy. For example, see Ref. [5]. We also study the positive and negative energy states obtained from the corresponding Hamiltonian and the form that this equation takes for the propagation of a photon in a magnetodielectric medium.

As in Ref. [6], we begin with a derivation of a Dirac‐like equation for a photon starting from Dirac equation for a massless particle in free motion, so that we postulate an equation of the form:

Since Eq. (1) is linear in the time derivative, it seems natural to construct a Hamiltonian operator also linear in the spatial derivatives. This is compatible with the energy‐momentum relation for the photon

where

As in Ref. [7], a possible election for

where

In Section 3, it is shown that the above equations lead to the relation

There exist a variety of Dirac‐like formulations of Maxwell's equations and alternative ways for choosing the wave function ψ. Considering that photons have only energy and no other scalar quantities such as mass or charge, it is convenient to choose ψ so that its modulus squared correspond to energy density not of probability density for localization as is the case of a particle with mass like the electron. Therefore, we choose as wave function the following column vector of dimension

In this last equation, the components of

This last expression corresponds to the density of energy in the electromagnetic field.

For an electromagnetic wave propagating in a linear magnetodielectric and nonconducting medium, there is an induced polarization and magnetization classically represented by the polarization and magnetization vectors

where

where

From Eq. (8), a continuity equation for the wave function

Using Eqs. (6)–(9) and the definition of the column matrix

## 3. Positive and negative energy states for a photon

The energy eigenvalues are obtained by looking for stationary state solutions of the Dirac‐like equation. A plane wave solution of this equation has the form

for

The corresponding solutions for the energy are *z* direction with momentum

It is important to consider that Dirac's equation for a particle like an electron only gives positive and negative energy solutions.

## 4. Interaction between photons and the medium

### 4.1. Introduction

When light passes through a material, there is an electromagnetic interaction with the particles of the medium. This interaction is macroscopically manifested by two main effects: absorption of energy from the incident beam and scattering.

Considering that every particle has electric charge that acquires a motion due to the electric field associated to the incident electromagnetic wave, the absorption of energy may be understood using a phenomenological model of electric dipoles with negative charges whose positions oscillate with respect to the centre of positive charges, with a frequency corresponding to that of the incident light. This oscillatory motion has a damping associated to the dielectric losses.

The scattering of light may be thought of as the redirection that takes place when an electromagnetic wave encounters an obstacle or non‐homogeneity. The accelerated motion of the charges gives rise to radiation of electromagnetic energy in all directions producing secondary waves, process known as scattering.

### 4.2. Classical model for the interaction between light and matter

A classical model for representing the optical response of a polarizable medium through which travels a monochromatic electromagnetic wave of frequency *ω* is the Drude model where it is imagined that due to the electric field vector associated to this wave, each electron bound to the nucleus of an atom performs harmonic oscillations. Therefore, oscillating electric dipoles are formed, and considering that the electric dipolar moment is related to the electric field through the atomic polarizability

In this last equation,

At a macroscopic scale, the formation of electric dipoles in a dielectric material subjected to an applied electric field is described by means of the polarization vector

In this last equation,

If *Ze* electrons, the polarization vector is rewritten as

Therefore, we obtain the Drude model for the dielectric function:

where

From Eq. (16), it can be seen that a lossy dielectric medium has a complex refraction index

For the magnetic permeability, as in Ref. [10], a Drude‐Lorentz model similar to that given by Eq. (16) can be used:

where

When both

### 4.3. Quantum treatment

When an incident pulse enters into a dielectric medium, it undergoes modifications due to dispersion and absorption and in the case of a dielectric slab due to reflections from its surfaces. These modifications give rise to a distortion of the transmitted pulse in comparison with the incident pulse. Furthermore, the transmission of the pulse may be affected by thermal emission from the slab at elevated temperatures.

In classical electrodynamics, the interaction of light with matter is performed in two stages. First, an explicit model of the medium is assumed and its response to an electromagnetic field is calculated. The interaction is represented by the dielectric function, which embodies the optical properties of the material. In the second stage, this dielectric function is used for studying the propagation of the electromagnetic wave through the medium, determining effects as energy absorption and velocity of propagation. Nevertheless, for a finite number of photons, there are effects that cannot be described by a classical approach such as zero average electric field between two conducting plates and electric force between the plates even if the number of photons is zero (the so‐called vacuum fluctuation and Casimir forces, respectively). In addition, for a nonclassical pulse propagating in an absorbing and dispersive medium, there are modifications in the correlation properties that can only be described by a quantum theory of the photon.

As in Refs. [12, 13], the formalism for electromagnetic field quantization in a dispersive and absorbing dielectric, in general, includes the following steps:

Express Maxwell’s equations in terms of transverse electric and magnetic vector operators obtained from a vector potential operator.

Express the above vector potential operator as a function of the complex refraction index and of a current operator associated with noise sources coupled with the electromagnetic field in presence of lossy dielectrics.

Incorporate boson‐type operators and commutation relations between the electromagnetic field operators.

In what follows, we illustrate the application of this procedure for the case of light propagating in the

As shown in Ref. [12] for a state with *N* photons, the quantum field‐field correlations between two points placed over the *x* axis at positions

We consider the following definition of the spatial first‐order coherence function for two points placed over the *x* axis and separated by a distance *s* (e.g., see Refs. [14, 15])

Therefore, from Eqs. (18–20), we obtain

As an example, we make a comparison of this last result with that obtained calculating the classical coherence function. For that purpose, we consider a beam of light produced by excitation of two linearly polarized waves with frequencies

The Fourier transform of this field is

From the definition of the classical spatial coherence function:

we obtain that the classical spatial first order coherence function calculated for this case is

Due to the exponential factor of this last equation, the classical coherence function for the considered case has a value lower than 1 meaning that measurements of the electric field at two separated points are partially correlated in a medium with absorption and for

Comparing Eqs. (21) and (25), it can be seen that if absorption is neglected the classical model predicts total coherence, while the quantum treatment in this case gives an oscillatory behaviour of the spatial coherence function with respect to the distance between the points considered, with null partial coherence for some values of this distance. It is worth to note that this result is also valid for a left‐handed medium if absorption may be neglected.

## 5. Propagation of photons through a waveguide

When the propagation of an electromagnetic wave of a given frequency

This is illustrated considering the propagation of a transverse electric (TE) mode in a rectangular waveguide with transversal section having dimensions

where

For a relativistic fermion of mass *m*, Klein‐Gordon's equation can be obtained from Dirac's equation. For example, see Ref. [18]. This makes natural to wonder what is the form that the Dirac‐like equation takes for a photon moving along a waveguide. This can be determined writing the Dirac‐like equation in the following form:

with

For a photon moving along the waveguide, let us write the wave functions as

From Eqs. (27) and (29), we get

This last equation may be recast as

For propagation along the *z* direction, as it is considered in this case:

Equation (31) is the Dirac‐like equation for photons moving along the waveguide and it can be shown that leads to the Klein‐Gordon equation applying the operator

## 6. Some technological applications

On the last years, several articles about the properties of a special kind of materials, known as metamaterials, have been published. These materials can exhibit negative values on their permittivity or permeability. They are also named “left handed materials” and can have negative refraction index, which leads to interesting phenomena for the wave propagation. The effect of negative refraction was predicted in 1968 by Veselago in Ref. [19], principle that has led to many technological applications. For example, see Ref. [20]. Nowadays is possible to build artificial metamaterials with different geometries. These arrays can achieve negative values of permittivity or permeability, achieving either single negative material (SNG), where *ϵ* or *μ* are negative, or double negative material (DNG), where both *ϵ* and *μ* are negative. On the other hand, natural medias, such as plasmas, can behave as an SNG media, depending on its physical characteristics.

Metamaterials constructed by circuit arrays are based on a group of elements organized periodically, and designed in order to respond to an impinging electromagnetic field. The size and spacing of each element of the array must be much lower than the wavelength of the wave interacting with the array. This will allows that the impinging wave interacts with the artificial material as a homogeneous material with certain

There are different type of structures to obtain negative permittivity and permeability. To obtain a negative permittivity, periodic structures based on wire arrays that are based on the Drude‐Lorentz model for dielectric constant are used. For example, see Refs. [21–23]. On the other hand, to obtain negative permeability values, split ring resonators (SRR) and the induced current on wire structures are used, as can be seen in Refs. [21, 24, 25]. Currently, we can find two kind of artificial metamaterials that can exhibit negative refraction: photonic crystals, as shown in Refs. [26, 27], and composite materials as shown in Ref. [28]. Composite materials exhibit simultaneously negative permittivity and permeability within a certain frequency range. This immediately leads to a negative index of refraction. Dielectric photonic crystals are composed of materials with positive

The possibility of having SNG or DNG metamaterials opens a huge number of new applications that these physical characteristics can offer. Some examples of these applications are invisible materials or cloaking, as shown in Refs. [30, 31], phase control of propagating modes on waveguides, as shown in Refs. [32, 33], antenna miniaturization as shown in Refs. [34, 35], and superlens, as shown in Refs. [36, 37].

For each application, the design parameters of the metamaterials are an important issue. If the parameters that determine the permittivity and permeability of the material are known beforehand, it is possible to predict the behaviour of the electromagnetic wave that propagates on the media. The following study is intended to analyse the possible variations of the wavenumber of metamaterials, depending on their design parameters.

### 6.1. Wavenumber on SNG media depending on metamaterial parameters

The first analysis to be developed is the variation of the wavenumber in terms of the material properties. The complex relative permittivity and permeability of the metamaterials (either real or artificial) can be modelled by adopting a simplified Drude's model, which uses the following expressions based in Ref. [38]:

where

One way to describe the behaviour of the wave is to describe the permittivity o the permeability as function of a parameter ratio, and use this relation on the wavenumber k. Let us take for example the permittivity described in Eq. (48) and express it on terms of

**Figure 1** shows the real and imaginary part of the relative permittivity as function of

From **Figure 1(a)** and **(b)**, we can notice the necessary ratios of

After the analysis of the permittivity or permeability on a media than can achieve SNG characteristics, we can express the complex wavenumber *k* of a wave propagating on an infinite media where its permittivity can have negative values, normalized by the wavenumber on free space **Figure 2**). By knowing the complex wavenumber of the wave propagating on a media, we can determine for example the cut‐off frequency of a waveguide filled with the material.

If we consider the wavenumber

When

Equation (35) may be related with the description in terms of massive photons propagating in a plasma: As shown in Ref. [16], the presence of the plasma decreases the rate of electromagnetic energy flow, reaching a zero speed when

### 6.2. Wavenumber on DNG media depending on metamaterial parameters

It is possible to do a further analysis considering now a variation of the permittivity and permeability where both takes negatives values (DNG media). In this case, we will consider a material with its permittivity an permeability following the Drude's model expressed in Eqs. (20) and (21), and its parameters changes equally in terms of the ratios **Figure 3** shows the real and complex values of the wavenumber *k*'’ can result on negative values (remember that here both permittivity and permeability can be negative). As this is not possible, because the imaginary part must be positive due to the conservation of energy, as shown in Ref. [39], there is a change of sign for those values. For these calculations, some values of *k*'’ can result on negative values (remember that here both permittivity and permeability can be negative).

In a more practical way, both ratios are not necessarily modified in the same way. Normally, the periodical structures can present with different parameters of design its negative permittivity or permeability. This means that

Expressing the complex permittivity and complex permeability in terms of **Figure 4** shows the real and imaginary part of the relative permittivity as function of

From **Figure 4**, we can observe when the material have negative or positive values of **Figure 5**, we can trace the wavenumber in function of the different variations of

From this analysis, we can notice that the values of the attenuation and of the phase constant will vary depending of the design parameters of

In conclusion, for the design of metamaterials, it is important to know how the parameters that are described on the Drude's model vary, in order to predict the behaviour of the wave that propagates on the media.

## 7. Discussion

From the point of view of unification of electromagnetic fields and relativistic quantum theory, it is useful to study the dynamics of a photon in a form comparable with the case of a particle like an electron. As a first stage towards this unification, an important result is that the Dirac‐like equation allows to write Maxwell's equations in a compact form and that for light propagating in an homogeneous medium this equation has energy solutions similar to those obtained by Dirac’s equation for fermions in the limit of zero mass, except that in this case there is no a solution with zero energy. Nevertheless, this must be considered only as a formal analogy since photons have spin 1 and Dirac's equation is applicable for particles of spin ½.

Among other issues related with the behaviour of photons that have been a matter of discussion in several publications, it is worth to mention those concerned with localizability, Zitterbewegung and its relation with spin.

The localizability of massless photons was first examined in Ref. [40] by Newton and Wigner and later by Wightman in Ref. [41], showing that there is no position operator for a massless particle with spin higher than ½ leading many authors to conclude that it is not possible to define a wave function for a photon, which has zero mass and spin (or helicity) 1. Wightman has proved that the only localizable massless elementary system has spin zero and that a free photon is not localizable.

However, Bialynicki‐Birula in Refs. [42, 43] and, independently, Sipe in Ref. [44] introduced a function of the position and time coordinates that completely describes the quantum state of a photon. Such function may be referred to as the photon wave function. The wave equation for this function can be derived from the Einstein kinematics for a particle with spin 1 and zero mass in the same way that the Dirac equation is obtained for a massive particle with spin ½. For example, see Refs. [45, 46]. A strong argument in favour of this photon wave function formulation is that the corresponding wave equation is completely equivalent to the Maxwell equations in vacuum. In addition in Ref. [1], a position operator has been postulated whose eigenvectors form bases of localized states.

The concept of spin of a photon and its relation with Zitterbewegung is still a matter of discussion and deserves further research. In Ref. [47], it has been postulated that the spin of the photon can be considered as a consequence of the orbital angular momentum due to the photon's Zitterbewegung. This postulate is based on a Schrödinger‐like equation, having a velocity operator that undergoes oscillations in a direction orthogonal to its momentum, effect known as Zitterbewegung, with a spatial amplitude equal to the classical wavelength. The spin of the photon would be the orbital angular momentum due to the Zitterbewegung. Nevertheless, up to now, this is a theoretical postulate which results from an equation of the form of Schrödinger's equation which is known, was formulated for a nonrelativistic particle. In this context it seems formally more suitable to use a Dirac‐like equation.

## 8. Conclusion

We have presented a conceptual frame for understanding the propagation of light through a dispersive and absorptive medium, considering both the classical description based on electromagnetic waves and a quantum description considering photons as elementary excitations of the quantized electromagnetic field.

A semiclassical description of the dynamics of a photon propagating freely in an unbounded medium has been presented using a Dirac‐like equation, discussing the solutions for the energy and comparing with those corresponding to a free electron as given by the Dirac's equation.

The interaction of light with the medium of propagation has been analysed using both a classical and a quantum treatment. In particular and as a specific example, the first‐order field‐field spatial correlation for a beam of light produced by the excitation of two linearly polarized waves has been calculated, comparing results between the classical and quantum model. For this specific case, it is concluded that in absence of absorption, the classical model predicts total coherence at all points, while the quantum treatment predicts that for some distances between the considered points there is null coherence.

The propagation of light along a rectangular waveguide has been studied showing how the Dirac‐like equation previously studied is modified due to the bounding conditions in the propagation imposed by the conducting walls verifying the result obtained in several publications that in this condition photons appear to acquire an effective mass.

As an application to communication engineering, we have analysed conditions for which the dielectric permittivity and magnetic permeability of a medium filled with plasma behaves as a metamaterial.