Corrected Energy-Momentum Conservation in Inhomogeneous Media

In the classical description of electromagnetism, the conservation of energy and momentum in vacuum is expressed through elegant, symmetric equations. When a field propagates in a uniform medium, the same formalism can be adapted with only minor adjustments. However, once the medium becomes inhomogeneous—its permittivity (\epsilon(\mathbf{r})) and permeability (\mu(\mathbf{r})) vary with position—the underlying symmetries that guarantee simple conservation laws are broken. The result is a richer, more intricate set of equations that must capture how the electromagnetic field exchanges energy and momentum with the material background.

This article presents a concise, yet comprehensive, derivation of the corrected energy–momentum conservation laws for linear, lossless, inhomogeneous media. It highlights the role of spatial gradients of material parameters, explains the physical origin of the additional terms, and illustrates their significance with a familiar experimental context: optical tweezers.

Revisiting Poynting’s Theorem in a Spatially Varying Medium

Constitutive Relations

For a linear, isotropic, non‑dissipative medium, the displacement and magnetic induction fields are defined as

[
\mathbf{D}(\mathbf{r},t)=\epsilon(\mathbf{r}),\mathbf{E}(\mathbf{r},t),\qquad
\mathbf{B}(\mathbf{r},t)=\mu(\mathbf{r}),\mathbf{H}(\mathbf{r},t).
]

Unlike in a homogeneous medium, (\epsilon) and (\mu) now depend on the spatial coordinate (\mathbf{r}).

Differential Form of Energy Conservation

Starting from Maxwell’s equations and the standard vector identity, one obtains

[
\nabla!\cdot(\mathbf{E}\times\mathbf{H})=
-\mathbf{E}!\cdot!\frac{\partial\mathbf{D}}{\partial t}
-\mathbf{H}!\cdot!\frac{\partial\mathbf{B}}{\partial t}.
]

This relation remains formally valid even when (\epsilon) and (\mu) vary. The subtlety lies in the definition of the local energy density:

[
u(\mathbf{r},t)=\frac{1}{2}\bigl(\mathbf{E}!\cdot!\mathbf{D}+\mathbf{H}!\cdot!\mathbf{B}\bigr)
=\frac{1}{2}\bigl(\epsilon(\mathbf{r})E^{2}+\mu(\mathbf{r})H^{2}\bigr).
]

Because (\epsilon) and (\mu) are functions of position, the time derivative of (u) contains extra terms:

[
\frac{\partial u}{\partial t}
=\mathbf{E}!\cdot!\frac{\partial\mathbf{D}}{\partial t}
+\mathbf{H}!\cdot!\frac{\partial\mathbf{B}}{\partial t}
+\frac{1}{2}E^{2}\frac{\partial\epsilon}{\partial t}
+\frac{1}{2}H^{2}\frac{\partial\mu}{\partial t}.
]

In a static medium ((\partial\epsilon/\partial t=\partial\mu/\partial t=0)), the last two terms vanish, and the usual Poynting theorem follows. When the material parameters vary in space, the gradient of the energy density acquires a contribution from the spatial variation of (\epsilon) and (\mu). Consequently, the flux of electromagnetic energy—captured by the Poynting vector (\mathbf{S}=\mathbf{E}\times\mathbf{H})—is no longer simply the transport of a scalar quantity; it is also reshaped by the medium’s internal structure.

Momentum Conservation and the Maxwell Stress Tensor

From Symmetry to Force Density

In a homogeneous medium, the momentum density of the field is (\mathbf{g}=\mathbf{S}/c^{2}), and the force density acting on matter is given by the divergence of the Maxwell stress tensor (\mathbf{T}):

[
\mathbf{f}=\nabla!\cdot!\mathbf{T}-\frac{\partial\mathbf{g}}{\partial t}.
]

The tensor (\mathbf{T}) in vacuum is

[
T_{ij}= \epsilon_{0}!\left(E_{i}E_{j}-\frac{1}{2}\delta_{ij}E^{2}\right)
+\frac{1}{\mu_{0}}!\left(H_{i}H_{j}-\frac{1}{2}\delta_{ij}H^{2}\right).
]

In an inhomogeneous medium, the same expression can be used, but the material dependence of (\epsilon) and (\mu) introduces additional terms when taking the divergence.

Gradient Force in Linear Media

Consider a non‑magnetic, isotropic dielectric ((\mu=\mu_{0})). The force density simplifies to

