History of Topics in Special Relativity/Stress-energy tensor (electromagnetic)

Overview
The electromagnetic stress–energy tensor is the contribution to the stress–energy tensor due to the electromagnetic field in terms of the electromagnetic field tensor $$F^{\mu\nu}$$, Poynting vector $$\mathbf{S}$$, Maxwell stress tensor $$\sigma_{ij}$$:


 * $$\begin{matrix}T^{\mu\nu}=\begin{bmatrix}\frac{1}{2}\left(\epsilon_{0}E^{2}+\frac{1}{\mu_{0}}B^{2}\right) & S_{\text{x}}/c & S_{\text{y}}/c & S_{\text{z}}/c\\

S_{\text{x}}/c & -\sigma_{\text{xx}} & -\sigma_{\text{xy}} & -\sigma_{\text{xz}}\\ S_{\text{y}}/c & -\sigma_{\text{yx}} & -\sigma_{\text{yy}} & -\sigma_{\text{yz}}\\ S_{\text{z}}/c & -\sigma_{\text{zx}} & -\sigma_{\text{zy}} & -\sigma_{\text{zz}} \end{bmatrix}=\frac{1}{\mu_{0}}\left[F^{\mu\alpha}F^{\nu}{}_{\alpha}-\frac{1}{4}\eta^{\mu\nu}F_{\alpha\beta}F^{\alpha\beta}\right]\\ \left[\mathbf{S}=\frac{1}{\mu_{0}}\mathbf{E}\times\mathbf{B},\ \sigma_{ij}=\epsilon_{0}E_{i}E_{j}+\frac{1}{\mu_{0}}B_{i}B_{j}-\frac{1}{2}\left(\epsilon_{0}E^{2}+\frac{1}{\mu_{0}}B^{2}\right)\delta_{ij}\right] \end{matrix}$$

This tensor in the form written above is universally accepted when it comes to electromagnetic processes in vacuum. It was introduced by in matrix notation and used in different notation by, ,. Alternative representations have been used by (dyadics) and  (quaternions). However, several competing tensors have been proposed when it comes to electromagnetic processes in physical media, most notably the asymmetric one by and the symmetric one by  (Abraham-Minkowski controversy). Abraham's form was used by and, while Minkowski's was used by , with  summarizing and combining both approaches.

Minkowski (1907/8)
In his lecture from December 1907 (published 1908), Hermann Minkowski used the electromagnetic tensor in vacuum $$f$$ and in matter $$F$$, as well as their duals $$f^{*},F^{*}$$, in order to define the (dielectric) electromagnetic energy-momentum tensor $$S$$ using $$\mathfrak{e},\mathfrak{m}$$ (i.e. $$\mathbf{E},\mathbf{B}$$) and $$\mathfrak{E},\mathfrak{M}$$ (= $$\mathbf{D},\mathbf{H}$$):


 * $$\begin{matrix}S=\left|\begin{array}{cccc}

S_{11}, & S_{12}, & S_{13}, & S_{14}\\ S_{21}, & S_{22}, & S_{23}, & S_{23}\\ S_{31}, & S_{32}, & S_{33}, & S_{34}\\ S_{41}, & S_{42}, & S_{43}, & S_{44} \end{array}\right|=\left|\begin{array}{cccc} X_{x}, & Y_{x}, & Z_{x}, & -iT_{x}\\ X_{y}, & Y_{y}, & Z_{y}, & -iT_{y}\\ X_{z}, & Y_{z}, & Z_{z}, & -iT_{z}\\ -iX_{t}, & -iY_{t}, & -iZ_{t}, & T_{t} \end{array}\right|\\ \hline {\scriptstyle fF=S-L=\left|\begin{array}{llll} S_{11}-L, & S_{12}, & S_{13}, & S_{14}\\ S_{21}, & S_{22}-L, & S_{23}, & S_{23}\\ S_{31}, & S_{32}, & S_{33}-L, & S_{34}\\ S_{41}, & S_{42}, & S_{43}, & S_{44}-L \end{array}\right|,\quad F^{*}f^{*}=-S-L=\left|\begin{array}{llll} -S_{11}-L, & -S_{12}, & -S_{13}, & -S_{14}\\ -S_{21}, & -S_{22}-L, & -S_{23}, & -S_{23}\\ -S_{31}, & -S_{32}, & -S_{33}-L, & -S_{34}\\ -S_{41}, & -S_{42}, & -S_{43}, & -S_{44}-L \end{array}\right|}\\ \hline L=\frac{1}{2}\left(f_{23}F_{23}+f_{31}F_{31}+f_{12}F_{12}+f_{14}F_{14}+f_{24}F_{24}+f_{34}F_{34}\right)\\ S_{11}=\frac{1}{2}\left(f_{23}F_{23}+f_{34}F_{34}+f_{42}F_{42}-f_{12}F_{12}-f_{13}F_{13}-f_{14}F_{14}\right)\\ S_{12}=f_{13}F_{32}+f_{14}F_{42},etc.\\ \hline X_{y}=\mathfrak{m}_{x}\mathfrak{M}_{y}+\mathfrak{e}_{y}\mathfrak{E}_{x},\quad Y_{x}=\mathfrak{m}_{y}\mathfrak{M}_{x}+\mathfrak{e}_{x}\mathfrak{E}_{y},\ {\rm etc}.\\ X_{t}=\mathfrak{e}_{y}\mathfrak{M}_{z}-\mathfrak{e}_{z}\mathfrak{M}_{y},\\ T_{x}=\mathfrak{m}_{z}\mathfrak{E}_{y}-\mathfrak{m}_{y}\mathfrak{E}_{z},\ {\rm etc}.\\ T_{t}=\tfrac{1}{2}(\mathfrak{m}_{x}\mathfrak{M}_{x}+\mathfrak{m}_{y}\mathfrak{M}_{y}+\mathfrak{m}_{z}\mathfrak{M}_{z}+\mathfrak{e}_{x}\mathfrak{E}_{x}+\mathfrak{e}_{y}\mathfrak{E}_{y}+\mathfrak{e}_{z}\mathfrak{E}_{z})\\ L=\tfrac{1}{2}(\mathfrak{m}_{x}\mathfrak{M}_{x}+\mathfrak{m}_{y}\mathfrak{M}_{y}+\mathfrak{m}_{z}\mathfrak{M}_{z}-\mathfrak{e}_{x}\mathfrak{E}_{x}-\mathfrak{e}_{y}\mathfrak{E}_{y}-\mathfrak{e}_{z}\mathfrak{E}_{z}) \end{matrix}$$

where $$X_{x},Y_{x},Z_{x},\dots$$ are the Maxwell stresses, $$T_{x},T_{y},T_{z}$$ the Poynting vector, $$T_{t}$$ the energy density and $$L$$ the Lagrangian. While Minkowski's tensor is symmetrical in vacuum, it is evidently asymmetric in a dielectric because $$\left(X_{t},Y_{t},Z_{t}\right)=\mathbf{E}\times\mathbf{H}$$ while their counterparts have the form $$\left(T_{x},T_{y},T_{z}\right)=\mathbf{B}\times\mathbf{D}$$. Alternatively, symmetrically defined the tensor even in a dielectric, which started the Abraham-Minkowski controversy.

Using that tensor, Minkowski went on to derive the four-force density using four-velocity $$w$$:


 * $$\begin{matrix}K=\text{lor }S\Rightarrow K+(w\overline{K})w\\

\left[\text{lor }=\left|\frac{\partial}{\partial x_{1}},\ \frac{\partial}{\partial x_{2}},\ \frac{\partial}{\partial x_{3}},\ \frac{\partial}{\partial x_{4}}\right|\right] \end{matrix}$$

Abraham (1909-10)
Max Abraham (1909) expressed the stress energy tensor components in terms of arbitrary dielectric constant and permeability. He defined the fields $$\mathfrak{E},\mathfrak{H}$$ (= $$\mathbf{E},\mathbf{H}$$) and $$\mathfrak{D},\mathfrak{B}$$ (i.e. $$\mathbf{D},\mathbf{B}$$) for resting media, and $$\mathfrak{E}',\mathfrak{H}'$$ as the fields acting on moving electric and magnetic unit poles. He started by relating the six components of the spatial part $$X_{x}\dots$$ of the stress-energy tensor to the “relative stresses” $$X_{x}^{\prime}\dots$$ (= the ordinary elastic stresses seen in the moving frame) using momentum density $$\mathfrak{g}$$ and velocity $$\mathfrak{w}$$:


 * $$\begin{matrix}X_{x},\ Y_{y},\ Z_{z},\ X_{y}=Y_{x},\ Y_{z}=Z_{y},\ Z_{x}=X_{z}\\

