The Dual Field Tensor

Two of Maxwell's equations have sources on the right. Two do not. Written as vectors, that looks like an accident of history — Gauss and Ampère got charges and currents, Faraday and no-monopoles got zero. Written in spacetime, it stops being an accident and becomes a statement. Contract the field tensor with the four-dimensional Levi-Civita symbol and you build a second antisymmetric tensor, the dual, Gμν=12 εμναβFαβG^{\mu\nu} = \tfrac12\,\varepsilon^{\mu\nu\alpha\beta}F_{\alpha\beta}. Component by component, GG is FF with E⃗→cB⃗\vec E \to c\vec B and B⃗→−E⃗/c\vec B \to -\vec E/c — an electric-magnetic rotation performed by index gymnastics. And the two sourceless Maxwell equations, ∇⋅B⃗=0\nabla\cdot\vec B = 0 and Faraday's law, collapse into one line:

∂μGμν=0.\partial_\mu G^{\mu\nu} = 0.

Now put the two equations side by side. ∂μFμν=μ0Jν\partial_\mu F^{\mu\nu} = \mu_0 J^\nu says electric charge sources the field. ∂μGμν=0\partial_\mu G^{\mu\nu} = 0 says there is nothing on the right — and the thing that would go there has a name. If magnetic charge existed the equation would read ∂μGμν=Jmν/c\partial_\mu G^{\mu\nu} = J_m^\nu/c, with JmνJ_m^\nu a magnetic four-current. The zero is not decoration. It is a measurement. Every time anyone has looked for a magnetic monopole — in cosmic rays, in moon rock, at the LHC — the answer has come back zero. That is why this article's portrait is the magnetic monopole: the dual tensor makes the monopole's absence into a physical law rather than an omission.

Key takeaways
  • Signature (+,−,−,−)(+,-,-,-), η=diag(+1,−1,−1,−1)\eta = \mathrm{diag}(+1,-1,-1,-1), p⋅p=+m2c2p\cdot p = +m^2c^2; ε0123=+1\varepsilon^{0123} = +1; mm is the invariant rest mass.
  • The dual: Gμν≡12 εμναβFαβ\boxed{G^{\mu\nu} \equiv \tfrac12\,\varepsilon^{\mu\nu\alpha\beta}F_{\alpha\beta}} — the 12\tfrac12 undoes the double-counting of the antisymmetric pair (αβ)(\alpha\beta). It implements E⃗→cB⃗\vec E \to c\vec B,   B⃗→−E⃗/c\;\vec B \to -\vec E/c: a 90°90° rotation in the (E⃗/c,B⃗)(\vec E/c, \vec B) plane. Twice is 180°180°, hence ∗∗Fμν=−Fμν\boxed{{}^{**}F^{\mu\nu} = -F^{\mu\nu}}.
  • The homogeneous pair, in one equation: ∂μGμν=0\boxed{\partial_\mu G^{\mu\nu} = 0}. ν=0\nu = 0 gives ∇⋅B⃗=0\nabla\cdot\vec B = 0; ν=i\nu = i gives ∇×E⃗=−∂tB⃗\nabla\times\vec E = -\partial_t\vec B. It is equivalent to the Bianchi identity ∂[λFμν]=0\partial_{[\lambda}F_{\mu\nu]} = 0 — automatic once F=dAF = \mathrm{d}A, because partial derivatives commute. Not dynamics: the price of having a potential.
  • Contracting FF with GG gives the pseudoscalar FμνGμν=−4c E⃗⋅B⃗\boxed{F_{\mu\nu}G^{\mu\nu} = -\tfrac{4}{c}\,\vec E\cdot\vec B} — the invariant that FF alone cannot see.
  • Owned result: the monopole source term, with its coefficient derived, not guessed. ∂μGμν=1cJmν\boxed{\partial_\mu G^{\mu\nu} = \tfrac1c J_m^\nu} with Jmν=(cρm,J⃗m)J_m^\nu = (c\rho_m, \vec J_m) reproduces ∇⋅B⃗=ρm\nabla\cdot\vec B = \rho_m and ∇×E⃗=−∂tB⃗−J⃗m\nabla\times\vec E = -\partial_t\vec B - \vec J_m. The AllQFT base prints μ0c Jmν\mu_0 c\,J_m^\nu; that coefficient is dimensionally impossible and is refused here.
  • Key idea: duality maps electric charge qq to magnetic charge g=μ0qcg = \mu_0 q c. Do that to an electron and you get 6.04×10−17 Wb6.04\times10^{-17}\ \mathrm{Wb} — but Dirac's condition says the smallest legal monopole carries h/e=4.14×10−15 Wbh/e = 4.14\times10^{-15}\ \mathrm{Wb}, exactly 1/2α=68.51/2\alpha = 68.5 times larger. Duality is a symmetry of the classical equations; quantum mechanics sets the price of admission.
The core insight

The Levi-Civita symbol is the only Lorentz-invariant object that knows about orientation, so contracting with it is the only way to build a new tensor out of FF without differentiating it. That single construction does three jobs at once: it rotates E⃗\vec E into B⃗\vec B, it packages the two source-free Maxwell equations into a divergence, and it exposes the pseudoscalar E⃗⋅B⃗\vec E\cdot\vec B, which no contraction of FF with itself can ever produce. The dual is not a notational convenience. It is the second half of the electromagnetic field.

Where this lives today

