canonical setting

Strongly log-concave + log-smooth

π∝e−V\pi\propto e^{-V}, αI⪯∇2V⪯βI\alpha I\preceq\nabla^2V\preceq\beta I, κ=β/α\kappa=\beta/\alpha.

Upper and lower bounds

Comparison table

13 scoped results. Tildes suppress logarithmic factors.

Edit table
Best-known and comparison results for Strongly log-concave + log-smooth
ResultAlgorithm or modelComplexityGuaranteeOracle / startAssumptions and notesReviewSources
best upperGaussian cloud samplerO~(κd1/5 polylog⁡(1/ε))\widetilde O(\kappa d^{1/5}\,\operatorname{polylog}(1/\varepsilon))∥μN−π∥TV≤ε\lVert\mu_N-\pi\rVert_{\rm TV}\le\varepsilon∇V\nabla V (expected queries)Given xrefx_{\rm ref} with ∥∇V(xref)∥≤αd\lVert\nabla V(x_{\rm ref})\rVert\le\sqrt{\alpha d}; no warm start0<ε≤1/40<\varepsilon\le1/4; only needs β\beta-smoothness and an LSI with constant 1/α1/\alpha.Exact bound: O(κ{d+log⁡(κd/ε)}1/5log⁡13/2(κd/ε))O\bigl(\kappa\{d+\log(\kappa d/\varepsilon)\}^{1/5}\log^{13/2}(\kappa d/\varepsilon)\bigr), via the proximal sampler around a FORS-corrected smoothed Picard HMC. Ignoring logarithms, the best checked rate for κ≪d1/10\kappa\ll d^{1/10} and for d1/5≪κ≪d2/5d^{1/5}\ll\kappa\ll d^{2/5}; Picard warm start + Proximal BPS is better in between, and exact ULD / FORS for κ≫d2/5\kappa\gg d^{2/5}.Checkedpreprint
best upperPicard warm start + Proximal BPSO~((κ7/6d1/6+κ1/2d1/4) polylog⁡(1/ε))\widetilde O\bigl((\kappa^{7/6}d^{1/6}+\kappa^{1/2}d^{1/4})\,\operatorname{polylog}(1/\varepsilon)\bigr)∥μN−π∥TV≤ε\lVert\mu_N-\pi\rVert_{\rm TV}\le\varepsilon∇V\nabla V (expected queries)Given xrefx_{\rm ref} with ∥∇V(xref)∥≤αd\lVert\nabla V(x_{\rm ref})\rVert\le\sqrt{\alpha d}0<ε<1/40<\varepsilon<1/4; standard curvature bounds.Includes warm start generation; trades condition-number dependence for dimension dependence versus FORS.Checkedpreprint
best upperSmoothed Picard HMC + recursive warm startO~(κ7/6d1/6/ε1/3)\widetilde O(\kappa^{7/6}d^{1/6}/\varepsilon^{1/3})KL⁡(μN∥π)≤ε2\operatorname{KL}(\mu_N\Vert\pi)\le\varepsilon^2∇V\nabla V (expected queries)Given xrefx_{\rm ref} with V(xref)−inf⁡V≤dV(x_{\rm ref})-\inf V\le d0<ε≤1/20<\varepsilon\le1/2; standard curvature bounds.Generates its own initialization. The additive κ\kappa term is dominated; there is no κ2\kappa^2 term. Improves on FORS when d≳κ3/ε2d\gtrsim\kappa^3/\varepsilon^2, ignoring logarithms.Checkedpreprint
provisional upperProximal sampler + Picard warm start + Proximal BPSO~(κd3/16 polylog⁡(1/ε))\widetilde O(\kappa d^{3/16}\,\operatorname{polylog}(1/\varepsilon))∥μN−π∥TV≤ε\lVert\mu_N-\pi\rVert_{\rm TV}\le\varepsilon∇V\nabla V (expected queries)Known mode / Gaussian cold start; inner reference points computed by convex optimizationκ≳d1/8\kappa\gtrsim d^{1/8}; approximate RGO error budgeted across outer steps.Derived bound, not a stated theorem: inner condition number t≍d1/8t\asymp d^{1/8} gives (κ/t) (t7/6d1/6+t1/2d1/4)=O(κd3/16)(\kappa/t)\,\bigl(t^{7/6}d^{1/6}+t^{1/2}d^{1/4}\bigr)=O(\kappa d^{3/16}). Improves the listed high-accuracy rates for d1/8≪κ≪d7/16d^{1/8}\ll\kappa\ll d^{7/16}, ignoring logarithms; no matching general lower bound.AI-only noteno primary source
best upperExact ULD / FORSO~(κ2/3d1/3 polylog⁡(1/ε))\widetilde O(\kappa^{2/3}d^{1/3}\,\operatorname{polylog}(1/\varepsilon))Rq≤ε2\mathcal R_q\le\varepsilon^2V,∇VV,\nabla VCold start under the paper's setupExact path-space rejection implementation; fixed Rényi order.Complementary high-accuracy frontier; theorem-level transcription still needs human checking.Citedpreprint
best lowerGeneral first-order oracle lower boundΩ~(min⁡{κlog⁡d,d})\widetilde\Omega(\min\{\sqrt\kappa\log d,d\})Constant total-variation accuracyV,∇VV,\nabla V queriesAny randomized algorithm in the stated modelGaussian subclass; parameter ranges in the theorem.Fixed d≥2d\ge2 also has the tight Θd(log⁡κ)\Theta_d(\log\kappa) dependence.Checkedpublished
upper/lowerLMC / ULAUpper: O~(κd/ε)\widetilde O(\kappa\sqrt d/\varepsilon); lower: Ω~(κ+κd/ε)\widetilde\Omega\bigl(\kappa+\sqrt{\kappa d}/\varepsilon\bigr)α W2(μN,π)≤ε\sqrt\alpha\,W_2(\mu_N,\pi)\le\varepsilon∇V\nabla V; fixed-step LMCAny finite-second-moment law; initial Wasserstein distance enters logarithmically0<ε≤10<\varepsilon\le1; standard curvature bounds only.Lower bound: anisotropic Gaussian targets with one curvature-α\alpha direction and d−1d-1 curvature-β\beta directions, for d≥2d\ge2 and κ≥2\kappa\ge2, under worst-case initialization. Stability and slow-direction relaxation require Ω(κ)\Omega(\kappa) steps; stationary bias and relaxation give Ω~(κd/ε)\widetilde\Omega(\sqrt{\kappa d}/\varepsilon).Checkedpreprint
upper/lowerMALA (cold start)Θ~(κd)\widetilde\Theta(\kappa d)Constant TV accuracy; upper bound extends to high accuracy with polylogarithmic overheadV,∇VV,\nabla V; fixed-step MALAUpper: N(x⋆,β−1I)\mathcal N(x_\star,\beta^{-1}I); lower: a cold law with χ2(μ0∥π)≤ed\chi^2(\mu_0\Vert\pi)\le e^dKnown mode for the upper bound; lower bound is worst-case over cold initial laws.Matching dependence on dimension and condition number in the cold-start analysis; the lower bound does not assert hardness for the particular Gaussian upper-bound initialization.Checkedmonograph
upperULMC warm start + MALAO~(κd polylog⁡(1/ε))\widetilde O(\kappa\sqrt d\,\operatorname{polylog}(1/\varepsilon))∥μN−π∥TV≤ε\lVert\mu_N-\pi\rVert_{\rm TV}\le\varepsilon∇V\nabla V for ULMC; V,∇VV,\nabla V for MALAKnown mode x⋆x_\star; Gaussian cold start with R2(μ0∥π)=O~(d)\mathcal R_2(\mu_0\Vert\pi)=\widetilde O(d)Standard strong convexity and smoothness only.ULMC generates R2=O(1)\mathcal R_2=O(1), hence χ2=O(1)\chi^2=O(1), in O~(κd)\widetilde O(\kappa\sqrt d) queries. The displayed cost includes initialization and the subsequent MALA run.Checkedmonograph
upper/lowerRandomized / Poisson midpoint ULMCUpper: RM O~(κ5/6d1/3/ε2/3)\widetilde O(\kappa^{5/6}d^{1/3}/\varepsilon^{2/3}); Poisson (p=3p=3) O~(κ13/10d11/30/ε1/3)\widetilde O(\kappa^{13/10}d^{11/30}/\varepsilon^{1/3}). Lower: Ω(d1/3/ε1/3)\Omega(d^{1/3}/\varepsilon^{1/3}) iterationsRM upper: ∥μN−π∥TV≤ε\lVert\mu_N-\pi\rVert_{\rm TV}\le\varepsilon; Poisson upper and lower bound: α W2(μN,π)≤ε\sqrt\alpha\,W_2(\mu_N,\pi)\le\varepsilon∇V\nabla V (Poisson: expected queries); one random quadrature time shared across coordinatesUpper: known mode x⋆x_\star; position N(x⋆,β−1I)\mathcal N(x_\star,\beta^{-1}I) and standard Gaussian momentum. Lower: position N(ed,I)\mathcal N(e_d,I), standard Gaussian momentumPoisson upper: 0<ε≤cmin⁡{d/κ,d/κ}0<\varepsilon\le c\min\{\sqrt{d/\kappa},d/\sqrt\kappa\}; Corollary 2 with p=3p=3. Lower: standard Gaussian target (κ=1\kappa=1), d≥2d\ge2, 0<ε≤ε00<\varepsilon\le\varepsilon_0, Poisson grid size k≥2k\ge2.Poisson accuracy is rescaled from W22≤ε2d/αW_2^2\le\varepsilon^2d/\alpha in the paper. Its Wasserstein guarantee and the RM TV guarantee are distinct; larger pp gives further tradeoffs. Lower bound: the shared quadrature time makes every coordinate's variance fluctuate together, giving stationary bias ≳min⁡{dh3,d h3/2}\gtrsim\min\{dh^3,\sqrt d\,h^{3/2}\}. It matches Poisson ULMC for fixed p>3p>3 at sufficiently small ε\varepsilon, and also holds for overdamped RLMC and PLMC; it does not bound the RM TV guarantee.Checkedmonograph
upper/lowerHeun (deterministic midpoint LMC)Upper: O~(κ4/3d1/3/ε2/3)\widetilde O(\kappa^{4/3}d^{1/3}/\varepsilon^{2/3}); lower: Ω~(κmin⁡{d1/3/ε2/3, d1/2/ε1/2})\widetilde\Omega\bigl(\kappa\min\{d^{1/3}/\varepsilon^{2/3},\,d^{1/2}/\varepsilon^{1/2}\}\bigr)α W2(μN,π)≤ε\sqrt\alpha\,W_2(\mu_N,\pi)\le\varepsilon∇V\nabla V; fixed-step Heun, two gradients per stepUpper: known mode, X0=x⋆X_0=x_\star; lower: a law with E[V(X0)−V(x⋆)]≤Cd\mathbb E[V(X_0)-V(x_\star)]\le CdUpper: 0<ε≤d/20<\varepsilon\le\sqrt d/2; lower: d≥d0d\ge d_0, κ≥κ0\kappa\ge\kappa_0, 0<ε≤ε00<\varepsilon\le\varepsilon_0. No third-derivative bounds.Lower bound: step-size-dependent oscillatory targets built from cosine-perturbed blocks coupled to a curvature-α\alpha coordinate. For d−1≤ε≤ε0d^{-1}\le\varepsilon\le\varepsilon_0 it reads Ω~(κd1/3/ε2/3)\widetilde\Omega(\kappa d^{1/3}/\varepsilon^{2/3}), matching the dd and ε\varepsilon powers with a κ1/3\kappa^{1/3} gap. Dominated by smoothed Picard HMC, whose KL guarantee gives α W2≤2 ε\sqrt\alpha\,W_2\le\sqrt2\,\varepsilon.Checkedpreprint
upper/lowerUnderdamped exponential midpointUpper: O~(κ5/4d1/4/ε1/2)\widetilde O(\kappa^{5/4}d^{1/4}/\varepsilon^{1/2}); lower: Ω~(κ+κ3/4d1/4/ε1/2)\widetilde\Omega\bigl(\kappa+\kappa^{3/4}d^{1/4}/\varepsilon^{1/2}\bigr)α W2(μN,π)≤ε\sqrt\alpha\,W_2(\mu_N,\pi)\le\varepsilon (position marginal)∇V\nabla V; fixed step, friction γ=2β\gamma=2\sqrt\beta, two gradients per stepUpper: known mode, X0=x⋆X_0=x_\star, P0∼N(0,I)P_0\sim\mathcal N(0,I); lower: X0=α−1/2edX_0=\alpha^{-1/2}e_d, P0∼N(0,I)P_0\sim\mathcal N(0,I)Upper: 0<ε≤d/20<\varepsilon\le\sqrt d/2; lower: d≥2d\ge2, κ≥2\kappa\ge2, 0<ε≤1/20<\varepsilon\le1/2. No third-derivative bounds.Lower bound: anisotropic Gaussian with d−1d-1 curvature-β\beta directions and one curvature-α\alpha direction; matches the dd and ε\varepsilon powers with a κ1/2\kappa^{1/2} gap. At bounded κ\kappa, beats the d1/3/ε1/3d^{1/3}/\varepsilon^{1/3} lower bound for every randomized midpoint method when ε≫d−1/2\varepsilon\gg d^{-1/2}. Dominated by smoothed Picard HMC, whose KL guarantee gives α W2≤2 ε\sqrt\alpha\,W_2\le\sqrt2\,\varepsilon.Checkedpreprint
upperBlock-Krylov Gaussian samplerO((d∧κ)log⁡(d/ε2))O((d\wedge\sqrt\kappa)\log(d/\varepsilon^2))KL⁡≤ε\sqrt{\operatorname{KL}}\le\varepsilonGaussian matrix-vector / gradientCentered Gaussian modelGaussian subclass.Nearly tight for Gaussians, not a general-target upper bound.Checkedmonograph

Scope

The default oracle returns V(x)V(x) and ∇V(x)\nabla V(x). Rows are not directly comparable unless the error metric, initialization, and additional smoothness assumptions match.

Conventions

O~\widetilde O and Ω~\widetilde\Omega suppress logarithmic factors. Rows use an explicit initialization or include the cost of generating one. Algorithm-specific matching bounds do not imply a general oracle lower bound.