Before the discussion of trajectory analysis, we need to define and clarify some nomenclatures, which will be used extensively below. In this study, we abbreviate the vibronic density matrix of the system, which takes into account both electronic and nuclear DoFs, as density matrix. If the system wavefunction is known, density matrix can be obtained directly. For a system composed of nuclear and electronic DoFs, the system wavefunction can be generally expressed as