Electromagnetic wave equation

From Wikipedia, the free encyclopedia

The electromagnetic wave equation is a second-order partial differential equation that describes the propagation of electromagnetic waves through a medium or in a vacuum. The homogeneous form of the equation, written in terms of either the electric field E or the magnetic field B, takes the form:

  \left( \nabla^2 - { 1 \over {c}^2 } {\partial^2 \over \partial t^2} \right) \mathbf{E} \ \ = \ \ 0
  \left( \nabla^2 - { 1 \over {c}^2 } {\partial^2 \over \partial t^2} \right) \mathbf{B} \ \ = \ \ 0

where c is the speed of light in the medium. In a vacuum, c = c0 = 299,792,458 meters per second, which is the speed of light in free space.

The electromagnetic wave equation derives from Maxwell's equations.

It should also be noted that in most older literature, B is called the "magnetic flux density" or "magnetic induction".

Contents

Speed of propagation

In vacuum

If the wave propagation is in vacuum, then

c = c_o = { 1 \over \sqrt{ \mu_o \varepsilon_o } } = 2.99792458 \times 10^8
\ \mu_0 c_0 \   299\ 792\ 458   \ \varepsilon_0  8.854\ 187\ 817...\times 10^{-12} \begin{matrix}\frac {1}{\mu_0 {c_0}^2}\end{matrix}\  \mu_0 \  4 \pi \times 10^{-7}  \ \Z_0   376.730\ 313\ 461...
Symbol Name Numerical Value SI Unit of Measure Type
speed of light in vacuum meters per second defined
electric constant farads per meter
magnetic constant henries per meter defined
characteristic impedance of vacuum ohms derived; μ0c0

In a material medium

The speed of light in a linear, isotropic, and non-dispersive material medium is

c = { c_0 \over n } =  { 1 \over \sqrt{ \mu \varepsilon } }

where

 n = \sqrt{ \mu \varepsilon \over  \mu_0 \varepsilon_0  }
\ \mu

The origin of the electromagnetic wave equation

Conservation of charge

Conservation of charge requires that the time rate of change of the total charge enclosed within a volume V must equal the net current flowing into the surface S enclosing the volume:

 \oint \limits_S \mathbf{j} \cdot d \mathbf{A}  = - {d \over d t} \int \limits_V \rho \cdot dV

where j is the current density (in Amperes per square meter) flowing through the surface and ρ is the charge density (in coulombs per cubic meter) at each point in the volume.

From the divergence theorem, this relationship can be converted from integral form to differential form:

 \nabla \cdot \mathbf{j} = - { \partial \rho \over \partial t}

Ampère's circuital law prior to Maxwell's correction

In its original form, Ampère's circuital law relates the magnetic field B to the current density j:

 \oint \limits_C \mathbf{B} \cdot d \mathbf{l} =  \iint \limits_S \mu \mathbf{j} \cdot d \mathbf{A}

where S is an open surface terminated in the curve C. This integral form can be converted to differential form, using Stokes' theorem:

 \nabla \times \mathbf{B} = \mu_0 \mathbf{j}

Inconsistency between Ampère's circuital law and the law of conservation of charge

Taking the divergence of both sides of Ampère's circuital law gives:

 \nabla \cdot  ( \nabla \times \mathbf{B} ) = \nabla \cdot \mu_0 \mathbf{j}

The divergence of the curl of any vector field, including the magnetic field B, is always equal to zero:

 \nabla \cdot  ( \nabla \times \mathbf{B}) = 0

Combining these two equations implies that

\nabla \cdot \mu_0 \mathbf{j} = 0
 \ \mu_0
\nabla \cdot \mathbf{j} = 0

However, the law of conservation of charge tells that

 \nabla \cdot \mathbf{j} = - { \partial \rho \over \partial t }

Hence, as in the case of Kirchhoff's circuit laws, Ampère's circuital law would appear only to hold in situations involving constant charge density. This would rule out the situation that occurs in the plates of a charging or a discharging capacitor.

Maxwell's correction to Ampère's circuital law

Gauss's law in integral form states:

 \oint \limits_S \mathbf{E} \cdot d \mathbf{A}  = \frac{1}{\varepsilon_0} \int \limits_V \rho \cdot dV \ ,

where S is a closed surface enclosing the volume V. This integral form can be converted to differential form using the divergence theorem:

 \nabla \cdot \varepsilon_0 \mathbf{E}  =  \rho

Taking the time derivative of both sides and reversing the order of differentiation on the left-hand side gives:

 \nabla \cdot   \varepsilon_0   {\partial  \mathbf{E}   \over \partial t }     = { \partial \rho \over \partial t}

