The energy momentum tensor is a symmetric, rank-2 contravariant tensor $T$ which determines the right-hand side of the Einstein field equations.

If FieldType = "DiracWeyl", then the additional arguments for EnergyMomentumTensor are: a solder form (compatible with the metric $g$), a rank 1 covariant spinor $\mathrm{\ψ}$, and the complex conjugate $\overline{\mathrm{\ψ}}$.

If FieldType = "Dust", then the additional arguments for EnergyMomentumTensor are: a vector field $U$, a scalar $\mathrm{\μ}$ (energy density).

If FieldType = "Electromagnetic", then the additional arguments are either: a 1-form $A$ (the electromagnetic 4-potential), or a skew-symmetric rank 2 tensor $F$ (the field strength tensor).

If FieldType = "PerfectFluid", then the additional arguments for EnergyMomentumTensor are: a vector field $U$, and scalars $\mathrm{\μ}$ (energy density) and $p$ (pressure).

If FieldType = "Scalar", then the additional argument for EnergyMomentumTensor is a scalar $\mathrm{\φ}$.

If FieldType = "NMCScalar", then the additional argument for EnergyMomentumTensor is a non-minimally coupled scalar $\mathrm{\φ}$.

See the Details help page for the explicit formulas used to calculate the various energy-momentum tensors, the matter field equations and the divergence identities.

These commands are part of the DifferentialGeometry:-Tensor: package, and so can be used in the form EnergyMomentumTensor(...), MatterFieldEquations(...), DivergenceIdentities(...) only after executing the commands with(DifferentialGeometry); with(Tensor); in that order. They can always be used in the long form DifferentialGeometry:-Tensor:-EnergyMomentumTensor, DifferentialGeometry:-Tensor-MatterFieldEquations, DifferentialGeometry:-Tensor:-DivergenceIdentities.

$\mathrm{with}\left(\mathrm{DifferentialGeometry}\right)\:$$\mathrm{with}\left(\mathrm{Tools}\right)\:$$\mathrm{with}\left(\mathrm{Tensor}\right)\:$

$\mathrm{DGsetup}\left(\left[t\,x\,y\,z\right]\,\left[\mathrm{z1}\,\mathrm{z2}\,\mathrm{w1}\,\mathrm{w2}\right]\,N\right)$

$\mathrm{frame\; name:\; N}$

$\mathrm{g1}≔\mathrm{evalDG}\left(x^4\mathrm{dt}\&t\mathrm{dt}-\mathrm{dx}\&t\mathrm{dx}-\mathrm{dy}\&t\mathrm{dy}-\mathrm{dz}\&t\mathrm{dz}\right)$

$\mathrm{g1}≔x^4\mathrm{dt}\mathrm{dt}-\mathrm{dx}\mathrm{dx}-\mathrm{dy}\mathrm{dy}-\mathrm{dz}\mathrm{dz}$

$\mathrm{OTetrad}≔\mathrm{evalDG}\left(\left[\frac{1}{x^2}\mathrm{D\_t}\,\mathrm{D\_x}\,\mathrm{D\_y}\,\mathrm{D\_z}\right]\right)$

$\mathrm{OTetrad}≔\left[\frac{1}{x^2}\mathrm{D\_t}\,\mathrm{D\_x}\,\mathrm{D\_y}\,\mathrm{D\_z}\right]$

$\mathrm{\σ1}≔\mathrm{SolderForm}\left(\mathrm{OTetrad}\right)$

$\mathrm{\σ1}≔\frac{x^2\sqrt{2}}{2}\mathrm{dt}\mathrm{D\_z1}\mathrm{D\_w1}+\frac{x^2\sqrt{2}}{2}\mathrm{dt}\mathrm{D\_z2}\mathrm{D\_w2}+\frac{\sqrt{2}}{2}\mathrm{dx}\mathrm{D\_z1}\mathrm{D\_w2}+\frac{\sqrt{2}}{2}\mathrm{dx}\mathrm{D\_z2}\mathrm{D\_w1}-\frac{\mathrm{I}}{2}\sqrt{2}\mathrm{dy}\mathrm{D\_z1}\mathrm{D\_w2}+\frac{\mathrm{I}}{2}\sqrt{2}\mathrm{dy}\mathrm{D\_z2}\mathrm{D\_w1}+\frac{\sqrt{2}}{2}\mathrm{dz}\mathrm{D\_z1}\mathrm{D\_w1}-\frac{\sqrt{2}}{2}\mathrm{dz}\mathrm{D\_z2}\mathrm{D\_w2}$

