  We consider the spreading of a local operator $A$ in Euclidean time in
one-dimensional many-body systems with Hamiltonian $H$ by calculating the
$k$-fold commutator $[H,[H,[...,[H,A]]]]$. We derive general bounds for the
operator norm of this commutator in free and interacting fermionic systems with
and without disorder. We show, in particular, that in a localized system the
norm does grow at most exponentially and that the contributions of operators to
the total norm are exponentially suppressed with their length. We support our
general results by considering one specific example, the XXZ chain with random
magnetic fields. We solve the operator spreading in the XX case without
disorder exactly. For the Anderson and Aubry-Andr\'e models we provide strict
upper bounds. We support our results by symbolic calculations of the commutator
up to high orders. For the XXX case with random magnetic fields, these symbolic
calculations show a growth of the operator norm faster than exponential and
consistent with the general bound for a non-localized system. Also, there is no
exponential decay of the contribution of operators as function of their length.
We conclude that there is no indication for a many-body localization
transition. Finally, we also discuss the differences between the interacting
and non-interacting cases when trying to perturbatively transform the
microscopic to an effective Hamiltonian of local conserved charges by
consecutive Schrieffer-Wolff transformations. We find that such an approach is
not well-defined in the interacting case because the transformation generates
$\sim 4^\ell$ terms connecting sites a distance $\ell$ apart which can
overwhelm the exponential decay with $\ell$ of the amplitude of each individual
term.
