  We present an algorithm for the efficient simulation of the half-filled
spinless $t$-$V$ model on bipartite lattices, which combines the stochastic
series expansion method with determinantal quantum Monte Carlo techniques
widely used in fermionic simulations. The algorithm scales linearly in the
inverse temperature, cubically with the system size and is free from the
time-discretization error. We use it to map out the finite temperature phase
diagram of the spinless $t$-$V$ model on the honeycomb lattice and observe a
suppression of the critical temperature of the charge density wave phase in the
vicinity of a fermionic quantum critical point.
