Efficient Simulation of Inhomogeneously Correlated Systems Using Block Interaction Product States

The strength of the density matrix renormalization group (DMRG) in handling strongly correlated systems lies in its unbiased and simultaneous treatment of identical sites that are both energetically degenerate and spatially similar, as typically encountered in physical models. However, this very feature becomes a drawback when DMRG is applied to quantum chemistry calculations for large, realistic correlated systems. This is because entangled orbitals often span broad ranges in both energy and space, with their interactions being notably inhomogeneous. In this study, we suggest addressing the strong intrafragment correlations and weak interfragment correlations separately, utilizing a large-scale multiconfigurational calculation framework grounded in the block interaction product state formulation. The strong intrafragment correlation can be encapsulated in several electronic states located on fragments, which are obtained by considering the entanglement between fragments and their environments. Moreover, we incorporate non-Abelian spin-SU(2) symmetry in our work to target the desired states we interested with well-defined particle number and spin, providing deeper insights into the corresponding chemical processes. The described method has been examined in various chemical systems and demonstrates high efficiency in addressing the inhomogeneous effects in strong correlation quantum chemistry.