if you have a Schrödinger operator on a sphere ( $\mathbb{S}^2$) $-\Delta_{\theta,\phi} \psi(\theta,\phi) + V(\theta) \psi(\theta,\phi) = E\psi(\theta,\phi),$ where the potential does not depend on the azimuthal coordinate explicitely, then the substitution $\psi(\theta,\phi) = \Theta(\theta)e^{i n \phi}$ leaves you with the differential equation
$$- \Theta''(\theta)- \cot(\theta) \Theta'(\theta) + \frac{n^2}{\sin^2(\theta)} \Theta(\theta)+ V(\theta) \Theta(\theta) = E \Theta(\theta).$$
For $n=0$ this simplifies to $$- \Theta''(\theta)- \cot(\theta) \Theta'(\theta) + V(\theta) \Theta(\theta) = E \Theta(\theta).$$
Has this equation ever been studied? I mean if $V \in C^{\infty}$, then it should have a discrete spectrum, as the initial PDE has a discrete spectrum. But is there a further literature on this kind of problem?
I mean, by just googling Schrödinger operator on sphere, I got this result for example on arxiv, but if you are aware of better overview articles, please do not hesitate to give them to me. But apparently, the special case, where you can separate that easily is not (as far as I see) that well-studied or am I wrong?