This last result, along with Ampère's circuital law and the conservation of charge equation, suggests that there are actually two origins of the magnetic field: the current density j, as Ampère had already established, and the so-called displacement current:

   {\partial  \mathbf{D}   \over \partial t }   =  \varepsilon_0   {\partial  \mathbf{E}   \over \partial t }

So the corrected form of Ampère's circuital law becomes:

 \nabla \times \mathbf{B} = \mu_0 \mathbf{j} + \mu_0 \varepsilon_0   {\partial  \mathbf{E}   \over \partial t }

Maxwell's hypothesis that light is an electromagnetic wave

A postcard from Maxwell to Peter Tait. A postcard from Maxwell to Peter Tait.

In his 1864 paper entitled A Dynamical Theory of the Electromagnetic Field, Maxwell utilized the correction to Ampère's circuital law that he had made in part III of his 1861 paper On Physical Lines of Force. In PART VI of his 1864 paper which is entitled 'ELECTROMAGNETIC THEORY OF LIGHT', Maxwell combined displacement current with some of the other equations of electromagnetism and he obtained a wave equation with a speed equal to the speed of light. He commented:

The agreement of the results seems to show that light and magnetism are affections of the same substance, and that light is an electromagnetic disturbance propagated through the field according to electromagnetic laws.

Maxwell's derivation of the electromagnetic wave equation has been replaced in modern physics by a much less cumbersome method involving combining the corrected version of Ampère's circuital law with Faraday's law of induction.

To obtain the electromagnetic wave equation in a vacuum using the modern method, we begin with the modern 'Heaviside' form of Maxwell's equations. In a vacuum, these equations are:

 \nabla \cdot \mathbf{E} = \frac {\rho} {\epsilon_0}
 \nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}} {\partial t}
 \nabla \cdot \mathbf{B} = 0
 \nabla \times \mathbf{B} =\mu_0 \varepsilon_0 \frac{ \partial \mathbf{E}} {\partial t}

Taking the curl of the curl equations gives:

 \nabla \times \nabla \times \mathbf{E} = -\frac{\partial } {\partial t} \nabla \times \mathbf{B} = -\mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{E} }  {\partial t^2}
 \nabla \times \nabla \times \mathbf{B} = \mu_0 \varepsilon_0 \frac{\partial } {\partial t} \nabla \times \mathbf{E} = -\mu_o \varepsilon_o \frac{\partial^2 \mathbf{B}}{\partial t^2}

By using the vector identity

\nabla \times \left( \nabla \times \mathbf{V} \right) = \nabla \left( \nabla \cdot \mathbf{V} \right) - \nabla^2 \mathbf{V}
 \mathbf{V}
 {\partial^2 \mathbf{E} \over \partial t^2} \ - \  {c_0}^2 \cdot \nabla^2 \mathbf{E}  \ \ = \ \ 0
 {\partial^2 \mathbf{B} \over \partial t^2} \ - \  {c_0}^2 \cdot \nabla^2 \mathbf{B}  \ \ = \ \ 0

where

c_0 = { 1 \over \sqrt{ \mu_0 \varepsilon_0 } } = 2.99792458 \times 10^8

is the speed of light in free space.

Covariant form of the homogeneous wave equation

Time dilation in transversal motion. The requirement that the speed of light is constant in every inertial reference frame leads to the theory of Special RelativityTime dilation in transversal motion. The requirement that the speed of light is constant in every inertial reference frame leads to the theory of Special Relativity

These relativistic equations can be written in covariant form as

\ \Box A^{\mu} = 0

where the electromagnetic four-potential is

A^{\mu}=(\varphi, \mathbf{A} c)

with the Lorenz gauge condition:

\partial_{\mu} A^{\mu} = 0\,

Here

\Box = \nabla^2 - { 1 \over c^2} \frac{   \partial^2} { \partial t^2}

Homogeneous wave equation in curved spacetime

Main article: Maxwell's equations in curved spacetime

The electromagnetic wave equation is modified in two ways, the derivative is replaced with the covariant derivative and a new term that depends on the curvature appears.

 - {A^{\alpha ; \beta}}_{; \beta} + {R^{\alpha}}_{\beta} A^{\beta} = 0

where

  {R^{\alpha}}_{\beta}

is the Ricci curvature tensor and the semicolon indicates covariant differentiation.

The generalization of the Lorenz gauge condition in curved spacetime is assumed:

  {A^{\mu}}_{ ; \mu} =0

