greypilgrim said:
That indeed sounds highly nontrivial. Is QM essential to understand mass, even in classical physics?
My 2 cents: no, this already plays at the classical level.
The action of a point particle with embedding coordinates ##x^0(\tau),x^i(\tau)## and wordlineparameter ##\tau## reads
<br />
S = \frac{m}{2} \int \frac{\dot{x}^i \dot{x}_i}{\dot{x}^0} d \tau<br />
This action is invariant under the infinitesimal Galilei transformations
<br />
\delta_H x^0 = \zeta^0, \ \ \delta_P x^i = \zeta^i, \ \ \delta_J x^i = \lambda^i{}_j x^i, \ \ \delta_G x^i = v^i x^0<br />
From this you can calculate the commutators. In particular, you'll find
<br />
[\delta_P, \delta_G] x^0 = [\delta_P, \delta_G] x^i = 0 \ \ \rightarrow [\delta_P, \delta_G] = 0<br />
But now comes the subtlety: the action is invariant, but the Lagrangian ##L## is only
quasi-invariant under boosts, i.e. it transforms to a total derivative as you can check:
<br />
\delta_G L = \frac{d}{d\tau}\Bigl( m x_i v^i \Bigr)<br />
This has an important consequence: the Noether charge ##Q_G## of the boosts needs to be adjusted to this such that it remains conserved:
<br />
Q_G = p_i v^i x^0 - m x^i v_i<br />
If we now calculate the corresponding Poisson-brackets between the Noether charges ##Q_P = p_i \zeta^i## and ##Q_G##, you will find due to this (signs depend on your conventions for these brackets)
<br />
\{Q_G, Q_P \}_{poisson} \sim m \zeta_i v^i<br />
where ##m## is the mass appearing in the action. This suggests that the Lie algebra of the Galilei group allows for a central extension (i.e. an extra generator ##M## which commutes with all the other generators), playing the role of mass:
<br />
[G,P] = M<br />
You can also check by the Jacobi-identities that this works out. Algebraically, in the non-relativistic limit ##M## replaces the Casimir operator ##P_{\mu}P^{\mu}## of the Poincare algebra which plays the role of mass. This centrally-extended Galilei algebra is called the Bargmann algebra. Interestingly, one can gauge this algebra and show that the corresponding gauge field of the generator ##M## roughly plays the role of the Newton potential. I like this result, because usually people tend to associate central extensions to quantization, but they already play an important role at the classical level.
By the way, this should also answer the question of the OP: in varying the action to obtain the EOM or symmetries, we treat the mass parameter ##m## as a constant. If you want to promote it to a function depending on ##x^0## or ##x^i##, you should also vary it. I'm not sure if this makes sense, though, but it should be clear that the Galilei-symmetries are not apparent anymore.