For more details about the implementation of TSH in the CPMD package, see Refs. Recently, some of us implemented rRESPA in the CPMD package for BOMD. The important details of this implementation are summarized in section II B 1, while in section II B 2, we suggest an extension of the MTS algorithm to enable non-adiabatic dynamics in the context of TSH-MD.