Efficient calculation of imaginary time displaced correlation functions in the projector auxiliary field quantum Monte-Carlo algorithm

The calculation of imaginary time displaced correlation functions with the auxiliary field projector quantum Monte-Carlo algorithm provides valuable insight (such as spin and charge gaps) in the model under consideration. One of the authors and M. Imada [F.F. Assaad and M. Imada, J. Phys. Soc. Jpn. 65 189 (1996).] have proposed a numerically stable method to compute those quantities. Although precise this method is expensive in CPU time. Here, we present an alternative approach which is an order of magnitude quicker, just as precise, and very simple to implement. The method is based on the observation that for a given auxiliary field the equal time Green function matrix, $G$, is a projector: $G^2=G$.