Ab initio study of the radiation pressure on dielectric and magnetic media

The Maxwell stress tensor and the distributed Lorentz force are applied to calculate forces on lossless media and are shown to be in excellent agreement. From the Maxwell stress tensor, we derive analytical formulae for the forces on both a half-space and a slab under plane wave incidence. It is shown that a normally incident plane wave pushes the slab in the wave propagation direction, while it pulls the half-space toward the incoming wave. Zero tangential force is derived at a boundary between two lossless media, regardless of incident angle. The distributed Lorentz force is applied to the slab in a direct way, while the half-space is dealt with by introducing a finite conductivity. In this regard, we show that the ohmic losses have to be properly accounted for, otherwise differing results are obtained. This contribution, together with a generalization of the formulation to magnetic materials, establishes the method on solid theoretical grounds. Agreement between the two methods is also demonstrated for the case of a 2-D circular dielectric particle. © 2005 Optical Society of America OCIS codes: (260.2110) Electromagnetic Theory; (140.7010) Trapping, (290.5850) Particle Scattering References and links 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. M. Mansuripur, A. R. Zakharian, and J. V. Moloney, ”Radiation pressure on a dielectric wedge,” Opt. Express 13, 2064-2074 (2005), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-6-2064. 3. M. Mansuripur, ”Radiation pressure and the linear momentum of light in dispersive dielectric media,” Opt. Express 13, 2245-2250 (2005), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-6-2245. 4. A. R. Zakharian, M. Mansuripur, and J. V. Moloney, ”Radiation pressure and the distribution of electromagnetic force in a dielectric media,” Opt. Express 13, 2321-2336 (2005), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-7-2321. 5. R. Loudon, “Theory of radiation pressure on dielectric surfaces,” J. Mod. Opt. 49, 821-838 (2002). 6. R. Loudon, S. M. Barnett, and C. Baxter, “Radiation pressure and momentum transfer in dielectrics: the photon drag effect,” Phys. Rev. A 71, 063802 (2005). 7. J. A. Stratton, Electromagnetic Theory (McGraw-Hill, 1941), ISBN 0-07-062150-0. 8. J. A. Kong, Electromagnetic Wave Theory (EMW, 2005), ISBN 0-9668143-9-8. 9. L. Tsang, J. A. Kong, and K. Ding, Scattering of Electromagnetic Waves: Theories and Applications (Wiley, 2000), ISBN 0-471-38799-7.


Introduction
Recently, the radiation pressure exerted by an electromagnetic field impinging on a medium was derived by the direct application of the Lorentz law [1]. This method applies the Lorentz force to bound currents distributed throughout the medium and bound charges at the surface of the medium. The approach allows for the computation of force at any point inside a dielectric [1,2,3,4] and has been shown applicable to numerical methods, such as the finite-differencetime-domain (FDTD) [4]. A similar approach was taken to study the radiation pressure on a dielectric surface [5] and a semiconductor exhibiting the photon drag effect [6]. Comparison with the established stress tensor approach has not been done previously. The purpose of this paper is to compare this approach with the application of the Maxwell stress tensor [7,8] to lossless media. In this paper, two important generalizations of the method proposed in [1] are introduced that allow for the calculation of forces on lossless media. First, the force on a magnetic polarization vector and a bound magnetic charge density are introduced to the distributed Lorentz force to model magnetic materials. The second generalization allows for the discrimination of the force on free carriers in a slightly conducting medium, which may not contribute to the total force on the bulk medium. In addition to generalizing the Lorentz method, we also derive closed-form expressions from the Maxwell stress tensor for the force on a lossless slab and a half-space. At normal incidence, the force on a slab is in the wave propagation direction, while the force on a half-space pulls it toward the incoming wave. At any incident angle, the tangential force at a boundary is shown to be zero. The force on a 2D dielectric cylinder is calculated by both methods and shown to be in agreement. Furthermore, we contrast the two methods in their relative advantages and disadvantages.

Force calculation methods
The time-average forces of electromagnetic waves incident on dielectric and magnetic bodies are calculated from the total complex field vectors due to both the incident waves and the scattered waves from the bodies. The Maxwell stress tensor approach to compute these forces is considered first. Second, the direct application of the Lorentz force is generalized from [1,4].

