Time-domain (TD) methods for the solution of Maxwell's equations are particularly appealing for their ability to provide the overall characteristics of an electrical system in a single simulation run. In many situations, such TD methods require computing the system's impulse response and using it in a convolution-based solver. In this work, we propose the evaluation of the scattering-parameters-type impulse response of partial element equivalent circuit (PEEC) models by firstly computing the scattering parameters pertaining to a unit-step excitation via the Numerical Inversion of Laplace Transform (NILT) technique, followed by recovering the corresponding impulse response. The accuracy and effectiveness of the advocated approach is validated by means of numerical experiments comparing its performance with that of more standard methods.@en