Duality is the entry point to a great deal of modern theory, and the dual tensor is where it starts. The invariant FμνGμνF_{\mu\nu}G^{\mu\nu} is a total derivative in vacuum, so adding it to the Lagrangian changes nothing classically — but multiply it by a field and it becomes the axion coupling, a leading candidate both for solving the strong-CP problem and for the dark matter of the universe. Axion haloscopes such as ADMX are, quite literally, experiments to detect a term built out of GμνG^{\mu\nu}. Duality also organises strong coupling: Montonen–Olive and Seiberg–Witten theory exploit the fact that a theory with electric charges at strong coupling can be re-described as one with magnetic charges at weak coupling. The factor 68.5=1/2α68.5 = 1/2\alpha computed above is exactly why: duality maps weak coupling to strong. And the null result keeps its own importance. Every bound on the monopole flux — MoEDAL, cosmic-ray detectors, the "Parker bound" from the survival of the galactic magnetic field — is a bound on the right-hand side of ∂μGμν=0\partial_\mu G^{\mu\nu} = 0. Be clear about the logic: the equation is not a theorem, it is a summary of the data. Were a monopole found tomorrow the Bianchi identity would not be violated — the potential would be, and electromagnetism would acquire a topological sector. A far more interesting failure mode than a broken equation.

The physical picture

FμνF^{\mu\nu} is an antisymmetric 4×44\times4 table with six independent entries: three time-space slots holding E⃗/c\vec E/c, three space-space slots holding B⃗\vec B. That split is not sacred — a boost reshuffles the entries — but the shape is: six numbers, arranged antisymmetrically. Ask what else you can build from those six numbers using only Lorentz-invariant machinery. Raising and lowering with η\eta gets you FμνF_{\mu\nu} and FμνF^\mu{}_\nu — nothing new. Multiplying FF by itself is quadratic. The only linear operation that produces a genuinely new antisymmetric tensor is contraction with εμναβ\varepsilon^{\mu\nu\alpha\beta}, and it produces exactly one object: the dual. It exchanges the two blocks — the entries that held E⃗/c\vec E/c now hold B⃗\vec B, and those that held B⃗\vec B now hold −E⃗/c-\vec E/c: a rotation by ninety degrees in the (E⃗/c,B⃗)(\vec E/c, \vec B) plane. Twice is 180°180°, which is why ∗∗F=−F{}^{**}F = -F. The payoff is that Maxwell's equations become symmetric in appearance. ∂μFμν=μ0Jν\partial_\mu F^{\mu\nu} = \mu_0 J^\nu and ∂μGμν=0\partial_\mu G^{\mu\nu} = 0 have the identical form; the only difference is what stands on the right. Electricity has a source. Magnetism has a zero. The dual turns that asymmetry from a fact about how the equations are written into a fact about the world.

Figure
THE DUAL IS A 90° ROTATION IN THE (E/c, B) PLANE. F holds E/c in its time-space slots and B in its space-space slots; G = ½ε F swaps them, sending E → cB and B → −E/c. Do it twice and you have rotated 180°: **F = −F. The two Maxwell equations then have IDENTICAL form — ∂_μF^{μν} = μ₀J^ν and ∂_μG^{μν} = 0 — and the ONLY difference is what stands on the right. Electric charge exists. Magnetic charge, so far as anyone has ever measured, does not.
G^μν = ½ ε^μναβ F_αβ · ε⁰¹²³ = +1 · signature (+,−,−,−) F^μν — E in the time-space slots, B in the space-space slots ⎡ 0 −Ex/c −Ey/c −Ez/c ⎤ ⎢ Ex/c 0 −Bz By ⎥ ⎢ Ey/c Bz 0 −Bx ⎥ ⎣ Ez/c −By Bx 0 ⎦ G^μν — the SAME table with E → cB and B → −E/c ⎡ 0 −Bx −By −Bz ⎤ ⎢ Bx 0 Ez/c −Ey/c ⎥ ⎢ By −Ez/c 0 Ex/c ⎥ ⎣ Bz Ey/c −Ex/c 0 ⎦ B E/c dual = 90° E/c → B , B → −E/c twice = 180° ⇒ **F = −F INHOMOGENEOUS — has a source ∂_μ F^μν = μ₀ J^ν ν=0 → Gauss · ν=i → Ampère–Maxwell electric charge EXISTS HOMOGENEOUS — the source is ZERO ∂_μ G^μν = 0 (would be: J_m^ν / c) ν=0 → ∇·B = 0 · ν=i → Faraday magnetic charge: NEVER OBSERVED

The core equations

The construction, and the combinatorial factor that everyone forgets:

  Gμν  ≡  12 εμναβ Fαβ   ,ε0123=+1(1)\boxed{\;G^{\mu\nu} \;\equiv\; \tfrac12\,\varepsilon^{\mu\nu\alpha\beta}\,F_{\alpha\beta}\;}\ ,\qquad \varepsilon^{0123} = +1 \qquad (1)

where: εμναβ\varepsilon^{\mu\nu\alpha\beta} — the totally antisymmetric Levi-Civita symbol; Fαβ=ηαμηβνFμνF_{\alpha\beta} = \eta_{\alpha\mu}\eta_{\beta\nu}F^{\mu\nu}. Physical reading: the sum over (α,β)(\alpha,\beta) counts every independent pair twice — swapping α↔β\alpha\leftrightarrow\beta flips the sign of ε\varepsilon and of FF, so the contributions reinforce. The 12\tfrac12 divides that out, giving G01=F23G^{01} = F_{23}. Two conventions are in play — signature and the sign of ε0123\varepsilon^{0123} — and flipping either flips GG. Written out, GG is FF with the fields rotated into one another:

