, and thus the constraint in Eq. ( ) is automatically maintained. In particular, the positivity violation problem of trajectory-analyzed populations is also circumvented. And we expect this approach will improve the accuracy of coherence, as the populations based on active states in surface hopping dynamics are much more reliable compared with those by electronic wavefunctions. As the phase of coherence can only be extracted from the electronic wavefunction, we use