We use a two-site variant of this algorithm, which allows the bond dimension to grow with the evolution time at the expense of forgoing strict energy conservation owing to the truncation step involved ...