Skip to content
Preprint

An FPRAS for Antiferromagnetic Ising Models on Random Regular Bipartite Graphs

Aug 2026 · 0 citations · 22 references
Computer Science Mathematics

Abstract

We design randomized approximation schemes for the partition function of antiferromagnetic Ising models with uniform external field on random regular bipartite graphs. Our algorithm generalizes the approach of Kocurek, Oveis Gharan and Tjowasi (arXiv, 2026) for hard-core models on the same random graph model beyond the uniqueness threshold. We show that, as long as $\lambda$ is upper bounded by a constant and $\lambda(1 - \beta) \lesssim \Delta^{-1/2}$, an efficient randomized algorithm approximates the partition function with high probability. The algorithm first truncates configurations that are large on either side of the bipartition and then samples from Gibbs distributions conditioned on fixed sizes on one or both sides. To choose an optimal truncation bound, we establish concentration properties of the Gibbs distribution on random regular bipartite graphs. Then we apply high-dimensional expansion and prove trickle-down theorems to obtain fast samplers for the conditioned distributions.

View source

Similar papers

Preprint Aug 2026

Scaling Limits for Ising Models on Inhomogeneous Random Graphs and Applications

In this paper, we derive quenched scaling limits for linear functionals and the empirical spin field of Ising models on inhomogeneous random graphs generated by a graphon (encompassing both dense and sparse graphs), in the high-temperature regime. We first prove a joint central limit theorem (CLT) for finite collections of linear statistics of the spin configurations, where the limiting covariance is characterized by the resolvent of the associated graphon integral operator. Building on this result, we establish functional CLTs for the average magnetization and for the spin field indexed by suitable classes of regular test functions. We further prove convergence of the full empirical spin field, viewed as a random generalized function in negative Sobolev spaces. These scaling limits provide applications to both Bayesian neural networks and causal inference. Specifically, for the former, we derive infinite-width Gaussian-process limits for two-layer Bayesian neural networks with Ising-dependent output-layer signs, while for the latter, we establish the asymptotic normality of H\'{a}jek estimators for average treatment effects under network interference.

Sanchayan Bhowal, Anirban Chatterjee, Somabha Mukherjee · 0 citations
Preprint Jul 2026

Marked vertex search on disordered graphs with Rosenzweig-Porter phases

These results establish a direct and quantitative link between random matrix disorder on graphs and the performance of continuous-time quantum walk search, and suggest that disorder, rather than being merely an obstacle, can be exploited as a tunable parameter in quantum search protocols.

Sabyasachi Chakraborty, T. Čadež, Sonjoy Majumder et al. · 0 citations
Preprint Aug 2026

The critical probability for percolation on finite graphs

We determine the critical probability for Bernoulli bond percolation on essentially any finite graph. Namely, letting $\lambda(G)$ denote the spectral radius (maximum eigenvalue) of $G$, we prove that the critical probability is at $1/\lambda(G)$: above this probability there is typically a component of order $\Omega(\lambda(G))$, whereas below it all components are of order at most $O(\sqrt{|G|})$. These results in particular confirm a conjecture of Krivelevich and Samotij about percolation on graphs of a given average degree, and vastly extend theorems of Bollob\'as, Borgs, Chayes, and Riordan, who proved analogous results but only for dense graphs. Our theorems are optimal in many regimes, and also demonstrate that percolation has an unexpectedly subtle behaviour on graphs whose spectral radius is roughly the square root of their maximum degree.

Micha Christoph, Patryk Morawski, Yuval Wigderson · 0 citations
Conference Jul 2026

Mixing Time Bounds for Asymmetric Random Walks on Hypercubes, Cycles, and Tori Using Coupling