Gμν=(0−Bx−By−BzBx0Ez/c−Ey/cBy−Ez/c0Ex/cBzEy/c−Ex/c0)⟺  E⃗→cB⃗ ,B⃗→− E⃗c  (2)G^{\mu\nu} = \begin{pmatrix} 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{pmatrix} \qquad\Longleftrightarrow\qquad \boxed{\;\vec E \to c\vec B\ ,\quad \vec B \to -\,\frac{\vec E}{c}\;} \qquad (2)

where: compare entry by entry with FμνF^{\mu\nu}, whose 0i0i slots hold −Ei/c-E_i/c and whose ijij slots hold −εijkBk-\varepsilon^{ijk}B_k. Physical reading: GG is not a new field. It is the same six numbers, re-slotted — which is why duality is a symmetry of the vacuum equations and fails once you put a charge on the right. The two source-free Maxwell equations, in one line — and its equivalent cyclic (Bianchi) form, which is an identity:

  ∂μGμν=0    ⟺  ∇⋅B⃗=0 (ν=0),  ∇×E⃗=−∂B⃗∂t (ν=i)  ⟺    ∂λFμν+∂μFνλ+∂νFλμ=0  (3)\boxed{\;\partial_\mu G^{\mu\nu} = 0\;} \iff \nabla\cdot\vec B = 0\ (\nu=0),\ \ \nabla\times\vec E = -\frac{\partial\vec B}{\partial t}\ (\nu=i) \iff \boxed{\;\partial_\lambda F_{\mu\nu} + \partial_\mu F_{\nu\lambda} + \partial_\nu F_{\lambda\mu} = 0\;} \qquad (3)

where: ∂μ=(1c∂t,∇)\partial_\mu = (\tfrac1c\partial_t, \nabla); AμA_\mu is the four-potential. Physical reading: one four-vector equation, four components, two familiar laws — and not a dynamical law. Substitute Fμν=∂μAν−∂νAμF_{\mu\nu} = \partial_\mu A_\nu - \partial_\nu A_\mu into the cyclic sum and every term cancels against another, because ∂μ∂ν=∂ν∂μ\partial_\mu\partial_\nu = \partial_\nu\partial_\mu. The homogeneous equations carry no content beyond "a potential exists". And the equation nobody has ever needed — the article's owned result, with its coefficient derived:

  ∂μGμν=1c Jmν   ,Jmν=(cρm, J⃗m)⟺∇⋅B⃗=ρm ,∇×E⃗=− ∂B⃗∂t−J⃗m(4)\boxed{\;\partial_\mu G^{\mu\nu} = \frac{1}{c}\,J_m^\nu\;}\ ,\quad J_m^\nu = (c\rho_m,\ \vec J_m) \qquad\Longleftrightarrow\qquad \nabla\cdot\vec B = \rho_m\ ,\quad \nabla\times\vec E = -\,\frac{\partial\vec B}{\partial t} - \vec J_m \qquad (4)

where: ρm\rho_m — magnetic charge density in Wb m−3\mathrm{Wb\,m^{-3}}; J⃗m\vec J_m — magnetic current density in V m−2\mathrm{V\,m^{-2}}. Physical reading: the coefficient is not free. It is fixed the moment you choose units for ρm\rho_m, and drops out of ∂μGμ0=∇⋅B⃗\partial_\mu G^{\mu 0} = \nabla\cdot\vec B by inspection. Observation says Jmν=0J_m^\nu = 0, so (4) collapses to (3) — an experimental result, not a definition.

Building the dual, and reading Maxwell off it

Derivation AThe dual tensor, its factor of one half, and what it does to E and B

The goal is to construct, from FF alone, a second antisymmetric tensor that is not just FF in disguise. There is essentially one way: the Levi-Civita symbol is the unique totally antisymmetric invariant of the proper orthochronous Lorentz group, so contracting with it is the only linear operation available. Everything else — raising, lowering, tracing — either gives back FF or gives zero. The only freedom is the prefactor, and it is fixed by counting.

Step 1. Write Gμν=k εμναβFαβG^{\mu\nu} = k\,\varepsilon^{\mu\nu\alpha\beta}F_{\alpha\beta} and count the double-counting. For fixed (μ,ν)(\mu,\nu) the sum runs over both (α,β)(\alpha,\beta) and (β,α)(\beta,\alpha), and swapping gives εμνβαFβα=(−εμναβ)(−Fαβ)=εμναβFαβ\varepsilon^{\mu\nu\beta\alpha}F_{\beta\alpha} = (-\varepsilon^{\mu\nu\alpha\beta})(-F_{\alpha\beta}) = \varepsilon^{\mu\nu\alpha\beta}F_{\alpha\beta} — the two sign flips reinforce, so each independent pair contributes twice with the same sign. Hence k=12k = \tfrac12. Why: two antisymmetric objects contracted on two indices always double-count. Drop the 12\tfrac12 and every subsequent factor is wrong by two.

Step 2. Compute the entries. With ε0123=+1\varepsilon^{0123} = +1 and the derived components F0i=+Ei/cF_{0i} = +E_i/c, Fij=−εijkBkF_{ij} = -\varepsilon_{ijk}B_k: G^{01} = \tfrac12\,\varepsilon^{01\alpha\beta}F_{\alpha\beta} = \varepsilon^{0123}F_{23} = F_{23} = -B_x,\qquad G^{12} = \varepsilon^{1230}F_{30} = -F_{03} = E_z/c. Why: ε1230\varepsilon^{1230} is three transpositions from ε0123\varepsilon^{0123}, so it is −1-1. Getting that sign right is the entire difficulty of the calculation.

