  The Stochastic Series Expansion method (SSE) is a Quantum Monte Carlo (QMC)
technique working directly in the imaginary time continuum and thus avoiding
"Trotter discretization" errors. Using a non-local "operator-loop update" it
allows treating large quantum mechanical systems of many thousand sites. In
this paper we first give a comprehensive review on SSE and present benchmark
calculations of SSE's scaling behavior with system size and inverse
temperature, and compare it to the loop algorithm, whose scaling is known to be
one of the best of all QMC methods. Finally we introduce a new and efficient
algorithm to measure Green's functions and thus dynamical properties within
SSE.
