Articles
Here is a commented bibliography on unbiased MCMC.
For a recent review on the topic, see Atchadé and Jacob (2024), Unbiased Markov Chain Monte Carlo: what, why and how, available on arXiv, a version of which will hopefully be a chapter in the upcoming 2nd edition of the Handbook of Markov chain Monte Carlo.
Core ideas of unbiased estimators using couplings
Glynn and Rhee (2014) : first paper presenting the idea that contractive couplings of Markov chains, combined with a telescopic sum argument, lead to implementable unbiased estimators of expectations with respect to the stationary distribution of these Markov chains.
Agapiou, Roberts, and Vollmer (2018) : paper applying these ideas to the setting of Markov chain Monte Carlo methods, with contractive couplings.
Vihola (2018) : presents a general unbiased scheme with an explicit expression for its variance.
Pierre E. Jacob, O’Leary, and Atchadé (2020) and its discussion : paper applying these ideas with couplings that induce exact meetings, and with modifications that make the unbiased estimators nearly as efficient as the original MCMC estimators when tuned properly.
Convergence diagnostics
Biswas, Jacob, and Vanetti (2019) : develops an idea presented at the end of Pierre E. Jacob, O’Leary, and Atchadé (2020) to obtain bounds on the total variation distance and the 1-Wasserstein distance between a Markov chain at fixed time \(t\) and its limiting distribution.
Johnson (1996), Johnson (1998) : propose convergence diagnostics based on couplings, including couplings that induce exact meetings for random walk Metropolis-Rosenblut-Teller-Hastings algorithms.
Corenflos and Dau (2025) : proposes computable bounds on the f-divergence (e.g. KL, TV or \(\chi^2\)) between the target distribution \(\pi\) and the marginal distribution of the chain at time \(t\). The method involves running \(N\) pairs of chains up to time \(t\), coupled such that pairs can meet exactly at each iteration. Each chain is weighted and the weights are harmonized as the pairs coalesce. The matching between chains is randomized at each time step, so that all chains eventually interact with one another.
Normalizing constant estimation
- Rischard, Jacob, and Pillai (2018) : starts from thermodynamic integration, with a path of distributions \((\pi_\lambda)\) linking a tractable \(\pi_0\) and the target \(\pi\). Then randomizes \(\lambda\) and runs unbiased MCMC for \(\pi_\lambda\), to ultimately unbiased estimators of log ratios of normalizing constants.
Estimating nonlinear functions of expectations
- Wang and Wang (2022)
Specific classes of MCMC algorithms
Metropolis-Rosenbluth-Teller-Hastings
Wang, O’Leary, and Jacob (2021), O’Leary and Wang (2024) : construction of one-step maximal coupling of general MRTH kernels, and characterization of the set of all couplings of these kernels.
Papp and Sherlock (2024) : propose couplings of random walk proposal MRTH that are optimal for a class of targets in high-dimension, improving upon the reflection couplings employed in Pierre E. Jacob, O’Leary, and Atchadé (2020).
Deligiannidis et al. (2025) : the case of independent proposals, with a possibly unbounded ratio target/proposal. The common draws coupling of independent proposal MRTH is maximal. Connects with the idea of debiasing self-normalised importance sampling, as initiated in Middleton et al. (2019). Also looks at interconnection between unbiased MCMC and robust mean estimation.
Pseudo-marginal methods and conditional particle filters
Gradient-based MCMC
Gibbs sampling
Bacallado et al. (2021)
Nguyen, Trippe, and Broderick (2022)
Biswas et al. (2022)
Ceriani and Zanella (2024)
Ascolani and Zanella (2025) : Gibbs sampler high-dimensional probit regression. The paper obtains results about the mixing time of the sampler, not using couplings; but coupling-based total variation upper bounds are obtained to verify the theory in numerical experiments.
Du and He (2024)
Piecewise deterministic Markov chain Monte Carlo
- Corenflos, Sutton, and Chopin (2023)