Magnetic Field Boundary Conditions in the Monopole Viewpoint

In classical electrodynamics, Maxwell's equations are traditionally formulated under the assumption that magnetic monopoles do not exist. Consequently, the magnetic flux through any closed surface is invariably zero, mathematically expressed as $\nabla \cdot \mathbf{B} = 0$. However, from the perspective of theoretical symmetry, introducing the concept of "magnetic charge" reveals a profound elegance in electromagnetic theory. By adopting the monopole viewpoint, the magnetic field is no longer viewed solely as a product of electric currents; it can also arise from the motion of magnetic charges (magnetic currents) or their static distribution. This paradigm shift not only deepens our understanding of the fundamental nature of electromagnetic fields but also provides a unified mathematical framework for deriving boundary conditions at interfaces between different media.

To construct a theoretical system based on magnetic charges, we must extend Maxwell's equations. We hypothesize the existence of a magnetic charge density $\rho_m$ (measured in $\text{Wb}/\text{m}^3$) and a magnetic current density $\mathbf{J}_m$ (measured in $\text{A}/\text{m}$ or $\text{Wb}/(\text{m}^2\cdot\text{s})$). The expanded Maxwell's equations, presented in both integral and differential forms, are as follows:

  1. Gauss's Law for Magnetism (Extended):
    $$\nabla \cdot \mathbf{B} = \mu_0 \rho_m$$
    This equation signifies that the divergence of the magnetic field is no longer zero but is determined by the local magnetic charge density.

  2. Faraday's Law of Induction (Extended):
    $$\nabla \times \mathbf{E} = -\mu_0 \mathbf{J}_m - \frac{\partial \mathbf{B}}{\partial t}$$
    Here, the magnetic current $\mathbf{J}_m$ contributes to the curl of the electric field, illustrating the induction of an electric field generated by moving magnetic charges.

  3. Gauss's Law for Electricity:
    $$\nabla \cdot \mathbf{D} = \rho_e$$

  4. Ampère-Maxwell Law:
    $$\nabla \times \mathbf{H} = \mathbf{J}_e + \frac{\partial \mathbf{D}}{\partial t}$$

Through this extension, the mathematical forms of the electric and magnetic fields achieve a high degree of symmetry, treating electric and magnetic sources on an equal footing.

Derivation of Magnetic Field Boundary Conditions

When a magnetic field transitions from one medium (Medium 1) to another (Medium 2), specific continuity conditions must be satisfied at the interface. By applying integral theorems—specifically Gauss's theorem and Stokes' theorem—to the expanded equation set, we can derive these conditions rigorously.

1. Normal Component Boundary Condition for Magnetic Flux Density $\mathbf{B}$

Consider an infinitesimal "pillbox" straddling the interface, with a base area $A$ and height $\Delta h \to 0$. According to the integral form of the extended Gauss's law for magnetism:
$$\oint_S \mathbf{B} \cdot d\mathbf{A} = \mu_0 Q_{m, \text{encl}}$$

As the height approaches zero, the contribution from the side walls becomes negligible, leaving only the flux through the top and bottom faces:
$$(B_{n1} - B_{n2}) A = \mu_0 (\sigma_m A)$$
where $\sigma_m$ represents the magnetic surface charge density distributed on the interface. Dividing by the area $A$, we obtain:
$$B_{n1} - B_{n2} = \mu_0 \sigma_m$$

Conclusion: The discontinuity in the normal component of the magnetic flux density $\mathbf{B}$ is directly proportional to the magnetic surface charge density $\sigma_m$ at the interface. If no magnetic charge is present ($\sigma_m = 0$), then $B_{n1} = B_{n2}$, maintaining continuity.

2. Tangential Component Boundary Condition for Magnetic Field Intensity $\mathbf{H}$

