Approximate conservation laws in perturbed integrable lattice models

We develop a numerical algorithm for identifying approximately conserved quantities in models perturbed away from integrability. In the long-time regime, these quantities fully determine correlation functions of local observables. Applying the algorithm to the perturbed XXZ model, we find that the main effect of perturbation consists in expanding the support of conserved quantities. This expansion follows quadratic dependence on the strength of perturbation. The latter result, together with correlation functions of conserved quantities obtained from the memory function analysis, confirms the feasibility of the perturbation theory.