Calculating Green functions from finite systems

In calculating Green functions for interacting quantum systems numerically one often has to resort to finite systems which introduces a finite size level spacing. In order to describe the limit of system size going to infinity correctly, one has to introduce an artificial broadening larger than the finite size level discretization. In this work we compare various discretization schemes for impurity problems, i.e. a small system coupled to leads. Starting from a naive linear discretization we will then discuss the logarithmic discretization of the Wilson NRG, compare it to damped boundary conditions and arbitrary discretization in energy space. We then discuss the importance of choosing the right single particle basis when calculating bulk spectral functions. Finally we show the influence of damped boundary conditions on the time evolution of wave packets leading to a NRG-tsunami.