Step 3. Read off the pattern. GG's time-space slots hold −Bi-B_i where FF's held −Ei/c-E_i/c; GG's space-space slots hold +εijkEk/c+\varepsilon^{ijk}E_k/c where FF's held −εijkBk-\varepsilon^{ijk}B_k. Matching GG to FF's shape gives \boxed{\vec E \to c\vec B\ ,\qquad \vec B \to -\vec E/c}. Why: this is the substitution rule and the whole reason the dual is useful — a 90°90° rotation in the (E⃗/c,B⃗)(\vec E/c, \vec B) plane. Apply it twice and you have ∗∗F=−F{}^{**}F = -F.

Step 4. Note what GG is not: a tensor. ε\varepsilon carries a fixed orientation of spacetime, so under parity (det⁡Λ=−1\det\Lambda = -1) GG picks up an extra minus sign relative to a genuine tensor. GG is a pseudotensor. Why: this is the index-notation statement that B⃗\vec B is axial while E⃗\vec E is polar. The pseudo-ness is physics, not bookkeeping. *Verified in verify_aq-sr-dual-field-tensor.py: FμνF^{\mu\nu} is first rebuilt from Aμ=(ϕ/c,A⃗)A^\mu = (\phi/c, \vec A) via Fμν=∂μAν−∂νAμF^{\mu\nu} = \partial^\mu A^\nu - \partial^\nu A^\mu, so the component dictionary is derived rather than assumed; then G=12εFG = \tfrac12\varepsilon F is contracted explicitly and matched entry-by-entry against the matrix printed in equation (2); and GG is shown to equal FF under E⃗→cB⃗\vec E\to c\vec B, B⃗→−E⃗/c\vec B\to-\vec E/c.*

Go deeper · Derivation B · 2 minThe vanishing divergence IS the homogeneous pair, and the Bianchi identity makes it automatic
Derivation BThe vanishing divergence IS the homogeneous pair, and the Bianchi identity makes it automatic

Two things must be shown, and they differ in kind. First, that ∂μGμν=0\partial_\mu G^{\mu\nu} = 0 contains ∇⋅B⃗=0\nabla\cdot\vec B = 0 and Faraday's law — a computation. Second, that it is an identity rather than a law — a structural fact about potentials, and the more important of the two: it tells you the homogeneous equations are not really equations at all.

Step 1. Take ν=0\nu = 0. Only spatial derivatives survive (G00=0G^{00} = 0), and Gi0=+BiG^{i0} = +B_i from equation (2), so \partial_\mu G^{\mu 0} = \partial_i G^{i0} = \partial_i B_i = \nabla\cdot\vec B = 0. Why: <strong>no magnetic monopoles</strong>, straight out of the time slot. The vanishing right-hand side is doing the work.

Step 2. Take ν=1\nu = 1: \partial_\mu G^{\mu 1} = \tfrac1c\partial_t G^{01} + \partial_2 G^{21} + \partial_3 G^{31} = -\tfrac1c\Big[\partial_t B_x + \big(\partial_y E_z - \partial_z E_y\big)\Big] = -\tfrac1c\Big[\partial_t B_x + (\nabla\times\vec E)_x\Big] = 0, so ∇×E⃗=−∂tB⃗\nabla\times\vec E = -\partial_t\vec B. Why: <strong>Faraday&#39;s law of induction</strong>, from the spatial slots. The overall −1/c-1/c is common to all three components; the physics is in the bracket.

Step 3. Now show it is an identity. Insert the definition (with ε\varepsilon constant), then substitute Fαβ=∂αAβ−∂βAαF_{\alpha\beta} = \partial_\alpha A_\beta - \partial_\beta A_\alpha: \partial_\mu G^{\mu\nu} = \tfrac12\,\varepsilon^{\mu\nu\alpha\beta}\,\partial_\mu F_{\alpha\beta} = \tfrac12\,\varepsilon^{\mu\nu\alpha\beta}\big(\partial_\mu\partial_\alpha A_\beta - \partial_\mu\partial_\beta A_\alpha\big) = \varepsilon^{\mu\nu\alpha\beta}\,\partial_\mu\partial_\alpha A_\beta \;\equiv\; 0 , because ∂μ∂α\partial_\mu\partial_\alpha is symmetric in (μ,α)(\mu,\alpha) while εμναβ\varepsilon^{\mu\nu\alpha\beta} is antisymmetric in it, and symmetric contracted with antisymmetric on the same pair vanishes. Why: <strong>this is the whole story.</strong> No physics went in — only the commutativity of partial derivatives. The two homogeneous Maxwell equations are guaranteed the moment you write F=dAF = \mathrm{d}A.

Step 4. Recognise the equivalent form. Contracting the cyclic sum ∂λFμν+∂μFνλ+∂νFλμ\partial_\lambda F_{\mu\nu} + \partial_\mu F_{\nu\lambda} + \partial_\nu F_{\lambda\mu} with ε\varepsilon reproduces 3×3\times the divergence above, and ε\varepsilon is non-degenerate, so ∂μGμν=0  ⟺  ∂[λFμν]=0\partial_\mu G^{\mu\nu} = 0 \iff \partial_{[\lambda}F_{\mu\nu]} = 0. Why: the cyclic form is the <strong>Bianchi identity</strong>, and it survives to curved spacetime — the connection terms cancel by antisymmetry, which is why electromagnetism couples to gravity painlessly.

