arxiv: v3 [physics.optics] 14 May 2012

Transcription

1 Elecromagnec energy-momenum n dspersve meda T. G. Phlbn School of Physcs and Asronomy, Unversy of S Andrews, Norh Haugh, S Andrews, Ffe, KY16 9SS, Scoland, Uned Kngdom arxv: v3 [physcs.opcs] 1 May 01 The sandard dervaons of elecromagnec energy and momenum n meda ake Maxwell s equaons as he sarng pon. I s well known ha for dspersve meda hs approach does no drecly yeld exac expressons for he energy and momenum denses. Alhough Maxwell s equaons fully descrbe elecromagnec felds, he general approach o conserved quanes n feld heory s no based on he feld equaons, bu raher on he acon. Here an acon prncple for macroscopc elecromagnesm n dspersve, lossless meda s used o derve he exac conserved energy-momenum ensor. The me-averaged energy densy reduces o Brlloun s smple formula when he felds are monochromac. The momenum densy s no gven by he famlar Mnkowsk expresson D B, even for me-averaged monochromac felds. The resuls are unaffeced by he debae over momenum balance n lgh-maer neracons. PACS numbers: De,.5.Bs,.50.Wk [Noe: Ths arxv verson ncorporaes he Erraum o he publshed verson [1].] I. INTRODUCTION Momenum balance n lgh-maer neracons s a suble subjec ha has been debaed for over a cenury. Dscussons of he ssues nvolved, ogeher wh references o much of he relevan leraure, can be found n [ 6]. Ths paper makes no aemp o conrbue o he momenum-balance debae, n ha here wll be no nvesgaon of he ransfer of momenum (or energy) from lgh o maer or vce versa. Insead we solae a relaed problem, bu one ha s mahemacally wellposed and herefore capable of an exac and unambguous soluon. The problem solved here s: wha s he conserved energy-momenum ensor of lgh propagang n a homogeneous, dspersve, lossless magneodelecrc medum? Ths s a well-posed problem of macroscopc elecromagnesm, bu appears no o have been fully addressed n he exsng leraure on elecromagnec energy-momenum, where he prmary focus has been on momenum ransfer beween lgh and maer. Alhough he soluon s an example of Noeher s heorem, he presence of dsperson means ha he well-known expresson for he energy-momenum ensor of a feld heory n erms of s Lagrangan densy s no vald n hs case. In fac he resul canno be derved n any sysemac way from Maxwell s equaons or a Lagrangan densy. Dspersve macroscopc elecromagnesm mus be formulaed n erms of an acon prncple, and he Noeher heorem for spaceme ranslaon symmery mus be derved from he acon, leadng o an exac conserved energy-momenum ensor. The problem n fac presens a sgnfcan calculaonal challenge, whch may explan s absence from he momenum-balance leraure. By solang hs problem from he conundrum of momenum ransfer, we show ha he ssue of elecromagnec energy-momenum n dspersve meda s amenable o some mahemacally exac and unambguous saemens. The resrcon o lossless meda n hs paper requres commen on he physcal sgnfcance of he resuls. A complee reamen of macroscopc elecromagnesm a all frequences mus nclude he absorpon n he maerals ha s necessarly srong n some frequency ranges. Mahemacally, hs dsspaon s a consequence of he resrcon o rearded soluons of Maxwell s equaons, whch leads o he Kramers-Krong relaons [7, 8]. As s well known [8 10], quanfyng elecromagnec energy-momenum n dsspave meda s problemac, and n any case such energy-momenum wll of course no be conserved. In opcs, on he oher hand, one ofen deals wh lmed frequency ranges where losses are neglgble, bu where maeral dsperson canno be gnored. An example of he umos expermenal and praccal mporance s fber opcs [11], where n many crcumsances absorpon can be gnored no only n he enre vsble range, bu also no he nfrared and ulravole. In such a large frequency range dsperson plays a crucal role n lgh propagaon and canno be negleced [11]. The queson of how much conserved energy-momenum s beng ranspored by lgh has an unambguous answer n hese crcumsances, and he resul s expermenally sgnfcan. For example, a lgh beam encounerng an nense pulse can be frequency shfed hrough nonlneary of he medum, and he amoun of elecromagnec energy ha s frequency shfed depends, among oher facors, on he energy-momenum ensor of he lgh beam. One example of such an expermen s descrbed n [1], where he lgh beam encounerng he nense pulse was monochromac; hs represens he smples case, bu for a lgh beam wh a complcaed specrum he general expresson for he energy-momenum ensor s requred o fully descrbe he behavour of he beam. Ths exac energy-momenum ensor s derved here. The resuls obaned below generalze a well-known approxmaon for elecromagnec energy n dspersve meda. As s famlar from he exbook reamens [7, 8], he Poynng heorem ha follows from Maxwell s equaons does no drecly lead o an expresson for he elecromagnec energy densy n he case of lossless, dspersve meda. The sandard procedure [7, 8] s o make a resrcon o quasmonochromac felds, n whch case one can exrac an approxmae expresson for he me-averaged energy densy, due o Brlloun (Eq. (1) below). For general felds one s reduced o makng a formal me negraon of he Poynng

2 heorem [9]. Bu as noed above, n he dspersve case one canno derve Noeher s heorem from a sysemac manpulaon of he dynamcal equaons of he heory, as s done n he Poynng heorem; nsead, one mus sar wh he acon ha underles he equaons. Noeher s heorem guaranees he exsence of an exac energy densy for arbrary felds n lossless dspersve meda and hs energy densy s derved here. In he case of elecromagnec momenum densy, he dspersonless resul s he famlar Mnkowsk expresson [ 6] D B. We noe agan ha hs saemen refers o he conserved elecromagnec energy-momenum ensor: he Mnkowsk momenum densy D B unambguously gves he conserved elecromagnec momenum of lgh propagang n a non-dspersve, lossless, homogeneous medum [ 6]. I s shown below ha n dspersve meda he momenum densy s more complcaed, and s no gven by Mnkowsk s formula. II. THE ACTION The resuls follow sraghforwardly once an acon prncple s wren for elecromagnesm n dspersve, lossless meda. For hs purpose he dynamcal varables mus be aken o be he scalar poenal φ and vecor poenal A, defned by E = φ A, B = A. (1) The relaonshp beween felds n he me and frequency domans (he laer denoed by a lde) s, for he example of he elecrc feld, E(r, ) = 1 π = 1 π 0 ) dω (Ẽ(r, ω)e ω + c.c dω Ẽ(r, ω)e ω, where he realy of E(r, ) mples Ẽ(r, ω) = Ẽ (r, ω). In he frequency doman we have () D = ε 0 ε(r, ω)ẽ, H = κ0 κ(r, ω) B, (3) where ε(r, ω) s he relave permvy of he (n general nhomogeneous) medum, κ 0 = µ 1 0, and he relave permeably s µ(r, ω) = κ(r, ω) 1. As he medum can be assumed o be lossless n he frequency range of neres, ε(r, ω) and κ(r, ω) are real and even funcons of ω [8], and so hey have he seres expansons ε(r, ω) = ε n (r) ω n, κ(r, ω) = κ n (r) ω n, () n=0 n=0 where hese equaons defne he coeffcens ε n (r) and κ n (r). In pracce, he seres () wll represen a f o he dsperson daa of he maeral for he frequency range of neres; we ake hese seres o be nfne, bu he subsequen resuls also hold when hey are fne seres. Equaons () () show ha n he me doman he Maxwell equaons (wh no free charges or currens) are ε 0 [ε(r, )E] = 0, (5) κ 0 [κ(r, )B] = ε 0 ε(r, ) E, (6) whle he oher wo Maxwell equaons are auomacally sasfed because of (1). Noe ha seres expansons of he form () are sandard n reang dsperson n frequency ranges where losses are neglgble, for example n lnear and nonlnear fber opcs [11]. Moreover he ransformaon from he frequency o he me doman apparen n (5) (6), and he reverse ransformaon, are a sandard par of numercal soluon procedures for he propagaon of wave packes wh specra whn he frequency range where expansons () are vald, for example n he spl-sep mehod [11, 13]. The acon S[φ, A] for he poenals φ and A ha gves he dynamcal equaons (5) (6) s S = d x κ 0 { } 1 c E [ε(r, )E] B [κ(r, )B]. (7) Varaon of φ n (7) gves he Maxwell equaon (5), whle varaon of A gves (6). The fac ha only even-order me dervaves occur n (5) (6) s essenal o beng able o wre an acon prncple for hese equaons. For dsspave meda, erms wh odd-order me dervaves occur n (5) (6), arsng from he magnary pars of ε(r, ω) and κ(r, ω) whch are odd funcons of ω [8]; bu he acon (7) would no generae hese erms. In he acon (7), erms of he form E n+1 E, for example, would no conrbue o he dynamcal equaons because her varaon gves zero afer negraons by pars: δ(e n+1 E) = (δe) n+1 E + E n+1 δe, and he second erm becomes mnus he frs upon negraons by pars. I mus of course be mpossble o wre an acon n he dsspave case ha s a funconal only of φ and A, snce hs would mply he exsence of a conserved elecromagnec energy n lossy meda (see below). III. ENERGY The acon (7) s nvaran under acve me ranslaons of he dynamcal felds φ and A, and hs nvarance mples, hrough Noeher s heorem, he exsence of a conserved quany, he energy. The exracon of he conservaon law from he acon s a sandard echnque of feld heory (see [1], for example). Even f one s unfamlar wh Noeher s heorem, one can of course verfy usng he feld equaons (5) (6) ha he resulng conservaon law does n fac hold. The heorem [1] shows ha f we make an acve nfnesmal me ranslaon φ(r, ) φ(r, + ζ(r, )), A(r, ) A(r, + ζ(r, )), bu ake he ranslaon ζ(r, ) o vary n space and me, hen he change n he acon can be reduced o he form δs = d x (ρ ζ + S ζ), (8)

3 3 where ρ s he energy densy and S s he energy flux, obeyng he conservaon law ρ + S = 0. (9) The varaon δs wll clearly be lnear n he nfnesmal ranslaon ζ(r, ), bu varous numbers of negraons by pars are requred o acheve he form (8) (surface erms produced by he negraons are o be dropped). Use mus be made of he followng denes, whch hold for arbrary funcons Y and Z, up o surface erms: [ d x Y ε(r, )(ζz) = d x ζy ε(r, )Z ( 1) n+m ε n (r) m 1 Y n m Z ζ [ d x Y ε(r, ) (ζz) = d x ζy ε(r, ) Z + n=0 m=0 ( 1) n+m ε n (r) m Y n m Z ζ ] ], (10). (11) These denes were found by frs workng wh small values of n, where all he erms can be checked by hand, and hen verfyng he general expressons (10) (11) usng Mahemaca. Ideny (10) wh ε(r, ) replaced by κ(r, ) mus also be used. In hs way one aans he form (8) wh ρ = κ 0 1 c + { 1 c ( φ A) [ε(r, )E] + B [κ(r, )B] ( 1) n+m ε n (r) m 1 E n m+1 E ( 1) n+m κ n (r) m 1 B n m+1 B }, (1) S = ε 0 φ ε(r, )E κ 0 A [κ(r, )B]. (13) I s sraghforward o verfy ha (1) and (13) obey he conservaon law (9) when he felds sasfy Maxwell s equaons (5) (6). As n he case of dspersonless meda, and ndeed vacuum, he energy densy (1) and flux (13) ha drecly emerge from Noeher s heorem are no gauge nvaran [7, 1]. They are however equvalen o gauge-nvaran quanes because hey fal o be gauge nvaran up o erms ha dencally sasfy he conservaon law (9). Specfcally, he quanes f dencally sasfy := ε 0 φ ε(r, )E =: f, (1) f j := κ 0 φ κ(r, )( A j j A ), (15) f + ( f + j f j ) = 0. (16) Comparng (16) wh (9), we see ha f f s added o ρ, and f + j f j s added o S, hen he conservaon law (9) wll sll hold. Moreover, wh use of Maxwell s equaons (5) (6), he energy densy and flux ha resul from hese addons are gauge nvaran and are gven by ρ = κ 0 1 c + { 1 c E [ε(r, )E] + B [κ(r, )B] ( 1) n+m ε n (r) m 1 E n m+1 E ( 1) n+m κ n (r) m 1 B n m+1 B }, (17) S =κ 0 E [κ(r, )B]. (18) The conservaon law (9) for he fnal expressons (17) and (18) can also be verfed usng Maxwell s equaons (5) (6). In (18) we see ha he energy flux s gven by he Poynng vecor E H n he me doman, he obvous generalzaon of he non-dspersve resul. The energy densy (17), however, has no smple relaon o he non-dspersve resul, whch only conans he frs wo erms n (17). I s a smple maer o express (17) and (18) n he frequency doman. We can use he exac energy densy (17) o derve he sandard exbook resul [7, 8] for he me-averaged energy densy of a monochromac wave n a dspersve medum. Such a wave has an elecrc feld E(r, ) = 1 ( E0 (r)e ω0 + c.c ), (19) and he B feld s of he same form. When he monochromac E and B felds are subsued no (17) and a me average s aken, all -dependen erms vansh. The seres n (17) conanng he E feld s hen ε 0 ε n (r)ω n 0 E 0 = ε 0 n=1 n ε n (r)ω n 0 E 0 = ε 0 ω dε(r, ω 0 ) 0 E 0. (0) The seres n (17) conanng he B feld undergoes a smlar smplfcaon and he me-averaged monochromac energy densy ρ mono can be wren ρ mono = ε 0 d[ω 0 ε(r, ω 0 )] E 0 + µ 0 d[ω 0 µ(r, ω 0 )] H 0, (1) whch s Brlloun s formula [7, 8]. For homogeneous meda he monochromac wave s a plane wave and s easy o show, by dvdng he me-averaged energy flow S by ρ, ha he elecromagnec energy of he plane wave moves hrough he medum a he group velocy c/[d(ω 0 ε(ω0 )µ(ω 0 ))/ ]. I s of course expeced ha he energy of a monochromac wave n a lossless, homogeneous medum should move a he group velocy, bu here we have derved hs fac for elecromagnesm n lossless meda wh arbrary dsperson.

4 E permvy s neglgble; he real par s Ω(ω ω ε(ω) =1 0) (ω ω0 ) + γ ω =1 + C n ω n, ω < ω 0, () C n = Ω n=0 ω0 ωn m=0 n ( ) ( ) ( 1) m n + m γ m. (3) m ω 0 Ρ FIG. 1: The emporal profle of he elecrc feld (op) and energy densy (boom) of a gaussan pulse n a medum wh relave permvy () and µ = 1. Noe ha he energy densy s no zero a all he local mnma; hs feaure s due solely o he complcaed elecrc-feld seres n (17). In he case of felds ha are nearly monochromac, Brlloun s expresson (1) can serve as an approxmaon of he me-averaged energy densy, bu only f ε(r, ω) and µ(r, ω) do no vary sgnfcanly over he range of frequency componens n he felds [7]. As was recenly poned ou n [10], Brlloun s formula s also exac for me-averaged felds whose frequency componens are uncorrelaed, he mporan example beng hermal radaon. Bu he exac energy densy for all felds n dspersve, lossless meda s gven by (17). As menoned n Secon II, he seres represenaons () on whch our resuls are based wll n pracce be numercal fs o he measured dsperson daa of he maeral n queson, n he frequency range of neres. Clearly, a fne number of erms n he seres wll be suffcen for an accurae reamen of he dsperson, and he forgong resuls hold n hs case where all seres ermnae. As a heory exercse, however, we now consder an example of an nfne seres expanson of he form () for he permvy, and verfy ha he expresson (17) for he energy densy s a well defned quany. We ake he sandard, damped harmonc oscllaor model of a homogeneous permvy [7]. A frequences well below he resonan frequency ω 0 of he oscllaor, he magnary par of he relave The seres expanson () converges for ω < ω 0 and we consder a lgh beam n hs maeral wh a specrum lyng far enough below ω 0 for losses o be neglgble. The evoluon of he lgh beam n hs dspersve medum can of course only be compued numercally, bu we consder an nal, npu gaussan pulse wh cenral frequency ω 0 / and calculae s nal energy dsrbuon. As s clear from (17), he effec of dsperson on he nal energy densy wll be seen n he emporal profle of he pulse, raher han he spaal profle, so we consder a fxed pon r = 0, whch can be vewed as a boundary hrough whch he pulse eners he medum. We ake uns wh c = 1 and choose he values ω 0 = 10, Ω = 30 and γ = 1 n he relave permvy (). The npu elecrc feld of he gaussan pulse cenered on ω = ω 0 / a he boundary r = 0 s aken as E() r=0 = e / cos(ω 0 /) () and s ploed n Fg. 1. We can numercally compue he npu magnec feld B() r=0 by ransformng o he frequency doman and use of he dsperson relaon k = εω/c and he Maxwell equaon ω B = Ẽ. To avod nessenal complcaons we gnore he ransverse spaal profle of he pulse and consder a one-dmensonal propagaon (an acual one-dmensonal propagaon n a wavegude such as an opcal fber wll nvolve an effecve dsperson dfferen from (), bu our consderaons here are purely for demonsraon purposes). Usng E() r=0 and B() r=0 we compue he npu emporal profle ρ() r=0 of he energy densy (17) (wh µ = 1/κ = 1) and he resul s shown n Fg. 1. Noe ha no cycle-averagng has been performed; he exac energy densy (17) s no a cycle-averaged quany. As well as he quanave change n he energy densy caused by he dsperson, here s a qualave dfference, vsble n Fg. 1, compared o he same pulse n a non-dspersve medum. The local mnma of he energy densy n Fg. 1 occur a nodes of he elecrc (and magnec) feld, bu s clearly seen ha he energy densy does no drop o zero a many of hese local mnma. Ths s purely an effec of he dsperson; n a non-dspersve medum hese local mnma are zeros of he energy densy. In fac hs effec of dsperson s due solely o he complcaed elecrc-feld seres n (17). The famlar E D erm n (17) does no conrbue o hs feaure, whch can only be seen usng he exac resul (17).

5 5 IV. MOMENTUM Turnng o he elecromagnec momenum, we noe ha momenum wll be conserved only f he medum s homogeneous, so ha Noeher s heorem apples o spaal ranslaons. I s nsrucve however o rean ε and κ as funcons of poson and hereby see how he momenum conservaon law fals o hold n he nhomogeneous case. To exrac he conservaon law assocaed wh spaal ranslaon nvarance we make he acve nfnesmal ranslaon φ(r, ) φ(r + η(r, ), ), A(r, ) A(r + η(r, ), ) n he acon (7). Noeher s heorem [1] shows ha when he resulng change n he acon s wren n he form ( δs = d x P η + σ j j η ), (5) hen (n homogeneous meda) P and σ j obey he conservaon law P + j σ j = 0. (6) In hs way we fnd he elecromagnec momenum densy P and sress ensor σ j. We wll carry ou hs procedure wh ε and κ varyng n space; hen (6) wll fal o hold because exra erms wll appear n hs equaon ha conan spaal dervaves of ε and κ. Agan, use mus be made of he denes (10) and (11) o acheve he form (5) and he resul s P = κ [ 0 c A j ε(r, )E j σ j + 1 c ( 1) n+m ε n (r) m 1 E j n m E j ] ( 1) n+m κ n (r) m 1 B j n m B j, (7) = L δ j + ε 0 φ ε(r, )E j + κ 0 A k κ(r, )[ j A k k A j ], (8) where L n (8) s he Lagrangan densy,.e. he negrand n he acon (7). To oban gauge-nvaran expressons for he momenum densy and sress ensor, we noe ha he quanes f j dencally sasfy := ε 0 A ε(r, )E j =: f j, (9) f k j := κ 0 A κ(r, )( j A k k A j ), (30) j f j + j ( f j + k f k j ) = 0. (31) Thus, addon of j f j o P, and of f j + kf k j o σ j, does no affec he momenum conservaon law (6). Afer hese addons and use of Maxwell s equaons (5) (6), he momenum densy and sress ensor are gauge nvaran and are gven by P = ε 0 ɛ jk [ε(r, )E j ]B k σ j + ε 0 κ 0 ( 1) n+m ε n (r) m 1 E j n m E j ( 1) n+m κ n (r) m 1 B j n m B j, = ε 0 E ε(r, )E j κ 0 [κ(r, )B ]B j + 1 δ j [ ε0 E k ε(r, )E k + κ 0 B k κ(r, )B k], (3) (33) where ɛ jk s he (compleely ansymmerc) Lev-Cva ensor. I s sraghforward o verfy ha, when he felds sasfy Maxwell s equaons (5) (6), he momenum densy (3) and sress ensor (33) sasfy P + j σ j = ε 0 E j[ ε(r, )]E j κ 0 B j[ κ(r, )]B j, (3) so ha he momenum conservaon law (6) ndeed holds for homogeneous meda. The dervaon of (3) n he nondspersve case s famlar from he exbooks [15]. Noe ha he elecromagnec sress ensor (33) s here defned so ha s he spaal par of he elecromagnec energymomenum ensor; a wdespread convenon for he sress ensor dffers from hs by a mnus sgn [7]. The sress ensor (33) s he obvous generalzaon o he dspersve case of he non-dspersve resul; he fac ha hs expresson descrbes he sress ensor n dspersve meda s deduced n [8] from compleely dfferen consderaons. For he monochromac wave (19), he seres n he momenum densy (3) conanng he E feld s, afer a me averagng, ε 0 8 ( ) ε n ω0 n 1 E 0j E j 0 E 0j E j 0 = ε 0 dε(ω 0 ) Im(E 0j E j 0 ). (35) There s a smlar resul for he seres n (3) conanng he magnec feld, and he complee resul for he meaveraged monochromac momenum densy P n a homogeneous medum s P = ε 0 ε(ω 0)Re(E 0 B 0) ε 0 + κ 0 dε(ω 0 ) Im(E 0 E0 ) dκ(ω 0 ) Im(B 0 B0 ). (36) The Mnkowsk momenum densy D B would gve only he frs erm n (36). Wh an addonal resrcon o a sngle plane wave of frequency ω 0 n he homogeneous medum,.e. E 0 (r) = E 0 e k r, B 0 (r) = 1 k E 0 (r), (37) ω 0 k = ω 0 ε(ω0 )µ(ω 0 ), (38) c

6 6 P FIG. : The emporal profle of he momenum densy of a gaussan pulse n a dspersve medum. The pulse and dsperson are he same as hose descrbed a he end of Secon III and used for he energy densy n Fg. 1. Noe ha he momenum densy s no zero a local mnma; hs feaure s due solely o he complcaed elecrcfeld seres n (3). he expresson (36) reduces o P = ε 0 c ε(ω 0)n g (ω 0 )ˆk E 0, (39) n g (ω 0 ) = n p (ω 0 ) + ω 0 dn p (ω 0 ), (0) where n g (ω 0 ) s he group ndex, n p (ω 0 ) = ε(ω 0 )µ(ω 0 ) s he phase ndex and ˆk = k/ k. Usng he quanum prescrpon (for µ = 1) n [16] corresponds o akng E 0 ω 0 /[ε 0 n p (ω 0 )n g (ω 0 )] for a sngle phoon; hs gves from (39) a sngle phoon momenum ha s n p (ω 0 ) mes he vacuum value, n agreemen wh he canoncal momenum n [6]. We noe ha he Mnkowsk expresson (he frs erm n (36)) would gve a correcon facor of n p(ω 0 )/n g (ω 0 ) for he sngle phoon [16], whereas he facor n p (ω 0 ) has been obaned here from he general momenum densy (3). A he end of Secon III we showed an example of he elecromagnec energy densy for a gaussan pulse n a dspersve medum. For he same pulse and medum we use (3) o compue he npu emporal profle of he momenum densy, where he fnal seres n he magnec feld n (3) s absen for hs example. The spaal dervave of he elecrc feld n (3) s compued by ransformng o he frequency doman and usng he dsperson relaon k = εω/c (as dscussed n Secon III, we assume one-dmensonal propagaon). The profle for he momenum densy s shown n Fg. and s seen o have he same shape as he energy densy n Fg. 1. Smlar o wha was found for he energy densy, he fac ha he momenum densy does no fall o zero a local mnma s due solely o he complcaed elecrc-feld seres n (3); hs feaure can only be seen wh he exac expresson (3) and s absen from he Mnkowsk conrbuon D B n (3). Elecromagnec angular momenum n dspersve meda can be derved n a smlar manner, hough hs wll no be done here. The angular momenum densy s no o be consruced from he energy-momenum ensor bu raher derved from he nvarance of he acon (7) under acve roaons (for roaonally symmerc meda). V. CONCLUSIONS We have derved he exac conserved energy-momenum ensor of lgh propagang n lossless, dspersve, homogeneous meda. Elecromagnec energy s also conserved when he medum s nhomogeneous and hs case has been ncluded. The energy flux and sress ensor n dspersve meda have he same general form as n he non-dspersve case, wh he permvy and permeably becomng dervave operaors n he me doman. In conras, he energy densy and momenum densy for arbrary felds have no smple relaon o he non-dspersve resuls. For me-averaged monochromac waves he Brlloun energy densy s recovered, and he momenum densy also akes a smple form (equaon (36)). As sressed n he Inroducon, he resuls do no address he queson of momenum balance n lgh-maer neracons, whch requres an analyss of momenum ransfer beween lgh and maer [ 6]. The conserved energy-momenum ensor of lgh n lossless meda s a well-defned quany wh s own expermenal sgnfcance [1]; has been derved here n complee generaly whn he framework of macroscopc elecromagnesm. Acknowledgemens I hank Olver Allanson for ponng ou an error n he orgnal monochromac momenum resul. Ths research s suppored by he Scosh Governmen and he Royal Socey of Ednburgh. [1] T. G. Phlbn, Phys. Rev. A 83, (011). Erraum: Phys. Rev. A 85, 05990(E) (01). [] I. Brevk, Phys. Rep. 5, 133 (1979). [3] U. Leonhard, Naure (London), 83 (006). [] R. N. C. Pfefer, T. A. Nemnen, N. R. Heckenberg, and H. Rubnszen-Dunlop, Rev. Mod. Phys. 79, 1197 (007). [5] S. M. Barne and R. Loudon, Phl. Trans. R. Soc. A 368, 97 (010). [6] S. M. Barne, Phys. Rev. Le. 10, (010). [7] J. D. Jackson, Classcal Elecrodynamcs, 3rd ed. (Wley, New York, 1999). [8] L. D. Landau, E. M. Lfshz, and L. P. Paevsk, Elecrodynamcs of Connuous Meda, nd ed. (Buerworh-Henemann, Oxford, 198). [9] Y. S. Barash and V. L. Gnsburg, Sov. Phys. Usp. 19, 63 (1976).

7 7 [10] F. S. S. Rosa, D. A. R. Dalv, and P. W. Mlonn, Phys. Rev. A 81, (010). [11] G. Agrawal, Nonlnear Fber Opcs, h ed. (Academc Press, San Dego, 006). [1] T. G. Phlbn, C. Kuklewcz, S. Roberson, S. Hll, F. Köng, and U. Leonhard, Scence (008). [13] W. H. Reeves, D. V. Skryabn, F. Bancalana, J. C. Kngh, P. S. J. Russell, F. G. Omeneo, A. Efmov, and A. J. Taylor, Naure 511 (003). [1] S. Wenberg, The Quanum Theory of Felds, Volume I (Cambrdge Unversy Press, Cambrdge, 1995). [15] J. A. Sraon, Elecromagnec Theory, (McGraw-Hll, New York, 191). [16] J. C. Garrson and R. Y. Chao, Phys. Rev. A 70, (00).