Dear ALPS Users and Developers,
I have been observing some unexpected behavior when using translation symmetry with the ALPS fulldiag application. The detailed description (PDF), as well as a simple example and the output in question (.zip), are attached.
More recently, I have discovered another issue: If I consider a single spin-1/2 in the Sz basis subjected to magnetic field pointing along the y-direction, the eigenvectors and energies are calculated correctly but the expectation value of e.g. the identity-operator is zero. I have tracked down the problem to a surplus complex conjugation in the calculation of the scalar product of complex vectors and the attached patch seems to fix it. Moreover, the originally encountered problem (cf. first paragraph) seems to be no more than a different manifestation of the same bug.
I hope this information will be useful.
Best regards,
Marek Pikulski Laboratory for Solid State Physics ETH Zurich