Abstract We exploit a suite of large N-body simulations (up to N = 40963) performed with abacus, scale-free models range spectral indices n, better understand and quantify convergence the matter power spectrum. Using self-similarity identify converged regions, we show that maximal wavenumber resolved at given level accuracy increases monotonically as function time. At 1 per cent it starts early...