- 24,488
- 15,057
Charles Link said:The magnetic surface currents are present even in the static case, and occur also in a permanent magnet. The magnetic field for any magnet can be computed from Biot-Savart (or Ampere's law) if the magnetization ## \vec{M} ## is known, so that the surface currents can be computed. (In complete detail, the magnetic current density ## \vec{J}_m=\frac{\nabla \times \vec{M}}{\mu_o} ##. The result is that there are surface currents at boundaries given by surface current per unit length ## \vec{K}_m=\frac{\vec{M} \times \hat{n}}{\mu_o} ##, but there can also be magnetic current densities in the bulk material in regions of non-uniform ## \vec{M} ##). ## \\ ## For a somewhat detailed example of the magnetic surface currents, see https://www.physicsforums.com/threads/magnetic-field-of-a-ferromagnetic-cylinder.863066/
[Later note:] There seems to be varying conventions in the physics and electrical-engineering literature: I've just checked Simonyi's comprehensive textbook on "theoretical electric engineering" (the title of my German edition is "Theoretische Elektrotechnik", I don't know, whether this amazing book is available in English; it's simply great!) for another question preparing my lecture notes. There I stumbled over the surprising fact that he defines the magnetization as part of ##\vec{B}##. As I write below, for me as a physicist used to interpreting everything from the point of view of the fundamental principles, the magnetization should clearly belong to ##\vec{H}##, because it's an external source for the magnetic field ##\vec{B}##. There seems to be more historical balast in E&M than I ever dreamt of ;-))). Of course, at the end both conventions lead to the correct result, the only difference is that ##\vec{M}_{\text{engineers}}=\mu_0 \vec{M}_{\text{physicists}}##, i.e., the constitutive equation reads
##\vec{B}=\mu_0 \vec{H} + \vec{M}_{\text{engineers}} \equiv \mu_0 (\vec{H}+\vec{M}_{\text{physicists}}## ;-)).
Just a hint. I'm teaching E&M right now to high school-teacher students, and there I've of course to stick to the SI units. So I'm a bit thinking how to explain to the students the somewhat idosyncratic terminology of permittivities and permeabilities occurring in opposite ways (and even occurring for the vacuum just because of the SI units), because in the early days of electromagnetism the physicists somehow mixed up ##\vec{H}## and ##\vec{B}##. It's very important not to confuse where the ##\mu_0## (or ##\mu=\mu_0 \mu_{\text{r}}## in para- and diamagnetic substances) occur.
For hard ferromagnetics you only give the magnetization (magnetic-dipole density per unit volume) and assume full saturation. The magnetization belongs to the "excitation quantities", i.e., to the "magnetic excitation" ##\vec{H}## rather than the "magnetic field" ##\vec{B}##, where I use the modern naming scheme following from the physical meaning and relativity, where ##\vec{E}/c## and ##\vec{B}## together form an antisymmetric four-tensor as well as ##\vec{D}## and ##\vec{H}##. Physically ##\vec{E}## and ##\vec{B}## are the fields occurring in the Lorentz force law, while ##\vec{D}## and ##\vec{H}## are auxilliary quantities to separate the external (or "free") sources from the sources due to the presence of matter from bound charges (i.e. the electric and magnetic dipole moments induced by the perturbation through the fields produced by the free charges and current densities).
Thus for hard ferromagnetics the basic static equations read
$$\vec{\nabla} \times \vec{H}=\vec{j}_f, \quad \vec{\nabla} \cdot \vec{D}=\rho_f, \quad \vec{\nabla} \times \vec{E}=0, \quad \vec{\nabla} \cdot \vec{B}=0$$
and the constitutive equations
$$\vec{D}=\epsilon_0 \epsilon_r \vec{E}, \quad \vec{B}=\mu_0 (\vec{H}+\vec{M}),$$
where ##\vec{M}## is the magnetization unchanged by external influence due to the presence of hard ferromagnets. The electric and magnetic fields decouple completely, and the electric part is just as in electrostatics. So we can concentrate on the magnetic fields.
In this most general case for static fields, i.e., with currents in addition to permanent magnetization present, the most general way to solve the equations is to introduce the vector potential for ##\vec{B}##. For the case without currents, ##\vec{H}## is curl-free and you can simplify the task somewhat by introducing a magnetostatic potential for ##\vec{H}##.
The first method, which applies everywhere in the most general static case, you get formally vacuum equations for ##\vec{B}## with an effective current composed by the free current and the curl of the magnetization. This becomes clear by taking the curl of the 2nd constitutive equation and using Ampere's Law for ##\vec{H}##:
$$\vec{\nabla} \times \vec{B}=\mu_0 \vec{\nabla} \times (\vec{H}+\vec{M}) = \mu_0 (\vec{j}_e + \vec{j}_m), \quad \vec{j}_m=\vec{\nabla} \times \vec{M}.$$
Here ##\vec{j}_e## and ##\vec{j}_m## are usually considered to be given. If one assumes ##\vec{M}## to be constant within the magnet and (of course) jumps instantaneously to 0 at the boundary of the magnet, its curl leads to a singular ##\delta##-function like ##\vec{j}_m## which is equivalent to a surface current along the boundary surface of the magnet. This is in some sense in accordance with Ampere's old molecular-current hypothesis, which of course today must be substituted by the proper quantum-mechanical treatment. For the phenomenology it's just a calculational trick to substitute the approximate analysis with homogeneous magnetization within the magnet by the equivalent surface current given also in this case by ##\vec{\nabla} \times \vec{M}##, leading however to a ##\delta##-distribution due to the jump at the boundary.
In any case the solution lies in the introduction of the vector potential. Due to gauge freedom we can also impose a gauge constraint, which in this static case is most conveniently chosen as the Coulomb gauge, i.e., we set
$$\vec{B}=\vec{\nabla} \times \vec{A}, \quad \vec{\nabla} \cdot \vec{A}=0,$$
and for Cartesian components we can thus write
$$\vec{\nabla} \times \vec{B}=\vec{\nabla} \times (\vec{\nabla} \times \vec{A})=-\Delta \vec{A}=\mu_0 \vec{j}, \quad \vec{j}=\vec{j}_f + \vec{\nabla} \times \vec{M}$$
with the Biot-Savart-like solution
$$\vec{A}(\vec{r}) = \mu_0 \int_{\mathbb{R}^3} \mathrm{d}^3 r' \frac{\vec{j}(\vec{r}')}{|\vec{r}-\vec{r}'|},$$
or for the field
$$\vec{B}(\vec{r}) =\nabla \times \vec{A}(\vec{r})=\mu_0 \int_{\mathbb{R}^3} \mathrm{d^3} r' \vec{j}(\vec{r}') \times \frac{\vec{r}-\vec{r}'}{|\vec{r}-\vec{r}'|^3}.$$
This is Biot-Savart's Law in its usual form but with the total current density including the effective magnetization current.
If somebody is interested in my German manuscript for these lectures (unfinished, but quickly growing ;-)), here they are:
https://th.physik.uni-frankfurt.de/~hees/publ/theo2-l3.pdf
Last edited: