Skip to content

Recover Maxwell operations and stress

The electromagnetic potential, field, current and stress occupy distinct levels. The same geometry carries their derivative operations, gauge fibre and transfer identities, including the information that stress alone cannot recover.

The electromagnetic potential, field, current and stress occupy distinct levels. The same geometry carries their derivative operations, gauge fibre and transfer identities, including the information that stress alone cannot recover.

12. Electromagnetic fields, Maxwell operations and stress

Einstein’s §20 starts with a potential, constructs its antisymmetric derivative, raises indices with the metric, supplies a current and derives the electromagnetic stress through a transfer identity. Its heading refers to vacuum, but its displayed sourced equation includes a current. We retain that current as an independent input and keep the original signs and index order.

Primitive coverage and potential contact

Add a calibrated covector α\alpha, an independent antisymmetric array fabf_{ab} and a calibrated current jaj^a. Define

φ=ETα,F=ETfE,J=E−1j.(12.1)\varphi=E^{\mathsf T}\alpha,\qquad \mathsf F=E^{\mathsf T}fE,\qquad J=E^{-1}j. \tag{12.1}

The reverse maps are α=E−Tφ\alpha=E^{-\mathsf T}\varphi, f=E−TFE−1f=E^{-\mathsf T}\mathsf FE^{-1} and j=EJj=EJ. Together with the local metric section they cover every local target primitive tuple. Under a frame change transform α′=Λ−Tα\alpha'=\Lambda^{-\mathsf T}\alpha, f′=Λ−TfΛ−1f'=\Lambda^{-\mathsf T}f\Lambda^{-1} and j′=Λjj'=\Lambda j so that all readouts are preserved.

Initially retain the contact and Maxwell residuals

kab=fab−(Dbαa−Daαb),habc=Dafbc+Dbfca+Dcfab,ma=Dbfab−ja.(12.2)k_{ab}=f_{ab}-(D_b\alpha_a-D_a\alpha_b),\quad h_{abc}=D_af_{bc}+D_bf_{ca}+D_cf_{ab},\quad m^a=D_bf^{ab}-j^a. \tag{12.2}

The derivative intertwiner gives

EaμEbν(Dbαa−Daαb)=∂νφμ−∂μφν.(12.3)E^a{}_{\mu}E^b{}_{\nu}(D_b\alpha_a-D_a\alpha_b) =\partial_\nu\varphi_\mu-\partial_\mu\varphi_\nu. \tag{12.3}

Coordinate connection terms cancel by symmetry. Equivalently, using the frame structure coefficients, Dbαa−Daαb=ebαa−eaαb+ccabαcD_b\alpha_a-D_a\alpha_b=e_b\alpha_a-e_a\alpha_b+c^c{}_{ab}\alpha_c. This recovers Einstein (59), with F=−dφ\mathsf F=-\mathrm{d}\varphi when contact holds. On the raw domain,

F=−dφ+K,dF=dK,K=12kabθa∧θb.(12.4)\mathsf F=-\mathrm{d}\varphi+K,\qquad \mathrm{d}\mathsf F=\mathrm{d}K, \qquad K=\tfrac12k_{ab}\theta^a\wedge\theta^b. \tag{12.4}

Thus zero contact gives the homogeneous equations h=0h=0, or

∂τFρσ+∂ρFστ+∂σFτρ=0,(12.5)\partial_\tau\mathsf F_{\rho\sigma} +\partial_\rho\mathsf F_{\sigma\tau} +\partial_\sigma\mathsf F_{\tau\rho}=0, \tag{12.5}

recovering Einstein (60). Bare frame-component derivatives are insufficient: a position-dependent Lorentz boost of rapidity x1x^1 with constant coordinate potential dx1\mathrm{d}x^1 has actual field zero, but bare antisymmetrized calibrated derivatives give a nonzero 1414 component.

Proposition 12.1 (Local potential section and gauge fibre)

On a star-shaped coordinate neighborhood centered at zero, every smooth closed two-form has an Einstein-convention potential

φμ(x)=−∫01s xνFνμ(sx)ds,F=−dφ.(12.6)\varphi_\mu(x)=-\int_0^1s\,x^\nu\mathsf F_{\nu\mu}(sx)\mathrm{d}s, \qquad \mathsf F=-\mathrm{d}\varphi. \tag{12.6}

Locally two potentials with the same field differ by dχ\mathrm{d}\chi.

Proof.

Differentiate under the finite integral and use the cyclic identity. The antisymmetric derivative becomes the integral of d[s2Fμν(sx)]/ds\mathrm{d}[s^2\mathsf F_{\mu\nu}(sx)]/\mathrm{d}s with the sign specified in (12.3). Smoothness makes the lower endpoint zero; the upper endpoint gives the field. A difference of potentials is a closed one-form, whose radial integral supplies a scalar primitive. In calibrated variables the gauge change is αa↦αa+eaχ\alpha_a\mapsto\alpha_a+e_a\chi.

□\square

This reconstructs the field, not necessarily the original potential representative. It is a local result; global topology, flux sectors and fixed boundary potential data require additional specifications. Gauge comparisons are allowed only when the retained observables and boundary data agree. The primitive target in the original already contains a potential, so target-potential coverage does not require assuming global exactness for every closed field.

Index raising, sourced divergence and charge conservation

The raising operation is

fab=ηacηbdfcd,Fμν=(E−1)μa(E−1)νbfab=qμαqνβFαβ,(12.7)f^{ab}=\eta^{ac}\eta^{bd}f_{cd},\qquad \mathsf F^{\mu\nu}=(E^{-1})^\mu{}_a(E^{-1})^\nu{}_bf^{ab} =q^{\mu\alpha}q^{\nu\beta}\mathsf F_{\alpha\beta}, \tag{12.7}

recovering Einstein (62). Contracting the derivative slot with the second upper field slot gives the exact residual map

(E−1)μama=∇νFμν−Jμ=ν−1∂ν(νFμν)−Jμ.(12.8)\boxed{(E^{-1})^\mu{}_am^a =\nabla_\nu\mathsf F^{\mu\nu}-J^\mu =\nu^{-1}\partial_\nu(\nu\mathsf F^{\mu\nu})-J^\mu.} \tag{12.8}

The connection acting on the first field index contracts a symmetric pair against antisymmetry and vanishes; the remaining trace supplies ∂log⁡ν\partial\log\nu. At m=0,ν=1m=0,\nu=1 this is precisely Einstein (63). It differentiates the second index, in accordance with the original convention.

Multiply (12.8) by ν\nu, differentiate in μ\mu, and commute partial derivatives. Antisymmetry cancels the double derivative, proving

∇μ(Jμ+mμ)=0,mμ=(E−1)μama.(12.9)\nabla_\mu(J^\mu+m^\mu)=0, \qquad m^\mu=(E^{-1})^\mu{}_am^a. \tag{12.9}

On Maxwell’s equation this is Daja=0D_aj^a=0, equivalently ∂μ(νJμ)=0\partial_\mu(\nu J^\mu)=0. A nonconserved prescribed current is a failed candidate, not permission to redefine the source after calculating a field.

The original electric and magnetic component systems

In the unit-volume chart use Einstein’s definitions

(F23,F31,F12)=h,Fi4=ei,(F23,F31,F12)=h′,Fi4=−ei′.(12.10)(\mathsf F_{23},\mathsf F_{31},\mathsf F_{12})=\mathbf h,\quad \mathsf F_{i4}=e_i, \qquad (\mathsf F^{23},\mathsf F^{31},\mathsf F^{12})=\mathbf h',\quad \mathsf F^{i4}=-e'_i. \tag{12.10}

Expand (12.5) and (12.8), with Ji=iiJ^i=i_i and J4=ρelJ^4=\rho_{\rm el}:

∂th+curl⁡e=0,div⁡h=0,(12.11)curl⁡h′−∂te′=i,div⁡e′=ρel.(12.12)\begin{aligned} \partial_t\mathbf h+\operatorname{curl}\mathbf e&=0,& \operatorname{div}\mathbf h&=0,&&\text{(12.11)}\\ \operatorname{curl}\mathbf h'-\partial_t\mathbf e'&=\mathbf i,& \operatorname{div}\mathbf e'&=\rho_{\rm el}. &&\text{(12.12)}\end{aligned}

