Learning the Inverse Temperature of Ising Models under Hard Constraints using One Sample
Rohan Chauhan, Ioannis Panageas
TL;DR
The paper develops a near-linear-time method to estimate the inverse temperature $\beta$ of a truncated Ising model from a single sample, where the truncation set $S$ is the solution space of a bounded-degree $k$-SAT formula. It relies on a maximum pseudolikelihood estimator, and its consistency is proven by bounding the first and second derivatives of the pseudo-likelihood: a controlled gradient at $\beta^*$ and a strongly convex objective with parameter $\Omega(n/\Delta^3)$, shown via an independent-set conditioning argument and a Lovász Local Lemma construction that guarantees many one-flip neighbors within $S$. The analysis combines exchangeable-pair techniques for concentration with intricate combinatorial arguments (LLL, disjoint neighborhoods, and magnetization bounds) to bound the Hessian from below, enabling a projected gradient-descent algorithm running in $\mathcal{O}(\Delta^3 n \log n)$ time that achieves $|\hat{\beta}-\beta^*| \le c\Delta^3/\sqrt{n}$ with high probability, for $k \ge \Omega(\Delta^3(1+\log(d^2k+1)))$. These results extend single-sample learning to hard-constrained Ising models and provide a quantitative bridge between truncation structure and statistical-optimization guarantees, with implications for parameter inference in complex, non-product graphical models shaped by logical constraints.
Abstract
We consider the problem of estimating inverse temperature parameter $β$ of an $n$-dimensional truncated Ising model using a single sample. Given a graph $G = (V,E)$ with $n$ vertices, a truncated Ising model is a probability distribution over the $n$-dimensional hypercube $\{-1,1\}^n$ where each configuration $\mathbfσ$ is constrained to lie in a truncation set $S \subseteq \{-1,1\}^n$ and has probability $\Pr(\mathbfσ) \propto \exp(β\mathbfσ^\top A\mathbfσ)$ with $A$ being the adjacency matrix of $G$. We adopt the recent setting of [Galanis et al. SODA'24], where the truncation set $S$ can be expressed as the set of satisfying assignments of a $k$-SAT formula. Given a single sample $\mathbfσ$ from a truncated Ising model, with inverse parameter $β^*$, underlying graph $G$ of bounded degree $Δ$ and $S$ being expressed as the set of satisfying assignments of a $k$-SAT formula, we design in nearly $O(n)$ time an estimator $\hatβ$ that is $O(Δ^3/\sqrt{n})$-consistent with the true parameter $β^*$ for $k \gtrsim \log(d^2k)Δ^3.$ Our estimator is based on the maximization of the pseudolikelihood, a notion that has received extensive analysis for various probabilistic models without [Chatterjee, Annals of Statistics '07] or with truncation [Galanis et al. SODA '24]. Our approach generalizes recent techniques from [Daskalakis et al. STOC '19, Galanis et al. SODA '24], to confront the more challenging setting of the truncated Ising model.
