“…Trajectory-based semi-classical non-adiabatic dynamics simulations approaches with different schemes were extensively employed to simulate a series of ultrafast excited-state relaxation processes in chemical and biological systems, and materials [ 87 , 95 , 96 , 112 , 113 , 114 , 115 , 116 , 117 , 118 , 119 , 120 , 121 , 122 , 123 , 124 , 125 , 126 ]. In the present work, trajectory-based surface hopping dynamics simulations involving the lowest two singlet states were carried out with the recently developed method proposed by Zhu et al [ 127 ].…”