The Maxwell stress tensor
The momentum conservation theorem is derived from the Maxwell equations and the Lorentz force and is given by [8] wherer and t refer to position and time, respectively,Ḡ(r,t) =D(r,t)×B(r,t) is the momentum density vector,T (r,t) is the time-domain Maxwell stress tensor, and f (r,t) is the force density in [N/m 3 ]. The momentum density vector is fundamentally defined in terms of the electric flux densityD(r,t) and the magnetic flux densityB(r,t). By integrating over a volume V enclosed by a surface S and using the divergence theorem, the total force can be written as We are generally interested in the time average force, which is since the time average of the first term on the right-hand side of equation (2) is zero. The complex Maxwell stress tensor for lossless media is whereĪ is the 3 × 3 identity matrix and ( * ) denotes the complex conjugate. Eqs. (3) and (4) are applied later to calculate the force on media.

Lorentz force on a host medium
The Lorentz force is applied to both bound currents due to the polarization of a dielectric and bound charges at the boundaries resulting from discontinuity ofn · ε 0Ē [1]. The time-average Lorentz force in [N/m 3 ] is found to bē where ( * ) denotes a complex conjugate. In Eq. (5), bound magnetic currentM and bound magnetic charge ρ m have been added to the formulation of [1] to account for permeable media. The magnetic fluxB is modeled by a magnetic polarization vectorP m such thatB = µ r µ 0H = µ 0H −P m and it follows thatP By application of Gauss's law, the bound magnetic charges are given by ρ m = ∇ ·P m = ∇ · µ 0H . Invoking charge continuity leads to an expression for the bound magnetic current density, which can be used directly in equation (5) M = −iωP m = iω µ 0 (µ r − 1)H.
By duality a similar expression for the bound electric current density is obtained as in [1] where ε r can be complex to account for material losses [4]. The charge distributions in equation (5) are found at a medium boundary by considering the discontinuity in the normal fields.
For example, at a boundary between two media referenced by the subscripts 0 and 1, wheren is a unit vector normal to the surface pointing from region 0 to region 1. When applying Eq. (9), the average of the normal field vectors across the boundary should be used [1]. To get the total force on a material body from Eq. (5), the contribution from distributed current densitiesJ andM are integrated over the volume of the medium and the effect of bound charge densities ρ e and ρ m are included at the medium boundaries.

