Abstract:We consider the computational task of sampling a bit string $x$ from a distribution $\pi(x)=|\langle x|\psi\rangle|^2$, where $\psi$ is the unique ground state of a local Hamiltonian $H$. Our main result describes a direct link between the inverse spectral gap of $H$ and the mixing time of an associated continuous-time Markov Chain with steady state $\pi$. The Markov Chain can be implemented efficiently whenever ratios of ground state amplitudes $\langle y|\psi\rangle/\langle x|\psi\rangle$ are efficiently computable, the spectral gap of $H$ is at least inverse polynomial in the system size, and the starting state of the chain satisfies a mild technical condition that can be efficiently checked. This extends a previously known relationship between sign-problem free Hamiltonians and Markov chains. The tool which enables this generalization is the so-called fixed-node Hamiltonian construction, previously used in Quantum Monte Carlo simulations to address the fermionic sign problem. We implement the proposed sampling algorithm numerically and use it to sample from the ground state of Haldane-Shastry Hamiltonian with up to 56 qubits. We observe empirically that our Markov chain based on the fixed-node Hamiltonian mixes more rapidly than the standard Metropolis-Hastings Markov chain.

What problem does this paper attempt to address?

### Problems Addressed by the Paper This paper attempts to address the problem of generating samples from the ground state of quantum many-body systems. Specifically, given a local Hamiltonian $ H $ whose ground state $ \psi $ is unique and has an energy gap $ \gamma $ from the excited states, the goal is to sample bit strings $ x $ from the distribution $ \pi(x) = |\langle x | \psi \rangle|^2 $. ### Main Contributions 1. **Rapidly Mixing Markov Chain**: - The authors propose a continuous-time Markov chain (CTMC) whose stationary distribution is exactly the desired $ \pi(x) $. - By designing this Markov chain, the authors prove that its mixing time is related to the spectral gap $ \gamma $ of the Hamiltonian $ H $. Specifically, the mixing time $ T $ satisfies $ T \leq \tilde{O}(n k s / \gamma) $, where $ n $ is the number of qubits, $ k $ is the size of the local interactions, and $ s $ is a sensitivity parameter. 2. **Fixed-Node Hamiltonian**: - The authors introduce a construction method called the "fixed-node Hamiltonian" $ F $, which maps any local Hamiltonian $ H $ to a new Hamiltonian $ F $ such that $ F $ and $ H $ share the same ground state, and $ F $ is sign-problem-free. - This construction allows the use of relative phase information of the ground state amplitudes when computing probability ratios, thereby improving sampling efficiency. 3. **Algorithm Implementation**: - The authors implement the aforementioned continuous-time Markov chain using Gillespie’s algorithm and numerically validate its performance on the Haldane-Shastry Hamiltonian. - Experimental results show that the Markov chain based on the fixed-node Hamiltonian mixes faster than the standard Metropolis-Hastings Markov chain. ### Key Techniques - **Fixed-Node Hamiltonian**: By converting the original Hamiltonian $ H $ to a fixed-node Hamiltonian $ F $, the sign problem is resolved, making the sampling process more efficient. - **Continuous-Time Markov Chain**: Compared to discrete-time Markov chains, continuous-time Markov chains have advantages in handling large matrix elements, effectively simulating the sampling process. - **Gillespie’s Algorithm**: Used to implement the continuous-time Markov chain, ensuring the accuracy and efficiency of the sampling. ### Conclusion By introducing the fixed-node Hamiltonian and continuous-time Markov chain, the paper provides an efficient method for generating samples from the ground state of quantum many-body systems. This method is not only theoretically well-founded but also demonstrates good performance in practical applications.

A rapidly mixing Markov chain from any gapped quantum many-body system

Mixing time of quantum Gibbs sampling for random sparse Hamiltonians

Faster quantum mixing for slowly evolving sequences of Markov chains

Slow Mixing of Quantum Gibbs Samplers

Quasi-optimal sampling from Gibbs states via non-commutative optimal transport metrics

For every quantum walk there is a (classical) lifted Markov chain with faster mixing time

An efficient and exact noncommutative quantum Gibbs sampler

On Mixing of Markov Chains: Coupling, Spectral Independence, and Entropy Factorization

Quantum generalizations of Glauber and Metropolis dynamics

Quantum-enhanced Markov chain Monte Carlo

Rapidly mixing loop representation quantum Monte Carlo for Heisenberg models on star-like bipartite graphs

Bounding speedup of quantum-enhanced Markov chain Monte Carlo

Quantum Sampling Algorithms, Phase Transitions, and Computational Complexity

Quantum Sampling Algorithms for Near-Term Devices

Enabling Quantum Speedup of Markov Chains using a Multi-level Approach

Boson Sampling with Gaussian input states: efficient scaling and certification

Quantum-enhanced Markov Chain Monte Carlo for systems larger than your Quantum Computer

Stochastic Quantum Sampling for Non-Logconcave Distributions and Estimating Partition Functions

Locally Stationary Distributions: A Framework for Analyzing Slow-Mixing Markov Chains

Dissipative Quantum Gibbs Sampling

Quantum Metropolis Sampling via Weak Measurement