Implementation of bond weighting method for the Grassmann tensor renormalization group

Abstract: We demonstrate the efficiency of the bond weighting method for the Grassmann tensor renormalization group (TRG). Benchmarking with the two-dimensional Gross-Neveu model with the Wilson fermion at finite density, we show that the bond weighting method improves the accuracy of the original Grassmann TRG. We also provide a sample code of the bond-weighted TRG that can be applied to the two-dimensional models including fermions on a square lattice.