{"id":"arxiv/2609.38710","url":"https://arxiv.org/abs/2609.38710","tags":["math.ST","cs.DS","math.NA"],"score":0,"title":"High-accuracy simulation of Picard HMC, part I: Gaussian cloud correction","domain":"math.ST","source":"arxiv","authors":["Fan Chen","Sinho Chewi","Jianfeng Lu","Matthew S. Zhang"],"reviews":0,"subject":"math.ST","summary":"We study the problem of sampling from a continuous density $π\\propto \\exp(-V)$ on $\\mathbb R^d$, where $V\\in C^2(\\mathbb R^d)$ has a $β$-Lipschitz gradient and $π$ satisfies a logarithmic Sobolev inequality with constant $α^{-1}$, and write $κ= β/α$. We introduce the Gaussian cloud sampler, which achieves total variation accuracy $\\varepsilon$ using $\\widetilde O(κd^{1/5}\\,\\text{polylog}(1/\\varepsilon))$ gradient queries in expectation. The algorithm uses first-order rejection sampling (FORS) to correct the law of smoothed Picard HMC trajectories. To do so, we represent the iterates of the…","category":"research","comments":0,"disputed":false,"division":"math","mentions":0,"discussed":0,"published":"2026-09-30","reactions":{}}