Various linear and nonlinear vibrational and electronic spectroscopy experiments in liquids are usually analyzed within the second-cumulant approximation, and therefore the fundamental quantity of interest is the equilibrium time-correlation function of the fluctuating transition frequency. In the usual approach the ''bath'' variables responsible for the fluctuating frequency are treated classically, leading to a classical time-correlation function. Alternatively, sometimes a quantum correction appropriate for relatively high temperatures is included, which adds an imaginary part to the classical time-correlation function. This approach, although appealing, does not satisfy detailed balance. One can consider a similar correction, but where detailed balance is satisfied, by using the harmonic quantum correction factor. In this article, we compare these approaches for a model system and two realistic examples. Our conclusion is that for linear spectroscopy the classical result is usually adequate, whereas for nonlinear spectroscopy it can be more important to include quantum corrections.dynamics ͉ liquids S pectroscopy is a powerful experimental tool for obtaining information about molecular dynamics in condensed phases, since the frequency at which a molecule absorbs light (be it an electronic or vibrational transition) is perturbed by its local environment, and as this environment changes in time the frequency of the molecule fluctuates accordingly. In the limit that these fluctuations are relatively fast (the homogeneous limit), then one can extract some information about their dynamics by simply measuring the width of the absorption line shape (1, 2). On the other hand, when these frequency fluctuations occur slowly (the inhomogeneous limit), then the line shape contains no dynamical information, and in fact is simply proportional to the distribution of frequencies (2). In this limit experiments such as the three-pulse echo are particularly useful, as they allow one to recover dynamical information absent in the line shape (3). Even within the two-level approximation (the molecule is considered to have only a ground state and a single vibrational or electronic excited state), the interpretation of experiments such as the three-pulse echo is difficult. The theoretical expression that describes the echo response is quite complicated and involves averages of exponentials of integrals of the fluctuating frequency operator (3, 4). To make progress, one typically uses a truncated cumulant expansion (3-5). Although this approximation is uncontrolled [which is sometimes accurate and sometimes not (6-9)], it allows one to express the echo response as a function only of the equilibrium frequency time-correlation function (FTCF). One of the main goals of echo experiments, then, is to extract this FTCF.In general, the FTCF should be calculated quantum mechanically, and as such it is, of course, complex (10). To simplify matters, one often assumes that the system is classical enough that the real part of the quantum FTCF c...