Matrix-product-state GPLHR for non-Hermitian many-body eigenproblems

Computing multiple eigenstates of non-Hermitian many-body operators plays an important role in investigating excitation spectra, spectral properties, and relaxation dynamics in interacting and open quantum systems. However, this task is challenging because the eigenvectors are generally nonorthogonal and the Hilbert space dimension grows exponentially with system size. Here, we formulate the generalized preconditioned locally harmonic residual (GPLHR) method in a matrix product state (MPS) framework to efficiently compute multiple targeted right eigenstates whose eigenvalues lie near a specified shift in the complex plane. Successive randomized compression (SRC) is employed to accelerate operator applications. We test the method on transcorrelated and complex Hubbard chains, an interacting non-Hermitian Su-Schrieffer-Heeger (SSH) chain, and a Liouvillian model. The results demonstrate that GPLHR combined with MPS provides an efficient and accurate approach for computing multiple eigenstates of non-Hermitian many-body operators.