Radiation pressure and the linear momentum of the electromagnetic field Masud mansuripur Optical Sciences Center. The University of Arona, Tucson, Arisona 85721 Abstract: We derive the force of the electromagnetic radiation on material objects by a direct application of the Lorentz law of classical electro- dynamics. The derivation is straightforward in the case of solid metals and solid dielectrics, where the mass density and the optical constants of the media are assumed to remain unchanged under internal and external pressures, and where material flow and deformation can be ignored. For metallic mirrors, we separate the contribution to the radiation pressure of the electrical charge density from that of the current density of the conduction electrons. In the case of dielectric media. we examine the forces experienced by bound charges and currents, and determine the contribution of each to the radiation pressure. These analyses reveal the existence of a lateral radiation pressure inside the dielectric media, one that is exerted at and around the edges of a finite-diameter light beam. The lateral pressure turns out to be compressive for s-polarized light and expansive for p- polarized light. Along the way, we derive an expression for the momentum density of the light field inside dielectric media, one that has equal contributions from the traditional minkowski and abraham forms this ew expression for the momentum density, which contains both electromagnetic and mechanical terms, is used to explain the behavior of light pulses and individual photons upon entering and exiting a dielectric slab. In all the cases considered, the net forces and torques experienced by material bodies are consistent with the relevant conservation laws. Our method of calculating the radiation pressure can be used in conjunction with numerical simulations to yield the distribution of fields and forces in diverse systems of practical interest C 2004 Optical Society of America OCIS codes:(2602110)Electromagnetic theory; (1407010)Trapping;(0207010) Trapping, (3106860)Thin films, optical properties References 1. H Minkowski, Nachr. Ges. Wiss. Gottingen 53(1908) 2. H Minkowski, Math. Annalon 68, 472(1910) 3. M. Abraham, R C Circ. Mat. Palermo 28, 1(1909) 4. M. Abraham, R C. Circ. Mat. Palermo 30, 33(1910) 5. J. P. Gordon, "Radiation forces and momenta in dielectric media, "Phys. Rev. A8, 14-21(1973) 6. R Loudon. "Radiation Pressure and Momentum in Dielectrics. "De Martini lecture. in Fortschritte der Physik(2004) 7. R. Loudon, "Theory of the radiati re on dielectric surfaces, J Mod. Opt. 49, 821-838(2002 8. L.Landau, E. Lifshitz, Electrodynamics of Continous Media, Pergamon, New York, 1960 9. J.D. Jackson, Classical Electrodynamics, 2 edition, wiley, New York, 1975 10. M. Planck, The Theory of Heat Radiation, translated by M. Masius form the German edition of 1914, Dover Publications, New York(1959) I. R. V Jones and J. C S. Richards, Proc. Roy. Soc. A 221, 480(1954) #5025-S1500US Received 10 August 2004; revised 13 October 2004; accepted 20 October 2004 (C)2004OSA November 2004/Vol 12. No 22/OPTICS EXPRESS 5375
Radiation pressure and the linear momentum of the electromagnetic field Masud Mansuripur Optical Sciences Center, The University of Arizona, Tucson, Arizona 85721 masud@u.arizona.edu Abstract: We derive the force of the electromagnetic radiation on material objects by a direct application of the Lorentz law of classical electrodynamics. The derivation is straightforward in the case of solid metals and solid dielectrics, where the mass density and the optical constants of the media are assumed to remain unchanged under internal and external pressures, and where material flow and deformation can be ignored. For metallic mirrors, we separate the contribution to the radiation pressure of the electrical charge density from that of the current density of the conduction electrons. In the case of dielectric media, we examine the forces experienced by bound charges and currents, and determine the contribution of each to the radiation pressure. These analyses reveal the existence of a lateral radiation pressure inside the dielectric media, one that is exerted at and around the edges of a finite-diameter light beam. The lateral pressure turns out to be compressive for s-polarized light and expansive for ppolarized light. Along the way, we derive an expression for the momentum density of the light field inside dielectric media, one that has equal contributions from the traditional Minkowski and Abraham forms. This new expression for the momentum density, which contains both electromagnetic and mechanical terms, is used to explain the behavior of light pulses and individual photons upon entering and exiting a dielectric slab. In all the cases considered, the net forces and torques experienced by material bodies are consistent with the relevant conservation laws. Our method of calculating the radiation pressure can be used in conjunction with numerical simulations to yield the distribution of fields and forces in diverse systems of practical interest. © 2004 Optical Society of America OCIS codes: (260.2110) Electromagnetic theory; (140.7010) Trapping; (020.7010) Trapping; (310.6860) Thin films, optical properties References 1. H. Minkowski, Nachr. Ges. Wiss. Gottingen 53 (1908). 2. H. Minkowski, Math. Annalon 68, 472 (1910). 3. M. Abraham, R. C. Circ. Mat. Palermo 28, 1 (1909). 4. M. Abraham, R. C. Circ. Mat. Palermo 30, 33 (1910). 5. J. P. Gordon, “Radiation forces and momenta in dielectric media,” Phys. Rev. A 8, 14-21 (1973). 6. R. Loudon, “Radiation Pressure and Momentum in Dielectrics,” De Martini lecture, to appear in Fortschritte der Physik (2004). 7. R. Loudon, “Theory of the radiation pressure on dielectric surfaces,” J. Mod. Opt. 49, 821-838 (2002). 8. L. Landau, E. Lifshitz, Electrodynamics of Continuous Media, Pergamon, New York, 1960. 9. J. D. Jackson, Classical Electrodynamics, 2nd edition, Wiley, New York, 1975. 10. M. Planck, The Theory of Heat Radiation, translated by M. Masius form the German edition of 1914, Dover Publications, New York (1959). 11. R. V. Jones and J. C. S. Richards, Proc. Roy. Soc. A 221, 480 (1954). (C) 2004 OSA 1 November 2004 / Vol. 12, No. 22 / OPTICS EXPRESS 5375 #5025- $15.00 US Received 10 August 2004; revised 13 October 2004; accepted 20 October 2004
12. A. Ashkin, J. M. Dziedzic, J. E Bjorkholm and S. Chu, "Observation of a single-beam gradient force optical trap for dielectric particles, Opt. Lett. 11, 288-290(1986) 13. A. Ashkin and J M. Dziedzic, "Optical trapping and manipulation of viruses and bacteria, "Science 235,1517. 1520(1987) 14. A Rohrbach and E Stelzer, "Trapping force onstants, and potential depths for dielectric sphere presence of spherical aberrations, " AppL. Opt. 41, 2494(2002) 15. Y N Obukhov and F WHehl, "Electromagnetic energy-momentum and forces in matter, "Phys.Lett. A,311 16. G.Barlow, Proc. Roy. Soc. Lond. A 87, 1-16(1912) 17. R. Loudon, "Theory of the forces exerted by Laguerre-Gaussian light beams on dielectrics, "Phys. Rev. A68 013806(2003) 1. Introduction It is well-known that the electromagnetic radiation carries both energy and momentum, that the energy flux is given by the Poynting vector S of the classical electrodynamics, and that, in e-space, the momentum density p(i. e, momentum per unit volume)is given by p=S/c where c is the speed of light in vacuum. What has been a matter of controversy for quite some time now is the proper form for the momentum of the electromagnetic waves in dielectric media. The question is whether the momentum density in a material medium has the form P=DXB, due to Minkowski [1, 2], or p=Ex H/c2, due to Abraham [3, 4). J. P. Gordon [51 attributes the following comment to E. I Blount: "The argument has not, it is true, been carried out at high volume, but the list of disputants is very distinguished. For a historical perspective on the subject and a summary of the relevant experimental results see R. Loudon 6,7] and J. P. Gordon [5 Traditionally, the electromagnetic stress tensor has been used to derive the mechanical force exerted by the radiation field on ponderable media [8, 9]. This approach, while having the advantage of generality, tends to obscure behind complicated mathematics the physical origin of the forces. It is possible, however, to calculate the force of the electromagnetic radiation on various media by direct invocation of the Lorentz law of force. The derivation is especially straightforward in the case of solid metals and solid dielectrics, where the mass tensity and the optical constants of the media may be assumed to remain constant under internal and external pressures, and where material flow and deformation can be ignored Loudon [7] has emphasized "the simplicity and safety of calculations based on the Lorentz force and the dangers of calculations based on derived expressions involving elements of the Maxwell stress tensor, whose contributions may vanish in some situations but not in others In this paper we use the Lorentz law to derive the force of electromagnetic radiation on isotropic solid media in several simple situations. In the case of metallic mirrors, we separate following Planck [101 the contribution to the radiation pressure of the electrical charge tensity from that of the current density(both due to conduction electrons ). In the case of dielectric media, we examine the force experienced by bound charges and currents, an determine the contribution of each to the radiation pressure. Along the way, we derive a new expression for the momentum density of the light field inside dielectric media, one that has equal contributions from the aforementioned Minkowski and Abraham forms. This new expression for the momentum density, which contains both mechanical and electromagnetic terms, is subsequently used to elucidate the behavior of individual photons upon entering and exiting a dielectric slab. With the exception of the semi-quantitative results of Section 13, all the results obtained in this paper are exact, in the sense that no approximations or implifications have been introduced; all derivations are based directly on the Lorentz law of force in conjunction with Maxwells equations, using the standard constitutive relations for homogeneous, isotropic, linear, non-magnetic, and non-dispersive media The organization of the paper is as follows. In Section 2 we describe the notation and define the various parameters used throughout the paper. Section 3 considers the reflection of #5025-S1500US Received 10 August 2004; revised 13 October 2004; accepted 20 October 2004 (C)2004OSA November 2004/Vol 12. No 22/OPTICS EXPRESS 5376
12. A. Ashkin, J. M. Dziedzic, J. E. Bjorkholm and S. Chu, “Observation of a single-beam gradient force optical trap for dielectric particles,” Opt. Lett. 11, 288-290 (1986). 13. A. Ashkin and J. M. Dziedzic, “Optical trapping and manipulation of viruses and bacteria,” Science 235, 1517- 1520 (1987). 14. A. Rohrbach and E. Stelzer, “Trapping forces, force constants, and potential depths for dielectric spheres in the presence of spherical aberrations,” Appl. Opt. 41, 2494 (2002). 15. Y. N. Obukhov and F. W.Hehl, “Electromagnetic energy-momentum and forces in matter,” Phys. Lett. A, 311, 277-284 (2003). 16. G. Barlow, Proc. Roy. Soc. Lond. A 87, 1-16 (1912). 17. R. Loudon, “Theory of the forces exerted by Laguerre-Gaussian light beams on dielectrics,” Phys. Rev. A 68, 013806 (2003). 1. Introduction It is well-known that the electromagnetic radiation carries both energy and momentum, that the energy flux is given by the Poynting vector S of the classical electrodynamics, and that, in free-space, the momentum density p (i.e., momentum per unit volume) is given by p = S/c 2 , where c is the speed of light in vacuum. What has been a matter of controversy for quite some time now is the proper form for the momentum of the electromagnetic waves in dielectric media. The question is whether the momentum density in a material medium has the form p = D × B, due to Minkowski [1,2], or p = E × H/c 2 , due to Abraham [3,4]. J. P. Gordon [5] attributes the following comment to E. I. Blount: “The argument has not, it is true, been carried out at high volume, but the list of disputants is very distinguished.” For a historical perspective on the subject and a summary of the relevant experimental results see R. Loudon [6,7] and J. P. Gordon [5]. Traditionally, the electromagnetic stress tensor has been used to derive the mechanical force exerted by the radiation field on ponderable media [8,9]. This approach, while having the advantage of generality, tends to obscure behind complicated mathematics the physical origin of the forces. It is possible, however, to calculate the force of the electromagnetic radiation on various media by direct invocation of the Lorentz law of force. The derivation is especially straightforward in the case of solid metals and solid dielectrics, where the mass density and the optical constants of the media may be assumed to remain constant under internal and external pressures, and where material flow and deformation can be ignored. Loudon [7] has emphasized “the simplicity and safety of calculations based on the Lorentz force and the dangers of calculations based on derived expressions involving elements of the Maxwell stress tensor, whose contributions may vanish in some situations but not in others.” In this paper we use the Lorentz law to derive the force of electromagnetic radiation on isotropic solid media in several simple situations. In the case of metallic mirrors, we separate, following Planck [10], the contribution to the radiation pressure of the electrical charge density from that of the current density (both due to conduction electrons). In the case of dielectric media, we examine the force experienced by bound charges and currents, and determine the contribution of each to the radiation pressure. Along the way, we derive a new expression for the momentum density of the light field inside dielectric media, one that has equal contributions from the aforementioned Minkowski and Abraham forms. This new expression for the momentum density, which contains both mechanical and electromagnetic terms, is subsequently used to elucidate the behavior of individual photons upon entering and exiting a dielectric slab. With the exception of the semi-quantitative results of Section 13, all the results obtained in this paper are exact, in the sense that no approximations or simplifications have been introduced; all derivations are based directly on the Lorentz law of force in conjunction with Maxwell’s equations, using the standard constitutive relations for homogeneous, isotropic, linear, non-magnetic, and non-dispersive media. The organization of the paper is as follows. In Section 2 we describe the notation and define the various parameters used throughout the paper. Section 3 considers the reflection of (C) 2004 OSA 1 November 2004 / Vol. 12, No. 22 / OPTICS EXPRESS 5376 #5025- $15.00 US Received 10 August 2004; revised 13 October 2004; accepted 20 October 2004
a plane electromagnetic wave from a perfect conductor, and shows that both the Lorentz law of force and the momentum of the electromagnetic field in free-space can consistently account for the radiation pressure exerted on the mirror surface. In Section 4 we use the Lorentz law to determine the radiation pressure on the surface of a semi-infinite dielectric medium at normal incidence. Here we derive an expression for the momentum of the field inside the dielectrics. The results of Section 4 are then extended to cover the case of oblique incidence on a semi-infinite dielectric, first with s-polarized light in Section 5, then with p- polarized light in Section 6. These analyses lead to the discovery of a lateral radiation pressure inside the dielectric medium, exerted at and around the edges of a finite-diameter lane wave. This lateral presond expansive for p-light. To the author's best knowledge, the ressure. while having the same magnitude in both cases turns out to be compressive for s-light expansive lateral force on the dielectric host of p-polarized light has not been discussed in the existing literature, making it a novel prediction that requires experimental verification In Section 7 we examine the torque experienced by a dielectric slab, illuminated at the Brewsters angle by a p-polarized plane wave. The torque is calculated directly from the Lorentz law applied to the induced(bound) charges at the surfaces of the slab, then shown to be consistent with the change in the angular momentum of the incident light. The case of an anti-reflection coated, semi-infinite dielectric medium is taken up in Section 8, where the increase in the momentum of the incident beam upon transmission into the dielectric medium is shown to result in a net force on the anti-reflection coating layer that tends to peel the layer away from its substrate, another new result that requires experimental verification. In Section 9 we analyze the case of a dielectric slab of finite thickness, and show that optical interference within the slab is responsible for the (longitudinal) stress induced by the electromagnetic radiation For a different perspective on the lateral pressure at the edges of a finite-diameter beam a dielectric, Section 10 is devoted to an analysis of the one-dimensional Gaussian beam inside a dielectric medium. Depending on the direction of the E-field, we show that the lateral pressure on the medium can be compressive or expansive, and that the magnitude and direction of this radiation force are in complete accord with the results of Sections 5 and 6 The generality of this lateral pressure(and the dependence of its direction on the state of polarization) are brought to the fore in Section 1l, where the simple fringes produced by the interference between two plane waves are shown to exhibit the same phenomena In Section 12 we extend our method of calculation of force and momentum to(finite- duration) light pulses, where the leading and trailing edges of the pulse are shown to play an important role in exchanging the electromagnetic momentum of the light with the mechanical momentum of the medium of propagation. The physical basis for the designation(in Section 4)of a fraction of the photon's momentum as"mechanical"is clarified in Section 12 The classical experiments pertaining to metallic mirrors immersed in liquid dielectrics [11] are discussed in Section 13, where they are shown to be in complete agreement with our theoretical calculations. It is well known, e.g, from optical tweezers experiments[12-14], that a focused laser beam tends to attract small dielectric beads toward the center of the focused spot. At first glance, this observation might seem at odds with the presence of an expansive lateral pressure inside the dielectric medium of the bead. To resolve this apparent discrepancy, Section 14 offers a semi-quantitative analysis of a simplified model of the optical tweezers experiment. Final remarks and a summary of our important results appear in Section 15 The MKSA system of units is used throughout the paper. Time harmonic fields are written as E(x,,=,0)=E(x,y, =)exp(@r), where @=2f is the angular frequency. For brevity, we omit the explicit dependence of the fields on x, y, = t. To specify their magnitude and phase, #5025-S1500US Received 10 August 2004; revised 13 October 2004; accepted 20 October 2004 (C)2004OSA November 2004/Vol 12. No 22/OPTICS EXPRESS 5377
a plane electromagnetic wave from a perfect conductor, and shows that both the Lorentz law of force and the momentum of the electromagnetic field in free-space can consistently account for the radiation pressure exerted on the mirror surface. In Section 4 we use the Lorentz law to determine the radiation pressure on the surface of a semi-infinite dielectric medium at normal incidence. Here we derive an expression for the momentum of the field inside the dielectrics. The results of Section 4 are then extended to cover the case of oblique incidence on a semi-infinite dielectric, first with s-polarized light in Section 5, then with ppolarized light in Section 6. These analyses lead to the discovery of a lateral radiation pressure inside the dielectric medium, exerted at and around the edges of a finite-diameter plane wave. This lateral pressure, while having the same magnitude in both cases, turns out to be compressive for s-light and expansive for p-light. To the author’s best knowledge, the expansive lateral force on the dielectric host of p-polarized light has not been discussed in the existing literature, making it a novel prediction that requires experimental verification. In Section 7 we examine the torque experienced by a dielectric slab, illuminated at the Brewster’s angle by a p-polarized plane wave. The torque is calculated directly from the Lorentz law applied to the induced (bound) charges at the surfaces of the slab, then shown to be consistent with the change in the angular momentum of the incident light. The case of an anti-reflection coated, semi-infinite dielectric medium is taken up in Section 8, where the increase in the momentum of the incident beam upon transmission into the dielectric medium is shown to result in a net force on the anti-reflection coating layer that tends to peel the layer away from its substrate; another new result that requires experimental verification. In Section 9 we analyze the case of a dielectric slab of finite thickness, and show that optical interference within the slab is responsible for the (longitudinal) stress induced by the electromagnetic radiation. For a different perspective on the lateral pressure at the edges of a finite-diameter beam in a dielectric, Section 10 is devoted to an analysis of the one-dimensional Gaussian beam inside a dielectric medium. Depending on the direction of the E-field, we show that the lateral pressure on the medium can be compressive or expansive, and that the magnitude and direction of this radiation force are in complete accord with the results of Sections 5 and 6. The generality of this lateral pressure (and the dependence of its direction on the state of polarization) are brought to the fore in Section 11, where the simple fringes produced by the interference between two plane waves are shown to exhibit the same phenomena. In Section 12 we extend our method of calculation of force and momentum to (finiteduration) light pulses, where the leading and trailing edges of the pulse are shown to play an important role in exchanging the electromagnetic momentum of the light with the mechanical momentum of the medium of propagation. The physical basis for the designation (in Section 4) of a fraction of the photon’s momentum as “mechanical” is clarified in Section 12. The classical experiments pertaining to metallic mirrors immersed in liquid dielectrics [11] are discussed in Section 13, where they are shown to be in complete agreement with our theoretical calculations. It is well known, e.g., from optical tweezers experiments [12-14], that a focused laser beam tends to attract small dielectric beads toward the center of the focused spot. At first glance, this observation might seem at odds with the presence of an expansive lateral pressure inside the dielectric medium of the bead. To resolve this apparent discrepancy, Section 14 offers a semi-quantitative analysis of a simplified model of the optical tweezers experiment. Final remarks and a summary of our important results appear in Section 15. 2. Notation and basic definitions The MKSA system of units is used throughout the paper. Time harmonic fields are written as E (x, y, z, t) = E (x, y, z) exp(−iω t), where ω = 2πf is the angular frequency. For brevity, we omit the explicit dependence of the fields on x, y, z, t. To specify their magnitude and phase, (C) 2004 OSA 1 November 2004 / Vol. 12, No. 22 / OPTICS EXPRESS 5377 #5025- $15.00 US Received 10 August 2004; revised 13 October 2004; accepted 20 October 2004
complex amplitudes such as E are expressed as elexp(idE). Time-averaged products of two fields, say, Real (A)=| I cos(@t-oA) and Real(B)=B I cos(@t-oB), given by v2) AB I cos(OA-oB), may also be written /Real(AB") To compute the force exerted by the electromagnetic field on a given medium, we use Maxwell's equations to determine the distributions of the E-and H-fields both inside and outside the medium. We then apply the lorentz law F=q(E+vx B), which gives the electromagnetic force on a particle of charge q and velocity V. The magnetic induction B is assumed to be related to the H-field via B=loH, where lo=4I x 10 henrys/meter is the permeability of free space, in other words, any magnetic moments that might exist in the medium and their interactions with the radiation field at optical frequencies are ignored Typically, there are no free charges in the system, so V D=0, where D=EE+ P is the electric displacement vector, E =8.82 x 10 farads/meter is the permittivity of free space, and P is the local polarization density within the medium. In linear media, D=Eoe E, where E is the mediums relative permittivity; hence, P=E(E-lE. We ignore the frequency- dependence of e throughout the paper and treat the media as non-dispersive When V D=0, the density of bound charges P,=-V P may be expressed as Pb=EV- E Inside a homogeneous and isotropic medium, E being proportional to D and V.D=0 imply that PB=0; no bound charges, therefore, can exist inside such media. However, at the interface between two different media, the component of d perpendicular to the interface, Di, must be continuous. The implication is that El is discontinuous and therefore, bound charges can exist at such interfaces; the interfacial bound charges will thus have an areal density o=E(E21- Eln. Under the influence of the local E-field, these charges give rise to an electric Lorentz force F=VReal(oE), where F is the force per unit area of there is no ambiguity as to which field must be used in conjunction with the lorentz law s the interface. Since the tangential E-field, El, is generally continuous across the interi for the perpendicular component, the average E across the boundary, 2(E11+ E21), must be used in calculating the interfacial force. (The use of the av erage El in this context is not a matter of choice; it is the only way to get the calculated force at the boundary to agree with the time rate of change of the momentum that passes through the interface. From a physical andpoint, the interfacial charges produce a local Ei that has the same magnitude but opposite directions on the two sides of the interface. It is this locally-generated Ei that is responsible for the E-fields discontinuity. Averaging El across the interface eliminates the local field, as it should, since the charge cannot exert a force on itself. Since V. B=0 and B=uoh, the perpendicular H-field, Hi, at the interface between adjacent media must remain continuous. The tangential H-field at such interfaces, however, may be discontinuous. This, in accordance with Maxwells equation VxH=J+ aD/a, gives rise to an interfacial current density J, =H21-HIl Such currents can exist on the surfaces of ood conductors, where Eyl is negligible, yet the high conductance of the medium permits the flow of the surface current. Elsewhere, the only source of electrical currents are bound charges, with the bound current density being Jp=ap/dt=E(e-1)dE/ar. Assuming time harmonic fields with the time-dependence factor exp(ian), we can write J=-ioE(E-lE The H-field of the electromagnetic wave then exerts a force on the bound current according to the Lorentz law, namely, F=Real ( bx B), where F is force per unit volume Note: For time-harmonic fields, the contribution of conduction electrons to current density may be combined with that of bound electrons. Since J=o E, where O is the conductivity of the medium, the net current density Jc Jb may be attributed to an effective dielectric constant E+i(oc/E@). In general, since E is complex-valued, there is no need to distinguish conduction electrons from bound electrons, and e may be treated as an effective dielectric constant that contains both contributions. An exception will be made in Section 3 in the case #5025-S1500US Received 10 August 2004; revised 13 October 2004; accepted 20 October 2004 (C)2004OSA November 2004/Vol 12. No 22/OPTICS EXPRESS 5378
complex amplitudes such as E are expressed as | E | exp(iφE). Time-averaged products of two fields, say, Real (A ) = | A | cos(ω t – φA) and Real (B ) = |B | cos(ω t – φB), given by ½| AB | cos(φA – φB), may also be written ½Real (AB*). To compute the force exerted by the electromagnetic field on a given medium, we use Maxwell’s equations to determine the distributions of the E- and H-fields both inside and outside the medium. We then apply the Lorentz law F = q (E + V × B), which gives the electromagnetic force on a particle of charge q and velocity V. The magnetic induction B is assumed to be related to the H-field via B = µ oH, where µ o = 4π × 10−7 henrys/meter is the permeability of free space; in other words, any magnetic moments that might exist in the medium and their interactions with the radiation field at optical frequencies are ignored. Typically, there are no free charges in the system, so ∇ · D = 0, where D = εoE + P is the electric displacement vector, εo = 8.82 × 10−12 farads/meter is the permittivity of free space, and P is the local polarization density within the medium. In linear media, D =εoε E, where ε is the medium’s relative permittivity; hence, P =εo(ε – 1)E. We ignore the frequencydependence of ε throughout the paper and treat the media as non-dispersive. When ∇ · D = 0, the density of bound charges ρb = −∇ · P may be expressed as ρb =εo∇ · E. Inside a homogeneous and isotropic medium, E being proportional to D and ∇ · D = 0 imply that ρb = 0; no bound charges, therefore, can exist inside such media. However, at the interface between two different media, the component of D perpendicular to the interface, D⊥, must be continuous. The implication is that E⊥ is discontinuous and, therefore, bound charges can exist at such interfaces; the interfacial bound charges will thus have an areal density σ =εo(E2⊥ − E1 ⊥). Under the influence of the local E-field, these charges give rise to an electric Lorentz force F = ½Real(σ E*), where F is the force per unit area of the interface. Since the tangential E-field, E| | , is generally continuous across the interface, there is no ambiguity as to which field must be used in conjunction with the Lorentz law. As for the perpendicular component, the average E across the boundary, ½(E1 ⊥ + E2 ⊥), must be used in calculating the interfacial force. (The use of the average E⊥ in this context is not a matter of choice; it is the only way to get the calculated force at the boundary to agree with the time rate of change of the momentum that passes through the interface. From a physical standpoint, the interfacial charges produce a local E⊥ that has the same magnitude but opposite directions on the two sides of the interface. It is this locally-generated E⊥ that is responsible for the E-field’s discontinuity. Averaging E⊥ across the interface eliminates the local field, as it should, since the charge cannot exert a force on itself.) Since ∇ · B = 0 and B = µ oH, the perpendicular H-field, H⊥ , at the interface between adjacent media must remain continuous. The tangential H-field at such interfaces, however, may be discontinuous. This, in accordance with Maxwell’s equation ∇ × H = J + ∂D/∂t, gives rise to an interfacial current density Js = H2 | | − H1 | | . Such currents can exist on the surfaces of good conductors, where E| | is negligible, yet the high conductance of the medium permits the flow of the surface current. Elsewhere, the only source of electrical currents are bound charges, with the bound current density being Jb = ∂P/∂t = εo(ε – 1)∂E/∂t. Assuming timeharmonic fields with the time-dependence factor exp(−iω t), we can write Jb = −iω εo(ε – 1)E. The H-field of the electromagnetic wave then exerts a force on the bound current according to the Lorentz law, namely, F = ½Real (Jb × B*), where F is force per unit volume. Note: For time-harmonic fields, the contribution of conduction electrons to current density may be combined with that of bound electrons. Since Jc = σcE, where σc is the conductivity of the medium, the net current density Jc + Jb may be attributed to an effective dielectric constant ε + i(σc /εoω ). In general, since ε is complex-valued, there is no need to distinguish conduction electrons from bound electrons, and ε may be treated as an effective dielectric constant that contains both contributions. An exception will be made in Section 3 in the case (C) 2004 OSA 1 November 2004 / Vol. 12, No. 22 / OPTICS EXPRESS 5378 #5025- $15.00 US Received 10 August 2004; revised 13 October 2004; accepted 20 October 2004
of perfect conductors, where Oc -oo. Here both E-and H-fields inside the medium tend to zero, and the contributions of bound charges/currents become negligible. The effect of the radiation in this case will be the creation of a surface current density Js and a surface charge density or, both of which may be attributed in their entirety to the conduction electrons. The conservation of charge then requires that V J+dg/dt=0 3. Reflection of plane wave from a perfect conductor The material in this section is not new, but the line of reasoning and the methodology will be needed to build the necessary arguments in subsequent sections. The case of normal incidence on perfect conductors is well-known, but the two cases of oblique incidence, originally ublished in [10], have all but vanished from modern textbooks H (a) ++-+++-+++--+++ Js cose J Fig. 1. A linearly-polarized plane wave is reflected from a perfectly conducting mirror Whereas the parallel component of the E-field at the mirror surface is zero, the parallel nt of the H-field is at its maximum. The surface current J, is equal in magnitude cular in direction to the tic field at the surface. (a) Normal incidence que incidence with s-polarization. (c) oblique incidence with p-polarizatie In Fig. 1(a) a plane wave of wavelength no, having E-field amplitude Eo (units=V/m) and H-field amplitude Ho=EZ(units=A/m), where Z=vudEo-37722 is the free-space impedance, is incident on a perfect conductor. The Poynting vector is S=rEal(E xH") and the momentum density (per unit volume)is p=S/c(vacuum speed of light c= I/ueo In unit time, the incoming momentum over a unit area of the reflector is that contained Imn of base A= 1.0 m2 and height c. The same momentum returns to the source after ng reflected from the mirror, so the net rate of change of the field momentum over a unit #5025-S1500US Received 10 August 2004; revised 13 October 2004; accepted 20 October 2004 (C)2004OSA November 2004/Vol 12. No 22/OPTICS EXPRESS 5379
of perfect conductors, where σc→ ∞. Here both E- and H-fields inside the medium tend to zero, and the contributions of bound charges/currents become negligible. The effect of the radiation in this case will be the creation of a surface current density Js and a surface charge density σ, both of which may be attributed in their entirety to the conduction electrons. The conservation of charge then requires that ∇ · Js + ∂σ /∂t = 0. 3. Reflection of plane wave from a perfect conductor The material in this section is not new, but the line of reasoning and the methodology will be needed to build the necessary arguments in subsequent sections. The case of normal incidence on perfect conductors is well-known, but the two cases of oblique incidence, originally published in [10], have all but vanished from modern textbooks. Fig. 1. A linearly-polarized plane wave is reflected from a perfectly conducting mirror. Whereas the parallel component of the E-field at the mirror surface is zero, the parallel component of the H-field is at its maximum. The surface current Js is equal in magnitude and perpendicular in direction to the magnetic field at the surface. (a) Normal incidence. (b) Oblique incidence with s-polarization. (c) Oblique incidence with p-polarization. In Fig. 1(a) a plane wave of wavelength λo, having E-field amplitude Eo (units = V/m) and H-field amplitude Ho = Eo/Zo (units = A/m), where Zo = √µo/εo ~ 377Ω is the free-space impedance, is incident on a perfect conductor. The Poynting vector is S = ½Real (E × H*), and the momentum density (per unit volume) is p = S/c 2 (vacuum speed of light c = 1/√µoεo). In unit time, the incoming momentum over a unit area of the reflector is that contained in a column of base A = 1.0 m2 and height c. The same momentum returns to the source after being reflected from the mirror, so the net rate of change of the field momentum over a unit H E Js H E Js cosθ θ θ (a) (b) H E Js θ θ +++ --- +++ --- +++ --- +++ (c) (C) 2004 OSA 1 November 2004 / Vol. 12, No. 22 / OPTICS EXPRESS 5379 #5025- $15.00 US Received 10 August 2004; revised 13 October 2004; accepted 20 October 2004