Along the QM/MM trajectory, vertical excitation energies were computed in 2 fs intervals using time-dependent density-functional theory (TDDFT) as implemented in the TeraChem code, with the same basis set and functional as for the ground state dynamics. A total of 10,000 individual vertical excitation energies were then used to compute classical correlation functions C