Next, consider an infinitesimal "loop" straddling the interface with length $L$ and width $w \to 0$. While the extended Faraday's law governs the electric field, symmetry dictates that the behavior of $\mathbf{H}$ follows a form analogous to Ampère's law. In the monopole viewpoint, a magnetic surface current density $\mathbf{K}m$ on the interface induces a jump in the tangential component of $\mathbf{H}$:
$$\oint_C \mathbf{H} \cdot d\mathbf{l} = I{m, \text{encl}}$$

As the width of the loop vanishes, the line integral simplifies to:
$$(H_{t1} - H_{t2}) L = K_m L$$
which yields the vector form:
$$\mathbf{n} \times (\mathbf{H}_1 - \mathbf{H}_2) = \mathbf{K}_m$$

Conclusion: The discontinuity in the tangential component of the magnetic field intensity $\mathbf{H}$ is equal to the magnetic surface current density $\mathbf{K}_m$. In the absence of magnetic surface currents ($\mathbf{K}_m = 0$), the tangential component of $\mathbf{H}$ remains continuous across the boundary.

Summary: Symmetric Comparison of Electromagnetic Boundary Conditions

By adopting the monopole viewpoint, the boundary conditions for electric and magnetic fields can be organized into a symmetric table. This structure is invaluable for grasping the physical essence and facilitating quick recall:

Physical Quantity Component Boundary Condition (With Sources) Boundary Condition (Source-Free)
Electric Field $\mathbf{E}$ Tangential $\mathbf{E}_t$ $\mathbf{n} \times (\mathbf{E}_1 - \mathbf{E}_2) = \mathbf{J}_e \cdot \mathbf{n} \times \mathbf{l}$ (Generalized) $E_{t1} = E_{t2}$
Electric Field $\mathbf{E}$ Normal $\mathbf{E}_n$ $D_{n1} - D_{n2} = \sigma_e$ $D_{n1} - D_{n2} = 0$
Magnetic Field $\mathbf{H}$ Tangential $\mathbf{H}_t$ $\mathbf{n} \times (\mathbf{H}_1 - \mathbf{H}_2) = \mathbf{K}_m$ $H_{t1} = H_{t2}$
Magnetic Field $\mathbf{B}$ Normal $\mathbf{B}_n$ $B_{n1} - B_{n2} = \mu_0 \sigma_m$ $B_{n1} = B_{n2}$

Practical Applications and Implications

Although macroscopic engineering applications typically assume the non-existence of magnetic monopoles (effectively setting $\sigma_m = 0$ and $\mathbf{K}_m = 0$), the monopole viewpoint holds significant theoretical value in several contexts:

  1. Equivalent Modeling of Magnetic Materials: In materials with strong anisotropy or complex domain structures, "effective magnetic charges" are often employed to simplify the calculation of magnetic field distributions.
  2. Theoretical Physics Research: In Quantum Electrodynamics (QED) or the study of topological insulators, the excitation of magnetic monopoles is a critical phenomenon. Accurate modeling in these domains necessitates the use of boundary conditions incorporating magnetic charge terms.
  3. Numerical Simulation Verification: When developing electromagnetic simulation software (such as Finite Element Analysis tools), utilizing the symmetric Maxwell equations with monopole terms enhances the robustness and generality of the algorithms.

Example Problem:
Consider a planar interface dividing space into Region 1 and Region 2, covered by a layer of magnetic surface charge density $\sigma_m$. If the magnetic flux density in Region 1 is known to be $\mathbf{B}_1 = 2\mathbf{a}z$, determine the normal component of the magnetic flux density $B{n2}$ in Region 2.

Solution:
Based on the normal boundary condition derived for the monopole viewpoint:
$$B_{n1} - B_{n2} = \mu_0 \sigma_m$$
Substituting the known values:
$$2 - B_{n2} = \mu_0 \sigma_m \implies B_{n2} = 2 - \mu_0 \sigma_m$$
This result demonstrates that the presence of magnetic charge directly causes a discontinuity in the normal component of the magnetic flux density, fundamentally altering the field distribution across the interface.