We present a theory of the four-wave-mixing (FWM) spectroscopy for the mult
ilevel system coupled to the bath with arbitrary time scale, in contrast to
most previous work that assumes the fast bath dynamic limit (motionally na
rrowed limit). Exact quantum mechanical expressions for the relevant nonlin
ear response tensor and the nonlinear polarization vector are obtained for
a model Hamiltonian. Effects of molecular rotation and the polarization dir
ections of the input radiation fields are also fully taken into account. Th
e results show how the cross correlation as well as the autocorrelation fun
ctions of the fluctuations in the energy-gaps or transition frequencies bet
ween eigenstates contribute to the FWM spectroscopy. The result is applied
to three-pulse echo spectroscopy. We examine the results for a variety of m
odel spectral density functions, and recover the results of the previously
reported theories as limiting cases. (C) 2001 American Institute of Physics
.