MPO for Anyon Hubbard model

The graph you’re showing is obtained using Open Boundary Condition (OBC), while you are using PBC, could that be the problem? Otherwise, can you show how you are taking the expectation values?