Electromagnetism

Table of Contents

1. Charge

A fundamental property matter can have.

1.1. Surface Charge Density

Charge per unit area

1.2. Volume Charge Density

Charge per unit volume

2. Coulomb's Law

The force \(\mathbf{F}\) between two point charge \(q_1\) and \(q_2\) is given by: \[ \mathbf{F} = \frac{1}{4\pi\varepsilon_0} \frac{q_1q_2}{r^2} \vu{r}, \] where \(\varepsilon_0\) is the permittivity of free space, and \(\mathbf{r}\) is the separation vector between two charge.

2.1. Electric Field

  • \(\mathbf{E}\), 전기장, 전장, 전계

Electric field is the mediating field that is generated by a stationary charge: \[ \mathbf{E} = \frac{1}{4\pi\varepsilon_0} \frac{q_1}{r^2}\vu{r}, \] and exerts Coulomb force on another charge: \[ \mathbf{F} = q_2\mathbf{E}. \]

2.2. Electric Potential

  • \(\phi\), \(V\)

Irrotational electric field in electrostatics admits a scalar potential function \(\phi\) such that: \[ \mathbf{E} = -\grad \phi. \]

When there is changing magnetic field, Faraday's law can be applied

\begin{align*} &\curl{\mathbf{E}} = - \pdv{t} \curl{\mathbf{A}} \\ \implies & \curl(\mathbf{E} + \pdv{\mathbf{A}}{t}) = \vb{0}, \end{align*}

and the definition of \(\phi\) is modified, and becomes dependent on the gauge: \[ \mathbf{E} + \pdv{\mathbf{A}}{t} = - \grad{\phi}. \] Rearranging, we get: \[ \mathbf{E} = - \grad{\phi} - \pdv{\mathbf{A}}{t}. \]

3. Current

  • \(I\), \(i\)

Amount of charge passing through a given surface per unit time.

3.1. Surface Current Density

  • \(\mathbf{K}\)

Amount of charge passing through, per unit time, a unit length perpendicular to the direction of the current.

\[ \mathbf{K} = \sigma \mathbf{v}. \]

3.2. Volume Current Density

  • \(\mathbf{J}\)
  • Current Density

Amount of charge passing through, per unit time, a unit surface area normal to the direction of the current.

\[ \mathbf{J} = \rho\mathbf{v}. \]

3.3. Conductivity

  • \(\sigma\)

In most continuous media, a charged particle gains drift velocity \(\mathbf{v}_{\rm d}\) under electric field. The proportional constant is called the mobility \(\mu\): \[ \mathbf{v}_{d} = \mu \mathbf{E}. \]

The drift velocity directly relates to current by charge density: \[ \mathbf{J} = \rho \mu \mathbf{E} = \sigma \mathbf{E}. \] The combined constant \(\sigma\) is called the conductivity.

In general, the proportionality constant between force per unit charge \(\mathbf{f}\) and current density \(\mathbf{J}\) is given by the conductivity: \[ \mathbf{J} = \sigma \mathbf{f}. \]

3.4. Resistivity

Resistivity \(\rho\) is defined as the inverse of conductivity \[ \rho = \frac{1}{\sigma}. \]

3.5. 4-Current

  • \(J^{\mu}\), \(j^{\mu}\)

Form of Number Flux. \[ J^{\mu} := \rho u^{\mu} \] where \(\rho\) is the rest charge density, and \(u^{\mu}\) is the 4-velocity of charge.

4. Ohm's Law

\[ V = IR \] where \(V\) is the electric potential difference across a systrem, or voltage, \(I\) is the current flowing through a system, \(R\) is the resistance that combines the resistivity of each material within the system into a single number.

5. Lorentz Force

\[ \mathbf{F} = q(\mathbf{E} + \mathbf{v} \times \mathbf{B}). \]

5.1. Magnetic Field

  • \(\mathbf{B}\), Magnetic Flux Density, Magnetic Induction

Magnetic field is the mediating field that is generated by a moving charge, and exerts force on another charge, according to Lorentz force.

5.1.1. Biot-Savart Law

The magnetic field generated by steady infinitesimal current element. By steady current: \[ \pdv{\rho}{t} = 0, \quad \pdv{\mathbf{J}}{t} = \mathbf{0} \] is meant.

\[ \mathbf{B}(\mathbf{r}) = \frac{\mu_0}{4\pi} \frac{\mathbf{J}(\mathbf{r}')\times (\widehat{\mathbf{r} - \mathbf{r}'})}{|\mathbf{r} - \mathbf{r}'|^2} \dd[3]{\mathbf{r}'}. \]

5.2. Magnetic Vector Potential

  • \(\mathbf{A}\)
  • \(V = -q\mathbf{v}\cdot\mathbf{A}\)

Solenoidal magnetic field admits a vector potential function \(A\) such that: \[ \mathbf{B} = \curl{\mathbf{A}}. \]

5.2.1. Gauge Invariance.

There is scalar field degrees of freedom when defining \(\mathbf{A}\). Particular choice of \(\mathbf{A}\) is called the gauge, and the law of physics should be invariant under different gauge.

5.2.2. Coulomb Gauge

  • Transverse Gauge, Radiation Gauge

\[ \div{\mathbf{A}} = 0. \] Such gauge exists in all physical circumstances, even in general relativity.

