Mixing times of a Burnside process Markov chain on set partitions
J. E. Paguyo
TL;DR
This work analyzes the Burnside process on $[k]^n$ with $S_k$ acting by permuting coordinate values, showing that the induced chain on set partitions $\\\\ $ mixes rapidly in two regimes: $k\ge n$ and $1\le k<n$. Using coupling and minorization, the authors establish explicit mixing-time bounds: in the first regime, $t_{mix}(\\varepsilon)\le\\lceil 2k\\log(n/\\varepsilon)\\rceil$, and in the second, $t_{mix}(\\varepsilon)\le\\lceil (k-1)! \\log(1/\\varepsilon)\\rceil$, with independence from $n$ for fixed $k$. They also derive explicit transition probabilities both for the original process on $[k]^n$ and its lumped chain on partitions, and obtain eigenvalue bounds for the lumped and original chains, linking spectral properties to mixing via relaxation times. These results provide a rigorous foundation for using the Burnside process to uniformly sample set partitions and offer insights into the chain’s spectral structure, while highlighting open problems in diagonalizing the transition operators and proving sharp lower bounds. The work thus contributes practical mixing-time guarantees for an orbit-based sampler and connects combinatorial structure with Markov-chain techniques such as coupling, minorization, and geometric bounds.
Abstract
Let $X$ be a finite set and let $G$ be a finite group acting on $X$. The group action splits $X$ into disjoint orbits. The Burnside process is a Markov chain on $X$ which has a uniform stationary distribution when the chain is lumped to orbits. We consider the case where $X = [k]^n$ with $k \geq n$ and $G = S_k$ is the symmetric group on $[k]$, such that $G$ acts on $X$ by permuting the value of each coordinate. The resulting Burnside process gives a novel algorithm for sampling a set partition of $[n]$ uniformly at random. We obtain bounds on the mixing time and show that the chain is rapidly mixing. For the case $k < n$, the algorithm corresponds to sampling a set partition of $[n]$ with at most $k$ blocks, and we obtain a mixing time bound which is independent of $n$. Along the way, we obtain explicit formulas for the transition probabilities and bounds on the second largest eigenvalue for both the original process and the lumped chain.
