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.
The easiest way I could think of would be to use the TEBD algorithm, which would be easy enough in the C++ or Julia versions. Have you looked into that?
You can look here: