PHYSICAL REVIEW E VOLUME 55. NUMBER 1 JANUA Propagation of electromagnetic energy and momentum through an absorbing dielectric R. Loudon and L. allen Department of Physics, Essex University, Colchester CO4 3S0, England D. F Nelson Department of Physics, Worcester Polytechnic Institute, Worcester, Massachusetts 01609 We calculate the energy and momentum densities and currents associated with electromagnetic wave propa- gation through an absorbing and dispersive diatomic dielectric, which is modeled by a single- resonance Lor- entz oscillator. The relative and center-of-mass coordinates of the dielectric sublattices and the electromagnetic field vectors are treated as dynamical variables, while the dielectric loss is modeled by a phenomenological damping force. The characteristics of the energy propagation agree with previous work, including the form of the energy velocity. The treatment of momentum propagation extends previous work to lossy media, and it is found that the damping plays an important role in the transfer of momentum from the electromagnetic field to the center of mass of the dielectric. We discuss the significances of the momentum, the pseudomomentum, and their sum, the wave momentum. For each of these quantities we derive the density, the current density, and the appropriate conservation or continuity equation. The general expressions are illustrated by applications to a steady-state monochromatic wave and to an excitation in the form of a localized Gaussian pulse. The velocities associated with propagation of the various kinds of momentum are derived and discussed PACS number(s):4225Bs,03.40.-t,03.50.De,41.20.Jb . INTRODUCTION where the momentum density is a vector quantity propor- tional to the energy current, G=S/c, and the momentum The nature of electromagnetic energy and the characteris- current density is a second rank tensor, or 3X3 matrix, re- tics of its propagation through dielectric media have been lated to the Maxwell stress tensor [1, 2, 9). For the electro- studied since the early years of electromagnetic theory. For magnetic momentum in material media, it is necessary to propagation through the simplest kind of linear, isotropic, take account of contributions from both the electromagnetic and homogeneous medium, the energy density W and energy field and the dielectric medium. The momentum current in a current density, or Poynting vector S are routinely treated in lossless dielectric was obtained by this approach as a modi standard texts [1, 2]. The forms of these energy densities and fied form of the Maxwell stress tensor. In addition, the na- their conservation law have also been evaluated for much tures of the momentumlike quantities that have been defined more general dielectric media [3]. For propagation through for the coupled system of electromagnetic field and dielectric absorbing or scattering materials, the classic treatment of material, including the densities proposed by Abraham and electromagnetic wave propagation, and particularly the iden- Minkowski, were identified [6] tification of the several distinct velocities that are associated The controversy has always revolved around a linear light with an optical pulse, was provided by Sommerfeld and Bril- wave for which deformation of the dielectric medium is ir- louin [4]. The detailed theory for lossy dielectrics is quite relevant, but a key ingredient of its recent resolution is the complicated, but the main features of the energy density and inclusion of deformation of the medium. This necessitates current, and of energy propagation, are correctly predicted by the use of both spatial(Eulerian) and material (Lagrangian) a simple calculation [5], based on the standard model of coordinates, and it allows the deduction of conservation laws electromagnetic waves in a lorentzian dielectric with a from Noether's theorem Thus the momentum conservation single resonance. The essential feature of this theory is the law follows from invariance to displacements of the spatial inclusion of contributions to the total energy density W and coordinates(homogeneity of free space), and the pseudomo- energy current density S of the optical excitation from both mentum conservation law follows from invariance to dis the electromagnetic field and the dielectric medium placements of the material coordinates(homogeneity of the interesting to determine whether there is an analo- material medium). This approach [6] found the electromag gous theory for electromagnetic momentum propagation in a netic momentum density Gm to be EoEXB, close to, but in lossy medium, and this is the primary purpose of the present general different from, the Abraham form EoloEXH. It also paper. Such an inquiry is particularly topical because recent found the pseudomomentum density Gsm to be PXB plus a rk [6]on di but lossless dielectrics has G and a resolution of the long-standing Minkowski-Abraham con- pseudomomentum densities, which we call the wave momen- troversy concerning the correct expressions for the densities tum, is the generalization of the Minkowski momentum of electromagnetic momentum and current, denoted here by DXB to include dispersion. However, the Minkowski mo- G and T, respectively(see [7, 8]for reviews). These densities mentum was proposed as being the ordinary momentum are well understood for electromagnetic fields in free space, while this derivation shows instead that it is the sum of or- 1063-651X/97/55(1)/1071(l5)10.00 c 1997 The American Physical Society
Propagation of electromagnetic energy and momentum through an absorbing dielectric R. Loudon and L. Allen Department of Physics, Essex University, Colchester CO4 3SQ, England D. F. Nelson Department of Physics, Worcester Polytechnic Institute, Worcester, Massachusetts 01609 ~Received 22 August 1996! We calculate the energy and momentum densities and currents associated with electromagnetic wave propagation through an absorbing and dispersive diatomic dielectric, which is modeled by a single-resonance Lorentz oscillator. The relative and center-of-mass coordinates of the dielectric sublattices and the electromagnetic field vectors are treated as dynamical variables, while the dielectric loss is modeled by a phenomenological damping force. The characteristics of the energy propagation agree with previous work, including the form of the energy velocity. The treatment of momentum propagation extends previous work to lossy media, and it is found that the damping plays an important role in the transfer of momentum from the electromagnetic field to the center of mass of the dielectric. We discuss the significances of the momentum, the pseudomomentum, and their sum, the wave momentum. For each of these quantities we derive the density, the current density, and the appropriate conservation or continuity equation. The general expressions are illustrated by applications to a steady-state monochromatic wave and to an excitation in the form of a localized Gaussian pulse. The velocities associated with propagation of the various kinds of momentum are derived and discussed. @S1063-651X~97!04901-5# PACS number~s!: 42.25.Bs, 03.40.2t, 03.50.De, 41.20.Jb I. INTRODUCTION The nature of electromagnetic energy and the characteristics of its propagation through dielectric media have been studied since the early years of electromagnetic theory. For propagation through the simplest kind of linear, isotropic, and homogeneous medium, the energy density W and energy current density, or Poynting vector S are routinely treated in standard texts @1,2#. The forms of these energy densities and their conservation law have also been evaluated for much more general dielectric media @3#. For propagation through absorbing or scattering materials, the classic treatment of electromagnetic wave propagation, and particularly the identification of the several distinct velocities that are associated with an optical pulse, was provided by Sommerfeld and Brillouin @4#. The detailed theory for lossy dielectrics is quite complicated, but the main features of the energy density and current, and of energy propagation, are correctly predicted by a simple calculation @5#, based on the standard model of electromagnetic waves in a Lorentzian dielectric with a single resonance. The essential feature of this theory is the inclusion of contributions to the total energy density W and energy current density S of the optical excitation from both the electromagnetic field and the dielectric medium. It is interesting to determine whether there is an analogous theory for electromagnetic momentum propagation in a lossy medium, and this is the primary purpose of the present paper. Such an inquiry is particularly topical because recent work @6# on dispersive, but lossless, dielectrics has proposed a resolution of the long-standing Minkowski-Abraham controversy concerning the correct expressions for the densities of electromagnetic momentum and current, denoted here by G and T, respectively ~see @7,8# for reviews!. These densities are well understood for electromagnetic fields in free space, where the momentum density is a vector quantity proportional to the energy current, G5S/c2 , and the momentum current density is a second rank tensor, or 333 matrix, related to the Maxwell stress tensor @1,2,9#. For the electromagnetic momentum in material media, it is necessary to take account of contributions from both the electromagnetic field and the dielectric medium. The momentum current in a lossless dielectric was obtained by this approach as a modi- fied form of the Maxwell stress tensor. In addition, the natures of the momentumlike quantities that have been defined for the coupled system of electromagnetic field and dielectric material, including the densities proposed by Abraham and Minkowski, were identified @6#. The controversy has always revolved around a linear light wave for which deformation of the dielectric medium is irrelevant, but a key ingredient of its recent resolution is the inclusion of deformation of the medium. This necessitates the use of both spatial ~Eulerian! and material ~Lagrangian! coordinates, and it allows the deduction of conservation laws from Noether’s theorem. Thus the momentum conservation law follows from invariance to displacements of the spatial coordinates ~homogeneity of free space!, and the pseudomomentum conservation law follows from invariance to displacements of the material coordinates ~homogeneity of the material medium!. This approach @6# found the electromagnetic momentum density Gm to be «0E3B, close to, but in general different from, the Abraham form «0m0E3H. It also found the pseudomomentum density Gpsm to be P3B plus a dispersive term. Thus the sum G of the momentum and pseudomomentum densities, which we call the wave momentum, is the generalization of the Minkowski momentum D3B to include dispersion. However, the Minkowski momentum was proposed as being the ordinary momentum, while this derivation shows instead that it is the sum of orPHYSICAL REVIEW E VOLUME 55, NUMBER 1 JANUARY 1997 1063-651X/97/55~1!/1071~15!/$10.00 1071 © 1997 The American Physical Society 55
R LOUDON. L. ALLEN. AND D. F. NELSON dinary momentum and pseudomomentum. The name"wave nomentum'' was introduced for this reason v×E While the inclusion of material deformation has played ar essential role in the clarification of what is momentum and what is pseudomomentum, it also acts as a barrier to simple V×B=j+ (24) physical understanding. The aim of the present paper is thus to find simplified versions of the conservation laws for mo- mentum and pseudomomentum, even after generalization of where the fields are functions of position and time previous work to include loss. This is achieved by the sim- E=E(r, t), and so on. The bound charge and current densi with cubic isotropy, essentially the single-resonance Lorentz time; they can be expressed in terms of the dielectric po i plification of the dielectric to a nonmagnetic diatomic crystal ties, p and j, respectively, are also functions of position model. The ions are assumed to be coupled to the electro- Ization P as magnetic field only by an electric-dipole interaction Before proceeding to the main calculations, we present a V. P in Sec. Il. An improved version of previous calculations of and the energy continuity equations and the velocity of energy propagation [5] leads to essentially the same results as be- (26) fore, but our method facilitates parallel discussions of the momentum propagation characteristics. It is found, however that the propagation of momentum involves both the center The electric displacement is defined in the usual way, of mass and relative coordinates of the diatomic dielectric whose proper treatment requires a Lagrangian formulation D=EoE+P current obey a conservation law when the center-of-mass and the magnetic field is given by momentum is included, but that the pseudomomentum, and hence the wave momentum, suffer dissipation on account of (28) he dielectric loss. The various electromagnetic densities de- Note that in the view implicit in these equations, E and b are rived in Secs. Il and Ill are evaluated for a steady-state the fundamental electromagnetic fields, while P describes the monochromatic wave in Sec. IV. where the velocities of propagation of energy and wave momentum are derived, and response of the matter, and Eqs. (2.7)and(2. 8)are constitu- for an optical pulse in Sec. V. The results are discussed in tive equations for D and H We consider a polar, diatomic, cubic, crystal lattice in which the relative spatial displacement field of the two ions in the unit cell is denoted s=s(r, t). The long-wavelength Il. SIMPLE THEORY OF ELECTROMAGNETIC ENERGY optic modes of vibration have a basic threefold degeneracy AND MOMENTUM PROPAGATION which is lifted by the long-range electrical forces to form a wofold-degenerate transverse mode and a nondegenerate The present section is devoted to a derivation of some longitudinal mode[10]. Then, if the frequency of the trans basic results for electromagnetic fields in a dielectric material verse mode is denoted or and its damping rate is denoted T, reated in the Lorentz model. We present a simple derivation the standard form of the lorentz equation for the ith Carte- of the equations that describe the propagation of energy, and sian component of the internal coordinate of the ionic motion show that the corresponding description of momentum IS propagation cannot be obtained by so simple a theory. The detailed derivations of the equations that describe momen- ms,+mrs+mo tum propagation and the identification of the different char- acters of the momentumlike contributions are given careful Here m is the reduced mass density of the two ions, of consideration in Sec. Ill masses Mi and M2, in the primitive unit cell of volume n2 A. Basic equations M1M2 1=(M1+M2) (2.10) The fundamental energy and momentum properties electromagnetic fields in matter are governed by maxwells nd the chars equations and by the equations of motion for the matter. we is given by ity s associated with the internal motion consider a nonmagnetic dielectric material that has no free charges or currents. The Maxwell-Lorentz forms of the equa s=e/① (2.11) tions in conventional notation and Systeme International (SI nits are then where e and -e are the charges on the two kinds of ion. The VE=p/Eo, (2. 1 polarization is expressed in terms of the internal coordinate V.B=0, (2.12)
dinary momentum and pseudomomentum. The name ‘‘wave momentum’’ was introduced for this reason. While the inclusion of material deformation has played an essential role in the clarification of what is momentum and what is pseudomomentum, it also acts as a barrier to simple physical understanding. The aim of the present paper is thus to find simplified versions of the conservation laws for momentum and pseudomomentum, even after generalization of previous work to include loss. This is achieved by the simplification of the dielectric to a nonmagnetic diatomic crystal with cubic isotropy, essentially the single-resonance Lorentz model. The ions are assumed to be coupled to the electromagnetic field only by an electric-dipole interaction. Before proceeding to the main calculations, we present a simplified discussion of energy and momentum propagation in Sec. II. An improved version of previous calculations of the energy continuity equations and the velocity of energy propagation @5# leads to essentially the same results as before, but our method facilitates parallel discussions of the momentum propagation characteristics. It is found, however, that the propagation of momentum involves both the center of mass and relative coordinates of the diatomic dielectric, whose proper treatment requires a Lagrangian formulation. Thus it is shown in Sec. III that the momentum density and current obey a conservation law when the center-of-mass momentum is included, but that the pseudomomentum, and hence the wave momentum, suffer dissipation on account of the dielectric loss. The various electromagnetic densities derived in Secs. II and III are evaluated for a steady-state monochromatic wave in Sec. IV, where the velocities of propagation of energy and wave momentum are derived, and for an optical pulse in Sec. V. The results are discussed in Sec. VI. II. SIMPLE THEORY OF ELECTROMAGNETIC ENERGY AND MOMENTUM PROPAGATION The present section is devoted to a derivation of some basic results for electromagnetic fields in a dielectric material treated in the Lorentz model. We present a simple derivation of the equations that describe the propagation of energy, and show that the corresponding description of momentum propagation cannot be obtained by so simple a theory. The detailed derivations of the equations that describe momentum propagation and the identification of the different characters of the momentumlike contributions are given careful consideration in Sec. III. A. Basic equations The fundamental energy and momentum properties of electromagnetic fields in matter are governed by Maxwell’s equations and by the equations of motion for the matter. We consider a nonmagnetic dielectric material that has no free charges or currents. The Maxwell-Lorentz forms of the equations in conventional notation and Syste`me International ~SI! units are then “•E5r/«0 , ~2.1! “•B50, ~2.2! “3E52 ]B ]t , ~2.3! 1 m0 “3B5j1«0 ]E ]t , ~2.4! where the fields are functions of position and time, E[E~r,t!, and so on. The bound charge and current densities, r and j, respectively, are also functions of position and time; they can be expressed in terms of the dielectric polarization P as r52“•P ~2.5! and j5 ]P ]t . ~2.6! The electric displacement is defined in the usual way, D5«0E1P, ~2.7! and the magnetic field is given by H5B/m0 . ~2.8! Note that in the view implicit in these equations, E and B are the fundamental electromagnetic fields, while P describes the response of the matter, and Eqs. ~2.7! and ~2.8! are constitutive equations for D and H. We consider a polar, diatomic, cubic, crystal lattice in which the relative spatial displacement field of the two ions in the unit cell is denoted s[s~r,t!. The long-wavelength optic modes of vibration have a basic threefold degeneracy which is lifted by the long-range electrical forces to form a twofold-degenerate transverse mode and a nondegenerate longitudinal mode @10#. Then, if the frequency of the transverse mode is denoted vT and its damping rate is denoted G, the standard form of the Lorentz equation for the ith Cartesian component of the internal coordinate of the ionic motion is ms¨i1mGs˙i1mvT 2 si5§Ei . ~2.9! Here m is the reduced mass density of the two ions, of masses M1 and M2 , in the primitive unit cell of volume V, m5 M1M2 V~M11M2! ~2.10! and the charge density § associated with the internal motion is given by §5e/V, ~2.11! where e and 2e are the charges on the two kinds of ion. The polarization is expressed in terms of the internal coordinate by P5§s. ~2.12! 1072 R. LOUDON, L. ALLEN, AND D. F. NELSON 55
PROPAGATION OF ELECTROMAGNETIC ENERGY AN 1073 B. Energy propagation B. B The flow of electromagnetic energy through the dielectric (Tem)ji=-EoE E,+iEoE-Sjn or is determined by the energy current density, or Poynting vec- This quantity is usually identified as the negative of the Max (2.13) well stress tensor [1, 2], but occasionally as the Maxwell stress tensor [9]. The momentum continuity equation for the for the cubic isotropic material assumed here. It is straight- electromagnetic field is obtained from Maxwell's equations forward to show with the use of Eqs. (2.3)and(2.4)that products of EoE with Eq. (2.1), B/Ho Eq.(2.3), and B with Eq(2. 4), and (S-m)计+Wcm=-Ej (2.14) then adding the equations. The result after use of stan- dard vector operator identities is here the repeated index i is summed over the Cartesian coordinates x,y, and : and dr: (Tem)yi+ at gem) --PEj-GXB=-Fy, is the usual electromagnetic energy density. Equation(2.14)where sses the continuity of the electromagnetic energy, and the term on the right represents the rate of loss of energy Gcm=E0E×B from the field by transfer to the dielectric. Multiplication of (2.9)by s, gives is the electromagnetic momentum density. Equation (2.21) expresses continuity of electromagnetic momentum. The ms,+mIs:+moTS S = sE 5S, =Ej, ( 2.16) terms on the right represent the rate of loss of momentum from the field by transfer to the dielectric, in the form of where Eqs.(2.6)and(2. 12) have been used, similar to a minus the usual Lorentz force density, denoted Fi calculation in [4]. The rate of loss on the right of the elec- The transfer of momentum from the electromagnetic field tromagnetic energy continuity equation(2.14)is thus bal- implies that the dielectric as a whole is set into motion.The anced by the rate of gain of energy represented by the term internal relative displacement field s is itself invariant under on the right of Eq(2. 16)for the dielectric lattice The a uniform displacement of the crystal, and cannot therefore sum of Eqs. (2. 14)and(2. 16)can be written in the form of carry momentum. The dielectric momentum is carried by the an energy continuity equation for the coupled electromag- motion of the spatial displacement field R=R(r, t) defined by netic field and dielectric lattice he position of the center of mass of the two ions in the unit cell. A treatment of the propagation of momentum through as an (2. 17) cludes both the relative and center-of-mass coordinates s and R; this is provided by the Lagrangian formalism presented in here the total-energy current density Sec. ll S=Scm=E×B/0 The effect of the dissipation term in the internal equation (2.18) of motion(2.9) is to remove energy from the optic modes of is the same as the electromagnetic current density(2.13),but vibration. The sink for this energy is provided by a reservoir, the total energy density is whose nature is determined by the microscopic mechanism of the dissipation. For example, anharmonic forces in the W=IE0E+uoH2+ms2+mas2).(2.19) lattice transfer the optic-mode energy into continuous distri butions of other vibrational modes which are not directly The excitation of the dielectric lattice. that is of the Lorent- coupled to the electromagnetic field. Thus an initial excita- zian oscillator or optic mode, thus makes no explicit contri- tion of the coupled electromagnetic field and optic modes bution to the energy current density. The lattice does, decays to a steady state in which all of the energy is trans- though, have an implicit effect via the scaling of the ratio of ferred to the reservoir. This transfer has implications for both the magnetic and electric fields by the complex refractive the momentum and the kinetic energy associated with the index of the medium [see 0)]. However, the energy motion of the dielectric cryst: density explicitly contains the kinetic and potential energies Suppose that the initial excitation has N quanta of wave of the optic vibrational mode in addition to the electromag- vector k and frequency w per unit volume. The magnitude of netic energy density(2.15). The term on the right of Eq. the momentum density acquired by the dielectric crystal as a coupled field-lattice system by the optic mode sity from the whole, when all of the energy has been transferred to the (2. 17)represents the rate of loss of energy de damping reservoIr. IS C Momentum propagation MR= Nhk=Nho/c 2.23 The flow of electromagnetic momentum is determined by where M is the dielectric mass dens the momentum current density, whose components are given by[12,6,9,10] M=(M1+M2)/, (224)
B. Energy propagation The flow of electromagnetic energy through the dielectric is determined by the energy current density, or Poynting vector, given by Sem5E3B/m0 ~2.13! for the cubic isotropic material assumed here. It is straightforward to show with the use of Eqs. ~2.3! and ~2.4! that ] ]ri ~Sem!i1 ] ]t Wem52E•j, ~2.14! where the repeated index i is summed over the Cartesian coordinates x, y, and z, and Wem5 1 2 «0E21 1 2 m0H2 ~2.15! is the usual electromagnetic energy density. Equation ~2.14! expresses the continuity of the electromagnetic energy, and the term on the right represents the rate of loss of energy from the field by transfer to the dielectric. Multiplication of ~2.9! by s˙i gives ms¨is˙i1mGs˙i 2 1mvT 2 sis˙i5§Eis˙i5E•j, ~2.16! where Eqs. ~2.6! and ~2.12! have been used, similar to a calculation in @4#. The rate of loss on the right of the electromagnetic energy continuity equation ~2.14! is thus balanced by the rate of gain of energy represented by the term on the right of Eq. ~2.16! for the dielectric lattice mode. The sum of Eqs. ~2.14! and ~2.16! can be written in the form of an energy continuity equation for the coupled electromagnetic field and dielectric lattice, ]Si ]ri 1 ]W ]t 52mGs˙ 2, ~2.17! where the total-energy current density S5Sem5E3B/m0 ~2.18! is the same as the electromagnetic current density ~2.13!, but the total energy density is W5 1 2 $«0E21m0H21ms˙ 21mvT 2 s2 %. ~2.19! The excitation of the dielectric lattice, that is of the Lorentzian oscillator or optic mode, thus makes no explicit contribution to the energy current density. The lattice does, though, have an implicit effect via the scaling of the ratio of the magnetic and electric fields by the complex refractive index of the medium @see Eq. ~4.10!#. However, the energy density explicitly contains the kinetic and potential energies of the optic vibrational mode in addition to the electromagnetic energy density ~2.15!. The term on the right of Eq. ~2.17! represents the rate of loss of energy density from the coupled field-lattice system by the optic mode damping. C. Momentum propagation The flow of electromagnetic momentum is determined by the momentum current density, whose components are given by @1,2,6,9,10# ~Tem!ji52«0EjEi1 1 2 «0E2d ji2 BjBi m0 1 B2 2m0 d ji . ~2.20! This quantity is usually identified as the negative of the Maxwell stress tensor @1,2#, but occasionally as the Maxwell stress tensor @9#. The momentum continuity equation for the electromagnetic field is obtained from Maxwell’s equations by forming the vector products of «0E with Eq. ~2.1!, B/m0 with Eq. ~2.2!, E with Eq. ~2.3!, and B with Eq. ~2.4!, and then adding the four equations. The result after use of standard vector operator identities is ] ]ri ~Tem!ji1 ] ]t ~Gem!j52rEj2~j3B
1074 R LOUDON. L. ALLEN. AND D. F. NELSON and the relation between the frequency and wave vector has been taken in its free-space form, for the purpose of an order- of-magnitude estimate. Clearly it is important to include the where the charge and current densities are related to the di center-of-mass momentum of the crystal in any theory of electric polarization by Eqs. (2.5)and(2.6).However, when nomentum propagation through an absorbing dielectric the center-of-mass motion is included, the latter expression The transfer of momentum to the dielectric must be ac- should be augmented by inclusion of the Rontgen current ompanied by a growth in its kinetic energy density, whose [11], to give a total current density value for the momentum density given by Eq (2.23)is aP MR +V×(P×R (36) where r is again the continuum center-of-mass coordinate The rest-mass energy density Mc2 of the crystal is The interaction Lagrangian density(3. 5)can be converted ery much larger than the initial energy density Ni with the use of this expression to rystal kinetic energy is thus completely negligible cor to Nho. This justifies the neglect of center-of-mass motion C1=P(E+R×B) (3.7) in the theory of energy propagation given in Sec. II B, de spite its importance in the theory of momentum propagation. where some perfect space and time derivative terms, which make no contribution to the Lagrange equations of motion, IIL LAGRANGIAN THEORY OF ELECTROMAGNETIC have been discarded [10]. The material Lagrangian density MOMENTUM PROPAGATION for a rigid body is This section is devoted to a rigorous derivation of the CM=MR2+显ms2-显mos2 (38) various momentum densities associated with the propagation of electromagnetic waves through absorbing dielectrics. The where the dielectric parameters are as defined in Sec. IL. The basic dielectric model is the same as that used in Sec. II, but theory also needs to include a term that allows for damping it is necessary to generalize the model to include center-of- of the internal motion at a rate proportional to T. This mass motion in order to describe momentum propagation. It conveniently implemented by a Rayleigh dissipation func- is also necessary to distinguish the contributions of momen- tion of the form tum and pseudomomentum. The continuum mechanics back ground to the calculations is described in detail in Ref [10] R=m2, (3.9) It is assumed throughout that the dielectric material fills all of space the effects of crystal boundaries are excluded from which is incorporated into the Euler-Lagrange equations by the calculations an appropriate additional term [12] The equations of motion for the electromagnetic and ma- A. Lagrangian formulation terial field variables are obtained by the standard Lagrangian The system of dielectric material (M) and electron procedures. Thus the Maxwell-Lorentz equations(2. 1)and netic field (F) coupled by an electric-dipole interaction(n) (2. 4)are rederived straightforwardly, while Eqs.(2.2)and described by a Lagrangian density (2.3)are satisfied automatically from the definitions (3.3) and (3. 4)of the fields in terms of the potentials. It should L=LMt Lte (3. 1) however be noted that the Rontgen term in the current den- sity(3. 6)causes a generalization of relation(2.8)between where the Lagrangian itself is formed by integration over the magnetic field and magnetic induction to [see, for example, Lagrangian density in the usual way Eq(76. 11)of Ref [21 The Lagrangian density of the electromagnetic field is (3.10) The new term is a function of both the internal relative- where the electric and magnetic fields are determined by the displacement coordinate and the center-of-mass coordinate calar potential and the vector potential A in the usual of the dielectric material For the dielectric spatial displacement variables, the equa- tion of motion for the relative position of the two ions in the ⅴd-at (3.3) unit cell is obtained with the use of Eqs. (3. 7)-(3.9)as 5+mIs;+mo3s;=s(E1+(R×B)),(3,1) which is identical to Eq.(2.9)except for the addition of the B=V×A (3.4) term proportional to the center-of-mass velocity R. The equation of motion for the continuum center-of-mass coordi- The interaction Lagrangian density is is obtained similarly
and the relation between the frequency and wave vector has been taken in its free-space form, for the purpose of an orderof-magnitude estimate. Clearly it is important to include the center-of-mass momentum of the crystal in any theory of momentum propagation through an absorbing dielectric. The transfer of momentum to the dielectric must be accompanied by a growth in its kinetic energy density, whose value for the momentum density given by Eq. ~2.23! is MR˙ 2 2 5N\v N\v 2M c2 . ~2.25! The rest-mass energy density M c2 of the crystal is always very much larger than the initial energy density N\v. The crystal kinetic energy is thus completely negligible compared to N\v. This justifies the neglect of center-of-mass motion in the theory of energy propagation given in Sec. II B, despite its importance in the theory of momentum propagation. III. LAGRANGIAN THEORY OF ELECTROMAGNETIC MOMENTUM PROPAGATION This section is devoted to a rigorous derivation of the various momentum densities associated with the propagation of electromagnetic waves through absorbing dielectrics. The basic dielectric model is the same as that used in Sec. II, but it is necessary to generalize the model to include center-ofmass motion in order to describe momentum propagation. It is also necessary to distinguish the contributions of momentum and pseudomomentum. The continuum mechanics background to the calculations is described in detail in Ref. @10#. It is assumed throughout that the dielectric material fills all of space; the effects of crystal boundaries are excluded from the calculations. A. Lagrangian formulation The system of dielectric material (M) and electromagnetic field (F) coupled by an electric-dipole interaction (I) is described by a Lagrangian density L5LM1LI1LF , ~3.1! where the Lagrangian itself is formed by integration over the Lagrangian density in the usual way. The Lagrangian density of the electromagnetic field is LF5«0 2 E22 1 2m0 B2, ~3.2! where the electric and magnetic fields are determined by the scalar potential f and the vector potential A in the usual way, E52“f2 ]A ]t ~3.3! and B5“3A. ~3.4! The interaction Lagrangian density is LI5j•A2rf, ~3.5! where the charge and current densities are related to the dielectric polarization by Eqs. ~2.5! and ~2.6!. However, when the center-of-mass motion is included, the latter expression should be augmented by inclusion of the Ro¨ntgen current @11#, to give a total current density j5 ]P ]t 1“3~P3R˙ !, ~3.6! where R is again the continuum center-of-mass coordinate. The interaction Lagrangian density ~3.5! can be converted with the use of this expression to Ll5P•~E1R˙ 3B!, ~3.7! where some perfect space and time derivative terms, which make no contribution to the Lagrange equations of motion, have been discarded @10#. The material Lagrangian density for a rigid body is LM5 1 2 MR˙ 21 1 2 ms˙22 1 2 mvT 2 s 2, ~3.8! where the dielectric parameters are as defined in Sec. II. The theory also needs to include a term that allows for damping of the internal motion at a rate proportional to G. This is conveniently implemented by a Rayleigh dissipation function of the form R5 1 2 mGs˙2, ~3.9! which is incorporated into the Euler-Lagrange equations by an appropriate additional term @12#. The equations of motion for the electromagnetic and material field variables are obtained by the standard Lagrangian procedures. Thus the Maxwell-Lorentz equations ~2.1! and ~2.4! are rederived straightforwardly, while Eqs. ~2.2! and ~2.3! are satisfied automatically from the definitions ~3.3! and ~3.4! of the fields in terms of the potentials. It should however be noted that the Ro¨ntgen term in the current density ~3.6! causes a generalization of relation ~2.8! between magnetic field and magnetic induction to @see, for example, Eq. ~76.11! of Ref. @2## H5 B m0 2P3R˙ . ~3.10! The new term is a function of both the internal relativedisplacement coordinate and the center-of-mass coordinate of the dielectric material. For the dielectric spatial displacement variables, the equation of motion for the relative position of the two ions in the unit cell is obtained with the use of Eqs. ~3.7!–~3.9! as ms¨i1mGs˙i1mvT 2 si5§„Ei1~R˙ 3B!i…, ~3.11! which is identical to Eq. ~2.9! except for the addition of the term proportional to the center-of-mass velocity R˙ . The equation of motion for the continuum center-of-mass coordinate is obtained similarly as 1074 R. LOUDON, L. ALLEN, AND D. F. NELSON 55
PROPAGATION OF ELECTROMAGNETIC ENERGY AN 1075 with the same magnitude as the center-of-mass kinetic en- MR,-sS, dr:(E, +(RXB))-s dt (sx b))=0 rgy has been omitted from the momentum current density (3.12)(3.15)in accordance with the discussion that follows Eq A more convenient form of this equation is found after con- It should be noted that the equation of motion(3. 11)for siderable manipulation [10], using Eqs. 2.2),(2.5)and(3.6), the internal coordinate plays no role in the above derivation, to be on account of the inability of the relative coordinate to carry momentum [6]. As the dissipation described by Eq (3.9)acts s(E+(R×B)s}=pE+(j×B)=F only on the internal coordinate, the absence of this coordi nate from the derivation of Eqs. (3.14)-(3. 16)accounts for (3.13) the lack of damping terms in these results. The fulfilment of momentum conservation in the presence of energy dissipa- The significances of the terms on the left are clarified by tion, as described by the term on the right-hand side of ec their contributions to the momentum densities defined in Sec. (2.17), is of course a common occurrence in mechanics III B below As discussed in connection with Eq. (2.25), the kinetic contribution mr from the center-of-mass motion and a con- energy of the dielectric is negligible compared to the energy tribution EoEXB from the electromagnetic field. The latter densities associated with the electromagnetic field and the differs from the abraham form EoloEXH on account of the optic mode vibration of the lattice. The contribution of the generalized relation(3. 10) between the magnetic field and center-of-mass motion to the Lorentz force on the right-hand the induction. It is seen that, in contrast to the magnetic side of (3. 11)is also small, as RB is of order RE/c, which is induction B, which is purely a property of the electromag certainly much smaller than E. With terms in R removed, the netic field, the magnetic field H contains a contribution that Lagrangian theory reproduces the results of Sec. IB as Eq. depends on the material internal coordinate contained in P (3. 11)reduces to Eq.(2.9), and the Maxwell equations are The occurrence of b rather than h ensures that the electro. unchanged. The energy continuity equation (2.17)and the magnetic momentum density is properly independent of any energy densities given in Eqs.(2.18)and (2.19)thus con- material variables The total momentum is obtained by integration of Gn (r, t) subtle problems and distinctions associated with the propa- over all space, and this quantity is conserved for a closed sorbing dielectric Thus integration of Eq. (3. 14)over the effectively-infinite dielectric material gives B. Momentum conservation The law of momentum conservation is a consequence of (3.17) the invariance of the laws of physics to arbitrary infinitesimal dt dr gm(r, 0=0, displacements of the spatial coordinates. The momentum is thus defined with respect to the Maxwell equations(2.1)- provided that Tm(r, t)vanishes at r=oo. The total momentum (2.4)and the center-of-mass equation (3. 13). The continuty electromagnetic field and the dielectric center-of-mass mo- equation(2.21)for the electromagnetic momentum, which is based entirely on Maxwell's equations, therefore remains tion changes with time as, for example, in the propagation of valid. The Lorentz force densities F, on the right-hand sides a pulse of excitation through the crystal Although the material motion makes an important contri of Eqs.(2.21)and(3. 13)are equal and opposite, demonstrat- bution to the momentum density, its contribution to the mo- magnetic and material parts of the coupled system. Addition mentum current density(3. 15)is generally less important of these equations, using definition(2. 12)of the polarization, E, and, with the corresponding term removed, Eq. (3.15) reduces to (3.14) (Tm=-E,P+(T which is the conservation law for the momentum of the field- C. Pseudomomentum aterial system. Here Pseudomomentum has been a much neglected quantity in (Tm)ji=-(E+(RX B)))P, +(Tem)ji (3. 15) continuum mechanics, and a regularly misin erpreted quar tity in quantum mechanics. Quantum-mechanical treatments is the momentum current density, and of excitations in solids have often called hk the pseudomo- mentum of an excitation quantum. This was shown to be Gm=Mr+ gem (3. 16) wrong [6] on the basis of an unambiguous definition of the is the momentum density of the coupled field and material. served by virtue of the homogeneity of the material dog a pseudomomentum as the momentumlike quantity that is ce Expressions for the electromagnetic contributions to the Noether's theorem [10] can be used with a Lagrangian mentum current density and the momentum density are given formulation to obtain a rigorous derivation of the pseudomo- in Eqs. (2.20)and(2.22), respectively. A tensor contribution mentum conservation law for a homogeneous body. This is a
MR¨ j2§si ] ]rj „Ei1~R˙ 3B!i…2§ d dt ~s3B!j50. ~3.12! A more convenient form of this equation is found after considerable manipulation @10#, using Eqs.~ 2.2!, ~2.5! and ~3.6!, to be MR¨ j2 ] ]ri $§„Ej1~R˙ 3B!j…si%5rEj1~j3B!j5Fj . ~3.13! The significances of the terms on the left are clarified by their contributions to the momentum densities defined in Sec. III B below. As discussed in connection with Eq. ~2.25!, the kinetic energy of the dielectric is negligible compared to the energy densities associated with the electromagnetic field and the optic mode vibration of the lattice. The contribution of the center-of-mass motion to the Lorentz force on the right-hand side of ~3.11! is also small, as R ˙ B is of order R ˙ E/c, which is certainly much smaller than E. With terms in R˙ removed, the Lagrangian theory reproduces the results of Sec. II B as Eq. ~3.11! reduces to Eq. ~2.9!, and the Maxwell equations are unchanged. The energy continuity equation ~2.17! and the energy densities given in Eqs. ~2.18! and ~2.19! thus continue to hold. The derivations that follow consider the more subtle problems and distinctions associated with the propagation of momentum and pseudomomentum through the absorbing dielectric. B. Momentum conservation The law of momentum conservation is a consequence of the invariance of the laws of physics to arbitrary infinitesimal displacements of the spatial coordinates. The momentum is thus defined with respect to the Maxwell equations ~2.1!– ~2.4! and the center-of-mass equation ~3.13!. The continuity equation ~2.21! for the electromagnetic momentum, which is based entirely on Maxwell’s equations, therefore remains valid. The Lorentz force densities Fj on the right-hand sides of Eqs. ~2.21! and ~3.13! are equal and opposite, demonstrating the action and reaction of the forces between the electromagnetic and material parts of the coupled system. Addition of these equations, using definition ~2.12! of the polarization, gives ] ]ri ~Tm!ji1 ] ]t ~Gm!j50, ~3.14! which is the conservation law for the momentum of the fieldmaterial system. Here ~Tm!ji52„Ej1~R˙ 3B!j…Pi1~Tem!ji ~3.15! is the momentum current density, and Gm5MR˙ 1Gem ~3.16! is the momentum density of the coupled field and material. Expressions for the electromagnetic contributions to the momentum current density and the momentum density are given in Eqs. ~2.20! and ~2.22!, respectively. A tensor contribution with the same magnitude as the center-of-mass kinetic energy has been omitted from the momentum current density ~3.15! in accordance with the discussion that follows Eq. ~2.25!. It should be noted that the equation of motion ~3.11! for the internal coordinate plays no role in the above derivation, on account of the inability of the relative coordinate to carry momentum @6#. As the dissipation described by Eq. ~3.9! acts only on the internal coordinate, the absence of this coordinate from the derivation of Eqs. ~3.14!–~3.16! accounts for the lack of damping terms in these results. The fulfilment of momentum conservation in the presence of energy dissipation, as described by the term on the right-hand side of Eq. ~2.17!, is of course a common occurrence in mechanics. The momentum density ~3.16! is clearly separated into a contribution MR˙ from the center-of-mass motion and a contribution «0E3B from the electromagnetic field. The latter differs from the Abraham form «0m0E3H on account of the generalized relation ~3.10! between the magnetic field and the induction. It is seen that, in contrast to the magnetic induction B, which is purely a property of the electromagnetic field, the magnetic field H contains a contribution that depends on the material internal coordinate contained in P. The occurrence of B rather than H ensures that the electromagnetic momentum density is properly independent of any material variables. The total momentum is obtained by integration of Gm~r,t! over all space, and this quantity is conserved for a closed system with no flow of momentum through its boundaries. Thus integration of Eq. ~3.14! over the effectively-infinite dielectric material gives ] ]t E dr Gm~r,t!50, ~3.17! provided that Tm~r,t! vanishes at r5`. The total momentum is therefore conserved, and only its division between the electromagnetic field and the dielectric center-of-mass motion changes with time as, for example, in the propagation of a pulse of excitation through the crystal. Although the material motion makes an important contribution to the momentum density, its contribution to the momentum current density ~3.15! is generally less important. Thus R ˙ B is again of order R ˙ E/c, which is much smaller than E, and, with the corresponding term removed, Eq. ~3.15! reduces to ~Tm!ji52EjPi1~Tem!ji . ~3.18! C. Pseudomomentum Pseudomomentum has been a much neglected quantity in continuum mechanics, and a regularly misinterpreted quantity in quantum mechanics. Quantum-mechanical treatments of excitations in solids have often called \k the pseudomomentum of an excitation quantum. This was shown to be wrong @6# on the basis of an unambiguous definition of the pseudomomentum as the momentumlike quantity that is conserved by virtue of the homogeneity of the material body. Noether’s theorem @10# can be used with a Lagrangian formulation to obtain a rigorous derivation of the pseudomomentum conservation law for a homogeneous body. This is a 55 PROPAGATION OF ELECTROMAGNETIC ENERGY AND... 1075