A general formalism for introducing nuclear quantum effects in the expression of the quantum time correlation function of an operator in a multi-level electronic system is presented in the adiabatic limit. The final formula includes the nuclear quantum time correlation functions of the operator matrix elements, of the energy gap, and their cross terms.These quantities can be inferred and evaluated from their classical analogs obtained by mixed quantum-classical molecular dynamics simulations. The formalism is applied to the absorption spectrum of a hydrated electron, expressed in terms of the time correlation function of the dipole operator in the ground electronic state. We find that both static and dynamic nuclear quantum effects distinctly influence the shape of the absorption spectrum, especially its high-energy tail related to transitions to delocalized electron states. Their inclusion does improve significantly the agreement between theory and experiment for both the low and high frequency edges of the spectrum. It does not appear sufficient, however, to d