When you step through the procedure of deriving a phonon dispersion relationship for a given crystal structure (i.e. small oscillations from equilibrium, harmonic approximation, collective coordinate transformation, normal mode orthogonal transformation), you can rewrite your original interacting Hamiltonian in second quantization language like so, \begin{equation} \mathcal{H}=\sum_{q,\lambda}\hbar\omega_{q,\lambda}(b^\dagger_{q,\lambda}b_{q,\lambda}+\frac{1}{2}), \end{equation} where $q$ and $\lambda$ are the crystal wave number and branch respectively. The above procedure diagonalizes the Hamiltonian, and the space of states that is acted upon by the Hamiltonian operator is spanned by simple harmonic oscillator wave functions with respect to the collective coordinates. i.e. \begin{align} \psi_n(Q)&=(\frac{m\omega}{\pi\hbar})^{\frac{1}{4}}\frac{1}{\sqrt{2^nn!}}H_n(\sqrt{\frac{m\omega}{\hbar}}Q)e^{-\frac{m\omega}{2\hbar}Q^2},\\ \\ Q_q&=\sum_\alpha\hat{\eta} e^{-iq\cdot R_\alpha}, \end{align} meaning we can think of our phonons as decoupled oscillations in our new collective coordinates.
I think my question is pretty simple, but I can't find an answer anywhere: What is the mass of this oscillator?
My naïve guess would be that it has something to do with a derivative of the dispersion relationship similar to the effective mass of Bloch electrons, but I haven't seen a definitive answer.