Accurate Simulation of the Hubbard Model with Finite Fermionic Projected Entangled Pair States

We demonstrate the use of finite-size fermionic projected entangled pair states, in conjunction with variational Monte Carlo, to perform accurate simulations of the ground-state of the 2D Hubbard model. Using bond dimensions of up to D=28, we show that we can surpass state-of-the-art DMRG energies that use up to m=32000 SU(2) multiplets on 8-leg ladders. We further apply our methodology to 10x16, 12x16 and 16x16 lattices at 1/8 hole doping, and observe the dimensional crossover between stripe orientations. Our work shows the power of finite-size fermionic tensor networks to resolve the physics of the 2D Hubbard model and related problems.