Fast, Accurate Simulation of Polaron Dynamics and Multidimensional Spectroscopy by Multiple Davydov Trial States.