Higher order tensor renormalization group for relativistic fermion systems

We apply the higher order tensor renormalization group to two and three dimensional relativistic fermion systems on the lattice. In order to perform a coarse-graining of tensor networks including Grassmann variables, we introduce Grassmann higher order tensor renormalization group. We test the validity of the new algorithm by comparing its results with those of exact or previous methods.