Ingeneral relativity, thevan Stockum dust or theLanczos–van Stockum dust is an exact solution of theEinstein field equations where the gravitational field is generated bydust rotating about an axis of cylindrical symmetry. Since the density of the dust isincreasing with distance from this axis, the solution is rather artificial, but as one of the simplest known solutions in general relativity, it stands as apedagogically important example.

This solution is named afterWillem Jacob van Stockum, who rediscovered it in 1938 independently of a much earlier discovery byCornelius Lanczos in 1924.

One way of obtaining this solution is to look for a cylindrically symmetric perfectfluid solution in which the fluid exhibitsrigid rotation. That is, we demand that the world lines of the fluid particles form a timelike congruence having nonzerovorticity but vanishingexpansion and shear. (In fact, since dust particles feel no forces, this will turn out to be a timelikegeodesic congruence, but this needs not be assumed in advance.)

A simpleansatz corresponding to this demand is expressed by the followingframe field, which contains two undetermined functions of${\displaystyle r}$:

${\displaystyle {\overrightarrow{e}}_0={\mathrm{\partial }}_t,{\overrightarrow{e}}_1=f(r){\mathrm{\partial }}_z,{\overrightarrow{e}}_2=f(r){\mathrm{\partial }}_r,{\overrightarrow{e}}_3=\frac{1}{r}{\mathrm{\partial }}_{\phi }-h(r){\mathrm{\partial }}_t}$

To prevent misunderstanding, we should emphasize that taking thedual coframe

${\displaystyle {\sigma }^0=dt+h(r)rd\phi ,{\sigma }^1=\frac{1}{f(r)}dz,{\sigma }^2=\frac{1}{f(r)}dr,{\sigma }^3=rd\phi }$

gives the metric tensor in terms of the same two undetermined functions:

${\displaystyle g=-{\sigma }^0\otimes {\sigma }^0+{\sigma }^1\otimes {\sigma }^1+{\sigma }^2\otimes {\sigma }^2+{\sigma }^3\otimes {\sigma }^3}$

Multiplying out gives

${\displaystyle d{s}^2=-d{t}^2-2h(r)rdtd\phi +(1-h(r{)}^2)r^{2}d{\phi }^2+\frac{d{z}^2+d{r}^2}{f(r{)}^2}}$

We compute the Einstein tensor with respect to this frame, in terms of the two undetermined functions, and demand that the result have the form appropriate for a perfect fluid solution with the timelike unit vector${\displaystyle {\overrightarrow{e}}_0}$ everywhere tangent to the world line of a fluid particle. That is, we demand that

${\displaystyle G^{\hat{m}\hat{n}}=8\pi \mu \mathrm{diag}(1,0,0,0)+8\pi p\mathrm{diag}(0,1,1,1)}$

This gives the conditions

${\displaystyle f^{\mathrm{\prime }\mathrm{\prime }}=\frac{(f^{\mathrm{\prime }}{)}^2}{f}+\frac{f^{\mathrm{\prime }}}{r},(h^{\mathrm{\prime }}{)}^2+\frac{2h^{\mathrm{\prime }}h}{r}+\frac{h^2}{r^2}=\frac{4f^{\mathrm{\prime }}}{rf}}$

Solving for${\displaystyle f}$ and then for${\displaystyle h}$ gives the desired frame defining the van Stockum solution:

${\displaystyle {\overrightarrow{e}}_0={\mathrm{\partial }}_t,{\overrightarrow{e}}_1=\exp (a^2{r}^2/2){\mathrm{\partial }}_z,{\overrightarrow{e}}_2=\exp (a^2{r}^2/2){\mathrm{\partial }}_r,{\overrightarrow{e}}_3=\frac{1}{r}{\mathrm{\partial }}_{\varphi }-ar{\mathrm{\partial }}_t}$

Note that this frame is only defined on${\displaystyle r>0}$.

Computing the Einstein tensor with respect to our frame shows that in factthe pressure vanishes, so we have adust solution. The mass density of the dust turns out to be

${\displaystyle \mu =\frac{a^2}{2\pi }\exp (a^2{r}^2)}$

Happily, this is finite on the axis of symmetry${\displaystyle r=0}$, but the densityincreases with radius, a feature which unfortunately severely limits possible astrophysical applications.

Solving theKilling equations shows that this spacetime admits a three-dimensionalabelian Lie algebra ofKilling vector fields, generated by

${\displaystyle {\overrightarrow{\xi }}_1={\mathrm{\partial }}_t,{\overrightarrow{\xi }}_2={\mathrm{\partial }}_z,{\overrightarrow{\xi }}_3={\mathrm{\partial }}_{\varphi }}$

Here,${\displaystyle {\overrightarrow{\xi }}_1}$ has nonzero vorticity, so we have astationary spacetime invariant under translation along the world lines of the dust particles, and also under translation along the axis of cylindrical symmetry and rotation about that axis.

Note that unlike theGödel dust solution, in the van Stockum dust the dust particles are rotating about ageometrically distinguished axis.

As promised, the expansion and shear of the timelike geodesic congruence${\displaystyle {\overrightarrow{e}}_0}$ vanishes, but the vorticity vector is

${\displaystyle \overrightarrow{\mathrm{\Omega }}=-a\exp (a^2{r}^2/2){\overrightarrow{e}}_1}$

