Source author record

Philippe Corboz

Philippe Corboz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2026arXiv

Itinerant magnetism in the triangular lattice Hubbard model at half-doping: Application to twisted transition-metal dichalcogenides

We use unrestricted Hartree-Fock, density matrix renormalization group, and variational projected entangled pair state calculations to investigate the ground state phase diagram of the triangular lattice Hubbard model at "half doping" relative to single occupancy, i.e. at a filling of $(1\pm \frac{1}{2})$ electrons per site. The electron-doped case has a nested Fermi surface in the non-interacting limit, and hence a weak-coupling instability towards density-wave orders whose wavevectors are determined by Fermi surface nesting conditions. We find that at moderate to strong interaction strengths other spatially-modulated orders arise, with wavevectors distinct from the nesting vectors. In particular, we identify a series closely-competing itinerant long-wavelength magnetically ordered states, yielding to uniform ferromagnetic order at the largest interaction strengths. For half-hole doping and a similar range of interaction strengths, our data indicate that magnetic orders are most likely absent.

preprint2022arXiv

Variational methods for contracting projected entangled-pair states

The norms or expectation values of infinite projected entangled-pair states (PEPS) cannot be computed exactly, and approximation algorithms have to be applied. In the last years, many efficient algorithms have been devised -- the corner transfer matrix renormalization group (CTMRG) and variational uniform matrix product state (VUMPS) algorithm are the most common -- but it remains unclear whether they always lead to the same results. In this paper, we identify a subclass of PEPS for which we can reformulate the contraction as a variational problem that is algorithm independent. We use this variational feature to assess and compare the accuracy of CTMRG and VUMPS contractions. Moreover, we devise a new variational contraction scheme, which we can extend to compute general N-point correlation functions.

preprint2021arXiv

Automatic differentiation applied to excitations with Projected Entangled Pair States

The excitation ansatz for tensor networks is a powerful tool for simulating the low-lying quasiparticle excitations above ground states of strongly correlated quantum many-body systems. Recently, the two-dimensional tensor network class of infinite entangled pair states gained new ground state optimization methods based on automatic differentiation, which are at the same time highly accurate and simple to implement. Naturally, the question arises whether these new ideas can also be used to optimize the excitation ansatz, which has recently been implemented in two dimensions as well. In this paper, we describe a straightforward way to reimplement the framework for excitations using automatic differentiation, and demonstrate its performance for the Hubbard model at half filling.

preprint2021arXiv

Phase diagram of the Shastry-Sutherland Compound SrCu2(BO3)2 under extreme combined conditions of field and pressure

Motivated by the intriguing properties of the Shastry-Sutherland compound SrCu2(BO3)2 under pressure, with a still debated intermediate plaquette phase appearing at around 20 kbar and a possible deconfined critical point at higher pressure upon entering the antiferromagnetic phase, we have investigated its high-field properties in this pressure range using tunnel diode oscillator (TDO) measurements. The two main new phases revealed by these measurements are fully consistent with those identified by infinite Projected Entangled Pair states (iPEPS) calculations of the Shastry-Sutherland model, a 1/5 plateau and a 10 x 2 supersolid. Remarkably, these phases are descendants of the full-plaquette phase, the prominent candidate for the intermediate phase of SrCu2(BO3)2. The emerging picture for SrCu2(BO3)2 is shown to be that of a system dominated by a tendency to an orthorhombic distortion at intermediate pressure, an important constraint on any realistic description of the transition into the antiferromagnetic phase.

preprint2021arXiv

Quantum Monte Carlo simulations of highly frustrated magnets in a cluster basis: The two-dimensional Shastry-Sutherland model

Quantum Monte Carlo (QMC) simulations constitute nowadays one of the most powerful methods to study strongly correlated quantum systems, provided that no "sign problem" arises. However, many systems of interest, including highly frustrated magnets, suffer from an average sign that is close to zero in standard QMC simulations. Nevertheless, a possible sign problem depends on the simulation basis, and here we demonstrate how a suitable choice of cluster basis can be used to eliminate or at least reduce the sign problem in highly frustrated magnets that were so far inaccessible to efficient QMC simulations. We focus in particular on the application of a two-spin (dimer)-based QMC method to the thermodynamics of the spin-1/2 Shastry-Sutherland model for SrCu$_2$(BO$_3$)$_2$.

preprint2021arXiv

Tensor network study of the $m=1/2$ magnetization plateau in the Shastry-Sutherland model at finite temperature

The two-dimensional infinite projected entangled pair states tensor network is evolved in imaginary time with the full update (FU) algorithm to simulate the Shastry-Sutherland model in a magnetic field at finite temperature directly in the thermodynamic limit. We focus on the phase transition into the $m=1/2$ magnetization plateau, which was observed in experiments on SrCu$_2$(BO$_3$)$_2$. For the largest simulated bond dimension, the early evolution in the high-temperature regime is simulated with the simple update (SU) scheme and then, as the correlation length increases, continued with the FU scheme towards the critical regime. We apply a small symmetry-breaking bias field and then extrapolate towards zero bias using a simple scaling theory in the bias field. The combined SU + FU scheme provides an accurate estimate of the critical temperature, even though the results could not be fully converged in the bond dimension in the vicinity of the transition. The critical temperature estimate is improved with a generalized scaling theory that combines two divergent length scales: One due to the bias, and the other due to the finite bond dimension. The obtained results are consistent with the transition being in the universality class of the two-dimensional classical Ising model. The estimated critical temperature is $3.5(2)$ K, which is well above the temperature $2.1$ K used in the experiments.

preprint2021arXiv

The Hubbard model: A computational perspective

The Hubbard model is the simplest model of interacting fermions on a lattice and is of similar importance to correlated electron physics as the Ising model is to statistical mechanics or the fruit fly to biomedical science. Despite its simplicity, the model exhibits an incredible wealth of phases, phase transitions, and exotic correlation phenomena. While analytical methods have provided a qualitative description of the model in certain limits, numerical tools have shown impressive progress in achieving quantitative accurate results over the last years. This article gives an introduction to the model, motivates common questions, and illustrates the progress that has been achieved over the last years in revealing various aspects of the correlation physics of the model.

preprint2020arXiv

Excitations with projected entangled pair states using the corner transfer matrix method

We present an extension of a framework for simulating single quasiparticle or collective excitations on top of strongly correlated quantum many-body ground states using infinite projected entangled pair states, a tensor network ansatz for two-dimensional wave functions in the thermodynamic limit. Our approach performs a systematic summation of locally perturbed states in order to obtain excited eigenstates localized in momentum space, using the corner transfer matrix method, and generalizes the framework to arbitrary unit cell sizes, the implementation of global Abelian symmetries and fermionic systems. Results for several test cases are presented, including the transverse Ising model, the spin-$\frac{1}{2}$ Heisenberg model and a free fermionic model, to demonstrate the capability of the method to accurately capture dispersions. We also provide insight into the nature of excitations at the $k=(π,0)$ point of the Heisenberg model.

preprint2016arXiv

Gradient methods for variational optimization of projected entangled-pair states

We present a conjugate-gradient method for the ground-state optimization of projected entangled-pair states (PEPS) in the thermodynamic limit, as a direct implementation of the variational principle within the PEPS manifold. Our optimization is based on an efficient and accurate evaluation of the gradient of the global energy functional by using effective corner environments, and is robust with respect to the initial starting points. It has the additional advantage that physical and virtual symmetries can be straightforwardly implemented. We provide the tools to compute static structure factors directly in momentum space, as well as the variance of the Hamiltonian. We benchmark our method on Ising and Heisenberg models, and show a significant improvement on the energies and order parameters as compared to algorithms based on imaginary-time evolution.

preprint2016arXiv

Improved energy extrapolation with infinite projected entangled-pair states applied to the 2D Hubbard model

An infinite projected entangled-pair state (iPEPS) is a variational tensor network ansatz for 2D wave functions in the thermodynamic limit where the accuracy can be systematically controlled by the bond dimension $D$. We show that for the doped Hubbard model in the strongly correlated regime ($U/t=8$, $n=0.875$) iPEPS yields lower variational energies than state-of-the-art variational methods in the large 2D limit, which demonstrates the competitiveness of the method. In order to obtain an accurate estimate of the energy in the exact infinite $D$ limit we introduce and test an extrapolation technique based on a truncation error computed in the iPEPS imaginary time evolution algorithm. The extrapolated energies are compared with accurate quantum Monte Carlo results at half filling and with various other methods in the doped, strongly correlated regime.

preprint2016arXiv

SU(6) Heisenberg model on the honeycomb lattice: competition between plaquette and chiral order

We revisit the SU(6) Heisenberg model on the honeycomb lattice, which has been predicted to be a chiral spin liquid by mean-field theory [G. Szirmai et al., Phys. Rev. A 84, 011611 (2011)]. Using exact diagonalizations of finite clusters, infinite projected entangled pair states simulations, and variational Monte Carlo simulations based on Gutzwiller projected wave functions, we provide strong evidence in favour of the competing plaquette state, which was reported to be higher but close by in energy according to mean-field theory. This is further confirmed by the investigation of the model with a ring exchange term, which shows that there is a transition between the plaquette state and the chiral state at a finite value of the ring exchange term.

preprint2015arXiv

Efficient Continuous-time Quantum Monte Carlo Method for the Ground State of Correlated Fermions

We present the ground state extension of the efficient quantum Monte Carlo algorithm for lattice fermions of arXiv:1411.0683. Based on continuous-time expansion of imaginary-time projection operator, the algorithm is free of systematic error and scales \emph{linearly} with projection time and interaction strength. Compared to the conventional quantum Monte Carlo methods for lattice fermions, this approach has greater flexibility and is easier to combine with powerful machinery such as histogram reweighting and extended ensemble simulation techniques. We discuss the implementation of the continuous-time projection in detail using the spinless $t-V$ model as an example and compare the numerical results with exact diagonalization, density-matrix-renormalization-group and infinite projected entangled-pair states calculations. Finally we use the method to study the fermionic quantum critical point of spinless fermions on a honeycomb lattice and confirm previous results concerning its critical exponents.

preprint2015arXiv

The iPEPS algorithm, improved: fast full update and gauge fixing

The infinite Projected Entangled Pair States (iPEPS) algorithm [J. Jordan et al, PRL 101, 250602 (2008)] has become a useful tool in the calculation of ground state properties of 2d quantum lattice systems in the thermodynamic limit. Despite its many successful implementations, the method has some limitations in its present formulation which hinder its application to some highly-entangled systems. The purpose of this paper is to unravel some of these issues, in turn enhancing the stability and efficiency of iPEPS methods. For this, we first introduce the fast full update scheme, where effective environment and iPEPS tensors are both simultaneously updated (or evolved) throughout time. As we shall show, this implies two crucial advantages: (i) dramatic computational savings, and (ii) improved overall stability. Besides, we extend the application of the \emph{local gauge fixing}, successfully implemented for finite-size PEPS [M. Lubasch, J. Ignacio Cirac, M.-C. Bañuls, PRB 90, 064425 (2014)], to the iPEPS algorithm. We see that the gauge fixing not only further improves the stability of the method, but also accelerates the convergence of the alternating least squares sweeping in the (either "full" or "fast full") tensor update scheme. The improvement in terms of computational cost and stability of the resulting "improved" iPEPS algorithm is benchmarked by studying the ground state properties of the quantum Heisenberg and transverse-field Ising models on an infinite square lattice.

preprint2014arXiv

Competing states in the t-J model: uniform d-wave state versus stripe state

Variational studies of the t-J model on the square lattice based on infinite projected-entangled pair states (iPEPS) confirm an extremely close competition between a uniform d-wave superconducting state and different stripe states. The site-centered stripe with an in-phase d-wave order has an equal or only slightly lower energy than the stripe with anti-phase d-wave order. The optimal stripe filling is not constant but increases with J/t. A nematic anisotropy reduces the pairing amplitude and the energies of stripe phases are lowered relative to the uniform state with increasing nematicity.

preprint2014arXiv

Crystals of bound states in the magnetization plateaus of the Shastry-Sutherland model

Using infinite projected entangled-pair states (iPEPS), we show that the Shastry-Sutherland model in an external magnetic field has low-magnetization plateaus which, in contrast to previous predictions, correspond to crystals of bound states of triplets, and {\it not} to crystals of triplets. The first sizable plateaus appear at magnetization 1/8, 2/15 and 1/6, in agreement with experiments on the orthogonal-dimer antiferromagnet SrCu$_2$(BO$_3$)$_2$, and they can be naturally understood as regular patterns of bound states, including the intriguing 2/15 one. We also show that, even in a confined geometry, two triplets bind into a {\it localized} bound state with $S_z=2$. Finally, we discuss the role of competing domain-wall and supersolid phases as well as that of additional anisotropic interactions.

preprint2014arXiv

Fermionic Quantum Critical Point of Spinless Fermions on a Honeycomb Lattice

Spinless fermions on a honeycomb lattice provide a minimal realization of lattice Dirac fermions. Repulsive interactions between nearest neighbors drive a quantum phase transition from a Dirac semimetal to a charge-density-wave state through a fermionic quantum critical point, where the coupling of Ising order parameter to the Dirac fermions at low energy drastically affects the quantum critical behavior. Encouraged by a recently discovery of absence of the fermion sign problem in this model, we study the fermionic quantum critical point using the continuous time quantum Monte Carlo method with worm sampling technique. We estimate the transition point $V/t= 1.356(1)$ with the critical exponents $ν=0.80(3)$ and $η=0.302(7)$. Compatible results for the transition point are also obtained with infinite projected entangled-pair states.

preprint2014arXiv

Probing the stability of the spin liquid phases in the Kitaev-Heisenberg model using tensor network algorithms

We study the extent of the spin liquid phases in the Kitaev-Heisenberg model using infinite Projected Entangled-Pair States tensor network ansatz wave functions directly in the thermodynamic limit. To assess the accuracy of the ansatz wave functions we perform benchmarks against exact results for the Kitaev model and find very good agreement for various observables. In the case of the Kitaev-Heisenberg model we confirm the existence of 6 different phases: Néel, stripy, ferromagnetic, zigzag and two spin liquid phases. We find finite extents for both spin liquid phases and discontinuous phase transitions connecting them to symmetry-broken phases.

preprint2014arXiv

RVB superconductors with fermionic projected entangled pair states

We construct a family of simple fermionic projected entangled pair states (fPEPS) on the square lattice with bond dimension $D=3$ which are exactly hole-doped resonating valence bond (RVB) wavefunctions with short-range singlet bonds. Under doping the insulating RVB spin liquid evolves immediately into a superconductor with mixed $d+is$ pairing symmetry whose pair amplitude grows as the square-root of the doping. The relative weight between $s$-wave and $d$-wave components can be controlled by a single variational parameter $c$. We optimize our ansatz w.r.t. $c$ for the frustrated $t-J_1-J_2$ model (including both nearest and next-nearest neighbor antiferromagnetic interactions $J_1$ and $J_2$, respectively) for $J_2\simeq J_1/2$ and obtain an energy very close to the infinite-PEPS state (using full update optimization and same bond dimension). The orbital symmetry of the optimized RVB superconductor has predominant d-wave character, although we argue a residual (complex s-wave) time reversal symmetry breaking component should always be present. Connections of the results to the physics of superconducting cuprates and pnictides are outlined.

preprint2013arXiv

Competing states in the SU(3) Heisenberg model on the honeycomb lattice: Plaquette valence-bond crystal versus dimerized color-ordered state

Conflicting predictions have been made for the ground state of the SU(3) Heisenberg model on the honeycomb lattice: Tensor network simulations found a plaquette order [Zhao et al, Phys. Rev. B 85, 134416 (2012)], where singlets are formed on hexagons, while linear flavor-wave theory (LFWT) suggested a dimerized, color ordered state [Lee and Yang, Phys. Rev. B 85, 100402 (2012)]. In this work we show that the former state is the true ground state by a systematic study with infinite projected-entangled pair states (iPEPS), for which the accuracy can be systematically controlled by the so-called bond dimension $D$. Both competing states can be reproduced with iPEPS by using different unit cell sizes. For small $D$ the dimer state has a lower variational energy than the plaquette state, however, for large $D$ it is the latter which becomes energetically favorable. The plaquette formation is also confirmed by exact diagonalizations and variational Monte Carlo studies, according to which both the dimerized and plaquette states are non-chiral flux states.

preprint2013arXiv

Competition between three-sublattice order and superfluidity in the quantum 3-state Potts model of ultracold bosons and fermions on a square optical lattice

