Link auxiliary field method in the Extended Hubbard model

The aim of this work is to generalize the method of Hubbard fields in fermion Quantum Monte Carlo simulation to the case of link fields. The introduced Hubbard link fields play a role of the interaction fields responsible for the attraction and repulsion of electronic excitations at the nodes. Such improvements leads to a number of advantages concerning thermalization time, autocorrelations, metastable states and probability distributions of the quantities. It is demonstrated that computations using five fields are more stable and reliable in case of energy observables.