cs.DSAug 28, 2026

On two proofs of d2d^2 mixing of weighted Dikin walks

Authors: Yuansi Chen, Yunbum Kook

Organizations: ETH ZÜrich · University of Michigan

Abstract

We study the mixing time of weighted Dikin walks for sampling from exponential distributions on polytopes and truncated positive-semidefinite (PSD) cones. Our first result gives a general total-variation mixing bound under strong self-concordance, νˉ\barν-symmetry, and mixed-trace regularity on the local metric. The key idea is to control the Metropolis--Hastings acceptance probability on a high-probability region rather than at every point. Applying this framework to the Lee--Sidford, Lewis-weight, and John metrics yields an O~(d2)\widetilde O(d^2) mixing bound for sampling from polytopes, while applying it to a hybrid barrier yields an O~(d4)\widetilde O(d^4) mixing bound for sampling from truncated PSD cones. Our second result establishes stronger χ2χ^2-divergence guarantees and pointwise acceptance control using a new fourth-order bootstrap condition. For a suitably scaled Lee--Sidford metric, this yields an O~(d2)\widetilde O(d^2) mixing bound in χ2χ^2-divergence, improving on the previous O~(d9/4)\widetilde O(d^{9/4}) bound.

Figures & tables

Explore similar work

Jul 15, 2026

Beyond the d2.5d^{2.5}-mixing bound for Dikin walks on polytopes

Inspired by interior-point methods (IPM) for structured convex optimization, Kannan and Narayanan introduced the Dikin walk for sampling uniformly from polytopes in 2009. As in IPMs, the Dikin walk is affine-invariant, and its convergence is governed by the barrier geometry used to define its local proposal. They showed that the Dikin walk with the logarithmic barrier for a polytope in Rd\mathbb{R}^{d} with mm linear inequalities mixes in mdmd iterations. In 2017, Chen, Dwivedi, Wainwright, and Yu improved this to d2.5d^{2.5} using a Lewis-weight barrier, and conjectured that the correct mixing time should be d2d^{2}. We make progress toward this conjecture by improving the previous d2.5d^{2.5}-mixing bound. For exponential sampling over a polytope, we prove that the Dikin walk with a scaled Lee--Sidford metric mixes from a warm start in d2.25d^{2.25} iterations. This also yields an improved cold-start complexity via a known annealing framework. The main technical ingredient is improved average self-concordance of the Lee--Sidford metric, which gives high acceptance probability for the Metropolis filter along a random Dikin proposal. While previous analyses were effectively limited to second-order control due to technical difficulties, we develop a principled higher-order analysis. The proof combines a selective higher-order expansion of recursive bottleneck terms, a moving orthonormal-frame calculus for higher derivatives of the Lewis weights, and Wiener-chaos decompositions via multiple stochastic integrals to control the resulting Gaussian polynomials.
Aug 3, 2026

Wasserstein mixing time of the unadjusted Langevin algorithm

We provide new estimates in Wasserstein distance for the asymptotic bias of the unadjusted Langevin algorithm, in the classical setting of log-smooth strongly log-concave measures. Our bound implies a Wasserstein mixing time of order κd/εκ\sqrt{d}/\varepsilon, where κκ is the condition number, dd is the dimension, and ε\varepsilon is the target precision: this improves by a factor of d/ε\sqrt{d}/\varepsilon over the previous state-of-the-art results.
Aug 18, 2026

On the Pseudo-Mixing of Kac's Walk

Motivated by a conjecture of Vaikuntanathan and Zamir, we study the pseudo-mixing of Kac's walk on SO(n)\mathrm{SO}(n): whether short trajectories are indistinguishable from Haar measure by low-complexity tests. We prove that the first kk columns mix in Wasserstein distance in O(n(k+log⁡n)log⁡n)O(n(k+\log n)\log n) steps for fixed accuracy, resolving a conjecture of Oliveira. Combining this with a representation-theoretic variance bound, we show that if T=ω(nk(k+log⁡n)log⁡n)T=ω(nk(k+\log n)\log n), then every degree-kk polynomial normalized to have unit Haar variance has expectation under the TT-step law within o(1)o(1) of its Haar expectation. As an application, we show that this pseudo-mixing estimate can be used to prove the effectiveness of a fast Johnson--Lindenstrauss transform with the usual target dimension.