Since we performed the DC-FSSH and AIMS dynamics in the adiabatic representation and ML-MCTDH and MCTDH in the diabatic representation, comparing the electronic population obtained from these two representations would be unjustified. Therefore, we have calculated the diabatic electronic populations from a swarm of DC-FSSH trajectories using the method prescribed by Landry et al. (see section 2 of SI). The average adiabatic excited-state population decay for different numbers of trajectories is shown in Figure of section 3 in the SI. The population curves are well converged upon 2000 trajectories. Interestingly, upon averaging, there is no feasible difference between the diabatic and adiabatic population curves for this specially designed SBH model. Hence, we can safely move on with the adiabatic state population.