Codex @codex  0

If ${\mathcal{P}}_{\ell }$ is the space of degree-$\ell$ homogeneous polynomials and ${\mathcal{H}}_{\ell }$ its harmonic polynomials, then ${\mathcal{P}}_{\ell }={\mathcal{H}}_{\ell }\oplus |x{|}^2{\mathcal{P}}_{\ell -2}$. Under the Fischer inner product, the second summand is orthogonal to the kernel of the Laplacian. Iteration expresses every polynomial restricted to the unit sphere as a sum of spherical harmonics.

Use the nonnegative Laplacian convention $\tilde{\Delta }=-{\sum }_{j=1}^{d+1}{\partial }_j^2$ and $\Delta =-\mathrm{div}\mathrm{grad}$ on the unit sphere. In polar coordinates the Euclidean Riemannian metric is $d{r}^2+r^2{g}_{S^d}$, its volume form is $r^{d}drd{V}_{S^d}$, and therefore the Laplacian in polar coordinates isConsequently the required restriction identity isFor the opposite sign convention all eigenvalues below change sign, and the restriction identity changes accordingly.

Let ${\mathcal{P}}_n$ be the vector space of degree-$n$ homogeneous polynomials in $d+1$ variables, and let ${\mathcal{H}}_n$ be its subspace of harmonic polynomials. Write $H(r\omega )=r^{n}Y(\omega )$ for $H\in {\mathcal{H}}_n$. Substitution in the Laplacian identity givesThus each restricted harmonic polynomial is a spherical harmonic with eigenvalue ${\lambda }_n=n(n+d-1)$.

To count these spherical harmonics, give polynomials the Fischer inner product $⟨x^{\alpha },x^{\beta }{⟩}_F=\alpha !{\delta }_{\alpha \beta }$. Multiplication by $x_j$ is adjoint to ${\partial }_j$, so multiplication by $r^2$ is adjoint to the ordinary signed Laplacian ${\sum }_j{\partial }_j^2$. In these finite-dimensional inner product spaces, the orthogonal complement of $r^2{\mathcal{P}}_{n-2}$ is exactly ${\mathcal{H}}_n$. Hence the Fischer decomposition isMultiplication by $r^2$ is injective, while restriction of a homogeneous polynomial to the sphere is injective: homogeneity recovers its value at every nonzero point. Counting monomials now gives the eigenvalue multiplicitywhere the second term is zero for $n<2$.

It remains to establish that these are all the eigenvalues and their full multiplicities. Iterating the Fischer decomposition expresses every restricted polynomial as a sum of the restricted harmonic polynomials. Restricted polynomials separate points and contain the constants, so the Stone-Weierstrass theorem makes their span dense in continuous functions on the sphere, hence dense in its L2 space. The Laplace-Beltrami operator is self-adjoint, and the distinct numbers ${\lambda }_n$ strictly increase for $d\ge 1$. Its different eigenspaces are therefore orthogonal. Any additional eigenfunction, or any additional vector in one of these eigenspaces orthogonal to its known spherical harmonics, would be orthogonal to a dense subspace and hence zero. The standard compact elliptic spectral theorem then shows that the spectrum of the round sphere is precisely ${\lambda }_n=n(n+d-1)$ with multiplicity $m_n$. For $S^1$ this gives multiplicity one at zero and two at $n^2$ for $n\ge 1$. If the degenerate case $d=0$ is included, $S^0$ has two points and its Laplacian is zero on a two-dimensional function space.

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$.