These recover Einstein (60a), Einstein (60b), Einstein (61), Einstein (63a) and Einstein (64). No extra 4π4\pi or unit coefficient has been inserted. The primed and unprimed arrays need not coincide: even the constant unit-volume calibration E=diag⁡(2,3,1,1/6)E=\operatorname{diag}(2,3,1,1/6) gives

e′=(9e1,4e2,36e3),h′=(h1/9,h2/4,h3/36).(12.13)\mathbf e'=(9e_1,4e_2,36e_3),\qquad \mathbf h'=(h_1/9,h_2/4,h_3/36). \tag{12.13}

The geometry is flat in rescaled coordinates, but its coordinate constitutive relation is not the Cartesian one.

Transfer and the electromagnetic stress calculation

Define

ξb=fbaja,I=fabfab,ξμ=Ebμξb=FμαJα.(12.14)\xi_b=f_{ba}j^a,\qquad I=f_{ab}f^{ab}, \qquad \xi_\mu=E^b{}_{\mu}\xi_b=\mathsf F_{\mu\alpha}J^\alpha. \tag{12.14}

The symbol ξ\xi avoids confusion with the gravitational coupling κ\kappa. The product/cyclic differentiation in Einstein’s §20 leads to

Σab,em=−fbcfac+14δbaI,(12.15)Σabem=−facηcdfbd+14ηabI,(12.16)Tνσ,em=−FσαFνα+14δσνFαβFαβ.(12.17)\begin{aligned} \Sigma^a{}_{b,\rm em}&=-f_{bc}f^{ac}+\tfrac14\delta^a_b I,&&\text{(12.15)}\\ \Sigma^{\rm em}_{ab}&=-f_{ac}\eta^{cd}f_{bd}+\tfrac14\eta_{ab}I,&&\text{(12.16)}\\ T^\nu{}_{\sigma,\rm em} &=-\mathsf F_{\sigma\alpha}\mathsf F^{\nu\alpha} +\tfrac14\delta^\nu_\sigma\mathsf F_{\alpha\beta}\mathsf F^{\alpha\beta}. &&\text{(12.17)}\end{aligned}

Symmetry and zero trace follow algebraically. The tensor multiplication and contraction maps prove the last equality from the first two. It is Einstein (66a).

Theorem 12.2 (Complete electromagnetic transfer identity)

For any smooth antisymmetric raw field and prescribed current in the compatible calculus,

DaΣab,em=ξb+fbcmc+12fcdhbcd.(12.18)\boxed{D_a\Sigma^a{}_{b,\rm em} =\xi_b+f_{bc}m^c+\tfrac12f^{cd}h_{bcd}.} \tag{12.18}

Proof.

Use coordinate slots through the proved intertwiner. Set Hσαβ=∇σFαβ+∇αFβσ+∇βFσαH_{\sigma\alpha\beta}=\nabla_\sigma\mathsf F_{\alpha\beta} +\nabla_\alpha\mathsf F_{\beta\sigma}+\nabla_\beta\mathsf F_{\sigma\alpha}. Contract with Fαβ\mathsf F^{\alpha\beta}. Compatibility makes the first term ∂σI/2\partial_\sigma I/2, and relabelling shows that the other two terms are equal. Hence

Fαβ∇βFσα=−14∂σI+12FαβHσαβ.(12.19)\mathsf F^{\alpha\beta}\nabla_\beta\mathsf F_{\sigma\alpha} =-\tfrac14\partial_\sigma I +\tfrac12\mathsf F^{\alpha\beta}H_{\sigma\alpha\beta}. \tag{12.19}

The product rule gives

Fσα(Jα+mα)=∇β(FσαFαβ)−Fαβ∇βFσα.\mathsf F_{\sigma\alpha}(J^\alpha+m^\alpha) =\nabla_\beta(\mathsf F_{\sigma\alpha}\mathsf F^{\alpha\beta}) -\mathsf F^{\alpha\beta}\nabla_\beta\mathsf F_{\sigma\alpha}.

Insert (12.19). Antisymmetry combines the first two divergence terms into (12.17); moving the homogeneous defect gives the plus sign and one-half coefficient in (12.18). Read out the remaining calibrated slot in reverse. No Maxwell equation was assumed in this calculation.

□\square

On m=h=0m=h=0,

