Accurate calculations of the spectral density in a strongly correlated quantum many body system are of fundamental importance to study many-particle dynamics in the linear response regime. Typical examples are the calculation of inclusive and semi-exclusive scattering cross sections in atomic nuclei and transport properties of nuclear and neutron star matter.
Integral transform techniques have played an important role in accessing the spectral density in a variety of nuclear systems. However, their accuracy is in practice limited by the need to perform a numerical inversion which is often ill-conditioned.
In order to circumvent this problem, a quantum algorithm based on an appropriate expansion in Chebyshev polynomials was recently proposed. In the present work we build on this idea. We show how to perform controllable reconstructions of the spectral density over a finite energy resolution with rigours error estimates while allowing for efficient simulations on classical computers. We apply our idea to simple model response functions and comment on the applicability of the method to study realistic systems using scalable nuclear many-body methods.