w:Wilhelm Killing discussed Newtonian mechanics in non-Euclidean spaces by expressing position, velocity, acceleration and force in terms of four components. He expressed the four components of force R=(P,X,Y,Z) , its norm, its inner product with coordinates (p,x,y,z) , and the equations of motion as the product of mass with acceleration as follows:[ M 1]
R
⇒
P
,
X
,
Y
,
Z
R
2
=
P
2
k
2
+
X
2
+
Y
2
+
Z
2
p
P
+
x
X
+
y
Y
+
z
Z
=
0
m
k
2
d
2
p
d
t
2
=
P
,
m
d
2
x
d
t
2
=
X
,
m
d
2
y
d
t
2
=
Y
,
m
d
2
z
d
t
2
=
Z
{\displaystyle {\begin{matrix}R\Rightarrow P,X,Y,Z\\\hline R^{2}={\frac {P^{2}}{k^{2}}}+X^{2}+Y^{2}+Z^{2}\\pP+xX+yY+zZ=0\\\hline mk^{2}{\frac {d^{2}p}{dt^{2}}}=P,\ m{\frac {d^{2}x}{dt^{2}}}=X,\ m{\frac {d^{2}y}{dt^{2}}}=Y,\ m{\frac {d^{2}z}{dt^{2}}}=Z\end{matrix}}}
If the Gaussian curvature
1
/
k
2
{\displaystyle 1/k^{2}}
(with k as radius of curvature) is negative the force becomes related to the hyperboloid model of hyperbolic space, which at first sight becomes similar to the relativistic four-force in Minkowski space by setting
k
2
=
−
c
2
{\displaystyle k^{2}=-c^{2}}
with c as speed of light. However, Killing obtained his results by differentiation with respect to Newtonian time t , not relativistic proper time, so his expressions aren't relativistic four-vectors in the first place, in particular they don't involve a limiting speed.
w:Henri Poincaré (July 1905, published January 1906) argued, that the Lorentz transformation not only applies to electrodynamics, but to all other phenomena as well including mechanics. For instance, he explicitly defined gravitation as non-electromagnetic in origin and applied the following expression of four-force to his Lorentz invariant model of gravitation:[ R 1]
(
k
0
X
1
,
k
0
Y
1
,
k
0
Z
1
,
k
0
T
1
)
{\displaystyle \left(k_{0}X_{1},\quad k_{0}Y_{1},\quad k_{0}Z_{1},\quad k_{0}T_{1}\right)}
with
T
1
=
Σ
X
1
ξ
{\displaystyle T_{1}=\Sigma X_{1}\xi }
and
k
0
=
1
1
−
ϵ
2
{\displaystyle k_{0}={\tfrac {1}{\sqrt {1-\epsilon ^{2}}}}}
equivalent to (b) because
ϵ
=
v
c
,
(
X
1
,
Y
1
,
Z
1
)
=
f
,
Σ
X
1
ξ
=
f
⋅
v
{\displaystyle \epsilon ={\frac {v}{c}},\ \left(X_{1},\ Y_{1},\ Z_{1}\right)=\mathbf {f} ,\ \Sigma X_{1}\xi =\mathbf {f} \cdot \mathbf {v} }
.
In an appendix to his lecture from December 1907 (published 1908), w:Hermann Minkowski extended the postulate of relativity to mechanics in general, defining four-force density
K
+
(
K
w
¯
)
w
{\displaystyle K+(K{\overline {w}})w}
with w as four-velocity and S as stress energy tensor[ R 2]
K
+
(
K
w
¯
)
w
K
=
lor
S
⇒
K
h
=
∂
S
1
h
∂
x
1
+
∂
S
2
h
∂
x
2
+
∂
S
3
h
∂
x
3
+
∂
S
4
h
∂
x
4
[
lor
=
|
∂
∂
x
1
,
∂
∂
x
2
,
∂
∂
x
3
,
∂
x
4
|
]
{\displaystyle {\begin{matrix}K+(K{\overline {w}})w\\\hline K={\text{lor }}S\Rightarrow K_{h}={\frac {\partial S_{1h}}{\partial x_{1}}}+{\frac {\partial S_{2h}}{\partial x_{2}}}+{\frac {\partial S_{3h}}{\partial x_{3}}}+{\frac {\partial S_{4h}}{\partial x_{4}}}\\\left[{\text{lor }}=\left|{\frac {\partial }{\partial x_{1}}},\ {\frac {\partial }{\partial x_{2}}},\ {\frac {\partial }{\partial x_{3}}},\ {\frac {\partial }{x_{4}}}\right|\right]\end{matrix}}}
equivalent to (d1). He went on to show that these relations can be also used to define the equations of motion of mechanics in terms of constant rest mass density
ν
{\displaystyle \nu }
:[ R 3]
∂
ν
w
h
w
1
∂
x
1
+
∂
ν
w
h
w
2
∂
x
2
+
∂
ν
w
h
w
3
∂
x
3
+
∂
ν
w
h
w
4
∂
x
4
ν
d
w
h
d
τ
⇒
ν
d
d
τ
(
d
x
h
d
τ
)
⇒
K
+
(
K
w
¯
)
w
=
(
X
,
Y
,
Z
,
i
T
)
ν
d
d
τ
d
x
d
τ
=
X
,
ν
d
d
τ
d
y
d
τ
=
Y
,
ν
d
d
τ
d
z
d
τ
=
Z
,
ν
d
d
τ
d
t
d
τ
=
T
{\displaystyle {\begin{matrix}{\frac {\partial \nu \,w_{h}w_{1}}{\partial x_{1}}}+{\frac {\partial \nu \,w_{h}w_{2}}{\partial x_{2}}}+{\frac {\partial \nu \,w_{h}w_{3}}{\partial x_{3}}}+{\frac {\partial \nu \,w_{h}w_{4}}{\partial x_{4}}}\\\nu {\frac {dw_{h}}{d\tau }}\Rightarrow \nu {\frac {d}{d\tau }}\left({\frac {dx_{h}}{d\tau }}\right)\Rightarrow K+(K\,{\overline {w}})w=(X,Y,Z,iT)\\\hline \nu {\frac {d}{d\tau }}{\frac {dx}{d\tau }}=X,\ \nu {\frac {d}{d\tau }}{\frac {dy}{d\tau }}=Y,\ \nu {\frac {d}{d\tau }}{\frac {dz}{d\tau }}=Z,\ \nu {\frac {d}{d\tau }}{\frac {dt}{d\tau }}=T\end{matrix}}}
equivalent to (b1, c1, d1) as well as (f) since the first line includes the dust stress-energy tensor. Eventually he defined a “moving force” as the product of constant rest mass and four-acceleration[ R 4]
m
d
d
τ
d
x
d
τ
=
R
x
,
m
d
d
τ
d
y
d
τ
=
R
y
,
m
d
d
τ
d
z
d
τ
=
R
z
,
m
d
d
τ
d
t
d
τ
=
R
t
{\displaystyle m{\frac {d}{d\tau }}{\frac {dx}{d\tau }}=R_{x},\quad m{\frac {d}{d\tau }}{\frac {dy}{d\tau }}=R_{y},\quad m{\frac {d}{d\tau }}{\frac {dz}{d\tau }}=R_{z},\quad m{\frac {d}{d\tau }}{\frac {dt}{d\tau }}=R_{t}}
equivalent to (c). Minkowski's assumption of constant rest mass was later challenged by Abraham (see next section).
In 1909, w:Max Abraham pointed out that the relativity principle requires that the mechanical forces must transform like the electromagnetic ones, so there must be a four-dimensional tensor for mechanics (i.e. mechanical stress energy tensor) in analogy to the electromagnetic one, and that the relation
c
2
g
=
S
{\displaystyle c^{2}{\mathfrak {g}}={\mathfrak {S}}}
can alternatively be interpreted as relation between mechanical momentum and energy density:[ R 5]
K
x
=
∂
X
x
∂
x
+
∂
X
y
∂
y
+
∂
X
z
∂
z
+
∂
X
u
∂
u
K
y
=
∂
Y
x
∂
x
+
∂
Y
y
∂
y
+
∂
Y
z
∂
z
+
∂
Y
u
∂
u
K
z
=
∂
Z
x
∂
x
+
∂
Z
y
∂
y
+
∂
Z
z
∂
z
+
∂
Z
u
∂
u
K
u
=
∂
U
x
∂
x
+
∂
U
y
∂
y
+
∂
U
z
∂
z
+
∂
U
u
∂
u
[
u
=
i
c
t
]
{\displaystyle {\begin{matrix}{\begin{aligned}{\mathfrak {K}}_{x}&={\frac {\partial X_{x}}{\partial x}}+{\frac {\partial X_{y}}{\partial y}}+{\frac {\partial X_{z}}{\partial z}}+{\frac {\partial X_{u}}{\partial u}}\\{\mathfrak {K}}_{y}&={\frac {\partial Y_{x}}{\partial x}}+{\frac {\partial Y_{y}}{\partial y}}+{\frac {\partial Y_{z}}{\partial z}}+{\frac {\partial Y_{u}}{\partial u}}\\{\mathfrak {K}}_{z}&={\frac {\partial Z_{x}}{\partial x}}+{\frac {\partial Z_{y}}{\partial y}}+{\frac {\partial Z_{z}}{\partial z}}+{\frac {\partial Z_{u}}{\partial u}}\\{\mathfrak {K}}_{u}&={\frac {\partial U_{x}}{\partial x}}+{\frac {\partial U_{y}}{\partial y}}+{\frac {\partial U_{z}}{\partial z}}+{\frac {\partial U_{u}}{\partial u}}\end{aligned}}\\\left[u=ict\right]\end{matrix}}}
equivalent to (a) when interpreted as mechanical force. While Minkowski (1907) assumed constant rest mass density, Abraham (1909) held that mass-energy equivalence, according to which mass depends on its energy content, would suggest a variable rest mass
m
0
{\displaystyle m_{0}}
.
In 1912, Abraham introduced the expression “world tensor of motion”
T
∗
{\displaystyle T^{\ast }}
(equivalent to the dust tensor) while formulating his first theory of gravitation. It has ten components representing kinetic stresses, energy flux
S
∗
{\displaystyle {\mathfrak {S}}^{\ast }}
and momentum
g
∗
{\displaystyle {\mathfrak {g}}^{\ast }}
of matter in terms of rest mass density
ν
{\displaystyle \nu }
, which he combined with the world tensor
T
{\displaystyle T}
(representing the electromagnetic-, gravitational-, and stress field) in order to formulate the momentum and energy conservation theorems:[ R 6]
∂
(
c
g
x
+
c
g
x
∗
)
c
∂
t
=
∂
(
X
x
+
X
x
∗
)
∂
x
+
∂
(
X
y
+
X
y
∗
)
∂
y
+
∂
(
X
z
+
X
z
∗
)
∂
z
∂
(
c
g
y
+
c
g
y
∗
)
c
∂
t
=
∂
(
Y
x
+
Y
x
∗
)
∂
x
+
∂
(
Y
y
+
Y
y
∗
)
∂
y
+
∂
(
Y
z
+
Y
z
∗
)
∂
z
∂
(
c
g
z
+
c
g
z
∗
)
c
∂
t
=
∂
(
Z
x
+
Z
x
∗
)
∂
x
+
∂
(
Z
y
+
Z
y
∗
)
∂
y
+
∂
(
Z
z
+
Z
z
∗
)
∂
z
∂
(
ε
+
ε
∗
)
∂
t
=
−
c
d
i
v
{
S
+
S
∗
c
}
{\displaystyle {\begin{aligned}{\frac {\partial \left(c{\mathfrak {g}}_{x}+c{\mathfrak {g}}_{x}^{\ast }\right)}{c\partial t}}&={\frac {\partial \left(X_{x}+X_{x}^{\ast }\right)}{\partial x}}+{\frac {\partial \left(X_{y}+X_{y}^{\ast }\right)}{\partial y}}+{\frac {\partial \left(X_{z}+X_{z}^{\ast }\right)}{\partial z}}\\{\frac {\partial \left(c{\mathfrak {g}}_{y}+c{\mathfrak {g}}_{y}^{\ast }\right)}{c\partial t}}&={\frac {\partial \left(Y_{x}+Y_{x}^{\ast }\right)}{\partial x}}+{\frac {\partial \left(Y_{y}+Y_{y}^{\ast }\right)}{\partial y}}+{\frac {\partial \left(Y_{z}+Y_{z}^{\ast }\right)}{\partial z}}\\{\frac {\partial \left(c{\mathfrak {g}}_{z}+c{\mathfrak {g}}_{z}^{\ast }\right)}{c\partial t}}&={\frac {\partial \left(Z_{x}+Z_{x}^{\ast }\right)}{\partial x}}+{\frac {\partial \left(Z_{y}+Z_{y}^{\ast }\right)}{\partial y}}+{\frac {\partial \left(Z_{z}+Z_{z}^{\ast }\right)}{\partial z}}\\{\frac {\partial \left(\varepsilon +\varepsilon ^{\ast }\right)}{\partial t}}&=-c{\rm {div}}\left\{{\frac {{\mathfrak {S}}+{\mathfrak {S}}^{\ast }}{c}}\right\}\end{aligned}}}
equivalent to (f) when only the dust tensor
T
∗
{\displaystyle T^{\ast }}
is considered.
The force definitions of both Abraham (variable rest mass) and Minkowski (constant rest mass) were elaborated by w:Gunnar Nordström (1910), who defined two variants of four-force density using a “four-dimensional tensor” (i.e. dust solution)
γ
a
m
a
n
{\displaystyle \gamma {\mathfrak {a}}_{m}{\mathfrak {a}}_{n}}
consisting of rest mass density
γ
{\displaystyle \gamma }
and four-velocity
a
{\displaystyle {\mathfrak {a}}}
. The first formulation was based on Abraham's assumption of variable rest mass density:[ R 7]
K
x
′
=
∂
∂
x
γ
a
x
2
+
∂
∂
y
γ
a
x
a
y
+
∂
∂
z
γ
a
x
a
z
+
∂
∂
u
γ
a
x
a
u
K
y
′
=
∂
∂
x
γ
a
y
a
x
+
∂
∂
y
γ
a
y
2
+
∂
∂
z
γ
a
y
a
z
+
∂
∂
u
γ
a
y
a
u
K
z
′
=
∂
∂
x
γ
a
z
a
x
+
∂
∂
y
γ
a
z
a
y
+
∂
∂
z
γ
a
z
2
+
∂
∂
u
γ
a
z
a
u
K
u
′
=
∂
∂
x
γ
a
u
a
x
+
∂
∂
y
γ
a
u
a
y
+
∂
∂
z
γ
a
u
a
z
+
∂
∂
u
γ
a
u
2
(
u
=
i
c
t
)
{\displaystyle {\begin{aligned}{\mathfrak {K}}_{x}^{\prime }&={\frac {\partial }{\partial x}}\gamma {\mathfrak {a}}_{x}^{2}+{\frac {\partial }{\partial y}}\gamma {\mathfrak {a}}_{x}{\mathfrak {a}}_{y}+{\frac {\partial }{\partial z}}\gamma {\mathfrak {a}}_{x}{\mathfrak {a}}_{z}+{\frac {\partial }{\partial u}}\gamma {\mathfrak {a}}_{x}{\mathfrak {a}}_{u}\\{\mathfrak {K}}_{y}^{\prime }&={\frac {\partial }{\partial x}}\gamma {\mathfrak {a}}_{y}{\mathfrak {a}}_{x}+{\frac {\partial }{\partial y}}\gamma {\mathfrak {a}}_{y}^{2}+{\frac {\partial }{\partial z}}\gamma {\mathfrak {a}}_{y}{\mathfrak {a}}_{z}+{\frac {\partial }{\partial u}}\gamma {\mathfrak {a}}_{y}{\mathfrak {a}}_{u}\\{\mathfrak {K}}_{z}^{\prime }&={\frac {\partial }{\partial x}}\gamma {\mathfrak {a}}_{z}{\mathfrak {a}}_{x}+{\frac {\partial }{\partial y}}\gamma {\mathfrak {a}}_{z}{\mathfrak {a}}_{y}+{\frac {\partial }{\partial z}}\gamma {\mathfrak {a}}_{z}^{2}+{\frac {\partial }{\partial u}}\gamma {\mathfrak {a}}_{z}{\mathfrak {a}}_{u}\\{\mathfrak {K}}_{u}^{\prime }&={\frac {\partial }{\partial x}}\gamma {\mathfrak {a}}_{u}{\mathfrak {a}}_{x}+{\frac {\partial }{\partial y}}\gamma {\mathfrak {a}}_{u}{\mathfrak {a}}_{y}+{\frac {\partial }{\partial z}}\gamma {\mathfrak {a}}_{u}{\mathfrak {a}}_{z}+{\frac {\partial }{\partial u}}\gamma {\mathfrak {a}}_{u}^{2}\\&(u=ict)\end{aligned}}}
equivalent to (f), and the second one on Minkowski's assumption of constant rest mass density:[ R 8]
K
x
′
′
=
γ
a
x
∂
a
x
∂
x
+
γ
a
y
∂
a
x
∂
y
+
γ
a
z
∂
a
x
∂
z
+
γ
a
u
∂
a
x
∂
u
K
y
′
′
=
γ
a
x
∂
a
y
∂
x
+
γ
a
y
∂
a
y
∂
y
+
γ
a
z
∂
a
y
∂
z
+
γ
a
u
∂
a
y
∂
u
K
z
′
′
=
γ
a
x
∂
a
z
∂
x
+
γ
a
y
∂
a
z
∂
y
+
γ
a
z
∂
a
z
∂
z
+
γ
a
u
∂
a
z
∂
u
K
u
′
′
=
γ
a
x
∂
a
u
∂
x
+
γ
a
y
∂
a
u
∂
y
+
γ
a
z
∂
a
u
∂
z
+
γ
a
u
∂
a
u
∂
u
{\displaystyle {\begin{aligned}{\mathfrak {K}}_{x}^{\prime \prime }&=\gamma {\mathfrak {a}}_{x}{\frac {\partial {\mathfrak {a}}_{x}}{\partial x}}+\gamma {\mathfrak {a}}_{y}{\frac {\partial {\mathfrak {a}}_{x}}{\partial y}}+\gamma {\mathfrak {a}}_{z}{\frac {\partial {\mathfrak {a}}_{x}}{\partial z}}+\gamma {\mathfrak {a}}_{u}{\frac {\partial {\mathfrak {a}}_{x}}{\partial u}}\\{\mathfrak {K}}_{y}^{\prime \prime }&=\gamma {\mathfrak {a}}_{x}{\frac {\partial {\mathfrak {a}}_{y}}{\partial x}}+\gamma {\mathfrak {a}}_{y}{\frac {\partial {\mathfrak {a}}_{y}}{\partial y}}+\gamma {\mathfrak {a}}_{z}{\frac {\partial {\mathfrak {a}}_{y}}{\partial z}}+\gamma {\mathfrak {a}}_{u}{\frac {\partial {\mathfrak {a}}_{y}}{\partial u}}\\{\mathfrak {K}}_{z}^{\prime \prime }&=\gamma {\mathfrak {a}}_{x}{\frac {\partial {\mathfrak {a}}_{z}}{\partial x}}+\gamma {\mathfrak {a}}_{y}{\frac {\partial {\mathfrak {a}}_{z}}{\partial y}}+\gamma {\mathfrak {a}}_{z}{\frac {\partial {\mathfrak {a}}_{z}}{\partial z}}+\gamma {\mathfrak {a}}_{u}{\frac {\partial {\mathfrak {a}}_{z}}{\partial u}}\\{\mathfrak {K}}_{u}^{\prime \prime }&=\gamma {\mathfrak {a}}_{x}{\frac {\partial {\mathfrak {a}}_{u}}{\partial x}}+\gamma {\mathfrak {a}}_{y}{\frac {\partial {\mathfrak {a}}_{u}}{\partial y}}+\gamma {\mathfrak {a}}_{z}{\frac {\partial {\mathfrak {a}}_{u}}{\partial z}}+\gamma {\mathfrak {a}}_{u}{\frac {\partial {\mathfrak {a}}_{u}}{\partial u}}\end{aligned}}}
equivalent to (f). In 1911, Nordström only used the first variant with variable rest mass density and considered pressure in a “material fluid” as well.[ R 9]
In 1913, he added an “elastic stress tensor” p in order to reformulated Laue's symmetrical four dimensional tensor T representing spatial stresses and mechanical momentum and energy density, which he used to add an elastic component
K
e
{\displaystyle {\mathfrak {K}}^{e}}
to the four-force-density
K
{\displaystyle {\mathfrak {K}}}
to give the equation of motion:[ R 10]
K
x
+
K
x
e
=
∂
∂
x
ν
B
x
2
+
∂
∂
y
ν
B
x
B
y
+
∂
∂
z
ν
B
x
B
z
+
∂
∂
u
ν
B
x
B
u
K
y
+
K
y
e
=
∂
∂
x
ν
B
y
B
x
+
∂
∂
y
ν
B
y
2
+
∂
∂
z
ν
B
y
B
z
+
∂
∂
u
ν
B
y
B
u
K
z
+
K
z
e
=
∂
∂
x
ν
B
z
B
x
+
∂
∂
y
ν
B
z
B
y
+
∂
∂
z
ν
B
z
2
+
∂
∂
u
ν
B
z
B
u
K
u
+
K
u
e
=
∂
∂
x
ν
B
u
B
x
+
∂
∂
y
ν
B
u
B
y
+
∂
∂
z
ν
B
u
B
z
+
∂
∂
u
ν
B
u
2
{\displaystyle {\begin{aligned}{\mathfrak {K}}_{x}+{\mathfrak {K}}_{x}^{e}&={\frac {\partial }{\partial x}}\nu {\mathfrak {B}}_{x}^{2}+{\frac {\partial }{\partial y}}\nu {\mathfrak {B}}_{x}{\mathfrak {B}}_{y}+{\frac {\partial }{\partial z}}\nu {\mathfrak {B}}_{x}{\mathfrak {B}}_{z}+{\frac {\partial }{\partial u}}\nu {\mathfrak {B}}_{x}{\mathfrak {B}}_{u}\\{\mathfrak {K}}_{y}+{\mathfrak {K}}_{y}^{e}&={\frac {\partial }{\partial x}}\nu {\mathfrak {B}}_{y}{\mathfrak {B}}_{x}+{\frac {\partial }{\partial y}}\nu {\mathfrak {B}}_{y}^{2}+{\frac {\partial }{\partial z}}\nu {\mathfrak {B}}_{y}{\mathfrak {B}}_{z}+{\frac {\partial }{\partial u}}\nu {\mathfrak {B}}_{y}{\mathfrak {B}}_{u}\\{\mathfrak {K}}_{z}+{\mathfrak {K}}_{z}^{e}&={\frac {\partial }{\partial x}}\nu {\mathfrak {B}}_{z}{\mathfrak {B}}_{x}+{\frac {\partial }{\partial y}}\nu {\mathfrak {B}}_{z}{\mathfrak {B}}_{y}+{\frac {\partial }{\partial z}}\nu {\mathfrak {B}}_{z}^{2}+{\frac {\partial }{\partial u}}\nu {\mathfrak {B}}_{z}{\mathfrak {B}}_{u}\\{\mathfrak {K}}_{u}+{\mathfrak {K}}_{u}^{e}&={\frac {\partial }{\partial x}}\nu {\mathfrak {B}}_{u}{\mathfrak {B}}_{x}+{\frac {\partial }{\partial y}}\nu {\mathfrak {B}}_{u}{\mathfrak {B}}_{y}+{\frac {\partial }{\partial z}}\nu {\mathfrak {B}}_{u}{\mathfrak {B}}_{z}+{\frac {\partial }{\partial u}}\nu {\mathfrak {B}}_{u}^{2}\end{aligned}}}
equivalent to (f). He went on to employ this notion in his theory of gravitation.
w:Wladimir Ignatowski derived the hydrodynamic four-force density in terms of mass density
ϱ
{\displaystyle \varrho }
and pressure
π
{\displaystyle \pi }
in case of perfect fluids:[ R 11]
(
K
1
,
n
v
K
1
+
1
−
n
v
2
D
)
=
(
−
Δ
π
,
n
∂
π
∂
t
)
=
(
−
Δ
π
,
n
v
Δ
π
+
n
d
π
d
t
)
[
K
1
=
ϱ
−
d
n
v
1
−
n
v
2
d
t
+
v
1
−
n
v
2
D
=
−
Δ
π
,
D
=
n
1
−
n
v
2
d
π
d
t
=
d
ϱ
d
t
+
ϱ
d
i
v
v
,
n
=
1
c
2
,
Δ
=
g
r
a
d
]
{\displaystyle {\begin{matrix}\left({\mathfrak {K}}_{1},\ n{\mathfrak {v}}{\mathfrak {K}}_{1}+{\sqrt {1-n{\mathfrak {v}}^{2}}}D\right)=\left(-\Delta \pi ,\ n{\frac {\partial \pi }{\partial t}}\right)=\left(-\Delta \pi ,\ n{\mathfrak {v}}\Delta \pi +n{\frac {d\pi }{dt}}\right)\\\hline \left[{\mathfrak {K}}_{1}=\varrho -{\frac {d{\frac {n{\mathfrak {v}}}{\sqrt {1-n{\mathfrak {v}}^{2}}}}}{dt}}+{\frac {\mathfrak {v}}{\sqrt {1-n{\mathfrak {v}}^{2}}}}D=-\Delta \pi ,\ D={\frac {n}{\sqrt {1-n{\mathfrak {v}}^{2}}}}{\frac {d\pi }{dt}}={\frac {d\varrho }{dt}}+\varrho \mathrm {div} {\mathfrak {v}},\ n={\frac {1}{c^{2}}},\ \Delta =\mathrm {grad} \right]\end{matrix}}}
which corresponds to the four-force of perfect fluid (e).
In the first textbook on relativity (1911), w:Max von Laue defined the mechanical ponderomotive force F based on world tensor T (i.e. mechanical stress-energy tensor), implying the complete reduction of mechanical inertia to energy and stresses:[ R 12]
F
=
−
Δ
i
v
T
{\displaystyle F=-\varDelta iv\,T}
equivalent to (c).
In the second edition (1912, published 1913), he followed #Herglotz (1911) and #Lamla (1911/12) in defining the four-force K in order to produce the Lagrangian (1) and Eulerian (2) fundamental equations of hydrodynamics, as well as for the case of least compressibility (3), where
δ
m
n
{\displaystyle \delta m_{n}}
is the normal rest mass, p the pressure, Y the four-velocity:[ R 13]
(
1
)
d
d
τ
{
δ
m
n
(
1
+
P
c
2
)
Y
}
+
δ
V
0
c
Γ
ϱ
α
δ
p
=
K
(
2
)
Y
x
∂
∂
x
(
ϰ
0
Y
m
)
+
Y
y
∂
∂
y
(
ϰ
0
Y
m
)
+
Y
z
∂
∂
z
(
ϰ
0
Y
m
)
+
Y
l
∂
∂
l
(
ϰ
0
Y
m
)
+
1
c
2
Γ
ϱ
α
δ
m
P
=
K
m
c
δ
m
n
(
3
)
k
0
(
1
c
⋅
d
Y
d
τ
−
Y
D
i
v
Y
)
+
Γ
ϱ
α
δ
k
0
=
k
n
0
K
c
δ
m
n
(
l
=
i
c
t
,
ϰ
0
=
1
+
P
c
2
,
P
=
∫
p
n
p
d
p
k
0
,
k
0
=
δ
m
n
δ
V
0
,
Γ
ϱ
α
δ
=
four-gradient
)
{\displaystyle {\begin{matrix}(1)&{\frac {d}{d\tau }}\left\{\delta m_{n}\left(1+{\frac {P}{c^{2}}}\right)Y\right\}+{\frac {\delta V^{0}}{c}}\Gamma \varrho \alpha \delta \ p=K\\(2)&Y_{x}{\frac {\partial }{\partial x}}\left(\varkappa ^{0}Y_{m}\right)+Y_{y}{\frac {\partial }{\partial y}}\left(\varkappa ^{0}Y_{m}\right)+Y_{z}{\frac {\partial }{\partial z}}\left(\varkappa ^{0}Y_{m}\right)+Y_{l}{\frac {\partial }{\partial l}}\left(\varkappa ^{0}Y_{m}\right)+{\frac {1}{c^{2}}}\Gamma \varrho \alpha \delta _{m}P={\frac {K_{m}}{c\delta m_{n}}}\\(3)&k^{0}\left({\frac {1}{c}}\cdot {\frac {dY}{d\tau }}-Y\ Div\ Y\right)+\Gamma \varrho \alpha \delta \ k^{0}={\frac {k_{n}^{0}K}{c\delta m_{n}}}\\&\left(l=ict,\ \varkappa ^{0}=1+{\frac {P}{c^{2}}},\ P=\int _{p_{n}}^{p}{\frac {dp}{k^{0}}},\ k^{0}={\frac {\delta m_{n}}{\delta V^{0}}},\ \Gamma \varrho \alpha \delta ={\text{four-gradient}}\right)\end{matrix}}}
equivalent to (c).
In the fourth edition (1921), he defined the four-force density using the kinetic stress energy tensor, with rest energy density
W
0
{\displaystyle W^{0}}
, rest energy
E
0
{\displaystyle E^{0}}
, rest volume
V
0
{\displaystyle V^{0}}
:[ R 14]
T
j
k
=
W
0
Y
j
Y
k
Δ
i
v
x
T
=
Y
x
D
i
v
(
W
0
Y
)
+
W
0
(
Y
x
∂
Y
x
∂
x
+
Y
y
∂
Y
x
∂
x
+
Y
z
∂
Y
x
∂
x
+
Y
l
∂
Y
x
∂
x
)
Δ
i
v
x
T
=
W
0
c
d
Y
d
τ
=
−
F
∫
T
i
k
δ
V
0
=
E
0
Y
i
Y
k
(
l
=
i
c
t
,
d
i
v
=
four-divergence of four-vector
Δ
i
v
=
four-divergence of six-vector
)
{\displaystyle {\begin{matrix}T_{jk}=W^{0}Y_{j}Y_{k}\\\varDelta iv_{x}T=Y_{x}Div\left(W^{0}Y\right)+W^{0}\left(Y_{x}{\frac {\partial Y_{x}}{\partial x}}+Y_{y}{\frac {\partial Y_{x}}{\partial x}}+Y_{z}{\frac {\partial Y_{x}}{\partial x}}+Y_{l}{\frac {\partial Y_{x}}{\partial x}}\right)\\\varDelta iv_{x}T={\frac {W^{0}}{c}}{\frac {dY}{d\tau }}=-F\\\int T_{ik}\delta V^{0}=E^{0}Y_{i}Y_{k}\\\left(l=ict,\ {\begin{matrix}div={\text{four-divergence of four-vector}}\\\varDelta iv={\text{four-divergence of six-vector}}\end{matrix}}\right)\end{matrix}}}
equivalent to (b,c) and the dust solution (f).
w:Gustav Herglotz gave a complete theory of elasticity in special relativity, including equations of motion in different forms, which he defined using coordinates
x
…
{\displaystyle x\dots }
after deformation,
ξ
…
{\displaystyle \xi \dots }
and
ξ
0
…
{\displaystyle \xi ^{0}\dots }
before deformation, from which he derived the deformation quantities
a
i
j
{\displaystyle a_{ij}}
and
A
i
j
{\displaystyle A_{ij}}
, together with the kinetic potential
Φ
{\displaystyle \varPhi }
. He gave the Lagrangian equation of motion in terms of four-force density:[ R 15]
1
a
44
Ξ
,
1
a
44
H
,
1
a
44
Z
,
1
a
44
T
Ξ
=
∂
Φ
11
∂
ξ
+
∂
Φ
12
∂
η
+
∂
Φ
13
∂
ζ
+
∂
Φ
14
∂
τ
H
=
∂
Φ
21
∂
ξ
+
∂
Φ
22
∂
η
+
∂
Φ
23
∂
ζ
+
∂
Φ
24
∂
τ
Z
=
∂
Φ
31
∂
ξ
+
∂
Φ
32
∂
η
+
∂
Φ
33
∂
ζ
+
∂
Φ
34
∂
τ
T
=
∂
Φ
41
∂
ξ
+
∂
Φ
42
∂
η
+
∂
Φ
43
∂
ζ
+
∂
Φ
44
∂
τ
=
−
(
u
Ξ
+
v
H
+
w
Z
+
ϑ
∂
ε
∂
τ
)
E
=
−
∂
Φ
∂
ε
=
a
44
ϑ
{\displaystyle {\begin{matrix}{\frac {1}{a_{44}}}\Xi ,\ {\frac {1}{a_{44}}}\mathrm {H} ,\ {\frac {1}{a_{44}}}\mathrm {Z} ,\ {\frac {1}{a_{44}}}\mathrm {T} \\\hline {\begin{aligned}{\begin{aligned}\Xi &={\frac {\partial \varPhi _{11}}{\partial \xi }}+{\frac {\partial \varPhi _{12}}{\partial \eta }}+{\frac {\partial \varPhi _{13}}{\partial \zeta }}+{\frac {\partial \varPhi _{14}}{\partial \tau }}\\\mathrm {H} &={\frac {\partial \varPhi _{21}}{\partial \xi }}+{\frac {\partial \varPhi _{22}}{\partial \eta }}+{\frac {\partial \varPhi _{23}}{\partial \zeta }}+{\frac {\partial \varPhi _{24}}{\partial \tau }}\\\mathrm {Z} &={\frac {\partial \varPhi _{31}}{\partial \xi }}+{\frac {\partial \varPhi _{32}}{\partial \eta }}+{\frac {\partial \varPhi _{33}}{\partial \zeta }}+{\frac {\partial \varPhi _{34}}{\partial \tau }}\\\mathrm {T} &={\frac {\partial \varPhi _{41}}{\partial \xi }}+{\frac {\partial \varPhi _{42}}{\partial \eta }}+{\frac {\partial \varPhi _{43}}{\partial \zeta }}+{\frac {\partial \varPhi _{44}}{\partial \tau }}\\&=-\left(u\Xi +v\mathrm {H} +w\mathrm {Z} +\vartheta {\frac {\partial \varepsilon }{\partial \tau }}\right)\\\mathrm {E} &=-{\frac {\partial \varPhi }{\partial \varepsilon }}=a_{44}\vartheta \end{aligned}}\end{aligned}}\end{matrix}}}
and the Euler equations of motion by defining stress-energy tensor
F
i
j
{\displaystyle F_{ij}}
, whose components can be related to momentum density
X
,
Y
,
Z
{\displaystyle {\mathfrak {X}},{\mathfrak {Y}},{\mathfrak {Z}}}
, energy density
E
{\displaystyle {\mathfrak {E}}}
, velocity u,v,w :[ R 16]
X
=
∂
F
11
∂
x
+
∂
F
12
∂
y
+
∂
F
13
∂
z
+
∂
F
14
∂
t
Y
=
∂
F
21
∂
x
+
∂
F
22
∂
y
+
∂
F
23
∂
z
+
∂
F
24
∂
t
Z
=
∂
F
31
∂
x
+
∂
F
32
∂
y
+
∂
F
33
∂
z
+
∂
F
34
∂
t
T
=
∂
F
41
∂
x
+
∂
F
42
∂
y
+
∂
F
43
∂
z
+
∂
F
44
∂
t
=
−
(
u
X
+
v
Y
+
w
Z
+
Q
)
ϑ
=
Δ
1
−
s
2
∂
F
∂
ε
=
a
44
ϑ
{\displaystyle {\begin{aligned}{\begin{aligned}X&={\frac {\partial F_{11}}{\partial x}}+{\frac {\partial F_{12}}{\partial y}}+{\frac {\partial F_{13}}{\partial z}}+{\frac {\partial F_{14}}{\partial t}}\\Y&={\frac {\partial F_{21}}{\partial x}}+{\frac {\partial F_{22}}{\partial y}}+{\frac {\partial F_{23}}{\partial z}}+{\frac {\partial F_{24}}{\partial t}}\\Z&={\frac {\partial F_{31}}{\partial x}}+{\frac {\partial F_{32}}{\partial y}}+{\frac {\partial F_{33}}{\partial z}}+{\frac {\partial F_{34}}{\partial t}}\\T&={\frac {\partial F_{41}}{\partial x}}+{\frac {\partial F_{42}}{\partial y}}+{\frac {\partial F_{43}}{\partial z}}+{\frac {\partial F_{44}}{\partial t}}\\&=-\left(uX+vY+wZ+Q\right)\\\vartheta &=\varDelta {\sqrt {1-s^{2}}}{\frac {\partial F}{\partial \varepsilon }}=a_{44}\vartheta \end{aligned}}\end{aligned}}}
Then he formulated a third kind of equations of motion by introducing “relative” stresses
S
i
j
{\displaystyle S_{ij}}
into the Euler equations:[ R 17]
X
=
D
t
X
+
∂
S
11
∂
x
+
∂
S
12
∂
y
+
∂
S
13
∂
z
Y
=
D
t
Y
+
∂
S
21
∂
x
+
∂
S
22
∂
y
+
∂
S
23
∂
z
Z
=
D
t
Z
+
∂
S
31
∂
x
+
∂
S
32
∂
y
+
∂
S
33
∂
z
−
T
=
D
t
E
+
∂
(
X
−
u
E
)
∂
x
+
∂
(
Y
−
v
E
)
∂
y
+
∂
(
Z
−
w
E
)
∂
z
[
D
t
(
f
)
=
∂
f
∂
t
+
∂
(
u
f
)
∂
x
+
∂
(
v
f
)
∂
y
+
∂
(
w
f
)
∂
z
]
{\displaystyle {\begin{matrix}{\begin{aligned}{\begin{aligned}X&=D_{t}{\mathfrak {X}}+{\frac {\partial S_{11}}{\partial x}}+{\frac {\partial S_{12}}{\partial y}}+{\frac {\partial S_{13}}{\partial z}}\\Y&=D_{t}{\mathfrak {Y}}+{\frac {\partial S_{21}}{\partial x}}+{\frac {\partial S_{22}}{\partial y}}+{\frac {\partial S_{23}}{\partial z}}\\Z&=D_{t}{\mathfrak {Z}}+{\frac {\partial S_{31}}{\partial x}}+{\frac {\partial S_{32}}{\partial y}}+{\frac {\partial S_{33}}{\partial z}}\\-T&=D_{t}{\mathfrak {E}}+{\frac {\partial \left({\mathfrak {X}}-u{\mathfrak {E}}\right)}{\partial x}}+{\frac {\partial \left({\mathfrak {Y}}-v{\mathfrak {E}}\right)}{\partial y}}+{\frac {\partial \left({\mathfrak {Z}}-w{\mathfrak {E}}\right)}{\partial z}}\end{aligned}}\end{aligned}}\\\left[D_{t}(f)={\frac {\partial f}{\partial t}}+{\frac {\partial (uf)}{\partial x}}+{\frac {\partial (vf)}{\partial y}}+{\frac {\partial (wf)}{\partial z}}\right]\end{matrix}}}
He finally showed how to modify
X
,
Y
,
Z
,
E
,
S
i
j
{\displaystyle {\mathfrak {X,Y,Z,E}},S_{ij}}
using mass density m and pressure p , so that previous equations become the equations of motion of a perfect fluid:[ R 18]
D
t
(
m
u
)
+
∂
p
∂
x
=
X
D
t
(
m
v
)
+
∂
p
∂
y
=
Y
D
t
(
m
w
)
+
∂
p
∂
z
=
Z
−
D
t
(
m
)
+
∂
p
∂
t
=
T
[
m
=
F
+
p
1
−
s
2
,
s
=
u
2
+
v
2
+
w
2
]
{\displaystyle {\begin{matrix}D_{t}(mu)+{\frac {\partial p}{\partial x}}=X\\D_{t}(mv)+{\frac {\partial p}{\partial y}}=Y\\D_{t}(mw)+{\frac {\partial p}{\partial z}}=Z\\-D_{t}(m)+{\frac {\partial p}{\partial t}}=T\\\left[m={\frac {F+p}{1-s^{2}}},\ s={\sqrt {u^{2}+v^{2}+w^{2}}}\right]\end{matrix}}}
which corresponds to the four-force of perfect fluid (e).
Ernst Lamla (1911, published 1912) derived the equation of motion of hydrodynamics independently of Herglotz. Using pressure p and rest density g , he gave the Lagrangian form:[ R 19]
c
⋅
d
d
τ
{
w
(
1
+
P
c
2
)
}
+
l
o
r
P
−
R
=
0
[
R
1
=
X
1
−
q
2
c
2
,
…
R
4
=
i
c
1
−
q
2
c
2
(
X
x
˙
+
Y
y
˙
+
Z
z
˙
)
]
d
d
t
{
x
˙
1
−
q
2
c
2
(
1
+
P
c
2
)
}
+
1
−
q
2
c
2
⋅
δ
P
δ
x
−
X
=
0
d
d
t
{
y
˙
1
−
q
2
c
2
(
1
+
P
c
2
)
}
+
1
−
q
2
c
2
⋅
δ
P
δ
y
−
Y
=
0
d
d
t
{
z
˙
1
−
q
2
c
2
(
1
+
P
c
2
)
}
+
1
−
q
2
c
2
⋅
δ
P
δ
z
−
Z
=
0
d
d
t
{
c
2
1
−
q
2
c
2
(
1
+
P
c
2
)
}
−
1
−
q
2
c
2
δ
P
δ
z
−
(
X
x
˙
+
Y
y
˙
+
Z
z
˙
)
=
0
[
g
=
k
1
−
q
2
c
2
,
P
(
g
)
=
∫
p
=
0
d
p
g
,
l
o
r
P
=
|
∂
P
∂
x
1
,
∂
P
∂
x
2
,
∂
P
∂
x
3
,
∂
P
∂
x
4
|
]
{\displaystyle {\begin{matrix}c\cdot {\frac {d}{d\tau }}\left\{w\left(1+{\frac {P}{c^{2}}}\right)\right\}+\mathrm {lor} \,P-R=0\\\left[R_{1}={\frac {X}{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}},\dots R_{4}={\frac {i}{c{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}}}\left(X{\dot {x}}+Y{\dot {y}}+Z{\dot {z}}\right)\right]\\\hline {\begin{aligned}{\begin{aligned}{\frac {d}{dt}}\left\{{\frac {\dot {x}}{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}}\left(1+{\frac {P}{c^{2}}}\right)\right\}+{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}\cdot {\frac {\delta P}{\delta x}}-X&=0\\{\frac {d}{dt}}\left\{{\frac {\dot {y}}{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}}\left(1+{\frac {P}{c^{2}}}\right)\right\}+{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}\cdot {\frac {\delta P}{\delta y}}-Y&=0\\{\frac {d}{dt}}\left\{{\frac {\dot {z}}{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}}\left(1+{\frac {P}{c^{2}}}\right)\right\}+{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}\cdot {\frac {\delta P}{\delta z}}-Z&=0\\{\frac {d}{dt}}\left\{{\frac {c^{2}}{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}}\left(1+{\frac {P}{c^{2}}}\right)\right\}-{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}{\frac {\delta P}{\delta z}}-\left(X{\dot {x}}+Y{\dot {y}}+Z{\dot {z}}\right)&=0\end{aligned}}\end{aligned}}\\\hline \left[g=k{\sqrt {1-{\frac {q^{2}}{c^{2}}}}},\ P(g)=\int _{p=0}{\frac {dp}{g}},\ \mathrm {lor} \,P=\left|{\frac {\partial P}{\partial x_{1}}},\ {\frac {\partial P}{\partial x_{2}}},\ {\frac {\partial P}{\partial x_{3}}},\ {\frac {\partial P}{\partial x_{4}}}\right|\right]\end{matrix}}}
which corresponds to the four-force of perfect fluid (e), as well as Euler's equation of motion using four-velocity w :[ R 20]
(
w
l
o
r
¯
)
[
w
(
c
2
+
P
)
]
+
l
o
r
P
−
R
=
0
δ
(
u
ε
)
δ
t
+
u
δ
(
u
ε
)
δ
x
+
v
δ
(
u
ε
)
δ
y
+
w
δ
(
u
ε
)
δ
z
+
1
−
q
2
c
2
δ
P
δ
x
−
X
=
0
δ
(
v
ε
)
δ
t
+
u
δ
(
v
ε
)
δ
x
+
v
δ
(
v
ε
)
δ
y
+
w
δ
(
v
ε
)
δ
z
+
1
−
q
2
c
2
δ
P
δ
y
−
Y
=
0
δ
(
w
ε
)
δ
t
+
u
δ
(
w
ε
)
δ
x
+
v
δ
(
w
ε
)
δ
y
+
w
δ
(
w
ε
)
δ
z
+
1
−
q
2
c
2
δ
P
δ
z
−
Z
=
0
c
2
(
δ
ε
δ
t
+
u
δ
ε
δ
x
+
v
δ
ε
δ
y
+
w
δ
ε
δ
z
)
−
1
k
δ
p
δ
t
−
(
X
u
+
Y
v
+
Z
w
)
=
0
(
ε
=
1
+
P
c
2
1
−
q
2
c
2
,
(
w
l
o
r
¯
)
=
w
1
δ
δ
x
1
+
w
2
δ
δ
x
2
+
w
3
δ
δ
x
3
+
w
4
δ
δ
x
4
)
{\displaystyle {\begin{matrix}\left(w{\overline {\mathrm {lor} }}\right)\left[w\left(c^{2}+P\right)\right]+\mathrm {lor} \,P-R=0\\\hline {\begin{aligned}{\frac {\delta (u\varepsilon )}{\delta t}}+u{\frac {\delta (u\varepsilon )}{\delta x}}+v{\frac {\delta (u\varepsilon )}{\delta y}}+w{\frac {\delta (u\varepsilon )}{\delta z}}+{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}{\frac {\delta P}{\delta x}}-X&=0\\{\frac {\delta (v\varepsilon )}{\delta t}}+u{\frac {\delta (v\varepsilon )}{\delta x}}+v{\frac {\delta (v\varepsilon )}{\delta y}}+w{\frac {\delta (v\varepsilon )}{\delta z}}+{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}{\frac {\delta P}{\delta y}}-Y&=0\\{\frac {\delta (w\varepsilon )}{\delta t}}+u{\frac {\delta (w\varepsilon )}{\delta x}}+v{\frac {\delta (w\varepsilon )}{\delta y}}+w{\frac {\delta (w\varepsilon )}{\delta z}}+{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}{\frac {\delta P}{\delta z}}-Z&=0\\c^{2}\left({\frac {\delta \varepsilon }{\delta t}}+u{\frac {\delta \varepsilon }{\delta x}}+v{\frac {\delta \varepsilon }{\delta y}}+w{\frac {\delta \varepsilon }{\delta z}}\right)-{\frac {1}{k}}{\frac {\delta p}{\delta t}}-\left(Xu+Yv+Zw\right)&=0\end{aligned}}\\\hline \left(\varepsilon ={\frac {1+{\frac {P}{c^{2}}}}{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}},\ \left(w{\overline {\mathrm {lor} }}\right)=w_{1}{\frac {\delta }{\delta x_{1}}}+w_{2}{\frac {\delta }{\delta x_{2}}}+w_{3}{\frac {\delta }{\delta x_{3}}}+w_{4}{\frac {\delta }{\delta x_{4}}}\right)\end{matrix}}}
and also gave the Euler equations of motion in the case of substances of least compressibility[ R 21]
l
o
r
P
=
g
(
w
(
w
l
o
r
¯
)
−
1
c
d
w
d
τ
)
+
g
0
c
2
R
(
1
−
q
2
c
2
)
δ
g
δ
x
=
−
g
c
2
{
δ
u
δ
t
+
u
δ
u
δ
x
+
v
δ
u
δ
y
+
w
δ
u
δ
z
−
u
d
i
v
q
}
+
X
g
0
c
2
1
−
q
2
c
2
(
1
−
q
2
c
2
)
δ
g
δ
y
=
−
g
c
2
{
δ
v
δ
t
+
u
δ
v
δ
x
+
v
δ
v
δ
y
+
w
δ
v
δ
z
−
v
d
i
v
q
}
+
Y
g
0
c
2
1
−
q
2
c
2
(
1
−
q
2
c
2
)
δ
g
δ
z
=
−
g
c
2
{
δ
w
δ
t
+
u
δ
w
δ
x
+
v
δ
w
δ
y
+
w
δ
w
δ
z
−
w
d
i
v
q
}
+
Z
g
0
c
2
1
−
q
2
c
2
(
1
−
q
2
c
2
)
δ
g
δ
t
=
−
g
d
i
v
q
−
g
0
c
2
1
−
q
2
c
2
(
X
u
+
Y
v
+
Z
w
)
(
g
=
g
0
1
+
2
p
g
0
c
2
)
{\displaystyle {\begin{matrix}\mathrm {lor} \,P=g\left(w\left(w{\overline {\mathrm {lor} }}\right)-{\frac {1}{c}}{\frac {dw}{d\tau }}\right)+{\frac {g_{0}}{c^{2}}}R\\\hline {\begin{aligned}\left(1-{\frac {q^{2}}{c^{2}}}\right){\frac {\delta g}{\delta x}}&=-{\frac {g}{c^{2}}}\left\{{\frac {\delta u}{\delta t}}+u{\frac {\delta u}{\delta x}}+v{\frac {\delta u}{\delta y}}+w{\frac {\delta u}{\delta z}}-u\,\mathrm {div} \,q\right\}+X{\frac {g_{0}}{c^{2}}}{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}\\\left(1-{\frac {q^{2}}{c^{2}}}\right){\frac {\delta g}{\delta y}}&=-{\frac {g}{c^{2}}}\left\{{\frac {\delta v}{\delta t}}+u{\frac {\delta v}{\delta x}}+v{\frac {\delta v}{\delta y}}+w{\frac {\delta v}{\delta z}}-v\,\mathrm {div} \,q\right\}+Y{\frac {g_{0}}{c^{2}}}{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}\\\left(1-{\frac {q^{2}}{c^{2}}}\right){\frac {\delta g}{\delta z}}&=-{\frac {g}{c^{2}}}\left\{{\frac {\delta w}{\delta t}}+u{\frac {\delta w}{\delta x}}+v{\frac {\delta w}{\delta y}}+w{\frac {\delta w}{\delta z}}-w\,\mathrm {div} \,q\right\}+Z{\frac {g_{0}}{c^{2}}}{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}\\\left(1-{\frac {q^{2}}{c^{2}}}\right){\frac {\delta g}{\delta t}}&=-g\,\mathrm {div} \,q-{\frac {g_{0}}{c^{2}}}{\sqrt {1-{\frac {q^{2}}{c^{2}}}}}\left(Xu+Yv+Zw\right)\end{aligned}}\\\hline \left(g=g_{0}{\sqrt {1+{\frac {2p}{g_{0}c^{2}}}}}\right)\end{matrix}}}
Lewis & Wilson (1912)
edit
w:Edwin Bidwell Wilson and w:Gilbert Newton Lewis devised an alternative 4D vector calculus based on w:Dyadics which, however, never gained widespread support. They defined the dyadic
μ
0
w
w
{\displaystyle \mu _{0}\mathbf {ww} }
using four-velocity
w
{\displaystyle \mathbf {w} }
and rest mass density
μ
0
{\displaystyle \mu _{0}}
in order to formulate the fundamental equation of hydrodynamics:[ R 22]
◊
⋅
(
μ
0
w
w
)
=
0
[
◊
=
k
1
∂
∂
x
1
+
k
2
∂
∂
x
2
+
k
3
∂
∂
x
3
−
k
4
∂
∂
x
4
]
{\displaystyle {\begin{matrix}\lozenge \cdot (\mu _{0}\mathbf {ww} )=0\\\left[\lozenge =\mathbf {k} _{1}{\frac {\partial }{\partial x_{1}}}+\mathbf {k} _{2}{\frac {\partial }{\partial x_{2}}}+\mathbf {k} _{3}{\frac {\partial }{\partial x_{3}}}-\mathbf {k} _{4}{\frac {\partial }{\partial x_{4}}}\right]\end{matrix}}}
equivalent to the (f).
In 1913, in the context of his Entwurf theory (a precursor of general relativity), w:Albert Einstein defined the equation for incoherent matter using the dust stress energy tensor
Θ
μ
ν
{\displaystyle \Theta _{\mu \nu }}
:[ R 23]
∑
ν
∂
T
σ
ν
∂
x
ν
=
1
2
∑
μ
ν
τ
∂
g
μ
ν
x
σ
γ
μ
τ
T
τ
ν
[
T
σ
ν
=
∑
μ
−
g
g
σ
μ
Θ
μ
ν
]
{\displaystyle {\begin{matrix}\sum _{\nu }{\frac {\partial {\mathfrak {T}}_{\sigma \nu }}{\partial x_{\nu }}}={\frac {1}{2}}\sum _{\mu \nu \tau }{\frac {\partial g_{\mu \nu }}{x_{\sigma }}}\gamma _{\mu \tau }{\mathfrak {T}}_{\tau \nu }\\\left[{\mathfrak {T}}_{\sigma \nu }=\sum _{\mu }{\sqrt {-g}}g_{\sigma \mu }\Theta _{\mu \nu }\right]\end{matrix}}}
equivalent to (f), and in 1916 he used the perfect fluid tensor
T
σ
α
{\displaystyle \left.T_{\sigma }\right.^{\alpha }}
:[ R 24]
∂
T
σ
α
∂
x
α
+
1
2
∂
g
μ
ν
x
σ
T
μ
ν
=
0
{\displaystyle {\frac {\partial \left.T_{\sigma }\right.^{\alpha }}{\partial x_{\alpha }}}+{\frac {1}{2}}{\frac {\partial g^{\mu \nu }}{x_{\sigma }}}T_{\mu \nu }=0}
equivalent to (e) in the case of
g
μ
ν
{\displaystyle g_{\mu \nu }}
being the Minkowski tensor.
w:Friedrich Kottler derived the equations of motion from the “Nordström tensor” (i.e. dust tensor) in terms of rest mass density
ν
0
{\displaystyle \nu _{0}}
:[ R 25]
∑
k
=
1
4
∂
∂
x
(
k
)
(
ν
d
x
(
h
)
d
t
d
x
(
k
)
d
t
)
=
ν
d
2
x
(
h
)
d
t
2
{\displaystyle \sum _{k=1}^{4}{\frac {\partial }{\partial x^{(k)}}}\left(\nu {\frac {dx^{(h)}}{dt}}{\frac {dx^{(k)}}{dt}}\right)=\nu {\frac {d^{2}x^{(h)}}{dt^{2}}}}
equivalent to (f), which he then related to the action of a constant external electromagnetic field.