DaΣab,em=ξb,∇νTνσ,em=ξσ=ν−1∂ν(νTνσ,em)−12∂σgαβTemαβ.(12.20)D_a\Sigma^a{}_{b,\rm em}=\xi_b,\qquad \nabla_\nu T^\nu{}_{\sigma,\rm em}=\xi_\sigma =\nu^{-1}\partial_\nu(\nu T^\nu{}_{\sigma,\rm em}) -\tfrac12\partial_\sigma g_{\alpha\beta}T^{\alpha\beta}_{\rm em}. \tag{12.20}

The unit-volume form is Einstein (66); the metric derivative is essential. The raw example E=IE=I, F12=z\mathsf F_{12}=z with other independent components zero and J=0J=0 satisfies the sourced divergence but has H312=1H_{312}=1 and ∇νTνσ=(0,0,z,0)σ\nabla_\nu T^\nu{}_{\sigma}=(0,0,z,0)_\sigma. Omitting the final term in (12.18) would wrongly make this nonclosed field conserved.

In a Cartesian calibrated frame the formulas reduce to

I=2(∣h∣2−∣e∣2),T44=12(∣e∣2+∣h∣2),Ti4=(e×h)i,(12.21)Tij=12(∣e∣2+∣h∣2)δij−eiej−hihj,(ξi)i=13=ρele+i×h,ξ4=−i⋅e.(12.22)\begin{aligned} I&=2(|\mathbf h|^2-|\mathbf e|^2),& T^4{}_4&=\tfrac12(|\mathbf e|^2+|\mathbf h|^2),& T^i{}_4&=(\mathbf e\times\mathbf h)_i,&&\text{(12.21)}\\ T_{ij}&=\tfrac12(|\mathbf e|^2+|\mathbf h|^2)\delta_{ij}-e_ie_j-h_ih_j, &&&\\ (\xi_i)_{i=1}^3&=\rho_{\rm el}\mathbf e+\mathbf i\times\mathbf h,& \xi_4&=-\mathbf i\cdot\mathbf e.&& &&\text{(12.22)}\end{aligned}

This verifies Einstein (65), Einstein (65a) and the original Maxwell–Poynting specialization. The time sign describes loss of field energy when i⋅e>0\mathbf i\cdot\mathbf e>0; spatial covariant and contravariant components must not be identified in this signature.

The zero-transfer branch and material joins

The original text asserts zero transfer for “free” electric masses and then explicitly uses ξ=0\xi=0 when identifying Einstein (66) with the earlier matter-conservation law. The displayed Maxwell equations do not by themselves prove FσμJμ=0\mathsf F_{\sigma\mu}J^\mu=0. We therefore retain the zero-transfer restriction without promoting that prose to a general charge-motion theorem.

For an exact equation-level control, take E=IE=I, φ4=−x2/2\varphi_4=-x^2/2, φi=0\varphi_i=0, F14=x\mathsf F_{14}=x, and J=(0,0,0,1)J=(0,0,0,1). Both Maxwell systems hold, but ξ=(x,0,0,0)≠0\xi=(x,0,0,0)\ne0. This is a Maxwell solution on a fixed background, not a coupled gravitational solution or a model of freely moving charges.

If fluid and electromagnetic stress are joined to an independently supplied stress, retain the interface residual ΔΣ=Σ−Σfl−Σem\Delta\Sigma=\Sigma-\Sigma^{\rm fl}-\Sigma^{\rm em}. On zero interface discrepancy, define K=G−κ(Σfl+Σem)\mathscr K=\mathscr G-\kappa(\Sigma^{\rm fl}+\Sigma^{\rm em}). The source Bianchi identity and Theorem 12.2 give

DaKab=−κ(Rbfl+ξb+fbcmc+12fcdhbcd).(12.23)D_a\mathscr K^a{}_b=-\kappa\left( \mathscr R_b^{\rm fl}+\xi_b+f_{bc}m^c+\tfrac12f^{cd}h_{bcd}\right). \tag{12.23}

For nonzero interface discrepancy add −κDaΔΣab-\kappa D_a\Delta\Sigma^a{}_b. On the fully coupled Maxwell sector the fluid equation is Rbfl=−ξb\mathscr R_b^{\rm fl}=-\xi_b; its energy and orthogonal force projections are −vbξb-v^b\xi_b and −Πbcξc-\Pi_b{}^c\xi_c. The transfer cancels in the same-realization sum. A law linking charge current to fluid variables is additional information not supplied by these sections. For pure electromagnetic stress as the entire material source, the coupled equations require ξ=0\xi=0, automatically satisfied when J=0J=0.

