So, I want to ask a follow-up question about this function. Is it possible to create a thermal state using that method ? that mean, if dt is the same and I am just changing the number of the steps I am taking.
site=siteinds("Qubit",n)
psi= randomITensor(site)
dt=0.1
num=t/dt
for j=1:num
psi= tdvp(H,psi,-dt)
normalize!(psi)
end
I know from reading that this method is usually good for getting the ground state in the limit where n dt \rightarrow \infty, but I am asking if in ITensor this thing can go around, and get your thermal state for ndt =t