We study a quantum version of the three-state Potts model that includes as special cases the effective models of bosons and fermions on the square lattice in the Mott insulating limit. It can be viewed as a model of quantum permutations with amplitudes J_parallel and J_perp for identical and different colors, respectively. For J_parallel=J_perp>0, it is equivalent to the SU(3) Heisenberg model, which describes the Mott insulating phase of 3-color fermions, while the parameter range J_perp<min(0,-J_parallel) can be realized in the Mott insulating phase of 3-color bosonic atoms. Using linear flavor wave theory, infinite projected entangled-pair states (iPEPS), and continuous-time quantum Monte-Carlo simulations, we construct the full T=0 phase diagram, and we explore the T>0 properties for J_perp<0. For dominant antiferromagnetic J_parallel interactions, a three-sublattice long-range ordered stripe state is selected out of the ground state manifold of the antiferromagnetic Potts model by quantum fluctuations. Upon increasing |J_perp|, this state is replaced by a uniform superfluid for J_perp<0, and by an exotic three-sublattice superfluid followed by a two-sublattice superfluid for J_perp>0. The transition out of the uniform superfluid (that can be realized with bosons) is shown to be a peculiar type of Kosterlitz-Thouless transition with three types of elementary vortices.

preprint2013arXiv

Phase diagram of $^{4}$He on graphene

The low temperature phase diagram of $^4$He adsorbed on a single graphene sheet is studied by computer simulation of a system comprising nearly thousand helium atoms. In the first layer, two commensurate solid phases are observed, with fillings 1/3 and 7/16 respectively, separated by a domain wall phase, as well as an incommensurate crystal at higher coverage. No evidence of a thermodynamically stable superfliuid phase is found for the first adlayer. Second layer promotion occurs at a coverage of 0.111(4) $Å^{-2}$. In the second layer two phases are observed, namely a superfluid and an incommensurate solid, with no commensurate solid intervening between these two phases. The computed phase diagram closely resembles that predicted for helium on graphite.

preprint2012arXiv

Comment on "Topological quantum phase transitions of attractive spinless fermions in a honeycomb lattice" by Poletti D. et al

In a recent letter [D. Poletti et al., EPL 93, 37008 (2011)] a model of attractive spinless fermions on the honeycomb lattice at half filling has been studied by mean-field theory, where distinct homogenous phases at rather large attraction strength $V>3.36$, separated by (topological) phase transitions, have been predicted. In this comment we argue that without additional interactions the ground states in these phases are not stable against phase separation. We determine the onset of phase separation at half filling $V_{ps}\approx 1.7$ by means of infinite projected entangled-pair states (iPEPS) and exact diagonalization.

preprint2012arXiv

Simplex solids in SU(N) Heisenberg models on the kagome and checkerboard lattices

We present a numerical study of the SU(N) Heisenberg model with the fundamental representation at each site for the kagome lattice (for N=3) and the checkerboard lattice (for N=4), which are the line graphs of the honeycomb and square lattices and thus belong to the class of bisimplex lattices. Using infinite projected entangled-pair states (iPEPS) and exact diagonalizations, we show that in both cases the ground state is a simplex solid state with a two-fold ground state degeneracy, in which the N spins belonging to a simplex (i.e. a complete graph) form a singlet. Theses states can be seen as generalizations of valence bond solid states known to be stabilized in certain SU(2) spin models.

preprint2012arXiv

Three-sublattice order in the SU(3) Heisenberg model on the square and triangular lattice

We present a numerical study of the SU(3) Heisenberg model of three-flavor fermions on the triangular and square lattice by means of the density-matrix renormalization group (DMRG) and infinite projected entangled-pair states (iPEPS). For the triangular lattice we confirm that the ground state has a three-sublattice order with a finite ordered moment which is compatible with the result from linear flavor wave theory (LFWT). The same type of order has recently been predicted also for the square lattice [PRL 105, 265301 (2010)] from LFWT and exact diagonalization. However, for this case the ordered moment cannot be computed based on LFWT due to divergent fluctuations. Our numerical study clearly supports this three-sublattice order, with an ordered moment of m=0.2-0.4 in the thermodynamic limit.

preprint2011arXiv

Simulation of fermionic lattice models in two dimensions with Projected Entangled-Pair States: Next-nearest neighbor Hamiltonians

