A quantum Langevin equation approach for two-dimensional electronic spectra of coupled vibrational and electronic dynamics

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

The influence of a Hamiltonian vibration vs a bath vibration on the 2D electronic spectra of a homodimer

We elucidate the influence of the system–bath boundary placement within an open quantum system, with emphasis on the two-dimensional electronic spectra, through the application of the hierarchical equations of motion formalism for an exciton system. We apply two different models, the Hamiltonian vibration model (HVM) and bath vibration model (BVM), to a monomer and a homodimer. In the HVM, we specifically include the vibronic states in the Hamiltonian capturing vibronic quenching, whereas in the BVM, all vibrational details are contained within the bath and described by an underdamped spectral density. The resultant spectra are analyzed in terms of energetic peak position and thermodynamic broadening precision in order to evaluate the efficacy of the two models. The HVM produces 2D spectra with accurate peak positional information, while the BVM is well suited to modeling dynamic peak broadening. For the monomer, both models produce equivalent spectra in the limit where additional damping associated with the underdamped vibration in the BVM approaches zero. This is supported by analytical results. However, for the homodimer, the BVM spectra are redshifted with respect to the HVM due to an absence of vibronic quenching in the BVM. The computational efficiency of the two models is also discussed in order to inform us of the most appropriate use of each method.