Lossless slab
First, the problem of TE incidence on a lossless slab is considered, as shown in Fig. 1. The forces from a TM polarized wave are directly obtained from the duality principle. The forces are calculated by using the two methods which are referred to as stress tensor and Lorentz, from section 2.1 and section 2.2, respectively. The force on an isotropic slab characterized by µ 1 = µ r µ 0 and ε 1 = ε r ε 0 is evaluated in free space (µ 2 = µ 0 = 4π · 10 −7 H/m, ε 2 = ε 0 = 8.85 · 10 −12 F/m) due to an incident TE plane waveĒ where E i is the incident field magnitude (see Fig. 1). The total fields in the three regions are found by application of the boundary conditions [8]. The Maxwell stress tensor is applied through Eq. (3) by selecting the integration path shown in Fig. 1. This integration path is chosen such that the fields are evaluated at the boundary as δ z → 0. The force per unit area on the slab of thickness d is, therefore, given bȳ whereT (z = z 0 ) is the Maxwell stress tensor of equation (4) evaluated at the point z = z 0 . The contributions from the fields inside the slab are restricted to the termŝ By substitution of the fields in the slab, it can be shown that and the radiation pressure on the slab reduces tō Therefore, the force on a lossless slab can be computed solely from the knowledge of the fields outside the slab. This is to be expected and can be generalized to media of arbitrary geometry since the divergence of the stress tensor applied to continuous, lossless media is where θ 0 is the incident angle and R slab and T slab are the slab reflection and transmission coefficients, respectively [8]. This analytic expression shows that thex-component of the force is zero (F x = 0), regardless of incident angle. At normal incidence, equation (15) can be written simply as the summation of force components due to momentum conservation, whereF i ,F r , andF t are the forces due to the incident, reflected, and transmitted wave momentums, respectively. Equation (5) is used to calculate the Lorentz force on the slab, and the results are compared with Eq. (15). Figure 2 shows excellent agreement between the two methods applied to compute the force on a lossless dielectric slab (µ r = 1, ε r = 4). The maxima and minima in the force are due to the periodic dependence of the reflection coefficient R slab and the transmission coefficient T slab on the slab thickness. Minimum force is observed for slab thicknesses equal to multiples of half wavelength n λ 1 2 and maximum force is observed for slab thicknesses equal to odd multiples of quarter wavelength (2n + 1) λ 1 4 , where λ 1 is the wavelength of the electromagnetic wave inside the slab and n ∈ [0, 1, 2, . . .]. For comparison, note that is calculated by the Lorentz method and f = 3.184 [pN/m 2 ] is calculated by the application of the Maxwell stress tensor, which is in reasonable agreement with f = 3.188 [pN/m 2 ] previously reported [4].  Force density from a normal incident wave onto a lossless slab as a function of slab thickness d. The free space wavelength is λ 0 = 640nm and ε r = 4, µ r = 1, E i = 1. Shown are the forces calculated from the distributed Lorentz force (circles) and the Maxwell stress tensor (line). The background medium (region 0 and region 2) is free space. Figure 3 shows the force on a permeable slab (µ r = 4, ε r = 1) as a function of incident angle θ 0 , measured from the surface normal. This example requires the calculation of both magnetic currents from Eq. (7) and magnetic surface charge from Eq. (9), which are used to compute the Lorentz force. As seen in Fig. 3, there is excellent agreement between the two force calculation methods at all incident angles. The Brewster angle for total transmission of a TE incident wave is given by [8] and is calculated for this example to be θ b = 63.4 • . This is evident in Fig. 3 by the zero force at this particular angle. The force is shown as a function of dielectric constant ε r for a lossless dielectric slab in Fig. 4. Again, there is excellent agreement between the two force calculation methods. It is seen that the force goes to zero as expected when the slab is impedance matched to free space ε r = 1. For all positive permittivities, including the region 0 ≤ ε r ≤ 1, the force is in the positivê z-direction, indicating that the force is pushing the slab.