Step 5. Draw the conclusion. The homogeneous equations have no source because they cannot have one, so long as FF derives from a potential. To give magnetism a source you must give up the global existence of AμA^\mu — exactly what Dirac did in 1931, with a string singularity. That is why the monopole is a statement about topology rather than about an extra Lagrangian term. Why: this is the deep reason the two Maxwell pairs are not symmetric, and no amount of pretty notation can make them so. *Verified in verify_aq-sr-dual-field-tensor.py: ∂μGμ0=∇⋅B⃗\partial_\mu G^{\mu 0} = \nabla\cdot\vec B and ∂μGμi=−1c[(∇×E⃗)i+∂tBi]\partial_\mu G^{\mu i} = -\tfrac1c[(\nabla\times\vec E)_i + \partial_t B_i], symbolically, with E⃗\vec E and B⃗\vec B arbitrary functions of (t,x⃗)(t,\vec x); and εμναβ∂μ∂αAβ=0\varepsilon^{\mu\nu\alpha\beta}\partial_\mu\partial_\alpha A_\beta = 0 by explicit summation over all 444^4 index combinations for an arbitrary Aμ(t,x⃗)A_\mu(t,\vec x).*

Go deeper · Derivation C · 2 minThe double dual, the pseudoscalar, and the monopole source term whose coefficient the base gets wrong
Derivation CThe double dual, the pseudoscalar, and the monopole source term whose coefficient the base gets wrong

Three loose ends, the third a repair. The double dual pins down the sign that distinguishes Minkowski from Euclidean geometry; the pseudoscalar FμνGμνF_{\mu\nu}G^{\mu\nu} is the invariant only the dual can see; and the monopole equation's coefficient is a number the base article simply asserts — wrongly.

Step 1. Dualise twice, using εμναβεμνρσ=−2(δραδσβ−δσαδρβ)\varepsilon^{\mu\nu\alpha\beta}\varepsilon_{\mu\nu\rho\sigma} = -2\big(\delta^\alpha_\rho\delta^\beta_\sigma - \delta^\alpha_\sigma\delta^\beta_\rho\big): {}^{}F^{\rho\sigma} = \tfrac14\,\varepsilon^{\rho\sigma\mu\nu}\varepsilon_{\mu\nu\alpha\beta}F^{\alpha\beta} = -\tfrac12\big(F^{\rho\sigma} - F^{\sigma\rho}\big) = \boxed{-F^{\rho\sigma}} . *Why: the minus traces to det⁡η=−1\det\eta = -1 — a signature** effect, not a convention. In Euclidean four-space the Hodge star squares to +1+1 and real self-dual fields exist; in Minkowski space it squares to −1-1 and they do not. That is why instantons live in Euclidean signature.*

Step 2. Contract FF with its dual. Every surviving product pairs an EE-component with the same BB-component: \boxed{\;F_{\mu\nu}G^{\mu\nu} = -\,\frac{4}{c}\,\vec E\cdot\vec B\;}. Why: <strong>this invariant is invisible to FF alone.</strong> No contraction of FF with itself and the metric can produce E⃗⋅B⃗\vec E\cdot\vec B — you need the orientation ε\varepsilon supplies. It is a <strong>pseudoscalar</strong>: unchanged by boosts, sign-flipped by a mirror.

Step 3. Now the monopole. Write ∂μGμν=κ Jmν\partial_\mu G^{\mu\nu} = \kappa\,J_m^\nu with Jmν=(cρm,J⃗m)J_m^\nu = (c\rho_m, \vec J_m) and κ\kappa unknown, then fix it from the time component. The previous derivation gave ∂μGμ0=∇⋅B⃗\partial_\mu G^{\mu 0} = \nabla\cdot\vec B exactly. If magnetic charge is measured in webers, so the monopole Gauss law reads ∇⋅B⃗=ρm\nabla\cdot\vec B = \rho_m with ρm\rho_m in Wb m−3\mathrm{Wb\,m^{-3}}, then κ (cρm)=ρm\kappa\,(c\rho_m) = \rho_m and \boxed{\kappa = \tfrac1c}. Why: no fitting, no guessing. The left side is already computed; the right side is a choice of units; κ\kappa is whatever makes them agree. It is <strong>not</strong> a free parameter.

Step 4. Check the spatial components, which are now a prediction, not an input. With κ=1/c\kappa = 1/c the ν=i\nu = i components read −1c[(∇×E⃗)i+∂tBi]=1cJm,i-\tfrac1c[(\nabla\times\vec E)_i + \partial_t B_i] = \tfrac1c J_{m,i}, i.e. \nabla\times\vec E = -\frac{\partial\vec B}{\partial t} - \vec J_m , exactly the magnetic-current Faraday law of the standard duality-symmetric Maxwell equations. Why: κ\kappa was fixed by ν=0\nu = 0 alone, and the three ν=i\nu = i components then came out right on their own. That is the check, and it is what distinguishes a derived coefficient from a guessed one.

