In this article, we intend to address this question using the spin-Boson Hamiltonian (SBH) model, which provides a simple, low-computational cost mathematical representation of the nonadiabatic process. Indeed, studies showcased the efficiency of trajectory surface hopping, AIMS, and MCTDH on the SBH model and its variants in the ultrafast timescale. Nevertheless, we aim to push those methodologies to the edge, into the realm of tens of picoseconds. To demonstrate how hard it has been to accomplish this task, even with the simple analytical SBH model, some of the simulations we reported here have been running for over two years.