Semi-infinite half-space
Calculation of the Lorentz force on a half-space is more involved than what might be expected. For a lossless half-space, the fields propagate inside the medium without attenuation. This poses a problem in finding the radiation pressure by integrating over the Lorentz force from z = 0 to z → ∞. Previously, this issue has been sidestepped by introducing a small amount of loss in the medium, applying the distributed Lorentz force, and allowing the losses to approach zero after integration [1,5]. The problem with this approach is that not all of the force on free charges can be attributed to force on the bulk medium. Some of the wave energy may be lost, such as ohmic losses in a conducting medium, which must be considered when calculating the total force on a material body. This problem will be addressed subsequently.
First, we derive the exact solution to the half-space problem of an incident TE plane wave by applying the Maxwell stress tensor to the boundary as described in section 2.1. The path of integration for the Maxwell stress tensor is the same as that shown in Fig. 1 except that there is only one boundary at z = 0. Again we allow the integration to just enclose the boundary such that δ z → 0. The force per unit area on the half space medium is given bȳ The contributions on the two sides of the interface arê and the force tangential to the boundary is seen to be Upon substitution of the fields on both sides of the boundary, the tangential force is simplified to where R hs and T hs are the half-space reflection and transmission coefficients, respectively and p 01 is given by [8] p 01 = µ 0 k 1z µ 1 k 0z .
The simplification to Eq. (21) is a direct result of phase matching, which requires k x to be continuous across the boundary. Applying the formulas T hs = 1 + R hs and p 01 T hs = 1 − R hs obtained from the boundary conditions, it immediately follows that so that the tangential force is zero (F x = 0). This result is a direct consequence of the boundary conditions and phase matching at a single interface, and it can be generalized to multiple interfaces, including the slab of section 3, the details of which were omitted for brevity. This differs from previous reports of derived nonzero tangential force at a planar boundary [1,4]. The radiation pressure on a lossless half-space medium is found from Eq. (18) to bē where ε r is the relative permittivity of the half-space and θ 1 is the transmitted angle found from Snell's law. Note that this equation also accounts for variations in relative permeability µ r , which appears in the reflection and transmission coefficients. Again we can see the contribution from the incident, reflected, and transmitted wave momentums on the force at the medium interface.
Equation (24) is applied to evaluate the radiation pressure on the half-space medium due to a normal incident plane wave. By taking the limit ε r → ∞, it is simple to show that the radiation pressure approaches lim This limit is seen in Fig. 5 where for very large ε r . Thus, the force of a normally incident plane wave on a lossless dielectric half-space is pulling toward the incoming wave, while in the slab case the plane wave pushes the medium in the wave propagating direction. Likewise, the analytical limit of ε r → 0 gives which is also seen in Fig. 5. Equation (24) differs from the TE incident results F z = ε 0 2 E 2 i 1 + R 2 hs cos 2 θ 0 and F x = ε 0 2 E 2 i 1 − R 2 hs sin θ 0 cos θ 0 previously reported from the distribution of Lorentz force in a lossless media [1]. The previously reported results produce different values from Eq. (24) for three cases in particular. First, for an impedance matched interface such as a fictive interface between two free space media, R hs = 0, which would yield a nonzero value of F z =ẑ ε 0 2 E 2 i . Since such a fictive boundary could be placed anywhere, it would indicate that forces in free space are present everywhere. Second, a nonzero tangential force F x = 0 at oblique incidence [1,4] disagrees with Eq. (24). Third, the predicted normal force F z is always positive for a lossless dielectric, which contradicts the known theory that the force of a normally incident plane wave on a lossless half-space is pulling toward the incident wave [8].
The differing results obtained from the distribution of Lorentz force and the Maxwell stress tensor do not imply that either result is incorrect. Instead, we interpret the force on a semiinfinite half-space obtained from the method of [1] to be the total Lorentz force on all bound and free charges and currents with the assumption that the fields attenuate to zero as z → ∞ due to some finite loss, which holds once the force density is integrated, even if the losses approach zero. To illustrate this point, the force on a lossless half-space due to a normally incident wave can be computed from the distributed Lorentz force by considering the ohmic loss due to the  Fig. 5. Force density from a normal incident wave on a lossless half-space medium as calculated from the Maxwell stress tensor. Theẑ-directed force is shown as a function of the relative permittivity ε r . The free space wavelength is λ 0 = 640nm and µ r = 1, E i = 1. The incident media (region 0) is free space. The dotted lines denote ±ε 0 E 2 i = 8.85 pN/m 2 conductivity σ in a slightly lossy medium. Some of the energy transferred from the wave to the conduction currentJ c is dissipated as ohmic loss. The time average power dissipated is [8] which must be subtracted from the Lorentz force calculated on the conduction currentJ c . Therefore, the total radiation pressure on a dielectric half-space due to a normally incident plane wave isF whereJ is the current from Eq. (8),J c = σĒ is the conduction current, and v e is the energy velocity. After integration, theẑ-directed force reduces to where k I is the imaginary part of the wavenumber(k I = Im{k 1 }). As σ → 0, the energy velocity and attenuation constant k I can be approximated by [8] v e ≈ 1 √ Thus, the total force on the dielectric half-space reduces to where the first term is a result of the calculated Lorentz force [1] and the second term is derived from the ohmic losses due to the finite conductivity. This expression is in agreement with Eq. (24) at normal incidence θ 0 = 0, which is derived from the Maxwell stress tensor without approximation.