$\mathrm{\ψ1}≔\mathrm{evalDG}\left(h\left(x\right)\mathrm{dz1}-f\left(x\right)\mathrm{dz2}\right)$

$\mathrm{\ψ1}≔h\left(x\right)\mathrm{dz1}-f\left(x\right)\mathrm{dz2}$

$\mathrm{barpsi1}≔\mathrm{evalDG}\left(h\left(x\right)\mathrm{dw1}-f\left(x\right)\mathrm{dw2}\right)$

$\mathrm{barpsi1}≔h\left(x\right)\mathrm{dw1}-f\left(x\right)\mathrm{dw2}$

$\mathrm{T1}≔\mathrm{EnergyMomentumTensor}\left(''DiracWeyl''\,\mathrm{g1}\,\mathrm{\σ1}\,\mathrm{\ψ1}\,\mathrm{barpsi1}\right)$

$\mathrm{T1}≔-\frac{\sqrt{2}\left({f\left(x\right)}^2-{h\left(x\right)}^2\right)}{x^3}\mathrm{D\_t}\mathrm{D\_y}+\sqrt{2}\left(-h\left(x\right)f\prime \left(x\right)+f\left(x\right)h\prime \left(x\right)\right)\mathrm{D\_x}\mathrm{D\_y}-\frac{\sqrt{2}\left({f\left(x\right)}^2-{h\left(x\right)}^2\right)}{x^3}\mathrm{D\_y}\mathrm{D\_t}+\sqrt{2}\left(-h\left(x\right)f\prime \left(x\right)+f\left(x\right)h\prime \left(x\right)\right)\mathrm{D\_y}\mathrm{D\_x}$

$\mathrm{E1}≔\mathrm{MatterFieldEquations}\left(''DiracWeyl''\,\mathrm{g1}\,\mathrm{\σ1}\,\mathrm{\ψ1}\,\mathrm{barpsi1}\right)$

$\mathrm{E1}≔\frac{\frac{\mathrm{I}}{2}\sqrt{2}\left(f\prime \left(x\right)x+f\left(x\right)\right)}{x}\mathrm{D\_w1}-\frac{\frac{\mathrm{I}}{2}\sqrt{2}\left(h\prime \left(x\right)x+h\left(x\right)\right)}{x}\mathrm{D\_w2},-\frac{\frac{\mathrm{I}}{2}\sqrt{2}\left(f\prime \left(x\right)x+f\left(x\right)\right)}{x}\mathrm{D\_z1}+\frac{\frac{\mathrm{I}}{2}\sqrt{2}\left(h\prime \left(x\right)x+h\left(x\right)\right)}{x}\mathrm{D\_z2}$

$\mathrm{Div1},\mathrm{RHS1}≔\mathrm{DivergenceIdentities}\left(''DiracWeyl''\,\mathrm{g1}\,\mathrm{\σ1}\,\mathrm{\ψ1}\,\mathrm{barpsi1}\,\mathrm{T1}\,\mathrm{E1}\right)$

$\mathrm{Div1},\mathrm{RHS1}≔\frac{\sqrt{2}\left(-xh\left(x\right)f″\left(x\right)+xf\left(x\right)h″\left(x\right)+2f\left(x\right)h\prime \left(x\right)-2h\left(x\right)f\prime \left(x\right)\right)}{x}\mathrm{D\_y},\frac{\sqrt{2}\left(-xh\left(x\right)f″\left(x\right)+xf\left(x\right)h″\left(x\right)+2f\left(x\right)h\prime \left(x\right)-2h\left(x\right)f\prime \left(x\right)\right)}{x}\mathrm{D\_y}$

$\mathrm{Div1}-\mathrm{RHS1}$

$0$

$\mathrm{map}\left(\mathrm{DGsimplify}\,\mathrm{eval}\left(\left[\mathrm{E1}\right]\,\left[f\left(x\right)=\frac{1}{x}\,h\left(x\right)=\frac{1}{x}\right]\right)\right)$

$\left[0\mathrm{D\_z1}\,0\mathrm{D\_z1}\right]$

$\mathrm{DGsimplify}\left(\mathrm{eval}\left(\mathrm{Div1}\,\left[f\left(x\right)=\frac{1}{x}\,h\left(x\right)=\frac{1}{x}\right]\right)\right)$

