Continuous Time Quantum Monte Carlo method for fermions

We present numerically exact continuous-time Quantum Monte Carlo algorithm for fermions with a general non-local in space-time interaction. The new determinantal grand-canonical scheme is based on a stochastic series expansion for the partition function in the interaction representation. The method is particularly applicable for multi-band time-dependent correlations since it does not invoke the Hubbard-Stratonovich transformation. The test calculations for exactly solvable models as well results for the Green function and for the time-dependent susceptibility of the multi-band super-symmetric model with a spin-flip interaction are discussed.