Efficient calculation of the polarization induced by N coherent laser pulses

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.

Plenty of Room on the Top: Pathways and Spectroscopic Signatures of Singlet Fission from Upper Singlet States

We investigate dynamic signatures of the singlet fission (SF) process triggered by the excitation of a molecular system to an upper singlet state S_N (N > 1) and develop a computational methodology for the simulation of nonlinear spectroscopic signals revealing the S_N → TT₁ SF in real time. We demonstrate that SF can proceed directly from the upper state S_N, bypassing the lowest excited state, S₁. We determine the main S_N → TT₁ reaction pathways and show by computer simulation and spectroscopic measurements that the S_N-initiated SF can be faster and more efficient than the traditionally studied S₁ → TT₁ SF. We claim that the S_N → TT₁ SF offers novel promising opportunities for engineering SF systems and enhancing SF yields.

Time-dependent simulation of photocurrent-detected two-dimensional spectroscopy of open systems

A new time-domain simulation protocol of two-dimensional electronic spectra with photocurrent detection is presented. Time-dependent density functional theory for open systems at finite temperature is applied to evaluate the photocurrent response to four laser pulses, and a non-perturbative phase-matching approach is implemented to extract the fourth-order photocurrent signal with a desired phase-matching condition. Simulations for an open three-level model indicates that transition dipoles interact resonantly with the incident pulses and that different sample-electrode couplings may be identified by appearance of different peaks/valleys in photocurrent spectra from different electrodes. Moreover, qualitative reproduction of experimental spectra of a PbS quantum dot photocell [Karki et al., Nat. Commun. 5(1), 5869 (2014)] reveals the stimulated electron dynamics.

Simulation of Time- and Frequency-Resolved Four-Wave-Mixing Signals at Finite Temperatures: A Thermo-Field Dynamics Approach

We propose a new approach to simulate four-wave-mixing signals of molecular systems at finite temperatures by combining the multiconfigurational Ehrenfest method with the thermo-field dynamics theory. In our approach, the four-time correlation functions at finite temperatures are mapped onto those at zero temperature in an enlarged Hilbert space with twice the vibrational degrees of freedom. As an illustration, we have simulated three multidimensional spectroscopic signals, time- and frequency-resolved fluorescence spectra, transient-absorption pump-probe spectra, and electronic two-dimensional (2D) spectra at finite temperatures, for a conical intersection-mediated singlet fission model of a rubrene crystal. It is shown that a detailed dynamical picture of the singlet fission process can be extracted from the three spectroscopic signals. An increasing temperature leads to lower intensities of the signals and broadened vibrational peaks, which can be attributed to faster singlet-triplet population transfer and stronger bath-induced electronic dephasing at higher temperatures.

Efficient numerical method for predicting nonlinear optical spectroscopies of open systems

Nonlinear optical spectroscopies are powerful tools for probing quantum dynamics in molecular and nanoscale systems. While intuition about ultrafast spectroscopies is often built by considering impulsive optical pulses, actual experiments have finite-duration pulses, which can be important for interpreting and predicting experimental results. We present a new freely available open source method for spectroscopic modeling, called Ultrafast Ultrafast (UF²) spectroscopy, which enables computationally efficient and convenient prediction of nonlinear spectra, such as treatment of arbitrary finite duration pulse shapes. UF² is a Fourier-based method that requires diagonalization of the Liouvillian propagator of the system density matrix. We also present a Runge-Kutta-Euler (RKE) direct propagation method. We include open system dynamics in the secular Redfield, full Redfield, and Lindblad formalisms with Markovian baths. For non-Markovian systems, the degrees of freedom corresponding to memory effects are brought into the system and treated nonperturbatively. We analyze the computational complexity of the algorithms and demonstrate numerically that, including the cost of diagonalizing the propagator, UF² is 20-200 times faster than the direct propagation method for secular Redfield models with arbitrary Hilbert space dimension; it is similarly faster for full Redfield models at least up to system dimensions where the propagator requires more than 20 GB to store; and for Lindblad models, it is faster up to Hilbert space dimension near 100 with speedups for small systems by factors of over 500. UF² and RKE are part of a larger open source Ultrafast Software Suite, which includes tools for automatic generation and calculation of Feynman diagrams.

STM Imaging of Electron Migration in Real Space and Time: A Simulation Study

Using a simulation protocol that mimics ultrafast scanning tunneling microscopy (STM) experiments, we demonstrate how pump-probe ultrafast STM may be used to image electron migration in molecules. Two pulses are applied to a model system, and the time-integrated current through the tip is calculated versus the delay time and tip position to generate STM images. With suitable pump and probe parameters, the images can track charge migration with atomistic spatial and femtosecond temporal resolutions.

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.

Femtosecond stimulated Raman spectroscopy as a tool to detect molecular vibrations in ground and excited electronic states

We give a detailed theoretical analysis of the simplest variant of femtosecond stimulated Raman spectroscopy, where a picosecond Raman pump pulse and a femtosecond Raman probe pulse are applied resonantly to a chromophore in thermal equilibrium in the ground electronic state. We demonstrate that this technique is capable of the detection of dephasing-free Raman-like lines revealing vibrational modes not only in the electronic ground state but also in the excited electronic state of the chromophore. The analytical results obtained with simplifying assumptions for the shape of the laser pulses are substantiated by numerical simulations with realistic laser pulses, employing the equation-of-motion phase-matching approach.

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.

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.