---
title: Spectral Convergence of Graph Laplacians
url: https://www.emergentmind.com/topics/spectral-convergence-of-graph-laplacians
type: topic
---

# Spectral Convergence of Graph Laplacians

Searching arXiv for recent and foundational papers on spectral convergence of graph Laplacians to ground the article in the literature.
Spectral convergence of graph Laplacians is the study of how spectral data of discrete graph operators—eigenvalues, eigenvectors, eigenspaces, resolvents, spectral projections, heat semigroups, or even empirical spectral distributions—approach the corresponding objects of a limiting operator under graph refinement, random sampling, bandwidth shrinkage, or structural degeneration. In the literature represented here, the limiting object may be the Laplace–Beltrami operator on a manifold, a weighted Laplacian under non-uniform sampling, a connection Laplacian on a vector bundle, a Laplacian on a fractal, an effective Laplacian on a coarsened graph, or another discrete graph Laplacian related by perturbation or contraction [1507.00151][1704.00064][1306.1587][2601.18057][2412.15937].

## 1. Operator settings and model classes

The subject encompasses several distinct approximation regimes. One major regime begins with i.i.d. samples from a compact manifold and forms a weighted graph using an $\epsilon$-graph, a $k$-NN graph, a Gaussian kernel, or a symmetrized normalization. Another begins with an already discrete graph and varies its potential, self-loops, or intra-cluster weights. A third treats graphs as approximants of non-manifold spaces such as post-critically finite fractals. A fourth studies vector-valued analogues, where the scalar Laplacian is replaced by a connection Laplacian acting on sections of an associated vector bundle [1910.13476][2110.06988][1704.00064][1306.1587].

| Setting | Discrete object | Limiting or comparison operator |
|---|---|---|
| Random samples on manifolds | $\epsilon$-graphs, $k$-NN graphs, Gaussian-kernel graph Laplacians | Laplace–Beltrami or weighted Laplacian [1910.13476][2101.09875] |
| Manifolds with boundary | Symmetrized or truncated graph Laplacian | Neumann or Dirichlet Laplace–Beltrami operator [2110.06988] |
| Vector-valued data on bundles | Connection graph Laplacian | Connection Laplacian on an associated bundle [1306.1587] |
| pcf fractals | Finite-dimensional weighted graph Laplacians | Fractal Laplacian via quasi-unitary equivalence [1704.00064] |
| General discrete graphs | $L_0$ and $L_c$ on $\ell^2(X,m)$ | Spectral comparison within the discrete category [2412.15937] |
| Tightly connected clusters | $L_\beta$ with diverging intra-cluster weights | Effective Laplacian on a coarsened graph [2601.18057] |

This breadth is important because “spectral convergence of graph Laplacians” does not designate a single theorem. It designates a family of approximation theories whose hypotheses, normalizations, and limit operators differ substantially. In particular, the limit need not be the ordinary Laplace–Beltrami operator: non-uniform sampling produces weighted operators, bundle structure produces connection Laplacians, and purely discrete perturbations can yield exact comparison identities rather than continuum limits [1507.00151][1306.1587][2412.15937].

## 2. Modes of convergence and analytic mechanisms

The strongest results in this literature are not restricted to eigenvalue-by-eigenvalue convergence. For pcf fractals, the graph energy forms and fractal energy form are $\delta_m$-quasi-unitarily equivalent, with $\delta_m \to 0$ exponentially fast, yielding norm resolvent convergence,
\[
\left\| (\Delta + 1)^{-1} J - J (\Delta_m + 1)^{-1} \right\| \leq \delta_m,
\]
as well as operator-norm convergence of heat operators, spectral projections, and other functions of the Laplacian; eigenfunctions converge in energy norm [1704.00064].

