hardware algorithms

Determinant Quantum-Quantum Monte Carlo: Coherent Auxiliary-Field Sampling

Curator's Take

AI Commentary

This article shows how a core workhorse of condensed‑matter simulation—determinant quantum Monte Carlo—can be lifted onto a quantum processor by encoding the single‑particle action matrix and using quantum singular‑value transforms to generate Hubbard‑Stratonovich weights without ever enumerating them. By eliminating the classical Markov chain and its autocorrelation time, the approach promises dramatically lower scaling of circuit depth with system size for fermionic models such as the Hubbard Hamiltonian, complementing recent quantum‑algorithmic advances that target many‑body ground states. The main caveat is a post‑selection cost that grows exponentially for sharply peaked weight distributions, so practical advantage will hinge on problem instances where the probability landscape remains smooth.

— Mark Eatherly

Summary

We introduce determinant quantum-quantum Monte Carlo (DQ$^2$MC), a quantum algorithm that lifts the auxiliary-field sampling and averaging at the operational core of determinant quantum Monte Carlo onto a quantum computer. A determinant oracle synthesizes the DQMC amplitudes directly from a block encoding of the single-particle action matrix via quantum singular value transformations, so that the exponentially many Hubbard-Stratonovich weights are never enumerated, precomputed, or stored. Since the fermions are free for fixed auxiliary fields, the construction operates entirely at the single-particle level, requiring $O(\log N_{\mathrm{st}})$ system qubits and no Jordan-Wigner or Bravyi-Kitaev encoding, where $N_{\mathrm{st}}$ is the space-time volume. A full-quantum protocol makes observables interference amplitudes, eliminating the Markov chain and its autocorrelation time altogether; a hybrid quantum-classical protocol retains a constant-size active block of qubits and replaces the Metropolis-Hastings acceptance step with an exact heat-bath draw, so that cluster updates of any size are rejection-free, and passes only classical information between updates, admitting parallel tempering and distributed execution across quantum processors. The circuit-depth scales more favorably with spatial volume than classical DQMC, at the price of a post-selection overhead determined exactly by the largest target probability --- polynomial for smooth distributions, exponential for sharply peaked ones. Finally, the reweighting estimator underlying the fermion sign problem maps exactly onto a quantum weak value, placing the exponential cost of sign-problematic DQMC in precise correspondence with the post-selection overhead of weak-value extraction.