Estoy en una clase de mecánica cuántica, y se da en el libro que los operadores $\hat{L^{2}}$ y $\hat{H}$ para el Oscilador Armónico 3D, pero no se da ninguna prueba matemática definitiva, y yo mismo estoy teniendo dificultades para demostrarlo, y conceptualizar por qué debe ser cierto.
He intentado utilizar la coordenada esférica para demostrarlo, y sé que en coordenadas esféricas $$\hat{L^{2}}= \frac{-h^{2}}{4\pi^{2}}(\frac{1}{\sin\theta}\frac{d}{d\theta}(\sin\theta \frac{d}{d\theta})+\frac{1}{\sin^{2}\theta}\frac{d^{2}}{d\phi^{2}})$$
Y $$\hat{H}=\frac{-h^{2}}{4\pi^{2}(2m)}\Delta+\frac{1}{2}kr^{2}.$$
He estado tratando de probarlo usando el muy básico $[\hat{L^{2}},\hat{H}]f= \hat{L^{2}}\hat{H}f-\hat{H}\hat{L^{2}}f$ método para mostrar la relación de conmutación. Mi primer pensamiento fue que al aplicar $\hat{L^{2}}$ a $\hat{H}$ todos los términos que dependen de r o $\frac{r}{dr}$ desaparecería, pero si alguna función arbitraria f tuviera términos cruzados, esto no es necesariamente cierto, y el álgebra se vuelve bastante confusa después. ¿Hay una forma mejor de demostrar esto?