# Physics:Gauss's law for magnetism

Short description: Foundational law of classical magnetism

In physics, Gauss's law for magnetism is one of the four Maxwell's equations that underlie classical electrodynamics. It states that the magnetic field B has divergence equal to zero,[1] in other words, that it is a solenoidal vector field. It is equivalent to the statement that magnetic monopoles do not exist.[2] Rather than "magnetic charges", the basic entity for magnetism is the magnetic dipole. (If monopoles were ever found, the law would have to be modified, as elaborated below.)

Gauss's law for magnetism can be written in two forms, a differential form and an integral form. These forms are equivalent due to the divergence theorem.

The name "Gauss's law for magnetism"[1] is not universally used. The law is also called "Absence of free magnetic poles".[2] It is also referred to as the "transversality requirement"[3] because for plane waves it requires that the polarization be transverse to the direction of propagation.

## Differential form

The differential form for Gauss's law for magnetism is:

$\displaystyle{ \nabla\cdot\mathbf{B} = 0 }$

where ∇ · denotes divergence, and B is the magnetic field.

## Integral form

Definition of a closed surface.
Left: Some examples of closed surfaces include the surface of a sphere, surface of a torus, and surface of a cube. The magnetic flux through any of these surfaces is zero.
Right: Some examples of non-closed surfaces include the disk surface, square surface, or hemisphere surface. They all have boundaries (red lines) and they do not fully enclose a 3D volume. The magnetic flux through these surfaces is not necessarily zero.

The integral form of Gauss's law for magnetism states:

$\displaystyle{ \scriptstyle S }$ $\displaystyle{ \mathbf{B} \cdot \mathrm{d}\mathbf{S} = 0 }$

where S is any closed surface (see image right), and dS is a vector, whose magnitude is the area of an infinitesimal piece of the surface S, and whose direction is the outward-pointing surface normal (see surface integral for more details).

The left-hand side of this equation is called the net flux of the magnetic field out of the surface, and Gauss's law for magnetism states that it is always zero.

The integral and differential forms of Gauss's law for magnetism are mathematically equivalent, due to the divergence theorem. That said, one or the other might be more convenient to use in a particular computation.

The law in this form states that for each volume element in space, there are exactly the same number of "magnetic field lines" entering and exiting the volume. No total "magnetic charge" can build up in any point in space. For example, the south pole of the magnet is exactly as strong as the north pole, and free-floating south poles without accompanying north poles (magnetic monopoles) are not allowed. In contrast, this is not true for other fields such as electric fields or gravitational fields, where total electric charge or mass can build up in a volume of space.

## Vector potential

Main page: Physics:Magnetic vector potential

Due to the Helmholtz decomposition theorem, Gauss's law for magnetism is equivalent to the following statement:[4][5]

There exists a vector field A such that $\displaystyle{ \mathbf{B} = \nabla\times\mathbf{A}. }$

The vector field A is called the magnetic vector potential.

Note that there is more than one possible A which satisfies this equation for a given B field. In fact, there are infinitely many: any field of the form ϕ can be added onto A to get an alternative choice for A, by the identity (see Vector calculus identities): $\displaystyle{ \nabla\times \mathbf{A} = \nabla\times(\mathbf{A} + \nabla \phi) }$ since the curl of a gradient is the zero vector field: $\displaystyle{ \nabla\times \nabla \phi=\boldsymbol{0} }$

This arbitrariness in A is called gauge freedom.

## Field lines

Main page: Field line

The magnetic field B can be depicted via field lines (also called flux lines) – that is, a set of curves whose direction corresponds to the direction of B, and whose areal density is proportional to the magnitude of B. Gauss's law for magnetism is equivalent to the statement that the field lines have neither a beginning nor an end: Each one either forms a closed loop, winds around forever without ever quite joining back up to itself exactly, or extends to infinity.

## Incorporating magnetic monopoles

Main page: Physics:Magnetic monopole

If magnetic monopoles were to be discovered, then Gauss's law for magnetism would state the divergence of B would be proportional to the magnetic charge density ρm, analogous to Gauss's law for electric field. For zero net magnetic charge density (ρm = 0), the original form of Gauss's magnetism law is the result.

The modified formula for use with the SI is not standard and depends on the choice of defining equation for the magnetic charge and current; in one variation, magnetic charge has units of webers, in another it has units of ampere-meters.

System Equation
SI (weber convention)[6] $\displaystyle{ \nabla\cdot\mathbf{B} = \rho_{\mathrm m} }$
SI (ampere-meter convention)[7] $\displaystyle{ \nabla\cdot\mathbf{B} = \mu_0\rho_{\mathrm m} }$
CGS-Gaussian[8] $\displaystyle{ \nabla\cdot\mathbf{B} = 4\pi\rho_{\mathrm m} }$

where μ0 is the vacuum permeability.

So far, examples of magnetic monopoles are disputed in extensive search,[9] although certain papers report examples matching that behavior. [10]

## History

This idea of the nonexistence of the magnetic monopoles originated in 1269 by Petrus Peregrinus de Maricourt. His work heavily influenced William Gilbert, whose 1600 work De Magnete spread the idea further. In the early 1800s Michael Faraday reintroduced this law, and it subsequently made its way into James Clerk Maxwell's electromagnetic field equations.

## Numerical computation

In numerical computation, the numerical solution may not satisfy Gauss's law for magnetism due to the discretization errors of the numerical methods. However, in many cases, e.g., for magnetohydrodynamics, it is important to preserve Gauss's law for magnetism precisely (up to the machine precision). Violation of Gauss's law for magnetism on the discrete level will introduce a strong non-physical force. In view of energy conservation, violation of this condition leads to a non-conservative energy integral, and the error is proportional to the divergence of the magnetic field.[11]

There are various ways to preserve Gauss's law for magnetism in numerical methods, including the divergence-cleaning techniques,[12] the constrained transport method,[13] potential-based formulations[14] and de Rham complex based finite element methods[15][16] where stable and structure-preserving algorithms are constructed on unstructured meshes with finite element differential forms.

## References

1. Chow, Tai L. (2006). Electromagnetic Theory: A modern perspective. Jones and Bartlett. p. 134. ISBN 0-7637-3827-1.
2. Jackson, John David (1999). Classical Electrodynamics (3rd ed.). John Wiley & Sons. p. 237. ISBN 0-471-30932-X.
3. Joannopoulos, John D.; Johnson, Steve G.; Winn, Joshua N.; Meade, Robert D. (2008). Photonic Crystals: Molding the Flow of Light (2nd ed.). Princeton University Press. p. 9. ISBN 978-0-691-12456-8.
4. Schilders, W. H. A. (2005). Handbook of Numerical Analysis. p. 13. ISBN 978-0-444-51375-5.
5. Jackson, John David (1999). Classical Electrodynamics (3rd ed.). John Wiley & Sons. p. 180. ISBN 0-471-30932-X.
6. Jackson, John David (1999). Classical Electrodynamics (3rd ed.). John Wiley & Sons. p. 273, eq. 6.150.
7. See for example equation 4 in Nowakowski, M.; Kelkar, N. G. (2005). "Faraday's law in the presence of magnetic monopoles". Europhysics Letters 71 (3): 346. doi:10.1209/epl/i2004-10545-2. Bibcode2005EL.....71..346N.
8. Moulin, F. (2001). "Magnetic monopoles and Lorentz force". Il Nuovo Cimento B 116 (8): 869–877. Bibcode2001NCimB.116..869M.
9. Magnetic Monopoles, report from Particle data group, updated August 2015 by D. Milstead and E.J. Weinberg. "To date there have been no confirmed observations of exotic particles possessing magnetic charge."
10. Castelnovo, C.; Moessner, R.; Sondhi, S. L. (January 3, 2008). "Magnetic monopoles in spin ice". Nature. 451 (7174): 42–45. arXiv:0710.5515. Bibcode:2008Natur.451...42C. doi:10.1038/nature06433. PMID 18172493. S2CID 2399316.
11. Brackbill, J.U; Barnes, D.C (May 1980). "The Effect of Nonzero ∇ · B on the numerical solution of the magnetohydrodynamic equations". Journal of Computational Physics 35 (3): 426–430. doi:10.1016/0021-9991(80)90079-0. Bibcode1980JCoPh..35..426B.
12. Tóth, Gábor (1 July 2000). "The ∇·B=0 Constraint in Shock-Capturing Magnetohydrodynamics Codes". Journal of Computational Physics 161 (2): 605–652. doi:10.1006/jcph.2000.6519. ISSN 0021-9991. Bibcode2000JCoPh.161..605T.
13. Hernquist, Lars; Vogelsberger, Mark; Mocz, Philip (21 July 2014). "A constrained transport scheme for MHD on unstructured static and moving meshes" (in en). Monthly Notices of the Royal Astronomical Society 442 (1): 43–55. doi:10.1093/mnras/stu865. ISSN 0035-8711. Bibcode2014MNRAS.442...43M.
14. Jardin, Stephen (2010). Computational Methods in Plasma Physics (1st ed.). Boca Raton: CRC Press. ISBN 9780429075537.
15. Hu, Kaibo; Ma, Yicong; Xu, Jinchao (1 February 2017). "Stable finite element methods preserving ∇·B=0 exactly for MHD models" (in en). Numerische Mathematik 135 (2): 371–396. doi:10.1007/s00211-016-0803-4. ISSN 0945-3245.
16. Ma, Yicong; Hu, Kaibo; Hu, Xiaozhe; Xu, Jinchao (July 2016). "Robust preconditioners for incompressible MHD models". Journal of Computational Physics 316: 721–746. doi:10.1016/j.jcp.2016.04.019. Bibcode2016JCoPh.316..721M.