Sorry, I haven't been able to respond for several days.
First, the argument I gave in post #17, although it was in response to
tiny-tim's example in post #11 of a rotating wheel, works perfectly well for any system of particles with arbitrary velocities. I made no assumption about a "rigid body". So my argument also shows that a hot object has more invariant mass than an "identical" cold object (the thermal energy being, of course, kinetic energy of constituent molecules).
The equations I used are the same equations (rearranged and in different notation) as appear in the Wikipedia article
Invariant mass.
What I did do was choose a frame in which the total momentum \textbf{P} = \Sigma \textbf{p} is zero. (As invariant mass is
invariant, I can choose to calculate via any frame I wish.) This frame is called the
centre of momentum frame (a.k.a. centre of mass frame) (I'm British, that's how we spell "centre").
1effect has asked, how do you measure the overall velocity of a system of particles? The answer is, it is the velocity of the centre of momentum frame (relative to the observer asking the question).
For those readers who understand 4-vectors, all the above can be expressed in terms of 4-momentum \left( E/c, \textbf{p} \right). It can be shown this is indeed a 4-vector (you use the Lorentz transform to change inertial frames), and it's additive (by the laws of conservation of energy and momemtum). Invariant mass is just the norm ("length") of this 4-vector, which is automatically an invariant (frame-independent) quantity. To find the centre of momentum frame, calculate the total 4-momentum (the sum of the individual particle 4-momenta), which will be a timelike 4-vector, and choose a coordinate system in which this vector is straight up the "
t" axis.
In post #31,
1effect claimed that the formula (with
c = 1)
M^2 = m_1^2 + m_2^2 + 2\left(E_1 E_2 - \textbf{p}_1 \cdot \textbf{p}_2 \right)
for the invariant mass of a two-particle system was "not invariant" due to the presence of the "variant" term \textbf{p}_1 \cdot \textbf{p}_2. However the term E_1 E_2 isn't invariant either and the "variants" cancel out. In fact the expression E_1 E_2 - \textbf{p}_1 \cdot \textbf{p}_2 is (when
c = 1) the "inner product" of two 4-momentum vectors and is therefore invariant.