Global fermionic mode optimization via swap gates

Abstract

An optimal choice of single-particle modes leads to an optimal tensor network-based compression of the many-body wave function for systems of indistinguishable particles, which compression can be essential for large, close-to-critical systems. The proposed mode optimization relies on the minimization of the block entropy area, a global quantity that measures entanglement of the wave function along the one-dimensional chain of modes. Extension of the two-site DMRG algorithm by nearest-neighbor mode optimizations is straightforward, but it performs poorly without systematic reordering of modes. Moreover, finding a stationary optimal point in the combination of the fixed-rank MPS manifold and the unitary group of mode rotations requires optimization for every generator, i.e., for every two-mode pair. Here, a systematic joint optimization protocol over the MPS manifold and the full unitary group is presented in which every two-mode pair is accessed by a carefully constructed sequence of two-qubit operations, including swap-gate controlled permutations. Large-scale DMRG simulations of strongly correlated two-dimensional fermionic lattice models and the multireference Fe${}_4$S${}_4$ transition metal cluster demonstrate rapid convergence towards the optimum, achieving low energy and low entanglement.

6 pages, 4 figures + Supplemental Material (7 pages, 9 figures)