This means that even though in our comoving chart the world lines of the dust particles appear as vertical lines, in fact they are twisting about one another as the dust particles swirl about the axis of symmetry. In other words, if we follow the evolution of a small ball of dust, we find that it rotates about its own axis (parallel to${\displaystyle r=0}$), but does not shear or expand; the latter properties define what we mean byrigid rotation. Notice that on the axis itself, the magnitude of the vorticity vector becomes simply${\displaystyle a}$.

The tidal tensor is

${\displaystyle E_{\hat{m}\hat{n}}=a^2\exp (a^2{r}^2)\mathrm{diag}(0,1,1)}$

which shows that observers riding on the dust particles experience isotropic tidal tension in the plane of rotation. The magnetogravitic tensor is

Consider athought experiment in which an observer riding on a dust particle sitting on the axis of symmetry looks out at dust particles with positive radial coordinate. Does he see them to berotating, or not?

Since the top array of null geodesics is obtained simply by translating upwards the lower array, and since the three world lines are all vertical (invariant undertime translation), it might seem that the answer is "no". However, while the frame given above is aninertial frame, computing thecovariant derivatives

${\displaystyle {\mathrm{\nabla }}_{{\overrightarrow{e}}_0}{\overrightarrow{e}}_1,{\mathrm{\nabla }}_{{\overrightarrow{e}}_0}{\overrightarrow{e}}_2,{\mathrm{\nabla }}_{{\overrightarrow{e}}_0}{\overrightarrow{e}}_3}$

shows that only the first vanishes identically. In other words, the remaining spatial vectors arespinning about${\displaystyle {\overrightarrow{e}}_1}$ (i.e. about an axis parallel to the axis of cylindrical symmetry of this spacetime).

Thus, to obtain anonspinning inertial frame we need to spin up our original frame, like this:

${\displaystyle {\overrightarrow{f}}_0={\overrightarrow{e}}_0,{\overrightarrow{f}}_1={\overrightarrow{e}}_1,{\overrightarrow{f}}_2=\cos (\theta ){\overrightarrow{e}}_2+\sin (\theta ){\overrightarrow{e}}_3,{\overrightarrow{f}}_3=-\sin (\theta ){\overrightarrow{e}}_2+\cos (\theta ){\overrightarrow{e}}_3}$

where${\displaystyle \theta =tq(r)}$ where q is a new undetermined function of r. Plugging in the requirement that the covariant derivatives vanish, we obtain

${\displaystyle \theta =at\exp (a^2{r}^2/2)}$

The new frame appears, in our comoving coordinate chart, to be spinning, but in fact it is gyrostabilized. In particular, since our observer with the green world line in the figure is presumably riding anonspinning dust particle (otherwisespin-spin forces would be apparent in the dynamics of the dust), he in fact observes nearby radially separated dust particles to be rotating clockwise about his location with angular velocity a. This explains the physical meaning of the parameter which we found in our earlier derivation of the first frame.

(Pedantic note: alert readers will have noticed that we ignored the fact that neither of our frame fields is well defined on the axis. However, we can define a frame for an on-axis observer by an appropriate one-sided limit; this gives a discontinuous frame field, but we only need to define a framealong the world line of our on-axis observer in order to pursue the thought experiment considered in this section.)

It is worth remarking that the null geodesicsspiral inwards in the above figure. This means that our on-axis observer sees the other dust particles attime-lagged locations, which is of course just what we would expect. The fact that the null geodesics appear "bent" in this chart is of course an artifact of our choice ofcomoving coordinates in which the world lines of the dust particles appear as vertical coordinate lines.

Let us draw the light cones for some typical events in the van Stockum dust, to see how their appearance (in our comoving cylindrical chart) depends on the radial coordinate:

As the figure^[which?] shows, at${\displaystyle r=a^{-1}}$, the cones become tangent to the coordinate plane${\displaystyle t=t_0}$, and we obtain a closed null curve (the red circle). Note that this isnot a null geodesic.

As we move further outward, we can see that horizontal circles with larger radii areclosed timelike curves. The paradoxical nature of these CTCs was apparently first pointed out by van Stockum: observers whose world lines form a closed timelike curve can apparently revisit or affect their own past. Even worse, there is apparently nothing to prevent such an observer from deciding, on his third lifetime, say, to stop accelerating, which would give him multiple biographies.

These closed timelike curves arenot timelike geodesics, so these paradoxical observers mustaccelerate to experience these effects. Indeed, as we would expect, the required accelerationdiverges as these timelike circles approach the null circles lying in the critical cylinder${\displaystyle r=a^{-1}}$.

Closed timelike curves turn out to exist in many other exact solutions in general relativity, and their common appearance is one of the most troubling theoretical objections to this theory. However, very few physicists refuse to use general relativity at all on the basis of such objections; rather most take the pragmatic attitude that using general relativity makes sense whenever one can get away with it, because of the relative simplicity and well established reliability of this theory in many astrophysical situations. This is not unlike the fact that many physicists use Newtonian mechanics every day, even though they are well aware that Galilean kinematics has been "overthrown" by relativistic kinematics.

• Dust solution
• Gödel dust solution

• Exact solutions in general relativity
• Thought experiments in physics

• Articles with short description
• Short description is different from Wikidata
• Pages using sidebar with the child parameter
• All articles with specifically marked weasel-worded phrases
• Articles with specifically marked weasel-worded phrases from May 2022