Radiation pressure and the linear momentum of light in dispersive dielectric media Masud mansuripur Optical Sciences Center, The University of Arona, Tucson, Arisona 85721 Abstract: We derive an exact expression for the radiation pressure of a quasi-monochromatic plane wave incident from the free space onto the flat surface of a semi-infinite dielectric medium In order to account for the total optical momentum (incident plus reflected) that is transferred to the dielectric, the mechanical momentum acquired by the medium must be added to the rate of flow of the electromagnetic momentum( the so-called Abraham momentum) inside the dielectric. We confirm that the lectromagnetic momentum travels with the group velocity of light inside the medium. The photon drag effect in which the photons captured in a emiconductor appear to have the minkowski momentum is explained by analyzing a model system consisting of a thin absorptive layer embedded in I transparent dielectric OCIS codes:(2602110)Electromagnetic theory, (1407010)Trapping. References M. Mansuripur, "Radiation pressure and the linear momentum of the electromagnetic field, "Opt. Express 12 5375-5401(2004),http://www.opticsexpress.org/abstract.cfm?uri=opex-12-22-5375 2. J. P. Gordon, Radiation forces and momenta in dielectric media, "Phys. Rev. A8, 14-21(1973) 3. J.D. Jackson, Classical Electrod namics, 2d edition( Wiley, New York, 1975) 4. R Loudon, "Theory of the radiation pressure on dielectric surfaces, J Mod. Opt. 49, 821-838(2002) 5. Y.N. Obukhov and F. w. Hehl, "Electromagnetic energy-momentum and forces in matter, " Phys. Lett. A 311 277-284(2003) 6. M. Mansuripur, A R. Zakharian, and J. V Moloney, "Radiation pressure on a dielectric wedge, accepted for 7. A.F. Gibson, M. F Kimmitt, and A C. Walker, "Photon drag in Germanium, Appl. Phys. Lett. 17, 75-77 (1970) 8. R Loudon, S M. Barnett, and C. Baxter, "Theory of radiation pressure and momentum transfer in dielectrics: the photon drag effect, "to appear in 2005 1 Introduction In a previous paper [1] we showed that the momentum density p of a plane electro-magnetic wave inside a dispersionless dielectric medium may be expressed as the average of the Minkowski and Abraham momentum densities [2], namely, p=v Real(ExH")c2+ 4 Real(D x B*) Here the complex amplitudes of the electric and magnetic fields within the medium of refractive index n are denoted by E and H, respectively; B=4H, and D= E+P=EEE, where P is the polarization density induced in the medium by the local E-field. E is the permittivity and uo the permeability of free space; E=n is the relative permittivity of the dielectric material [3]. We derived Eq (la) by a direct application of the Lorentz law of force to bound charges and bound currents within the medium- a method that has been the subject of other recent studies as well [4, 5]. It was concluded that the light carries its own electro 6629-$1500US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005 (C)2005OSA 21 March 2005/ Vol 13. No 6/ OPTICS EXPRESS 2245
Radiation pressure and the linear momentum of light in dispersive dielectric media Masud Mansuripur Optical Sciences Center, The University of Arizona, Tucson, Arizona 85721 masud@u.arizona.edu Abstract: We derive an exact expression for the radiation pressure of a quasi- monochromatic plane wave incident from the free space onto the flat surface of a semi-infinite dielectric medium. In order to account for the total optical momentum (incident plus reflected) that is transferred to the dielectric, the mechanical momentum acquired by the medium must be added to the rate of flow of the electromagnetic momentum (the so-called Abraham momentum) inside the dielectric. We confirm that the electromagnetic momentum travels with the group velocity of light inside the medium. The photon drag effect in which the photons captured in a semiconductor appear to have the Minkowski momentum is explained by analyzing a model system consisting of a thin absorptive layer embedded in a transparent dielectric. © 2005 Optical Society of America OCIS codes: (260.2110) Electromagnetic theory; (140.7010) Trapping. References 1. M. Mansuripur, “Radiation pressure and the linear momentum of the electromagnetic field,” Opt. Express 12, 5375-5401 (2004), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-12-22-5375. 2. J. P. Gordon, “Radiation forces and momenta in dielectric media,” Phys. Rev. A 8, 14-21 (1973). 3. J. D. Jackson, Classical Electrodynamics, 2nd edition (Wiley, New York, 1975). 4. R. Loudon, “Theory of the radiation pressure on dielectric surfaces,” J. Mod. Opt. 49, 821-838 (2002). 5. Y. N. Obukhov and F. W. Hehl, “Electromagnetic energy-momentum and forces in matter,” Phys. Lett. A 311, 277-284 (2003). 6. M. Mansuripur, A. R. Zakharian, and J. V. Moloney, “Radiation pressure on a dielectric wedge,” accepted for publication, Opt. Express, 2005. 7. A. F. Gibson, M. F. Kimmitt, and A. C. Walker, “Photon drag in Germanium,” Appl. Phys. Lett. 17, 75-77 (1970). 8. R. Loudon, S. M. Barnett, and C. Baxter, “Theory of radiation pressure and momentum transfer in dielectrics: the photon drag effect,” to appear in 2005. 1. Introduction In a previous paper [1] we showed that the momentum density p of a plane electro-magnetic wave inside a dispersionless dielectric medium may be expressed as the average of the Minkowski and Abraham momentum densities [2], namely, p = ¼ Real (E × H*)/c2 + ¼ Real (D × B*). (1a) Here the complex amplitudes of the electric and magnetic fields within the medium of refractive index n are denoted by E and H, respectively; B = µoH, and D = εoE + P = εoε E, where P is the polarization density induced in the medium by the local E-field. εo is the permittivity and µo the permeability of free space; ε = n2 is the relative permittivity of the dielectric material [3]. We derived Eq.(1a) by a direct application of the Lorentz law of force to bound charges and bound currents within the medium – a method that has been the subject of other recent studies as well [4,5]. It was concluded that the light carries its own electro- (C) 2005 OSA 21 March 2005 / Vol. 13, No. 6 / OPTICS EXPRESS 2245 #6629 - $15.00 US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005
magnetic momentum inside the dielectric, while an additional momentum is transferred to the medium in the form of mechanical force. Equation(la)may be rewritten as follows p=Rea(E×H)a2+%Real(P×B*) (1b) In the above equation, the first term is the Abraham momentum density of the field, while the second term is the mechanical momentum density imparted to the medium. (If the coefficient of the second term were 2 instead of / 4, the total momentum density p would have been equal to the Minkowski momentum. )The electromagnetic and mechanical momenta of the light inside the dielectric are not decoupled from each other. This fact is better appreciated if one observes, for instance, that the same beam of light, upon emerging into the free-space at the exit facet of a dielectric slab, recovers its total initial momentum (i.e, the momentum it possessed before entering the slab) by re-converting the mechanical momentum(manifested in the motion of the medium)to electromagnetic momentum [1]. If, for simplicity' s sake, we assume that the entrance and exit facets of the dielectric slab are anti-reflection coated then upon transmission, the emerging beams momentum will be identical to the momentum it possessed before entering the slab; in other words, the(partial)conversion of the beam's momentum into mechanical form that takes place while the beam passes through the slab, is fully reversed when the beam leaves the slab and returns to the free space. Another example of the"connectedness"of the electromagnetic and mechanical momenta was provided in [6] where the radiation pressure on a dielectric wedge and its surrounding liquid was found to arise from the total momentum of the beam as opposed to, say, from one or the other of its constituents The present paper extends the results of our previous work to the case of light beams that propagate in dispersive dielectrics. We show that the earlier results obtained for non- dispersive media remain valid if the abraham momentum is assumed to travel with the group velocity Vg=c/n+nf) inside the dielectric. (f is the optical frequency, n=dn/df is the derivative of the refractive index n )Our new expression for the mechanical momentum density reverts to the old expression, 74 Real(Px B ) in the limit of n0, i.e., when the medium becomes dispersionless. Finally, to resolve the discrepancy between the theory and certain experiments in which the light appears to possess the Minkowski momentum, we propose a model system for analyzing the photon drag effect observed in certain bulk semiconductors 2. Superposition of two plane waves in free space Figure I shows a beam of light consisting of two equal-amplitude plane-waves of differing frequencies f i and i, incident on a semi-infinite dielectric medium of refractive index n(f) The beam is linearly polarized, having its E-field along the x-axis and H-field along the y axis. The field amplitudes in free space are given by Ex=,t=Eo sin(2T fi[(=c)-t1-Eo sin(2Tf2((/)-JF H(=,D)=(EZ)sin{2f[(/c)-l}-(E/Z)sin{2f[(=lc)-l]}(2b) Here Zo=NudE is the free-space impedance, and c= 1/lE is the speed of light in vacuum The traveling wave is readily seen to be a sinusoid of frequency f=7(fi+/2), modulated with another(envelope)sinusoid of frequency Af=fi-fi, exhibiting a beat period T=1 and traveling(in free-space)with the speed c The Poynting vector S=ExH has a component only along the z-axis, S(=,n) Ere, OH, n). For a fixed value of = if S(=, t) is averaged over the time interval T, and if terms of order Af and higher are neglected, the time-averaged Poynting vector will bece independent of the coordinate : and will be given by 6629-$1500US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005 (C)2005OSA 21 March 2005/Vol 13. No 6/ OPTICS EXPRESS 2246
magnetic momentum inside the dielectric, while an additional momentum is transferred to the medium in the form of mechanical force. Equation (1a) may be rewritten as follows: p = ½Real (E × H*)/c2 + ¼ Real (P × B*). (1b) In the above equation, the first term is the Abraham momentum density of the field, while the second term is the mechanical momentum density imparted to the medium. (If the coefficient of the second term were ½ instead of ¼, the total momentum density p would have been equal to the Minkowski momentum.) The electromagnetic and mechanical momenta of the light inside the dielectric are not decoupled from each other. This fact is better appreciated if one observes, for instance, that the same beam of light, upon emerging into the free-space at the exit facet of a dielectric slab, recovers its total initial momentum (i.e., the momentum it possessed before entering the slab) by re-converting the mechanical momentum (manifested in the motion of the medium) to electromagnetic momentum [1]. If, for simplicity’s sake, we assume that the entrance and exit facets of the dielectric slab are anti-reflection coated, then, upon transmission, the emerging beam’s momentum will be identical to the momentum it possessed before entering the slab; in other words, the (partial) conversion of the beam’s momentum into mechanical form that takes place while the beam passes through the slab, is fully reversed when the beam leaves the slab and returns to the free space. Another example of the “connectedness” of the electromagnetic and mechanical momenta was provided in [6], where the radiation pressure on a dielectric wedge and its surrounding liquid was found to arise from the total momentum of the beam as opposed to, say, from one or the other of its constituents. The present paper extends the results of our previous work to the case of light beams that propagate in dispersive dielectrics. We show that the earlier results obtained for nondispersive media remain valid if the Abraham momentum is assumed to travel with the group velocity Vg = c/(n+ n′ f ) inside the dielectric. ( f is the optical frequency; n′ = dn/df is the derivative of the refractive index n.) Our new expression for the mechanical momentum density reverts to the old expression, ¼Real(P×B*), in the limit of n′→0, i.e., when the medium becomes dispersionless. Finally, to resolve the discrepancy between the theory and certain experiments in which the light appears to possess the Minkowski momentum, we propose a model system for analyzing the photon drag effect observed in certain bulk semiconductors. 2. Superposition of two plane waves in free space Figure 1 shows a beam of light consisting of two equal-amplitude plane-waves of differing frequencies f1 and f2, incident on a semi-infinite dielectric medium of refractive index n( f ). The beam is linearly polarized, having its E-field along the x-axis and H-field along the yaxis. The field amplitudes in free space are given by Ex(z, t) = Eo sin{2π f1 [(z /c) – t]} – Eo sin{2π f2 [(z /c) – t]} (2a) Hy(z, t) = (Eo/Zo)sin{2π f1 [(z /c) – t]} – (Eo/Zo)sin{2π f2 [(z /c) – t]} (2b) Here Zo = √µo/εo is the free-space impedance, and c = 1/√µoεo is the speed of light in vacuum. The traveling wave is readily seen to be a sinusoid of frequency f = ½( f1 + f2), modulated with another (envelope) sinusoid of frequency ∆f = f2 – f1, exhibiting a beat period T= 1/∆f, and traveling (in free-space) with the speed c. The Poynting vector S = E×H has a component only along the z-axis, Sz(z, t) = Ex(z, t)Hy(z, t). For a fixed value of z, if Sz(z, t) is averaged over the time interval T, and if terms of order ∆f and higher are neglected, the time-averaged Poynting vector will become independent of the coordinate z, and will be given by (C) 2005 OSA 21 March 2005 / Vol. 13, No. 6 / OPTICS EXPRESS 2246 #6629 - $15.00 US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005
S(,1)>=E21Z乙 (The above equation will be exact, i.e., terms in Af and higher order will be absent, iff=NAf, e, if the center frequency happens to be an integer-multiple of Af) Each participating frequency thus contributes its own(time-averaged) Poynting vector = yEo Zo to the energy flux of the beam. The time-averaged momentum density of the beam (i.e, momentum per unit volume), P:=/c, is thus uniform throughout the free-space region, having equal contributions from the two frequency components of the beam Reflected Refractive index=以(/ Www Incident beam light beam consisting of two equal-amplitude plane-waves of slightly is linearly polarized, having its E-field along the x-axis and H-field along action of the beam is reflected at the surface. the remainder enters the ng it at the group velocity g=c/(n+nf). Heref=7(fi+f2)is the center frequency, and n=dn/df; both n and n are evaluated at the center frequency Suppose now that the beam reaches the interface with a dielectric medium of refractive index n(f), as shown in Fig. 1. Each frequency component gets reflected at the interface, with an(amplitude)reflection coefficient p=(1-n)(1+ n). The time rate of change of the optica momentum g on the vacuum side of the interface is then equal to the rate of arrival of the incident momentum, cp:=Eo/(cZ=Eo Eo plus the rate of departure of the reflected optical momentum,namely 1+(1-n1)(1+n2+k(1-m)(1+n2)2}eE2 In the limit -f2, the refractive indices n1, n2 in Eq (4a)are nearly identical, and the above expression simplifies to dg /dr=2(n+1)(n+1)jeE where n=7(n1 + n2). Momentum conservation requires the above dqs /dt to be balanced by the force exerted on the dielectric medium plus the time rate of change of any momentum taken by the transmitted beam into the dielectric 3. Superposition of two plane waves in a dielectric medium Inside the medium, each frequency component arrives with an amplitude transmission coefficient(for the E-field) given by T=2/(1+n). The transmitted E- and H-fields may thus be written as follows: E(=,D)=[2E(n1+1)sin{2xf[(m11)-l}-[2EM(m2+1)sin{2r/[(n2=lc)-l}(5a) 6629-$1500US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005 (C)2005OSA 21 March 2005/ Vol 13. No 6/ OPTICS EXPRESS 2247
= Eo 2 /Zo (3) (The above equation will be exact, i.e., terms in ∆f and higher order will be absent, if f = N ∆f, i.e., if the center frequency happens to be an integer-multiple of ∆f.) Each participating frequency thus contributes its own (time-averaged) Poynting vector = ½Eo 2 /Zo to the energy flux of the beam. The time-averaged momentum density of the beam (i.e., momentum per unit volume), pz = /c 2 , is thus uniform throughout the free-space region, having equal contributions from the two frequency components of the beam. Fig. 1. A light beam consisting of two equal-amplitude plane-waves of slightly differing frequencies, f1 and f2, is normally incident on a semi-infinite dielectric of refractive index n( f ). The beam is linearly polarized, having its E-field along the x-axis and H-field along the y-axis. While a fraction of the beam is reflected at the surface, the remainder enters the dielectric, penetrating it at the group velocity Vg = c/(n + n′ f ). Here f = ½( f1 + f2) is the center frequency, and n′ = dn/d f ; both n and n′ are evaluated at the center frequency. Suppose now that the beam reaches the interface with a dielectric medium of refractive index n( f ), as shown in Fig. 1. Each frequency component gets reflected at the interface, with an (amplitude) reflection coefficient ρ = (1 – n)/(1 + n). The time rate of change of the optical momentum q on the vacuum side of the interface is then equal to the rate of arrival of the incident momentum, c pz = Eo 2 /(cZo) = εoEo 2 , plus the rate of departure of the reflected optical momentum, namely, d qz/dt = {1 + ½|(1 – n1)/(1 + n1)| 2 + ½|(1 – n2)/(1 + n2)| 2 }εo Eo 2 . (4a) In the limit f1 → f2, the refractive indices n1, n2 in Eq.(4a) are nearly identical, and the above expression simplifies to d qz /dt = 2[(n2 + 1)/(n + 1)2 ]εoEo 2 , (4b) where n = ½(n1 + n2). Momentum conservation requires the above d qz/dt to be balanced by the force exerted on the dielectric medium plus the time rate of change of any momentum taken by the transmitted beam into the dielectric. 3. Superposition of two plane waves in a dielectric medium Inside the medium, each frequency component arrives with an amplitude transmission coefficient (for the E-field) given by τ = 2/(1+ n). The transmitted E- and H-fields may thus be written as follows: Ex(z, t)= [2Eo/(n1 + 1)]sin{2π f1 [(n1 z /c) – t]}– [2Eo/(n2 + 1)]sin{2π f2 [(n2 z /c) – t]} (5a) Incident beam Transmitted X Reflected Refractive index = n( f ) Z (C) 2005 OSA 21 March 2005 / Vol. 13, No. 6 / OPTICS EXPRESS 2247 #6629 - $15.00 US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005
,1={2mE[Ln1+1)}sin{2rf(n/c)-l]} -{2nE[Zm2+1)]}sin{2xf2[(n2z/c)-l} As before, the Poynting vector may be calculated and time-averaged over one beat period T Once again, ignoring the terms in A and higher order(or assuming= NAf so that these terms would automatically vanish), we find the average rate of flow of energy(per unit area per unit time)within the dielectric medium to be independent of z and given by =%{4n(m1+1)+[4m(m2+1)}E21Z乙 The electromagnetic(Abraham) momentum inside the dielectric whose volumetric density P=/c must propagate with the group velocity V,=c(2-fiV(2f2-nifi)=c/(n+n) Here n=7(m+n2),f=2(1+ 2), and n=(n2-n)/(fi-fi is the derivative of n with respect to frequency. In the limit i-2, the refractive indices n, and n2 appearing in Eq(6) are nearly identical, and the time rate of flow of the electromagnetic mor medium is given by dq /dt=p k=4nE Eo/[(n+nf)(n+1)] (7) Clearly, this is not equal to the momentum flux across the vacuum-dielectric interface given by Eq (4b). The reason for the discrepancy is that some of the incoming momentum has been converted into mechanical force exerted on the dielectric medium. To compute this force w use the Lorentz law, taking into account the fact that the polarization density within the medium is P(=,t)=E(E-lE(, n), where E=n is a function of the frequency f The bound current density is thus given by JME, t=aP(a, t/ar=-4Ifi(n-1)Eo Eo cos 2I fi [(n=/c)-1 +4If2(n -1) Eo cos(2f(n =/c)-1) The Lorentz force density F=, t=JxXHo Hy is obtained by multiplying Eqs. 8)and(5b) At time t the leading edge of the beat"waveform"(see Fig. 1) has penetrated a distance =o=Vgt into the dielectric. The integrated force density from ==0 to =o thus yields the force per unit cross-sectional area, F:(0, exerted at time t on the medium. We find F(1)(26E2)={[n2-(mh-m2f)(m22-n/[(m+1)n2+1)}{1-cos2m(f-fn]} -{(m2-m)f2-f)/[(n1+1)n2+1)mnf+n2f2)}cos[4an2-n)t/m2-n1f +{[lnn2-(mf2+n2fn1f+n2/2)/[(n1+1)n2+1)}cos{2m(f1+f) [(n1-1)(n1+1)]cos(4xft)-%(n2-1)(n2+1)cos(4f1 Averaging F:(n) over the duration of a single beat waveform, T=1/(f2-fi), we find that second, third, fourth, and fifth terms in the above expression contribute very little to average force, and that the only significant contribution arises from the first term, namely, =(1/T)/F0dr=2eE2[mn2-(n-n2f(n2压-m)/[(m+1)m2+1)](10) In the limit f2-fi the above expression simplifies to yield the net(average) force per unit cross-sectional area exerted on the dielectric medium as follow {n2-[(n-nn+nf)]}E。n+1) 6629-$1500US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005 (C)2005OSA 21 March 2005/VoL 13. No 6/ OPTICS EXPRESS 2248
Hy(z, t) = {2n1Eo/[Zo(n1 + 1)]} sin{2πf1 [(n1 z/c) – t ]} – {2n2Eo/[Zo(n2 + 1)]}sin{2π f2 [(n2 z/c) – t ]} (5b) As before, the Poynting vector may be calculated and time-averaged over one beat period T. Once again, ignoring the terms in ∆f and higher order (or assuming f = N∆f so that these terms would automatically vanish), we find the average rate of flow of energy (per unit area per unit time) within the dielectric medium to be independent of z and given by = ½{[4n1/(n1 + 1)2 ]+[4n2 /(n2 + 1)2 ]}Eo 2 /Zo. (6) The electromagnetic (Abraham) momentum inside the dielectric whose volumetric density pz = /c 2 must propagate with the group velocity Vg = c( f2 – f1)/(n2 f2 – n1 f1) = c/(n + n′ f ). Here n = ½(n1 + n2), f = ½( f1 + f2), and n′ = (n2 – n1)/( f2 – f1) is the derivative of n with respect to frequency. In the limit f1 → f2, the refractive indices n1 and n2 appearing in Eq. (6) are nearly identical, and the time rate of flow of the electromagnetic momentum through the medium is given by d qz /dt = pzVg = 4nεoEo 2 /[(n + n′ f )(n + 1)2 ]. (7) Clearly, this is not equal to the momentum flux across the vacuum-dielectric interface given by Eq.(4b). The reason for the discrepancy is that some of the incoming momentum has been converted into mechanical force exerted on the dielectric medium. To compute this force we use the Lorentz law, taking into account the fact that the polarization density within the medium is P(z, t) = εo(ε – 1)E(z, t), where ε = n2 is a function of the frequency f. The bound current density is thus given by Jx(z, t) = ∂Px(z, t)/∂t = – 4πf1 (n1 – 1)εoEo cos{2π f1 [(n1z /c) – t]} + 4π f2 (n2 – 1)εoEo cos{2π f2 [(n2 z /c) – t]}. (8) The Lorentz force density Fz(z, t)= Jx ×µoHy is obtained by multiplying Eqs.(8) and (5b). At time t the leading edge of the beat “waveform” (see Fig. 1) has penetrated a distance zo = Vg t into the dielectric. The integrated force density from z = 0 to zo thus yields the force per unit cross-sectional area, Fz(t), exerted at time t on the medium. We find Fz(t)/(2εoEo 2 ) = {[n1n2 – (n1 f2 – n2 f1)/(n2 f2 – n1 f1)]/[(n1+ 1)(n2 + 1)]}{1 – cos[2π( f2 – f1)t]} – {(n2 – n1)( f2 – f1)/[(n1+ 1)(n2 + 1)(n1 f1+ n2 f2)]}cos[4π(n2 – n1) f1 f2 t /(n2 f2 – n1 f1)] + {[n1n2 – (n1 f2+ n2 f1)/(n1 f1+ n2 f2)]/[(n1 + 1)(n2+ 1)]}cos[2π( f1 +f2)t] –½[(n1 – 1)/(n1 + 1)] cos(4π f1 t) –½[(n2 – 1)/(n2 + 1)] cos(4π f2 t). (9) Averaging Fz(t) over the duration of a single beat waveform, T= 1/( f2 – f1), we find that the second, third, fourth, and fifth terms in the above expression contribute very little to the average force, and that the only significant contribution arises from the first term, namely, = (1/T )∫ Fz(t) dt = 2εoEo 2 [n1 n2 – (n1 f2 – n2 f1)/(n2 f2 – n1 f1)]/[(n1 + 1)(n2 + 1)]. (10) In the limit f2 → f1 the above expression simplifies to yield the net (average) force per unit cross-sectional area exerted on the dielectric medium as follows: = 2{n2 – [(n – n′ f )/(n + n′ f )]}εoEo 2 /(n+1)2 . (11) 0 T (C) 2005 OSA 21 March 2005 / Vol. 13, No. 6 / OPTICS EXPRESS 2248 #6629 - $15.00 US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005
In averaging F(n)of Eq (9)over the time interval T, we neg but the first term. This is easily justified for the third, fourth, and fifth terms, which exhibit rapid oscillations. However, the second term is harder to neglect, especially when n=0(i.e for a nearly dispersionless medium), because the cosine term under these circumstances is not rapidly oscillating. However, the magnitude of the second term is proportional to AnAf, which reduces the term's significance when the cosine is weakly oscillating. All in all, the time average of the second term in Eq (9)turns out to be negligible even when the medium is dispersionless(or nearly so) Adding of Eq (11)to the rate of flow of the electromagnetic (i.e, Abraham) momentum inside the dielectric given by Eq (7) yields PV+=2[(m2+1)(n+1)]eE This is identical to the total momentum per unit area per unit time imparted to the dielectric as given by Eq (4b). We have thus confirmed the conservation of momentum in the system of Fig. I by deriving the expression for the radiation pressure, Eq (11), and by requiring propagation at the group velocity Vg for the electromagnetic momentum inside the medium; see Eq (7) If mechanical momentum is assumed to travel through the dielectric with the group velocity Vs a mechanical momentum density P(ech)=/Vs can be defined which with the electromagnetic momentum density p total optical momentum residing in the medium. In the limit when n0, for each frequency component of the beam (i.e, fi,f2), p mech -V(E-1)/c2, in agreement with the result obtained in [1] for dispersionless media 4. The photon drag effect An intriguing experimental observation in certain (weakly absorbing) semiconductors notably Si and Ge, is the photon drag effect [7]. When a photon of energy hf from a monochromatic beam(vacuum wavelength Mo=c/f) is absorbed within a semiconductor of refractive index n(i. e, the real part of the complex refractive index n+ ik, whose imaginary part is the absorption coefficient K), the excited charge carrier acquires a momentum equal to thf/c, the so-called Minkowski momentum of the photon. In contrast, the photons Abraham momentum-seen from the preceding sections discussions to be hf/[(n+nf)cl-is clearly different from the Minkowski momentum. Combining the photons electromagnetic and mechanical momenta does not resolve the discrepancy either, as the total photon momentum v[n+(I/n)] hf/c, differs from the Minkowski momentum as well Loudon et al [ 8] have given a comprehensive theory of the photon drag effect, arguing that the"transparent part"of the semiconductor(associated with the real part of the complex refractive index)takes up the difference between the photon's momentum and the Minkowski momentum, when the latter is transferred to the "absorbing part"of the material (i.e, the part associated with the imaginary component of the complex refractive index). We present a similar(though by no means identical)explanation of the photon drag effect by showing that the momentum picked up by a thin absorbing layer embedded in a transparent dielectric is equal to the minkowski momentum of the incident photon With reference to Fig. 2. the reflection and transmission coefficients of g lay of thickness d and complex index n+ix, in the limit of d<<no(where A is the vacuum wavelength of the incident beam), can be shown to be p=-[1+(k/2n)](2xxdD) (13a) 1-(2兀kd)+i[2n-(x2n)] (13b) 6629-$1500US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005 (C)2005OSA 21 March 2005/ Vol 13. No 6/ OPTICS EXPRESS 2249
In averaging Fz(t) of Eq.(9) over the time interval T, we neglected the contributions of all but the first term. This is easily justified for the third, fourth, and fifth terms, which exhibit rapid oscillations. However, the second term is harder to neglect, especially when n′ ≈ 0 (i.e., for a nearly dispersionless medium), because the cosine term under these circumstances is not rapidly oscillating. However, the magnitude of the second term is proportional to ∆n∆f, which reduces the term’s significance when the cosine is weakly oscillating. All in all, the time average of the second term in Eq.(9) turns out to be negligible even when the medium is dispersionless (or nearly so). Adding of Eq.(11) to the rate of flow of the electromagnetic (i.e., Abraham) momentum inside the dielectric given by Eq.(7) yields: pzVg + = 2[(n2 + 1)/(n + 1)2 ]εoEo 2 . (12) This is identical to the total momentum per unit area per unit time imparted to the dielectric as given by Eq.(4b). We have thus confirmed the conservation of momentum in the system of Fig. 1 by deriving the expression for the radiation pressure, Eq.(11), and by requiring propagation at the group velocity Vg for the electromagnetic momentum inside the medium; see Eq.(7). If mechanical momentum is assumed to travel through the dielectric with the group velocity Vg, a mechanical momentum density pz (mech) = /Vg can be defined which, combined with the electromagnetic momentum density pz (Abraham) = /c 2 , accounts for the total optical momentum residing in the medium. In the limit when n′ →0, for each frequency component of the beam (i.e., f1, f2), pz (mech) →½(ε – 1) /c 2 , in agreement with the result obtained in [1] for dispersionless media. 4. The photon drag effect An intriguing experimental observation in certain (weakly absorbing) semiconductors, notably Si and Ge, is the photon drag effect [7]. When a photon of energy hf from a monochromatic beam (vacuum wavelength λo = c/f ) is absorbed within a semiconductor of refractive index n (i.e., the real part of the complex refractive index n + iκ, whose imaginary part is the absorption coefficient κ ), the excited charge carrier acquires a momentum equal to nhf /c, the so-called Minkowski momentum of the photon. In contrast, the photon’s Abraham momentum – seen from the preceding section’s discussions to be hf/[(n+ n′ f ) c] – is clearly different from the Minkowski momentum. Combining the photon’s electromagnetic and mechanical momenta does not resolve the discrepancy either, as the total photon momentum, ½[n + (1/n)] hf /c, differs from the Minkowski momentum as well. Loudon et al [8] have given a comprehensive theory of the photon drag effect, arguing that the “transparent part” of the semiconductor (associated with the real part of the complex refractive index) takes up the difference between the photon’s momentum and the Minkowski momentum, when the latter is transferred to the “absorbing part” of the material (i.e., the part associated with the imaginary component of the complex refractive index). We present a similar (though by no means identical) explanation of the photon drag effect by showing that the momentum picked up by a thin absorbing layer embedded in a transparent dielectric is equal to the Minkowski momentum of the incident photon. With reference to Fig. 2, the reflection and transmission coefficients of an absorbing layer of thickness d and complex index n + iκ, in the limit of d << λo (where λo is the vacuum wavelength of the incident beam), can be shown to be ρ = − [1 + i(κ /2n)](2πκ d/λo), (13a) τ = 1 − (2πκ d/λo) + i[2n – (κ2 /n)] πd/λo. (13b) (C) 2005 OSA 21 March 2005 / Vol. 13, No. 6 / OPTICS EXPRESS 2249 #6629 - $15.00 US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005
The above equations are obtained by assuming the existence of a pair of counter-propagating plane-waves in the absorbing layer, then determining the various unknown amplitudes by matching the boundary conditions for the E-and H-fields while ignoring second-and higher order terms in d/no. (The algebra is straightforward but tedious. The absorbed optical power y within the layer is thus given by 2(1-|p2-|r)nE21Z=(2nKd)E21z X E TE E0/Z0 Fig. 2. A thin absorptive layer of thickness d and complex refractive index n+ in a transparent, homogeneous dielectric medium of refractive index n(same n of the complex index of the absorptive layer ). A monochromatic plane-wave, wavelength Mo=c/f E-field amplitude Eo, and H-field amplitude Ho=nEo coefficients are p and I, respectively. Each absorbed of energy hf equivalent of its Minkowski momentum nhf le to the absorbing layer 6. The force experienced by the absorbing layer may be derived from the Lorentz law, lowing the procedure outlined in [1], then integrated over the film thickness d to yield =(2Tn'xd/oJEEO Thus =(n/c)% namely, the force experienced by the absorbing layer is (n/c) times the captured optical power. Consequently, the momentum transferred to the layer in a given time interval At must also be(n/c)times the energy absorbed by the layer during the same time interval. For a captured photon of energy hf the momentum transfer is thus equal to nhflc,i.e the Minkowski momentum of the photon in its dielectric environment. Since this is greater than the total photon momentum prior to being absorbed, the host medium(i.e, the dielectric) must experience a recoil equal to the difference between the incident photon,s initial momentum and the Minkowski value picked up by the excited charge carrier Acknowledgments The author is grateful to Ewan Wright, Armis Zakharian, Pavel Polynkin, and Walter Hoyer for many helpful discussions. This work has been supported by the AFOSR contract F49620 02-1-0380 with the Joint Technology Office, by the Office of Naval Research mURi grant No.NO0014-03-1-0793, , and he National Science Foundation STC Pro agreement DMR-0120967 6629-$1500US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005 (C)2005OSA 21 March 2005/ Vol 13. No 6/ OPTICS EXPRESS 2250
The above equations are obtained by assuming the existence of a pair of counter-propagating plane-waves in the absorbing layer, then determining the various unknown amplitudes by matching the boundary conditions for the E- and H-fields while ignoring second- and higherorder terms in d/λo. (The algebra is straightforward but tedious.) The absorbed optical power γ within the layer is thus given by γ = ½(1 – |ρ |2 − |τ |2 ) nEo 2 /Zo = (2π nκ d /λo)Eo 2 /Zo. (14) Fig. 2. A thin absorptive layer of thickness d and complex refractive index n + iκ is embedded in a transparent, homogeneous dielectric medium of refractive index n (same n as the real part of the complex index of the absorptive layer). A monochromatic plane-wave, having vacuum wavelength λo = c/f, E-field amplitude Eo, and H-field amplitude Ho = nEo/Zo, is normally incident on the absorbing layer. The layer’s (amplitude) reflection and transmission coefficients are ρ and τ, respectively. Each absorbed photon of energy hf transfers the equivalent of its Minkowski momentum nhf /c to the absorbing layer. The force experienced by the absorbing layer may be derived from the Lorentz law, following the procedure outlined in [1], then integrated over the film thickness d to yield = (2πn2 κ d/λo)εoEo 2 . (15) Thus = (n/c)γ, namely, the force experienced by the absorbing layer is (n/c) times the captured optical power. Consequently, the momentum transferred to the layer in a given time interval ∆t must also be (n/c) times the energy absorbed by the layer during the same time interval. For a captured photon of energy hf the momentum transfer is thus equal to nhf /c, i.e., the Minkowski momentum of the photon in its dielectric environment. Since this is greater than the total photon momentum prior to being absorbed, the host medium (i.e., the dielectric) must experience a recoil equal to the difference between the incident photon’s initial momentum and the Minkowski value picked up by the excited charge carrier. Acknowledgments The author is grateful to Ewan Wright, Armis Zakharian, Pavel Polynkin, and Walter Hoyer for many helpful discussions. This work has been supported by the AFOSR contract F49620- 02-1-0380 with the Joint Technology Office, by the Office of Naval Research MURI grant No.N00014-03-1-0793, and by the National Science Foundation STC Program under agreement DMR-0120967. d n ρEo n +iκ Eo τEo Z X Ho = nEo /Zo · · × (C) 2005 OSA 21 March 2005 / Vol. 13, No. 6 / OPTICS EXPRESS 2250 #6629 - $15.00 US Received 18 February 2005; revised 14 March 2005; accepted 15 March 2005