The quantum Fokker-Planck equation is derived for a system nonlinearly coupled to a harmonic oscillator bath. The system-bath interaction is assumed to be linear in the bath coordinates but quadratic in the system coordinate. The relaxation induced dynamics of a harmonic system are investigated by simulating the higher-order correlation functions of the Raman polarizability and the dipole moment, which represent the nonlinear optical responses of Raman or infrared spectroscopy. The 5th-order Raman response shows that, in addition to the frequency fluctuations induced by the bath, higher-order energy transfer between the system and bath plays a role. The nonlinearity of the system-bath interaction yields also an interesting feature in the 7th-order Raman echo or the 3rd-order infrared photon echo response: The calculations predict a finite signal for the case of a harmonic potential and a linear coordinate dependence of the polarizability or dipole while for linear system-bath coupling this response vanishes completely due to destructive interference of different Liouville space pathways.KEYWORDS: 2D Raman, 2D IR, quantum Fokker-Planck, nonlinear system-bath interaction monly assumes a certain functional form for the Hamiltonian of the medium, which allows for the derivation of closed form expressions that are fit to the experimental results. Such model approach, which is useful to understand the complex molecular interactions involved in optical processes, is complementary to more "realistic" approaches such as Molecular dynamics simulations. A typical example for this strategy is the Brownian oscillator model, which can be used to calculate, e.g., the first-order resonant IR response or the third-order offresonant Raman response. 5,8,9) Since the dipole moment or the polarizability is usually expanded in the power of molecular coordinate Q as μ(Q) = μ 1 Q + μ 2 Q 2 · · · or α(Q) = α 1 Q + α 2 Q 2 · · ·, the first-order IR or the third-order Raman onse in lowest order of Q is nothing but the two-time correlation functions of Brownian resp 3115 * E-mail: tanimura@ims.ac.jp rules for the transitions involved.To calculate the (non)linear optical response one com- §1. IntroductionFemtosecond nonlinear optical spectroscopies have proven to be valuable and versatile tools to obtain information on the dynamic characteristics of condensed phase systems.1-6) Usually, these experiments are described by a response function formalism, 5, 7) which is based on a perturbative expansion of the optical polarization in powers of the applied electric fields. For resonant spectroscopy, which is carried by the Infrared (IR) laser for molecular vibrational modes, the laser interaction is described by μE(t), where μ is the dipole moment of the system. Then, for example, the lowestorder signal is expressed by the dipole correlation function μ(t )μ(t) . For off-resonant Raman spectroscopy, in which resonance arises from a pair of laser pulses through Raman excitation processes, the laser interaction is described by αE 2 (t), where α is ...