# Discrete point particles stress energy tensor

I am trying to solve an exercise in Sean Carroll's GR book "Spacetime and Geometry". Basically we need to derive the stress-energy tensor of a perfect fluid (ie $T^{\mu\nu}=(\rho +p)U^{\mu}U^{\nu} + p\eta^{\mu\nu}$) from the stress-energy tensor of a discrete set of particles (ie $T^{\mu\nu}=\sum_a \frac{p^{\mu}_a p^{\nu}_a}{p^0_a}\delta^{(3)}(\mathbf x - \mathbf x^{(a)})$), under the hypothesis of isotropy.

I managed to get in for the $T^{00}$ component and the $T^{0i}$ components, by replacing $p^{\mu}$ by $p^0$ a trivial sum appears: energy density for the 00-component and momentum density the 0i-components (vanishing by isotropy). But I am still struggling with the pure spatial part, I was thinking of substituting the sum by an integral, then the non diagonal part vanishes by isotropy again. Could it be something like: $\sum_a = \int d^3x \rho(x)$? Because then I need a relation that relates the density distribution $\rho$ and $p^{\mu}$ to the pressure (more exactly the definition of the pressure from these two concepts).

There are two points I wish to highlight. 1) Simple substitution of the sum by an integral would not work and is not justified. However, one should switch from microscopical quantities to macroscopical by doing avaraging over 4-volume, throughout which interparticle distances and times can be considered small.

Macroscopic stress-energy tensor will be then: $${\bf T}^{\mu\nu}=\dfrac{1}{\Delta V_4}\int_{\Delta V_4}T^{\mu\nu}d V=\dfrac{1}{\sqrt{-g} d^3x^i dx^0}\int_{\Delta V_4}T^{\mu\nu} \sqrt{-g} d^3x^i dx^0$$ Then a)in $T^{\mu\nu}$ only delta-functions depend on x, b) metric determinant g is a macroscopic quantity, is constant over selected volume and can also be taken away from the integra. One arrives then at: $${\bf T}^{\mu\nu}=\dfrac{1}{d^3x^i dx^0}\sum_a\dfrac{p_a^\mu p_a^\nu}{p_a^0}\int_{\Delta V_4} \delta^{(3)}({\bf x}-{\bf x}^{(a)}) d^3 x^i dx^0 = \dfrac{1}{d^3x^i}\sum_a\dfrac{p_a^\mu p_a^\nu}{p_a^0}.$$ In the last expression the sum is taken over the particles which have world lines passing through $\Delta V_4$ (we ignore the fact that some particles could have left or entered the volume through its 3-boundary, as there are much less of them then the particles inside the volume).

Now the expression ${\bf T}^{\mu\nu}= \dfrac{1}{d^3x^i}\sum_a\dfrac{p_a^\mu p_a^\nu}{p_a^0}$ can be more comfortably treated.

2) The symmetry considerations themselves. Consider the conponenent of the macroscopic tensor:

${\bf T}^{0 0} = \dfrac{1}{d^3x^i}\sum_a p_a^0 \equiv \rho$

${\bf T}^{i 0} = \dfrac{1}{d^3x^i}\sum_a p_a^i$. As the sum $\sum_a p_a^i$ of 3-vectors is taken over a macrospoic volume, the result should result in a macroscopic 3-vector. However, if this vector was not zero, it would violate isotropy, which states that there exists no preferable direction. Hence ${\bf T}^{i 0} = 0$

${\bf T}^{i j} = \dfrac{1}{d^3x^i}\sum_a\dfrac{p_a^i p_a^j}{p_a^0}$. As just before, the sum should produce a symmetric macroscopic 3-tensor of second order. But all symmetric 3-tensors are defined by 3 eigenvectors. If eigenvalues are non-degenerate, then there exists 3 preferred directions (3 eigenvectors), if eigenvalues are single degenerate, then there are 2 preferred directions etc. No preferred directions correspond to the case when the matrix has all eigenvalues equal, that is when the matrix is proportional to kronecker delta. The coefficient of proportionality is the pressure: ${\bf T}^{i j} \equiv P \delta^{ij}$

Expressing $\delta^{ij}$ as $\eta^{ij}+U^0 U^0$ ($U^i$ is zero by symmetry considerations, and $U^0$ is hence equal to unity), and $T^{00}$ as $\rho U^0 U^0$, one arrives at the final expression for $T$.

Here is a more direct approach which is outlined in the solution of Problem 5.2 of Alar P. Lightman's book 'Problem book in relativity and gravitation'. I've slighted edited the book's presentation so that it looks similar with the accepted answer.

We assume the magnitude of velocity $$v$$ is fixed, so the particles have the same magnitude of velocity but they all point in different directions. But the result can be generalized to the case with a non-delta velocity magnitude distribution.

The energy-momentum tensor of point particles can be written as $$T^{\mu\nu} = \sum_{n} \frac{p_n^\mu p_n^\nu}{p_n^0} \delta^3(\vec{x} - \vec{x}_n(t)).$$

Average this over a small 3-volume; using (sum) = (count) $$\times$$ (average), $$T^{\mu\nu} = \frac{1}{\Delta V}\int_{\Delta V}\sum_{n} \frac{p_n^\mu p_n^\nu}{p_n^0} \delta^3(\vec{x} - \vec{x}_n(t)) d^3\vec{x} = \frac{\Delta N}{\Delta V}\left < \frac{p_n^\mu p_n^\nu}{p_n^0} \right >,$$

where $$\Delta N$$ is the number of particles within $$\Delta V$$. Now notice that $$p_n^\mu = \gamma m v_n^\mu$$, where $$v_n^\mu = (1, \vec{v})$$. We have four cases for $$\left < v_n^\mu v_n^\nu \right >$$:

• $$\left < v_n^0 v_n^0 \right > = 1$$.
• $$\left < v_n^0 v_n^i \right > = \left < v_n^i v_n^0 \right > = 0$$.
• $$\displaystyle \left < v_n^i v_n^j \right > = \frac{1}{3} \delta^{ij} v^2$$.

Here, in the last equality, we use that the velocities have isotropically distributed directions: $$\left = \left = \left = \frac{1}{3}\left( \left + \left + \left \right) = \frac{1}{3} v^2.$$

Finally, since the average velocity of the particles is zero in our frame, $$\Delta N/\Delta V = n$$. Hence $$T^{\mu\nu} = \begin{pmatrix} \gamma n m & & & \\ & \gamma n m \frac{1}{3}v^2 & & \\ & & \gamma n m \frac{1}{3}v^2 & \\ & & & \gamma n m \frac{1}{3}v^2 \end{pmatrix}.$$

Comparing this with the perfect fluid's formula we obtain $$\rho = \gamma n m$$ and $$p = \gamma n m \frac{1}{3}v^2$$.