I seem to be having trouble understanding how Wigner D-matrices rotate spherical harmonics. I asked this question on the Maths Stack Exchange but decided to cast my net a bit wider and ask the question here too.
Suppose that we want to understand the effect of a rotation $R$ on the normal (spin-zero) spherical harmonic $Y^l_m(\vec{n})$, defined by a unit direction vector $\vec{n}$. Every textbook (and Wikipedia) I've consulted thus far has claimed that the rotated harmonic $Y^l_m(R^{-1}\vec{n})$ can be expressed as a linear combination of harmonics of the same order $l$ $$Y^l_m(R^{-1}\vec{n}) = \sum_{m'}D^{(l)}_{mm'}[R]Y^{l}_{m'}(\vec{n}),$$ where $$D^{(l)}_{mm'}[R] \equiv \langle lm|R|lm'\rangle$$ are the Wigner D-matrices and $|lm\rangle$ are the usual eigen-vectors of the angular momentum operators $J_z$ and $\vec{J}^2$.
Now, I do not understand why the sum in the above equation only runs over $m'$. $l$ is associated with the inclination angle $\theta$ and $m$ is associated with the azimuthal angle $\phi$, both of which can be affected by a general rotation. I suppose another way of phrasing this is: why is the $l$ index of the spherical harmonic left unaffected by a general rotation?