Cylindrical particle
In this section, the force on a dielectric particle due to the superposition of three incident plane waves is calculated by the two methods of section 2. The fields are calculated by the Mie theory as in [9] and verified through comparison with results from the commercial package CST Microwave Studio ®. The force distribution [N/m 3 ] inside the particle is then calculated using the Lorentz method. The distributed Lorentz force is numerically integrated throughout the interior surface of the particle to obtain the total force [N/m]. This value is then compared with the force obtained from the Maxwell stress tensor by integrating the tensor around a circle completely enclosing the 2D particle. Three plane waves of free space wavelength λ 0 = 532 nm are incident in the (xy) plane with angles, {π/2, 7π/6, 11π/6} [rad]. The particle is a polystyrene cylinder (ε p = 2.56ε 0 ) of radius a = 0.3λ 0 with center position (x 0 , y 0 ) = (0, 100) [nm] , and the background medium is water (ε b = 1.69ε 0 ). The magnitude of the total incident electric field is shown in Fig. 6 with the position of the dielectric particle overlayed. It is seen that the particle is of the same size as a typical gradient trap in the Rayleigh regime (seen from the periodicity in the interference pattern), but poses no problem for the method used here. Figure 7 shows the total field magnitudes resulting from the incident fields and the scattered fields. The symmetry of the fields with respect to x produces zero net force in the x-direction, as we shall see.
The distributed Lorentz force is found from Eq. (5) to bē where the electric polarization vector isP e =ẑ(ε p − ε b )E z . The distributed Lorentz force components ( f x and f y ) inside the particle are shown in Fig. 8. It is seen that f x is symmetric about x = 0, as expected. However, the distribution of f y is more complex, and it is not immediately apparent from the distributed force if the total force F y is positive or negative. The total force is obtained from the Lorentz force of Fig. 8 by simple numerical integration and found to beF =ŷ2.11 · 10 −18 [N/m]. Subsequently, the Maxwell stress tensor of Eq. 4 was integrated along a circular path of radius 2a centered at (x 0 , y 0 ), resulting in a force of F =ŷ2.12 · 10 −18 [N/m]. Note that the exact path of integration for the Maxwell stress tensor is not important as long as it completely encloses the cylinder. The excellent agreement between the two force calculation methods for this example further establishes the method of [1] on a sound theoretical basis. However, there are some differences between the two methods. As shown in Fig. 8, the force distribution inside a medium can be easily calculated and displayed by the application of the Lorentz force. However, the stress tensor approach only yields the total force on the medium, thus it cannot be used to produce the plots in Fig. 8. The primary advantage of the stress tensor is the reduction of computation. For the 2D cylindrical particle, the surface integral used to calculate the total force from the distributed Lorentz force was reduced to a line integral for the application of the Maxwell stress tensor.

Conclusions
The distributed Lorentz force method of [1] was generalized to include contribution from magnetic material. This method was shown to provide results equivalent to those obtained from the application of the Maxwell stress tensor, which was used to derive simple closed form expressions for the force on a lossless slab in equation (15) and the force on a lossless half-space in equation (24). It was shown that in both cases, the force tangential to the interface is zero. A normally incident plane wave pushes a dielectric slab in the wave propagation direction and pulls a dielectric half-space toward the incident wave [8].
It was shown that the previous application of the distributed Lorentz force to the half-space problem lacks to account for the losses not manifested as force on the bulk medium. We argued that the difference between the force on a half-space medium that we calculated using the Maxwell stress tensor and the force derived from the Lorentz force previously reported [1,4,5]   is simply due to the loss of energy by mechanisms such as heat or absorption. As an example, Eq. (24) was derived for the normal incidence case from the Lorentz force with consideration given to ohmic energy loss in the conduction current. Although the two force calculation methods can be applied to calculate the total force on media, there are inherent advantages and disadvantages for each method. Through numerous examples, [4] demonstrated that the distributed Lorentz method is easily applied to finitedifference time domain (FDTD) computational and can be used to find the forces distributed throughout a host medium. In particular, the tension and compression forces due to a light beam inside a dielectric medium were calculated. The Maxwell stress tensor is not suitable for this particular task because it is generally applied to media boundaries. However, the advantage of the Maxwell stress tensor is that it reduces a volume integral in three dimensions to a surface integral, thus reducing the computation required to calculate the total force on a material object. We have demonstrated the validity of these results by comparing the Lorentz force on a 2D particle with the force obtained from the Maxwell stress tensor.