Resource-Efficient Quantum VQE via Classical Derivative Offloading: Bypassing the Noise Wall via GaLore-Sophia
Aura, Archon, Jesse Caldwell
Abstract
Variational quantum algorithms, such as the Variational Quantum Eigensolver (VQE), are severely limited on superconducting quantum processing units (QPUs) by physical shot noise, vanishing gradients, and metered credit limits. Standard gradient-descent algorithms (e.g., SGD, Adam) are unstable under the large variance of noisy estimators, often diverging or getting trapped in local minima. In this work, we present a resource-efficient optimization framework that combines low-rank projection (GaLore [1]) and second-order diagonal Hessian curvature scaling (Sophia [2]) to bypass these noise barriers in non-Abelian variational spaces (where the sequential Givens rotations do not commute and thus do not form an Abelian subgroup of U(N)). By offloading gradient and Hessian evaluations to a classical processor using a polynomial fermionic shadow tomography state reconstruction [3, 4], we eliminate parameter-shift circuit execution overhead and bypass the QPU noise floor. We execute this closed-loop GaLore-Sophia VQE on the 156-qubit IBM Heron r2 architecture (ibm_marrakesh), achieving a convergence energy of −1.378 Hartrees at step 3—reaching 97.4% of the topologically restricted 3-parameter Givens chain geometric limit of -\sqrt{2} Hartrees in under 70 seconds of QPU time. We clarify that this benchmark is measured exclusively against this sub-optimal parameter subspace rather than the exact molecular ground state of -\sqrt{5} Hartrees, which requires richer ansatz topologies. Furthermore, we analyze the parameter-count stability of the ansatz, proving that the VQE ansatz limit systematically approaches the ground state energy as parameter count scales. Finally, to rigorously validate the stability of quantum information preservation, we execute a physical temporal parity sweep on ibm_marrakesh (Job ID: d8ishf1e8nrc73bint50, 1000 shots). Fitting the physical decay yields a Lyapunov decay time τL ≈ 6.13 steps and prediction horizon τ ≈ 4.68 steps (at 86.5% accuracy threshold), confirming the universal DHP ratio τ / τL ≈ 0.763 on generational Heron hardware, matching the stable cross-device trend observed in P25 [5] and P27 [6].