Nonadiabatic molecular dynamics (NAMD), where the motion of nuclei and electrons cannot be uncoupled and the dynamics occur beyond the Born-Oppenheimer regime, occur in many key chemical processes, such as photo-initiated molecular dissociation, photocatalytic reactions that use electronic excited states to drive chemical reactions, and photoinduced charge separation dynamics in photovoltaics. However, such dynamics are challenging to model with a fully quantum treatment of nuclei and electrons, especially for complex and condensed phase systems.