Inhomogeneous electromagnetic wave equation

Main article: Inhomogeneous electromagnetic wave equation

Localized time-varying charge and current densities can act as sources of electromagnetic waves in a vacuum. Maxwell's equations can be written in the form of a wave equation with sources. The addition of sources to the wave equations makes the partial differential equations inhomogeneous.

Solutions to the homogeneous electromagnetic wave equation

Main article: Wave equation

The general solution to the electromagnetic wave equation is a linear superposition of waves of the form

 \mathbf{E}( \mathbf{r}, t )  =  g(\phi( \mathbf{r}, t ))  =  g( \omega t  -  \mathbf{k} \cdot \mathbf{r}   )

and

 \mathbf{B}( \mathbf{r}, t )  =  g(\phi( \mathbf{r}, t ))  =  g( \omega t  -  \mathbf{k} \cdot \mathbf{r}   )

for virtually any well-behaved function g of dimensionless argument φ, where

 \ \omega  \mathbf{k} = ( k_x, k_y, k_z)

Although the function g can be and often is a monochromatic sine wave, it does not have to be sinusoidal, or even periodic. In practice, g cannot have infinite periodicity because any real electromagnetic wave must always have a finite extent in time and space. As a result, and based on the theory of Fourier decomposition, a real wave must consist of the superposition of an infinite set of sinusoidal frequencies.

In addition, for a valid solution, the wave vector and the angular frequency are not independent; they must adhere to the dispersion relation:

 k = | \mathbf{k} | = { \omega \over c } =  { 2 \pi \over \lambda }

where k is the wavenumber and λ is the wavelength.

Monochromatic, sinusoidal steady-state

The simplest set of solutions to the wave equation result from assuming sinusoidal waveforms of a single frequency in separable form:

\mathbf{E} ( \mathbf{r}, t ) = \mathrm {Re} \{ \mathbf{E} (\mathbf{r} )  e^{ j \omega t }  \}

where

Plane wave solutions

Main article: Sinusoidal plane-wave solutions of the electromagnetic wave equation

Consider a plane defined by a unit normal vector

 \mathbf{n} = { \mathbf{k} \over k }

Then planar traveling wave solutions of the wave equations are

 \mathbf{E}(\mathbf{r}) = E_0 e^{-j \mathbf{k} \cdot \mathbf{r} }

and

 \mathbf{B}(\mathbf{r}) = B_0 e^{-j \mathbf{k} \cdot \mathbf{r} }

where

 \mathbf{r} = (x, y, z)
 \mathbf{n}
    c  {\partial B \over \partial z} = {\partial E \over \partial t}

Because the divergence of the electric and magnetic fields are zero, there are no fields in the direction of propagation.

This solution is the linearly polarized solution of the wave equations. There are also circularly polarized solutions in which the fields rotate about the normal vector.

Spectral decomposition

Because of the linearity of Maxwell's equations in a vacuum, solutions can be decomposed into a superposition of sinusoids. This is the basis for the Fourier transform method for the solution of differential equations. The sinusoidal solution to the electromagnetic wave equation takes the form

Electromagnetic spectrum illustration. Electromagnetic spectrum illustration.
 \mathbf{E} ( \mathbf{r}, t ) = \mathbf{E}_0 \cos( \omega t  -  \mathbf{k} \cdot \mathbf{r} + \phi_0  )

and

 \mathbf{B} ( \mathbf{r}, t ) = \mathbf{B}_0 \cos(  \omega t  -  \mathbf{k} \cdot \mathbf{r} + \phi_0  )

where

 \ t  \ \omega  \mathbf{k} = ( k_x, k_y, k_z)  \phi_0 \,

The wave vector is related to the angular frequency by

 k = | \mathbf{k} | = { \omega \over c } =  { 2 \pi \over \lambda }

where k is the wavenumber and λ is the wavelength.

The electromagnetic spectrum is a plot of the field magnitudes (or energies) as a function of wavelength.

Other solutions

Spherically symmetric and cylindrically symmetric analytic solutions to the electromagnetic wave equations are also possible.

In cylindrical coordinates the wave equation can be written as follows:

 \mathbf{E} ( \mathbf{r}, t ) = {\mathbf{E}_0 \cos( \omega t  -  \mathbf{k} \cdot \mathbf{r} + \phi_0  )\over s}

and

 \mathbf{B} ( \mathbf{r}, t ) = {\mathbf{B}_0 \cos(  \omega t  -  \mathbf{k} \cdot \mathbf{r} + \phi_0  )\over s}