stat.MLJul 15, 2026

Spectral Concentration and Recovery in Sparse High-Dimensional Random Geometric Graphs

Authors: Manuel FernandezYizhe Zhu

Abstract

We study sparse random geometric graphs generated by connecting pairs of high-dimensional vectors whose inner product exceeds a threshold. The latent vectors are sampled either uniformly from the sphere or from a standard Gaussian distribution. Although every edge appears with probability pp, the edges are dependent through their shared latent vectors. For the spherical model, at the connectivity scale np=Ω(logn)np=Ω(\log n), we prove AEA=O(nplogn+npτ)\|A-\mathbb E A\|=O\left(\sqrt{np\log n}+npτ\right), with high probability, where ττ is the cap threshold. This sharpens the spectral norm bound of Liu, Mohanty, Schramm, and Yang (2023) under weaker assumptions. An analogous result holds for the Gaussian model after removing the fluctuations of the vector norms, yielding improved global synchronization guarantees for the homogeneous Kuramoto model. We then recover the latent geometry from the leading eigenspace. When nplognnp\gg\log n, both the latent vector and relative Gram matrix errors vanish provided dnplog(1/p)/lognd\ll np\log(1/p)/\log n. The required lower dimension is only dlog(1/p)d\gg\log(1/p) for the spherical model and dlog2(1/p)lognd\gg\log^2(1/p)\log n for the Gaussian model, improving the recovery guarantees of Li and Schramm (2023). Finally, we prove the first exact recovery result for the Gaussian mixture block model of Li and Schramm (2023). At the optimal connectivity scale np=Ω(logn)np=Ω(\log n), a polynomial-time semidefinite program exactly recovers all labels in a moderate-separation regime, whereas larger separation makes exact recovery impossible because isolated vertices appear with high probability. Our proofs combine orthogonal polynomial expansions, decoupling, and matrix concentration, avoiding the trace-moment arguments used in previous work.

Explore similar work

Jun 12, 2026math.ST

Recovery thresholds for hidden weighted sparse graphs

Recovering structural information from noisy high-dimensional data is a fundamental task in statistical inference. We investigate the recovery thresholds for a graph hidden in a randomly weighted complete graph. Specifically, an unknown graph HHnH^* \in H_n is chosen uniformly at random, and hidden in a complete graph of nn vertices as follows: the weight of an edge eHe \in H is distributed independently according to PnP_n; otherwise the weight is distributed independently according to QnQ_n. The goal is to recover almost all of HH from these edge weights. Assuming a local Lipschitzness of the Rényi divergence between distributions PnP_n and QnQ_n, and a mild density condition for the graphs HnH_n, we give a unified characterization of the information-theoretic limit for recovering almost all of HH (also known as almost exact recovery). Our characterization connects the KL divergence between PnP_n and QnQ_n to the logarithm of the first moment threshold of HH in the Erdős-Rényi random graph model G(n,p)G(n,p). Our lower bound also extends to the task of partial recovery, in which only a constant λλ-fraction of HH needs to be recovered. Last but not least, for certain Bernoulli and Exponential regimes, and for Gaussian distributions, we are able to show an All-or-Nothing (AoN) threshold phenomenon at the exponential scale.
Zhe Hou, Jingcheng Liu
Apr 22, 2026cs.LG

Spectral Embeddings Leak Graph Topology: Theory, Benchmark, and Adaptive Reconstruction

Graph Neural Networks (GNNs) excel on relational data, but standard benchmarks unrealistically assume the graph is centrally available. In practice, settings such as Federated Graph Learning, distributed systems, and privacy-sensitive applications involve graph data that are localized, fragmented, noisy, and privacy-leaking. We present a unified framework for this setting. We introduce LoGraB (Local Graph Benchmark), which decomposes standard datasets into fragmented benchmarks using three strategies and four controls: neighborhood radius dd, spectral quality kk, noise level σσ, and coverage ratio pp. LoGraB supports graph reconstruction, localized node classification, and inter-fragment link prediction, with Island Cohesion. We propose AFR (Adaptive Fidelity-driven Reconstruction), a method for noisy spectral fragments. AFR scores patch quality via a fidelity measure combining a gap-to-truncation stability ratio and structural entropy, then assembles fragments using RANSAC-Procrustes alignment, adaptive stitching, and Bundle Adjustment. Rather than forcing a single global graph, AFR recovers large faithful islands. We prove heat-kernel edge recovery under a separation condition, Davis--Kahan perturbation stability, and bounded alignment error. We establish a Spectral Leakage Proposition: under a spectral-gap assumption, polynomial-time Bayesian recovery is feasible once enough eigenvectors are shared, complementing AFR's deterministic guarantees. Experiments on nine benchmarks show that LoGraB reveals model strengths and weaknesses under fragmentation, AFR achieves the best F1 on 7/9 datasets, and under per-embedding (ε,δ)(ε,δ)-Gaussian differential privacy, AFR retains 75% of its undefended F1 at ε=2ε=2. Our anonymous code is available at https://anonymous.4open.science/r/JMLR_submission
Thinh Nguyen-Cong, Truong-Son Hy, Thang N. Dinh
Jul 20, 2026stat.ML

Mixing-Free and Signal-Optimal Learning of Gaussian Graphical Models from Glauber Dynamics

Gaussian graphical model selection is usually studied under independent sampling, but in many applications the data arise as a single trajectory of a dependent stochastic process. We study exact recovery of the graph from one trajectory of random-scan Gaussian Glauber dynamics. Existing techniques for this problem either inherit the mixing time of the chain, which can be super-polynomial in the dimension pp without strong assumptions, or are suboptimal in the minimum normalized edge strength κκ. We propose two algorithms that are mixing-free and attain the κ2κ^{-2} dependence of the information-theoretic lower bounds. Both instantiate a shared dueling-neighborhood search meta-algorithm with a local statistic built directly from the update sequence. For every fixed precision matrix and deterministic initialization, the first algorithm fits a least-squares regression at the updates of each node and has pointwise recovery horizon O~(pd2/κ2)\widetilde O(pd^{2}/κ^{2}), where dd is the maximum degree. Its horizon depends logarithmically on a local conditioning quantity and on the initialization potential. The second algorithm is based on counting occurences of a specific update pattern and requires O~(pd4/κ2)\widetilde O(pd^{4}/κ^{2}) updates, with no dependence on any condition number. The central technical challenge is that both statistics are built from dependent, non-stationary observations. Our analysis tackles this by demonstrating how to extract fresh Gaussian innovations from the update sequence, which yields mixing-free control of appropriate quantities. Neither the algorithms nor their analyses invoke stationarity, a spectral gap, or mixing conditions.
Vignesh Tirukkonda, Gautam Dasarathy