Electromagnetic stress tensor in ponderable media

doi:10.1364/OE.16.005193

« Show journal navigation

Electromagnetic stress tensor in ponderable media

Masud Mansuripur »View Author Affiliations

Optics Express, Vol. 16, Issue 8, pp. 5193-5198 (2008)
http://dx.doi.org/10.1364/OE.16.005193

View Full Text Article

Acrobat PDF (113 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Stress tensor in ponderable medium
3. Concluding remarks
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (6)
Cited By
Figures (3)
Metrics

Abstract

We derive an expression for the Maxwell stress tensor in a magnetic dielectric medium specified by its permittivity ε and permeability µ. The derivation proceeds from the generalized form of the Lorentz law, which specifies the force exerted by the electromagnetic E and H fields on the polarization P and magnetization M of a ponderable medium. Our stress tensor differs from the well-known tensors of Abraham and Minkowski, which have been at the center of a century-old controversy surrounding the momentum of the electromagnetic field in transparent materials.

1. Introduction

In the classical theory of electromagnetism, the Poynting vector S ( r , t)= E × H specifies the rate of flow of energy at a given location in space r and time t. The integral of S ( r , t) over a closed surface is equal to the rate of increase of stored energy, ∂ ε /∂t, within the volume enclosed by the surface, plus the rate of loss (or minus the rate of gain) of energy throughout that volume [1

J. D. Jackson, Classical Electrodynamics , 2^nd edition, Wiley, New York, 1975.

]. Similarly, the stress tensor T ( r , t) specifies the rate of flow of momentum (i.e., momentum crossing unit area per unit time) at a given point in space and time.

The electromagnetic stress tensor in the free space, known as Maxwell’s stress tensor, is readily derived from the Lorentz law of force in conjunction with the microscopic form of Maxwell’s equations [1

J. D. Jackson, Classical Electrodynamics , 2^nd edition, Wiley, New York, 1975.

]. However, in ponderable media where polarization density P and magnetization density M describe the electromagnetic properties of the material (and where the macroscopic version of Maxwell’s equations, incorporating P and M , are applicable), the form of the stress tensor has been the subject of debate and controversy for the past century [2

R. N. C. Pfeifer, T. A. Nieminen, N. R Heckenberg, and H. Rubinsztein-Dunlop, “Momentum of an electromagnetic wave in dielectric media,” Rev. Mod. Phys. 79, 1197–1216 (2007). [CrossRef]

]. In particular, the tensors of Minkowski and Abraham have been subjected to endless analysis; experiments have been designed to confirm or refute their corresponding predictions, and each tensor, while successful in some respects, has been found inadequate in others.

In a series of recent publications, we have shown that the pressure and momentum of the electromagnetic field can be obtained directly from the Lorentz law of force in conjunction with the macroscopic Maxwell equations, without the need to determine the stress tensor. However, like Poynting’s vector, stress tensor is a powerful tool that can be used to calculate forces and momenta by employing the knowledge of the fields at the boundaries, without having to pay attention to the details of what goes on inside a volume. In this paper we derive a general expression for the electromagnetic stress tensor inside ponderable media and show that, in cases that we have studied previously, the new tensor reproduces the old results.

2. Stress tensor in ponderable medium

In a recent publication [3

M. Mansuripur, “Radiation pressure and the linear momentum of the electro-magnetic field in magnetic media,” Opt. Express 15, 13502–13518 (2007). [CrossRef] [PubMed]

] we derived the following generalized expression for the Lorentz force density in a linear isotropic medium specified by its µ and ε parameters:

F (r, t) = (P \cdot \nabla) E + (M \cdot \nabla) H + (\frac{\partial P}{\partial t}) \times μ_{o} H - (\frac{\partial M}{\partial t}) \times ε_{o} E .

(1)

In conjunction with Eq. (1), Maxwell’s equations in the MKSA system of units are:

\nabla \cdot D = ρ_{free},

(2a)

\nabla \times H = J_{free} + \frac{\partial D}{\partial t},

(2b)

\nabla \times E = - \frac{\partial B}{\partial t},

(2c)

\nabla \cdot B = 0 .

(2d)

In what follows, the medium will be assumed to have neither free charges nor free currents, that is, ρ _free=0 and J _free=0. In the above equations, the electric displacement D and the magnetic induction B are related to the polarization density P and the magnetization density M as follows:

D = ε_{o} E + P = ε_{o} (1 + χ_{e}) E = ε_{o} ε E,

(3a)

B = μ_{o} H + M = μ_{o} (1 + χ_{m}) H = μ_{o} μ H .

(3b)

Using different arguments, Hansen and Yaghjian [4

T. B. Hansen and A. D. Yaghjian, Plane-Wave Theory of Time-Domain Fields: Near-Field Scanning Applications , IEEE Press, New York (1999). [CrossRef]

] have arrived at the same expression as Eq. (1) for the Lorentz force under quite general conditions. Also, Kemp et al, in their analysis of momentum in left-handed media [5

B. A. Kemp, J. A. Kong, and T. Grzegorczyk, “Reversal of wave momentum in isotropic left-handed media,” Phys. Rev. A 75, 053810 (2007). [CrossRef]

], use essentially the same force equation. Our derivation of Eq. (1) in [3

M. Mansuripur, “Radiation pressure and the linear momentum of the electro-magnetic field in magnetic media,” Opt. Express 15, 13502–13518 (2007). [CrossRef] [PubMed]

] started from the well-known Lorentz formula, F =q( E + V × B ), but it soon became apparent that magnetic dipoles cannot be treated as simple Amperian current loops; conservation of momentum demanded certain modifications of the original Lorentz law. In particular, a new term had to be introduced to account for the force experienced by magnetic dipoles. Rather than attempting to justify Eq. (1) on the basis of the original Lorentz law, we believe that one should simply accept it as a law of nature, on par with Maxwell’s equations. Not only are these five equations consistent among themselves, they also comply with the laws of energy and momentum conservation.

The first term on the right hand side of Eq. (1) may be rewritten using the identities

(P \cdot \nabla) E + (\nabla \cdot P) E = \frac{\partial (P_{x} E)}{\partial x} + \frac{\partial (P_{y} E)}{\partial y} + \frac{\partial (P_{z} E)}{\partial z},

(4a)

\nabla \cdot P = \nabla \cdot (D - ε_{o} E) = ρ_{free} - ε_{o} \nabla \cdot E = - ε_{o} \nabla \cdot E .

(4b)

A similar treatment can be applied to the second term in Eq. (1); here ∇· B is readily set to zero in accordance with Maxwell’s 4^th equation.

(M \cdot \nabla) H + (\nabla \cdot M) H = \frac{\partial (M_{x} H)}{\partial x} + \frac{\partial (M_{y} H)}{\partial y} + \frac{\partial (M_{z} H)}{\partial z} .

(5a)

\nabla \cdot M = \nabla \cdot (B - μ_{o} H) = \nabla . B - μ_{o} \nabla \cdot H = - μ_{o} \nabla \cdot H .

(5b)

The third term in Eq. (1) is rewritten by substituting for P in terms of D and E , then invoking Maxwell’s 2^nd equation. Similarly, the fourth term is rewritten by substituting for M in terms of B and H , then invoking Maxwell’s 3^rd equation. We find

(\frac{\partial P}{\partial t}) \times μ_{o} H = (\frac{\partial D}{\partial t} - \frac{ε_{o} \partial E}{\partial t}) \times μ_{o} H = μ_{o} (\nabla \times H) \times H - ε_{o} μ_{o} (\frac{\partial E}{\partial t}) \times H,

(6a)

(\frac{\partial M}{\partial t}) \times ε_{o} E = (\frac{\partial B}{\partial t} - \frac{μ_{o} \partial H}{\partial t}) \times ε_{o} E = {- ε}_{o} (\nabla \times E) \times E - ε_{o} μ_{o} (\frac{\partial H}{\partial t}) \times E .

(6b)

Substitution from Eqs. (4–6) into Eq.(1), followed by rearranging and combining the various terms yields,

F (r, t) = (\frac{\partial}{\partial x}) (P_{x} E + M_{x} H) + (\frac{\partial}{\partial y}) (P_{y} E + M_{y} H) + (\frac{\partial}{\partial z}) (P_{z} E + M_{z} H)

+ ε_{o} [(\nabla \cdot E) E + (\nabla \times E) \times E] + μ_{o} [(\nabla \cdot H) H + (\nabla \times H) \times H] - ε_{o} μ_{o} \partial \frac{(E \times H)}{\partial t} .

(7)

This equation can be further expanded and rearranged to yield,

F (r, t) + (\frac{\partial}{\partial t}) (E \times \frac{H}{c^{2}}) = (\frac{\partial}{\partial x}) {P_{x} E + M_{x} H + ε_{o} [\frac{1}{2} (E_{x}^{2} - E_{y}^{2} - E_{z}^{2}) \hat{x} + E_{x} E_{y} \hat{y} + E_{x} E_{z} \hat{z}]

+ μ_{o} [\frac{1}{2} (H_{x}^{2} - H_{y}^{2} - H_{z}^{2}) \hat{x} + H_{x} H_{y} \hat{y} + H_{x} H_{z} \hat{z}]}

+ (\frac{\partial}{\partial y}) {P_{y} E + M_{y} H + ε_{o} [E_{x} E_{y} \hat{x} + \frac{1}{2} + (E_{y}^{2} - E_{x}^{2} - E_{z}^{2}) \hat{y} + E_{y} E_{z} \hat{z}]

+ μ_{o} [H_{x} H_{z} \hat{x} + \frac{1}{2} (H_{y}^{2} - H_{x}^{2} - H_{z}^{2}) \hat{y} + H_{y} H_{z} \hat{z}]}

+ (\frac{\partial}{\partial z}) {P_{z} E + M_{z} H + ε_{o} [E_{x} E_{z} \hat{x} + E_{y} E_{z} \hat{y} + \frac{1}{2} + (E_{z}^{2} - E_{x}^{2} - E_{y}^{2}) \hat{z}]

+ μ_{o} [H_{x} H_{z} \hat{x} + H_{y} H_{z} \hat{y} + \frac{1}{2} (H_{z}^{2} - H_{x}^{2} - H_{y}^{2}) \hat{z}]} .

(8)

Fig. 1. A small cube of dimensions Δx×Δy×Δz within a magnetic dielectric is separated from the surrounding medium by a fictitious vacuum-filled gap; the medium is specified by its (ε, µ) parameters. Assuming the gap is sufficiently narrow (compared to the wavelength of the electromagnetic field), its presence should not affect the distribution of the fields throughout the medium. Within the gap, however, the various components of the electromagnetic field are determined by the standard boundary conditions derived from Maxwell’s equations. In general, the tangential components of E and H remain continuous across the gap, while, in the perpendicular direction, the components of D and B retain continuity.

Next, we integrate Eq. (8) over the small Δx×Δy×Δz cube depicted in Fig. 1, normalize the resultant by the cube’s volume, and consider the limit when (Δx, Δy, Δz)→0. The left-hand side of Eq. (8) thus remains intact, but several changes occur on the right-hand side. For instance, in the first term, integration over x yields the argument of ∂/∂x, evaluated in the gaps on the left- and right-hand sides of the cube, then subtracted from each other. In these gaps, P_x =0, M_x =0, ε_o E_x =D_x , and µ_o H_x =B_x , while the remaining components of E and H retain the values that they have in the adjacent material environment. (These gap fields are found by invoking standard boundary conditions, namely, the continuity of tangential E and H , as well as perpendicular D and B components.) Similar arguments apply to the second and third terms on the right-hand side of Eq. (8), provided that, in the case of the 2^nd (3^rd) term, the initial integration is carried over y (z). When the integrals are fully evaluated and the result is normalized by the volume of the cube, we find, in the limit of a vanishing cube,

F (r, t) + (\frac{\partial}{\partial t}) (E \times \frac{H}{c^{2}}) = (\frac{\partial}{\partial x}) {[\frac{1}{2} ({ε_{o}^{- 1} D}_{x}^{2} - {ε_{o} E}_{y}^{2} - {ε_{o} E}_{z}^{2}) \hat{x} + D_{x} E_{y} \hat{y} + D_{x} E_{z} \hat{z}]

+ [\frac{1}{2} ({μ_{o}^{- 1} B}_{x}^{2} - {μ_{o} H}_{y}^{2} - {μ_{o} H}_{z}^{2}) \hat{x} + B_{x} H_{y} \hat{y} + B_{x} H_{z} \hat{z}]}

+ \frac{\partial}{\partial y} {[E_{x} D_{z} \hat{x} + \frac{1}{2} (ε_{o}^{- 1} D_{y}^{2} - ε_{o} E_{x}^{2} - ε_{o} E_{z}^{2}) \hat{y} + D_{y} E_{z} \hat{z}] + [H_{x} B_{y} \hat{x} + \frac{1}{2} (μ_{o}^{- 1} B_{y}^{2} - μ_{o} H_{x}^{2} - μ_{o} H_{z}^{2}) \hat{y} + B_{y} H_{z} \hat{z}]}

+ \frac{\partial}{\partial z} {[E_{x} D_{z} \hat{x} + E_{y} D_{z} \hat{y} + \frac{1}{2} (ε_{o}^{- 1} D_{z}^{2} - ε_{o} E_{x}^{2} - ε_{o} E_{y}^{2}) \hat{z}] + [H_{x} B_{z} \hat{x} + H_{y} B_{z} \hat{y} + \frac{1}{2} (μ_{o}^{- 1} B_{z}^{2} - μ_{o} H_{x}^{2} - μ_{o} H_{y}^{2}) \hat{z}]} .

(9)

Equation (9) clearly identifies the Abraham momentum density E × H /c ² as the electromagnetic momentum density G ( r , t), and yields the following stress tensor T_ij (i.e., rate of flow of momentum per unit area per unit time) within the medium:

T_{xx} = \frac{1}{2} (ε_{o} E_{y}^{2} + ε_{o} E_{z}^{2} - ε_{o}^{- 1} D_{x}^{2}) + \frac{1}{2} (μ_{o} H_{y}^{2} + μ_{o} H_{z}^{2} - μ_{o}^{- 1} B_{x}^{2}),

(10a)

T_{yx} = - D_{x} E_{y} - B_{x} H_{y},

(10b)

T_{zx} = - D_{x} E_{z} - B_{x} H_{z},

(10c)

T_{xy} = - D_{x} E_{y} - B_{x} H_{y},

(10d)

T_{yy} = \frac{1}{2} (ε_{o} E_{x}^{2} + ε_{o} E_{z}^{2} - ε_{o}^{- 1} D_{y}^{2}) + \frac{1}{2} (μ_{o} H_{x}^{2} + μ_{o} H_{z}^{2} - μ_{o}^{- 1} B_{y}^{2}),

(10e)

T_{zy} = - D_{y} E_{z} - B_{y} H_{z},

(10f)

T_{xz} = - E_{x} D_{z} - H_{x} B_{z},

(10g)

T_{yz} = - E_{y} D_{z} - H_{y} B_{z},

(10h)

T_{zz} = \frac{1}{2} (ε_{o} E_{x}^{2} + ε_{o} E_{y}^{2} - ε_{o}^{- 1} D_{z}^{2}) + \frac{1}{2} (μ_{o} H_{x}^{2} + μ_{o} H_{y}^{2} - μ_{o}^{- 1} B_{z}^{2}) .

(10i)

Equation (9) may thus be written as the following streamlined expression of momentum conservation:

\nabla \cdot T + F (r, t) + \frac{\partial G (r, t)}{\partial t} = 0 .

(11)

In its specific combination of the various components of the E , D , H and B fields, the stress tensor of Eq. (10) differs from both Abraham and Minkowski tensors. A similar (although by no means identical) tensor has been derived by Yaghjian [6

A. D. Yaghjian, “Internal energy, Q-energy, Poynting’s theorem, and the stress dyadic in dispersive material,” IEEE Trans. Anten. Prop. 55, 1495–1505 (2007). [CrossRef]

], who has advocated methods of analysis that in many respects resemble our methods.

Example 1. A plane electromagnetic wave propagates along the z-axis inside a medium specified by its (ε, µ) parameters. The linearly polarized plane-wave has E-field amplitude E _o x̂ and H-field amplitude

H_{o} \hat{y} = Z_{o}^{- 1} \sqrt{\frac{ε}{μ}} E_{o} \hat{y}

, where

Z_{o} = \sqrt{\frac{μ_{o}}{ε_{o}}}

is the impedance of the free space. Assuming a monochromatic plane wave with angular frequency ω, the rate of flow of momentum (per unit area per unit time) along the z-axis will be given by Eq. (10i) as follows:

T_{zz} = \frac{1}{2} (ε_{o} E_{x}^{2} + μ_{o} H_{y}^{2}) + \frac{1}{2} ε_{o} [1 + (\frac{ε}{μ})] E_{o}^{2} \cos^{2} (ω t) .

(12)

The rate of flow of energy (per unit area per unit time) is the Poynting vector, S = E × H . Therefore,

S_{z} = E_{x} H_{y} = Z_{o}^{- 1} \sqrt{\frac{ε}{μ}} E_{o}^{2} \cos^{2} (ω t) .

(13)

Suppose a total of N photons cross the xy-plane at z=0 during the time interval [0, τ]. Since each photon has energy hf, we have <S_z >τ=Nhf (angled brackets denote time-averaging). The total momentum crossing the same plane during the same time interval will therefore be

< T_{zz} > τ = \frac{1}{2} N (\sqrt{\frac{ε}{μ}} + \sqrt{\frac{μ}{ε}}) \frac{hf}{c}

. A single photon’s momentum is thus

\frac{1}{2} (\sqrt{\frac{ε}{μ}} + \sqrt{\frac{μ}{ε}}) \frac{hf}{c}

, in agreement with our previous results [3

M. Mansuripur, “Radiation pressure and the linear momentum of the electro-magnetic field in magnetic media,” Opt. Express 15, 13502–13518 (2007). [CrossRef] [PubMed]

]. This photon momentum, which in non-magnetic dielectrics is equal to the arithmetic average of the Minkowski and Abraham momenta, is always greater than the photon momentum in free space, hf/c. In general, the photon momentum consists of an electromagnetic part and a mechanical part. In a non-dispersive medium where the group velocity of light equals its phase velocity, the electromagnetic momentum of a single photon is

\frac{hf}{(\sqrt{ε μ} c)}

under these circumstances, the photon’s mechanical momentum will be

[\frac{(ε + μ - 2)}{(2 \sqrt{ε μ})}] \frac{hf}{c}

Example 2. With reference to Fig. 2, consider a collimated, monochromatic beam of light propagating along the z-axis within a linear, isotropic, and homogeneous medium specified by its (ε,µ) parameters. The beam, which has a finite-diameter along the x-axis and an infinite diameter along y, is linearly polarized, having E-field amplitude E _o x̂ and H-field amplitude

H_{o} \hat{y} = Z_{o}^{- 1} \sqrt{\frac{ε}{μ}} E_{o} \hat{y}

at the center. In addition, there exists a weak E_z component of the field, which is an odd function of x and goes to zero at the center. At the central yz-plane, the time-averaged rate of flow of x-momentum along the x-axis is given by Eq. (10a), as follows:

< T_{xx} > = < - \frac{1}{2} ε_{o}^{- 1} D_{x}^{2} + \frac{1}{2} μ_{o} H_{y}^{2} > = \frac{1}{4} ε_{o} [(\frac{ε}{μ}) - ε^{2}] E_{o}^{2} .

(14)

Fig. 2. A collimated beam of light having a finite diameter along the x-axis, propagates along z in a homogeneous medium specified by its (ε,µ) parameters. The beam has transverse magnetic (TM) or p-polarization, that is, its electromagnetic field components are (E_x , E_z , H_y ). A narrow gap opened in the central region of the beam reveals the existence of a force on the adjacent layers of dipoles. Continuity of D _⊥ yields the E-field within the gap as ε E _o x̂ . The E-field acting on the negative charges of the upper layer of the dipoles (as well as that acting on the positive charges of the lower dipoles) is ½(ε+1)E _o x̂ , whereas the field acting on the positive charges of the upper dipoles (or negative charges of the lower dipoles) is E _o x̂ . These boundary dipole layers, therefore, experience an E-field gradient proportional to ½(ε-1)E _o x̂ . The net force of the E-field gradient exerted on the upper boundary layer is downward, while that on the lower boundary layer is upward. The two forces, being equal in magnitude, cancel each other out, but each force must be taken into account when considering the total force on the upper or lower halves of the medium. In addition to forces at the boundary layers, the sidewalls of the beam exert a force on the medium as well; the density of this force (per unit area of the sidewall) is denoted by F _sw .

This momentum flow is entirely converted to a force on the electric dipoles located just above the z-axis and a second force, F _sw , exerted on the medium by the upper sidewall of the beam. (The system being symmetric with respect to the yz-plane, identical forces, albeit in opposite directions, act on the lower half of the medium.) The force on the dipoles immediately above the z-axis is best understood if one introduces a gap in the middle of the beam as indicated in Fig. 2. The continuity of D _⊥ at this interface reveals the E-field within the gap as being equal to ε E _o x̂ . The average E-field at the interface is thus ½(ε+1)E _o, and the field gradient sensed by the interfacial dipole layer is proportional to ½(ε-1)E _o. The dipole density being P =ε _o(ε-1)E _o x̂ , we find a force density at the interface given by <F_x >=¼ε _o(ε-1)² E _o ². Adding this force density to <T_xx > of Eq. (14) yields <F_x ^(sw)>=¼ε _o[(ε/µ)-2ε+1]E _o ², which is consistent with the sidewall force density of finite-diameter beams found in [3

M. Mansuripur, “Radiation pressure and the linear momentum of the electro-magnetic field in magnetic media,” Opt. Express 15, 13502–13518 (2007). [CrossRef] [PubMed]

Example 3. Figure 3 shows a collimated, monochromatic beam of finite width propagating in a homogeneous medium specified by its (ε, µ) parameters. The propagation direction makes an angle θ with the z-axis in the xz-plane. The stress tensor of Eq. (10) gives the following time-averaged rate of flow of momentum (per unit area per unit time) across the xy-plane:

< T_{xz} \hat{x} + T_{zz} \hat{z} > = - < E_{x} D_{z} > \hat{x} + \frac{1}{2} < ε_{o} E_{x}^{2} - ε_{o}^{- 1} D_{z}^{2} + μ_{o} H_{y}^{2} > \hat{z}

= \frac{1}{2} ε_{o} E_{o}^{2} {ε \cos θ \sin θ \hat{x} + \frac{1}{2} [\cos^{2} θ - ε^{2} \sin^{2} θ + (\frac{ε}{μ})] \hat{z}} .

(15)

From Eq. (12), the rate of flow of momentum in a beam of cross-section cosθ should be

\frac{d < p >}{d t} = \frac{1}{4} ε_{o} [1 + (\frac{ε}{μ})] E_{o}^{2} \cos θ (\sin θ \hat{x} + \cos θ \hat{z}) .

(16)

From Example 2, we find the imbalance of the sidewall force on the lower sidewall depicted in Fig. 3 to be

< F_{sw} > = \frac{1}{4} ε_{o} [(\frac{ε}{μ}) - 2 ε + 1] E_{o}^{2} \sin θ (- \cos θ \hat{x} + \sin θ \hat{z}) .

(17)

A narrow gap, opened parallel to the x-axis, reveals the force exerted on the boundary electric dipoles due to the E_z discontinuity. The strength of this dipole layer is P =ε _o(ε-1)E _o sinθ $\hat{z}$ , and the effective E-field gradient acting on it is proportional to ½(ε-1)E _osinθ $\hat{z}$ . The effective force on the dipole layer (per unit area) is thus given by

< F > = - \frac{1}{4} ε_{o} {(ε - 1)}^{2} E_{o}^{2} \sin^{2} θ \hat{z} .

(18)

The combined forces in Eqs. (16–18) are identical with the momentum flow rate of Eq. (15). The stress tensor of Eq. (10) is thus seen to yield the correct rate of momentum flow along the propagation direction when the relevant boundary forces are properly taken into account.

Fig. 3. A collimated beam of finite diameter in the x-direction (and infinite diameter along y) propagates in a medium specified by its (ε, µ) parameters. The propagation direction makes an angle θ with the z-axis in the xz-plane. The beam’s foot-print on the x-axis has unit length, making the beam width equal to cosθ, as shown. The beam is transverse magnetic (TM) or p-polarized, that is, its electromagnetic field components are (E_x , E_z , H_y ). A narrow gap, opened parallel to the x-axis, reveals the force exerted on the boundary layer electric dipoles due to E_z discontinuity. The effective E-field gradient acting on the boundary dipole layer is proportional to ½(ε-1)E _osinθ $\hat{z}$ . There is also an imbalance between the forces acting at the beam’s upper and lower sidewalls, due to the extra length sinθ of the lower wall. When the force on the boundary dipole layer as well as the imbalance of the sidewall forces are taken into account, the stress tensor component T_xz x̂ +T_zz ẑ yields the rate of flow of momentum along the propagation direction.

3. Concluding remarks

We have derived the electromagnetic stress tensor of Eq. (10) for ponderable media in terms of the macroscopic fields E , D , H , and B . When the forces acting at the boundaries of a specified region of the medium are properly taken into account, we have shown that the stress tensor, in accordance with Eq. (11), yields the correct rate of flow of momentum, including both electromagnetic and mechanical momenta.

References and links

1.	J. D. Jackson, Classical Electrodynamics , 2^nd edition, Wiley, New York, 1975.
2.	R. N. C. Pfeifer, T. A. Nieminen, N. R Heckenberg, and H. Rubinsztein-Dunlop, “Momentum of an electromagnetic wave in dielectric media,” Rev. Mod. Phys. 79, 1197–1216 (2007). [CrossRef]
3.	M. Mansuripur, “Radiation pressure and the linear momentum of the electro-magnetic field in magnetic media,” Opt. Express 15, 13502–13518 (2007). [CrossRef] [PubMed]
4.	T. B. Hansen and A. D. Yaghjian, Plane-Wave Theory of Time-Domain Fields: Near-Field Scanning Applications , IEEE Press, New York (1999). [CrossRef]
5.	B. A. Kemp, J. A. Kong, and T. Grzegorczyk, “Reversal of wave momentum in isotropic left-handed media,” Phys. Rev. A 75, 053810 (2007). [CrossRef]
6.	A. D. Yaghjian, “Internal energy, Q-energy, Poynting’s theorem, and the stress dyadic in dispersive material,” IEEE Trans. Anten. Prop. 55, 1495–1505 (2007). [CrossRef]

OCIS Codes
(140.7010) Lasers and laser optics : Laser trapping
(260.2110) Physical optics : Electromagnetic optics

ToC Category:
Physical Optics

History
Original Manuscript: November 27, 2007
Revised Manuscript: March 25, 2008
Manuscript Accepted: March 27, 2008
Published: April 1, 2008

Citation
Masud Mansuripur, "Electromagnetic stress tensor in ponderable media," Opt. Express 16, 5193-5198 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-8-5193

Sort: Year | Journal | Reset

References

J. D. Jackson, Classical Electrodynamics, 2nd edition, Wiley, New York, 1975.
R. N. C. Pfeifer, T. A. Nieminen, N. R Heckenberg, and H. Rubinsztein-Dunlop, "Momentum of an electro-magnetic wave in dielectric media," Rev. Mod. Phys. 79, 1197-1216 (2007). [CrossRef]
M. Mansuripur, "Radiation pressure and the linear momentum of the electromagnetic field in magnetic media," Opt. Express 15, 13502-13518 (2007). [CrossRef] [PubMed]
T. B. Hansen and A. D. Yaghjian, Plane-Wave Theory of Time-Domain Fields: Near-Field Scanning Applications, IEEE Press, New York (1999). [CrossRef]
B. A. Kemp, J. A. Kong, and T. Grzegorczyk, "Reversal of wave momentum in isotropic left-handed media," Phys. Rev. A 75, 053810 (2007). [CrossRef]
A. D. Yaghjian, "Internal energy, Q-energy, Poynting’s theorem, and the stress dyadic in dispersive material," IEEE Trans. Anten. Prop. 55, 1495-1505 (2007). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper