  We generalize isometric tensor network states to fermionic systems, paving
the way for efficient adaptations of 1D tensor network algorithms to 2D
fermionic systems. As the first application of this formalism, we developed and
benchmarked a time-evolution block-decimation (TEBD) algorithm for real-time
and imaginary-time evolution. The imaginary-time evolution produces
ground-state energies for gapped systems, systems with a Dirac point, and
systems with gapless edge modes to good accuracy. The real-time TEBD captures
the scattering of two fermions and the chiral edge dynamics on the boundary of
a Chern insulator.