$0\mathrm{D\_t}$

$\mathrm{DGsetup}\left(\left[t\,x\,y\,z\right]\,M\right)$

$\mathrm{frame\; name:\; M}$

$\mathrm{g2}≔\mathrm{evalDG}\left(\mathrm{dt}\&t\mathrm{dt}-t^2\mathrm{dx}\&t\mathrm{dx}-\mathrm{dy}\&t\mathrm{dy}-\mathrm{dz}\&t\mathrm{dz}\right)$

$\mathrm{g2}≔\mathrm{dt}\mathrm{dt}-t^2\mathrm{dx}\mathrm{dx}-\mathrm{dy}\mathrm{dy}-\mathrm{dz}\mathrm{dz}$

$\mathrm{u2}≔\mathrm{evalDG}\left(\cosh \left(f\left(t\right)\right)\mathrm{D\_t}-\frac{\sinh \left(f\left(t\right)\right)}{t}\mathrm{D\_x}\right)$

$\mathrm{u2}≔\cosh \left(f\left(t\right)\right)\mathrm{D\_t}-\frac{\sinh \left(f\left(t\right)\right)}{t}\mathrm{D\_x}$

$\mathrm{TensorInnerProduct}\left(\mathrm{g2}\,\mathrm{u2}\,\mathrm{u2}\right)$

$1$

$\mathrm{\μ2}≔h\left(t\right)$

$\mathrm{\μ2}≔h\left(t\right)$

$\mathrm{T2}≔\mathrm{EnergyMomentumTensor}\left(''Dust''\,\mathrm{g2}\,\mathrm{u2}\,\mathrm{\μ2}\right)$

$\mathrm{T2}≔h\left(t\right){\cosh \left(f\left(t\right)\right)}^2\mathrm{D\_t}\mathrm{D\_t}-\frac{h\left(t\right)\cosh \left(f\left(t\right)\right)\sinh \left(f\left(t\right)\right)}{t}\mathrm{D\_t}\mathrm{D\_x}-\frac{h\left(t\right)\cosh \left(f\left(t\right)\right)\sinh \left(f\left(t\right)\right)}{t}\mathrm{D\_x}\mathrm{D\_t}+\frac{h\left(t\right)\left({\cosh \left(f\left(t\right)\right)}^2-1\right)}{t^2}\mathrm{D\_x}\mathrm{D\_x}$

$\mathrm{E2}≔\mathrm{MatterFieldEquations}\left(''Dust''\,\mathrm{g2}\,\mathrm{u2}\,\mathrm{\μ2}\right)$

$\mathrm{E2}≔\frac{h\left(t\right)\cosh \left(f\left(t\right)\right)+\stackrel{\mathbf{.}}{h}\left(t\right)\cosh \left(f\left(t\right)\right)t+h\left(t\right)\sinh \left(f\left(t\right)\right)\stackrel{\mathbf{.}}{f}\left(t\right)t}{t},\frac{{\cosh \left(f\left(t\right)\right)}^2-1+\cosh \left(f\left(t\right)\right)\sinh \left(f\left(t\right)\right)\stackrel{\mathbf{.}}{f}\left(t\right)t}{t}\mathrm{D\_t}-\frac{\cosh \left(f\left(t\right)\right)\left(\sinh \left(f\left(t\right)\right)+\cosh \left(f\left(t\right)\right)\stackrel{\mathbf{.}}{f}\left(t\right)t\right)}{t^2}\mathrm{D\_x}$

$\mathrm{Soln}≔\left[h\left(t\right)=\frac{\mathrm{\_C2}}{{\left(1+t^2{\mathrm{\_C1}}^2\right)}^{\frac{1}{2}}}\,f\left(t\right)=\mathrm{arcsinh}\left(\frac{1}{t\mathrm{\_C1}}\right)\right]$

$\mathrm{Soln}≔\left[h\left(t\right)=\frac{\mathrm{\_C2}}{\sqrt{1+t^2{\mathrm{\_C1}}^2}}\,f\left(t\right)=\mathrm{arcsinh}\left(\frac{1}{t\mathrm{\_C1}}\right)\right]$

$\mathrm{simplify}\left(\mathrm{eval}\left(\left[\mathrm{E2}\right]\,\mathrm{Soln}\right)\,\mathrm{symbolic}\right)$

