Worm algorithm for the O(2N) Gross-Neveu model

We study the lattice O(2N) Gross-Neveu model with Wilson fermions in the fermion loop formulation. Employing a worm algorithm for an open fermionic string, we simulate fluctuating topological boundary conditions and use them to tune the system to the critical point. We show how the worm algorithm can be extended to sample correlation functions of bound states involving an arbitrary number of Majorana fermions and present first results.