A second mechanism is resolvent convergence under graph contraction. When intra-cluster weights are scaled by a parameter $\beta \to \infty$, the corresponding Laplacians satisfy
\[
(L_\beta - zI)^{-1} \longrightarrow J^\uparrow (\underline{L} - zI)^{-1} J^\downarrow
\]
for $z \in \mathbb{C}\setminus[0,\infty)$, with strong resolvent convergence in the undirected setting and norm resolvent convergence for finite clusters; semigroup convergence follows as well [2601.18057].

A third mechanism is direct spectral comparison inside the discrete category. For general and possibly infinite discrete graphs with purely discrete spectrum, the perturbation identity
\[
\sum_{n=1}^\infty \left(\lambda_n(c)-\lambda_n(0)\right)=\sum_{x\in X} c(x)m(x)
\]
holds, with the right-hand side equal to infinity if $c\notin \ell^1(X,m)$. The same work proves the discrete local Weyl law
\[
\sum_{n=1}^\infty |f_n(x)|^2=\frac{1}{m(x)},
\]
using Mercer’s theorem and the heat kernel diagonal, without Tauberian arguments or the continuous Weyl law [2412.15937].

In periodic metric graphs, convergence and comparison are organized through spectral band bracketing. Spectral bands of the metric Laplacian are localized by Dirichlet and Neumann eigenvalues on a finite fundamental domain, and the passage between discrete and metric spectra is mediated by the Cattaneo correspondence, which maps discrete band edges through $-\cos z = \lambda$ [1406.7523].

These results show that spectral convergence is operator-theoretic before it is geometric. Depending on the setting, the decisive tool may be quasi-unitary equivalence, min–max comparison, RKHS methods, optimal transport, heat-kernel interpolation, quadratic-form monotonicity, or direct discrete heat-kernel expansions [1704.00064][2110.06988][1910.13476][2101.09875][2601.18057][2412.15937].

## 3. Quantitative convergence from random samples

The best-developed quantitative theory concerns low-lying eigenpairs of graph Laplacians built from i.i.d. samples on a manifold. A common structure is the balance of bias from kernel smoothing against variance from empirical approximation, with graph bandwidth chosen as a function of sample size and intrinsic dimension.

| Setting | Representative quantitative statement | Reference |
|---|---|---|
| Point Integral Method, non-uniform density | With probability at least $1-1/n$, \(\left| \lambda_i^{t,n} - \lambda_i \right| \leq C_1 \left( t^{1/2} + \frac{\log n + |\log t| + 1}{ t^{k+3}\sqrt{n} } \right)\) | [1507.00151] |
| Random geometric graphs | \(\frac{| \lambda_k(\Gamma)-\lambda_k(M)|}{\lambda_k(M)} = O\!\left(\sqrt{\frac{\log n^{p_m}}{n^{1/m}}}\right)\) for fixed low-lying spectrum | [1801.10108] |
| $\epsilon$-graphs and $k$-NN graphs | Optimal rate, up to log factors, is \(O(n^{-1/(m+4)})\) for both eigenvalues and eigenvectors | [1910.13476] |
| Gaussian kernelized graph Laplacian | If \(\epsilon \sim (\log N/N)^{1/(d/2+2)}\), eigenvalue rate is \(N^{-1/(d/2+2)}\) and eigenvector rate is \(N^{-1/(d+4)}\); if \(\epsilon \sim (\log N/N)^{1/(d/2+3)}\), both rates are \(N^{-1/(d/2+3)}\) | [2101.09875] |
| $L^\infty$ spectral convergence | Eigenvalue error \(O((\log n/n)^{3/(8d+26)})\) and eigenfunction sup-norm error \(O((\log n/n)^{1/(8d+16)})\), together with a heat-kernel reconstruction bound | [1912.05680] |

For the symmetric normalized graph Laplacian on a compact submanifold, one estimate is
\[
|\lambda_{n,h}-\lambda_h| \leq \frac{C'}{\sqrt{n}h^{5d+3}}, \qquad
\|a_nu_{n,h}-u_h\| \leq \frac{C''}{\sqrt{n}h^{4d+3}},
\]
and when the data are sampled from a tubular neighborhood rather than the manifold itself, the required scaling changes from intrinsic dimension $d$ to ambient dimension $D$; the accompanying numerical study indicates the necessity of a denoising step before applying spectral algorithms [1510.08110].

The rates are not uniform across constructions. They depend on normalization, kernel regularity, bandwidth scaling, the choice between $\epsilon$-graphs and $k$-NN graphs, the target norm for eigenvector comparison, and whether the claim concerns pointwise, $L^2$, $H^1$, or $L^\infty$ convergence. What is consistent across the literature is that low-lying spectral data converge under explicit sampling and bandwidth regimes, and that these regimes are controlled by intrinsic geometric dimension rather than ambient dimension in the noise-free setting [1507.00151][1801.10108][1910.13476][2101.09875][1912.05680][1510.08110].

## 4. Density, normalization, boundary conditions, and cluster structure

Non-uniform sampling changes the continuum limit. In the Point Integral Method, as $n\to\infty$ and $t\to 0$, the discrete spectrum converges to the weighted Neumann problem
\[
-\frac{1}{p(x)^2}\,\mathrm{div}\!\left(p(x)^2\nabla u(x)\right)=\lambda u(x),
\]
which reduces to the Laplace–Beltrami operator only when $p$ is uniform [1507.00151]. For Gaussian kernels, the density-corrected Laplacian is designed precisely so that, under non-uniform sampling, it recovers the Laplace–Beltrami operator rather than a Fokker–Planck-type limit [2101.09875].

Boundary effects are also spectrally visible. On manifolds with boundary, the symmetrized graph Laplacian converges to the Laplace–Beltrami operator with homogeneous Neumann boundary conditions, while the Dirichlet problem is recovered by a truncated graph Laplacian obtained by retaining only points sufficiently far from the boundary. In that setting, the eigenvalue error is
\[
\mathcal{O}\!\left(\frac{\sqrt{\log n}}{\epsilon^{d/2+1}\sqrt{n}}+\epsilon^{1/2}\right),
\]
and with optimized $\epsilon$ this becomes $\mathcal{O}((\log n/n)^{1/(2d+6)})$ [2110.06988].

Normalization is not a cosmetic choice. In the three-parameter family
\[
L u := -\frac{1}{\varrho^p}\operatorname{div}\!\Big(\varrho^q \nabla\Big[\frac{u}{\varrho^r}\Big]\Big),
\]
the nearly separated two-cluster analysis shows that a uniform spectral gap above the second eigenvalue appears only in the balanced case $q=p+r$; in unbalanced cases, only a ratio gap may remain [1909.06389]. Likewise, on unions of intersecting manifolds of different dimensions, the normalized graph Dirichlet energy converges to a dimension-adaptive limit on all manifolds simultaneously, whereas the unnormalized energy and its associated graph Laplacian asymptotically only see the variations within the manifold of the highest dimension [2509.24458].

A plausible implication is that several recurrent identifications in applied work—normalized versus unnormalized, density-corrected versus standard, Neumann versus Dirichlet, balanced versus unbalanced cluster normalization—cannot be treated as minor implementation details. The limit operator, gap structure, and even which geometric component remains visible in the spectrum depend on them [2101.09875][2110.06988][1909.06389][2509.24458].

## 5. Extensions beyond smooth closed manifolds

The scope of spectral convergence is wider than smooth, boundaryless manifolds in Euclidean space. On pcf self-similar fractals with arbitrary Borel regular probability measure of full support, a sequence of finite-dimensional weighted graph Laplacians approximates the fractal Laplacian in norm resolvent sense, and therefore heat semigroups, spectral projections, Hausdorff convergence of spectra, and eigenfunctions in energy norm follow as consequences [1704.00064].

For vector-valued problems, the principal-bundle framework extends the scalar theory to connection Laplacians on associated bundles. In this setting, eigenvalues and eigenvectors of discrete connection Laplacians converge, in the limit of infinitely many independent random samples, to the spectrum of the continuous connection Laplacian; the framework covers manifolds with and without boundary and non-uniform sampling, and subsumes Laplacian Eigenmaps, Diffusion Maps, Vector Diffusion Maps, and Orientable Diffusion Maps as special cases [1306.1587].

The 2025 Ricci-limit theory further weakens geometric regularity assumptions. Quantitative high-probability bounds are obtained for $\epsilon$-neighborhood graph Laplacians on closed Riemannian manifolds under a uniform lower Ricci curvature bound, a positive lower volume bound, and an upper diameter bound, without injectivity radius or upper sectional curvature bounds. The same framework extends to non-collapsed Ricci limit spaces, yielding spectral approximation of weighted Laplacians on manifolds with non-smooth points [2506.07427].

Periodic metric graphs supply a different extension. There the spectrum consists of absolutely continuous bands plus flat bands, and each spectral band is localized by Dirichlet and Neumann eigenvalues on a finite fundamental domain; the discrete-metric relation is explicit through the Cattaneo correspondence [1406.7523].

This collection suggests that spectral convergence is best understood as an approximation paradigm for Dirichlet forms and self-adjoint operators across manifold, fractal, metric-graph, and non-smooth metric-measure settings, rather than as a theorem tied to one ambient smooth category.

## 6. Discrete comparison principles, graph degeneration, and random-matrix limits

Not all spectral convergence problems pass through a continuum. In a general discrete-graph framework with countable vertex set $X$, symmetric edge weights $b$, vertex measure $m$, and non-negative potential $c$, the perturbation from $L_0$ to $L_c$ admits the exact comparison formula
\[
\sum_{n=1}^\infty \left(\lambda_n(c)-\lambda_n(0)\right)=\sum_{x\in X} c(x)m(x).
\]
If $c\in \ell^1(X,m)$, the eigenvalue differences form a summable null sequence, so $\lambda_n(c)-\lambda_n(0)\to 0$ as $n\to\infty$; if all eigenvalues coincide, then $c$ must vanish, yielding an Ambarzumian-type rigidity statement [2412.15937].

Graph degeneration can also force effective reduction. When intra-cluster edge weights tend to infinity, the Laplacian converges in resolvent sense to the Laplacian of a graph in which each tightly connected cluster is collapsed to a single node. In the undirected case, only the cluster partition and aggregated inter-cluster weights enter. In the directed case, the limit depends on the left and right kernel structure of the cluster Laplacian blocks, so the effective graph is sensitive to asymmetry inside the highly connected subgraphs [2601.18057].

Self-loops alter the discrete spectrum in a different way. For graphs with self-loops, the Laplacian decomposes as
\[
\mathcal{L}(G)=\mathcal{L}(G^o)+\sum_{(i,i)\in E} e_i e_i^T,
\]
and the lifting construction produces a self-loop-free graph $\hat G$ such that
\[
\sigma(\mathcal{L}(G)) \subseteq \sigma(\mathcal{L}(\hat G)) \cap [0,\,2d_{\max}(G^o)+1].
\]
If the graph is pseudo-connected, then the Laplacian is positive definite [1505.08133].

At the random-matrix scale, convergence may concern the empirical spectral distribution rather than individual low-lying eigenpairs. For centered scaled Laplacians of generalized Wigner matrices with a variance profile converging to a graphon, the empirical spectral distribution converges weakly in probability to a symmetric probability measure $\nu$, and the moments of $\nu$ admit graph-homomorphism formulas; spectral norm asymptotics are also derived [2011.07912].

Across these discrete regimes, the limit object depends on what is being varied: perturbing a potential yields an exact eigenvalue-shift identity, strengthening clusters yields a coarse-grained Laplacian, adding self-loops yields a lifted-graph inclusion principle, and increasing matrix size with graphon-convergent variance profile yields convergence of spectral distributions. The common theme is not a single limiting geometry, but the controlled transfer of spectral information under graph modification, approximation, or degeneration.

Source: https://www.emergentmind.com/topics/spectral-convergence-of-graph-laplacians