The scale spanned by (bio-)chemical processes ranges from the decomposition of small organic molecules due to attosecond laser pulses to the folding of proteins comprising hundreds of thousands of atoms over several micro-to milliseconds. Exploring the dynamics of those processes requires methods of varying precision and computational cost, where usually a trade-off between the former and the latter have to be found in order to describe chemically relevant time scales and system sizes. In particular, while the system sizes that can be afforded with accurate first-principles methods are constantly growing, the requirement of sampling a sufficiently long time scale often limits the system size to a few hundreds of atoms in Kohn-Sham density functional theory (KS-DFT) based molecular dynamics (MD) and to less than a few dozen for accurate wavefunction-based methods.