In a recent contribution [Phys. Rev. B 81, 165104 (2010)] fermionic Projected Entangled-Pair States (PEPS) were used to approximate the ground state of free and interacting spinless fermion models, as well as the $t$-$J$ model. This paper revisits these three models in the presence of an additional next-nearest hopping amplitude in the Hamiltonian. First we explain how to account for next-nearest neighbor Hamiltonian terms in the context of fermionic PEPS algorithms based on simulating time evolution. Then we present benchmark calculations for the three models of fermions, and compare our results against analytical, mean-field, and variational Monte Carlo results, respectively. Consistent with previous computations restricted to nearest-neighbor Hamiltonians, we systematically obtain more accurate (or better converged) results for gapped phases than for gapless ones.

preprint2011arXiv

Simultaneous dimerization and SU(4) symmetry breaking of 4-color fermions on the square lattice

Using infinite projected entangled pair states (iPEPS), exact diagonalization, and flavor-wave theory, we show that the SU(4) Heisenberg model undergoes a spontaneous dimerization on the square lattice, in contrast to its SU(2) and SU(3) counterparts, which develop Neel and three-sublattice stripe-like long-range order. Since the ground state of a dimer is not a singlet for SU(4) but a 6-dimensional irrep, this leaves the door open for further symmetry breaking. We provide evidence that, unlike in SU(4) ladders, where dimers pair up to form singlet plaquettes, here the SU(4) symmetry is additionally broken, leading to a gapless spectrum in spite of the broken translational symmetry.

preprint2010arXiv

Simulation of anyons with tensor network algorithms

Interacting systems of anyons pose a unique challenge to condensed matter simulations due to their non-trivial exchange statistics. These systems are of great interest as they have the potential for robust universal quantum computation, but numerical tools for studying them are as yet limited. We show how existing tensor network algorithms may be adapted for use with systems of anyons, and demonstrate this process for the 1-D Multi-scale Entanglement Renormalisation Ansatz (MERA). We apply the MERA to infinite chains of interacting Fibonacci anyons, computing their scaling dimensions and local scaling operators. The scaling dimensions obtained are seen to be in agreement with conformal field theory. The techniques developed are applicable to any tensor network algorithm, and the ability to adapt these ansaetze for use on anyonic systems opens the door for numerical simulation of large systems of free and interacting anyons in one and two dimensions.

preprint2010arXiv

Simulation of interacting fermions with entanglement renormalization

We propose and test an algorithm to simulate a lattice system of interacting fermions in two spatial dimensions. The approach is an extension of the entanglement renormalization technique [Phys. Rev. Lett. 99, 220405 (2007)] and the related multi-scale entanglement renormalization ansatz. Benchmark calculations for free and interacting fermions on a lattice of $6\times 6$ sites and periodic boundary conditions confirm the validity of this proposal.

preprint2010arXiv

Simulation of strongly correlated fermions in two spatial dimensions with fermionic Projected Entangled-Pair States

We explain how to implement, in the context of projected entangled-pair states (PEPS), the general procedure of fermionization of a tensor network introduced in [P. Corboz, G. Vidal, Phys. Rev. B 80, 165129 (2009)]. The resulting fermionic PEPS, similar to previous proposals, can be used to study the ground state of interacting fermions on a two-dimensional lattice. As in the bosonic case, the cost of simulations depends on the amount of entanglement in the ground state and not directly on the strength of interactions. The present formulation of fermionic PEPS leads to a straightforward numerical implementation that allowed us to recycle much of the code for bosonic PEPS. We demonstrate that fermionic PEPS are a useful variational ansatz for interacting fermion systems by computing approximations to the ground state of several models on an infinite lattice. For a model of interacting spinless fermions, ground state energies lower than Hartree-Fock results are obtained, shifting the boundary between the metal and charge-density wave phases. For the t-J model, energies comparable with those of a specialized Gutzwiller-projected ansatz are also obtained.

preprint2009arXiv

Fermionic multi-scale entanglement renormalization ansatz

In a recent contribution [arXiv:0904:4151] 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 multi-scale 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 \times 6$ and $162\times 162$ and with periodic boundary conditions. The present formulation of the approach applies to generic tensor network algorithms.