We establish the sharp logarithmic order $(\log N)^{-1/2}$ for the expected $p$-Wasserstein distance, induced by the supremum norm, between the empirical law of $N$ independent copies of a continuous It\^o process and their common path law. We only assume that the initial condition and the drift and diffusion integrands are controlled by a time-uniform random upper bound with a finite $\rho$-moment for some $\rho>p\geq1$. Under this assumption, we use an adaptive random time interval partition argument, which leads to a $(\log n)^{-1/2}$ functional quantization rate. A general transfer principle then converts the quantization estimate into a mean estimate and nonasymptotic deviation bounds for equal-weight empirical laws. Applications include empirical path-law estimates for path-dependent SDEs and a path-space propagation-of-chaos estimate for path-dependent McKean--Vlasov interacting particle systems.
We study the long-time behavior of the Wasserstein gradient flow of the squared Maximum Mean Discrepancy (MMD) between a probability measure $\rho$ and a target measure $\mu$, where the underlying kernel is given by a Coulomb potential. For $L^\infty$ target densities $\mu$, we establish the existence of global weak solutions starting from arbitrary Borel probability measures and prove that the density $\rho_t$ belongs to $L^\infty$ for any $t>0$. We also show that the H\"older norm can grow exponentially in time. On the flat torus ${\mathbb{T}}^\mathsf{d}$, we prove a global metric PL inequality for every finite-Coulomb-energy source and nearly uniform target. For general bounded, uniformly positive targets, we prove exponential decay of the squared MMD without requiring a lower bound on the initial data, using a defective PL inequality. We also prove that the usual PL inequality may fail when the target vanishes only at one point and that, when $\mathsf{d}\ge2$, no PL constant can hold uniformly over all targets satisfying a prescribed lower bound. On ${\mathbb{R}}^\mathsf{d}$, for $\mathsf{d}\ge2$, under radial symmetry, source-support inclusion, and target-positivity assumptions, we establish a PL inequality and exponential convergence. On the unrestricted whole-space class, neither a multiplicative squared-MMD decay modulus uniform over the initial datum nor a global PL inequality can hold. Finally, in every dimension and in both spatial settings, we prove that every Lagrangian critical point coincides with the target when $(\rho-\mu)^+$ is absolutely continuous. In dimension two, the energy supplies uniform tightness. This implies that if our constructed solutions have finite energy at some positive time, then they converge to the target narrowly and strongly in negative-order Sobolev spaces.
Antonin Chodron de Courcel, Matthew Rosenzweig· 0 citations
We consider a fractional Brownian motion $B$ with Hurst index $0<H<1/2$, and its maximiser $\tau$ on $[0,1]$. We show that the rescaled process $a^H(B_{\tau+\,\cdot\,/a}-B_\tau)$ converges in $C_{\mathrm{loc}}(\mathbb R)$ to a limiting tangent law that is $H$-self-similar, supported on nonpositive paths pinned at zero, and rerooting-rescaling invariant: rerooting the limit process at its maximum on any fixed compact interval separated from zero and rescaling again asymptotically reproduces the same law. We also identify the tangent law as the limit of two-sided finite-grid hard-wall laws as the mesh vanishes and both horizons diverge, which can informally be interpreted as conditioning fractional Brownian motion on a nonpositive path. As an application, we consider persistence probabilities for fractional Brownian motion: a tilted variant of $B$ yields a different tangent law with a finite left horizon and an infinite right horizon and we show that \[ \mathbb P(B_t\leq1\text{ for all }0\leq t\leq T) = \big(C+o(1)\big)\,T^{-(1-H)},\quad \text{as }T\to\infty, \] where the leading order coefficient $C\in(0,\infty)$ has an explicit representation in terms of the expected maximum and the tilted tangent law.
We prove absolute continuity of the running maximum $X^{\ast}_T=\sup_{0\leq s\leq T}X_s$ of one-dimensional time-inhomogeneous L\'evy--It\^o diffusions driven by a Brownian motion and an independent non-truncated pure-jump L\'evy process. Using Bismut's directional Malliavin calculus on the Wiener--Poisson space together with the running-maximum criteria of Song--Xie and Nakagawa--Suzuki, we reduce the problem to constructing an admissible direction $\Theta$. The key is to ensure that the directional derivative $D_\Theta X_t$ is strictly positive for all $t\in(0,T]$. We give explicit directions in two regimes. In the uniformly elliptic case with time-inhomogeneous coefficients, a purely Brownian perturbation yields an explicit positive integral representation for $D_\Theta X_t$, and hence $X^{\ast}_T$ admits a density without truncating the jump component. In a Brownian-degenerate pure-jump model with a bounded deterministic time-dependent jump weight $\kappa(t)$ that may vanish on subintervals, we prove absolute continuity under the minimal nondegeneracy-in-time condition $\int_0^t \kappa(s)^2\,ds>0$ for every $t>0$ and infinite activity of the L\'evy measure. A weighted Poisson positivity lemma is the key new input. Finally, we show that silent initial intervals can create atoms and we derive an explicit atom--density decomposition.
We study exact simulation of diffusions via rejection sampling on path space using unbiased estimators of the density ratio obtained from Girsanov's theorem. When applied to the underdamped Langevin diffusion, it yields an algorithm for sampling from a strongly log-concave and log-smooth distribution with condition number $\kappa$, in dimension $d$, to accuracy $\varepsilon$ in R\'enyi divergence, in $\widetilde O(\kappa^{2/3} d^{1/3}\,\mathrm{polylog}(1/\varepsilon))$ queries. Under a third derivative bound, the dimension dependence improves to $d^{1/5}$. This improves substantially over the prior state-of-the-art complexity of $\widetilde O(\kappa d^{1/2}\,\mathrm{polylog}(1/\varepsilon))$ for the Metropolis-adjusted Langevin algorithm, and over the $d^{1/4}$ dimension dependence of Metropolized Hamiltonian Monte Carlo under the same third derivative bound. We also present applications to the mirror Langevin diffusion, and for obtaining Fisher information bounds in the non-log-concave case.
Fan Chen, Sinho Chewi, Alexander Rakhlin et al.· 0 citations
We study a penalisation problem for two-dimensional Brownian motion. Starting from the Wiener measure, we consider a family of probability measures obtained by weighting paths by a nonnegative functional $F_t$ depending on $t \geq 0$, $F_t$ being measurable with respect to the $\sigma$-algebra generated by the path up to time $t$. Under suitable assumptions on the penalisation process, we establish the weak convergence of these measures when $t \rightarrow \infty$. The limiting law is identified explicitly in terms of a $\sigma-$finite measure $\mathbf{W}^{(2)}$, which admits a path decomposition involving the last hitting time of a circle. This decomposition plays a central role in the analysis and yields a martingale representation of the limiting measure. where ordering and local time techniques are no longer available. The proofs rely on Laplace transform methods and Tauberian theorems, which replace excursion-theoretic tools and allow a precise identification of the limiting measure and its structural properties.
Joseph Najnudel, Thammadol Tansrivorarat· 0 citations
This paper computes Poisson space noise functionals, $P^{\prime}(u)$, realised in a Gel'fand triple built from a L\'evy measure $\lambda_{\beta}(u)du$. We isolate three discretisation parameters: a small-amplitude cut-off, a Donsker delta truncation M, and a chaos order N. For a stable-type intensity $\lambda_{\beta}(u)=cu^{-1-\alpha}$ ($0<\alpha<2$), replacing discarded small amplitudes with matched Gaussian space noise improves the Wasserstein-1 error from $O(\epsilon^{1-\alpha/2})$ to $O(\epsilon)$. The residual is asymptotically normal at $O(\epsilon^{\alpha/2})$. This compensation reduces computational complexity from $O(\tau^{-2\alpha/(2-\alpha)})$ to $O(\tau^{-\alpha})$. We also evaluate the Gamma-type boundary ($\alpha=0$) and exponential tempering. Truncations converge algebraically (M) and super-geometrically (N). All predicted rates are tightly confirmed by deterministic numerical experiments via Gil-Pelaez inversion, eliminating Monte Carlo noise.