PAllen said:
[edit: perhaps all that’s needed is a statement that an accelerating particle or body is not a closed or isolated system; you must include the sources of acceleration]
That's indeed a very important point. Many misunderstandings (at least I had them, when I started to learn relativity) comes from not carefully keeping in mind that defining the total energy and momentum of extended systems is difficult. The only exception are indeed closed systems.
Math is, as always, your friend, and it's easy to see, where this problem comes from. First of all the natural description of relativistic physics are local laws for fields (classical and quantum, but I'll stay here within the realm of classical physics), i.e., you usually deal with local field equations of motion like, e.g., classical electrodynamics, where you describe matter naturally in terms of continuum mechanics, i.e., in terms of the energy-momentum tensor and electric charge and current densities.
Now, if you want to consider extended objects, you have to somehow integrate in a meaningful way the corresponding densities, well defined as relativistic field quantities. The mathematical difficulty is that naively taking spatial integrals in a given (inertial) reference frames in general do not lead to tensor quantities, but only when you have a conserved quantity in a closed system.
E.g. if you have electric charges and currents, it's described by the four-current density ##j^{\mu}##, which fulfills (from the Maxwell equations alone due to Noether's 2nd theorem/gauge invariance) the continuity equation,
$$\partial_{\mu} j^{\mu}=0.$$
Then the total charge of the system is given by
$$Q=\int_{\mathbb{R}^3} \mathrm{d}^3 r j^0(t,\vec{r}).$$
The equation of continuity immediately guarantees that this is a conserved quantity, i.e.,
$$\dot{Q}=\int_{\mathbb{R}^3} \mathrm{d}^3 r \partial_t j^0(t,\vec{r})=-\int_{\mathbb{R}^3} \mathrm{d}^3 r \vec{\nabla} \cdot \vec{j}(t,\vec{r}=0,$$
because for realistic situation the current goes to 0 at infinity sufficiently large so that the surface integral at infinity from Gauss's theorem vanishes.
Nothing in the above calculation, however ensures that ##Q## is a scalar. But that's also easy to prove from the
equation of continuity! Just use the 4D Gauss's theorem. As integration four-volume take a "cylindrical" volume with two time-like bottom and top non-overlapping hyper-surfaces and let the space-like "rim" of the cylinder go to infinity. Then the continuity equation shows
$$\int_{V^{(4)}} \mathrm{d}^4 x \partial_{\mu} j^{\mu}(t,\vec{x})=\int_{\partial V^{(4)}} \mathrm{d}^3 \sigma_{\mu} j^{\mu}(t,\vec{x})=0.$$
This shows that the integral over the top and bottom time-like hypersurface exactly cancel. As one of the hypersurfaces you can choose ##t=\text{t}_0=\text{const}## of the "lab frame" of the observer, which shows that the quantity which you get from integrating over either the top or the bottom is simply ##Q##, and this value doesn't depend on the chosen time-like hypersurface. That's why ##Q## is a Lorentz scalar, and it's only a Lorentz scalar, because the equation of continuity holds.
In the same way it's clear that energy and momentum are field-theoretically described by the (symmetric) energy-momentum tensor ##T^{\mu \nu}##, and the total momentum
$$P^{\mu}=\int_{\mathbb{R}^3} \mathrm{d}^3 r T^{\mu 0}(t,\vec{r})$$
is in general only a four-vector, if ##\partial_{\nu} T^{\mu \nu}=0##, i.e., if you consider the total energy-momentum tensor of a closed system. Then the total energy and momentum defined by the non-covariantly written integral is nevertheless a four-vector, and it's conserved: ##\dot{P}^{\mu}=0##.