Hello,
the calculation of MPSH states considers only elastic scattering. In principle, if you want to use it to analyse a IETS signal, you should extract the normal modes of the vibrating system, and calculate the MPSH when the geometry is displaced along a particular normal mode.
We could have a closer look at the reason why your calculation fails if you attach your atkpython script.
Regards,
Daniele.