I have to prove that $[A_j, H] = 0$, with; $$\vec{A} = \frac{1}{2Ze^{2}m}(\vec{L} \times \vec{P} - \vec{p} \times \vec{L}) + \frac{\vec{r}}{r}$$
$$H = \frac{p^2}{2m} - \frac{Ze^2}{r}$$
And, $Z, e, m$ are all constants.
Now, after expanding out the commutator bracket, I'm left with six terms, which I need to prove are all equal to zero.
Now, I can do this for $[\vec{L} \times \vec{p}, p^2]$, $[\vec{p} \times \vec{L}, p^2]$, and $[\frac{\vec{r}}{r}, \frac{1}{r}]$, but I'm having trouble with $[\vec{L} \times \vec{p}, \frac{1}{r}]$, $[\vec{p} \times \vec{L}, \frac{1}{r}]$, and $[\frac{\vec{r}}{r}, p^2]$.
I'm sure once I've $[\vec{L} \times \vec{p}, \frac{1}{r}]$, I can prove $[\vec{p} \times \vec{L}, \frac{1}{r}]$, but, I'm not really sure what to do.
With $[\vec{L} \times \vec{p}, \frac{1}{r}]$, I've done the following; $$[\vec{L} \times \vec{p}, \frac{1}{r}] = \varepsilon_{ijk}[L_{j}p_{k}, \frac{1}{r}]$$ $$ = \varepsilon_{ijk}(L_j[p_k, \frac{1}{r}] + [L_j, \frac{1}{r}]p_k)$$ $$ = \varepsilon_{ijk}L_{j}[p_k, \frac{1}{r}]$$ $$ = \frac{i \varepsilon_{ijk} L_j q_k}{r^3}$$ But from here, I'm not sure where to go. I was hoping to use the cyclic nature of the Levi-Civita symbol, but I can't really see how that helps here. Any input would be fantastic!!