For a transverse-field Ising chain with weak long-range interactions we develop a perturbative scheme, based on quantum kinetic equations, around the integrable nearest-neighbour model. We introduce, discuss, and benchmark several truncations of the time evolution equations up to eighth order in the Jordan-Wigner fermionic operators. The resulting set of differential equations can be solved for lattices with O(10 2 ) sites and facilitates the computation of spin expectation values and correlation functions to high accuracy, at least for moderate timescales. We use this scheme to study the relaxation dynamics of the model, involving prethermalisation and thermalisation. The techniques developed here can be generalised to other spin models with weak integrability-breaking terms.

No abstract