arXiv:2610.15387

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

Xingyao Wu, Zhangyu Yang, Yanbin Chen, Yu Wang, Feng Chen

Published Oct 11, 2026• 0 citations

Abstract

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.

DOI: 10.48550/arxiv.2610.15387

Related papers

Variational Matrix Product State Approach for Non-Hermitian System Based on a Companion Hermitian Hamiltonian

Zhen Guo, Zheng-Tao Xu, Meng Li +2 more

Non-Hermitian systems exhibiting topological properties are attracting growing interest. In this work, we propose an algorithm for solving the ground state of a non-Hermitian system in the matrix product state (MPS) formalism based on a companion Hermitian Hamiltonian. If the eigenvalues of the non-Hermitian system are known, the companion Hermitian Hamiltonian can be directly constructed and solved using Hermitian variational methods. When the eigenvalues are unknown, a gradient descent along with the companion Hermitian Hamiltonian yields both the ground state eigenenergy and the eigenstate. With the variational principle as a solid foundation, our algorithm ensures convergence and provides results in excellent agreement with the exact solutions of the non-Hermitian Su-Schrieffer-Heeger (nH-SSH) model as well as its interacting extension. The approach we present avoids solving any non-Hermitian matrix and overcomes numerical instabilities commonly encountered in large non-Hermitian systems.

arXiv:2210.14858Oct 26, 20221 citation

Density Matrix Renormalization Group Algorithm for non-Hermitian Systems

Peigeng Zhong, Wei Pan, Hai‐Qing Lin +2 more

A biorthonormal-block density-matrix renormalization group algorithm is proposed to accurately compute properties of large-scale non-Hermitian many-body systems, in which a renormalized-space partition of the non-Hermitian reduced density matrix is implemented to fulfill the prerequisite for the biorthonormality of the renormalization group (RG) transformation and to optimize the construction of saved Hilbert spaces. A redundancy in saved spaces of the reduced density matrix is exploited to reduce a condition number resulting from the nonunitarity of the left and right transformation matrices, in order to ensure the numerical stability of the RG procedure. The algorithm is successfully applied to an interacting fermionic Su-Schrieffer-Heeger model with nonreciprocal hoppings and staggered complex chemical potential, exhibiting novel many-body phenomena.

arXiv:2401.15000Physical Review LettersJul 31, 202516 citations

Density Matrix Renormalization Group for Transcorrelated Hamiltonians: Ground and Excited States in Molecules

Ke Liao, Huanchen Zhai, Evelin Martine Corvid Christlmaier +4 more

We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians and show applications in molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the Jastrow factor used for the similarity transformation. It also improves the accuracy of approximate methods such as coupled cluster singles and doubles (CCSD) as shown by recent studies. However, the non-Hermiticity of the TC Hamiltonians poses challenges for variational methods like DMRG. Imaginary-time evolution on the matrix product state (MPS) in the DMRG framework has been proposed to circumvent this problem, but this is currently limited to treating the ground state and has lower efficiency than the time-independent DMRG (TI-DMRG) due to the need to eliminate Trotter errors. In this work, we show that with minimal changes to the existing TI-DMRG algorithm, namely, replacing the original Davidson solver with the general Davidson solver to solve the non-Hermitian effective Hamiltonians at each site for a few low-lying right eigenstates, and following the rest of the original DMRG recipe, one can find the ground and excited states with improved efficiency compared to the original DMRG when extrapolating to the infinite bond dimension limit in the same basis set. An accelerated basis set convergence rate is also observed, as expected, within the TC framework.

Journal of Chemical Theory and ComputationMar 13, 202333 citations

Finding Matrix Product State Representations of Highly Excited Eigenstates of Many-Body Localized Hamiltonians

Xiongjie Yu, David Pekker, Bryan K. Clark

A key property of many-body localized Hamiltonians is the area law entanglement of even highly excited eigenstates. Matrix product states (MPS) can be used to efficiently represent low entanglement (area law) wave functions in one dimension. An important application of MPS is the widely used density matrix renormalization group (DMRG) algorithm for finding ground states of one-dimensional Hamiltonians. Here, we develop two algorithms, the shift-and-invert MPS (SIMPS) and excited state DMRG which find highly excited eigenstates of many-body localized Hamiltonians. Excited state DMRG uses a modified sweeping procedure to identify eigenstates, whereas SIMPS applies the inverse of the shifted Hamiltonian to a MPS multiple times to project out the targeted eigenstate. To demonstrate the power of these methods, we verify the breakdown of the eigenstate thermalization hypothesis in the many-body localized phase of the random field Heisenberg model, show the saturation of entanglement in the many-body localized phase, and generate local excitations.

arXiv:1509.01244Physical Review LettersJan 3, 2017121 citations

Efficient and systematic calculation of arbitrary observables for the matrix product state excitation ansatz

Jesse J. Osborne, Ian P. McCulloch

Numerical methods based on matrix product states (MPSs) are currently the standard for calculating the ground-state properties of (quasi-)one-dimensional quantum many-body systems. While the properties of the low-lying excitations in such systems are often studied in this MPS framework through by means of time-evolution simulations, we can also look at their by directly calculating eigenstates corresponding to these excitations. The so-called MPS excitation ansatz is a powerful method for finding such eigenstates with a single-particle character in the thermodynamic limit. Although this excitation ansatz has been used quite extensively, a general method for calculating expectation values for these states is lacking in the literature: We aim to fill this gap by presenting a recursive algorithm to calculate arbitrary observables expressed as matrix product operators. This method concisely encapsulates existing methods for—as well as extensions to—the excitation ansatz, such as excitations with a larger spatial support and multiparticle excitations, and is robust enough to handle further innovations. We demonstrate the versatility of our method by studying the low-lying excitations in the spin-1 Heisenberg chain and the one-dimensional Hubbard model, looking at how the excitations converge in the former, while in the latter, we present a refined method of targeting single-particle excitations inside a continuum by minimizing the energy rather than the energy itself. We hope that this technique will foster further advancements with the excitation ansatz.

arXiv:2408.17117Physical Review ResearchApr 7, 20253 citations

Simulating non-equilibrium dynamics and finite temperature physics: efficient representations for matrix product states

Lucas Kohn

Experimental advances have made it possible to realize and control quantum many-body systems, allowing the experimental study of non-equilibrium phenomena. Due to the Hilbert space growing exponentially with the number of particles, advanced numerical techniques are usually required to understand many-body effects at the theoretical level. In this thesis, we address the question of how to efficiently simulate the non-equilibrium dynamics of many-body systems using matrix product states (MPS). In particular, we consider quantum systems in interaction with macroscopic environments, as they appear in Anderson impurity problems, quantum thermodynamics, or open quantum systems. We develop a low-entanglement representation of the environment with only short-range interactions, perfectly suitable for simulations with MPS. We further show that an interleaved ordering of tight-binding chains can significantly reduce the creation of entanglement, as opposed to an intuitive implementation of the Hamiltonian's geometry. Our approach allows long-time simulations and an analysis of the environmental dynamics, in which, after a sudden quench, we find clear signatures of many-body effects. We employ our techniques to compute spectral functions of impurity models, with possible applications to dynamical mean-field theory impurity solvers, and to calculate dissipation in the non-equilibrium Anderson model with explicit time-dependence, as relevant for ongoing experiments with oscillating tip atomic force microscopes.

Sep 29, 20210 citations