Real-Frequency Correlation Functions from Neural Quantum States via Operator Lanczos Approach

Dynamical correlation functions are essential for revealing emergent phenomena in quantum matter, since they measure the frequency-resolved response of a system at equilibrium to a perturbation and thereby probe its full spectrum. While neural quantum states provide an accurate and flexible description of correlated ground states, accessing real-frequency dynamical correlators from them remains a major challenge. Here we introduce a general and systematically improvable framework to compute real-frequency linear-response functions directly from neural quantum state ground states. Our method performs the Lanczos subspace construction through an operator recursion and controls operator growth with a suited truncation scheme. We benchmark it on demanding quantum impurity models, where high precision is required close to the Fermi energy, and show that rich spectral structure and the self-energy can be resolved over several orders of magnitude. The approach generalizes straightforwardly, for example providing a direct path toward more broadly applicable dynamical mean-field theory solvers.