Challenging the Quantum Advantage Frontier with Large-Scale Classical Simulations of Annealing Dynamics

Abstract

Recent demonstrations of D-Wave's annealing-based quantum simulators have established new benchmarks for quantum computational advantage [arXiv:2403.00910]. However, the precise location of the classical-quantum computational frontier remains an open question, as classical simulation strategies continue to evolve. Here, we demonstrate that time-dependent variational Monte Carlo (t-VMC) with a physically motivated Jastrow-Feenberg wave function can efficiently simulate the quantum annealing of spin glasses up to system sizes previously thought to be intractable. Our approach achieves accuracy comparable to that of quantum processing units while requiring only polynomially scaling computational resources, in stark contrast to entangled-limited tensor network methods that scale exponentially. For systems up to 128 spins on a three-dimensional diamond lattice, we maintain correlation errors below 7%, which match or exceed the precision of existing quantum hardware. Rigorous assessments of residual energies and time-dependent variational principle errors establish clear performance benchmarks for classical simulations. These findings substantially shift the quantum advantage frontier and underscore that classical variational techniques, which are not fundamentally constrained by entanglement growth, remain competitive at larger system sizes than previously anticipated.