---
title: Sparse Random Geometric Graphs
url: https://www.emergentmind.com/topics/sparse-random-geometric-graphs-rggs
type: topic
---

# Sparse Random Geometric Graphs

Searching arXiv for recent and foundational work on sparse random geometric graphs to ground the article.
Sparse random geometric graphs (RGGs) are random graphs whose adjacency structure is induced by an underlying geometric embedding, with sparsity enforced by a small connection radius, a bounded mean degree, or an edge density that vanishes relative to graph size, depending on the model. Across the literature, the term covers several non-equivalent regimes: shrinking-radius binomial graphs on bounded domains, fixed-radius supercritical Gilbert graphs on Poisson point processes, and high-dimensional spherical models with edge probability \(p=\alpha/n\). In all of these settings, sparsity makes local geometry, fragmentation, and metric effects decisive, but it does so in model-specific ways rather than through a single universal definition [2007.02453] [2411.02046] [2111.11316].

## 1. Model families and notions of sparsity

The standard two-dimensional finite-window model studied in the spectral-topological literature is \(G(n,\ell)\): \(n\) vertices are sampled independently and uniformly on the unit square, and two vertices are joined whenever their Euclidean distance is at most \(\ell\in[0,\sqrt{2}]\). In this model, \(\ell=0\) yields only isolated vertices, \(\ell=\sqrt{2}\) yields the complete graph, and the average degree is explicitly
\[
k=(n-1)\left(\pi \ell^2-\frac{8}{3}\ell^3+\frac{1}{2}\ell^4\right),
\]
so that for \(\ell\ll 1\),
\[
k\approx (n-1)\pi \ell^2\sim n\pi \ell^2.
\]
This makes \(n\ell^2\) the natural sparse-regime control parameter in two dimensions [2007.02453].

A closely related adjacency-spectrum model uses the unit square with periodic boundary conditions, i.e. the unit torus, and varies a radius \(r\). There sparsity is described operationally by small \(r\), low mean degree \(\langle k\rangle\propto r^2\), low probability of full connectivity, and fragmentation into many disconnected components. In the numerical study with \(N=10^3\), the spacing-statistics regime \(r\in[0.03,0.4]\) spans a low-density fragmented phase through a connected phase toward the complete graph [1608.01154].

Another important notion of sparsity is the thermodynamic regime. For RGGs on the \(d\)-dimensional torus \(\mathbb T^d\), with adjacency rule \(A_{ij}=\mathbf 1\{\|x_i-x_j\|_p\le r_n,\ i\ne j\}\), the average degree is
\[
a_n=\theta^{(d)}nr_n^d.
\]
The thermodynamic regime fixes \(a_n\equiv \gamma<\infty\), equivalently \(r_n\asymp n^{-1/d}\). This is sparse in the graph-theoretic sense because \(|E|=\Theta(n)\) and degree does not diverge with \(n\) [1910.08869].

Sparse RGGs also arise in continuum percolation. In the Gilbert graph on a homogeneous Poisson point process \(\mathcal P_\lambda\subset\mathbb R^d\), edges join points within a fixed radius \(r\), and the expected degree is \(\lambda \upsilon_d r^d\). Here sparsity means fixed local connectivity scale rather than shrinking-radius finite-sample asymptotics. The regime of interest is supercritical, \(r>r_c(\lambda)\), so there exists a unique infinite connected component \(\mathcal H\) on which large-scale metric questions can be posed [2408.11615] [2411.02046].

High-dimensional sparse RGGs use yet another parameterization. In the spherical model \(\mathsf{Geo}_d(n,p)\), latent vectors are sampled uniformly on \(\mathbb S^{d-1}\), and edges are created by an inner-product threshold chosen to give marginal edge probability \(p\). The sparse case of main interest is
\[
p=\frac{\alpha}{n},
\]
so average degree remains \(O(1)\) while ambient dimension \(d\) varies [2111.11316]. A related spectral study uses the same spherical threshold model \(G(n,d,p)\) and investigates both \(p=\alpha/n\) and the more general sparse regime \(p\to 0\) with \(np\to\infty\) [2507.06556].

| Model | Sparse regime in the cited work | Structural hallmark |
|---|---|---|
| \(G(n,\ell)\) on the unit square | \(\ell\ll 1\), \(k\approx n\pi\ell^2\) | many isolated vertices and tiny components |
| Unit-torus adjacency RGG | small \(r\), \(\langle k\rangle\propto r^2\) | fragmented spectrum and motif effects |
| Thermodynamic torus RGG | \(a_n=\theta^{(d)}nr_n^d\equiv \gamma\) | constant expected degree |
| One-dimensional circular RGG | \(nr_n=o(1)\) with \(n^3r_n^2=\omega(1)\) | sparse but wedge/triangle counts still diverge |
| Spherical high-dimensional RGG | \(p=\alpha/n\) | sparse geometry versus Erdős–Rényi comparison |

## 2. Fragmentation, motifs, and clustering

In low-dimensional finite-window RGGs, the first structural effect of sparsity is fragmentation. In the unit-square model, small \(\ell\) produces low expected degree, many isolated vertices, many tiny disconnected components, and a subcritical or pre-connectivity regime. The number of non-isolated vertices,
\[
V_\times(G),
\]
is therefore a particularly informative observable: \(V_\times(G)=0\) at \(\ell=0\), \(V_\times(G)=n\) at \(\ell=\sqrt{2}\), and \(V_\times(G)\ll n\) in sparse graphs. The fraction \(V_\times(G)/n\) acts as an effective order parameter for departure from the fully isolated regime, and in the cited numerical study it almost completely controls the normalized Randić and harmonic indices as well as the averaged Shannon entropy of eigenvectors [2007.02453].

The torus adjacency-spectrum study makes the same fragmentation picture explicit from the spectral side: for small \(r\), the graph is highly likely to contain many disconnected components, the spectrum is close to a superposition of independent component spectra, and localization is naturally interpreted as concentration of eigenstates on disconnected components [1608.01154]. This is the sparse-RGG analogue of modular isolation.

A common extrapolation from non-geometric sparse random graph models does not hold uniformly for RGGs: sparsity need not force clustering to vanish. In the one-dimensional circular binomial RGG with adjacency
\[
A_{ij}=I[d(X_i,X_j)\le r_n],
\qquad
d(X_i,X_j)=\min\{|X_i-X_j|,1-|X_i-X_j|\},
\]
the sparse regime is defined by
\[
nr_n=o(1),
\qquad
n^3r_n^2=\omega(1).
\]
Here expected degree is of order \(nr_n\to 0\), but wedge and triangle counts are both of order \(n^3r_n^2\), and the global clustering coefficient
\[
\mathcal C_n=\frac{\sum_{i\neq j\neq k}A_{ij}A_{jk}A_{ki}}{\sum_{i\neq j\neq k}A_{ij}A_{jk}}
\]
still satisfies
\[
\mathcal C_n\to \frac34
\]
in probability. Moreover,
\[
\frac{8n\sqrt{n}r_n\sqrt{\mathbb E[f^2(X_1)]}}{3}\big(\mathcal C_n-\mu_n\big)\Rightarrow N(0,1)
\]
under the same sparse assumptions [2602.17896]. In this specific one-dimensional model, sparse geometry preserves transitivity rather than destroying it.

## 3. Spectral and eigenvector structure in sparse low-dimensional RGGs

Sparse RGGs display a characteristic spectral crossover from localized, near-Poisson behavior to delocalized, random-matrix behavior as connectivity increases. In the unit-square model \(G(n,\ell)\), the spectral analysis is performed not on the binary adjacency matrix but on a randomly weighted symmetric matrix \(\mathbf A\) whose off-diagonal entries are Gaussian on existing edges and whose diagonal entries are also random. This produces a Poisson ensemble (PE) limit at \(\ell=0\) and a GOE limit at \(\ell=\sqrt{2}\). The spacing-ratio reference values are
\[
\langle r\rangle_{PE}\approx 0.3863,
\qquad
\langle r\rangle_{GOE}\approx 0.5359,
\]
while eigenvector localization is measured through the participation ratio (“IPR” in the paper) and Shannon entropy \(S(G)\) [2007.02453].

