Efficient calculation of time- and frequency-resolved four-wave-mixing signals

Dissipative dynamics at conical intersections: simulations with the hierarchy equations of motion method

The effect of a dissipative environment on the ultrafast nonadiabatic dynamics at conical intersections is analyzed for a two-state two-mode model chosen to represent the S₂(ππ*)–S₁(nπ*) conical intersection in pyrazine (the system) which is bilinearly coupled to infinitely many harmonic oscillators in thermal equilibrium (the bath). The system–bath coupling is modeled by the Drude spectral function. The equation of motion for the reduced density matrix of the system is solved numerically exactly with the hierarchy equation of motion method using graphics-processor-unit (GPU) technology. The simulations are valid for arbitrary strength of the system–bath coupling and arbitrary bath memory relaxation time. The present computational studies overcome the limitations of weak system–bath coupling and short memory relaxation time inherent in previous simulations based on multi-level Redfield theory [A. Kühl and W. Domcke, J. Chem. Phys. 2002, 116, 263]. Time evolutions of electronic state populations and time-dependent reduced probability densities of the coupling and tuning modes of the conical intersection have been obtained. It is found that even weak coupling to the bath effectively suppresses the irregular fluctuations of the electronic populations of the isolated two-mode conical intersection. While the population of the upper adiabatic electronic state (S₂) is very efficiently quenched by the system–bath coupling, the population of the diabatic ππ* electronic state exhibits long-lived oscillations driven by coherent motion of the tuning mode. Counterintuitively, the coupling to the bath can lead to an enhanced lifetime of the coherence of the tuning mode as a result of effective damping of the highly excited coupling mode, which reduces the strong mode–mode coupling inherent to the conical intersection. The present results extend previous studies of the dissipative dynamics at conical intersections to the nonperturbative regime of system–bath coupling. They pave the way for future first-principles simulations of femtosecond time-resolved four-wave-mixing spectra of chromophores in condensed phases which are nonperturbative in the system dynamics, the system–bath coupling as well as the field-matter coupling.

Two-dimensional optical spectroscopy of homo- and heterodimers

We theoretically study the two-dimensional (2D) spectroscopy of molecular dimers.

Simulation of femtosecond two-dimensional electronic spectra of conical intersections

We have simulated femtosecond two-dimensional (2D) electronic spectra for an excited-state conical intersection using the wave-function version of the equation-of-motion phase-matching approach. We show that 2D spectra at fixed values of the waiting time provide information on the structure of the vibronic eigenstates of the conical intersection, while the evolution of the spectra with the waiting time reveals predominantly ground-state wave-packet dynamics. The results show that 2D spectra of conical intersection systems differ significantly from those obtained for chromophores with well separated excited-state potential-energy surfaces. The spectral signatures which can be attributed to conical intersections are discussed.

Quantum diffusion wave-function approach to two-dimensional vibronic spectroscopy

We apply the quantum diffusion wavefunction approach to calculate vibronic two-dimensional (2D) spectra. As an example, we use a system consisting of two electronic states with harmonic oscillator potentials which are coupled to a bath and interact with three time-delayed laser pulses. The first- and second-order perturbative wave functions which enter into the expression for the third-order polarization are determined for a sufficient number of stochastic runs. The wave-packet approach, besides being an alternative technique to calculate the spectra, offers an intuitive insight into the dissipation dynamics and its relation to the 2D vibronic spectra.

Polaron dynamics in two-dimensional photon-echo spectroscopy of molecular rings

We have developed a new approach to the computation of third-order spectroscopic signals of molecular rings, by incorporating the Davydov soliton theory into the nonlinear response function formalism. The Davydov D1 and D Ansätze have been employed to treat the interactions between the excitons and the primary phonons, allowing for a full description of arbitrary exciton-phonon coupling strengths. As an illustration, we have simulated a series of optical 2D spectra for two models of molecular rings.

Fifth-order three-dimensional electronic spectroscopy using a pump-probe configuration

We present the theoretical details and experimental demonstration of fifth-order three-dimensional (3D) electronic spectroscopy using a pump-probe beam geometry. This is achieved using a pulse shaper and appropriate phase cycling schemes. We show how 8-step and 27-step phase cycling schemes can measure purely absorptive 3D spectra as well as 3D spectra for the individual fifth-order processes that contribute to the purely absorptive spectrum. 3D spectra as a function of two separate controllable waiting time periods can be obtained. The peak shapes and positions of the peaks in the experimental measurement correspond well to theory.

Two-dimensional vibronic spectroscopy of coherent wave-packet motion

We theoretically study two-dimensional (2D) spectroscopic signals obtained from femtosecond pulse interactions with diatomic molecules. The vibrational wave-packet dynamics is monitored in the signals. During the motion in anharmonic potentials the wave packets exhibit vibrational revivals and fractional revivals which are associated with particular quantum phases. The time-dependent phase changes are identified by inspection of the complex-valued 2D spectra. We use the Na(2) molecule as a numerical example and discuss various pulse sequences which yield information about vibrational level structure and phase relationships in different electronic states.

Interference Effects in Vibronic 2D-Spectra of Diatomic Molecules

We theoretically study interference effects in two-dimensional (2D) vibronic spectra which arise from two electronically excited states taking part in the multi-photon process initiated by femtosecond laser pulses. Therefore, a model is employed which mimiques the situation encountered in many halogen and interhalogen molecules. There, upon excitation from the ground state, an excited bound state and a dissociative state exist which are close in energy. We demonstrate that the different pathways to final states which enter into the third-order polarization result in pronounced interference patterns in the 2D-spectra.

Optical N-wave-mixing spectroscopy with strong and temporally well-separated pulses: the doorway-window representation

We have extended the doorway-window representation of optical pump-probe spectroscopy with weak pulses toward N-wave-mixing spectroscopy with temporally well-separated pulses of arbitrary strength. The expressions for the signals in the strong-pulse doorway-window representation are derived in the framework of the nonperturbative theory of N-wave-mixing spectroscopy. The strong-pulse doorway-window representation is complementary to the equation-of-motion phase-matching approach. The latter fully accounts for pulse-overlap effects in signals induced by weak pulses but is computationally more expensive. The performance of the doorway-window approximation for temporally well-separated strong pulses is illustrated for an electronic two-level system with an underdamped Condon-active vibrational mode.