where V (R n ) is the state-dependent potential from Eq. 7 evaluated at nuclear configuration R n after n nuclear timesteps. Note that at t = 0, Û (0) = 1 N , which is the identity matrix in the electronic Hilbert space with N states. For spin-LSC, one could say that Û (t) = 1 N for all time t.