Fermionic multiscale entanglement renormalization ansatz

In a recent contribution [P. Corboz, G. Evenbly, F. Verstraete, and G. Vidal, arXiv: 0904.4151 (unpublished)] entanglement renormalization was generalized to fermionic lattice systems in two spatial dimensions. Entanglement renormalization is a real-space coarse-graining transformation for lattice systems that produces a variational ansatz, the multiscale entanglement renormalization ansatz (MERA), for the ground states of local Hamiltonians. In this paper we describe in detail the fermionic version of the MERA formalism and algorithm. Starting from the bosonic MERA, which can be regarded both as a quantum circuit or in relation to a coarse-graining transformation, we indicate how the scheme needs to be modified to simulate fermions. To confirm the validity of the approach, we present benchmark results for free and interacting fermions on a square lattice with sizes between 6 x 6 and 162 x 162 and with periodic boundary conditions. The present formulation of the approach applies to generic tensor network algorithms.