Step 5. Refuse the base's coefficient. The base prints ∂μGμν=μ0c Jmν\partial_\mu G^{\mu\nu} = \mu_0 c\,J_m^\nu. Dimensional analysis kills it: [∂G]=T m−1[\partial G] = \mathrm{T\,m^{-1}} and [μ0c]=T m2 A−1 s−1[\mu_0 c] = \mathrm{T\,m^2\,A^{-1}\,s^{-1}}, so JmνJ_m^\nu would carry units of A s m−3=C m−3\mathrm{A\,s\,m^{-3}} = \mathrm{C\,m^{-3}} — an electric charge density. Not a magnetic current at all. Why: the coefficient on the right of a field equation is the exchange rate between the units on the two sides, checkable in one line. (In the "ampere-metre" convention, ∇⋅B⃗=μ0ρm\nabla\cdot\vec B = \mu_0\rho_m, it would be κ=μ0/c\kappa = \mu_0/c — still not μ0c\mu_0 c.) *Verified in verify_aq-sr-dual-field-tensor.py: ∗∗F=−F{}^{**}F = -F by explicit double contraction; FμνGμν=−(4/c)E⃗⋅B⃗F_{\mu\nu}G^{\mu\nu} = -(4/c)\vec E\cdot\vec B symbolically; κ=1/c\kappa = 1/c derived by solving ∂μGμ0=κcρm\partial_\mu G^{\mu 0} = \kappa c\rho_m against ∇⋅B⃗=ρm\nabla\cdot\vec B = \rho_m, then independently confirmed against the ν=i\nu = i components; and the base's μ0c\mu_0 c refuted dimensionally.*

Limiting cases

CaseResultMeaning
B⃗=0\vec B = 0 (pure electric)G0i=0G^{0i} = 0,   Gij=εijkEk/c\;G^{ij} = \varepsilon^{ijk}E_k/cthe dual of a pure E⃗\vec E is a monopole field
plane wave (E=cBE = cB, E⃗⊥B⃗\vec E\perp\vec B)FμνFμν=0F_{\mu\nu}F^{\mu\nu} = 0,   FμνGμν=0\;F_{\mu\nu}G^{\mu\nu} = 0a null field: FF and GG carry the same information
F=dAF = \mathrm{d}A∂μGμν≡0\partial_\mu G^{\mu\nu} \equiv 0an identity. Homogeneous Maxwell is free
dual twice∗∗F=−F{}^{**}F = -Fa 180°180° duality rotation. The −- is det⁡η=−1\det\eta = -1
parity (x⃗→−x⃗\vec x\to-\vec x)G→− (tensor image of G)G \to -\,(\text{tensor image of } G)GG is a pseudotensor; B⃗\vec B is axial, E⃗\vec E is polar
magnetic charge present∂μGμν=Jmν/c\partial_\mu G^{\mu\nu} = J_m^\nu/cwhat nobody has ever seen. The zero is data
ε0123→−1\varepsilon^{0123} \to -1G→−GG \to -Gconventions matter: state yours, or compare nonsense

Worked examples

1 · Dualise an electron, and discover why it is not a monopole

Take a single electron. At r=1.00 A˚=1.00×10−10 mr = 1.00\ \mathrm{\AA} = 1.00\times10^{-10}\ \mathrm{m}, compute its Coulomb field, apply the duality rotation, interpret the result, and compare with Dirac's quantisation condition.

Step 1 — identify. A point charge at rest has B⃗=0\vec B = 0, so FF has entries only in its time-space slots and the dual only in its space-space slots. Equation (2) gives (E⃗G,B⃗G)=(cB⃗,−E⃗/c)=(0⃗, −E⃗/c)(\vec E_G, \vec B_G) = (c\vec B, -\vec E/c) = (\vec 0,\ -\vec E/c): purely magnetic, and radial.

Step 2 — the Coulomb field. With 1/4πε0=8.98755×109 N m2 C−21/4\pi\varepsilon_0 = 8.98755\times10^{9}\ \mathrm{N\,m^2\,C^{-2}}, E = \frac{e}{4\pi\varepsilon_0 r^2} = \frac{(8.98755\times10^{9})(1.602177\times10^{-19})}{(1.00\times10^{-10})^2} = \mathbf{1.43996\times10^{11}\ \mathrm{V\,m^{-1}}}.

Step 3 — the dual field. |\vec B_G| = \frac{E}{c} = \frac{1.43996\times10^{11}}{2.99792\times10^{8}} = \mathbf{480.32\ T}. A radial magnetic field of 480 tesla at one ångström, so ∇⋅B⃗G≠0\nabla\cdot\vec B_G \ne 0 at the origin. The dual of a point electric charge is a magnetic monopole.

Step 4 — its magnetic charge. Integrating ∇⋅B⃗G=ρm\nabla\cdot\vec B_G = \rho_m over a sphere gives g=μ0ecg = \mu_0 e c: g = (1.256637\times10^{-6})(1.602177\times10^{-19})(2.997925\times10^{8}) = \mathbf{6.0359\times10^{-17}\ Wb}.

Step 5 — the boundary number, and the punchline. Dirac's condition says the smallest magnetic charge consistent with quantum mechanics is g_D = \frac{h}{e} = \frac{6.62607\times10^{-34}}{1.602177\times10^{-19}} = \mathbf{4.13567\times10^{-15}\ Wb}, exactly twice the superconducting flux quantum Φ0=h/2e\Phi_0 = h/2e — and it exceeds our dual electron's charge by \frac{g_D}{\mu_0 e c} = \frac{h}{\mu_0 e^2 c} = \frac{1}{2\alpha} = \mathbf{68.518}. A dual electron is 68.5 times too weak to be a legal monopole. Duality is a symmetry of the classical equations; quantum mechanics then refuses most of its results. That gap — the fine-structure constant, appearing in the middle of a purely classical construction — is why monopoles, if they exist, are enormously strongly coupled.

2 · The invariants of the field and of its dual

A field has E⃗\vec E and B⃗\vec B <strong>parallel</strong>, with E=1.00×106 V m−1E = 1.00\times10^{6}\ \mathrm{V\,m^{-1}} and B=0.500 TB = 0.500\ \mathrm{T}. Compute both invariants for the field and for its dual, and say what a boost can and cannot do.

