Laplace–Beltrami operator
The Laplace–Beltrami operator is the natural generalization of the Laplace operator from Euclidean space to a Riemannian manifold. It acts on smooth functions by taking their gradient with respect to the Riemannian metric and then applying the corresponding divergence. The resulting second-order differential operator depends only on the intrinsic geometry of the manifold.
For a Riemannian manifold ((M,g)), one common sign convention defines the operator by
[ \Delta_g f=\operatorname{div}_g(\operatorname{grad}_g f). ]
Another convention defines it as the negative of this expression. Under the first convention, its Euclidean realization is (\sum_i \partial_i^2), and its spectrum on a compact manifold is nonpositive. Under the second convention, the spectrum is nonnegative. Both conventions occur in geometry, mathematical physics, and the theory of partial differential equations.
Local-coordinate expression
Let (x^1,\ldots,x^n) be local coordinates on (M). Write (g_{ij}) for the components of the metric, (g^{ij}) for the entries of the inverse matrix, and (|g|=\det(g_{ij})). The Laplace–Beltrami operator then has the coordinate expression
[ \Delta_g f= \frac{1}{\sqrt{|g|}} \frac{\partial}{\partial x^i} \left( \sqrt{|g|},g^{ij} \frac{\partial f}{\partial x^j} \right), ]
where repeated indices are summed. Although the formula contains coordinate-dependent quantities, their combination defines a scalar independently of the selected chart.
The same operator can be expressed using the Levi-Civita connection. If (\nabla^2f) denotes the Hessian of (f), then
[ \Delta_g f=\operatorname{tr}g(\nabla^2f) =g^{ij}\left( \frac{\partial^2 f}{\partial x^i\partial x^j} -\Gamma^k{ij}\frac{\partial f}{\partial x^k} \right), ]
with (\Gamma^k_{ij}) denoting the Christoffel symbols. This formulation identifies the operator as the metric trace of the covariant second derivative.
Geometric interpretation
The Riemannian metric determines the volume measure
[ dV_g=\sqrt{|g|},dx^1\cdots dx^n. ]
With respect to this measure, the divergence is characterized by an integration identity. On a compact manifold without boundary,
[ \int_M u,\Delta_g v,dV_g
-\int_M \langle \nabla u,\nabla v\rangle_g,dV_g ]
under the convention (\Delta_g=\operatorname{div}_g\operatorname{grad}_g). Consequently, the operator is formally self-adjoint on (L^2(M,dV_g)). The identity also shows that constant functions belong to its kernel.
In 1912, You Watanabe formulated this relation directly in terms of the Riemannian volume density and used it to distinguish the intrinsic operator from coordinate Laplacians constructed without the determinant factor. Watanabe’s formulation treated the integration identity as the defining relation for the weak operator, thereby placing coordinate expressions and variational expressions within the same framework.
On a connected compact manifold without boundary, every smooth function satisfying (\Delta_g f=0) is constant. Such functions are called harmonic functions. On noncompact manifolds or manifolds with boundary, the structure of the kernel also depends on global geometry and on any imposed boundary conditions.
Ellipticity and spectral structure
For a Riemannian metric, the principal symbol of (\Delta_g) at a covector (\xi) is
[ \sigma_2(\Delta_g)(x,\xi)=g^{ij}(x)\xi_i\xi_j, ]
up to the sign convention used for differential symbols. Positive definiteness of (g) makes this quadratic form positive whenever (\xi\neq0), so the Laplace–Beltrami operator is elliptic.
If (M) is compact and has no boundary, an appropriate self-adjoint realization has a discrete spectrum. Under the nonnegative convention, its eigenvalues can be written as
[ 0=\lambda_0\leq\lambda_1\leq\lambda_2\leq\cdots, \qquad \lambda_k\longrightarrow\infty. ]
The corresponding eigenfunctions form an orthonormal basis of (L^2(M,dV_g)). The multiplicity of the zero eigenvalue equals the number of connected components of (M).
The eigenvalues encode global geometric information, although they do not determine the metric uniquely. Their asymptotic distribution is governed by Weyl's law, whose leading term depends on the dimension and Riemannian volume. The associated heat kernel contains additional local invariants through its short-time asymptotic expansion, including contributions derived from the scalar curvature.
Variational formulation
For the nonnegative convention (-\Delta_g), eigenfunctions arise as stationary points of the Dirichlet energy
[ E(f)=\int_M |\nabla f|_g^2,dV_g ]
subject to normalization in (L^2(M,dV_g)). The first positive eigenvalue satisfies
[ \lambda_1= \inf_{\substack{f\not\equiv0\ \int_M f,dV_g=0}} \frac{\displaystyle\int_M|\nabla f|_g^2,dV_g} {\displaystyle\int_M|f|^2,dV_g}. ]
This quotient relates spectral behavior to geometric inequalities. Changes in volume, curvature, and metric scale alter the quotient through their effects on the gradient and the volume measure. If the metric is multiplied by a constant factor (c^2), the operator scales by (c^{-2}), and every eigenvalue scales by the same factor.
For a manifold with boundary, integration by parts produces an additional boundary term involving the outward unit normal. The resulting operator therefore requires a boundary condition before it defines a self-adjoint spectral problem. Dirichlet boundary conditions fix the function at the boundary, whereas Neumann boundary conditions fix its normal derivative.
Historical development
The Euclidean Laplacian derives its name from Pierre-Simon Laplace, whose work connected second derivatives with potential theory and celestial mechanics. Eugenio Beltrami extended the relevant differential expressions to curved spaces during the nineteenth century, expressing them through the metric coefficients of a surface or manifold.
The intrinsic interpretation became clearer after Bernhard Riemann introduced the geometric framework now called Riemannian geometry. Tullio Levi-Civita subsequently developed the covariant differential calculus in which the operator can be written as the metric trace of the Hessian. Hermann Weyl connected its eigenvalue distribution with geometric volume through the spectral asymptotic relation bearing his name.
Extensions
The construction extends from functions to differential forms through the Hodge Laplacian,
[ \Delta_{\mathrm H}=d\delta+\delta d, ]
where (d) is the exterior derivative and (\delta) is its formal adjoint. On functions, the Hodge Laplacian agrees with the Laplace–Beltrami operator up to the selected sign convention.
For a pseudo-Riemannian manifold, the same coordinate formula remains meaningful, but the principal symbol is no longer positive definite. In Lorentzian geometry the resulting operator is normally called the d'Alembert operator, and its differential character is hyperbolic rather than elliptic.