Theoretical simulation of nonadiabatic dynamics is essential for unraveling the intrinsic mechanisms of many important dynamical processes across physics, chemistry, biology and material science, e.g., proton transfer, charge transport, exciton relaxation and diffusion. In general, due to the structure of the Hilbert space, the computational cost of exact quantum dynamics scales exponentially as the number of degrees of freedom (DoFs) increases. In comparison, the computational cost in classical dynamics could scale linearly with the number of DoFs. As heavier particles often behave more classically, treating nuclear DoFs as classical trajectories is a more efficient strategy in nonadiabatic dynamics simulations of realistic systems.