Step 1 — identify. E⃗∥B⃗\vec E \parallel \vec B, so E⃗⋅B⃗=EB≠0\vec E\cdot\vec B = EB \ne 0: the pseudoscalar invariant is alive, and only the dual can see it.

Step 2 — the scalar. F_{\mu\nu}F^{\mu\nu} = 2\left(B^2 - \frac{E^2}{c^2}\right) = 2\left(0.250000 - \frac{1.00\times10^{12}}{8.98755\times10^{16}}\right) = 2(0.250000 - 1.11265\times10^{-5}) = \mathbf{0.499978\ T^2}. Positive: magnetically dominated, in every frame.

Step 3 — the pseudoscalar, which needs GG. E⃗⋅B⃗=(1.00×106)(0.500)=5.00×105 V T m−1\vec E\cdot\vec B = (1.00\times10^{6})(0.500) = 5.00\times10^{5}\ \mathrm{V\,T\,m^{-1}}, so F_{\mu\nu}G^{\mu\nu} = -\frac{4}{c}\,\vec E\cdot\vec B = -\frac{4 \times 5.00\times10^{5}}{2.99792\times10^{8}} = \mathbf{-6.6713\times10^{-3}\ T^2}. Non-zero.

Step 4 — the dual's invariants, and the boundary number. Applying E⃗→cB⃗\vec E\to c\vec B, B⃗→−E⃗/c\vec B\to -\vec E/c gives EG=cB=1.49896×108 V m−1E_G = cB = 1.49896\times10^{8}\ \mathrm{V\,m^{-1}} and BG=−E/c=−3.33564×10−3 TB_G = -E/c = -3.33564\times10^{-3}\ \mathrm{T}, whose scalar invariant is 2(BG2−EG2/c2)=−0.499978 T22(B_G^2 - E_G^2/c^2) = \mathbf{-0.499978\ T^2} — exactly minus the original, as a 90°90° rotation must give. And because FμνGμν≠0F_{\mu\nu}G^{\mu\nu} \ne 0 is frame-independent, no boost whatsoever can remove either field: a pure electric field has E⃗⋅B⃗=0\vec E\cdot\vec B = 0, a pure magnetic field has E⃗⋅B⃗=0\vec E\cdot\vec B = 0, and ours has 5.00×1055.00\times10^{5}. The best any observer can do is make E⃗\vec E and B⃗\vec B parallel — which they already are. The dual told us this before a single boost was attempted.

Physics of the portrait

The dual tensor exists because someone asked what the other two Maxwell equations are for. The answer — that they are an identity, and that the thing which would break them has never been found — is the story of the monopole.

Origins · 2 minMagnetic monopole and how this was discovered
MM
Magnetic monopole
1931–present · Dirac's quantisation · 't Hooft–Polyakov · Cabrera's Valentine's Day event · MoEDAL

The magnetic monopole is the particle that would put a source on the right of ∂μGμν=0\partial_\mu G^{\mu\nu} = 0, and its theoretical pedigree is far better than its experimental one. Paul Dirac's 1931 paper "Quantised Singularities in the Electromagnetic Field" (Proc. Roy. Soc. A 133, 60) showed something remarkable: the existence of even one monopole anywhere in the universe would force all electric charge to be quantised, in multiples of a unit fixed by eg=nℏc/2eg = n\hbar c/2. Charge quantisation is among the most precisely tested facts in physics and has no other explanation in ordinary electrodynamics. Dirac's argument turns on the fact that a monopole's vector potential cannot be defined globally: it must carry a singular "string", and demanding the string be unobservable to any quantum wavefunction is what quantises egeg. The minimum magnetic charge is gD=h/e=4.14×10−15 Wbg_D = h/e = 4.14\times10^{-15}\ \mathrm{Wb} — exactly two superconducting flux quanta — with a coupling 68.568.5 times the electron's, which makes a monopole spectacularly ionising. In 1974 Gerard 't Hooft and Alexander Polyakov independently showed that monopoles are not optional in grand unified theories: any GUT breaking a simple group to a U(1)U(1) predicts them near 1016 GeV10^{16}\ \mathrm{GeV} — a large part of why inflation was invented, to dilute them away. The experimental record is one long null result with a single famous asterisk: on the night of 14 February 1982, Blas Cabrera's superconducting loop at Stanford recorded exactly one event of the right size — the "Valentine's Day monopole". It never repeated. So far, the right-hand side of ∂μGμν=0\partial_\mu G^{\mu\nu} = 0 remains, stubbornly and informatively, zero.

History
1865
— Maxwell's "A Dynamical Theory of the Electromagnetic Field" presents the field equations as a system, two of them carrying sources and two not.
1873
— Maxwell's Treatise on Electricity and Magnetism codifies the asymmetry without explaining it.
1908
— Minkowski, in "Die Grundgleichungen für die elektromagnetischen Vorgänge in bewegten Körpern", casts Maxwell's equations in four-dimensional tensor form; his antisymmetric objects are called "six-vectors" for the next decade, and the dual is implicit in his formalism.
1931
— Dirac's "Quantised Singularities in the Electromagnetic Field" (Proc. Roy. Soc. A 133, 60) shows that a single magnetic monopole would force the quantisation of all electric charge, with eg = n ℏc/2.
1974
— 't Hooft and Polyakov independently show that grand unified theories necessarily contain magnetic monopoles, with masses near the unification scale.
1982
— On 14 February, Blas Cabrera's SQUID loop at Stanford records a single candidate monopole event, the "Valentine's Day monopole"; it is never repeated and remains unconfirmed.
2016
— MoEDAL at the LHC publishes its first monopole search, setting mass limits for magnetic charges up to several Dirac units; the result, like every result before it, is null.

