The simulation of realistic molecules using on-the-fly dynamics techniques has been the topic of much work over the last few decades. Performing on-the-fly dynamics requires two major components: (i) accurate electronic structure and (ii) propagation of the correlated electron-nuclear dynamics. Due to the size of the total Hilbert space needed to perform exact, full-dimensional quantum simulation, the community has resorted to various approximations. The most successful and computationally accessible methods are the mixed quantumclassical (MQC) approaches such as the fewest switches surface hopping (FSSH) and the mean-field Ehrenfest approach. These schemes rely on the output of electronic structure method to evolve the electronic subsystem quantum mechanically while treating the nuclear subsystem classically. Unfortunately, the MQC approximation introduces some known drawbacks, notably the breakdown of detailed balance , the artificial creation of electronic coherence, or incorrect chemical kinetics In response to these known drawbacks, nonadiabatic dynamics approaches continue to be developed in order to systematically improve the results while retaining a similar level of computational expense.