show the time evolution of 22 () t  and 2 12 () t  in the low-energy case with k0 = 16 au. It is apparent that BCSH-as gives the correct time evolution of population while the results of both FSSH-as and CT-MQC deviate significantly from the quantum reference for t > 1000 au. In comparison, the asymptotic population by CT-MQC is slightly better than that of FSSH-as. Concerning the indicator of coherence, despite being a coupled trajectory method, CT-MQC still cannot give accurate predictions. By employing our consistent construction method for the density matrix, even the trajectories generated by FSSH can roughly result in the correct asymptotic behavior of the indicator of coherence. After implementing phase and branching corrections, BCSH-cc achieves much better time-dependent behavior with only a small overestimation of the coherence.