\hline \left\{ \begin{array}{ccccc} X'_{x}=X_{x}+\mathfrak{w}_{x}\mathfrak{g}_{x}, & & X'_{y}=X_{y}+\mathfrak{w}_{y}\mathfrak{g}_{x}, &  & X'_{z}=X_{z}+\mathfrak{w}_{z}\mathfrak{g}_{x},\\ Y'_{x}=Y_{x}+\mathfrak{w}_{x}\mathfrak{g}_{y}, & & Y'_{y}=Y_{y}+\mathfrak{w}_{y}\mathfrak{g}_{y}, &  & Y'_{z}=Y_{z}+\mathfrak{w}_{z}\mathfrak{g}_{y},\\ Z'_{x}=Z_{x}+\mathfrak{w}_{x}\mathfrak{g}_{z}, & & Z'_{y}=Z_{y}+\mathfrak{w}_{y}\mathfrak{g}_{z}, &  & Z'_{z}=Z_{z}+\mathfrak{w}_{z}\mathfrak{g}_{z}, \end{array}\right.\\ \begin{cases} X'_{x}=\mathfrak{E}'_{x}\mathfrak{D}_{x}+\mathfrak{H}'_{x}\mathfrak{B}_{x}-\frac{1}{2}[\mathfrak{E'D+H'B}\}\\ X'_{y}=\mathfrak{E}'_{x}\mathfrak{D}_{y}+\mathfrak{H}'_{x}\mathfrak{B}_{y},\\ X'_{z}=\mathfrak{E}'_{x}\mathfrak{D}_{z}+\mathfrak{H}'_{x}\mathfrak{B}_{z};\\ Y'_{x}=\mathfrak{E}'_{y}\mathfrak{D}_{x}+\mathfrak{H}'_{y}\mathfrak{B}_{x},\\ Y'_{y}=\mathfrak{E}'_{y}\mathfrak{D}_{y}+\mathfrak{H}'_{y}\mathfrak{B}_{y}-\frac{1}{2}[\mathfrak{E'D+H'B}\},\\ Y'_{z}=\mathfrak{E}'_{y}\mathfrak{D}_{z}+\mathfrak{H}'_{y}\mathfrak{B}_{z};\\ Z'_{x}=\mathfrak{E}'_{z}\mathfrak{D}_{x}+\mathfrak{H}'_{z}\mathfrak{B}_{x},\\ Z'_{y}=\mathfrak{E}'_{z}\mathfrak{D}_{y}+\mathfrak{H}'_{z}\mathfrak{B}_{y},\\ Z'_{z}=\mathfrak{E}'_{z}\mathfrak{D}_{z}+\mathfrak{H}'_{z}\mathfrak{B}_{z}-\frac{1}{2}[\mathfrak{E'D+H'B}\}. \end{cases} \end{matrix}$$

and obtained the other components using energy density $$\psi$$ and Poynting vector $$\mathfrak{S}$$, by which he analyzed all electrodynamic theories of that time, including the following adaptation of Minkowski's model


 * $$\begin{matrix}\begin{cases}

\mathfrak{D}=\epsilon\mathfrak{E}'-[\mathfrak{qH}],\\ \mathfrak{B}=\mu\mathfrak{H}'+[\mathfrak{qE}]; \end{cases}\begin{cases} \mathfrak{E'}=\mathfrak{E}+[\mathfrak{qB}],\\ \mathfrak{H}'=\mathfrak{H}-[\mathfrak{qD}]. \end{cases}\\ \mathfrak{S}'=c[\mathfrak{E}'\mathfrak{H}'],\ \mathfrak{g=}\frac{\mathfrak{S}}{c^{2}}\\ \begin{align}\begin{align}\mathfrak{S} & =c[\mathfrak{E'H'}]+\mathfrak{w}\{\mathfrak{E'D}+\mathfrak{H'B}\}-\mathfrak{D(wE')-B(wH')+w(wg)}\\ c\mathfrak{g} & =\mathfrak{[E'H']+q(E'D)+q(H'B)-D(qE')-B(qH')+q(q}c\mathfrak{g)}\\ & =[\mathfrak{EH}]-\mathfrak{q(qW)}\\ & =[\mathfrak{EH}]\ (\mathfrak{q}=0)\\ \psi & =\frac{1}{2}\mathfrak{ED}+\frac{1}{2}\mathfrak{HB}+\mathfrak{qW} \end{align} \end{align} \\ \left[\begin{align}\mathfrak{W}_{x} & =k^{-2}\left\{ [\mathfrak{DB}]_{x}-[\mathfrak{EH}]_{x}\right\} \\ \mathfrak{W}_{y} & =[\mathfrak{DB}]_{y}-[\mathfrak{EH}]_{y}\\ \mathfrak{W}_{z} & =[\mathfrak{DB}]_{z}-[\mathfrak{EH}]_{z} \end{align} \ \left[k^{2}=1-|\mathfrak{q}|^{2},\ \mathfrak{q}=\frac{\mathfrak{w}}{c}\right]\right]\\ \hline c\mathfrak{g}_{x}=X_{t},\ c\mathfrak{g}_{y}=Y_{t},\ c\mathfrak{g}_{z}=Z_{t}\\ \mathfrak{S}_{x}=cT_{x},\ \mathfrak{S}_{y}=cT_{y},\ \mathfrak{S}_{z}=cT_{z},\ \psi=T_{t}\\ \hline \begin{align}X_{y} & =Y_{x} & Y_{z} & =Z_{y} & Z_{x} & =X_{z}\\ X_{t} & =T_{x} & Y_{t} & =T_{y} & Z_{t} & =T_{z} \end{align} \end{matrix}$$

He noted that the last relations exhibit a "remarkable" symmetry that doesn't occur in Minkowski's original formulation, that is, Abraham's tensor is symmetric because it satisfies Planck's (1907) relation $$\mathfrak{g=}\tfrac{\mathfrak{S}}{c^{2}}$$ even in a dielectric, which is the basis of the Abraham–Minkowski controversy. For instance, the above relations give in the case of a dielectric at rest $$\left(X_{t},Y_{t},Z_{t}\right)=\left(T_{x},T_{y},T_{z}\right)=\mathbf{E}\times\mathbf{H}$$ representing momentum density times the speed of light.

In a subsequent paper (1909) he denoted this system as a "four-dimensional tensor" or "tensor quadruple":

$$\begin{matrix}X_{x}=\frac{1}{2}\left(\mathfrak{E}_{x}^{2}-\mathfrak{E}_{y}^{2}-\mathfrak{E}_{z}^{2}\right)+\frac{1}{2}\left(\mathfrak{H}_{x}^{2}-\mathfrak{H}_{y}^{2}-\mathfrak{H}_{z}^{2}\right)\\ Y_{y}=\frac{1}{2}\left(\mathfrak{E}_{y}^{2}-\mathfrak{E}_{z}^{2}-\mathfrak{E}_{x}^{2}\right)+\frac{1}{2}\left(\mathfrak{H}_{y}^{2}-\mathfrak{H}_{z}^{2}-\mathfrak{H}_{x}^{2}\right)\\ Z_{z}=\frac{1}{2}\left(\mathfrak{E}_{z}^{2}-\mathfrak{E}_{x}^{2}-\mathfrak{E}_{y}^{2}\right)+\frac{1}{2}\left(\mathfrak{H}_{z}^{2}-\mathfrak{H}_{x}^{2}-\mathfrak{H}_{y}^{2}\right)\\ X_{y}=Y_{x}=\mathfrak{E}_{x}\mathfrak{E}_{y}+\mathfrak{H}_{x}\mathfrak{H}_{y}\\ Y_{z}=Z_{y}=\mathfrak{E}_{y}\mathfrak{E}_{z}+\mathfrak{H}_{y}\mathfrak{H}_{z}\\ Z_{x}=X_{z}=\mathfrak{E}_{z}\mathfrak{E}_{x}+\mathfrak{H}_{z}\mathfrak{H}_{x}\\ X_{u}=U_{x}=-i\mathfrak{s}_{x},\ Y_{u}=U_{y}=-i\mathfrak{s}_{y},\ Z_{u}=U_{z}=-i\mathfrak{s}_{z},\ U_{u}=\psi\\ \left[\mathfrak{s}=\mathfrak{S}/c=c\mathfrak{g}=[\mathfrak{EH}],\ \psi=\frac{1}{2}\mathfrak{E}^{2}+\frac{1}{2}\mathfrak{H}^{2},\ u=ict\right]\\ \begin{align}X_{y} & =Y_{x} & Y_{z} & =Z_{y} & Z_{x} & =X_{z}\\ \mathfrak{S}_{x} & =c^{2}\mathfrak{g}_{x} & \mathfrak{S}_{y} & =c^{2}\mathfrak{g}_{y} & \mathfrak{S}_{z} & =c^{2}\mathfrak{g}_{z} \end{align} \end{matrix}$$

and defined the corresponding ponderomotive four-force density. He again referred to the Abraham–Minkowski controversy by alluding to the result of this previous paper, that his tensor remains symmetrical even in the general case of arbitrary dielectric constant and permeability (as opposed to Minkowski's asymmetrical one). He further pointed out another difference to Minkowski's formulation: While Minkowski uses constant rest mass, Abraham showed that mass-energy equivalence requires variable rest mass because of the contribution of heat.

In 1910 he brought Minkowski's matrix formulation into vector form, writing the ten components of the "four-dimensional tensor" in a dielectric as:


 * $$\begin{matrix}\left\{ \begin{array}{ccccccc}

X_{x},\ Y_{y},\ Z_{z}, & & Y_{z}=Z_{y}, &  & Z_{x}=X_{z}, &  & X_{y}=Y_{x};\\ X_{u}=U_{x}, & & Y_{u}=U_{y}, &  & Z_{u}=U_{z}, &  & U_{u}, \end{array}\right.\\ \hline X_{u}=-ic\mathfrak{g}_{x},\ Y_{u}=-ic\mathfrak{g}_{y},\ Z_{u}=-ic\mathfrak{g}_{z}\\ U_{x}=-\frac{i}{c}\mathfrak{S}_{x},\ U_{y}=\frac{-i}{c}\mathfrak{S}_{y},\ U_{z}=\frac{-i}{c}\mathfrak{S}_{z},\ U_{u}=\psi;\\ \hline 2\mathfrak{f}=[\mathfrak{EH}]+[\mathfrak{DB}]-\mathfrak{W}-\mathfrak{q}(\mathfrak{qW})\\ 2\psi=\mathfrak{ED}+\mathfrak{HB}-2(\mathfrak{qW})\\ \left[\mathfrak{W}=(\epsilon\mu-1)\left\{ k^{-2}\mathfrak{f}'+k^{-4}\mathfrak{q}(\mathfrak{qf}')\right\} ,\ \mathfrak{f}=c\mathfrak{g}=\frac{1}{c}\mathfrak{S},\ k=\sqrt{1-\mathfrak{q}^{2}}\right] \end{matrix}$$

again pointing out that they satisfy symmetry conditions which are not satisfied by Minkowski's approach.

Sommerfeld (1910)
In an influential paper, Arnold Sommerfeld translated Minkowski's matrix formalism into a four-dimensional vector formalism involving four-vectors and six-vectors. Sommerfeld defined a tensor $$T=(ff)$$ as the product of six-vector $$f$$ (i.e. electromagnetic field tensor) and its dual $$f^{*}$$ using fields $$\mathfrak{E},\mathfrak{H}$$ (i.e. $$\mathbf{E},\mathbf{H}$$), with $$T$$ having 10 components:


 * $$\begin{matrix}T_{jh}=\left(\mathfrak{B}_{j}\mathfrak{B}_{h}\right)=\mathfrak{B}_{jx}\mathfrak{B}_{hx}+\mathfrak{B}_{jy}\mathfrak{B}_{hy}+\mathfrak{B}_{jz}\mathfrak{B}_{hz}\\

T_{jh}=T_{hj},\ T_{jj}\neq0\\ \hline \left(\mathfrak{B}_{j}\mathfrak{B}_{h}\right)\Rightarrow\left(f_{j}f_{h}\right)\\ \left(f_{j}f_{h}\right)=f_{jx}f_{hx}+f_{jy}f_{hy}+f_{jz}f_{hz}+f_{jl}f_{hl}\\ \left(f_{j}f_{h}\right)=\left(f_{h}f_{j}\right)\\ \hline j=h\Rightarrow\begin{array}{ll} \left(f_{x}f_{x}\right)= & \mathfrak{H}_{y}^{2}+\mathfrak{H}_{z}^{2}-\mathfrak{E}_{x}^{2}\\ \left(f_{y}f_{y}\right)= & \mathfrak{H}_{z}^{2}+\mathfrak{H}_{x}^{2}-\mathfrak{E}_{y}^{2}\\ \left(f_{z}f_{z}\right)= & \mathfrak{H}_{x}^{2}+\mathfrak{H}_{y}^{2}-\mathfrak{E}_{z}^{2}\\ \left(f_{l}f_{l}\right)=- & \mathfrak{E}_{x}^{2}-\mathfrak{E}_{y}^{2}-\mathfrak{E}_{z}^{2} \end{array}\\ j\ne h\Rightarrow\begin{array}{llr} \left(f_{x}f_{y}\right)= & \left(f_{y}f_{x}\right)= & -\left(\mathfrak{H}_{x}\mathfrak{H}_{y}+\mathfrak{E}_{x}\mathfrak{E}_{y}\right)\\ \left(f_{y}f_{z}\right)= & \left(f_{z}f_{y}\right)= & -\left(\mathfrak{H}_{y}\mathfrak{H}_{z}+\mathfrak{E}_{y}\mathfrak{E}_{u}\right)\\ \left(f_{z}f_{x}\right)= & \left(f_{x}f_{z}\right)= & -\left(\mathfrak{H}_{z}\mathfrak{H}_{x}+\mathfrak{E}_{z}\mathfrak{E}_{x}\right)\\ \left(f_{x}f_{l}\right)= & \left(f_{l}f_{x}\right)= & i\left(\mathfrak{E}_{y}\mathfrak{H}_{z}-\mathfrak{E}_{z}\mathfrak{H}_{y}\right)\\ \left(f_{y}f_{l}\right)= & \left(f_{l}f_{y}\right)= & i\left(\mathfrak{E}_{z}\mathfrak{H}_{x}-\mathfrak{E}_{x}\mathfrak{H}_{z}\right)\\ \left(f_{z}f_{l}\right)= & \left(f_{l}f_{z}\right)= & i\left(\mathfrak{E}_{x}\mathfrak{H}_{y}-\mathfrak{E}_{y}\mathfrak{H}_{x}\right) \end{array}\\ \hline \begin{array}{ll} \left(f_{x}^{*}f_{x}^{*}\right)= & -\mathfrak{E}_{y}^{2}-\mathfrak{E}_{z}^{2}+\mathfrak{H}_{x}^{2}\\ \left(f_{x}^{*}f_{y}^{*}\right)= & \left(f_{x}^{*}f_{y}^{*}\right)=+\left(\mathfrak{E}_{x}\mathfrak{E}_{y}+\mathfrak{E}_{x}\mathfrak{E}_{y}\right)\\ \left(f_{x}^{*}f_{l}^{*}\right)= & \left(f_{l}^{*}f_{x}^{*}\right)=i\left(\mathfrak{H}_{y}\mathfrak{E}_{z}-\mathfrak{H}_{z}\mathfrak{E}_{y}\right) \end{array} \end{matrix}$$

which he rewrote following a suggestion by Laue (see next section).

Laue (1910-11)
Max von Laue suggested to Sommerfeld in private conversation the following formulation of the stress energy tensor, published by Sommerfeld in 1910:

$$\begin{matrix}T=\frac{1}{2}\left((f\cdot f)-\left(f^{*}f^{*}\right)\right)\\ \hline \left\{ \begin{array}{ll} T_{xx} & =\frac{1}{2}\left(-\mathfrak{H}_{x}^{2}+\mathfrak{H}_{y}^{2}+\mathfrak{H}_{z}^{2}-\mathfrak{E}_{x}^{2}+\mathfrak{E}_{y}^{2}+\mathfrak{E}_{z}^{2}\right)\\ T_{yy} & =\frac{1}{2}\left(-\mathfrak{H}_{y}^{2}+\mathfrak{H}_{z}^{2}+\mathfrak{H}_{x}^{2}-\mathfrak{E}_{y}^{2}+\mathfrak{E}_{z}^{2}+\mathfrak{E}_{x}^{2}\right)\\ T_{zz} & =\frac{1}{2}\left(-\mathfrak{H}_{z}^{2}+\mathfrak{H}_{x}^{2}+\mathfrak{H}_{y}^{2}-\mathfrak{E}_{z}^{2}+\mathfrak{E}_{x}^{2}+\mathfrak{E}_{y}^{2}\right)\\ T_{ll} & =\frac{1}{2}\left(-\mathfrak{H}_{x}^{2}-\mathfrak{H}_{y}^{2}-\mathfrak{H}_{z}^{2}-\mathfrak{E}_{x}^{2}-\mathfrak{E}_{y}^{2}-\mathfrak{E}_{z}^{2}\right) \end{array}\right.\\ \left\{ \begin{array}{ll} T_{xy}=T_{yx} & =-\mathfrak{H}_{x}\mathfrak{H}_{y}+\mathfrak{E}_{x}\mathfrak{E}_{y}\\ T_{yz}=T_{zx} & =-\mathfrak{H}_{y}\mathfrak{H}_{z}+\mathfrak{E}_{y}\mathfrak{E}_{z}\\ T_{zx}=T_{xz} & =-\mathfrak{H}_{z}\mathfrak{H}_{x}+\mathfrak{E}_{z}\mathfrak{E}_{x} \end{array}\right.\\ \left\{ \begin{array}{ll} T_{xl}=T_{lx} & =i\left(\mathfrak{E}_{y}\mathfrak{H}_{z}-\mathfrak{E}_{z}\mathfrak{H}_{y}\right)\\ T_{yl}=T_{ly} & =i\left(\mathfrak{E}_{z}\mathfrak{H}_{x}-\mathfrak{E}_{x}\mathfrak{H}_{z}\right)\\ T_{zl}=T_{lz} & =i\left(\mathfrak{E}_{x}\mathfrak{H}_{y}-\mathfrak{E}_{y}\mathfrak{H}_{x}\right) \end{array}\right.\\ \hline T_{xx}+T_{yy}+T_{zz}+T_{ll}=0 \end{matrix}$$

In 1911, Laue himself published a paper on the dynamics of relativity, the contents of which he also published at almost the same time in his influential first textbook on relativity. He followed Abraham in defining a symmetric "world tensor" $$T$$ (i.e. electromagnetic stress energy tensor) as the tensor product of the electromagnetic field tensor $$\mathfrak{M}$$, including Poynting vector $$\mathfrak{S}$$, energy density $$W$$, momentum density $$\mathfrak{g}_{e}$$, and its Lorentz transformation:


 * $$\begin{matrix}\begin{matrix}T=\left[[\mathfrak{MM}]\right]=\begin{matrix}\mathbf{p}_{xx} & \mathbf{p}_{xy} & \mathbf{p}_{xz} & \frac{i}{c}\mathfrak{S}_{x}\\

\mathbf{p}_{yx} & \mathbf{p}_{yy} & \mathbf{p}_{yz} & \frac{i}{c}\mathfrak{S}_{y}\\ \mathbf{p}_{zx} & \mathbf{p}_{zy} & \mathbf{p}_{zz} & \frac{i}{c}\mathfrak{S}_{z}\\ \frac{i}{c}\mathfrak{S}_{x} & \frac{i}{c}\mathfrak{S}_{y} & \frac{i}{c}\mathfrak{S}_{z} & -W \end{matrix}\\ \hline \begin{align}T_{xx} & =\frac{1}{2}\left(\mathfrak{E}^{2}+\mathfrak{H}^{2}\right)-\mathfrak{E}_{x}^{2}-\mathfrak{H}_{x}^{2} & & =\mathbf{p}_{xx}\\ T_{xy} & =-\left(\mathfrak{E}_{x}\mathfrak{E}_{y}+\mathfrak{H}_{x}\mathfrak{H}_{x}\right) & & =\mathbf{p}_{xy}\\ T_{ll} & =-\frac{1}{2}\left(\mathfrak{E}^{2}+\mathfrak{H}^{2}\right) & & =-W\\ T_{xl} & =i\left(\mathfrak{E}_{y}\mathfrak{H}_{z}-\mathfrak{E}_{z}\mathfrak{H}_{y}\right) & & =\frac{i}{c}\mathfrak{S}_{x} \end{align} \\ \hline \mathfrak{g}_{e}=\frac{\mathfrak{S}}{c^{2}}=\frac{1}{c}[\mathfrak{EH}] \end{matrix}\\ \hline \\ \begin{matrix}\mathfrak{S}_{x}=\frac{\left(1+\beta^{2}\right)\mathfrak{S}_{x}^{\prime}+v\left(\mathbf{p}_{xx}^{\prime}+W^{\prime}\right)}{1-\beta^{2}},\\ \mathfrak{S}_{y}=\frac{\mathfrak{S}_{y}^{\prime}+v\mathbf{p}_{xy}^{\prime}}{\sqrt{1-\beta^{2}}},\ \mathfrak{S}_{z}=\frac{\mathfrak{S}_{z}^{\prime}+v\mathbf{p}_{xz}^{\prime}}{\sqrt{1-\beta^{2}}}\\ W=\frac{W'+\beta^{2}\mathbf{p}_{xx}^{\prime}+2\frac{v}{c^{2}}\mathfrak{S}_{x}^{\prime}}{1-\beta^{2}}\\ \mathbf{p}_{yy}=\mathbf{p}_{yy}^{\prime},\ \mathbf{p}_{zz}=\mathbf{p}_{zz}^{\prime},\ \mathbf{p}_{yz}=\mathbf{p}_{yz}^{\prime}\\ \mathbf{p}_{xx}=\frac{\mathbf{p}_{xx}^{\prime}+2\frac{v}{c^{2}}\mathfrak{S}_{x}^{\prime}+\beta^{2}W'}{c^{2}-q^{2}},\\ \mathbf{p}_{xy}=\frac{\mathbf{p}_{xy}^{\prime}+\frac{v}{c^{2}}\mathfrak{S}_{y}^{\prime}}{\sqrt{1-\beta^{2}}},\ \mathbf{p}_{xz}=\frac{\mathbf{p}_{xz}^{\prime}+\frac{v}{c^{2}}\mathfrak{S}_{z}^{\prime}}{\sqrt{1-\beta^{2}}}\\ \hline W=\frac{1}{1-\beta^{2}}\left\{ \left(W'+(\mathfrak{vg}')\right)+\frac{1}{c^{2}}\left(\mathfrak{v},\mathfrak{S}'+[\mathfrak{v}\mathbf{p}']\right)\right\} \\ \mathfrak{S}=\frac{1}{1-\beta^{2}}\left\{ \left(\mathfrak{S}'+[\mathfrak{v}\mathbf{p}']\right)+\mathfrak{v}\left(W'+(\mathfrak{vg}')\right)+\frac{1-\sqrt{1-\beta^{2}}}{v^{2}}\left[v,\left[\mathfrak{v},\mathfrak{S}'+[\mathfrak{v}\mathbf{p}']\right]\right]\right\} \end{matrix} \end{matrix}$$

From that he derived, among other things, the electromagnetic energy and momentum of a spherical electron


 * $$\begin{matrix}\mathfrak{G}_{e}=\frac{4E^{0}\mathfrak{q}}{3c\sqrt{c^{2}-q^{2}}},\ E=\frac{c^{2}+\frac{1}{3}q^{2}}{c\sqrt{c^{2}-q^{2}}}E^{0}\\

\left[\mathfrak{G}_{e}=\frac{1}{c^{2}}\int\mathfrak{S}\,dS,\ E=\int W\,dS\right] \end{matrix}$$

and its relation to four-force density $$\mathfrak{F}$$ and four-current $$P$$


 * $$\mathfrak{F}=-\varDelta iv\,T=-[P\mathfrak{M}]$$

Referring to Planck's (1907) "universal relation" relating momentum and energy density which he called "theorem of inertia of energy", he reproduced Abraham's expressions for energy and momentum in a dielectric at rest as well as in motion:


 * $$\begin{matrix}\mathfrak{F}^{0}=-\mathfrak{div}\mathbf{p}^{0}-\frac{1}{c}\frac{\partial}{\partial t}\left[\mathfrak{E}^{0}\mathfrak{H}^{0}\right]\\

\mathfrak{S}^{0}=c\left[\mathfrak{E}^{0}\mathfrak{H}^{0}\right],\ \mathfrak{g}_{e}^{0}=\frac{1}{c}\left[\mathfrak{E}^{0}\mathfrak{H}^{0}\right]\Rightarrow\mathfrak{g}_{e}^{0}=\frac{\mathfrak{S}^{0}}{c^{2}}\\ \hline \mathfrak{S}=c^{2}\mathfrak{g}_{e}=c[\mathfrak{EH}]+\mathfrak{q}\frac{(\mathfrak{q},[\mathfrak{EH}]-[\mathfrak{DB}])}{c^{2}-q^{2}}\\ W=\frac{1}{2}\left((\mathfrak{ED})+(\mathfrak{HB})\right)+\frac{c}{c^{2}-q^{2}}(\mathfrak{q},[\mathfrak{EH}]-[\mathfrak{DB}]) \end{matrix}$$

and wrote the relations also in terms of Abraham's "relative stresses" (i.e. the elastic stresses) as follows (Laue not only referred these equations to electrodynamics but to mechanics as well):


 * $$\begin{matrix}\mathbf{t}_{jk}=\mathbf{p}_{jk}-\mathfrak{g}_{j}\mathfrak{q}_{k}\\

c^{2}\mathfrak{g}=\mathfrak{S}=W\mathfrak{q}+[\mathfrak{q}\mathbf{t}]-\left[\mathfrak{q}[\mathfrak{qg}]\right]=W\mathfrak{q}+[\mathfrak{q}\mathbf{t}]-\mathfrak{q}[\mathfrak{qg}]+\mathfrak{g}q^{2} \end{matrix}$$

Ishiwara (1910-11)
Jun Ishiwara (1910) chose Minkowski's asymmetric formulation of the stress-energy tensor in a dielectric over Abraham's symmetric one, and also gave its Lorentz transformation:


 * $$\begin{matrix}\begin{align}X_{x} & =\frac{1}{4\pi}\left\{ \mathfrak{E}_{x}\mathfrak{D}_{x}+\mathfrak{H}_{x}\mathfrak{B}_{x}-\frac{1}{2}(\mathfrak{ED}+\mathfrak{HB})\right\} \\

X_{y} & =\frac{1}{4\pi}\left\{ \mathfrak{E}_{x}\mathfrak{D}_{y}+\mathfrak{H}_{x}\mathfrak{B}_{y}\right\} \\ X_{z} & =\frac{1}{4\pi}\left\{ \mathfrak{E}_{x}\mathfrak{D}_{z}+\mathfrak{H}_{x}\mathfrak{B}_{z}\right\} \\ Y_{x} & =\frac{1}{4\pi}\left\{ \mathfrak{E}_{y}\mathfrak{D}_{x}+\mathfrak{H}_{y}\mathfrak{B}_{x}\right\} \\ Y_{y} & =\frac{1}{4\pi}\left\{ \mathfrak{E}_{y}\mathfrak{D}_{y}+\mathfrak{H}_{y}\mathfrak{B}_{y}-\frac{1}{2}(\mathfrak{ED}+\mathfrak{HB})\right\} \\ Y_{z} & =\frac{1}{4\pi}\left\{ \mathfrak{E}_{y}\mathfrak{D}_{z}+\mathfrak{H}_{y}\mathfrak{B}_{z}\right\} \\ Z_{x} & =\frac{1}{4\pi}\left\{ \mathfrak{E}_{z}\mathfrak{D}_{x}+\mathfrak{H}_{z}\mathfrak{B}_{x}\right\} \\ Z_{y} & =\frac{1}{4\pi}\left\{ \mathfrak{E}_{z}\mathfrak{D}_{y}+\mathfrak{H}_{z}\mathfrak{B}_{y}\right\} \\ Z_{z} & =\frac{1}{4\pi}\left\{ \mathfrak{E}_{z}\mathfrak{D}_{z}+\mathfrak{H}_{z}\mathfrak{B}_{z}-\frac{1}{2}(\mathfrak{ED}+\mathfrak{HB})\right\} \\ \mathfrak{g} & =\frac{1}{4\pi c}[\mathfrak{DB}]\\ \mathfrak{S} & =\frac{c}{4\pi}[\mathfrak{EH}]\\ u & =\frac{1}{8\pi}\left\{ \mathfrak{E}\mathfrak{D}+\mathfrak{H}\mathfrak{B}\right\} \end{align} & &  & \left|\begin{align}\varkappa^{2}\frac{\mathfrak{S}_{x}^{\prime}}{c} & =\frac{\mathfrak{S}_{x}}{c}+\beta X_{x}-\beta u+\beta^{2}c\mathfrak{g}_{x}\\ \varkappa\frac{\mathfrak{S}_{y}^{\prime}}{c} & =\frac{\mathfrak{S}_{y}}{c}+\beta X_{y}\\ \varkappa\frac{\mathfrak{S}_{z}^{\prime}}{c} & =\frac{\mathfrak{S}_{z}}{c}+\beta X_{z}\\ \varkappa^{2}u^{\prime} & =u-\beta\frac{\mathfrak{S}_{x}}{c}-\beta c\mathfrak{g}_{x}-\beta^{2}X_{x}\\ \varkappa^{2}X_{x}^{\prime} & =X_{x}+\beta\frac{\mathfrak{S}_{x}}{c}+\beta c\mathfrak{g}_{x}-\beta^{2}u\\ \varkappa X_{y}^{\prime} & =X_{y}+\beta\frac{\mathfrak{S}_{y}}{c}\\ \varkappa X_{z}^{\prime} & =X_{z}+\beta\frac{\mathfrak{S}_{z}}{c}\\ \varkappa Y_{x}^{\prime} & =Y_{x}+\beta c\mathfrak{g}_{y}\\ Y_{y}^{\prime} & =Y_{y}\\ Y_{z}^{\prime} & =Y_{z}\\ \varkappa Z_{x}^{\prime} & =Z_{x}+\beta c\mathfrak{g}_{z}\\ Z_{y}^{\prime} & =Z_{y}\\ Z_{z}^{\prime} & =Z_{z}\\ \varkappa^{2}c\mathfrak{g}_{x}^{\prime} & =c\mathfrak{g}_{x}-\beta u+\beta X_{x}+\beta^{2}\frac{\mathfrak{S}_{x}}{c}\\ \varkappa c\mathfrak{g}_{y}^{\prime} & =c\mathfrak{g}_{y}+\beta Y_{x}\\ \varkappa c\mathfrak{g}_{z}^{\prime} & =c\mathfrak{g}_{z}+\beta Z_{x} \end{align} \right.\\ & &  & \left[\varkappa=\sqrt{1-\frac{\mathfrak{v}^{2}}{c^{2}}}\right] \end{matrix}$$

However, in order to derive the ponderomotive force from this tensor, he chose Abraham's approach of variable rest mass over Minkowski's assumption of constant rest mass.

In 1911 he wrote this tensor and its Lorentz transformation as follows (the indices $$\mathfrak{v}$$ and $$\mathfrak{\bar{v}}$$ indicate the components parallel and perpendicular to the velocity, respectively):

$$\begin{matrix}\begin{matrix}\begin{align}\mathbf{T}_{11} & =\mathfrak{T}_{xx} & \mathbf{T}_{12} & =\mathfrak{T}_{xy} & \mathbf{T}_{13} & =\mathfrak{T}_{xz} & \mathbf{T}_{14} & =-ic\mathfrak{g}_{x}\\ \mathbf{T}_{21} & =\mathfrak{T}_{yx} & \mathbf{T}_{22} & =\mathfrak{T}_{yy} & \mathbf{T}_{23} & =\mathfrak{T}_{yz} & \mathbf{T}_{24} & =-ic\mathfrak{g}_{y}\\ \mathbf{T}_{31} & =\mathfrak{T}_{zx} & \mathbf{T}_{32} & =\mathfrak{T}_{zy} & \mathbf{T}_{33} & =\mathfrak{T}_{zz} & \mathbf{T}_{34} & =-ic\mathfrak{g}_{z}\\ \mathbf{T}_{41} & =-\frac{i}{c}\mathfrak{S}_{x} & \mathbf{T}_{42} & =-\frac{i}{c}\mathfrak{S}_{y} & \mathbf{T}_{43} & =-\frac{i}{c}\mathfrak{S}_{z} & \mathbf{T}_{44} & =u \end{align} \\ \hline \begin{align}\mathfrak{T}_{jk} & =\frac{1}{4\pi}\left(\mathfrak{E}_{j}\mathfrak{D}_{k}+\mathfrak{H}_{j}\mathfrak{B}_{k}\right)-\frac{(\mathfrak{jk})}{8\pi}\left\{ (\mathfrak{E}\mathfrak{D})+(\mathfrak{H}\mathfrak{B})\right\} \\ \mathfrak{g} & =\frac{1}{4\pi c}[\mathfrak{DB}]\\ \mathfrak{S} & =\frac{c}{4\pi}[\mathfrak{EH}]\\ u & =\frac{1}{8\pi}\left\{ (\mathfrak{E}\mathfrak{D})+(\mathfrak{H}\mathfrak{B})\right\} \\ & \mathfrak{T}_{xx}+\mathfrak{T}_{yy}+\mathfrak{T}_{zz}+u=0 \end{align} \end{matrix}\\ \hline \\ \begin{align}\mathfrak{T_{vv}^{\prime}} & =\frac{1}{1-\frac{\mathfrak{v}^{2}}{c^{2}}}\left(\mathfrak{T_{vv}}+\mathfrak{v}\frac{\mathfrak{S_{v}}}{c^{2}}+\mathfrak{v\mathfrak{g}_{v}}-\frac{\mathfrak{v}^{2}}{c^{2}}u\right)\\ \mathfrak{T_{\bar{v}\bar{v}}^{\prime}} & =\mathfrak{T_{\bar{v}\bar{v}}}\\ \mathfrak{T_{v\bar{v}}^{\prime}} & =\frac{1}{\sqrt{1-\frac{\mathfrak{v}^{2}}{c^{2}}}}\left(\mathfrak{T_{v\bar{v}}}+\mathfrak{v}\frac{\mathfrak{S_{\bar{v}}}}{c^{2}}\right)\\ \mathfrak{T_{\bar{v}v}^{\prime}} & =\frac{1}{\sqrt{1-\frac{\mathfrak{v}^{2}}{c^{2}}}}\left(\mathfrak{T_{\bar{v}v}}+\mathfrak{vg_{\bar{v}}}\right)\\ \frac{\mathfrak{S_{v}^{\prime}}}{c} & =\frac{1}{c\left(1-\frac{\mathfrak{v}^{2}}{c^{2}}\right)}\left(\mathfrak{S_{v}}+\mathfrak{v}\mathfrak{T_{vv}}-\mathfrak{v}u-\mathfrak{v^{2}g_{v}}\right)\\ \frac{\mathfrak{S_{\bar{v}}^{\prime}}}{c} & =\frac{1}{c\sqrt{1-\frac{\mathfrak{v}^{2}}{c^{2}}}}\left(\mathfrak{S_{\bar{v}}}+\mathfrak{v}\mathfrak{T_{v\bar{v}}}\right)\\ c\mathfrak{g_{v}^{\prime}} & =\frac{1}{c\left(1-\frac{\mathfrak{v}^{2}}{c^{2}}\right)}\left(c^{2}\mathfrak{g_{v}}-\mathfrak{v}u+\mathfrak{v}\mathfrak{T_{vv}}+\mathfrak{v}^{2}\frac{\mathfrak{S_{v}}}{c^{2}}\right)\\ c\mathfrak{g_{\bar{v}}^{\prime}} & =\frac{1}{c\sqrt{1-\frac{\mathfrak{v}^{2}}{c^{2}}}}\left(c^{2}\mathfrak{g_{\bar{v}}}+\mathfrak{v}\mathfrak{T_{\bar{v}v}}\right)\\ u' & =\frac{1}{1-\frac{\mathfrak{v}^{2}}{c^{2}}}\left(u-\mathfrak{v}\frac{\mathfrak{S_{v}}}{c^{2}}-\mathfrak{v\mathfrak{g}_{v}}-\mathfrak{v}^{2}\frac{\mathfrak{\mathfrak{T_{vv}}}}{c^{2}}\right) \end{align} \end{matrix}$$

Einstein (1912-14)
In an unpublished manuscript on special relativity, written between 1912-1914, Albert Einstein followed Abraham and Laue in defining the symmetric stress energy tensors in vacuum (upper index 0) and matter (upper index $$e$$ electric polarization, $$m$$ magnetic polarization), in which the upper asterisk represent the dual, $$\mathfrak{g}$$ momentum density, $$\mathbf{s}$$ energy flux:


 * $$\begin{matrix}T_{\mu\nu}=\begin{matrix}p_{xx} & p_{xy} & p_{xz} & ic\mathfrak{g}_{x}\\

p_{yx} & p_{yy} & p_{yz} & ic\mathfrak{g}_{y}\\ p_{zx} & p_{zy} & p_{zz} & ic\mathfrak{g}_{z}\\ \frac{i}{c}\mathbf{s}_{x} & \frac{i}{c}\mathbf{s}_{y} & \frac{i}{c}\mathbf{s}_{z} & -\eta \end{matrix}\\ \left[p_{xy}=p_{yx}{\rm etc.,\ \mathfrak{g}=\frac{1}{c^{2}}\mathbf{s}}\right]\\ \hline \left(T_{\mu\nu}\right)=\left(T_{\mu\nu}^{0}\right)+\left(T_{\mu\nu}^{e}\right)+\left(T_{\mu\nu}^{m}\right)\\ \hline \begin{align}\left(T_{\mu\nu}^{0}\right) & =\frac{1}{2}\left\{ \left(\mathfrak{F}_{\mu\sigma}\right)\left(\mathfrak{F}_{\nu\sigma}\right)-\left(\mathfrak{F}_{\mu\sigma}^{\ast}\right)\left(\mathfrak{F}_{\nu\sigma}^{\ast}\right)\right\} \\ \left(T_{\mu\nu}^{e}\right) & =\frac{1}{\varepsilon-1}\left\{ \left(\mathfrak{P}_{\mu\sigma}\right)\left(\mathfrak{P}_{\nu\sigma}\right)-\frac{1}{4}\left(\delta_{\mu\nu}\right)\left(\mathfrak{P}_{\sigma\tau}\right)\left(\mathfrak{P}_{\sigma\tau}\right)\right\} \\ \left(T_{\mu\nu}^{m}\right) & =\frac{1}{\mu-1}\left\{ \left(\mathfrak{M}_{\mu\sigma}\right)\left(\mathfrak{M}_{\nu\sigma}\right)-\frac{1}{4}\left(\delta_{\mu\nu}\right)\left(\mathfrak{M}_{\sigma\tau}\right)\left(\mathfrak{M}_{\sigma\tau}\right)\right\} \end{align} \\ \hline \begin{align}p_{xx}^{0} & =-\mathfrak{e}_{x}^{2}-\mathfrak{h}_{x}^{2}+\frac{1}{2}\left(\mathfrak{e}^{2}+\mathfrak{h}^{2}\right) & p_{xx}^{e} & =\frac{1}{\varepsilon-1}\left\{ -\mathfrak{p}_{x}^{\ast2}-\mathfrak{p}_{x}^{\ast\ast2}+\frac{1}{2}\left(\mathfrak{p}^{\ast2}+\mathfrak{p}^{\ast\ast2}\right)\right\} \\ p_{xy}^{0} & =-\mathfrak{e}_{x}\mathfrak{e}_{y}-\mathfrak{h}_{x}\mathfrak{h}_{y} & p_{xy}^{e} & =\frac{1}{\varepsilon-1}\left\{ -\mathfrak{p}_{x}^{\ast}\mathfrak{p}_{y}^{\ast}-\mathfrak{p}_{x}^{\ast\ast}\mathfrak{p}_{y}^{\ast\ast}\right\} \\ & {\rm etc.} & & {\rm etc.}\\ \mathbf{s}^{0} & =c^{2}\mathfrak{g}^{0}=c[\mathfrak{e,h}] & \mathbf{s}^{e} & =c^{2}\mathfrak{g}^{e}=\frac{c}{\varepsilon-1}[\mathfrak{p^{\ast\ast},p}^{\ast}]\\ \eta^{0} & =\frac{\mathfrak{e}^{2}+\mathfrak{h}^{2}}{2} & \eta^{e} & =\frac{1}{\varepsilon-1}\frac{\mathfrak{p}^{\ast2}+\mathfrak{p}^{\ast\ast2}}{2} \end{align} \\ \hline T_{\mu\nu}^{\prime}=\sum_{\sigma\tau}\alpha_{\mu\sigma}\alpha_{\nu\tau}T_{\sigma\tau} \end{matrix}$$

Silberstein (1912-14)
Ludwik Silberstein used Quaternions to give an alternative representation of the components of the electromagnetic stress-energy tensor. He derived the Poynting vector $$\mathfrak{P}$$, energy density $$u$$ and the Maxwell stresses $$f_{n}$$ from the electromagnetic bivector (= Weber vector) $$\mathbf{G}$$ and its dual $$\mathbf{F}$$, and the Lorentz transformation of all these quantities:


 * $$\begin{matrix}\begin{align}\mathfrak{P} & =-\frac{\iota c}{2}V\mathbf{GF}=cV\mathbf{EM}\\

u & =\frac{1}{2}(\mathbf{GF})=\frac{1}{2}\left(E^{2}+M^{2}\right)\\ f_{n} & =f\mathbf{n}=\frac{1}{2}(\mathbf{GF})\mathbf{n}-\frac{1}{2}\mathbf{F}(\mathbf{Gn})-\frac{1}{2}\mathbf{G}(\mathbf{Fn})\\ & =\frac{1}{2}V\mathbf{G}n\mathbf{F} \end{align} \\ \left[\frac{1}{2}\mathbf{GF}=-u+\frac{\iota}{c}\mathfrak{P}\right]\\ \hline \begin{align}\frac{1}{\gamma^{2}}u & =u'+\frac{2}{c^{2}}(\mathfrak{P}'\mathbf{v})+\frac{1}{c^{2}}(\mathbf{v}f'\mathbf{v})\\ \frac{1}{\gamma^{2}}\epsilon\mathfrak{P} & =\mathfrak{P}'+\left[\frac{1}{c^{2}}(\mathbf{v}\mathfrak{P}')+u'+f'\right]\mathbf{v}\\ \frac{1}{\gamma^{2}}\epsilon f & =f'\frac{1}{\epsilon}+\frac{1}{c}\left[\mathfrak{P}'+u'\mathbf{v}\right]\left(\mathbf{v}+\frac{1}{c^{2}}\mathbf{v}\left(\frac{1}{\epsilon}\mathfrak{P}'\right.\right. \end{align} \end{matrix}$$

Silberstein used equivalent expressions also in his textbook on quaternionic special relativity in 1914, in which he also discussed Minkowski's electrodynamics of media by giving the electromagnetic bivectors of an isotropic medium:


 * $$\begin{matrix}\begin{align}u & =\frac{1}{2}\left(\mathbf{E}\mathfrak{E}+\mathbf{M}\mathfrak{M}\right)=\\

\mathfrak{P} & =cV\mathbf{EM}\\ \mathbf{g} & =\frac{1}{c}V\mathfrak{EM}\\ \mathbf{f}_{n} & =u\mathbf{n}-\mathbf{E}(\mathfrak{E}\mathbf{n})-\mathbf{M}(\mathfrak{M}\mathbf{n}) \end{align} \end{matrix}$$

Lewis/Wilson (1912)
Gilbert Newton Lewis and Edwin Bidwell Wilson devised an alternative 4D vector calculus based on Dyadics which, however, never gained widespread support. They gave the dyadic $$\Psi$$ (i.e. the stress-energy tensor) using unit dyadic $$\mathbf{I}$$ in terms of the electromagnetic 2-vector $$\mathbf{M}$$ (i.e. electromagnetic tensor) and its dual $$\mathbf{M}^{\ast}$$:


 * $$\begin{matrix}\Psi=\frac{1}{2}\left(\Phi+\Phi^{\ast}\right)\\

X_{x}\mathbf{k}_{1}\mathbf{k}_{1}+X_{y}\mathbf{k}_{1}\mathbf{k}_{2}+X_{z}\mathbf{k}_{1}\mathbf{k}_{3}-X_{t}\mathbf{k}_{1}\mathbf{k}_{4}\\ +Y_{x}\mathbf{k}_{2}\mathbf{k}_{1}+Y_{y}\mathbf{k}_{2}\mathbf{k}_{2}+Y_{z}\mathbf{k}_{2}\mathbf{k}_{3}-Y_{t}\mathbf{k}_{2}\mathbf{k}_{4}\\ +Z_{x}\mathbf{k}_{3}\mathbf{k}_{1}+Z_{y}\mathbf{k}_{3}\mathbf{k}_{2}+Z_{z}\mathbf{k}_{3}\mathbf{k}_{3}-Z_{t}\mathbf{k}_{3}\mathbf{k}_{4}\\ -T_{x}\mathbf{k}_{4}\mathbf{k}_{1}-T_{y}\mathbf{k}_{4}\mathbf{k}_{2}-T_{z}\mathbf{k}_{4}\mathbf{k}_{3}-T_{t}\mathbf{k}_{4}\mathbf{k}_{4}\\ \hline \Phi=(\mathbf{I}\cdot\mathbf{M})\cdot(\mathbf{I}\cdot\mathbf{M}),\ \Phi^{\ast}=(\mathbf{I}\cdot\mathbf{M}^{\ast})\cdot(\mathbf{I}\cdot\mathbf{M}^{\ast})\\ \hline \begin{align}\frac{1}{2}\left(\mathbf{k}_{1}\cdot\mathbf{M}\right)\cdot\mathbf{M}+\frac{1}{2}\left(\mathbf{k}_{1}\cdot\mathbf{M}^{\ast}\right)\cdot\mathbf{M}^{\ast} & =X_{x}\mathbf{k}_{1}+X_{y}\mathbf{k}_{2}+X_{z}\mathbf{k}_{3}-X_{t}\mathbf{k}_{4}\\ \frac{1}{2}\left(\mathbf{k}_{2}\cdot\mathbf{M}\right)\cdot\mathbf{M}+\frac{1}{2}\left(\mathbf{k}_{2}\cdot\mathbf{M}^{\ast}\right)\cdot\mathbf{M}^{\ast} & =Y_{x}\mathbf{k}_{1}+Y_{y}\mathbf{k}_{2}+Y_{z}\mathbf{k}_{3}-Y_{t}\mathbf{k}_{4}\\ \frac{1}{2}\left(\mathbf{k}_{3}\cdot\mathbf{M}\right)\cdot\mathbf{M}+\frac{1}{2}\left(\mathbf{k}_{3}\cdot\mathbf{M}^{\ast}\right)\cdot\mathbf{M}^{\ast} & =Z_{x}\mathbf{k}_{1}+Z_{y}\mathbf{k}_{2}+Z_{z}\mathbf{k}_{3}-Z_{t}\mathbf{k}_{4}\\ \frac{1}{2}\left(\mathbf{k}_{4}\cdot\mathbf{M}\right)\cdot\mathbf{M}+\frac{1}{2}\left(\mathbf{k}_{4}\cdot\mathbf{M}^{\ast}\right)\cdot\mathbf{M}^{\ast} & =T_{x}\mathbf{k}_{1}+T_{y}\mathbf{k}_{2}+T_{z}\mathbf{k}_{3}+T_{t}\mathbf{k}_{4} \end{align} \\ \hline \begin{align}X_{x} & =\frac{1}{2}\left(e_{1}^{2}-e_{2}^{2}-e_{3}^{2}+h_{1}^{2}-h_{2}^{2}-h_{3}^{2}\right)\\ Y_{y} & =\frac{1}{2}\left(e_{2}^{2}-e_{3}^{2}-e_{1}^{2}+h_{2}^{2}-h_{3}^{2}-h_{1}^{2}\right)\\ Z_{z} & =\frac{1}{2}\left(e_{3}^{2}-e_{1}^{2}-e_{2}^{2}+h_{3}^{2}-h_{1}^{2}-h_{2}^{2}\right)\\ T_{t} & =\frac{1}{2}\left(e_{1}^{2}+e_{2}^{2}+e_{3}^{2}+h_{1}^{2}+h_{2}^{2}+h_{3}^{2}\right)\\ X_{y} & =Y_{x}=e_{1}e_{2}+h_{1}h_{2}\ {\rm etc.}\\ T_{x} & =X_{t}=e_{2}h_{3}-e_{3}h_{2}\ {\rm etc.} \end{align} \end{matrix}$$

Henschke (1912/13)
Erich Henschke (in a thesis from 1912, published 1913) derived the stress-energy tensor from a variational principle, obtaining different formulations of four-force density on the basis of a symmetric tensor equivalent to the one of Abraham and Laue (also using their notation). He started with the case of vacuum:

$$\begin{matrix}\left\{ \begin{align}X_{x} & =\frac{1}{2}\left(\mathfrak{H}_{x}^{2}-\mathfrak{H}_{y}^{2}-\mathfrak{H}_{z}^{2}+\mathfrak{E}_{x}^{2}-\mathfrak{E}_{y}^{2}-\mathfrak{E}_{z}^{2}\right)\\ Y_{z} & =\frac{1}{2}\left(-\mathfrak{H}_{x}^{2}+\mathfrak{H}_{y}^{2}-\mathfrak{H}_{z}^{2}-\mathfrak{E}_{x}^{2}+\mathfrak{E}_{y}^{2}-\mathfrak{E}_{z}^{2}\right)\\ Z_{z} & =\frac{1}{2}\left(-\mathfrak{H}_{x}^{2}-\mathfrak{H}_{y}^{2}+\mathfrak{H}_{z}^{2}+\mathfrak{E}_{x}^{2}+\mathfrak{E}_{y}^{2}+\mathfrak{E}_{z}^{2}\right)\\ U_{u} & =\frac{1}{2}\left(\mathfrak{H}_{x}^{2}+\mathfrak{H}_{y}^{2}+\mathfrak{H}_{z}^{2}+\mathfrak{E}_{x}^{2}+\mathfrak{E}_{y}^{2}+\mathfrak{E}_{z}^{2}\right) \end{align} \right.\\ \left\{ \begin{align}X_{y}=Y_{x} & =\mathfrak{H}_{x}\mathfrak{H}_{y}+\mathfrak{E}_{x}\mathfrak{E}_{y}\\ X_{z}=Z_{x} & =\mathfrak{H}_{x}\mathfrak{H}_{z}+\mathfrak{E}_{x}\mathfrak{E}_{z}\\ Y_{z}=Z_{y} & =\mathfrak{H}_{y}\mathfrak{H}_{z}+\mathfrak{E}_{y}\mathfrak{E}_{z} \end{align} \right.\\ \left\{ \begin{align}U_{x}=X_{u} & =-i\left(\mathfrak{E}_{y}\mathfrak{H}_{x}-\mathfrak{E}_{z}\mathfrak{H}_{y}\right)\\ U_{y}=Y_{u} & =-i\left(\mathfrak{E}_{z}\mathfrak{H}_{x}-\mathfrak{E}_{z}\mathfrak{H}_{z}\right)\\ U_{x}=Z_{u} & =-i\left(\mathfrak{E}_{x}\mathfrak{H}_{y}-\mathfrak{E}_{y}\mathfrak{H}_{x}\right) \end{align} \right.\\ \hline T=\left[\left[\mathfrak{MM}\right]\right],\ \mathfrak{F}=-\varDelta iv\,T \end{matrix}$$

then in moving media


 * $$\begin{matrix}\left\{ \begin{align}X_{x} & =\frac{1}{2}\left(\mathfrak{E}_{x}\mathfrak{D}_{x}-\mathfrak{E}_{y}\mathfrak{D}_{y}-\mathfrak{E}_{z}\mathfrak{D}_{z}+\mathfrak{H}_{x}\mathfrak{B}_{x}-\mathfrak{H}_{y}\mathfrak{B}_{y}-\mathfrak{H}_{z}\mathfrak{B}_{z}+\frac{\mathfrak{v}_{x}}{c}\mathfrak{W}_{x}\right)\\

Y_{z} & =\frac{1}{2}\left(-\mathfrak{E}_{x}\mathfrak{D}_{x}+\mathfrak{E}_{y}\mathfrak{D}_{y}-\mathfrak{E}_{z}\mathfrak{D}_{z}-\mathfrak{H}_{x}\mathfrak{B}_{x}+\mathfrak{H}_{y}\mathfrak{B}_{y}-\mathfrak{H}_{z}\mathfrak{B}_{z}+\frac{\mathfrak{v}_{y}}{c}\mathfrak{W}_{y}\right)\\ Z_{z} & =\frac{1}{2}\left(-\mathfrak{E}_{x}\mathfrak{D}_{x}-\mathfrak{E}_{y}\mathfrak{D}_{y}+\mathfrak{E}_{z}\mathfrak{D}_{z}-\mathfrak{H}_{x}\mathfrak{B}_{x}-\mathfrak{H}_{y}\mathfrak{B}_{y}+\mathfrak{H}_{z}\mathfrak{B}_{z}+\frac{\mathfrak{v}_{z}}{c}\mathfrak{W}_{z}\right)\\ U_{u} & =\frac{1}{2}\left(\mathfrak{E}_{x}\mathfrak{D}_{x}+\mathfrak{E}_{y}\mathfrak{D}_{y}+\mathfrak{E}_{z}\mathfrak{D}_{z}+\mathfrak{H}_{x}\mathfrak{B}_{x}+\mathfrak{H}_{y}\mathfrak{B}_{y}+\mathfrak{H}_{z}\mathfrak{B}_{z}+\frac{\mathfrak{v}}{c}\mathfrak{W}\right) \end{align} \right.\\ \left\{ \begin{align}X_{y}=Y_{x} & =\frac{1}{2}\left(\mathfrak{E}_{x}\mathfrak{D}_{y}+\mathfrak{E}_{y}\mathfrak{D}_{x}+\mathfrak{H}_{x}\mathfrak{B}_{y}+\mathfrak{H}_{y}\mathfrak{B}_{x}+\frac{\mathfrak{v}_{x}}{c}\mathfrak{W}_{y}+\frac{\mathfrak{v}_{y}}{c}\mathfrak{W}_{x}\right)\\ X_{z}=Z_{x} & =\frac{1}{2}\left(\mathfrak{E}_{x}\mathfrak{D}_{z}+\mathfrak{E}_{z}\mathfrak{D}_{x}+\mathfrak{H}_{x}\mathfrak{B}_{z}+\mathfrak{H}_{z}\mathfrak{B}_{x}+\frac{\mathfrak{v}_{x}}{c}\mathfrak{W}_{z}+\frac{\mathfrak{v}_{z}}{c}\mathfrak{W}_{x}\right)\\ Y_{z}=Z_{y} & =\frac{1}{2}\left(\mathfrak{E}_{y}\mathfrak{D}_{z}+\mathfrak{E}_{z}\mathfrak{D}_{y}+\mathfrak{H}_{y}\mathfrak{B}_{z}+\mathfrak{H}_{z}\mathfrak{B}_{y}+\frac{\mathfrak{v}_{y}}{c}\mathfrak{W}_{z}+\frac{\mathfrak{v}_{z}}{c}\mathfrak{W}_{y}\right) \end{align} \right.\\ \left\{ \begin{align}U_{x}=X_{u} & =-i\frac{1}{2}\left([\mathfrak{EH}]_{x}+[\mathfrak{DB}]_{x}-\mathfrak{W}_{x}-\frac{\mathfrak{v}_{x}}{c}\left(\frac{\mathfrak{v}}{c}\mathfrak{W}\right)\right)\\ U_{y}=Y_{u} & =-i\frac{1}{2}\left([\mathfrak{EH}]_{y}+[\mathfrak{DB}]_{y}-\mathfrak{W}_{y}-\frac{\mathfrak{v}_{y}}{c}\left(\frac{\mathfrak{v}}{c}\mathfrak{W}\right)\right)\\ U_{x}=Z_{u} & =-i\frac{1}{2}\left([\mathfrak{EH}]_{z}+[\mathfrak{DB}]_{z}-\mathfrak{W}_{z}-\frac{\mathfrak{v}_{z}}{c}\left(\frac{\mathfrak{v}}{c}\mathfrak{W}\right)\right) \end{align} \right.\\ \hline \mathfrak{S}_{x}=c[\mathfrak{EH}]_{x}\ (\mathfrak{v}=0) \end{matrix}$$

Grammel (1913)
Richard Grammel developed three stress-energy tensors $$T', T, T'$$ in a dielectric satisfying the general four-force equation (using Laue's notation) $$F'=-\varDelta iv\,T',\dots$$. The first one is identical to Abraham's, the second one a modification of Minkowski's so that it becomes symmetrical in the case of rest, and the third one is related to the other ones by $$F'=2F'-F$$. He started with the case of rest


 * $$\begin{align}E & =\frac{1}{2}\left\{ (\mathfrak{ED})+(\mathfrak{HB})\right\} \\

\mathfrak{S} & =\mathfrak{g}=[\mathfrak{EH}]\\ \overline{X_{x}^{\prime}} & =\mathfrak{E}_{x}\mathfrak{D}_{x}-\frac{1}{2}(\mathfrak{ED})+\mathfrak{H}_{x}\mathfrak{B}_{x}-\frac{1}{2}(\mathfrak{HB})\\ \overline{Y_{y}^{\prime}} & =\overline{Y_{x}^{\prime}}=\frac{1}{2}\left\{ \mathfrak{E}_{x}\mathfrak{D}_{y}+\mathfrak{E}_{y}\mathfrak{D}_{x}+\mathfrak{H}_{x}\mathfrak{B}_{y}+\mathfrak{H}_{y}\mathfrak{B}_{x}\right\} \\ & \dots\\ \overline{X_{x}^{\prime\prime\prime}} & =\overline{X_{x}^{\prime\prime}}=\overline{X_{x}^{\prime}}\\ \overline{X_{y}^{\prime\prime\prime}} & =\overline{Y_{x}^{\prime\prime}}=\mathfrak{E}_{y}\mathfrak{D}_{x}+\mathfrak{H}_{y}\mathfrak{B}_{x}\\ \overline{Y_{x}^{\prime\prime\prime}} & =\overline{X_{y}^{\prime\prime}}=\mathfrak{E}_{x}\mathfrak{D}_{y}+\mathfrak{H}_{x}\mathfrak{B}_{y}\\ & \dots \end{align} $$

leading to the expressions for moving bodies:

$$\begin{matrix}\left\{ \begin{align}X_{x}^{\prime} & =\overline{X_{x}^{\prime}}+\mathfrak{v}_{x}\mathfrak{W}_{x},\ Y_{x}^{\prime}=\overline{Y_{x}^{\prime}}+\frac{1}{2}\left\{ \mathfrak{v}_{x}\mathfrak{W}_{y}+\mathfrak{v}_{y}\mathfrak{W}_{x}\right\} \\ X_{y}^{\prime} & =Y_{x}^{\prime}\dots,\\ \mathfrak{S}' & =\mathfrak{g}'=\frac{1}{2}\left\{ \mathfrak{U}+\mathfrak{S}-\mathfrak{W}-\mathfrak{v}(\mathfrak{vW})\right\} \\ E' & =E-(\mathfrak{vW}) \end{align} \right.\\ \left\{ \begin{align}X_{x}^{\prime\prime} & =\overline{X_{x}^{\prime\prime}}+\mathfrak{v}_{x}\mathfrak{W}_{x},\ Y_{x}^{\prime\prime}=\overline{Y_{x}^{\prime\prime}}+\mathfrak{v}_{y}\mathfrak{W}_{x},\\ X_{y}^{\prime\prime} & =X_{y}^{\prime\prime}+\mathfrak{v}_{y}\mathfrak{W}_{y}\dots\\ \mathfrak{S} & =\mathfrak{S}-\mathfrak{v}(\mathfrak{vW}),\ \mathfrak{g}=\mathfrak{U}-\mathfrak{W},\\ E'' & =E-(\mathfrak{vW}) \end{align} \right.\\ \left\{ \begin{align}X_{x}^{\prime\prime\prime} & =\overline{X_{x}^{\prime\prime}}+\mathfrak{v}_{x}\mathfrak{W}_{x},\ Y_{x}^{\prime\prime}=\overline{Y_{x}^{\prime\prime}}+\mathfrak{v}_{y}\mathfrak{W}_{x},\\ X_{y}^{\prime\prime\prime} & =Y_{x}^{\prime\prime}\dots\\ \mathfrak{S} & =\mathfrak{g},\ \mathfrak{g}=\mathfrak{S},\\ E' & =E \end{align} \right.\\ \left[\begin{matrix}\mathfrak{W}'=\varkappa^{2}\left\{ [\mathfrak{E'H'}]+\varkappa^{2}\mathfrak{v}(\mathfrak{v}[\mathfrak{E'H'}])\right\} \\ \mathfrak{W}=(\epsilon\mu-1)\mathfrak{W}',\ \mathfrak{U}=[\mathfrak{DB}],\\ \mathfrak{E}'=\mathfrak{E}+[\mathfrak{vB}],\ \mathfrak{H}'=\mathfrak{H}-[\mathfrak{vD}]\\ \varkappa=\frac{1}{\sqrt{1-\mathfrak{v}^{2}}} \end{matrix}\right] \end{matrix}$$