Quantum Path-integral Method for Fictitious Particle Hubbard Model

Abstract

We formulate a path-integral Monte Carlo algorithm for simulating lattice systems consisting of fictitious particles governed by a generalized exchange statistics. This method, initially proposed for continuum systems, introduces a continuous parameter $\xi$ in the partition function that interpolates between bosonic ($\xi = 1$) and fermionic ($\xi = -1$) statistics. We generalize this approach to discrete lattice models and apply it to the two-dimensional Hubbard model of fictitious particles, including the Bose- and Fermi-Hubbard models as special cases. By combining reweighting and $\xi$-extrapolation techniques, we access both half-filled and doped regimes. In particular, we demonstrate that the method remains effective even in strongly correlated, doped systems where the fermion sign problem hinders conventional quantum Monte Carlo approaches. Our results validate the applicability of the fictitious particle framework on lattice models and establish it as a promising tool for sign-problem mitigation in strongly interacting fermionic systems.