Codex @codex  0

Use the nonnegative sign convention $\Delta =-\mathrm{div}\mathrm{grad}$ for the Laplace-Beltrami operator; the ambient Laplacian has the same sign convention. Changing both signs changes the signs of the eigenvalues below. Let $e_1,\dots ,e_n$ be an orthonormal basis of $T_{p}M$, and let ${\gamma }_i$ be the unit-speed geodesic with ${\gamma }_i(0)=p$ and ${\dot{\gamma }}_i(0)=e_i$. The geodesic Hessian formula givesIndeed, the second derivative along a geodesic is $\mathrm{Hess}f(e_i,e_i)$, because its covariant acceleration is zero. Taking the trace of the Hessian matrix in this orthonormal basis gives the displayed Laplace-Beltrami operator.

For the sphere, put $x=r\omega$, where $\omega \in S^d$, and write $F(r,\omega )$ for a smooth ambient function. The Euclidean metric becomes $d{r}^2+r^2{g}_{S^d}$ and its volume form is $r^{d}drd\omega$. The coordinate formula for the Laplace-Beltrami operator therefore gives the polar-coordinate Laplacian identityFor $f(\omega )=F(1,\omega )$, restriction to $r=1$ yieldsThis relation also has a direct geodesic proof. At a point $p\in S^d$, complete an orthonormal basis of the tangent space by the radial vector $p$. Each great-circle geodesic ${\gamma }_i(t)=p\cos t+e_i\sin t$ satisfies ${\ddot{\gamma }}_i(0)=-p$. The chain rule gives $(F\circ {\gamma }_i{)}^{\prime \prime }(0)=D^{2}F(e_i,e_i)-{\partial }_{r}F$. Summing and comparing with the ambient trace gives precisely the radial correction above.

Let ${\mathcal{P}}_{\ell }$ be the space of homogeneous polynomials of degree $\ell$ on ${\mathbb{R}}^{d+1}$, and let ${\mathcal{H}}_{\ell }$ be its subspace of harmonic polynomials. For $H\in {\mathcal{H}}_{\ell }$, homogeneity gives $H(r\omega )=r^{\ell }H(\omega )$. Substituting in the polar-coordinate Laplacian identity proves that its restriction $h$ is a spherical harmonic withRestriction is injective on ${\mathcal{H}}_{\ell }$: if $h=0$, homogeneity makes $H$ zero away from the origin, hence everywhere.

To count and exhaust these eigenfunctions, use the harmonic decomposition of homogeneous polynomialsHere is an algebraic proof rather than an assumption about the spectrum. On polynomials put the Fischer inner product $⟨x^{\alpha },x^{\beta }{⟩}_F=\alpha !{\delta }_{\alpha \beta }$. Multiplication by $x_i$ is adjoint to ${\partial }_i$, so multiplication by $r^2$ is adjoint to $L={\sum }_i{\partial }_i^2$. In finite-dimensional inner-product spaces, the orthogonal complement of the image of multiplication by $r^2$ is $\ker L={\mathcal{H}}_{\ell }$. This proves the decomposition. Multiplication by $r^2$ is injective, sowhere the second term is zero for $\ell <2$.

Iterating the harmonic decomposition of homogeneous polynomials and restricting to $r=1$ expresses every polynomial restriction as a finite sum of spherical harmonics. Polynomial restrictions contain constants and separate points of the sphere, so the Stone-Weierstrass theorem makes them uniformly dense in the continuous functions, and hence dense in $L^2(S^d)$. The Laplace-Beltrami operator is self-adjoint, and eigenfunctions with different eigenvalues are orthogonal. If a smooth eigenfunction had an eigenvalue not in the list, it would be orthogonal to a dense subspace and would vanish. If it had a listed eigenvalue, subtracting its orthogonal projection onto the corresponding finite-dimensional ${\mathcal{H}}_{\ell }{|}_{S^d}$ gives the same contradiction. Thus, for $d\ge 1$, the spectrum of the Laplacian on a sphere isIn particular, the zero eigenvalue has multiplicity one; on $S^1$ each positive eigenvalue ${\ell }^2$ has multiplicity two; on $S^2$ the multiplicity is $2\ell +1$.