Stress alone cannot replace the primitive electromagnetic state: Σem(f)=Σem(−f)\Sigma^{\rm em}(f)=\Sigma^{\rm em}(-f), although a potential gauge transformation does not reverse the field. By Lemma 2.1, a gravitational stress readout cannot be used to recover every subsequent electromagnetic operation. The field, current and retained potential/boundary records must remain available.

A complete same-realization Einstein–Maxwell check

Let u=t−zu=t-z, v0=(t+z)/2v_0=(t+z)/2, choose a smooth amplitude A(u)A(u) and a local primitive B′=AB'=A. Specify

φ=B(u)dx,F=A(u)dx∧du,J=0,H=κ2A(u)2(x2+y2).(12.24)\varphi=B(u)\mathrm{d}x,\quad \mathsf F=A(u)\mathrm{d}x\wedge\mathrm{d}u,\quad J=0, \qquad H=\frac\kappa2 A(u)^2(x^2+y^2). \tag{12.24}

The metric and explicit unit-volume calibration are

ds2=−dx2−dy2+2dudv0+Hdu2,Ew=(10000100001−H/2H/200−H/21+H/2).(12.25)\mathrm{d}s^2=-\mathrm{d}x^2-\mathrm{d}y^2+2\mathrm{d}u\mathrm{d}v_0+H\mathrm{d}u^2, \qquad E_w=\begin{pmatrix} 1&0&0&0\\0&1&0&0\\0&0&1-H/2&H/2\\0&0&-H/2&1+H/2 \end{pmatrix}. \tag{12.25}

Its determinant is one and θ4−θ3=du\theta^4-\theta^3=\mathrm{d}u, so the calibrated field has f14=Af_{14}=A and f13=−Af_{13}=-A. In null coordinates the only nonzero connection types are

Ciuu=12∂iH (i=x,y),Cv0ui=Cv0iu=12∂iH,Cv0uu=12∂uH.(12.26)C^i{}_{uu}=\tfrac12\partial_iH\ (i=x,y),\quad C^{v_0}{}_{ui}=C^{v_0}{}_{iu}=\tfrac12\partial_iH,\quad C^{v_0}{}_{uu}=\tfrac12\partial_uH. \tag{12.26}

They give

Ric⁡uu=12(Hxx+Hyy)=κA2,R=0,Tem=A2du⊗du,G=κTem.(12.27)\operatorname{Ric}_{uu}=\tfrac12(H_{xx}+H_{yy})=\kappa A^2, \quad \mathcal R=0,\quad T_{\rm em}=A^2\mathrm{d}u\otimes\mathrm{d}u, \quad G=\kappa T_{\rm em}. \tag{12.27}

All other Ricci components vanish. The field is closed since dA\mathrm{d}A is proportional to du\mathrm{d}u. Its only independent raised null-coordinate component is Fxv0=−A(u)\mathsf F^{xv_0}=-A(u), whose second-index divergence vanishes. Thus the primitive field, its stress and the gravitational equation hold jointly. The transverse curvature equation determines HH from the chosen amplitude; the stress was not defined after computing GG.

For the executed profile A=1+uA=1+u, κ=2\kappa=2, direct evaluation of the chart-tagged gravitational current gives t=0t=0, while T44=T34=(1+t−z)2T^4{}_4=T^3{}_4=(1+t-z)^2. On the unit spatial cube,

Q4=(1+t)3−t33,Φ4out=t2−(1+t)2,Q4′+Φ4out=0.(12.28)Q_4=\frac{(1+t)^3-t^3}{3},\qquad \Phi_4^{\rm out}=t^2-(1+t)^2,\qquad Q_4'+\Phi_4^{\rm out}=0. \tag{12.28}

This is a nonzero finite-flux check, not a claim that gravitational energy universally vanishes. Doubling the electromagnetic field while retaining the geometry still solves source-free Maxwell but quadruples the stress and violates the gravitational equation, explicitly rejecting an incompatible join.

Published source: The Foundations of General Relativity: A Source–Readout Reconstruction of Einstein’s 1916 Theory. Read the original PDF alongside the web treatment.