Mixing time bounds for Markov chains play a central role in characterizing the sample complexity of learning and inference from correlated data. While the mixing behavior of symmetric random walks on standard graph structures such as cycles, tori, and hypercubes is well understood, the impact of transition asymmetry remains less explored. In this work, we study the mixing times of lazy, asymmetric random walks on cycles, tori, and hypercubes, motivated by their relevance in practical applications. For the $n$-cycle, we develop a novel coupling construction that yields an order-wise tight upper bound $O\left(\frac{n^{2}}{p+q}\right)$, explicitly capturing the dependence on asymmetric transition probabilities $p$ and $q$. Numerical results indicate that this dependence is highly accurate. Building on this result, we derive corresponding bounds for $d$-dimensional tori. For asymmetric random walks on $n$-dimensional hypercubes, motivated by applications, we consider the problem of estimating expectations of functions that depend only on a subset $\Delta \ll n$ of coordinates. We show that the effective sample complexity improves to $O(n \log \Delta)$, compared to $O(n \log n)$ for the full chain. In all cases, our bounds recover the tightest known results for symmetric walks as special cases.

Mrudula A Mahindrakar, Hrushikesh A Kant, Avhishek Chatterjee · 0 citations
Preprint Aug 2026

The Hard-Core Model on Bipartite Spectral Expanders: Counting and Sampling at All Fugacities

We study approximate counting and sampling algorithms for the hard-core model on $\Delta$-regular bipartite graphs under a spectral expansion condition. Let $M_G$ be the biadjacency matrix of $G$. For every fixed $\xi\in(0,1)$, we give an FPRAS for the hard-core partition function and an efficient approximate sampler whenever \[ \lambda\leq \frac{1-\xi}{\sigma_2(M_G)}. \] The main idea is to introduce a family of quadratic tilts in the left-right occupation imbalance and show that each tilted measure can be sampled efficiently using Glauber dynamics. A discrete Gaussian identity expresses the original hard-core model as an exact positive mixture of these tilted measures; truncation and simulated annealing then yield efficient counting and sampling algorithms. For the complementary high-fugacity regime, we refine the polymer-model approach and show that the required phase-dominance and cluster expansion conditions follow from the singular-spectrum bound alone. Combining the two regimes, we obtain efficient approximate counting and sampling at every fugacity $\lambda>0$ whenever \[ \sigma_2(M_G)\leq c\left(\frac{\Delta^2}{\log(\mathrm e\Delta)}\right)^{1/3} \] for an absolute constant $c>0$. In particular, this recovers all-fugacity algorithms for random $\Delta$-regular bipartite graphs for all sufficiently large $\Delta$, while providing an efficiently verifiable certificate of their success on a given instance.

Ijay Narang, Will Perkins · 0 citations
Preprint Jul 2026

Limiting spectral distribution for the adjacency matrix of the Watts-Strogatz random graph

The Watts-Strogatz random graph model on $n$ vertices with parameters $K$ (a positive even integer) and $p \in [0, 1]$ is constructed in two steps. First, one starts with a ring lattice on $n$ vertices, where each vertex is connected to its $K/2$ nearest neighbors on each side. Each edge in turn is then independently rewired with probability $p$ by replacing one endpoint with a uniformly chosen vertex not already adjacent to it. We study the empirical eigenvalue distribution of the adjacency matrix for this model, whose entries are highly dependent due to the rewiring construction. In the regime where both $K$ and $pK$ grow to infinity with the vertex size $n$, we show that, after appropriate scaling, the empirical eigenvalue distribution converges to the semicircle law. The proof is based on a novel coupling argument that approximates the adjacency matrix by a sum of two independent random matrices, one a sparse Wigner matrix and the other a random band matrix. In the case where $K$ and $p$ remain fixed, we propose conjectural formulas for the first five moments of the limiting eigenvalue distribution. These conjectures are supported by a convergence result relating the Watts-Strogatz model to another random graph model, together with numerical simulations.

Grégoire Meunier, Sean O’Rourke · 0 citations

We use cookies to run the site and, with your consent, for analytics and to show ads. See our Cookie Policy.