Exact Solution for Non-Hermitian Free Fermions: A Case Study of the XY Chain

We consider the non-Hermitian XY spin chain with open boundary conditions when the anisotropy parameter is extended to complex values. By analyzing the quasi-Hamiltonian matrix, we demonstrate that the free-fermion structure of the quasi-energy spectrum coincides with that of the Hermitian model and construct the corresponding biorthogonal fermionic basis away from exceptional points (EPs). We make use of an explicit Chebyshev-polynomial representation of the open-boundary eigenvectors in which the quasi-energy $\varepsilon$ is the natural spectral variable. This quasi-energy polynomial form is particularly useful at EPs, because EPs correspond to repeated roots of the same boundary polynomial, making the construction of generalized eigenvectors by $\varepsilon$-differentiation transparent. At EPs, where the quasi-Hamiltonian becomes defective, we derive the Jordan normal form and construct the associated generalized eigenvectors, which yields the correct counting of independent many-body eigenstates. We further show that EPs act as branch points in the complex anisotropy plane, leading to the characteristic permutation of eigenenergies and eigenstates upon encirclement. The branch-cut structure of the biorthogonal eigenstates provides direct evidence for the exchange of eigenstates when an EP is encircled. These results provide an analytically controlled many-body platform for studying EP physics and non-Hermitian topology beyond momentum-space descriptions.