Representation Independent Boundary Conditions for a Piecewise-Homogeneous Linear Magneto-dielectric Medium

At a boundary between two transparent, linear, isotropic, homogeneous materials, derivations of the electromagnetic boundary conditions and the Fresnel relations typically proceed from the Minkowski {E,B,D,H} representation of the macroscopic Maxwell equations. However, equations of motion for macroscopic fields in a transparent linear medium can be written using Ampere {E,B}, Chu {E,H}, Lorentz, Minkowski, Peierls, Einstein-Laub, and other formulations of continuum electrodynamics. We present a representation-independent derivation of electromagnetic boundary conditions and Fresnel relations for the propagation of monochromatic radiation through a piecewise-homogeneous, transparent, linear, magneto-dielectric medium. The electromagnetic boundary conditions and the Fresnel relations are derived from energy conservation coupled with the application of Stokes's theorem to the wave equation. Our representation-independent formalism guarantees the general applicability of the Fresnel relations. Specifically, the new derivation is necessary so that a valid derivation of the Fresnel equations exists for alternative, non-Minkowski formulations of the macroscopic Maxwell field equations.

In continuum dynamics, the usual procedure to derive boundary conditions is to apply conservation of energy and conservation of momentum across the surface of a control volume that contains the boundary. 1 However, such a procedure fails in continuum electrodynamics because the energy and momentum of light are not linearly independent as both depend quadratically on the macroscopic fields. For this case, rather than use conservation principles, the electromagnetic boundary conditions and the Fresnel relations are typically obtained by the application of Stokes's theorem and the divergence theorem to the Maxwell-Minkowski equations [2][3][4][5] ∇ × E + 1 c ∂B ∂t = 0 (1) The macroscopic Minkowski fields, E(r, t), D(r, t), B(r, t), and H(r, t), are functions of position r and time t.
Here, J f (r, t) is the free current, ρ f (r, t) is the free charge density, and ρ b (r) is the bound charge density. The macroscopic electric and magnetic fields are related by the constitutive relations where ε(r, t) is the electric permittivity and µ(r, t) is the magnetic permeability.
the material response is separated from the Chu electric field E c and the Chu magnetic field H c . The Chu polarization P c and Chu magnetization M c are treated as sources for the fields.
Although most of the different expressions of macroscopic electromagnetic theory appear to be justified, there have also been some questions about validity. Before the work of Penfield and Haus, 6 it was believed that optical forces predicted by the Chu formulation, Eqs. (7)-Eqs. (10), differed from other accepted theories. 7,8 Also, the validity of the Einstein-Laub formulation is still debated. 12,13 The problem that we address is that, in some of these representations, 10  The new derivation is necessary so that a valid derivation of the Fresnel relations exists for alternative, non-Minkowski formulations of the macroscopic Maxwell field equations. [6][7][8][9][10][11] As long as energy is conserved and the wave equation is valid in a specific formulation of continuum electrodynamics for a transparent linear medium, the Fresnel relations will also be valid.
We consider an arbitrarily long, nominally monochromatic pulse of light in the planewave limit to be propagating through a transparent, linear, homogeneous, magneto-dielectric medium with permittivity ε 1 and permeability µ 1 to be incident on a plane interface with a second homogeneous transparent linear medium with permittivity ε 2 and permeability µ 2 . The frequency of the field is assumed to be sufficiently far from material resonances that absorption can be neglected. The complete system that consists of a simple transparent linear magneto-dielectric medium plus finite radiation field isolated in free space is thermodynamically closed.
Acceleration of the material due to optically induced forces is viewed as negligible in the lowest order theory. The condition of a stationary, piecewise-homogeneous, linear medium with negligible absorption illuminated by an arbitrarily long, nominally monochromatic field (square/rectangular shape or "top-hat" finite pulse) in the plane-wave limit is a physical abstraction that is used in textbook derivations of the Fresnel relations from the Maxwell-Minkowski equations. [2][3][4][5] We use the same conditions in our derivation. Although "realworld" materials are more complicated than the theoretical models, appeal to complexity does not invalidate the theory for the idealized model.
For a transparent, linear, magneto-dielectric material, the energy formula depends on the representation that is used for the macroscopic fields. In the usual Maxwell-Minkowski representation, the total electromagnetic energy is conserved in a volume of space that contains all fields present by extending the region of integration to all-space σ. For a finite homogeneous medium, we can write the energy formula as where V is the integration volume and The electric refractive index n e and the magnetic refractive index n m can be related to the familiar permittivity, ε = n 2 e , and permeability, µ = n 2 m , in the Maxwell-Minkowski formulation of continuum electrodynamics for a simple linear medium. The use of n e and n m instead of √ ε and √ µ, respectively, is a matter of notational convenience.
For linearly polarized radiation normally incident on the boundary between two homogeneous linear media, we can write the vector potential of the incident (i), reflected (r), and refracted (R) waves in terms of the constant amplitudesÃ i ,Ã r , andÃ R of the rectangular/top-hat fields as 2-5 A r =ê xÃr e −i(ω d t+k 1 z) in the plane-wave limit. Here, ω d is the frequency of the field propagating in the direction of the z-axis andê x is a unit polarization vector that is perpendicular to the direction of propagation.
The total energy of a thermodynamically closed system is constant in time by virtue of being conserved. At time t 0 , the entire field is in medium 1 propagating toward the interface with medium 2. At time t 1 , the entire refracted field is in medium 2 and the reflected field is in medium 1. Both the refracted and reflected fields are propagating away from the material interface. The total energy at time t 0 , the incident energy U (t 0 ) = U i , is equal to the total energy at a later time t 1 , U (t 1 ) = U r + U R , which is the sum of the reflected energy U r (t 1 ) and the refracted energy U R (t 1 ). In terms of the incident, reflected, and refracted energy, upon cancelling common constant factors. Here, V 1 is the spatial volume of the incident/reflected field in the vacuum and V 2 is the volume of the refracted field in the medium.
In order to facilitate the integration of Eq. (19), we choose the incident pulse to be rectangular with a nominal width of w i . The refracted pulse in medium 2 has a width of n e1 n m1 w i /(n e2 n m2 ) due to the change in the velocity of light between the two media. Then, evaluating the integrals of Eq. (19) results in Grouping terms of like refractive index, the previous equation becomes n e1 n m1 The second-order equation, Eq. (21), can be written as two first-order equations, but the decomposition is not unique.
In order to derive boundary relations, we need a second linearly independent relation.
Substituting the relations between the vector potential and the fields, Eqs. (13) and (14), into the Maxwell-Ampère law, Eq. (2), we obtain the wave equation that describes the propagation of the electromagnetic field through an optically transparent linear magneto-dielectric medium.
Consider a thin right rectangular box or "Gaussian pillbox" that straddles the interface between the two mediums with the large surfaces parallel to the interface, Fig. 1. Then S is the surface of the pillbox, da is an element of area on the surface, andn is an outwardly directed unit vector that is normal to da. Integrating the wave equation, Eq. (22), we obtain There is no contribution to the surface integral from the large surfaces for a normally incident field in the plane-wave limit because both A and ∇ × (∇ × A) are orthogonal ton. The contributions from the smaller surfaces can be neglected as the box becomes arbitrarily thin.
Applying Stokes's theorem for an arbitrary vector v to Eq. (23), we have We choose the closed contour C in the form of a rectangular Stokesian loop with sides that bisect the two large surfaces and two of the small surfaces on opposite sides of the pillbox shown in Fig. 1. Here, dl is a directed line element that lies on the contour, C. Then C straddles the material interface. For normal incidence in the plane-wave limit, the field n −2 m ∇ × A can be oriented along the long sides of the contour C. Performing the contour integration in Eq. (25), the contribution from the short sides of the contour are neglected as the loop is made vanishingly thin and we obtain that corresponds to the usual electromagnetic boundary conditions 2-5 .
Subsequently, we eliminateÃ r to get For oblique incidence at an angle θ i with the normal to the interface, the width of the refracted pulse is w i (n 1 cos θ R )/(n 2 cos θ i ). The energy balance, Eq. (20) becomes n e1 n m1Ã for our nominally square/rectangular pulse. Then the usual Fresnel relations A r A i = n e1 n m2 cos θ i − n e2 n m1 cos θ R n e1 n m2 cos θ i + n e2 n m1 cos θ R (33) are obtained when we combine the energy balance, Eq. (32), with the relatioñ that is obtained by applying the Stokes's theorem to the wave equation for the case in which the vector potential is perpendicular to the plane of incidence. Likewise, one obtains A r A i = n e2 n m1 cos θ i − n e1 n m2 cos θ R n e2 n m1 cos θ i + n e1 n m2 cos θ R (36) A R A i = 2n e1 n m2 cos θ i n e2 n m1 cos θ i + n e1 n m2 cos θ R (37) by combining the energy balance, Eq. (32), with cos θ i (Ã i −Ã r ) = cos θ RÃR .
that is derived from the wave equation by using the Stokes's theorem in the case of a vector potential that is parallel to the plane of incidence.
In this communication, we obtained the electromagnetic boundary conditions and the