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.
Zhong et al. (Thu,) studied this question.