Related papers: Accelerating sampling via asymptotic relaxation en…
We propose quantum algorithms that provide provable speedups for Markov Chain Monte Carlo (MCMC) methods commonly used for sampling from probability distributions of the form $\pi \propto e^{-f}$, where $f$ is a potential function. Our…
This paper considers a sampling-based approach to stability verification for piecewise continuous nonlinear systems via Lyapunov functions. Depending on the system dynamics, the candidate Lyapunov function and the set of initial states of…
We consider the extreme value statistics of correlated random variables that arise from a Langevin equation. Recently, it was shown that the extreme values of the Ornstein-Uhlenbeck process follow a different distribution than those…
We propose a novel method for sampling from unnormalized Boltzmann densities based on a probability flow ordinary differential equation (ODE) derived from linear stochastic interpolants. The key innovation of our approach is the use of a…
We propose and study a new multilevel method for the numerical approximation of a Gibbs distribution $\pi$ on $\mathbb{R}^d$, based on (overdamped) Langevin diffusions. This method inspired by \cite{mainPPlangevin} and…
For the task of sampling from a density $\pi \propto \exp(-V)$ on $\mathbb{R}^d$, where $V$ is possibly non-convex but $L$-gradient Lipschitz, we prove that averaged Langevin Monte Carlo outputs a sample with $\varepsilon$-relative Fisher…
Asymptotics for Dickman's number theoretic function $\rho(u)$, as $u \rightarrow \infty$, were given de Bruijn and Alladi, and later in sharper form by Hildebrand and Tenenbaum. The perspective in these works is that of analytic number…
Exerting a nonequilibrium drive on an otherwise equilibrium Langevin process brings the dynamics out of equilibrium but can also speedup the approach to the Boltzmann steady-state. Transverse forces are a minimal framework to achieve…
We study the asymptotic behaviour of a real-valued diffusion whose non-regular drift is given as a sum of a dissipative term and a bounded measurable one. We prove that two trajectories of that diffusion converge a.s. to one another at an…
EEG time series are analyzed using the diffusion entropy method. The resulting EEG entropy manifests short-time scaling, asymptotic saturation and an attenuated alpha-rhythm modulation. These properties are faithfully modeled by a…
In this paper, we introduce a novel Extra-Gradient method with anchor term governed by general parameters. Our method is derived from an explicit discretization of a Tikhonov-regularized monotone flow in Hilbert space, which provides a…
Recent progress on the theory of variational hypocoercivity established that Randomized Hamiltonian Monte Carlo -- at criticality -- can achieve pronounced acceleration in its convergence and hence sampling performance over diffusive…
Sampling from probability distributions of the form $\sigma \propto e^{-\beta V}$, where $V$ is a continuous potential, is a fundamental task across physics, chemistry, biology, computer science, and statistics. However, when $V$ is…
We introduce an autonomous system with closed-loop damping for first-order convex optimization. While, to this day, optimal rates of convergence are almost exclusively achieved by non-autonomous methods via open-loop damping (e.g.,…
Efficient estimation of a non-Gaussian stable Levy process with drift and symmetric jumps observed at high frequency is considered. For this statistical experiment, the local asymptotic normality of the likelihood is proved with a…
The asymptotic expansion of digamma function is a starting point for the derivation of approximants for harmonic sums or Euler-Mascheroni constant. It is usual to derive such approximations as values of logarithmic function, which leads to…
We study the long time asymptotics of the relaxation dynamics of the totally asymmetric simple exclusion process on a ring. Evaluating the asymptotic amplitudes of the local currents by the algebraic Bethe ansatz method, we find the…
We study sampling from posterior distributions with nonsmooth composite potentials, a setting in which proximal-based Langevin methods are theoretically appealing but in practice limited to simple functions with closed-form proximal…
Suppose that A_1,\dots, A_N are independent random matrices whose atoms are iid copies of a random variable \xi of mean zero and variance one. It is known from the works of Newman et. al. in the late 80s that when \xi is gaussian then…
We show that the random point measures induced by vertices in the convex hull of a Poisson sample on the unit ball, when properly scaled and centered, converge to those of a mean zero Gaussian field. We establish limiting variance and…