(1) min and ϵ (2) min are also on the order of eV. The same challenge also exists in other propagation methods, such as using RK4 to solve Eq. 17, because when a large energy difference exists in a Hamiltonian, one would need a much smaller ∆t to fully capture the phase oscillations when propagating the coefficients.