Making an MPO from an MPS

Hi,

I have a question regarding creating an MPO rho from an MPS psi. The way I currently do this is rho = outer(psi’,psi) but this blows up when the bond dimension of psi is above 100. Is there anyway to do this more efficiently? Or should I just increase the truncation in outer to reduce the size of rho?

While I’m not sure that I know the perfect answer, it may depend on what you want to do after making the density matrix. What steps do you plan to take afterward? Will you be time evolving it using an open-system Lindbladian or something like that? Or are you just wanting to compute properties of the state?