Hi, I am trying to evaluate the infinite temperature expectation value of an operator at different times in a spin-1 chain, as in Fig. 1 of https://arxiv.org/abs/1904.04266.

The Hamiltonian contains three sites terms (Eq. 1 in the article), so I cannot use toExp() for the time evolution. Any idea of how I could overcome this problem?

I tried using the TDVP library but it time-evolves the MPS, it seems I need to use Heisenberg representation to time-evolve the MPO.

Best,

Dario