The NA-AIMD simulations are performed using Gaussian 16 interfaced with Newton-X. The 1000 structures for the ensemble absorption spectrum are used as initial conditions centered at the high energy absorption peaks as indicated by the blue lines in Fig. . The energy width of each excitation was chosen such that about 100 initial conditions were accepted.