Review of electromagnetic fields



1. Review of electromagnetic fields

1.1 Maxwell’s Equations

• Review the features of Maxwell’s Equations in differential forms.

• Discuss these equations in different representations: space-time (r, t), space-frequency (r, ω), wavevector-frequency (k, ω).

1.1.1 Maxwell’s Eqs. in the space-time representation

• The first equation, Faraday’s induction law, establishes that a varying magnetic field produces a rotating electric field,

[pic] (1.1)

• E(r, t) is the electric field (V/m) and B(r,t) the magnetic induction or magnetic flux density (Wb/m2).

• Ampere’s circuital law states a magnetic field can be generated by existing currents and by varying electric fields,

[pic] (1.2)

• H(r,t) is the magnetic field intensity (A/m), D(r,t) the electric induction or flux density (C/m2), Js(r,t) the source current density (A/m2), and Jc(r,t) the conduction current density (A/m2).

• Gauss’s law for the electric and magnetic fields:

[pic] (1.3)

• ρ(r, t) is the volume charge density (C/m3), and

[pic], (1.4)

• which establishes the non-existence of magnetic charge.

• In addition, the following constitutive relations apply:

[pic] (1.5)

• ε the dielectric permittivity, μ the magnetic permeability, and ( the conductivity

• Generally, ε, μ, and σ, are inhomogeneous (r-dependent) and anisotropic (direction-dependent, tensors). They can also be field-dependent for nonlinear media.

• For linear, homogeneous, and isotropic media, ε, μ, and σ are simple scalars. The relative permittivity ((r and permeability μr are defined with respect to the vacuum values, ε0=10-9/36π F/m and [pic]N/A2.

1.1.2 Boundary conditions

• To solve Maxwell’s equations for fields propagating between two media, boundary conditions for both the electric and magnetic fields are necessary.

[pic]

Figure 1. Fields at the boundary between two media

• 4 equations describe the relationship between the tangent components of E and H and the normal components of D and B in the two media (Fig. 1), as follows:

[pic] (1.6)

• n is the unit vector normal to the interface, Js is the surface current, and ρs is the surface charge density. In the absence of free charge and currents, the boundary conditions simply imply that the tangential fields (E, H) and normal inductions (D, B) are conserved across an interface.

1.1.3 Maxwell’s equations in the space-frequency representation (r, ()

• For broad-band fields, the spatial behavior of each temporal frequency is of interest. To obtain Maxwell’s equation in the space-frequency representation, we use the differentiation property of Fourier transforms,

[pic] (1.7)

• F is the Fourier transform operator (in time), F(t) is an arbitrary time-dependent vector, and F(ω) its Fourier transform.

• We will denote the Fourier transform of a function by the same symbol, with the understanding that F(t) and F(ω) are two distinct functions.

• Taking the Fourier transform of Eqs. 1-4 and using the differentiation property in Eq. 7, we can rewrite the Maxwell’s equations in the (r, ω) representation,

[pic] (1.8)

• The constitutive relations take the form

[pic] (1.9)

• Even for homogeneous media ε is a function of frequency, which establishes the dispersion relation associated with the medium. The boundary conditions apply for each individual frequency,

[pic] (1.10)

1.1.4 The Helmholtz equation

• Eliminating H, B, and D from Eqs. 1-5, we obtain an equation in E(r, t), the wave equation. The Helmholtz equation can be obtained in the (r, ω) representation, as follows.

• First let us consider linear, isotropic, and charge-free media ((=0). Eqs. 8a-d simplify to

[pic] (1.11)

• Applying the curl operator to Eq. 11a, we obtain

[pic] (1.12)

• Using the identity[pic] and the fact that [pic]for charge-free media, we obtain the Helmholtz equation:

[pic] (1.13)

• In Eq. 13a, the Laplace operator applied to vector E is defined as [pic].

• Equation 13b establishes the frequency dependence of the wavenumber ((ω), the dispersion relation of the medium. Typical biological structures are characterized by constant μ, the permittivity ε and conductivity σ can have strong dependence on frequency.

• In Cartesian coordinates, the vector Eq. 13a becomes three scalar equations.

[pic] (1.14)

• All the equations are of the same form, referred to as the scalar wave equation,

[pic].

• Assuming solutions of the form [pic], we obtain three 1D equations, which eventually lead to the plane wave solution of the form

[pic] (1.15)

• Decomposing the wavevector into its real and imaginary parts,

[pic], (1.16)

• We can rewrite Eq. 15a as

[pic] (1.17)

• In Eq. 17 the real and imaginary parts of k capture the refraction (phase term) and absorption/gain (amplitude term) of the medium.

• In lossless media, [pic]. In dielectric media, at optical frequencies[pic] and [pic]equals the refractive index, thus

[pic] (1.18)

• [pic]is the wavenumber in vacuum.

1.1.5 Maxwell’s equations in the (k, ω) representation

• Often we deal with optical fields with broad angular spectra. These can be decomposed in wavevectors k of different directions. The modulus of each k-vector is defined by the dispersion relation, depending on the material properties (ε, μ, σ) and optical frequency ω.

• The normal representation of the fields is in the k-ω space. A differentiation property similar to Eq. 7 holds for the ( operator,

[pic] (1.19)

• The (k, ω) representation of Maxwell’s Equations is obtained by the spatial Fourier transformation of Eqs. 11a-d.

• For media of no free charge (ρ=0) or currents (J=0), these equations are

[pic] (1.20)

• Equations 20a-d describe the propagation of frequency component ω and plane wave of wavevector k. Eqs. 20c-d establish that [pic] and [pic]. Generally μ is a scalar, i.e. B||H, but ε is a tensor, i.e. D is not necessarily parallel to E.

• Thus we see that[pic], such that

[pic] (1.21)

• Equations 20a-c show that H, D, and k are mutually orthogonal vectors.

[pic]

Figure 1-2. Mutually orthogonal set of vectors

• For isotropic media, D||E, such that H, E, and k are also mutually orthogonal.

• The characteristic impedance [pic] of the medium is defined as

[pic] (1.22)

• Using Eq. 15b to express k for media with non-zero conductivity, we obtain

[pic] (1.23)

• For non-conducting media (electric insulators),[pic]. For isotropic media, [pic] connects the moduli of E and H directly (from Eq. 19b),

[pic] (1.24)

1.1.6 Phase, group, and energy velocity

[pic]

Figure 3. Phase and group velocity.

• Consider the electric field associated with a light beam propagating along the +z direction with an average wavevector .

• As shown in Fig. 3, the temporal signal has a slow modulation (envelope) due to the superposition of different frequencies, and a fast sinusoidal modulation (carrier), at the average frequency ,

[pic] (1.25)

• Thus the phase delay of the field is

[pic] (1.26)

• The phase velocity is associated with the advancement of wave fronts, for which ((=constant. Differentiating Eq. 26, we obtain

[pic] (1.27)

• Thus, the phase velocity is

[pic] (1.28)

• Equation 28 represents the zeroth order approximation in the expansion:

[pic] (1.29)

• Equation 29 defines the group velocity, vg, at which the envelope of the light signal propagates. The group and phase velocities are equal in non-dispersive media ([pic]).

• We prove that the group velocity defines the energy velocity of the field. The energy velocity is the ratio between the Poynting vector S and electromagnetic volume density U,

[pic] (1.30)

• We now differentiate Eqs. 20a-b, as follows

[pic] (1.31)

• The meaning of δk and δω is that of a spread in k-vectors (directions) and optical frequency. Dot-multiply Eq. 31a by H and Eq. 31b by E,

[pic] (1.32)

• Adding the two Eqs. 32a-b, we obtain

[pic] (1.33)

• using the fact that the products [pic] and [pic] vanish.

• Finally, using the definition in Eq. 30a, we find

[pic] (1.34)

• Equation 34 establishes the result that the electromagnetic energy flows at group velocity; it is apparent from Eq. 29 that vg can exceed the speed of light in vacuum c, in special circumstances, where [pic] (anomalous dispersion).

• This does not pose a conflict with the postulate of the relativity theory, which states that the signal velocity, at which information can be transmitted via electromagnetic fields, is bounded by c.

1.1.7 The Fresnel Equations

• Consider light propagation at the interface between two media (Fig. 4). We will derive the expressions for the field reflection and transmission coefficients.

[pic]

Figure 4. Two media of refractive index n1 and n2 separated by the x-y plane. The subscripts i, t, r refer to incident, transmitted, and reflected.

• In Figure 4, the x-z plane is the plane of incidence (plane of the paper), i.e. the plane defined by the incident wavevector and normal ([pic]) at the interface.

• The tangential components of the fields and normal components of inductions are conserved in the absence of surface currents and charges (Eqs. 6a-d),

[pic] (1.35)

• For ρ=0 and J=0, Maxwell’s Eqs. in the k-ω representation yield (Eqs. 19)

[pic] (1.36)

• Expanding the cross products in Eqs. 36a-b, the problem breaks into two independent cases: a) transverse electric (TE) mode, when E is perpendicular to the plane of incidence (E||y) and b) transverse magnetic (TM) mode, when H||y.

i) TE mode (E||y)

• If E||y, the boundary conditions for the tangent E-field and normal B-field

[pic] (1.37)

• Using Eq. 36b to express Eq. 37b in terms of Ey components, we can rewrite the system of equations as

[pic] (1.38)

• Eqs. 38a-b must hold for any incident field Eyi, we obtain the following result

[pic] (1.39)

• Eq. 39 establishes Snell’s law,

[pic] (1.40)

• using that the wavevector k in a medium of refractive index n relates to the wavevector in vacuum as k=nk0.

• Eq. 40 implies kzi=-kzn.

• To obtain the field reflection coefficient, we use the continuity of tangent H components,

[pic] (1.41)

• which can be expressed in terms of E components (via Eq. 36b),

[pic] (1.42)

• Finally, combining Eqs. 42 and 37a to solve for the E field transmission and reflection coefficients, we obtain

[pic] (1.43)

ii) TM mode (H||y)

• Using the analog equations to the TE mode (Eqs. 37a-b), the conservation of Hy components and normal D components, we find that kx=const. for TM as well.

• To obtain the field reflection and transmission coefficients, we use the conservation of both the tangents fields components, Ex and Hy,

[pic] (1.44)

• The [pic] factor occurs due to the 1/ε factor in Eq. 36b (ε=n2). Thus the H field reflection and transmission coefficients for the TM mode are

[pic] (1.45)

• We expressed rTM and tTM in terms of H fields to emphasize the symmetry with respect to the TE case. Of course, the quantities can be further expressed in terms of E fields via [pic]. Conservation of energy is satisfied in both cases,

[pic] (1.46)

• Together, Eqs. 43a-b and 45a-b, the Fresnel equations, provide the reflected and transmitted fields for an arbitrary incident field. Because of the polarization dependence of the reflection and refraction coefficient, polarization properties of light can be modified via reflection and refraction. In the following we discuss two particular cases that follow from the Fresnel equations, where the transmission or reflection coefficients vanish.

1.1.8 Total internal reflection

• Setting tTE=tTM=0 yields the same condition for “no transmission” in both TE and TM modes, transmission coefficient vanishes, kzt=0. Thus,

[pic] (1.47)

• using kx=constant (Eq. 38), kxt=kxi. The transmission vanishes for

[pic] (1.48)

• θc is the critical angle at which total internal reflection takes place. Total internal reflection can occur for both TE and TM polarizations, the only restriction being that n2θc, the field reflection coefficient becomes

[pic] (1.49)

• [pic].

• For θi>θc, the reflection coefficient is purely imaginary, the power is 100% reflected, but the reflected field is shifted in phase by 2(TE.

• Similarly, for the TM mode we obtain

[pic] (1.50)

• Since (TM and (TE have different values, total internal reflection can be used to change the polarization state of optical fields.

• The transmitted plane wave has the form

[pic] (1.51)

• Equation 51 indicates that the field in medium 2 is decaying exponentially.

[pic]

Figure 1-5. Evanescent field decaying exponentially with depth z.

• Thus the field is significantly attenuated over a distance on the order of 1/kzt, i.e. the field does not propagate, or is evanescent.

1.1.9 Transmission at Brewster angle

• Another case of Fresnel’s equations is when the reflection coefficient vanishes.

• For the TE mode, we have

[pic] (1.52)

• For TE polarization, the only way to obtain maximum transmission through an interface is when there is no refractive index contrast between the two media, the trivial solution.

• For the TM mode the situation is very different:

[pic] (1.53)

• The condition is satisfied simultaneously with Snell’s law, such that we have

[pic] (1.54)

• Multiplying Eqs. 54a and 54b side by side, we obtain

[pic] (1.55)

• The angle of incidence at which rTM=0, referred to as the Brewster angle, is defined by combining Eq. 55 and Eq. 54a

[pic] (1.56)

• Unlike with total internal reflection, where the transmission can vanish for both polarizations, the reflection can only vanish in the TM mode.

• The absence of reflection at the Brewster angle for TM polarization can be understood by the absence of radiation by an (induced) dipole along its axis (Fig. 6). The concept of induced dipoles followed by re-radiation is essential for the Lorentz model of light-matter interaction, as detailed in the next section.

[pic]

Figure 1-6. Brewster angle incidence.

1.2 The Lorentz Model of Light Matter Interaction

• Review the main concepts in basic atom-field interactions. In particular the Lorentz model, a pre-quantum mechanics model, and its asymptotic case for metals, the Drude model.

• The Lorentz model explains much of classical optics via a physical picture borrowed from mechanics. The starting point is the “mass on a spring” description of electrons connected to nuclei. Thus, the incident electric field induces displacement to the electron that is under the influence of a spring-like restoring force due to the nucleus.

• The equation of motion for the electron can be expressed as

[pic] (1.57)

• ( is the damping constant, (0 is the resonant frequency, e is the electronic charge, m mass of the electron, and E the incident field. ( and (0 are characteristics of the material, the first describing the energy dissipation property of the medium and the second the ability of the medium to store energy.

• Since Eq. 57 is a linear differential equation, Fourier transforming both sides of the equation gives the frequency-domain solution.

• Using the Fourier property of the differential operator,

[pic] (1.58)

• we obtain Eq. 57 in the frequency domain

[pic]. (1.59)

• Thus, we find the solution for the charge displacement in the frequency domain,

[pic]. (1.60)

• To obtain the time domain solution, x(t), we need to Fourier-transform Eq. 60. However, we explore further the frequency domain solution. The induced dipole moment due to the charge displacement x(() is

[pic] (1.61)

• In Eqs. 59-60 we obtained microscopic quantities, the atomic-level response.

• The macroscopic behavior of the medium is obtained from the induced polarization P, which captures the contribution of all dipole moments within a certain volume,

[pic] (1.62)

• N is the volume concentration of dipoles (m-3) and the angular brackets denote ensemble average.

• Assuming that all induced dipoles are parallel within the volume, we obtain

[pic] (1.63)

• Generally, each atom has multiple resonances or dipole-active modes, such that Eq. 63 can be generalized to

[pic] (1.64)

• The summation is over all modes, characterized by different resonant frequencies and damping constants. The weight (i is called the oscillator strength and has the quantum mechanical meaning of a transition strength.

• For simplicity, we reverse to the single normal mode description, which captures the origin of absorption and refraction of materials. The induced polarization only captures the contribution of the medium itself, it excludes the vacuum contribution.

• Thus,

[pic] (1.65)

• ( is the dielectric susceptibility, which generally is a tensor quantity.

• However, for isotropic media, we obtain the complex scalar permittivity [pic]

[pic] (1.66)

• [pic] has units of frequency squared, [pic], and (p is the plasma frequency.

• From Eq. 66, we readily obtain the real and imaginary parts of (r,

[pic] (1.67)

• Figure 7 illustrates the main features of [pic] and [pic] vs. frequency. To gain further physical insight into Eqs. 67a-b, we discuss three different frequency regions, as follows.

1.2.1 Below the resonance, [pic]

• In this case, Eqs. 67a-b simplify to

[pic] (1.68)

• Since [pic], [pic], absorption is negligible, below the resonance the material is transparent. Further, [pic], which defines a region of normal dispersion.

• It can be seen that expanding the denominator in Eqs. 68a-b, we obtain [pic] and [pic]. In designing optics of imaging systems, the Sellmeier equation is very efficient for describing the refractive index vs. wavelength,

[pic], (1.68a)

• The summation is over several resonances, ai and bi are experimentally determined parameters, and c is the speed of light. It can be seen that the Sellmeier equation originates in the expression for [pic] in Eq. 1.68.

• As we approach resonance, this dependence becomes more complicated.

1.2.2 At resonance, [pic]

• For frequencies comparable to [pic], Eqs. 67a-b are well approximated by

[pic] (1.69)

• Under these conditions, the absorption is significant, [pic], and the absorption line has a characteristic shape, Lorentzian line. This shape has a central frequency [pic] and a full width half maximum of [pic].

• While [pic] has a clear physical significance of the frequency at which the system “resonates”, or absorbs strongly, the meaning of [pic] is somewhat more subtle.

• The damping constant [pic] represents the average frequency at which electrons collide with atoms, which generates loss of energy. Thus, [pic], with [pic] the average time between collisions.

• Finally, around resonance, [pic], which defines anomalous dispersion.

1.2.3 Above the resonance, [pic]

• Well above the resonance, the following equations apply:

[pic] (1.70)

• The absorption becomes less significant, as expected in a frequency range away from the resonance. The dispersion is normal again, [pic].

• This Lorentz oscillatory model provides great insight into the classical light-matter interaction. In the following section, we will investigate the particular situation of metals, when the charge moves freely within the material.

1.3 Drude model of light-metals interaction

• The optical properties of metals were first introduced by Drude in the context of conductivity. In highly conductive materials, the restoring force in Eq. 57, [pic], vanishes, establishing that the charge can move freely. Under these conditions, we obtain Drude’s model, in which Eqs. 67a-b reduce to [pic]

[pic] (1.71)

• Typically [pic], the frequency of collisions is much lower than that of optical frequencies.

• In this high frequency limit, [pic] and [pic]. From the Fresnel equations, we derive the reflectivity coefficient. For normal incidence, the intensity-based reflectivity is

[pic] (1.72)

• n is the (complex) refractive index. Since [pic], Eq. 72 becomes

[pic] (1.73)

• Figure 7 illustrates the frequency dependence of [pic], [pic], n’, n’’ and R for various values of [pic] and [pic].

[pic]

Figure 7. Frequency dependence of dielectric permeability and reflectivity around the plasma frequency ((p =10, (=1) .

• At [pic], [pic] vanishes. In this case, the real part of the refractive index, n’, can also vanish. This implies that the wavelength in the material is infinite, [pic].

• To gain a physical understanding of the plasma frequency [pic], consider a thin film of metal.

[pic]

Figure 1-9. Exiting surface plasmon resonance.

• The applied electric field induces a polarization [pic], with [pic].

• Tuning the frequency of the incident field to the plasma frequency, [pic], [pic], and [pic]. The induced polarization is the total charge times the displacement per unit volume,

[pic] (1.74)

• Therefore, the electric field is

[pic] (1.75)

• If we construct the electric force due to the charge displacement, [pic], we obtain

[pic] (1.76)

• In Eq. 76, we define ke as the “spring” constant of the restoring force. By definition, the system is characterized by a resonant frequency, [pic]. This is the plasma frequency associated with the thin film,

[pic]. (1.77)

• From Maxwell’s equations, we have that [pic]. At plasma frequency, [pic], and the magnetic field vanishes. This indicates that there is no bulk propagation of electromagnetic field.

• Since the E-field must be in the plane of the surface (plane of incidence), it is clear that only TM polarization can induce surface plasma resonance. The quantum of energy associated with the charge oscillations at plasma frequency is [pic] and the respective quantum particle is called plasmon.

• Light scattering on plasma oscillations, including in nanoparticles (Figure 9) can be described as a photon-plasmon interaction. These interactions are the basis for many modern sensors and lab-on-a-chip devices.

[pic]

Figure 1-10. Plasmon resonance in nanoparticles.

................
................

In order to avoid copyright disputes, this page is only a partial summary.

Google Online Preview   Download