Quantum Metropolis-Hastings algorithm with the target distribution calculated by quantum Monte Carlo integration

Abstract

The Markov chain Monte Carlo method (MCMC), especially the Metropolis-Hastings (MH) algorithm, is a widely used technique for sampling from a target probability distribution P on a state space and applied to various problems such as estimation of parameters in statistical models in the Bayesian approach. Quantum algorithms for MCMC have been proposed, yielding the quadratic speedup with respect to the spectral gap compered to classical counterparts. In this paper, we consider the quantum version of the MH algorithm in the case that calculating P is costly because the log-likelihood L for a state x∈ is obtained via computing the sum of many terms 1MΣi=0M-1 (i,x). We propose calculating L by quantum Monte Carlo integration and combine it with the existing method called quantum simulated annealing (QSA) to generate the quantum state that encodes P in amplitudes. We consider not only state generation but also finding a credible interval for a parameter, a common task in Bayesian inference. In the proposed method for credible interval calculation, the number of queries to the quantum circuit to compute scales on , the required accuracy ε and the standard deviation σ of as O(σ/ε23/2), in contrast to O(M/ε1/2) for QSA with L calculated exactly. Therefore, the proposed method is advantageous if σ scales on M sublinearly. As one such example, we consider parameter estimation in a gravitational wave experiment, where σ=O(M1/2).

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…