Doorway–window description of sequential three-pulse photon echo signals

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.

Ab Initio Quasiclassical Simulation of Femtosecond Time-Resolved Two-Dimensional Electronic Spectra of Pyrazine

Two-dimensional (2D) electronic spectroscopy is a powerful nonlinear technique which provides spectroscopic information on two frequency axes as well as dynamical information as a function of the so-called waiting time. Herein, an ab initio theoretical framework for the simulation of electronic 2D spectra has been developed. The method is based on the classical approximation to the doorway-window representation of three-pulse photon-echo signals and the description of nuclear motion by classical trajectories. Nonadiabatic effects are taken into account by a trajectory surface-hopping algorithm. 2D electronic spectra were simulated with ab initio on-the-fly trajectory calculations using the ADC(2) electronic-structure method for the pyrazine molecule, which is a benchmark system for ultrafast radiationless decay through conical intersections. It is demonstrated that 2D spectroscopy with subfemtosecond UV pulses can provide unprecedented detailed information on the ultrafast photodynamics of polyatomic molecules.

Quantum dynamics and spectroscopy of dihalogens in solid matrices. II. Theoretical aspects and G-MCTDH simulations of time-resolved coherent Raman spectra of Schrödinger cat states of the embedded I2Kr18 cluster

This study presents quantum dynamical simulations, using the Gaussian-based multiconfigurational time-dependent Hartree (G-MCTDH) method, of time-resolved coherent Raman four-wave-mixing spectroscopic experiments for the iodine molecule embedded in a cryogenic crystal krypton matrix [D. Picconi et al., J. Chem. Phys. 150, 064111 (2019)]. These experiments monitor the time-evolving vibrational coherence between two wave packets created in a quantum superposition (i.e., a "Schrödinger cat state") by a pair of pump pulses which induce electronic B Πu30⁺⟵XΣg+1 transitions. A theoretical description of the spectroscopic measurement is developed, which elucidates the connection between the nonlinear signals and the wave packet coherence. The analysis provides an effective means to simulate the spectra for several different optical conditions with a minimum number of quantum dynamical propagations. The G-MCTDH method is used to calculate and interpret the time-resolved coherent Raman spectra of two selected initial superpositions for a I₂Kr₁₈ cluster embedded in a frozen Kr cage. The time- and frequency-dependent signals carry information about the molecular mechanisms of dissipation and decoherence, which involve vibrational energy transfer to the stretching mode of the four "belt" Kr atoms. The details of these processes and the number of active solvent modes depend in a non-trivial way on the specific initial superposition.

Alternative view of two-dimensional spectroscopy

Femtosecond two-dimensional (2D) spectroscopy has become a widely employed method for the investigation of the dynamics of complex chemical and biological systems. In 2D spectroscopy, the sample is excited with three phase-locked femtosecond pulses, and the signal is heterodyned with the local oscillator field. The 2D spectrum is obtained by double Fourier transform with respect to the time delay between the first two pulses and the time delay between the third pulse and the local oscillator field. We show that 2D optical signals can alternatively be measured and computationally simulated as four-wave-mixing signals generated by two femtosecond pulses and two one-sided continuous-wave (CW) pulses. The first femtosecond pulse and one-sided CW pulse create the doorway state, while the second femtosecond pulse and one-sided CW pulse create the window state. This picture relates 2D spectroscopy to other mixed time-frequency-domain techniques, which is useful for the interpretation of the corresponding signals. Moreover, it allows a computationally efficient evaluation of 2D spectra.

Quantum process tomography of excitonic dimers from two-dimensional electronic spectroscopy. I. General theory and application to homodimers

Is it possible to infer the time evolving quantum state of a multichromophoric system from a sequence of two-dimensional electronic spectra (2D-ES) as a function of waiting time? Here we provide a positive answer for a tractable model system: a coupled dimer. After exhaustively enumerating the Liouville pathways associated to each peak in the 2D-ES, we argue that by judiciously combining the information from a series of experiments varying the polarization and frequency components of the pulses, detailed information at the amplitude level about the input and output quantum states at the waiting time can be obtained. This possibility yields a quantum process tomography (QPT) of the single-exciton manifold, which completely characterizes the open quantum system dynamics through the reconstruction of the process matrix. In this manuscript, we present the general theory as well as specific and numerical results for a homodimer, for which we prove that signals stemming from coherence to population transfer and vice versa vanish upon isotropic averaging, therefore, only allowing for a partial QPT in such case. However, this fact simplifies the spectra, and it follows that only two polarization controlled experiments (and no pulse-shaping requirements) suffice to yield the elements of the process matrix, which survive under isotropic averaging. Redundancies in the 2D-ES amplitudes allow for the angle between the two site transition dipole moments to be self-consistently obtained, hence simultaneously yielding structural and dynamical information of the dimer. Model calculations are presented, as well as an error analysis in terms of the angle between the dipoles and peak amplitude extraction. In the second article accompanying this study, we numerically exemplify the theory for heterodimers and carry out a detailed error analysis for such case. This investigation reveals an exciting quantum information processing (QIP) approach to spectroscopic experiments of excitonic systems, and hence, bridges an important gap between theoretical studies on excitation energy transfer from the QIP standpoint and experimental methods to study such systems in the chemical physics community.

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.