The same study finds that the topological sparse-to-dense crossover occurs at
\[
\ell^*=\mathcal C n^{-\gamma},
\qquad
\gamma\approx 0.5075\approx \frac12
\]
for \(V_\times(G)/n\), \(R(G)/(n/2)\), and \(H(G)/(n/2)\). Thus the principal topological crossover scale is \(\ell^*\asymp n^{-1/2}\), consistent with \(k\approx n\pi\ell^2\). For spectral quantities, the fitted exponents differ:
\[
\gamma=0.44193 \text{ for spacing ratio},\quad
\gamma=0.36272 \text{ for “IPR”},\quad
\gamma=0.44138 \text{ for Shannon entropy}.
\]
Most notably, the normalized average Shannon entropy is almost perfectly correlated with \(V_\times(G)/n\), so average eigenvector delocalization can be inferred numerically from a simple count of non-isolated vertices [2007.02453].

For the plain \(0/1\) adjacency spectrum on the two-dimensional torus, the low-radius regime is again Poisson-like, but now geometry generates an additional feature absent from many non-spatial models: sharp discrete peaks at integer eigenvalues caused by symmetric motifs, especially a large multiplicity at eigenvalue \(-1\). These peaks create a zero-spacing spike that must be removed before fitting the continuous part of the nearest-neighbor spacing distribution. Quantitatively, the Brody parameter rises from \(0.052\pm 0.005\) at \(r=0.03\) to \(0.958\pm 0.002\) at \(r=0.40\), with the main crossover around \(r=0.04\)–\(0.06\) for \(N=10^3\) [1608.01154].

