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.

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

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

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

Advantages of Unfair Quantum Ground-State Sampling

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

Sampling diverse near-optimal solutions via algorithmic quantum annealing

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

Optimization by Parallel Quasi-Quantum Annealing with Gradient-Based Sampling

Optimizing the Optimizer: Decomposition Techniques for Quantum Annealing

Scaling Advantage in Approximate Optimization with Quantum Annealing

Quantum Sampling Algorithms for Near-Term Devices

Landscape approximation of low energy solutions to binary optimization problems

borealis - A generalized global update algorithm for Boolean optimization problems

Optimization by quantum annealing: lessons from hard satisfiability problems

Reverse quantum annealing approach to portfolio optimization problems

Optimizing Schedules for Quantum Annealing

Quantum Sampling Algorithms, Phase Transitions, and Computational Complexity

Hybrid Quantum Annealing via Molecular Dynamics

Quantum-Enhanced Greedy Combinatorial Optimization Solver

Building an iterative heuristic solver for a quantum annealer

Non-Equilibrium Dynamics of Hybrid Continuous-Discrete Ground-State Sampling