Jonathan Scott said:
I'm just pointing out that there isn't anything available of the appropriate magnitude to replace the pressure integral in the static situation, which is exactly equal in magnitude to the potential energy.
Let's look at the actual Komar mass integral, which is the one that shows the issue we're discussing. For the case of a spherically symmetric, static mass distribution, that integral looks like this [note: edited to correct the factor under the square root]:
$$
M = \int 4 \pi r^2 dr \sqrt{g_{rr} g_{tt}} \left( \rho + 3 p \right)
$$
The exact form of ##g_{tt}## will depend on the exact form of ##\rho## and ##p## as functions of ##r##.
Notice that there is no "potential energy" term. There is, however, the extra factor under the square root sign, which is there because of the variation in "potential" with radius. It turns out that the effect of this factor, which reduces the value of the integral, exactly cancels the effect of adding the pressure, so the final integral is the same as if we just did
$$
M = \int 4 \pi r^2 dr \rho
$$
In other words, we get the same answer as we would if we just "naively" did the integration the way we would in Newtonian physics.
Now, what happens if, by some process, the pressure of this system vanishes? We'll defer (for a bit) the question of how fast it can vanish (it can't do so instantaneously, for the reasons I've given in previous posts). The point is, as long as the process is self-contained, i.e., no energy escapes to infinity, it can't change the value of ##M##. So how does ##M## stay the same if ##p## becomes zero?
The answer is that the above integral is not valid for a non-static system. The factor under the square root is derived from the norm of the timelike Killing vector field of the spacetime, and in a non-static system, such as a star that is collapsing because its internal pressure has been removed by the stoppage of nuclear reactions in its core, there is no timelike Killing vector field. (At the instant that the change begins, there is "almost" a timelike KVF--but at that instant, the pressure can't have instantaneously dropped to zero; what can change in an instant, as I've pointed out before, is that the constraint that is keeping the system static can be removed. But removing that constraint doesn't instantly remove the pressure; it just allows the system to start moving to remove the pressure. The pressure won't actually be zero until there has been enough motion to remove it, by which point the non-staticity will be significant.)
So how can we even evaluate ##M##? If the system is asymptotically flat, we can use the ADM mass (or the Bondi mass, but for the case of no radiation escaping to infinity, they're the same). The ADM mass does not look at the stress-energy tensor; it just looks at the metric coefficients, and it looks at them in a way that is not sensitive to the details of what is changing in the interior of the system; all it is really testing is what the metric looks like far away, to a test object in orbit about the system. By hypothesis, that won't change in the above scenario, and that is why we can say ##M## stays the same even though changes are happening in the interior of the system.
On this view, then, there really isn't a way to compute ##M## by adding up local contributions from all the individual pieces of a system if the system is not static. In other words, on this view, there is no way to say, once the system ceases being static, "where the energy is stored" locally; all you can do is say that ##M##, the total energy as viewed from far away, doesn't change. But many people find that unsatisfactory, so they try to find pseudotensors of one sort or another that capture "energy stored in the gravitational field". The key limitation of all of these, for someone who likes the fact that GR expresses everything in covariant form, is that defining any of these pseudotensors requires picking a particular frame, and only works in that frame. But for the case under discussion, that shouldn't be a serious limitation, since there is already a natural frame to pick: the one in which the system starts out static. If everything is self-contained, that frame should still work as the "center of mass" frame of the system, even when it is no longer static.
I have not seen an analysis of this specific case using the pseudotensor method, but I expect that such an analysis would be interesting.