$\left[0\,0\mathrm{D\_t}+0\mathrm{D\_x}\right]$

$\mathrm{Div2},\mathrm{RHS2}≔\mathrm{DivergenceIdentities}\left(''Dust''\,\mathrm{g2}\,\mathrm{u2}\,\mathrm{\μ2}\,\mathrm{T2}\,\mathrm{E2}\right)$

$\mathrm{Div2},\mathrm{RHS2}≔\frac{2h\left(t\right){\cosh \left(f\left(t\right)\right)}^2-h\left(t\right)+\stackrel{\mathbf{.}}{h}\left(t\right)t{\cosh \left(f\left(t\right)\right)}^2+2h\left(t\right)\cosh \left(f\left(t\right)\right)\sinh \left(f\left(t\right)\right)\stackrel{\mathbf{.}}{f}\left(t\right)t}{t}\mathrm{D\_t}-\frac{2h\left(t\right)\cosh \left(f\left(t\right)\right)\sinh \left(f\left(t\right)\right)+\stackrel{\mathbf{.}}{h}\left(t\right)\cosh \left(f\left(t\right)\right)\sinh \left(f\left(t\right)\right)t+2h\left(t\right){\cosh \left(f\left(t\right)\right)}^2\stackrel{\mathbf{.}}{f}\left(t\right)t-h\left(t\right)\stackrel{\mathbf{.}}{f}\left(t\right)t}{t^2}\mathrm{D\_x},\frac{2h\left(t\right){\cosh \left(f\left(t\right)\right)}^2-h\left(t\right)+\stackrel{\mathbf{.}}{h}\left(t\right)t{\cosh \left(f\left(t\right)\right)}^2+2h\left(t\right)\cosh \left(f\left(t\right)\right)\sinh \left(f\left(t\right)\right)\stackrel{\mathbf{.}}{f}\left(t\right)t}{t}\mathrm{D\_t}-\frac{2h\left(t\right)\cosh \left(f\left(t\right)\right)\sinh \left(f\left(t\right)\right)+\stackrel{\mathbf{.}}{h}\left(t\right)\cosh \left(f\left(t\right)\right)\sinh \left(f\left(t\right)\right)t+2h\left(t\right){\cosh \left(f\left(t\right)\right)}^2\stackrel{\mathbf{.}}{f}\left(t\right)t-h\left(t\right)\stackrel{\mathbf{.}}{f}\left(t\right)t}{t^2}\mathrm{D\_x}$

$\mathrm{Div2}\&minus\mathrm{RHS2}$

$0\mathrm{D\_t}$

$\mathrm{DGsetup}\left(\left[t\,x\,y\,z\right]\,M\right)$

$\mathrm{frame\; name:\; M}$

$\mathrm{g3}≔\mathrm{evalDG}\left(x^2\mathrm{dt}\&t\mathrm{dt}-\mathrm{dx}\&t\mathrm{dx}-\mathrm{dy}\&t\mathrm{dy}-\mathrm{dz}\&t\mathrm{dz}\right)$

$\mathrm{g3}≔x^2\mathrm{dt}\mathrm{dt}-\mathrm{dx}\mathrm{dx}-\mathrm{dy}\mathrm{dy}-\mathrm{dz}\mathrm{dz}$

$\mathrm{A3}≔\mathrm{evalDG}\left(\mathrm{f1}\left(x\right)\mathrm{dt}+\mathrm{f2}\left(x\right)\mathrm{dy}\right)$

$\mathrm{A3}≔\mathrm{f1}\left(x\right)\mathrm{dt}+\mathrm{f2}\left(x\right)\mathrm{dy}$

$\mathrm{T3}≔\mathrm{EnergyMomentumTensor}\left(''Electromagnetic''\,\mathrm{g3}\,\mathrm{A3}\right)$

$\mathrm{T3}≔-\frac{{\mathrm{f2}\prime \left(x\right)}^2{x}^2+{\mathrm{f1}\prime \left(x\right)}^2}{2x^4}\mathrm{D\_t}\mathrm{D\_t}+\frac{\mathrm{f1}\prime \left(x\right)\mathrm{f2}\prime \left(x\right)}{x^2}\mathrm{D\_t}\mathrm{D\_y}+\frac{-{\mathrm{f2}\prime \left(x\right)}^2{x}^2+{\mathrm{f1}\prime \left(x\right)}^2}{2x^2}\mathrm{D\_x}\mathrm{D\_x}+\frac{\mathrm{f1}\prime \left(x\right)\mathrm{f2}\prime \left(x\right)}{x^2}\mathrm{D\_y}\mathrm{D\_t}-\frac{{\mathrm{f2}\prime \left(x\right)}^2{x}^2+{\mathrm{f1}\prime \left(x\right)}^2}{2x^2}\mathrm{D\_y}\mathrm{D\_y}-\frac{-{\mathrm{f2}\prime \left(x\right)}^2{x}^2+{\mathrm{f1}\prime \left(x\right)}^2}{2x^2}\mathrm{D\_z}\mathrm{D\_z}$

$\mathrm{F3}≔\mathrm{ExteriorDerivative}\left(\mathrm{A3}\right)$

$\mathrm{F3}≔-\mathrm{f1}\prime \left(x\right)\mathrm{dt}\bigwedge \mathrm{dx}+\mathrm{f2}\prime \left(x\right)\mathrm{dx}\bigwedge \mathrm{dy}$

$\mathrm{EnergyMomentumTensor}\left(''Electromagnetic''\,\mathrm{g3}\,\mathrm{F3}\right)$

$-\frac{{\mathrm{f2}\prime \left(x\right)}^2{x}^2+{\mathrm{f1}\prime \left(x\right)}^2}{2x^4}\mathrm{D\_t}\mathrm{D\_t}+\frac{\mathrm{f1}\prime \left(x\right)\mathrm{f2}\prime \left(x\right)}{x^2}\mathrm{D\_t}\mathrm{D\_y}+\frac{-{\mathrm{f2}\prime \left(x\right)}^2{x}^2+{\mathrm{f1}\prime \left(x\right)}^2}{2x^2}\mathrm{D\_x}\mathrm{D\_x}+\frac{\mathrm{f1}\prime \left(x\right)\mathrm{f2}\prime \left(x\right)}{x^2}\mathrm{D\_y}\mathrm{D\_t}-\frac{{\mathrm{f2}\prime \left(x\right)}^2{x}^2+{\mathrm{f1}\prime \left(x\right)}^2}{2x^2}\mathrm{D\_y}\mathrm{D\_y}-\frac{-{\mathrm{f2}\prime \left(x\right)}^2{x}^2+{\mathrm{f1}\prime \left(x\right)}^2}{2x^2}\mathrm{D\_z}\mathrm{D\_z}$

$\mathrm{E3}≔\mathrm{MatterFieldEquations}\left(''Electromagnetic''\,\mathrm{g3}\,\mathrm{A3}\right)$

$\mathrm{E3}≔-\frac{\mathrm{f1}\prime \left(x\right)-\mathrm{f1}″\left(x\right)x}{x^3}\mathrm{D\_t}-\frac{\mathrm{f2}″\left(x\right)x+\mathrm{f2}\prime \left(x\right)}{x}\mathrm{D\_y},0\mathrm{dt}\bigwedge \mathrm{dx}\bigwedge \mathrm{dy}$

$\mathrm{MatterFieldEquations}\left(''Electromagnetic''\,\mathrm{g3}\,\mathrm{F3}\right)$

$-\frac{\mathrm{f1}\prime \left(x\right)-\mathrm{f1}″\left(x\right)x}{x^3}\mathrm{D\_t}-\frac{\mathrm{f2}″\left(x\right)x+\mathrm{f2}\prime \left(x\right)}{x}\mathrm{D\_y},0\mathrm{dt}\bigwedge \mathrm{dx}\bigwedge \mathrm{dy}$

$\mathrm{Div3},\mathrm{RHS3}≔\mathrm{DivergenceIdentities}\left(''Electromagnetic''\,\mathrm{g3}\,\mathrm{A3}\,\mathrm{T3}\,\mathrm{E3}\left[1\right]\right)$

$\mathrm{Div3},\mathrm{RHS3}≔-\frac{\mathrm{f2}\prime \left(x\right)x^3\mathrm{f2}″\left(x\right)-x\mathrm{f1}\prime \left(x\right)\mathrm{f1}″\left(x\right)+{\mathrm{f1}\prime \left(x\right)}^2+{\mathrm{f2}\prime \left(x\right)}^2{x}^2}{x^3}\mathrm{D\_x},-\frac{\mathrm{f2}\prime \left(x\right)x^3\mathrm{f2}″\left(x\right)-x\mathrm{f1}\prime \left(x\right)\mathrm{f1}″\left(x\right)+{\mathrm{f1}\prime \left(x\right)}^2+{\mathrm{f2}\prime \left(x\right)}^2{x}^2}{x^3}\mathrm{D\_x}$

$\mathrm{Div3}\&minus\mathrm{RHS3}$

$0\mathrm{D\_t}$

$\mathrm{DGsimplify}\left(\mathrm{eval}\left(\mathrm{E3}\,\left[\mathrm{f1}\left(x\right)=x^2\,\mathrm{f2}\left(x\right)=\ln \left(x\right)\right]\right)\right)$

$0\mathrm{D\_t}$

$\mathrm{DGsimplify}\left(\mathrm{eval}\left(\mathrm{Div3}\,\left[\mathrm{f1}\left(x\right)=x^2\,\mathrm{f2}\left(x\right)=\ln \left(x\right)\right]\right)\right)$

$0\mathrm{D\_t}$

$\mathrm{DGsetup}\left(\left[t\,x\,y\,z\right]\,M\right)$

$\mathrm{frame\; name:\; M}$

$\mathrm{g4}≔\mathrm{evalDG}\left(\mathrm{dt}\&t\mathrm{dt}-t^2\mathrm{dx}\&t\mathrm{dx}-\mathrm{dy}\&t\mathrm{dy}-\mathrm{dz}\&t\mathrm{dz}\right)$

$\mathrm{g4}≔\mathrm{dt}\mathrm{dt}-t^2\mathrm{dx}\mathrm{dx}-\mathrm{dy}\mathrm{dy}-\mathrm{dz}\mathrm{dz}$

$\mathrm{u4}≔\mathrm{evalDG}\left(2\mathrm{D\_t}+\frac{\mathrm{sqrt}\left(3\right)}{t}\mathrm{D\_x}\right)$

$\mathrm{u4}≔2\mathrm{D\_t}+\frac{\sqrt{3}}{t}\mathrm{D\_x}$

$\mathrm{TensorInnerProduct}\left(\mathrm{g4}\,\mathrm{u4}\,\mathrm{u4}\right)$

$1$

$\mathrm{\μ4}≔k\left(t\right)$

$\mathrm{\μ4}≔k\left(t\right)$

$\mathrm{p4}≔h\left(t\right)$

$\mathrm{p4}≔h\left(t\right)$

$\mathrm{T4}≔\mathrm{EnergyMomentumTensor}\left(''PerfectFluid''\,\mathrm{g4}\,\mathrm{u4}\,\mathrm{\μ4}\,\mathrm{p4}\right)$

$\mathrm{T4}≔\left(5h\left(t\right)+4k\left(t\right)\right)\mathrm{D\_t}\mathrm{D\_t}+\frac{2\left(k\left(t\right)+h\left(t\right)\right)\sqrt{3}}{t}\mathrm{D\_t}\mathrm{D\_x}+\frac{2\left(k\left(t\right)+h\left(t\right)\right)\sqrt{3}}{t}\mathrm{D\_x}\mathrm{D\_t}+\frac{2h\left(t\right)+3k\left(t\right)}{t^2}\mathrm{D\_x}\mathrm{D\_x}-h\left(t\right)\mathrm{D\_y}\mathrm{D\_y}-h\left(t\right)\mathrm{D\_z}\mathrm{D\_z}$

$\mathrm{E4}≔\mathrm{MatterFieldEquations}\left(''PerfectFluid''\,\mathrm{g4}\,\mathrm{u4}\,\mathrm{\μ4}\,\mathrm{p4}\right)$

$\mathrm{E4}≔\frac{7k\left(t\right)+7h\left(t\right)+5\stackrel{\mathbf{.}}{h}\left(t\right)t+4t\stackrel{\mathbf{.}}{k}\left(t\right)}{t}\mathrm{D\_t}+\frac{2\sqrt{3}\left(2k\left(t\right)+2h\left(t\right)+t\stackrel{\mathbf{.}}{k}\left(t\right)+\stackrel{\mathbf{.}}{h}\left(t\right)t\right)}{t^2}\mathrm{D\_x}$

$\mathrm{de}≔\mathrm{DGinfo}\left(\mathrm{E4}\,''CoefficientSet''\right)$

$\mathrm{de}≔\left\{\frac{7k\left(t\right)+7h\left(t\right)+5\stackrel{\mathbf{.}}{h}\left(t\right)t+4t\stackrel{\mathbf{.}}{k}\left(t\right)}{t}\,\frac{2\sqrt{3}\left(2k\left(t\right)+2h\left(t\right)+t\stackrel{\mathbf{.}}{k}\left(t\right)+\stackrel{\mathbf{.}}{h}\left(t\right)t\right)}{t^2}\right\}$

$\mathrm{dsolve}\left(\mathrm{de}\right)$

$\left\{h\left(t\right)=\mathrm{\_C1}+\frac{\mathrm{\_C2}}{t^2}\,k\left(t\right)=-\mathrm{\_C1}-\frac{3\mathrm{\_C2}}{t^2}\right\}$

$\mathrm{DGsetup}\left(\left[t\,x\,y\,z\right]\,M\right)$

$\mathrm{frame\; name:\; M}$

$\mathrm{g5}≔\mathrm{evalDG}\left(\mathrm{dt}\&t\mathrm{dt}-t^2\mathrm{dx}\&t\mathrm{dx}-\mathrm{dy}\&t\mathrm{dy}-\mathrm{dz}\&t\mathrm{dz}\right)$

$\mathrm{g5}≔\mathrm{dt}\mathrm{dt}-t^2\mathrm{dx}\mathrm{dx}-\mathrm{dy}\mathrm{dy}-\mathrm{dz}\mathrm{dz}$

$\mathrm{\φ5}≔f\left(t\right)$

$\mathrm{\φ5}≔f\left(t\right)$

$\mathrm{T5}≔\mathrm{EnergyMomentumTensor}\left(''Scalar''\,\mathrm{g5}\,\mathrm{\φ5}\right)$

$\mathrm{T5}≔-\left(-\frac{{\stackrel{\mathbf{.}}{f}\left(t\right)}^2}{2}+\frac{{\mathrm{\_m}}^2{f\left(t\right)}^2}{2}\right)\mathrm{D\_t}\mathrm{D\_t}+\frac{{\stackrel{\mathbf{.}}{f}\left(t\right)}^2+{\mathrm{\_m}}^2{f\left(t\right)}^2}{2t^2}\mathrm{D\_x}\mathrm{D\_x}+\left(\frac{{\stackrel{\mathbf{.}}{f}\left(t\right)}^2}{2}+\frac{{\mathrm{\_m}}^2{f\left(t\right)}^2}{2}\right)\mathrm{D\_y}\mathrm{D\_y}+\left(\frac{{\stackrel{\mathbf{.}}{f}\left(t\right)}^2}{2}+\frac{{\mathrm{\_m}}^2{f\left(t\right)}^2}{2}\right)\mathrm{D\_z}\mathrm{D\_z}$

$\mathrm{E5}≔\mathrm{MatterFieldEquations}\left(''Scalar''\,\mathrm{g5}\,\mathrm{\φ5}\right)$

$\mathrm{E5}≔\frac{\stackrel{\mathbf{.}}{f}\left(t\right)+\stackrel{\mathbf{..}}{f}\left(t\right)t}{t}-{\mathrm{\_m}}^{2}f\left(t\right)$

$\mathrm{Div5},\mathrm{RHS5}≔\mathrm{DivergenceIdentities}\left(''Scalar''\,\mathrm{g5}\,\mathrm{\φ5}\,\mathrm{T5}\,\mathrm{E5}\right)$

$\mathrm{Div5},\mathrm{RHS5}≔\frac{\stackrel{\mathbf{.}}{f}\left(t\right)\left(\stackrel{\mathbf{.}}{f}\left(t\right)+\stackrel{\mathbf{..}}{f}\left(t\right)t-t{\mathrm{\_m}}^{2}f\left(t\right)\right)}{t}\mathrm{D\_t},\frac{\stackrel{\mathbf{.}}{f}\left(t\right)\left(\stackrel{\mathbf{.}}{f}\left(t\right)+\stackrel{\mathbf{..}}{f}\left(t\right)t-t{\mathrm{\_m}}^{2}f\left(t\right)\right)}{t}\mathrm{D\_t}$

$\mathrm{Div5}\&minus\mathrm{RHS5}$

$0\mathrm{D\_t}$

$\mathrm{DGsetup}\left(\left[t\,x\,y\,z\right]\,M\right)$

$\mathrm{frame\; name:\; M}$

$\mathrm{g6}≔\mathrm{evalDG}\left(\mathrm{dt}\&t\mathrm{dt}-t^2\mathrm{dx}\&t\mathrm{dx}-\mathrm{dy}\&t\mathrm{dy}-\mathrm{dz}\&t\mathrm{dz}\right)$

$\mathrm{g6}≔\mathrm{dt}\mathrm{dt}-t^2\mathrm{dx}\mathrm{dx}-\mathrm{dy}\mathrm{dy}-\mathrm{dz}\mathrm{dz}$

$\mathrm{\φ6}≔f\left(t\right)$

$\mathrm{\φ6}≔f\left(t\right)$

$\mathrm{T6}≔\mathrm{EnergyMomentumTensor}\left(''NMCScalar''\,\mathrm{g6}\,\mathrm{\φ6}\right)$

$\mathrm{T6}≔\frac{-{\mathrm{\_m}}^2{f\left(t\right)}^{2}t+4\mathrm{\_\ξ}f\left(t\right)\stackrel{\mathbf{.}}{f}\left(t\right)+{\stackrel{\mathbf{.}}{f}\left(t\right)}^{2}t}{2t}\mathrm{D\_t}\mathrm{D\_t}-\frac{-{\mathrm{\_m}}^2{f\left(t\right)}^2+4\mathrm{\_\ξ}f\left(t\right)\stackrel{\mathbf{..}}{f}\left(t\right)+4{\stackrel{\mathbf{.}}{f}\left(t\right)}^2\mathrm{\_\ξ}-{\stackrel{\mathbf{.}}{f}\left(t\right)}^2}{2t^2}\mathrm{D\_x}\mathrm{D\_x}-\frac{-{\mathrm{\_m}}^2{f\left(t\right)}^{2}t+4\mathrm{\_\ξ}f\left(t\right)\stackrel{\mathbf{.}}{f}\left(t\right)+4\mathrm{\_\ξ}f\left(t\right)\stackrel{\mathbf{..}}{f}\left(t\right)t+4{\stackrel{\mathbf{.}}{f}\left(t\right)}^{2}t\mathrm{\_\ξ}-{\stackrel{\mathbf{.}}{f}\left(t\right)}^{2}t}{2t}\mathrm{D\_y}\mathrm{D\_y}-\frac{-{\mathrm{\_m}}^2{f\left(t\right)}^{2}t+4\mathrm{\_\ξ}f\left(t\right)\stackrel{\mathbf{.}}{f}\left(t\right)+4\mathrm{\_\ξ}f\left(t\right)\stackrel{\mathbf{..}}{f}\left(t\right)t+4{\stackrel{\mathbf{.}}{f}\left(t\right)}^{2}t\mathrm{\_\ξ}-{\stackrel{\mathbf{.}}{f}\left(t\right)}^{2}t}{2t}\mathrm{D\_z}\mathrm{D\_z}$

$\mathrm{E6}≔\mathrm{MatterFieldEquations}\left(''NMCScalar''\,\mathrm{g6}\,\mathrm{\φ6}\right)$

$\mathrm{E6}≔\frac{\stackrel{\mathbf{.}}{f}\left(t\right)+\stackrel{\mathbf{..}}{f}\left(t\right)t}{t}-{\mathrm{\_m}}^{2}f\left(t\right)$

$\mathrm{Div6},\mathrm{RHS6}≔\mathrm{DivergenceIdentities}\left(''Scalar''\,\mathrm{g6}\,\mathrm{\φ6}\,\mathrm{T6}\,\mathrm{E6}\right)$

$\mathrm{Div6},\mathrm{RHS6}≔\frac{\stackrel{\mathbf{.}}{f}\left(t\right)\left(\stackrel{\mathbf{.}}{f}\left(t\right)+\stackrel{\mathbf{..}}{f}\left(t\right)t-t{\mathrm{\_m}}^{2}f\left(t\right)\right)}{t}\mathrm{D\_t},\frac{\stackrel{\mathbf{.}}{f}\left(t\right)\left(\stackrel{\mathbf{.}}{f}\left(t\right)+\stackrel{\mathbf{..}}{f}\left(t\right)t-t{\mathrm{\_m}}^{2}f\left(t\right)\right)}{t}\mathrm{D\_t}$

$\mathrm{Div6}\&minus\mathrm{RHS6}$

$0\mathrm{D\_t}$