Derivation of the Integral Form of the Momentum Conservation Equation

In classical electrodynamics, the conservation of momentum stands as a foundational principle governing the intricate interplay between matter and electromagnetic fields. Unlike purely mechanical systems, an electromagnetic framework distributes momentum across two distinct domains: the tangible momentum harbored by moving charges (mechanical momentum) and the intrinsic momentum carried by the fields themselves (field momentum). To rigorously analyze these macroscopic interactions, physicists transition from local differential force densities to a global integral formulation of momentum conservation.

Before diving into the mathematical formulation, it is essential to define the primary physical quantities involved:

  • Mechanical Momentum Density ($\mathbf{g}_{mech}$): For a continuous medium with a charge density $\rho_m$ moving at a velocity $\mathbf{v}$, the mechanical momentum density is expressed as:
    $$\mathbf{g}_{mech} = \rho_m \mathbf{v}$$

  • Electromagnetic Field Momentum Density ($\mathbf{g}_{em}$): The momentum stored within the electromagnetic field is fundamentally tied to the Poynting vector $\mathbf{S} = \frac{1}{\mu_0} (\mathbf{E} \times \mathbf{B})$ and the speed of light, yielding:
    $$\mathbf{g}_{em} = \epsilon_0 (\mathbf{E} \times \mathbf{B})$$

  • Lorentz Force Density ($\mathbf{f}_{em}$): The net electromagnetic force exerted per unit volume on a distribution of charges and currents is given by:
    $$\mathbf{f}_{em} = \rho \mathbf{E} + \mathbf{J} \times \mathbf{B}$$
    where $\rho$ represents the free charge density and $\mathbf{J}$ is the current density.

  • Maxwell Stress Tensor ($\mathbf{T}$): A second-order symmetric tensor that encapsulates the transmission of momentum via electromagnetic fields through space. Its components $T_{ij}$ are defined as:
    $$T_{ij} = \epsilon_0 \left( E_i E_j - \frac{1}{2} \delta_{ij} E^2 \right) + \frac{1}{\mu_0} \left( B_i B_j - \frac{1}{2} \delta_{ij} B^2 \right)$$
    The core of the derivation relies on bridging the electromagnetic force density $\mathbf{f}_{em}$ with the divergence of the Maxwell stress tensor $\mathbf{T}$. By manipulating Maxwell's equations alongside vector calculus identities, one can rewrite the Lorentz force density as:

$$\mathbf{f}_{em} = \nabla \cdot \mathbf{T} - \epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t}$$

This differential relation carries profound physical implications:

  • The term $\nabla \cdot \mathbf{T}$ accounts for the net force transmitted across a surface by the electromagnetic field, functioning analogously to internal stress gradients in continuum mechanics.
  • The term $\epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t}$ reflects the temporal rate of change of the electromagnetic field momentum density (noting that $\epsilon_0 \mu_0 = 1/c^2$).

Step-by-Step Derivation of the Integral Form

We now evaluate this local differential relation over an arbitrary fixed control volume $V$ bounded by a closed surface $\partial V$.

Step 1: Connecting Force to Momentum Variation

Invoking the continuum formulation of Newton's second law, the total electromagnetic force acting on the matter enclosed within volume $V$ must equal the temporal rate of change of its total mechanical momentum:
$$\int_V \mathbf{f}{em} , dV = \frac{d}{dt} \int_V \mathbf{g}{mech} , dV \equiv \frac{d\mathbf{P}{mech}}{dt}$$
where $\mathbf{P}{mech}$ denotes the aggregate mechanical momentum inside $V$.

Step 2: Substituting the Stress Tensor Expansion

Next, we substitute the alternative expression for $\mathbf{f}{em}$ into the volume integral:
$$\frac{d\mathbf{P}{mech}}{dt} = \int_V \left( \nabla \cdot \mathbf{T} - \epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t} \right) dV$$

Using the linearity of integration, this separates into two distinct volume integrals:
$$\frac{d\mathbf{P}_{mech}}{dt} = \int_V (\nabla \cdot \mathbf{T}) , dV - \int_V \epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t} , dV$$

Step 3: Applying the Divergence Theorem

We apply the Divergence Theorem (Gauss's theorem) to the first term, converting the volume integral of a divergence into a closed surface integral over the bounding area $\partial V$:
$$\int_V (\nabla \cdot \mathbf{T}) , dV = \oint_{\partial V} \mathbf{T} \cdot d\mathbf{a}$$
where $d\mathbf{a}$ is the outward-pointing differential area vector.

For the second term, assuming a stationary control volume $V$, we can bring the time derivative outside the integral:
$$\int_V \epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t} , dV = \frac{d}{dt} \int_V \epsilon_0 \mu_0 \mathbf{S} , dV$$
Given that $\epsilon_0 \mu_0 \mathbf{S} = \mathbf{g}{em}$, this integral simply evaluates to the rate of change of the total field momentum $\mathbf{P}{em}$ within $V$:
$$\int_V \epsilon_0 \mu_0 \frac{\partial \mathbf{S}}{\partial t} , dV = \frac{d\mathbf{P}_{em}}{dt}$$

Step 4: Synthesizing the Final Equation

Combining these transformations yields:
$$\frac{d\mathbf{P}{mech}}{dt} = \oint{\partial V} \mathbf{T} \cdot d\mathbf{a} - \frac{d\mathbf{P}_{em}}{dt}$$

Rearranging the terms results in the definitive integral form of the momentum conservation equation:
$$\frac{d}{dt} (\mathbf{P}{mech} + \mathbf{P}{em}) = \oint_{\partial V} \mathbf{T} \cdot d\mathbf{a}$$

Physical Interpretation

The resulting balance equation, $\frac{d\mathbf{P}{total}}{dt} = \oint{\partial V} \mathbf{T} \cdot d\mathbf{a}$, offers vital insights into electrodynamic momentum dynamics:

  1. Total Momentum Accounting: The left-hand side monitors the time derivative of the combined system momentum—spanning both mechanical and field contributions.
  2. Momentum Flux: The right-hand surface integral signifies the net flux of momentum crossing the boundary $\partial V$. For an isolated system where no electromagnetic stresses cross the distant boundary ($\mathbf{T} \cdot d\mathbf{a} = 0$), the total momentum $\mathbf{P}_{total}$ remains strictly conserved.
  3. Dynamic Exchange: The equation emphasizes that variations in mechanical momentum are not solely driven by instantaneous local forces, but must also account for momentum stored and transported by the fields. This continuous exchange mechanism forms the bedrock for analyzing phenomena like radiation pressure and plasma confinement.

Practical Illustration: Radiation Pressure

Consider an intense laser beam striking a reflective mirror. The incident light carries a well-defined electromagnetic field momentum $\mathbf{P}{em}$. Upon reflection, the beam undergoes a momentum shift $\Delta \mathbf{P}{em}$. To satisfy overall momentum conservation as dictated by our integral equation, the mirror must absorb an equal and opposite change in mechanical momentum $\Delta \mathbf{P}_{mech}$. Evaluating the Maxwell stress tensor across the boundary enclosing the mirror allows physicists to calculate this photon-induced force, commonly known as radiation pressure, with high precision.