This gauge allows easy calculation of vector potential from current. \[ \laplacian \mathbf{A} = -\mu_0 \mathbf{J}, \] and therefore \[ \mathbf{A}(\mathbf{r}) = \frac{\mu_0}{4\pi} \frac{\mathbf{J}(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|} \dd[3]{\mathbf{r}'}. \]

5.3. Electromagnetic 4-Potential

The contravariant electromagnetic 4-potential is defined as: \[ A^{\mu} := \begin{bmatrix} \frac{1}{c}\phi \\ A_x \\ A_y \\ A_z \end{bmatrix} = \begin{bmatrix} \phi^\mathsf{G} \\ A_x^\mathsf{G} \\ A_y^\mathsf{G} \\ A_z^\mathsf{G} \end{bmatrix} \] where \(\phi\) is the electric potential, and \(A_i\) are the component of magnetic vector potential, and the superscript \(\mathsf{G}\) represent that the quantity is expressed in Gaussian units.

5.3.1. Gauge Tranformation

Arbitrary scalar field \(\lambda\) can be used for transformation:

\begin{align*} \mathbf{A}' &= \mathbf{A} + \grad{\lambda}, \\ \phi' &= \phi - \pdv{\lambda}{t}. \end{align*}

5.3.2. Lorenz Gauge

\[ \div{\mathbf{A}} + \frac{1}{c^2} \pdv{\phi}{t} = 0. \]

Using 4-gradient in spacetime coordinates: \[ \partial_{\mu}A^{\mu} = 0. \]

In this gauge, \[ \partial_{\mu}\partial^{\mu}A^{\nu} = \mu_0 J^{\nu} \] where \(J^{\nu}\) is the 4-current.

\(\partial_{\mu}\partial^{\mu}\) is called the d'Alembertian operator, and can be written as \(\partial^2\), \(\square\), \(\square^2\).

5.4. Potential Formulation

The potentials can be used instead of the fields: \[ \dv{\mathbf{p}}{t} = - q \left[ \pdv{\mathbf{A}}{t} + (\mathbf{v}\vdot \grad)\mathbf{A} + \grad(\phi - \mathbf{v}\vdot \mathbf{A}) \right]. \]

The first two terms are convective derivative, or total derivative, of \(\mathbf{A}\), and the equation rearranges to: \[ \dv{t} (\mathbf{p} + q\mathbf{A}) = - \grad[q(\phi - \mathbf{v}\vdot \mathbf{A})]. \] This motivates the potential energy term in Lagrangian, and its canonical momentum \((\partial_{\mathbf{v}}L)^{\mathsf{T}} := \mathbf{p} + q\mathbf{A}\): \[ L_{\rm em} = \frac{1}{2}mv^2 - q(\phi - \mathbf{v}\vdot \mathbf{A}). \]

6. Multipole Expansion

Charge distribution of arbitrary shape can be approximated by successive multipole moment:

\begin{align*} \phi(\mathbf{r}) &= \frac{1}{4\pi\varepsilon_0} \int \frac{\rho(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|} \dd[3]{\mathbf{r}'} \\ &= \frac{1}{4\pi\varepsilon_0} \sum_{n=0}^{\infty} \frac{1}{r^{n+1}} \int (r')^n P_n(\cos\alpha) \rho(\mathbf{r}')\dd[3]{\mathbf{r}'}, \end{align*}

where \(P_n\) is the Legendre polynomials, \(\alpha\) is the angle between \(\mathbf{r}'\) and \(\mathbf{r}\). The integral in each term corresponds to monopole, dipole, quadrupole, octopole, and so on.

Explicitly,

\begin{align*} \phi(\mathbf{r}) &= \frac{1}{4\pi\varepsilon_0} \Big[ \frac{1}{r}\int \rho(\mathbf{r}')\dd[3]{\mathbf{r}'} + \frac{1}{r^2} \int r' \cos\alpha \rho(\mathbf{r}')\dd[3]{\mathbf{r}} \\ &+ \frac{1}{r^3} \int (r')^2 \left( \frac{3}{2}\cos^2\alpha - \frac{1}{2}\right) \rho(\mathbf{r}') \dd[3]{\mathbf{r}} + \dots. \Big] \end{align*}

The expansion is derived from the Taylor expansion of \(1/|\mathbf{r} - \mathbf{r}'|\):

\begin{align*} |\mathbf{r} - \mathbf{r}'|^2 &= r^2 + (r')^2 - 2rr' \cos\alpha\\ \implies \frac{1}{|\mathbf{r} - \mathbf{r}'|} &= \frac{1}{r} \left( 1 + \frac{r'^2}{r^2} - \frac{2r'\cos\alpha }{r} \right)^{-1/2} \\ \end{align*}

6.1. Dipole Moment

6.1.1. Electric Dipole Moment

  • \(\mathbf{p}\)

The electric dipole memont when charge \(q\) is seperated with the displacement \(\mathbf{d}\) is: \[ \mathbf{p} := q\mathbf{d}. \]

In general, \[ \mathbf{p} := \int \mathbf{r} \rho \dd[3]{\mathbf{r}}. \]

The torque on the dipole \(\mathbf{p}\) under electric field is \[ \vb*{\tau} = \mathbf{p}\times \mathbf{E}. \] Notice the similarity to the torque on the magnetic dipole made of the circular circuit \[ \vb*{\tau} = I\mathbf{A}\times \mathbf{B} = \mathbf{m}\times \mathbf{B} \]

The force on the dipole due to electric field is \[ \mathbf{F} = \grad(\mathbf{p}\vdot\mathbf{E}) \]

The electric potential energy of the dipole \[ U = -\mathbf{p}\cdot \mathbf{E} \]

The electric potential caused by the dipole \[ \phi = \frac{1}{4\pi\varepsilon_0}\frac{\mathbf{p}\cdot \mathbf{\hat{r}}}{r^2} \]

6.1.2. Magnetic Dipole Moment

Magnetic dipole is a small current loop, which is called the Ampére dipole.

In general, \[ \mathbf{m} = \frac{1}{2} \int \mathbf{r}\cross \mathbf{J}\dd[3]{\mathbf{r}}. \] If the current is solenoidal, \[ \int r^i J^j \dd[3]{\mathbf{r}} = \varepsilon^{ij}{}_k m^k \] where identity metric is being used.

The torque on the dipole \(\mathbf{m}\) under magnetic field is \[ \vb*{\tau} = \mathbf{m} \cross \mathbf{B} \]

The force on the dipole due to magnetic field is \[ \mathbf{F} = \grad (\mathbf{m}\vdot\mathbf{B}) \]

The magnetic vector potential caused by the dipole is \[ \mathbf{A} = \frac{\mu_0}{4\pi}\frac{\mathbf{m}\cross \mathbf{\hat{r}}}{r^2}. \]

Under Lorentz boost by \(\mathbf{v}\),

6.1.2.1. Gilbert Dipole

Consists of two magnetic monopole.

6.1.3. Static Dipole Transformation

When the dipole moment is time-independent, the dipole moments transforms regularly under Lorentz transformation.

Under Lorentz boost by \(\mathbf{v}\),

\begin{align*} \tilde{\mathbf{p}}^{\parallel} &= \frac{\mathbf{p}^\parallel}{\gamma}, \\ \tilde{\mathbf{p}}^{\perp} &= \mathbf{p}^{\perp} - \frac{\mathbf{v} \cross \mathbf{m}}{c^2}, \\ \tilde{\mathbf{m}}^{\parallel} &= \frac{\mathbf{m}^\parallel}{\gamma}, \\ \tilde{\mathbf{m}}^{\perp} &= \frac{1 + \gamma^{-2}}{2}\mathbf{m}^{\perp} + \frac{1}{2}\mathbf{v} \cross \mathbf{p}. \end{align*}

This is reminiscent of the transformation of electromagnetic field.

7. Electromagnetism in Matter

7.1. Media

7.1.1. Conductor

Ideal conductor is defined to have

  • unlimited supply of free charges
  • electrons that are free to roam

Following properties can be derived:

  • \(\mathbf{E} = 0\) inside a conductor
  • \(\rho = 0\) inside a conductor
  • Net charge resides on the surface
  • Equipotential at all points
  • \(\mathbf{E}\) at the surface is normal to the surface.

7.1.2. Dielectric

  • Dielectric Medium

Linear dielectric is a ideal medium that has constant electric susceptibility for any range of electric field strength.

7.1.3. Paramagnet

Under external magnetic field, it is magnetized in the same direction.

7.1.4. Diamagnet

It is magnetized in the opposite direction.

7.1.5. Ferromagnet

It is magnetized in the same direction particularly strongly.

7.2. Polarization Density

  • \(\mathbf{P}\), Electric Polarization, Dielectric Polarization
  • Electric dipole moment density

7.2.1. Definition

Volume density of the electric dipole moment \[ \mathbf{P} = \operatorname*{equil\,lim}_{V\to 0} \frac{\mathbf{p}}{V} \] in locally equilibrated sense.

7.2.2. Properties

The electric potential due to the polarization field \(\mathbf{P}\) \[ V(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0}\int_\Omega \frac{\mathbf{P}(\mathbf{r}')\cdot \widehat{\mathbf{r}-\mathbf{r}'}}{|\mathbf{r}-\mathbf{r}'|^2}\,\mathrm{d}^{\wedge 3}\mathbf{r}' = \frac{1}{4\pi\varepsilon_0}\int_\Omega \mathbf{P}(\mathbf{r}')\cdot \nabla\big|^{\mathbf{r}'}\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right)\,\mathrm{d}^{\wedge 3}\mathbf{r}' \]

\begin{align*} V(\mathbf{r})\ \ =\quad\ & \frac{1}{4\pi\varepsilon_0}\left[\int_\Omega\nabla\big|^{\mathbf{r}'}\cdot\left(\frac{\mathbf{P}(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\right)\,\mathrm{d}^{\wedge 3}\mathbf{r}' -\int_\Omega \frac{1}{|\mathbf{r}-\mathbf{r}'|}\nabla\big|^{\mathbf{r}'}\mathbf{P}(\mathbf{r}')\,\mathrm{d}^{\wedge 3}\mathbf{r}' \right] \\[1em] \stackrel{\text{div. thm}}{=} \ & \frac{1}{4\pi\varepsilon_0}\left[\oint_{\partial\Omega}\frac{\mathbf{P}(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\cdot\mathrm{d}^{\wedge 2}\mathbf{r}' -\int_\Omega \frac{1}{|\mathbf{r}-\mathbf{r}'|}\nabla\big|^{\mathbf{r}'}\mathbf{P}(\mathbf{r}')\,\mathrm{d}^{\wedge 3}\mathbf{r}' \right] \\[1em] \ \ =\quad\ & \frac{1}{4\pi\varepsilon_0}\left[\oint_{\partial\Omega}\frac{\mathbf{P}(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\,\cdot{\mathbf{\hat{n}}}\,|\mathrm{d}^{\wedge 2}\mathbf{r}'| -\int_\Omega \frac{1}{|\mathbf{r}-\mathbf{r}'|}\nabla\big|^{\mathbf{r}'}\mathbf{P}(\mathbf{r}')\,\mathrm{d}^{\wedge 3}\mathbf{r}' \right] \\ \end{align*}

Using the bound charges and additionally letting \(\mathscr{r}'' := |\mathbf{r} - \mathbf{r}'|\), \[ V(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0}\oint_{\partial \Omega} \frac{\sigma_b}{r''}\,\mathrm{d}^{\wedge 2}\mathbf{r}' +\frac{1}{4\pi\varepsilon_0}\int_\Omega \frac{\rho_b}{r''}\,\mathrm{d}^{\wedge 3}\mathbf{r}' \]

7.3. Bound Charge

The bound surface charge density \( \sigma_b \) and the bound charge density \( \rho_b \) are defined as \[ \sigma_b := \mathbf{P}\cdot \mathbf{\hat{n}},\quad \rho_b := -\nabla\cdot \mathbf{P} \]

7.4. Displacement Field

  • \(\mathbf{D}\)

\[ \mathbf{D} = \varepsilon_0\mathbf{E} + \mathbf{P} \]

It equal to \(\mathbf{E}\) field outside of medium.

7.5. Molecular Polarizability

The induced dipole moment of individual molecule due to local electric field \(\mathbf{E}_{\rm local}\) is given by molecular polarizability: \[ \mathbf{p} = \varepsilon_0 \alpha \mathbf{E}_{\rm local}. \]

7.6. Electric Susceptibility

  • \(\chi_{\rm e}\)

Linear dielectric has linear relationship between imposed electric field and polarization density: \[ \mathbf{P} = \varepsilon_0 \chi_{\rm e} \mathbf{E}. \]

The permittivity of material can be defined \[ \varepsilon := \varepsilon_0 (1+\chi_{\rm e}) \implies \mathbf{D} = \varepsilon \mathbf{E} \] and the relative permittivity, or dielectric constant, is used to express them \[ \varepsilon_{\rm r} := \frac{\varepsilon}{\varepsilon_0} = 1 + \chi_{\rm e}. \]

7.6.1. Renormalization

Electric susceptibility can be renormalized from molecular polarizability using various models that describe particular cases.

Clausius-Mossotti Relation For weakly polarized, low density material, \[ \frac{\chi_{\rm e}{\chi_{\rm e} + 3} = \frac{n\alpha}{3\varepsilon_0} \] where \(n\) is the number density of the molecules.

7.7. Magnetization

  • \(\mathbf{M}\)
  • Magnetic dipole moment density

\[ \mathbf{M} = \operatorname*{equil\,lim}_{V\to 0}\frac{\mathbf{m}}{V} \] where \(\mathbf{m}\) is the magnetic moment.

7.8. Bound Current

The bound surface current density \(\mathbf{K}_b\) and the bound volume current density \(\mathbf{J}_b\) are defined as: \[ \mathbf{K}_b = \mathbf{M}\times \vu{n}, \quad\mathbf{J}_b = \nabla\times \mathbf{M} + \frac{\partial \mathbf{P}}{\partial t}. \]

The bound volume current density is sometimes devided into bound current \(\curl{\mathbf{M}}\) and polarization current \(\partial_t \mathbf{P}\).

7.8.1. Properties

\[ \div{\mathbf{J}_b} = 0 \]

7.9. Magnetizing Field

  • \(\mathbf{H}\), Magnetic Field Strength, H-Field, Magnetic Field, Auxiliary Field

\(\mathbf{H} = \mathbf{B}/\mu_0 - \mathbf{M}\)

It is the \(\mathbf{B}\) field outside, but without the solenoid field inside.

7.10. Magnetic Susceptibility

  • \(\chi_{\rm m}\), volume magnetic susceptibility

The relation between magnetization and magnetizing field is linear for most substances under moderate magnetic field: \[ \mathbf{M} = \chi_{\rm m} \mathbf{H} \]

It is non-linear in ferromagnets and and tensor-valued in general. Differential susceptibility is used for such cases: \[ \chi_{\mathrm{d}}{}^i{}_j = \pdv{M^i}{H^j}. \]

8. Boundary Conditions

8.1. Electriostatic

Only the normal component of electric field is discontinuous at a surface in the presense of surface charge:

\begin{align*} &\Delta E^{\perp} = \frac{\sigma}{\varepsilon_0} \\ &\Delta E^{\parallel} = 0. \end{align*}

They simplifies to: \[ \Delta \mathbf{E} = \frac{\sigma}{\varepsilon_0} \vu{n}. \]

The electric potential should be continuous across boundary, but its derivative can be discontinuous: \[ \Delta (\vu{n} \vdot \grad V) = - \frac{\sigma}{\varepsilon_0}. \]

For displacement field, the condition for electric field becomes:

\begin{align*} &\Delta D^{\perp} = \sigma_f \\ &\Delta \mathbf{D}^{\parallel} = \Delta \mathbf{P}^{\parallel} \end{align*}

8.2. Magnetostatic

Only tangential component of magnetic field is discontinuous at a surface in the presense of surface current:

\begin{align*} &\Delta B^{\perp} = 0\\ &\Delta B^{\parallel} = \mu_0 K \end{align*}

where the parallel component is in the direction given by the right hand rule. These conditions can be combined into \[ \Delta \mathbf{B} = \mu_0 (\mathbf{K} \times \vu{n}). \]

\(\mathbf{A}\) cannot be "streched" but can be "sheared" in Coulomb gauge. The vector potential must be continuous across boundary, but its derivative can be discontinuous: \[ \Delta (\vu{n} \vdot \grad \mathbf{A}) = - \mu_0 \mathbf{K}. \]

For magnetizing field, the condition for magnetic field becomes:

\begin{align*} &\Delta \mathbf{H}^{\perp} = - \Delta \mathbf{M}^{\perp} \\ &\Delta \mathbf{H}^{\parallel} = \mathbf{K}_f \cross \vu{n} \end{align*}

9. Electromotive Force

  • Electromotance, EMF

In a system, often a circuit, there exists two forces that drives the current, the external force per charge \(\mathbf{f}_s\), and electrostatic force generated within the system: \[ \mathbf{f} = \mathbf{f}_s + \mathbf{E}. \]

The electromotive force is the net effect of source around a desired path: \[ \mathcal{E} := \oint \mathbf{f} \vdot \dd{\mathbf{l}} = \oint \mathbf{f}_s \vdot \dd{\mathbf{l}}. \]

9.1. Wire Moving through Magnetic Field

moving-wire-export.svg

Two forces are applied to charge in vertical wire being pulled sideways. The electromotive force calculate the vertical component of magnetic force per unit charge \(\mathbf{f}_{\rm mag}\) at an instant: \[ \mathbf{f}_{\rm mag}\vdot \dd{\mathbf{l}} = vB\dd{\ell} \] where \(\dd{\mathbf{l}} = \mathbf{u} \dd{t}\). But the work is done by \(\mathbf{f}_{\rm pull}\): \[ \mathbf{f}_{\rm pull}\vdot \mathbf{w} \dd{t} = (uB) \frac{\dd{\ell}}{\cos\theta}\sin\theta = vB\dd{\ell} \]

9.2. Flux Rule

For a freely moving closed wire under magnetic field: \[ \mathcal{E} = - \dv{\Phi}{t}. \]

10. Faraday's Law

Due to relativity, moving magnetic field relative to the wire should also induce electromotive force. Electric field is generated by changing magnetic field to move stationary charge.

Even for changing magnetic field, the flux rule is true: \[ \mathcal{E} = \oint \mathbf{E} \vdot \dd{\mathbf{l}} = - \dv{\Phi}{t}. \] Lenz's law says that the electric field would be generated in a way that opposes the change in flux.

When \(\rho = 0\) thus \(\div{\mathbf{E}} = 0\), \[ \curl{\mathbf{E}} = - \pdv{\mathbf{B}}{t} \] suggests that \[ \mathbf{E}(\mathbf{r}) = - \frac{1}{4\pi} \int \frac{(\partial_t \mathbf{B}) \cross (\widehat{\mathbf{r} - \mathbf{r}'})}{|\mathbf{r} - \mathbf{r}'|^2} \dd[3]{\mathbf{r}}'. \]

11. Inductance

Inductance is proportionality constant \(M\) between the current and magnetic flux: \[ \Phi = MI. \] The same constant relates the electromotive force and change in current: \[ \mathcal{E} = - \dv{\Phi}{t} = - M \dv{I}{t}. \]

11.1. Mutual Inductance

The proportionality constant between the current in one loop and the magnetic flux through another loop: \[ \Phi_2 = M_{21}I_1. \]

11.1.1. Neumann Formula

Since \[ \Phi_2 = \oint \mathbf{A}_1 \vdot \dd{\mathbf{l}_2}, \] substituting in the formula of vector potential in terms of current gives: \[ M_{21} = \frac{\mu_0}{4\pi} \oint \oint \frac{\dd{\mathbf{l}_1}\dd{\mathbf{l}_2}}{|\mathbf{r}_2 - \mathbf{r}_1|}. \]

Mutual inductance is symmetric because the integral is the same regardless of the order of integration.

11.2. Self Inductance

The inductance on the loop itself.

12. Maxwell's Equations

\begin{align*} &\div{\mathbf{E}} = \frac{\rho}{\varepsilon_0} \\ &\div{\mathbf{B}} = 0 \\ &\curl{\mathbf{E}} = - \pdv{\mathbf{B}}{t} \\ &\curl{\mathbf{B}} = \mu_0 \mathbf{J} + \mu_0\varepsilon_0\pdv{\mathbf{E}}{t} \end{align*}

Maxwell called the extra term the displacement current: \[ \mathbf{J}_{\rm d} := \varepsilon_0 \pdv{\mathbf{E}}{t}. \]

12.1. Macroscopic Formulation

\begin{align*} &\div{\mathbf{D}} = \rho_f \\ &\div{\mathbf{B}} = 0 \\ &\curl{\mathbf{E}} = - \pdv{\mathbf{B}}{t} \\ &\curl{\mathbf{H}} = \mathbf{J}_f + \pdv{\mathbf{D}}{t} \end{align*}

The first equation:

\begin{align*} \varepsilon_0\div{\mathbf{E}} &= \rho \\ &= \rho_f + \rho_b \\ &= \rho_f - \div{\mathbf{P}} \\ \implies &\div{\mathbf{D}} = \rho_f \end{align*}

The fourth equation:

\begin{align*} \frac{1}{\mu_0} \curl{\mathbf{B}} &= \mathbf{J} + \varepsilon_0 \pdv{\mathbf{E}}{t} \\ &= \mathbf{J}_f + \mathbf{J}_b \\ &= \mathbf{J}_f + \curl{\mathbf{M}} + \pdv{\mathbf{P}}{t} \end{align*}

12.2. Tensor Formulation

Using Electromagnetic Tensor,

\begin{align*} &\partial_{\mu}G^{\mu\nu} = 0, \\ &\partial_{\mu}F^{\mu\nu} = \mu_0 J^{\nu} \end{align*}

The first equation is also given by: \[ \partial_{[\alpha} F_{\beta\gamma]} = \partial_{\alpha} F_{\beta\gamma} + \partial_{\beta}F_{\gamma\alpha} + \partial_{\gamma}F_{\alpha\beta} = 0. \]

12.3. Potential Formulation

\begin{align*} &\laplacian \phi + \pdv{t} (\div{\mathbf{A}}) = - \frac{\rho}{\varepsilon_0} \\ & \left( \laplacian \mathbf{A} - \frac{1}{c^2} \pdv[2]{\mathbf{A}}{t} \right) - \grad(\div{\mathbf{A}} + \frac{1}{c^2} \pdv{\phi}{t}) = - \mu_0\mathbf{J}. \end{align*}

Using tensors, \[ \partial_{\mu}\partial^{\mu}A^{\nu} - \partial_{\mu}\partial^{\nu}A^{\mu} = \mu_0J^{\nu}. \]

In Lorenz gauge, \(\partial_{\mu}A^{\mu} = 0\), simplifying it to \[ \partial_{\mu}\partial^{\mu}A^{\nu} = \mu_0 J^{\nu}, \] which expands as:

\begin{align*} \square^2 \phi &= -\frac{\rho}{\varepsilon_0},\\ \square^2 \mathbf{A} &= -\mu_0\mathbf{J}. \end{align*}

This is in the form of inhomogeneous wave equation, and \(\rho, \mathbf{J}\) become the source terms.

13. Continuity Equations

13.1. Charge Conservation

\[ \pdv{\rho}{t} = - \div{\mathbf{J}}. \]

Using 4-current, \[ \partial_{\mu} J^{\mu} = 0. \]

13.2. Energy Conservation

  • Poynting's Theorem

Work \(W\) done by electromagnetic field is given by:

\begin{align*} \dv{W}{t} &= \int_{\Omega} \mathbf{E} \vdot \mathbf{J} \dd[3]{\mathbf{r}} \\ &= - \dv{t} \int_{\Omega} \frac{1}{2} \left( \varepsilon_0 E^2 + \frac{1}{\mu_0} B^2\right) \dd[3]{\mathbf{r}} - \frac{1}{\mu_0} \oint_{\partial \Omega} (\mathbf{E}\cross \mathbf{B}) \vdot \dd[2]{\mathbf{r}} \end{align*}

The energy density \(u\) in electromagnetic field: \[ u := \frac{1}{2} \left( \varepsilon_0 E^2 + \frac{1}{\mu_0}B^2\right), \] and the Poynting vector \(\mathbf{S}\) that represent the energy flux (density) can be defined: \[ \mathbf{S} := \frac{1}{\mu_0} \mathbf{E}\cross \mathbf{B}. \]

When no work is done on the charges: \[ \pdv{u}{t} = - \div{\mathbf{S}} \] where \(u\) is the energy density in the first term.

13.3. Momentum Conservation

Force exerted by electromagnetic field is given by:

\begin{align*} \mathbf{F} &= \int \left(\rho \mathbf{E} + \mathbf{J}\cross \mathbf{B} \right) \dd[3]{\mathbf{r}}\\ &= \oint_{\partial\Omega} \left[ \varepsilon_0 \left( E^iE^j - \frac{1}{2}\delta^{ij} E^2\right) + \frac{1}{\mu_0} \left( B^iB^j - \frac{1}{2} \delta^{ij} B^2 \right)\right](\dd[2]{\mathbf{r}})_i \mathbf{e}_j - \varepsilon_0\mu_0 \dv{t} \int_{\Omega} \mathbf{S} \dd[3]{\mathbf{r}}. \end{align*}

The momentum density \(\mathbf{g}\): \[ \mathbf{g} := \mu_0\varepsilon_0 \mathbf{S} = \varepsilon_0 \mathbf{E}\cross \mathbf{B} \] and the electromagnetic stress or negative momentum flux (density) called Maxwell stress tensor \(\mathbf{T} = T^{ij}\mathbf{e}_i\otimes\mathbf{e}_j\) can be defined: \[ T^{ij} = \varepsilon_0 \left( E^iE^j - \frac{1}{2}\delta^{ij} E^2\right) + \frac{1}{\mu_0} \left( B^iB^j - \frac{1}{2} \delta^{ij} B^2 \right). \]

Total momentum of localized static field can be given in multiple equivalent forms:

\begin{align*} \mathbf{p} &= \varepsilon_0 \int \mathbf{E}\cross \mathbf{B}\dd[3]{\mathbf{r}} \\ &= \varepsilon_0\mu_0 \int \phi \mathbf{J} \dd[3]{\mathbf{r}} \\ &= \int \rho \mathbf{A}_{\rm C} \dd[3]{\mathbf{r}} \end{align*}

where \(\mathbf{A}_{\rm C}\) is the vector potential in Coulomb gauge.

Naturally, electromagnetic field has angular momentum: \[ \vb*{\ell} = \mathbf{r}\cross \mathbf{g}. \]

14. Electromagnetic Tensor

  • Electromagnetic Field Tensor, Field Strength Tensor, Faraday Tensor, Maxwell Bivector

Electromagnetic tensor is the exterior derivative of the covariant electromagnetic 4-potential: \[ F := \dd{A}. \] That is, \[ F_{\mu\nu} := \partial_{\mu}A_{\nu} - \partial_{\nu}A_{\mu} = \partial_{\mu}\eta_{\nu\alpha}A^{\alpha} - \partial_{\nu}\eta_{\mu\beta}A^{\beta} \] where \(\eta_{\mu\nu}\) is the Minkowski Metric.

To be explicit by employing mostly-minus Minkowski metric with signature (+,-,-,-),

\begin{alignat*}{2} F_{\mu\nu} &= \begin{bmatrix} 0 & E_x/c & E_y/c & E_z/c \\ -E_x/c & 0 & -B_z & B_y \\ -E_y/c & B_z & 0 & -B_x \\ -E_z/c & -B_y & B_x & 0 \end{bmatrix}, \quad &F^{\mu\nu} &= \begin{bmatrix} 0 & -E_x/c & -E_y/c & -E_z/c \\ E_x/c & 0 & -B_z & B_y \\ E_y/c & B_z & 0 & -B_x \\ E_z/c & -B_y & B_x & 0 \end{bmatrix} \end{alignat*}

The transformation rule can be derived from capacitor and solenoid(David J. Griffiths, 1989)

\begin{align*} \tilde{\mathbf{E}} &= \gamma (\mathbf{E} + \mathbf{v} \times \mathbf{B}) \\ \tilde{\mathbf{B}} &= \gamma\left(\mathbf{B} - \frac{1}{c^2} (\mathbf{v}\cross \mathbf{E})\right). \end{align*}

Note that Griffiths defines the electromagnetic tensor as \(F = - \dd{A}\) to match his mostly-plus convention, so that the definition of Minkowski force stay identical.

It has its Hodge dual: \[ G^{\mu\nu} = \frac{1}{2} \varepsilon^{\mu\nu\rho\sigma} F_{\rho\sigma} = \begin{bmatrix} 0 & -B_x & -B_y & -B_z \\ B_x & 0 & E_z/c & -E_y / c \\ B_y & -E_z / c & 0 & E_x /c \\ B_z & E_y /c & -E_x /c & 0 \end{bmatrix}. \]

14.1. Minkowski Force

\[ K^{\mu} = q u_{\nu} F^{\mu\nu} \] where \(u_{\mu}\) is the 4-velocity of the point charge.

The frame force of it is given by Lorentz force.

15. Electromagnetic Radiation

Radiation is energy emitted by charge and propagates infinitely far. The radiative power is given by: \[ P_{\rm rad}(t_0) = \lim_{r \to \infty} \oint \mathbf{S} \vdot \dd[2]{\mathbf{r}}. \]

Radiation occurs when there exists terms in electric and magnetic fields that falls off as inverse distance. For static system, the fields at most falls off as inverse square of distance.

15.1. Retarded Potential

In a dynamic system, the potential changes as the charge distribution changes:

\begin{align*} \phi(\mathbf{r}, t) &= \frac{1}{4\pi \varepsilon_0} \int \frac{\rho(\mathbf{r}', t - |\mathbf{r} - \mathbf{r}'|/c)}{|\mathbf{r} - \mathbf{r}'|} \dd[3]{\mathbf{r}'}, \\ \mathbf{A}(\mathbf{r}, t) &= \frac{\mu_0}{4\pi} \int \frac{\mathbf{J}(\mathbf{r}', t - |\mathbf{r} - \mathbf{r}'|/c)}{|\mathbf{r} - \mathbf{r}'|} \dd[3]{\mathbf{r}'}, \end{align*}

where \(t_{r} = t - |\mathbf{r} - \mathbf{r}'|/c\) is called the retarded time. These potentials satisfy the Lorenz gauge and its inhomogeneous wave equation.

The same formula does not apply for electric or magnetic field.

15.2. Jefimenko's Equations

The electric field and magnetic field can be obtained from the potential. When the retarded potential is substituted in, Jefimenko's equations are obtained, which help us calculate the fields from the charge and current, possibly easier:

\begin{align*} \mathbf{E}(\mathbf{r}, t) &= \frac{1}{4\pi\varepsilon_0} \int \left[ \frac{\rho(\mathbf{r}', t_r)}{|\mathbf{r} - \mathbf{r}'|^2}(\widehat{\mathbf{r} - \mathbf{r}'}) + \frac{\dot{\rho}(\mathbf{r}', t_r)}{c|\mathbf{r} - \mathbf{r}'|}(\widehat{\mathbf{r} - \mathbf{r}'}) - \frac{\dot{\mathbf{J}}(\mathbf{r}', t_r)}{c^2 |\mathbf{r} - \mathbf{r}'|} \right] \dd[3]{\mathbf{r}'} \\ \mathbf{B}(\mathbf{r}, t) &= \frac{\mu_0}{4\pi} \int \left[ \frac{\mathbf{J}(\mathbf{r}', t_r)}{|\mathbf{r} - \mathbf{r}'|^2} + \frac{\dot{\mathbf{J}}(\mathbf{r}', t_r)}{c|\mathbf{r} - \mathbf{r}'|} \right] \times (\widehat{\mathbf{r} - \mathbf{r}'}) \dd[3]{\mathbf{r}'} \\ \end{align*}

15.3. Liénard-Wiechert Potentials

This is the potential for a moving point charge. The charge appears to be different by the factor \((1 - (\widehat{\mathbf{r}- \mathbf{r}'})\vdot \mathbf{v} / c)^{-1}\).

The apparant elongation can be calculated from the difference in time it takes to arrive from the head and the tail of the object: \[ \frac{L'\cos\theta}{c} = \frac{L' - L}{v} \implies L' = \frac{L}{1 - v\cos\theta / c} \] where \(\theta\) is the angle between the velocity and the observer.

Given the trajectory of point charge \(\mathbf{w}(t)\),

\begin{align*} \phi(\mathbf{r}, t) &= \frac{1}{4\pi\varepsilon_0} \frac{q}{|\mathbf{r} - \mathbf{w}(t_r)| (1 -(\widehat{\mathbf{r}- \mathbf{w}(t_r)})\vdot \mathbf{v}(t_r) / c )}, \\ \mathbf{A}(\mathbf{r}, t) &= \frac{\mu_0}{4\pi} \frac{q\mathbf{v}(t_r)}{|\mathbf{r} - \mathbf{w}(t_r)| (1 -(\widehat{\mathbf{r}- \mathbf{w}(t_r)})\vdot \mathbf{v}(t_r) / c )} = \frac{\mathbf{v}(t_r)}{c^2}\phi(\mathbf{r}, t) \\ \end{align*}

where \(\mathbf{v} = \dot{\mathbf{w}}\), and \(t_r\) is given by the implicit equation \[ | \mathbf{r} - \mathbf{w}(t_r) | = c(t- t_r) \] with only one solution for each \(\mathbf{r}\), because the worldline of the point charge can only cross the light cone once.

From this potential, the fields can be calculated:

\begin{align*} \mathbf{E}(\mathbf{r}, t) &= \frac{1}{4\pi\varepsilon_0} \frac{q |\mathbf{r} - \mathbf{w}(t_r)|}{((\mathbf{r} - \mathbf{w}(t_r)) \vdot \mathbf{u}(t_r))^3} \left[ (c^2 - v(t_r)^2)\mathbf{u}(t_r) + (\mathbf{r} - \mathbf{w}(t_r)) \cross (\mathbf{u}(t_r) \cross \ddot{\mathbf{w}}(t_r)) \right] \\ \mathbf{B}(\mathbf{r}, t) &= \frac{1}{c} (\widehat{\mathbf{r} - \mathbf{w}(t_r)}) \times \mathbf{E}(\mathbf{r}, t). \end{align*}

where \(\mathbf{u} := c(\widehat{\mathbf{r} - \mathbf{w}(t_r)}) - \mathbf{v}(t_r)\).

The first term decay as inverse square of distance and does not depend on the acceleration of the charge, and it is called generalized Coulomb field, or velocity field. The second term decay as inverse of distance and is called radiation field, or acceleration field.

The radiative power of point charge with zero velocity and acceleration \(a\) can be calculated from radiative field term, and given in Larmor formula: \[ \expval{P} = \frac{\mu_0 q^2 a^2}{6\pi c}. \] In general, with point charge with velocity, the radiative power is given by Liénard's generalization: \[ P = \frac{\mu_0q^2\gamma^6}{6\pi c} \left( a^2 - \left| \frac{\mathbf{v}\times \mathbf{a}}{c} \right|^2 \right). \] where \(\gamma\) is the gamma factor in relativity.

15.4. Electric Dipole Radiation

Given a electric dipole with charge concentrated at both ends and oscillates: \[ \mathbf{p}(t) = q_0 d\cos(\omega t) \vu{z} = p_0 \cos\theta \vu{z}, \] the retarded potential is given by

\begin{align*} \phi(\mathbf{r}, t) &= \frac{1}{4\pi\varepsilon_0} \left[ \frac{q_0\cos(\omega (t - |\mathbf{r} - \mathbf{r}_+|/c))}{|\mathbf{r} - \mathbf{r}_+|} - \frac{q_0\cos(\omega(t - |\mathbf{r} - \mathbf{r}_-|))}{|\mathbf{r} - \mathbf{r}_-|} \right]\\ \mathbf{A}(\mathbf{r}, t) &= \frac{\mu_0}{4\pi} \int_{-d/2}^{d/2} \frac{-q_0\omega \sin[\omega(t - |\mathbf{r} - \mathbf{r}'|/c)]}{|\mathbf{r} - \mathbf{r}'|} \vu{z} \dd{z}. \end{align*}

Using the approximation at the radiative zone \(d \ll \lambda \ll r\):

\begin{align*} \phi(r, \theta, t) &= - \frac{p_0\omega}{4\pi\varepsilon_0 c} \left( \frac{\cos\theta}{r} \right) \sin [\omega(t - r/c)], \\ \mathbf{A}(r, \theta, t) &= - \frac{\mu_0 p_0\omega}{4\pi r} \sin [\omega(t - r/c)] \vu{z},\\ \end{align*}

and the fields becomes:

\begin{align*} \mathbf{E} &= - \frac{\mu_0p_0\omega^2}{4\pi} \left( \frac{\sin\theta}{r} \right)\cos[\omega(t - r/c)]\vu{\theta}, \\ \mathbf{B} &= - \frac{\mu_0p_0\omega^2}{4\pi c} \left( \frac{\sin\theta}{r} \right)\cos[\omega(t - r/c)]\vu{\phi}. \end{align*}

The average radiative power of dipole radiation is then: \[ \expval{P} = \frac{\mu_0p_0^2\omega^4}{12\pi c}. \]

15.5. Magnetic Dipole Radiation

Given a magnetic dipole of small circular circuit with oscillating current: \[ \mathbf{m}(t) = \pi b^2 I_0 \cos\theta \vu{z} = m_0 \cos\theta \vu{z}, \] the same procedure with radiative zone approximation \(d \ll \lambda \ll r\), almost identical result is obtained:

\begin{align*} \mathbf{E} &= - \frac{\mu_0m_0\omega^2}{4\pi c} \left( \frac{\sin\theta}{r} \right)\cos[\omega(t - r/c)]\vu{\phi}, \\ \mathbf{B} &= - \frac{\mu_0m_0\omega^2}{4\pi c^2} \left( \frac{\sin\theta}{r} \right)\cos[\omega(t - r/c)]\vu{\theta} \end{align*}

with radiative power of \[ \expval{P} = \frac{\mu_0m_0^2\omega^4}{12\pi c^3}. \]

The radiative power is tiny compared to electric dipole \[ \frac{P_{\rm magnetic}}{P_{\rm electric}} = \left( \frac{m_0}{p_0 c} \right)^2 = \left( \frac{\omega b}{c} \right)^2. \]

15.6. Radiation Reaction Force

The charge experience recoil force by the emitted field, From the average energy emitted by charge:

\begin{align*} &\int \mathbf{F}_{\rm rad} \vdot \mathbf{v} \dd{t} = - \frac{\mu_0q^2}{6\pi c} \int a^2 \dd{t} \\ \implies&\int \left( \mathbf{F}_{\rm rad} - \frac{\mu_0q^2}{6\pi c}\dot{\mathbf{a}} \right) \vdot \mathbf{v} \dd{t}. \end{align*}

Abraham-Lorentz formula is inspired by this: \[ \mathbf{F}_{\rm rad} = \frac{\mu_0q^2}{6 \pi c} \dot{\mathbf{a}}. \] This formula has to be used with care, because one solution to this is exponential growth.

Reaction force arises from the first order time-derivative of acceleration term of self-force, which charge exerts on itself retardedly.

16. Bloch Equation

  • Introduced by Felix Bloch in 1946.

Macroscopic equations that describe the change in magnetization as a magnetic field is imposed.

\[ \frac{d\mathbf{M}}{dt}=\gamma\mathbf{M}\times\mathbf{B}-\mathbf{T}(\mathbf{M}-M_0\hat{\mathbf{z}}) \] where \(\mathbf{T}\) is the relaxation time matrix with \(T_{11}=T_{22}=T_2\) and \(T_{33}=T_1\).

17. Gyromagnetic Ratio

17.1. Classical

\[ \gamma := \frac{q}{2m} \]

17.2. Isolated Electron

\[ \gamma_{\rm e} := \frac{-e}{2m_{\rm e}}|g_{\rm e}| = \frac{g_{\rm e}\mu_{\rm B}}{\hbar} \] It deviate from the classical result by the g-factor.

17.3. Nucleus

\[ \gamma = \frac{e}{2m_{\rm p}}g_{\rm n} = \frac{g_{\rm n}\mu_{\rm N}}{\hbar} \]

18. Scattering

19. References

Author: Jeemin Kim

Created: 2026-09-29 Tue 00:21