Studying non-adiabatic dynamics in quantum systems, particularly in condensed phase, is a central challenge for modern theoretical chemistry. In this context, we are interested in the quantum dynamics of a system with N electronic states coupled to nuclear degrees of freedom (DOFs) as follows Ĥ = TR + U 0 ( R) Î + Ve ( R)