Skip to content

High-Magnetization Sampling at Low Temperatures: Ising Models and Bayesian Sparse Linear Regression

Sep 2026 · 0 citations
Computer Science Mathematics

Abstract

Sparsity is a powerful structural resource in optimization and statistics. We develop frameworks for leveraging sparsity in sampling problems over the Hamming slice $\mathcal{X}_k^d:=\{\mathbf{x}\in\{\pm 1\}^d:|\{i:\mathbf{x}_i=1\}|=k\}$, in high-dimensional regimes where $k\ll d$ (i.e., where $\mathcal{X}_k^d$ is \emph{highly magnetized}). We use our frameworks to design improved samplers for canonical problems in the study of \emph{Ising models} and \emph{Bayesian sparse linear regression}. Our first main result considers the \emph{Sherrington--Kirkpatrick} (SK) model restricted to fixed-magnetization slices $\mathcal{X}_k^d$. We give a polynomial-time sampler for fixed-magnetization SK models at any inverse temperature $\beta>0$, under arbitrary external fields, provided that $k\le c_\beta d$ for an appropriate constant $c_\beta$. By combining this result with an annealing strategy for estimating normalizing constants, we obtain polynomial-time samplers for the SK model at arbitrarily low temperatures under a sufficiently strong external field of strength $h$. In the large-$\beta$ limit, our framework permits sampling at field strengths within constant factors of the \emph{Almeida--Thouless line} delineating the replica-symmetric and replica-symmetry-breaking regions ([dAT78]), improving polynomially over the field strength $h(\beta)$ required by the recent work of [BAR26]. Our second main result concerns the measurement complexity of polynomial-time Bayesian sparse linear regression. Recent work by [KSTZ25] shows how to sample from the canonical \emph{Gaussian spike-and-slab posterior} with expected sparsity $k$, at any signal-to-noise ratio, given $n\gtrsim k^3\log^3 d$ Gaussian measurements. We improve this requirement to $n\gtrsim k^{3/2}\log^2 d+k\log^3 d$, using a common sparsity-aware framework underlying both our results.

View source

Similar papers

Preprint Sep 2026

The anisotropic local law for sample covariance matrices under quadratic-form concentration

We study sample covariance matrices $K = \frac{1}{N} \sum_{i=1}^N \mathbf{x}_i \mathbf{x}_i^* \in \mathbb{R}^{n \times n}$ in the proportional regime $n \asymp N$. The columns $ \mathbf{x}_1, \ldots, \mathbf{x}_N \in \mathbb{R}^n$ are independent and centered, with common covariance $\mathbb{E} \mathbf{x}_i \mathbf{x}_...

Ren-Yuan Ma, Theodor Misiakiewicz · 0 citations
Preprint Oct 2026

Efficient Posterior Sampling for $\mathbb Z_2$ Synchronization

Consider the $\mathbb Z_2$ synchronization problem \[ \boldsymbol{Y} = \frac{\lambda}{\sqrt n} \theta\theta^{\top} + \boldsymbol{Z}, \] where $\theta$ is uniform on $\{-1,1\}^n$ and $\boldsymbol{Z}$ is an independent Gaussian Wigner matrix with off-diagonal variance one. We give a polynomial-time posterior sampling alg...

Zhang-Song Li · 0 citations
Preprint Aug 2026

Minimax Optimal Estimator and Improved Error Rate for the MLE in Logistic Regression with Gaussian Design

We study finite-sample parameter estimation in logistic regression with Gaussian design, where the goal is to estimate $\mathbf{\theta}^*\in \mathbb{R}^d$ with $R=\|\mathbf{\theta}^*\|_2\ge 1$ from i.i.d. samples $\{(\mathbf{x}_i,y_i)\}_{i=1}^n,$ $\mathbf{x}_i \sim N(0,\mathbf{I}_d)$, $y_i\mid \mathbf{x}_i \sim \mathrm...

Jun-Ren Chen, Arya Mazumdar · 0 citations
Preprint Sep 2026

A spectral gap for Metropolis-adjusted Langevin algorithm with a uniformly randomized step size

Let $\pi(\mathrm{d} x)\propto e^{-U(x)}\, \mathrm{d} x$ on $\mathbb{R}^d$, where $U$ is continuously differentiable and $m$-strongly convex with a globally $L$-Lipschitz gradient, $0<m\leq L<\infty$, and $\kappa=L/m$. Fixed-step Metropolis-adjusted Langevin algorithm (MALA) has known warm-start mixing-time upper bounds...

Qian Qin · 0 citations
Preprint Aug 2026

Uniform logarithmic Sobolev inequalities for the 2D Coulomb gas at the diffusive temperature scale

For $N\ge2$ and $\beta>0$, consider the canonical 2D Coulomb gas ensemble \begin{equation} \mathrm{d}\mathbb{P}_{N,\beta}(Z) =\mathsf{Z}_{N,\beta}^{-1}e^{-\beta|Z|^2/2} \prod_{i<j}|z_i-z_j|^{\beta/N}\,\mathrm{d}Z \qquad\text{on }\mathbb{C}^N. \end{equation} The factor $N^{-1}$ in the pair exponent places the ensemble a...

Matthew Rosenzweig · 1 citation
Preprint Aug 2026

Spectral phase transitions in Gaussian multi-index models

It is proved that the AMP-derived preprocessing of [Defilippis et al., 2025] is optimal among all bounded matrix-valued preprocessing maps of any fixed dimension and the general spectral conjecture of [Defilippis et al., 2025] is proved.

Florent Krzakala, Pierre Mergny, Vanessa Piccolo · 0 citations

Related blog posts

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