Dynamic correlation functions in quantum systems: A Monte Carlo algorithm

We propose a new algorithm to compute dynamic correlation functions in quantum systems at zero temperature using Monte Carlo techniques. The method, which avoids inverse-Laplace transforms, is based on a representation of the Dirac $\delta$ function as the limit of a Gaussian. We illustrate the method on the one-dimensional harmonic oscillator.