Research
Score-Repellent Monte Carlo: Toward Efficient Non-Markovian Sampler with Constant Memory in General State Spaces
arXiv:2604.22948v1 Announce Type: new Abstract: History-dependent sampling can reduce long-run Monte Carlo variance by discouraging redundant revisits, but existing schemes typically encode history th
arXiv:2604.22948v1 Announce Type: new Abstract: History-dependent sampling can reduce long-run Monte Carlo variance by discouraging redundant revisits, but existing schemes typically encode history through empirical measure on finite state spaces, which is infeasible in high-dimensional discrete configuration spaces or ill-posed in continuous domains. We propose Score-Repellent Monte Carlo (SRMC) framework that summarizes trajectory history by a running average of score evaluations in R^d, where d is the dimension of the score and state representation. This history is converted into a surrogate target through an exponential score tilt, indexed with alpha that represents the strength of repellence in controlling the magnitude of the history-based repulsion. The surrogate family is normalization-free in the standard MCMC sense, yielding a generic wrapper: at each iteration, any base kernel targeting pi can instead be run on the current surrogate pi_{heta_n} while the history is updated online. We analyze the coupled evolution of the history recursion and Monte Carlo estimators using stochastic approximation with controlled Markovian noise, establishing almost sure convergence and a joint central limit theorem. We further identify regimes in which the asymptotic covariance decreases as alpha increases, with scaling O(1/alpha), extending the near-zero-variance effect of finite-state history-dependent samplers to general state spaces with constant memory. Experiments on continuous targets and discrete energy-based models demonstrate improved estimator variance and mode coverage, while retaining O(d) memory usage and modest per-iteration overhead.
Related
- When Langevin Monte Carlo Meets Randomization: Non-asymptotic Error Bounds beyond Log-Concavity and Gradient Lipschitzness
- On two ways to use determinantal point processes for Monte Carlo integration
- A fast and Generic Energy-Shifting Transformer for Hybrid Monte Carlo Radiotherapy Calculation
- One-Step Score-Based Density Ratio Estimation
- Scalable Model-Based Clustering with Sequential Monte Carlo
Source: arXiv cs.LG | 2026-04-28