Going further

Texts. Jackson, Classical Electrodynamics, 3rd ed., §§11.9 and 6.11 · Griffiths, Introduction to Electrodynamics, 4th ed., §12.3.4 · Rindler, Introduction to Special Relativity, 2nd ed., ch. 7 · Dirac, Proc. Roy. Soc. A 133, 60 (1931) · Cabrera, Phys. Rev. Lett. 48, 1378 (1982).

Provenance. Signature (+,−,−,−)(+,-,-,-), ε0123=+1\varepsilon^{0123} = +1; p⋅p=+m2c2p\cdot p = +m^2c^2; mm is the invariant rest mass, and no velocity-dependent mass appears anywhere. Base-defect repairs (logged): (i) the base prints ∂μGμν=μ0cJmν\partial_\mu G^{\mu\nu} = \mu_0 c J_m^\nu; that coefficient is dimensionally impossible and is refused — the correct κ=1/c\kappa = 1/c is derived here and independently confirmed against the spatial components. (ii) The base's "compare with FμνF^{\mu\nu}" note mixes two incompatible sign conventions in one sentence, quoting F0i=+Ei/cF^{0i} = +E_i/c beside Fij=−εijkBkF^{ij} = -\varepsilon^{ijk}B_k; its own printed GG matrix does not follow from the components it states. (iii) Three of its ten equation cards are empty or contain the literal string `tex_note_placeholder`. Everything here is recomputed in verify_aq-sr-dual-field-tensor.py, which parses this article's own printed digits.

Scope. The construction of FμνF^{\mu\nu} belongs to The Electromagnetic Field Tensor; the covariant Maxwell system to Covariant Maxwell Equations; and the classification of fields by the two invariants belongs to Electromagnetic Field Invariants. This page owns the dual, its divergence, and the monopole source term.

Builds on. The Electromagnetic Field Tensor · Covariant Maxwell Equations · The Four-Current Density · The Four-Potential. Leads to. Electromagnetic Field Invariants · Transformation of Electric and Magnetic Fields · The Electromagnetic Stress-Energy Tensor.

Is G a different physical field from F?

No. GG contains exactly the same six numbers as FF, re-slotted; knowing FF everywhere means knowing GG everywhere and vice versa (up to the sign that ∗∗F=−F{}^{**}F = -F records). What GG buys is not new information but new contractions: FμνGμνF_{\mu\nu}G^{\mu\nu} is an invariant that no amount of contracting FF with itself and with η\eta can produce, because it needs the orientation ε\varepsilon carries.

If the homogeneous equations are just an identity, why bother writing them down?

Because "there exists a four-potential" is a physical assumption, not a tautology, and it is one a monopole would violate. The identity says: given AμA^\mu, you get ∇⋅B⃗=0\nabla\cdot\vec B = 0 and Faraday for free. The empirical content is the antecedent. In a spacetime containing a monopole, AμA^\mu cannot be defined globally — Dirac's string is the scar — and the "identity" fails, not because derivatives stopped commuting but because the object being differentiated does not exist.

Why does the pseudoscalar invariant matter if it is a total derivative?

Classically it does not — adding θ FμνGμν\theta\,F_{\mu\nu}G^{\mu\nu} to the Lagrangian changes no equation of motion. But "total derivative" means "sensitive only to boundary conditions and topology", and both matter quantum-mechanically. In QCD the analogous term produces observable CP violation (mysteriously absent — the strong-CP problem), and promoting θ\theta to a dynamical field is the axion.

What actually stops monopoles from existing?

Nothing in the classical equations. Duality-symmetric electrodynamics is perfectly consistent — the equations are more symmetric with monopoles than without. The obstructions are (i) experimental, nobody has found one, and (ii) structural: with a monopole present, FF is no longer dA\mathrm{d}A globally, so the potential formulation — and with it gauge theory, the path integral, and essentially all of QED as practised — needs surgery. Dirac's string, and Wu and Yang's later patchwork of overlapping gauge charts, are that surgery. The price is paid in topology, not in an extra term.

Does the sign of ε⁰¹²³ matter physically?

Not for any observable, but enormously for cross-checking formulas. Flipping ε0123\varepsilon^{0123} from +1+1 to −1-1 flips GG, hence FμνGμνF_{\mu\nu}G^{\mu\nu}, hence the sign in front of E⃗⋅B⃗\vec E\cdot\vec B. Two textbooks can print apparently contradictory results while agreeing completely on the physics. Always state both conventions — signature and ε0123\varepsilon^{0123} — before comparing. Here: (+,−,−,−)(+,-,-,-) and ε0123=+1\varepsilon^{0123} = +1.

How strongly would a monopole interact if it existed?

Ferociously. The Dirac charge corresponds to a coupling gD/(μ0ec)=1/2α=68.5g_D/(\mu_0 e c) = 1/2\alpha = 68.5 times the electron's, so the "magnetic fine-structure constant" is αm=1/(4α)≈34\alpha_m = 1/(4\alpha) \approx 34 — strongly coupled, and perturbation theory in it is worthless. A monopole passing through matter would ionise thousands of times more heavily than any known particle, which is precisely the signature MoEDAL's plastic track detectors are built to catch. If one goes through, you will not miss it.

Before this

Where this leads