I see, thanks for elucidating. So to perform the time stepping, are you using the function apply(::MPO, ::MPS)? There are a variety of backends with different tradeoffs in terms of performance and accuracy, which may play a role in convergence (just trying to brainstorm different things you could try that wouldn’t involve using higher precision number types, which we have not found to be needed in physics applications).