 Spherical inversion preserves angular derivatives

For $\overline{\mathbf{r}}=a^2\mathbf{r}/r^2$, the chain rule gives ${\nabla }_{\mathbf{r}}=(a^2/r^2)(I-2\hat{\mathbf{r}}{\hat{\mathbf{r}}}^T){\nabla }_{\overline{\mathbf{r}}}$. The radial term is killed by crossing with $\mathbf{r}$, giving the identity. Purely angular vector spherical harmonics are therefore unchanged as angular basis functions. This is distinct from pushing forward a physical vector field by the inversion Jacobian.