Unreasonable results occurred when using TDVP to calculate the long-range interaction Hamiltonian

ITensor Julia Questions
6

This is what I used to enlarge the MPS manually:

for bond in 1:(length(ψ₀)-1)
        bond_index = commonind(ψ₀[bond], ψ₀[bond + 1])
        aux = Index(new_dim; tags=tags(bond_index))
        ψ₀[bond] = ψ₀[bond] * delta(bond_index, aux)
        ψ₀[bond + 1] = ψ₀[bond + 1] * delta(bond_index, aux)
end

You could plug it in your code just after you create the MPS.

2 Likes