The magnetic field due to a steady current I I I in a wire element d l d\mathbf{l} d l :
d B = μ 0 I 4 π d l × r ^ r 2 d\mathbf{B} = \frac{\mu_0 I}{4\pi} \frac{d\mathbf{l} \times \hat{\mathbf{r}}}{r^2} d B = 4 π μ 0 I r 2 d l × r ^
For a complete circuit:
B ( r ) = μ 0 I 4 π ∮ d l × r ^ " ∣ r − r ′ ∣ 2 \mathbf{B}(\mathbf{r}) = \frac{\mu_0 I}{4\pi} \oint \frac{d\mathbf{l} \times \hat{\mathbf{r}}"}{|\mathbf{r} - \mathbf{r}'|^2} B ( r ) = 4 π μ 0 I ∮ ∣ r − r ′ ∣ 2 d l × r ^ "
For steady currents (∂ E / ∂ t = 0 \partial \mathbf{E} / \partial t = 0 ∂ E / ∂ t = 0 ):
∮ C B ⋅ d l = μ 0 I e n c \oint_C \mathbf{B} \cdot d\mathbf{l} = \mu_0 I_{\mathrm{enc}} ∮ C B ⋅ d l = μ 0 I enc
Example: Infinite straight wire carrying current I I I .
By cylindrical symmetry, B B B is constant on circles centred on the wire. Choose an Amperian loop of Radius r r r :
∮ B ⋅ d l = B ⋅ 2 π r = μ 0 I ⟹ B = μ 0 I 2 π r \oint \mathbf{B} \cdot d\mathbf{l} = B \cdot 2\pi r = \mu_0 I \implies B = \frac{\mu_0 I}{2\pi r} ∮ B ⋅ d l = B ⋅ 2 π r = μ 0 I ⟹ B = 2 π r μ 0 I
Example: Solenoid. For a long solenoid with n n n turns per unit length carrying current I I I :
B = μ 0 n I ( i n s i d e ) , B = 0 ( o u t s i d e ) B = \mu_0 n I \quad \mathrm{(inside)}, \quad B = 0 \quad \mathrm{(outside)} B = μ 0 n I ( inside ) , B = 0 ( outside )
Since ∇ ⋅ B = 0 \nabla \cdot \mathbf{B} = 0 ∇ ⋅ B = 0 We can write B = ∇ × A \mathbf{B} = \nabla \times \mathbf{A} B = ∇ × A Where A \mathbf{A} A is the magnetic vector potential .
In the Coulomb gauge (∇ ⋅ A = 0 \nabla \cdot \mathbf{A} = 0 ∇ ⋅ A = 0 ), the vector potential satisfies
∇ 2 A = − μ 0 J \nabla^2 \mathbf{A} = -\mu_0 \mathbf{J} ∇ 2 A = − μ 0 J
This is Poisson’s equation for each component of A \mathbf{A} A .
For a current loop, the solution is:
A ( r ) = μ 0 4 π ∫ J ( r ′ ) ∣ r − r ′ ∣ d 3 r ′ \mathbf{A}(\mathbf{r}) = \frac{\mu_0}{4\pi} \int \frac{\mathbf{J}(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|}\, d^3\mathbf{r}' A ( r ) = 4 π μ 0 ∫ ∣ r − r ′ ∣ J ( r ′ ) d 3 r ′
Example: Toroid. A toroid with N N N turns carrying current I I I has inner radius a a a and outer Radius b b b .
By symmetry, B \mathbf{B} B is tangential and constant on circular Amperian loops inside the Toroid. For a loop of radius r r r (a < r < b a \lt r \lt b a < r < b ):
B ⋅ 2 π r = μ 0 N I ⟹ B = μ 0 N I 2 π r B \cdot 2\pi r = \mu_0 N I \implies B = \frac{\mu_0 N I}{2\pi r} B ⋅ 2 π r = μ 0 N I ⟹ B = 2 π r μ 0 N I
For r < a r \lt a r < a or r > b r \gt b r > b : B = 0 B = 0 B = 0 (no enclosed current).
Unlike a solenoid, the field inside a toroid is not uniform --- it varies as 1 / r 1/r 1/ r .
Example: Infinite current sheet. A sheet in the x y xy x y -plane carries surface current density K = K x ^ \mathbf{K} = K\,\hat{\mathbf{x}} K = K x ^ .
By symmetry, B \mathbf{B} B is parallel to ± y ^ \pm\hat{\mathbf{y}} ± y ^ and depends only on z z z . Choose a rectangular Amperian loop straddling the sheet with sides parallel to y ^ \hat{\mathbf{y}} y ^ :
B ⋅ 2 L = μ 0 K L ⟹ B = μ 0 K 2 B \cdot 2L = \mu_0 K L \implies B = \frac{\mu_0 K}{2} B ⋅ 2 L = μ 0 K L ⟹ B = 2 μ 0 K
The field is uniform on each side, pointing in opposite directions:
B = { + μ 0 K 2 y ^ z > 0 − μ 0 K 2 y ^ z < 0 \mathbf{B} = \begin{cases} +\frac{\mu_0 K}{2}\,\hat{\mathbf{y}} & z \gt 0 \\[4pt] -\frac{\mu_0 K}{2}\,\hat{\mathbf{y}} & z \lt 0 \end{cases} B = { + 2 μ 0 K y ^ − 2 μ 0 K y ^ z > 0 z < 0
A current loop carrying current I I I enclosing area a \mathbf{a} a has magnetic dipole moment :
m = I a \mathbf{m} = I\mathbf{a} m = I a
For a planar loop of N N N turns: m = N I A n ^ \mathbf{m} = NIA\,\hat{\mathbf{n}} m = N I A n ^ Where A A A is the area And n ^ \hat{\mathbf{n}} n ^ is the unit normal given by the right-hand rule.
Field of a magnetic dipole (at position r \mathbf{r} r from the dipole):
B d i p ( r ) = μ 0 4 π [ 3 ( m ⋅ r ^ ) r ^ − m r 3 ] \mathbf{B}_{\mathrm{dip}(\mathbf{r}) = \frac{\mu_0}{4\pi}\left[\frac{3(\mathbf{m} \cdot \hat{\mathbf{r}})\hat{\mathbf{r}} - \mathbf{m}}{r^3}\right]} B dip ( r ) = 4 π μ 0 [ r 3 3 ( m ⋅ r ^ ) r ^ − m ]
This has the same angular structure as the electric dipole field.
Torque on a dipole in a uniform field:
τ = m × B \boldsymbol{\tau} = \mathbf{m} \times \mathbf{B} τ = m × B
Energy of a dipole in a field:
U = − m ⋅ B U = -\mathbf{m} \cdot \mathbf{B} U = − m ⋅ B
Force on a dipole in a non-uniform field:
F = ∇ ( m ⋅ B ) \mathbf{F} = \nabla(\mathbf{m} \cdot \mathbf{B}) F = ∇ ( m ⋅ B )
Example: Field on the axis of a circular loop A circular loop of radius R R R carries current I I I . On the axis at distance z z z from the centre, Every element d l d\mathbf{l} d l is perpendicular to r ^ \hat{\mathbf{r}} r ^ So:
d B = μ 0 I 4 π d l R 2 + z 2 d\mathbf{B} = \frac{\mu_0 I}{4\pi}\frac{dl}{R^2 + z^2} d B = 4 π μ 0 I R 2 + z 2 d l
The component perpendicular to the axis cancels by symmetry. The axial component is:
B z = ∮ d B sin α = μ 0 I 4 π ( R 2 + z 2 ) R R 2 + z 2 ∮ d l = μ 0 I R 2 2 ( R 2 + z 2 ) 3 / 2 B_z = \oint dB\,\sin\alpha = \frac{\mu_0 I}{4\pi(R^2+z^2)}\frac{R}{\sqrt{R^2+z^2}}\oint dl = \frac{\mu_0 I R^2}{2(R^2+z^2)^{3/2}} B z = ∮ d B sin α = 4 π ( R 2 + z 2 ) μ 0 I R 2 + z 2 R ∮ d l = 2 ( R 2 + z 2 ) 3/2 μ 0 I R 2
For z ≫ R z \gg R z ≫ R : B z ≈ μ 0 I R 2 2 z 3 = μ 0 4 π 2 m z 3 B_z \approx \frac{\mu_0 I R^2}{2z^3} = \frac{\mu_0}{4\pi}\frac{2\mathbf{m}}{z^3} B z ≈ 2 z 3 μ 0 I R 2 = 4 π μ 0 z 3 2 m Which matches the dipole formula with m = I π R 2 z ^ \mathbf{m} = I\pi R^2\,\hat{\mathbf{z}} m = I π R 2 z ^ . ■ \blacksquare ■
Starting from the Biot-Savart law and the identity r − r ′ ∣ r − r ′ ∣ 3 = − ∇ 1 ∣ r − r ′ ∣ \frac{\mathbf{r} - \mathbf{r}'}{|\mathbf{r}-\mathbf{r}'|^3} = -\nabla\frac{1}{|\mathbf{r}-\mathbf{r}'|} ∣ r − r ′ ∣ 3 r − r ′ = − ∇ ∣ r − r ′ ∣ 1 :
B ( r ) = μ 0 4 π ∫ J ( r ′ ) × ( r − r ′ ) ∣ r − r ′ ∣ 3 d 3 r ′ = − μ 0 4 π ∫ J ( r ′ ) × ∇ 1 ∣ r − r ′ ∣ d 3 r ′ \mathbf{B}(\mathbf{r}) = \frac{\mu_0}{4\pi}\int \mathbf{J}(\mathbf{r}') \times \frac{(\mathbf{r}-\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|^3}\,d^3\mathbf{r}' = -\frac{\mu_0}{4\pi}\int \mathbf{J}(\mathbf{r}') \times \nabla\frac{1}{|\mathbf{r}-\mathbf{r}'|}\,d^3\mathbf{r}' B ( r ) = 4 π μ 0 ∫ J ( r ′ ) × ∣ r − r ′ ∣ 3 ( r − r ′ ) d 3 r ′ = − 4 π μ 0 ∫ J ( r ′ ) × ∇ ∣ r − r ′ ∣ 1 d 3 r ′
Using the product rule J × ( ∇ f ) = ∇ × ( f J ) − f ( ∇ × J ) \mathbf{J} \times (\nabla f) = \nabla \times (f\mathbf{J}) - f(\nabla \times \mathbf{J}) J × ( ∇ f ) = ∇ × ( f J ) − f ( ∇ × J ) And noting that ∇ × J ( r ′ ) = 0 \nabla \times \mathbf{J}(\mathbf{r}') = 0 ∇ × J ( r ′ ) = 0 (since J \mathbf{J} J depends on r ′ \mathbf{r}' r ′ Not r \mathbf{r} r ):
B ( r ) = μ 0 4 π ∇ × ∫ J ( r ′ ) ∣ r − r ′ ∣ d 3 r ′ \mathbf{B}(\mathbf{r}) = \frac{\mu_0}{4\pi}\nabla \times \int \frac{\mathbf{J}(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\,d^3\mathbf{r}' B ( r ) = 4 π μ 0 ∇ × ∫ ∣ r − r ′ ∣ J ( r ′ ) d 3 r ′
Comparing with B = ∇ × A \mathbf{B} = \nabla \times \mathbf{A} B = ∇ × A :
A ( r ) = μ 0 4 π ∫ J ( r ′ ) ∣ r − r ′ ∣ d 3 r ′ \mathbf{A}(\mathbf{r}) = \frac{\mu_0}{4\pi}\int \frac{\mathbf{J}(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\,d^3\mathbf{r}' A ( r ) = 4 π μ 0 ∫ ∣ r − r ′ ∣ J ( r ′ ) d 3 r ′
This is the general solution for the vector potential in the Coulomb gauge. For a line current:
A ( r ) = μ 0 I 4 π ∮ d l ′ ∣ r − r ′ ∣ \mathbf{A}(\mathbf{r}) = \frac{\mu_0 I}{4\pi}\oint \frac{d\mathbf{l}'}{|\mathbf{r}-\mathbf{r}'|} A ( r ) = 4 π μ 0 I ∮ ∣ r − r ′ ∣ d l ′
Example: Vector potential of an infinite wire An infinite straight wire along the z z z -axis carries current I I I . In cylindrical coordinates ( s , ϕ , z ) (s, \phi, z) ( s , ϕ , z ) The vector potential can only depend on s s s by symmetry, and must point along z ^ \hat{\mathbf{z}} z ^ .
A ( s ) = μ 0 I 4 π ∫ − ∞ ∞ d z ′ s 2 + z ′ 2 z ^ \mathbf{A}(s) = \frac{\mu_0 I}{4\pi}\int_{-\infty}^{\infty}\frac{dz'}{\sqrt{s^2 + z'^2}}\,\hat{\mathbf{z}} A ( s ) = 4 π μ 0 I ∫ − ∞ ∞ s 2 + z ′2 d z ′ z ^
This integral diverges logarithmically. Introduce a cutoff at z ′ = ± L z' = \pm L z ′ = ± L :
A ( s ) ≈ μ 0 I 2 π ln ( 2 L s ) z ^ + c o n s t \mathbf{A}(s) \approx \frac{\mu_0 I}{2\pi}\ln\!\left(\frac{2L}{s}\right)\hat{\mathbf{z}} + \mathrm{const} A ( s ) ≈ 2 π μ 0 I ln ( s 2 L ) z ^ + const
Since A \mathbf{A} A is defined only up to a gauge transformation, we write:
A ( s ) = − μ 0 I 2 π ln ( s s 0 ) z ^ \mathbf{A}(s) = -\frac{\mu_0 I}{2\pi}\ln\!\left(\frac{s}{s_0}\right)\hat{\mathbf{z}} A ( s ) = − 2 π μ 0 I ln ( s 0 s ) z ^
Verify: B = ∇ × A = − ∂ A z ∂ s ϕ ^ = μ 0 I 2 π s ϕ ^ \mathbf{B} = \nabla \times \mathbf{A} = -\frac{\partial A_z}{\partial s}\,\hat{\boldsymbol{\phi}} = \frac{\mu_0 I}{2\pi s}\,\hat{\boldsymbol{\phi}} B = ∇ × A = − ∂ s ∂ A z ϕ ^ = 2 π s μ 0 I ϕ ^ . This matches the Ampere’s law result. ■ \blacksquare ■
Magnetization. The magnetization M \mathbf{M} M is the magnetic dipole moment per unit volume. It produces bound currents :
J b = ∇ × M , K b = M × n ^ \mathbf{J}_b = \nabla \times \mathbf{M}, \quad \mathbf{K}_b = \mathbf{M} \times \hat{\mathbf{n}} J b = ∇ × M , K b = M × n ^
The H field (magnetic field intensity) is defined as:
H = 1 μ 0 B − M \mathbf{H} = \frac{1}{\mu_0}\mathbf{B} - \mathbf{M} H = μ 0 1 B − M
Ampere’s law for H \mathbf{H} H :
∇ × H = J f \nabla \times \mathbf{H} = \mathbf{J}_f ∇ × H = J f
∮ H ⋅ d l = I f , e n c \oint \mathbf{H} \cdot d\mathbf{l} = I_{f,\mathrm{enc}} ∮ H ⋅ d l = I f , enc
This is simpler than Ampere’s law for B \mathbf{B} B because only free currents appear.
Linear magnetic materials. For isotropic linear materials:
M = χ m H , B = μ H \mathbf{M} = \chi_m \mathbf{H}, \quad \mathbf{B} = \mu \mathbf{H} M = χ m H , B = μ H
Where χ m \chi_m χ m is the magnetic susceptibility and μ = μ 0 ( 1 + χ m ) \mu = \mu_0(1 + \chi_m) μ = μ 0 ( 1 + χ m ) is the permeability. The relative permeability is μ r = 1 + χ m \mu_r = 1 + \chi_m μ r = 1 + χ m .
Diamagnetic materials (χ m < 0 \chi_m \lt 0 χ m < 0 , ∣ χ m ∣ ≪ 1 \lvert\chi_m\rvert \ll 1 ∣ χ m ∣ ≪ 1 ): Weakly repelled by Magnetic fields. The induced magnetization opposes the applied field (Lenz’s law at the Atomic level). Examples: bismuth, copper, water.
Paramagnetic materials (χ m > 0 \chi_m \gt 0 χ m > 0 , χ m ≪ 1 \chi_m \ll 1 χ m ≪ 1 ): Weakly attracted by magnetic fields. Atomic dipoles align partially with the applied field. Examples: aluminium, platinum, oxygen.
Ferromagnetic materials (χ m ≫ 1 \chi_m \gg 1 χ m ≫ 1 ): Strongly attracted by magnetic fields. Exhibit hysteresis : the magnetization depends on the history of the applied field.
The hysteresis loop traces B \mathbf{B} B vs H \mathbf{H} H as the external field cycles. Key Features:
Remanence B r B_r B r : the residual field when H = 0 H = 0 H = 0 .Coercivity H c H_c H c : the field required to demagnetize the material.Saturation : the maximum magnetization achievable.For soft ferromagnets (iron, nickel), H c H_c H c is small and the hysteresis loop is narrow. For hard ferromagnets (permanent magnets), H c H_c H c is large.
:::caution Common Pitfall The magnetic field B \mathbf{B} B is the fundamental quantity; H \mathbf{H} H is an auxiliary field Convenient for problems with free currents. The names “magnetic field” and “magnetic field Intensity” vary across textbooks --- always check which symbol a given text associates with Which name. In this document, B \mathbf{B} B is the magnetic field and H \mathbf{H} H is the Auxiliary H field.
:::