“…[18] To date, TDDFT have been employed to deal with a lot of highly nonlinear phenomena, such as multiphoton ionization, [19] high harmonic generation (HHG) in gases and solids, [20,21] as well as nonadiabatic molecular dynamics simulations. [18] Comparing with the other NAMD methods, such as Tully surface hopping, [22] multi configuration timedependent Hartree (MCTDHF), [23] directly solving TDSE, [24] TDDFT-MD scheme has several advantages: it has favorable system-size scaling, and its efficiency allows calculations on large systems for even picoseconds; besides, TDDFT provides an efficient way to describe the strong laser field assisting electron response (ionization, excitation), which is the key part in this work. This is the reason why we choose TDDFT-MD method in this work.…”