The problem was with using pm.Deterministic(). I chose a different solution, where the mass is not an observation, but a free variable. Samples of mass are then divided by volume samples to yield density distribution. Those density samples are then used for observation of density directly.
Without pm.Deterministic() this strange problem went away.