In terms of future prospects, we need to build a general and efficient decoherence scheme to obtain more accurate electronic wavefunction in surface hopping dynamics, instead of relying on active states of the trajectory ensemble. Such a self-consistent method will significantly extend the applicability of surface hopping to spectroscopy and related studies. Within the independent trajectory framework, however, it is not easy to get correct wavefunctions because WP branching events can only be depicted by more than one trajectory. Recently, we have proposed a new mean field method based on auxiliary trajectories on adiabatic PESs, which gives the correct spatial distribution of population. It is natural to adopt this strategy in surface hopping as adiabatic dynamics is intrinsically embedded. These studies are currently under way.