Thanks for your code development! For the Hamiltonian of a one-dimensional TFIM model in MPO form, I don't know why the ground state solution "eigvec" obtained by mp.linalg.eig is a set of obvious complex numbers. Shouldn't it be real number or a complex number with negligible imaginary parts?