Sparse Laplacian behavior is different but equally geometric. In the thermodynamic regime \(a_n\equiv \gamma\), the regularized normalized Laplacian
\[
\mathcal L = I-(D+\alpha I)^{-1/2}\left(A+\frac{\alpha}{n}\mathbf 1\mathbf 1^\top\right)(D+\alpha I)^{-1/2}
\]
is used to avoid singularities from zero degrees. For the \(\ell_\infty\) model on the torus, the low-end spectrum obeys
\[
F^{\mathcal L'}(x)\approx Cx^{d/2},
\qquad
\rho(x)\asymp x^{d/2-1},
\]
and the paper concludes that the spectral dimension satisfies
\[
d_s\approx d
\]
in the thermodynamic regime [1910.08869]. This indicates that sparse geometric randomness preserves the ambient diffusive dimension at low frequencies.

## 4. High-dimensional sparse RGGs and the disappearance of geometry

In high dimension, sparse RGGs can lose their geometric signature altogether. For the spherical model \(\mathsf{Geo}_d(n,p)\) in the constant-sparsity regime \(p=\alpha/n\), the total variation distance to Erdős–Rényi satisfies
\[
d_{\mathrm{TV}}\!\left(\mathsf{Geo}_d\!\left(n,\frac{\alpha}{n}\right),\mathsf G\!\left(n,\frac{\alpha}{n}\right)\right)\to 0
\]
whenever
\[
d=\Omega(\log^{36} n).
\]
The same paper recalls a signed-triangle lower bound showing detectability when \(d\ll \log^3 n\), so the sparse indistinguishability threshold is pinned down up to polylogarithmic factors [2111.11316].

Global spectral statistics can wash out geometry even earlier. In the spherical threshold model \(G(n,d,p)\), if
\[
p=\frac{\alpha}{n},
\qquad
d=\omega(\log^2 n),
\]
then the empirical spectral distribution of
\[
\frac{A}{\sqrt{\alpha}}
\]
converges to the same limiting law \(\nu_\alpha\) as \(G(n,\alpha/n)\). More generally, if \(np\to\infty\) and
\[
d=\omega\!\big(np\log^2(1/p)\big),
\]
then the empirical spectral distribution of
\[
\frac{A}{\sqrt{np(1-p)}}
\]
converges to the semicircle law [2507.06556]. This suggests that global spectral observables can become Erdős–Rényi-like before the full graph law is completely indistinguishable.

An entropy-based high-dimensional analysis sharpens the dependence on geometry and connection rule. For hard RGGs on the torus, and for soft RGGs in both cube and torus geometries, the graph law converges to \(G(n,p)\) as \(d\to\infty\). By contrast, for hard RGGs in the cube with coordinate kurtosis greater than \(1\), residual edge dependence survives, and the limiting maximum entropy is strictly below the Erdős–Rényi maximum. The same work derives an Edgeworth correction showing
\[
H(G)=H_d(G)+O(d^{-1/2})
\]
for the Shannon entropy of the graph distribution [2503.11418]. High-dimensional sparsity therefore does not by itself imply independent-edge behavior; the boundary geometry and the hard-versus-soft connection rule remain decisive.

## 5. Metric geometry, random walks, and algorithms

Sparse RGGs support a rich large-scale metric theory. In first-passage percolation (FPP) on the supercritical Gilbert graph built from a Poisson point process in \(\mathbb R^d\), with \(r>r_c(\lambda)\) and i.i.d. nonnegative edge weights, the growth set \(H_t\) satisfies a Euclidean-ball shape theorem:
\[
(1-\varepsilon)B(o,\varphi)\subseteq \frac1n H_n\subseteq (1+\varepsilon)B(o,\varphi)
\]
for all sufficiently large \(n\), almost surely, under the stated moment assumptions. Under exponential-moment assumptions, the quantitative shape theorem improves this to
\[
\left(1-c'\frac{\log t}{\sqrt t}\right)B(o,\varphi)\subseteq \frac1t H_t\subseteq \left(1+c'\frac{\log t}{\sqrt t}\right)B(o,\varphi),
\]
and passage times obey moderate deviations
\[
\mathbb P\left(\frac{|T(x)-\mathbb E[T(x)]|}{\sqrt{\|x\|}}>\ell\right)\le Ce^{-c\ell}
\]
for \(\ell\in[C'\log\|x\|,\sqrt{\|x\|}]\). Geodesics remain within Hausdorff distance \(\|y-x\|^{3/4+\epsilon}\) of the Euclidean segment with stretched-exponential tails [2408.11615] [2411.02046].

Random-walk metrics retain geometric information in sparse supercritical two-dimensional torus RGGs. For effective resistance \(E_{uv}\), the proposed and numerically validated model is
\[
E_{uv} \approx \beta_0 + \beta_1 \ln(d_{uv}) + \beta_2 d_{uv}^2 + \beta_3 d_{uv}^4\cos(4\theta_{uv}) + \beta_4 \left(\frac{1}{D_u} + \frac{1}{D_v}\right),
\]
with commute time given exactly by
\[
C_{uv}=\mathrm{Vol}(G)\,E_{uv}.
\]
The logarithmic term reflects the two-dimensional Laplacian Green’s function, the quadratic term the compact torus topology, and the quartic angular term the square anisotropy of the torus representation. This contrasts with the very dense regime, where commute times degenerate to an inverse-degree sum [2606.14895].

Algorithmically, sparse geometry can make exact distance problems easier than in general unit disk graphs. For two-dimensional RGGs with expected average degree \(\Theta(n^\delta)\), \(\delta\in(0,1)\), the diameter can be computed asymptotically almost surely in
\[
\tilde O\!\left(n^{\frac32(1+\delta)} + n^{2-\frac53\delta}\right)
\]
time on the square and
\[
\tilde O\!\left(n^{\frac32(1+\delta)} + n^{2-\frac13\delta}\right)
\]
time on the torus. The square bound is minimized at \(\delta=3/19\), yielding \(\tilde O(n^{33/19})\approx \tilde O(n^{1.737})\). The same paper shows that iFUB achieves a speedup of \(\tilde\Omega(n^{\delta/3})\) over naive \(O(nm)\) on square RGGs but still requires \(\Omega(nm)\) on torus RGGs, reflecting the algorithmic importance of boundary geometry [2603.16684].

## 6. Inference, transferability, and application-facing theory

Sparse RGGs have also become a setting for statistical inference from topology alone. For the toroidal \(L_\infty\)-RGG \(\mathcal G_n(m,r_n)\), a test of the latent dimension \(m\) is built from the contrast
\[
D_n=\sum_{i\neq j\neq k}A_{ij}A_{jk}A_{ki} -\left(\frac34\right)^{m_0}\sum_{i\neq j\neq k}A_{ij}A_{ik}.
\]
Under
\[
H_0:m=m_0,
\qquad
r_n=o(1),
\qquad
nr_n^m=\omega(1),
\]
the statistic
\[
\frac{\sqrt{2}D_n}{n^2\hat{\sigma}_{n2}}
\Rightarrow N(0,1),
\]
whereas under \(H_1:m\neq m_0\) it diverges at rate \(\Theta(n\sqrt{r_n^m})\) [2510.11844]. This is a sparse-but-not-ultra-sparse regime: edge density vanishes, but expected degree still diverges.

Learning-theoretic work uses sparse RGGs as wireless network models precisely because local interference yields bounded average degree rather than dense graphons. In the expanding-domain, fixed-density regime with connection radius \(r_c\), the expected number of neighbors is
\[
\pi r_c^2\rho.
\]
A transferability theory for GNNs is then obtained by comparing the RGG to a deterministic grid graph (DGG), assuming
\[
\mathbb E[\|S_n-S_{D_n}\|_2^2]=O(n^{-\alpha}),
\]
and using integral-Lipschitz graph filters. Under these assumptions, performance differences between small and large sparse geometric graphs can be bounded explicitly [2510.00896]. A closely related theory for wireless conflict graphs induced by sparse RGGs assumes an operator-norm bound
\[
\|\mathbf A-\tilde{\mathbf A}\|\le \varepsilon
\]
between RGG and DGG conflict-graph adjacencies and derives
\[
\|\Phi(\mathbf x,\tilde{\mathbf A};\mathbf H)-\Phi(\mathbf x,\mathbf A;\mathbf H)\| \le 2L\sqrt{\varepsilon C}\,\|\mathbf x\|
\]
for an \(L\)-layer GNN [2606.03794].

Sparse-RGG theory also transfers directly into sampling-based motion planning. Karaman and Frazzoli conjectured that a planner inherits a graph property if its underlying RGG has that property; the localization–tessellation framework proves this for several sparse roadmap models. In particular, PRM with
\[
r_n=\gamma\left(\frac{\log n}{n}\right)^{1/d},
\qquad
\gamma>\gamma^*=2(2d\theta_d)^{-1/d},
\]
is probabilistically complete, and Bluetooth-PRM, Soft-PRM, and Embedded-PRM inherit analogous completeness guarantees under their respective sparse-graph conditions. Standard PRM is furthermore asymptotically \(\zeta\)-optimal for some constant \(\zeta\), while Embedded-PRM is asymptotically optimal under
\[
p_n=\omega\left(\frac{\log^d n}{n}\right)
\]
[1602.05460].

Sparse RGGs are therefore not a single asymptotic object but a family of geometry-driven sparse graph models with distinct local, spectral, metric, and inferential behaviors. In low dimensions, sparsity foregrounds fragmentation, motif structure, and geometric random-walk effects; in high dimensions, it can render geometry statistically or spectrally invisible; and in applications, it supports both rigorous inference procedures and algorithmic transfer principles that depend essentially on locality rather than on dense-limit averaging.

Source: https://www.emergentmind.com/topics/sparse-random-geometric-graphs-rggs