Efficient computation of cumulant evolution and full counting statistics: application to infinite temperature quantum spin chains

We propose a numerical method to efficiently compute quantum generating functions (QGF) for a wide class of observables in one-dimensional quantum systems at high temperature. We obtain high-accuracy estimates for the cumulants and reconstruct full counting statistics from the QGF. We demonstrate its potential on spin $S=1/2$ anisotropic Heisenberg chain, where we can reach time scales hitherto inaccessible to state-of-the-art classical and quantum simulations. Our results challenge the conjecture of the Kardar--Parisi--Zhang universality for isotropic integrable quantum spin chains.