[
\mathbf{f}\approx
\rho,\mathbf{E}
+(\mathbf{P}!\cdot!\nabla)\mathbf{E}
-\frac{1}{2}E^{2}\nabla\epsilon.
]

The first two terms are familiar: the Lorentz force on free charges and the action of the polarization field. The third term is the new player in an inhomogeneous setting. It represents a gradient force that pulls the material toward regions of higher or lower permittivity, depending on the sign of (\nabla\epsilon). Physically, this arises because the electric field energy density (\frac{1}{2}\epsilon E^{2}) is locally higher where (\epsilon) is larger; the field therefore exerts a pressure that tends to push the material toward lower‑energy configurations.

Physical Interpretation of the Gradient Terms

  1. Energy Redistribution
    When a wave traverses a region where (\epsilon) changes, the stored electric energy per unit volume changes accordingly. To conserve the total energy, the field must do work on the medium, effectively transferring energy from high‑(\epsilon) zones to low‑(\epsilon) zones (or vice versa). This process is encoded in the (-\frac{1}{2}E^{2}\nabla\epsilon) term.

  2. Momentum Exchange
    The spatial variation of material parameters breaks translational symmetry. Noether’s theorem tells us that the corresponding conserved quantity—linear momentum—is no longer strictly conserved for the field alone. Instead, momentum is exchanged with the medium via the gradient force. The field’s momentum density (\mathbf{g}) is altered by the presence of (\nabla\epsilon) and (\nabla\mu).

  3. Directional Bias in Energy Flow
    In gradient‑index (GRIN) optics, the refractive index profile is engineered to steer light. The modified Poynting theorem shows that the direction of the energy flux can be bent by the spatial gradients of (\epsilon) and (\mu), leading to phenomena such as self‑focusing or beam steering without external lenses.

Case Study: Optical Tweezers

Setup

An optical tweezer system typically uses a tightly focused laser beam to trap a microscopic particle (e.g., a polystyrene bead or a living cell) suspended in a fluid. The particle’s refractive index differs from that of the surrounding medium, creating a sharp interface where (\epsilon) changes abruptly.

Gradient Force in Action

The intensity profile of a Gaussian beam falls off radially, producing a strong gradient in the electric field magnitude. Combined with the permittivity discontinuity at the particle surface, the gradient force

[
\mathbf{f}_{\text{grad}}\propto -\nabla\epsilon,E^{2}
]

directs the particle toward the beam’s focal point. This is the optical trapping effect: the particle experiences a restoring force that confines it in three dimensions.

Why the Gradient Term Matters

If one ignored the (\nabla\epsilon) contribution, the force would be predicted solely from the Lorentz term, which cannot explain the stable trapping observed experimentally. The gradient term provides the necessary stabilizing component that balances the scattering force pushing the particle along the beam direction. Thus, a correct description of momentum transfer in an inhomogeneous medium is essential for predicting and optimizing optical tweezer performance.

Implications for Emerging Photonic Technologies

  1. Metamaterials
    Engineered composites often exhibit spatially varying effective (\epsilon) and (\mu). The corrected conservation laws guide the design of negative‑index or hyperbolic metamaterials, ensuring that energy and momentum flow as intended.

  2. Non‑Linear Optics
    In media where the refractive index depends on the field intensity, the effective (\epsilon(\mathbf{r})) becomes a dynamic function. The gradient terms then couple to the field evolution, influencing phenomena like self‑phase modulation and soliton formation.

  3. Nanophotonics
    At sub‑wavelength scales, material inhomogeneities are ubiquitous. Accurate force calculations are critical for optical manipulation of nanoparticles, plasmonic trapping, and the development of optical MEMS devices.

Summary

  • Energy Conservation: Poynting’s theorem remains valid, but the energy density must include the spatial dependence of (\epsilon) and (\mu). The resulting flux of energy is modified by gradients of material parameters.
  • Momentum Conservation: The Maxwell stress tensor still governs the force density, yet an additional gradient force term (-\frac{1}{2}E^{2}\nabla\epsilon) appears, reflecting the exchange of momentum between field and medium.
  • Physical Consequences: These corrections explain how electromagnetic waves can exert forces on inhomogeneous media, enabling practical applications such as optical tweezers and the design of advanced photonic devices.
  • Broader Impact: Understanding and applying the corrected conservation laws is indispensable for the accurate modeling of light–matter interaction in any system where material properties vary in space.

By incorporating the spatial gradients of electromagnetic parameters into the conservation equations, we achieve a complete, physically faithful description of energy and momentum flow in inhomogeneous media—an essential foundation for both fundamental research and technological innovation in modern optics.