Abstract:Optimization problems pose challenges across various fields. In recent years, quantum annealers have emerged as a promising platform for tackling such challenges. To provide a new perspective, we develop a heuristic tensor-network-based algorithm to reveal the low-energy spectrum of Ising spin-glass systems with interaction graphs relevant to present-day quantum annealers. Our deterministic approach combines a branch-and-bound search strategy with an approximate calculation of marginals via tensor-network contractions. Its application to quasi-two-dimensional lattices with large unit cells of up to 24 spins, realized in current quantum annealing processors, requires a dedicated approach that utilizes sparse structures in the tensor network representation and GPU hardware acceleration. We benchmark our approach on random problems defined on Pegasus and Zephyr graphs with up to a few thousand spins, comparing it against the D-Wave Advantage quantum annealer and Simulated Bifurcation algorithm, with the latter representing an emerging class of classical Ising solvers. Apart from the quality of the best solutions, we compare the diversity of low-energy states sampled by all the solvers. For the biggest considered problems with over 5000 spins, the state-of-the-art tensor network approaches lead to solutions that are $0.1\%$ to $1\%$ worse than the best solutions obtained by Ising machines while being two orders of magnitude slower. We attribute those results to approximate contraction failures. While all three methods can output diverse low-energy solutions, e.g., differing by at least a quarter of spins with energy error below $1\%$, our deterministic branch-and-bound approach finds sets of a few such states at most. On the other hand, both Ising machines prove capable of sampling sets of thousands of such solutions.

Uncertain fate of fair sampling in quantum annealing

Exponentially Biased Ground-State Sampling of Quantum Annealing Machines with Transverse-Field Driving Hamiltonians.

Advantages of Unfair Quantum Ground-State Sampling

Fair Sampling of Ground-State Configurations of Binary Optimization Problems.

Sampling diverse near-optimal solutions via algorithmic quantum annealing

Biased Degenerate Ground-State Sampling of Small Ising Models with Converged QAOA

Best-case Performance of Quantum Annealers on Native Spin-Glass Benchmarks: How Chaos Can Affect Success Probabilities

Seeking Quantum Speedup Through Spin Glasses: the Good, the Bad, and the Ugly

Statistical Analysis of Quantum Annealing

Scaling Advantage in Approximate Optimization with Quantum Annealing

Quantum Annealing with chaotic driver Hamiltonians

Experimental demonstration of perturbative anticrossing mitigation using non-uniform driver Hamiltonians

Monte Carlo sampling from the quantum state space. II

Quantum sampling on a quantum annealer for large volumes in the strong coupling limit for gauge group U(3)

Nonstoquastic Hamiltonians and quantum annealing of an Ising spin glass

Can quantum Monte Carlo simulate quantum annealing?

Catastrophic failure of quantum annealing owing to non-stoquastic Hamiltonian and its avoidance by decoherence

Limitations of tensor network approaches for optimization and sampling: A comparison against quantum and classical Ising machines

Quantum versus classical annealing of Ising spin glasses

Quantum Sampling Algorithms, Phase Transitions, and Computational Complexity

Boltzmann sampling with quantum annealers via fast Stein correction