We use , whose spectrum on a compact manifold is nonpositive; write for the positive Laplace-Beltrami operator. This convention will also make the heat operator in Question 3 equal to .
Let be an orthonormal basis of , and let be the unit-speed geodesic with and . Since a geodesic has vanishing acceleration for the Levi-Civita connection, differentiation twice along it gives the Riemannian Hessian:
Taking its metric trace yields the geodesic trace formula for the Laplace-Beltrami operator
The positive-sign convention puts a minus sign in front of this sum.
For the unit sphere, take and tangent unit vectors orthogonal to . Its great circle geodesics are . For a smooth ambient function near the sphere, the ordinary chain rule gives
The ambient orthonormal basis is . Thus, writing and similarly for ,
This proves the ambient restriction formula for the spherical Laplacian and shows why ambient derivatives normal to the sphere must be subtracted. The equivalent formula in spherical coordinates is
Let be the vector space of homogeneous polynomials of degree in variables, and its harmonic polynomials. If , homogeneous polynomial scaling gives . Substitution into the spherical coordinates formula gives
The restrictions are the degree- spherical harmonics. A homogeneous polynomial vanishing on the sphere vanishes on every nonzero ray and hence identically, so restriction is injective on .
We prove the harmonic decomposition of homogeneous polynomials that both counts these eigenfunctions and establishes completeness. For an ambient harmonic homogeneous polynomial of degree and , the spherical coordinates formula gives
For its coefficient is nonzero. Inductively assume homogeneous polynomials of degree have been decomposed into sums of . Apply this decomposition to for , and lift each summand by multiplying by and dividing by the displayed nonzero coefficient. The resulting satisfies , so . The same coefficient formula shows that no nonzero element of can be harmonic. Starting at degrees zero and one therefore proves the direct decomposition
Since , the multiplicity is
where the second term is zero for .
On , iterating the decomposition writes every polynomial restriction as a finite sum of spherical harmonics. The Stone-Weierstrass theorem makes these polynomial restrictions dense in , since coordinate functions separate points and constants are included; they are consequently dense in . Integration by parts makes symmetric on smooth functions and gives orthogonality between distinct eigenvalues. Normalizing within each finite-dimensional eigenspace gives a complete orthonormal eigenbasis . To check the full operator spectrum, let . For smooth , integration by parts gives the coefficients of as . Therefore finite harmonic partial sums approximate both and in , hence approximate in the graph norm. Conversely every coefficient sequence with is the graph norm limit of its finite harmonic sums. Thus the closure of has exactly this operator domain and is the real diagonal self-adjoint operator with entries . If is not one of these entries, the inverse has entries and is bounded because . This proves that there are no missing spectral values. The spectrum of the Laplacian on a sphere is
with the above multiplicities; for , negate every eigenvalue. For example has eigenvalues , with multiplicities one at zero and two for , and has multiplicity . The zero-dimensional sphere consists of two isolated points and has only the zero eigenvalue, with multiplicity two.