A key distinction must be drawn between how the cumulants may be obtained in the context of realistic condensed-phase systems versus model Hamiltonians. For chromophores embedded in complex environments, the quantum correlation functions required to construct the cumulant expansion are generally inaccessible. Instead, cumulants are constructed through classical MD simulations in conjunction with calculations of the vertical excitation energy along the generated trajectory. From these vertical energies, classical correlation functions can be computed, and the approximate quantum correlation functions can be reconstructed with the help of quantum correction factors (QCFs). For the one-and two-time correlation functions necessary to construct the second and third order cumulant, the QCFs used in this work take the following form in the frequency domain: C qm (ω)