Unveiling Correlated Two-dimensional Topological Insulators through Fermionic Tensor Network States -- Classification, Edge Theories and Variational Wavefunctions

The study of topological band insulators has revealed fascinating phases characterized by band topology indices and anomalous boundary modes protected by global symmetries. In strongly correlated systems, where the traditional notion of electronic bands becomes obsolete, it has been established that topological insulator phases persist as stable phases, separate from the trivial insulators. However, due to the inability to express the ground states of such systems as Slater determinants, the formulation of generic variational wave functions for numerical simulations is highly desirable. In this paper, we tackle this challenge for two-dimensional topological insulators by developing a comprehensive framework for fermionic tensor network states. Starting from simple assumptions, we obtain possible sets of tensor equations for any given symmetry group, capturing consistent relations governing symmetry transformation rules on tensor legs. We then examine the connection between these tensor equations and