Hi @tsegev, I don’t have a direct answer to your question (that seems to refer to ITensor under Julia) but incidentally, the time-evolution of the density matrix of a system of qubits is what the Python/C++ package lindbladmpo solves - you can find it directly on github, and a link exists also in the Community Codes section (I am co-developer of this package). Not sure that it’s helpful.
Haggai.
about my original question what I eventually did is to trace out all the prime indices of the MPO and after that make the new object an MPS.
But I am not sure that this is a good approach (or even correct) to solve my original problem
Hi, no, this package is not available for Julia (but it can be used without writing a single C++ line, because everything can be done from Python - within the boundaries of the supported simulations).
Best.