where () n t  is the index of the active state for the n-th trajectory at time t. It is well known that population constructed by the occupation of active states is more reliable, but coherence cannot be directly obtained based on active states. In contrast, we can obtain both population and coherence distributions using the electronic wavefunctions along trajectories, which are denoted as ( ,