“…In the photoinduced CT processes in the solution phase, the triad is initially prepared on the ground (G) state, CPC 60 , in thermal equilibrium with the solvent, and then suddenly gets vertically photoexcited to the P-localized excitonic ππ* state, CP*C 60 . Then, there could occur nonradiative transitions to the excited P-to-C 60 CT state, CP + C 60 – , which is referred to as CT1, or to the excited C-to-C 60 charge-separated state, C + PC 60 – , which is referred to as CT2. , We previously found that the Marcus CT rate constant for ππ* → CT1 transition (picosecond) is larger than that for ππ* → CT2 transition (sub-microsecond), and the CT rate constants vary significantly with the triad conformations . Furthermore, LSC nonequilibrium Fermi’s golden rule (NE-FGR) ,, and its classical limit, that is, the instantaneous Marcus theory (IMT) , calculations, reveal that the triad has a significant nonequilibrium effect caused by the nonequilibrium initial nuclear state, which is manifested in the transient CT rate coefficient that can be enhanced to be about 40 times larger than the CT rate constant, dramatically changing the CT dynamics. , However, in previous studies, only a specific electronic transition can be simulated using LSC NE-FGR or IMT, that is, either pathway or pathway could be simulated on the all-atom level, where the population exchange with the third excited state is not allowed. , A general strategy for nonadiabatic dynamics on multiple electronic states with atomistic details is hence desirable.…”