where Â(0) and ρb (0) can be interpreted as the factorized initial electronic and nuclear configuration for the system at t = 0. In this work, we take Â and B to be operators in the N × N electronic sub-space. ρb is the nuclear density operator, and the trace is over both nuclear and electronic degrees of freedom (DOFs).