  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.
