  We present calculations on energy- and time-resolved two-photon photoemission
spectra of images states in Cu(100) and Cu(111) surfaces. The surface is
modeled by a 1D effective potential and the states are propagated within a
real-space, real-time method. To obtain the energy resolved spectra we employ a
geometrical approach based on a subdivision of space into two regions. We treat
electronic inelastic effects by taking into account the scattering rates
calculated within a GW scheme. To get further insight into the decaying
mechanism we have also studied the effect of the variation of the classical
Hartree potential during the excitation. This effect turns out to be small.
