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

Vibrational coherences in half-broadband 2D electronic spectroscopy: Spectral filtering to identify excited state displacements

Vibrational coherences in ultrafast pump-probe (PP) and 2D electronic spectroscopy (2DES) provide insights into the excited state dynamics of molecules. Femtosecond coherence spectra and 2D beat maps yield information about displacements of excited state surfaces for key vibrational modes. Half-broadband 2DES uses a PP configuration with a white light continuum probe to extend the detection range and resolve vibrational coherences in the excited state absorption (ESA). However, the interpretation of these spectra is difficult as they are strongly dependent on the spectrum of the pump laser and the relative displacement of the excited states along the vibrational coordinates. We demonstrate the impact of these convoluting factors for a model based upon cresyl violet. A careful consideration of the position of the pump spectrum can be a powerful tool in resolving the ESA coherences to gain insights into excited state displacements. This paper also highlights the need for caution in considering the spectral window of the pulse when interpreting these spectra.

Floquet non-equilibrium Green's function and Floquet quantum master equation for electronic transport: The role of electron-electron interactions and spin current with circular light

The non-equilibrium Green's function (NEGF) and quantum master equation (QME) are two main classes of approaches for electronic transport. We discuss various Floquet variances of these formalisms for transport properties of a quantum dot driven via interaction with an external periodic field. We first derived two versions of the Floquet NEGF. We also explore an ansatz of the Floquet NEGF formalism for the interacting systems. In addition, we derived two versions of Floquet QME in the weak interaction regime. With each method, we elaborate on the evaluation of the expectation values of the number and current operators. We examined these methods for transport through a two-level system that is subject to periodic driving. The numerical results of all four methods show good agreement for non-interacting systems in the weak regime. Furthermore, we have observed that circular light can introduce spin current. We expect these Floquet quantum transport methods to be useful in studying molecular junctions exposed to light.

Probing avoided crossings and conical intersections by two-pulse femtosecond stimulated Raman spectroscopy: Theoretical study

This study leverages two-pulse femtosecond stimulated Raman spectroscopy (2FSRS) to characterize molecular systems with avoided crossings (ACs) and conical intersections (CIs) in their low-lying excited electronic states. By simulating 2FSRS spectra of microscopically inspired ACs and CIs models, we demonstrate that 2FSRS not only delivers valuable information on the molecular parameters characterizing ACs and CIs but also helps distinguish between these two systems.

Nonlinear Molecular Electronic Spectroscopy via MCTDH Quantum Dynamics: From Exact to Approximate Expressions

We present an accurate and efficient approach to computing the linear and nonlinear optical spectroscopy of a closed quantum system subject to impulsive interactions with an incident electromagnetic field. It incorporates the effect of ultrafast nonadiabatic dynamics by means of explicit numerical propagation of the nuclear wave packet. The fundamental expressions for the evaluation of first- and higher-order response functions are recast in a general form that can be used with any quantum dynamics code capable of computing the overlap of nuclear wave packets evolving in different states. Here we present the evaluation of these expressions with the multiconfiguration time-dependent Hartree (MCTDH) method. Application is made to pyrene, excited to its lowest bright excited state S₂ which exhibits a sub-100-fs nonadiabatic decay to a dark state S₁. The system is described by a linear vibronic coupling Hamiltonian, parametrized with multiconfiguration electronic structure methods. We show that the ultrafast nonadiabatic dynamics can have a remarkable effect on the spectral line shapes that goes beyond simple lifetime broadening. Furthermore, a widely employed approximate expression based on the time scale separation of dephasing and population relaxation is recast in the same theoretical framework. Application to pyrene shows the range of validity of such approximations.

The hierarchy of Davydov's Ansätze: From guesswork to numerically "exact" many-body wave functions

This Perspective presents an overview of the development of the hierarchy of Davydov's Ansätze and a few of their applications in many-body problems in computational chemical physics. Davydov's solitons originated in the investigation of vibrational energy transport in proteins in the 1970s. Momentum-space projection of these solitary waves turned up to be accurate variational ground-state wave functions for the extended Holstein molecular crystal model, lending unambiguous evidence to the absence of formal quantum phase transitions in Holstein systems. The multiple Davydov Ansätze have been proposed, with increasing Ansatz multiplicity, as incremental improvements of their single-Ansatz parents. For a given Hamiltonian, the time-dependent variational formalism is utilized to extract accurate dynamic and spectroscopic properties using Davydov's Ansätze as its trial states. A quantity proven to disappear for large multiplicities, the Ansatz relative deviation is introduced to quantify how closely the Schrödinger equation is obeyed. Three finite-temperature extensions to the time-dependent variation scheme are elaborated, i.e., the Monte Carlo importance sampling, the method of thermofield dynamics, and the method of displaced number states. To demonstrate the versatility of the methodology, this Perspective provides applications of Davydov's Ansätze to the generalized Holstein Hamiltonian, variants of the spin-boson model, and systems of cavity-assisted singlet fission, where accurate dynamic and spectroscopic properties of the many-body systems are given by the Davydov trial states.

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.

Effect of mixed Frenkel and charge transfer states in time-gated fluorescence spectra of perylene bisimides H-aggregates: Hierarchical equations of motion approach

We theoretically investigated the effect of mixed Frenkel (F) and charge transfer (CT) states on the spectral properties of perylene bisimide (PBI) derivatives, focusing on the role of strong electron-phonon interactions. The model consists of a four-level system described by the Holstein Hamiltonian coupled to independent local heat-baths on each site, described by Brownian spectral distribution functions. We employ the reduced hierarchical equations of motion (HEOM) approach to calculate the time evolution of the system and compare it to the pure F exciton cases. We compute the absorption and time-gated fluorescence (TGF) spectra for different exciton transfer integrals and F-CT band gap conditions. The coherence length of excitons ($N_{coh}$) is evaluated employing two different definitions. We observe the presence of an excited hot state peak whose intensity is associated with the delocalization of the excited species and ultrafast dynamics that are solely dependent on the frequency of the local bath. The results indicate that the inclusion of CT states promotes localization of the excitons which is manifested in a decrease in the intensity of the hot state peak and the 0--1 peak, and an increase in the intensity of the 0--0 emission peak in TGF spectrum, leading to a decrease of $N_{coh}$.

Optical Cavity Manipulation and Nonlinear UV Molecular Spectroscopy of Conical Intersections in Pyrazine

Optical cavities provide a versatile platform for manipulating the excited-state dynamics of molecules via strong light-matter coupling. We employ optical absorption and two-multidimensional electronic spectroscopy simulations to investigate the effect of optical cavity coupling in the nonadiabatic dynamics of photoexcited pyrazine. We observe the emergence of a novel polaritonic conical intersection (PCI) between the electronic dark state and photonic surfaces as the cavity frequency is tuned. The PCI could significantly change the nonadiabatic dynamics of pyrazine by doubling the decay rate constant of the S₂ state population. Moreover, the absorption spectrum and excited-state dynamics could be systematically manipulated by tuning the strong light-matter interaction, e.g., the cavity frequency and cavity coupling strength. We propose that a tunable optical cavity-molecule system may provide promising approaches for manipulating the photophysical properties of molecules.

Ultrafast Internal Conversion Dynamics through the on-the-Fly Simulation of Transient Absorption Pump-Probe Spectra with Different Electronic Structure Methods

An on-the-fly surface-hopping simulation protocol is developed for the evaluation of transient absorption (TA) pump-probe (PP) signals of molecular systems exhibiting internal conversion to the electronic ground state. We study the nonadiabatic dynamics of azomethane and the associating TA PP spectra at three levels of the electronic-structure theory, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. The impact of these methods on the population dynamics and time-resolved TA PP signals is substantially different. This difference is attributed to the strong non-Condon effects that must be taken into account for the proper understanding and interpretation of time-resolved TA PP signals of nonadiabatic polyatomic systems. This shows that the combination of the dynamical and spectral simulations definitely provides more accurate and detailed information on the microscopic mechanisms of photophysical and photochemical processes. Hence the simulation of time-resolved spectroscopic signals provides another important dimension to examine the accuracy of quantum chemistry methods.

Signatures of coherent vibronic exciton dynamics and conformational control in two-dimensional electronic spectroscopy of conjugated polymers

Two-dimensional electronic spectroscopy (2DES) signals for homo-oligomer J-aggregates are computed, with a focus on the role of structural change induced by low-frequency torsional modes along with quasi-stationary trapping effects induced...

Simulation of Nonlinear Femtosecond Signals at Finite Temperature via a Thermo Field Dynamics-Tensor Train Method: General Theory and Application to Time- and Frequency-Resolved Fluorescence of the Fenna-Matthews-Olson Complex

Addressing needs of contemporary nonlinear femtosecond optical spectroscopy, we have developed a fully quantum, numerically accurate wave function-based approach for the calculation of third-order spectroscopic signals of polyatomic molecules and molecular aggregates at finite temperature. The systems are described by multimode nonadiabatic vibronic-coupling Hamiltonians, in which diagonal terms are treated in harmonic approximation, while off-diagonal interstate couplings are assumed to be coordinate independent. The approach is based on the Thermo Field Dynamics (TFD) representation of quantum mechanics and tensor-train (TT) machinery for efficient numerical simulation of quantum evolution of systems with many degrees of freedom. The developed TFD-TT approach is applied to the calculation of time- and frequency-resolved fluorescence spectra of the Fenna-Matthews-Olson (FMO) antenna complex at room temperature taking into account finite time-frequency resolution in fluorescence detection, orientational averaging, and static disorder.

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.

Theory of two-dimensional spectroscopy with intense laser fields

Two-dimensional vibrational and electronic spectroscopic observables of isotropically oriented molecular samples in solution are sensitive to laser field intensities and polarization. The third-order response function formalism predicts a signal that grows linearly with the field strength of each laser pulse, thus lacking a way of accounting for non-trivial intensity-dependent effects, such as saturation and finite bleaching. An analytical expression to describe the orientational part of the molecular response, which, in the weak-field limit, becomes equivalent to a four-point correlation function, is presented. Such an expression is evaluated for Liouville-space pathways accounting for diagonal and cross peaks for all-parallel and cross-polarized pulse sequences, in both the weak- and strong-field conditions, via truncation of a Taylor series expansion at different orders. The results obtained in the strong-field conditions suggest how a careful analysis of two-dimensional spectroscopic experimental data should include laser pulse intensity considerations when determining molecular internal coordinates.

Ab Initio Surface-Hopping Simulation of Femtosecond Transient-Absorption Pump-Probe Signals of Nonadiabatic Excited-State Dynamics Using the Doorway-Window Representation

An ab initio theoretical framework for the simulation of femtosecond time-resolved transient absorption (TA) pump-probe (PP) spectra with quasi-classical trajectories is presented. The simulations are based on the classical approximation to the doorway-window (DW) representation of third-order four-wave-mixing signals. The DW formula accounts for the finite duration and spectral shape of the pump and probe pulses. In the classical DW formalism, classical trajectories are stochastically sampled from a positive definite doorway distribution, and the signals are evaluated by averaging over a positive definite window distribution. Nonadiabatic excited-state dynamics is described by a stochastic surface-hopping algorithm. The method has been implemented for the pyrazine molecule with the second-order algebraic-diagrammatic construction (ADC(2)) ab initio electronic-structure method. The methodology is illustrated by ab initio simulations of the ground-state bleach, stimulated emission, and excited-state absorption contributions to the TA PP spectrum of gas-phase pyrazine.

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.

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.

Quantum dynamics and spectroscopy of dihalogens in solid matrices. I. Efficient simulation of the photodynamics of the embedded I2Kr18 cluster using the G-MCTDH method

The molecular dynamics following the electronic BΠu30⁺⟵XΣg+1 photoexcitation of the iodine molecule embedded in solid krypton are studied quantum mechanically using the Gaussian variant of the multiconfigurational time-dependent Hartree method (G-MCTDH). The accuracy of the Gaussian wave packet approximation is validated against numerically exact MCTDH simulations for a fully anharmonic seven-dimensional model of the I₂Kr₁₈ cluster in a crystal Kr cage. The linear absorption spectrum, time-evolving vibrational probability densities, and I₂ energy expectation value are accurately reproduced by the numerically efficient G-MCTDH approach. The reduced density matrix of the chromophore is analyzed in the coordinate, Wigner and energy representations, so as to obtain a multifaceted dynamical view of the guest-host interactions. Vibrational coherences extending over the bond distance range 2.7 Å < R_I-I < 4.0 Å are found to survive for several vibrational periods, despite extensive dissipation. The present results prepare the ground for the simulation of time-resolved coherent Raman spectroscopy of the I₂-krypton system addressed in Paper II.

One- to Two-Exciton Transitions in Perylene Bisimide Dimer Revealed by Two-Dimensional Electronic Spectroscopy

The excited-state energy levels of molecular dimers and aggregates play a critical role in their photophysical behavior and an understanding of the photodynamics in such structures is important for developing applications such as photovoltaics and optoelectronic devices. Here, exciton transitions in two different covalently bound PBI dimers are studied by two-dimensional electronic spectroscopy (2DES), a powerful spectroscopic method, providing the most complete picture of vibronic transitions in molecular systems. The data are accurately reproduced using the equation of motion-phase matching approach. The unambiguous presence of one-exciton to two-exciton transitions are captured in our results and described in terms of a molecular exciton energy level scheme based on the Kasha model. Furthermore, the results are supported by comparative measurements with the PBI monomer and another dimer in which the interchromophore distance is increased.

Spectral Filtering as a Tool for Two-Dimensional Spectroscopy: A Theoretical Model

Two-dimensional optical spectroscopy is a powerful technique for the probing of coherent quantum superpositions. Recently, the finite width of the laser spectrum has been employed to selectively tune experiments for the study of particular coherences. This involves the exclusion of certain transition frequencies, which results in the elimination of specific Liouville pathways. The rigorous analysis of such experiments requires the use of ever more sophisticated theoretical models for the optical spectroscopy of electronic and vibronic systems. Here we develop a nonimpulsive and non-Markovian model, which combines an explicit definition of the laser spectrum, via the equation of motion-phase matching approach (EOM-PMA), with the hierarchical equations of motion (HEOM). This theoretical framework is capable of simulating the 2D spectroscopy of vibronic systems with low frequency modes, coupled to environments of intermediate and slower time scales. In order to demonstrate the spectral filtering of vibronic coherences, we examine the elimination of lower energy peaks from the 2D spectra of a zinc porphyrin monomer upon blue-shifting the laser spectrum. The filtering of Liouville pathways is revealed through the disappearance of peaks from the amplitude spectra for a coupled vibrational mode.

Coherent two-dimensional electronic mass spectrometry

Coherent two-dimensional (2D) optical spectroscopy has revolutionized our ability to probe many types of couplings and ultrafast dynamics in complex quantum systems. The dynamics and function of any quantum system strongly depend on couplings to the environment. Thus, studying coherent interactions for different environments remains a topic of tremendous interest. Here we introduce coherent 2D electronic mass spectrometry that allows 2D measurements on effusive molecular beams and thus on quantum systems with minimum system–bath interaction and employ this to identify the major ionization pathway of 3d Rydberg states in NO₂. Furthermore, we present 2D spectra of multiphoton ionization, disclosing distinct differences in the nonlinear response functions leading to the ionization products. We also realize the equivalent of spectrally resolved transient-absorption measurements without the necessity for acquiring weak absorption changes. Using time-of-flight detection introduces cations as an observable, enabling the 2D spectroscopic study on isolated systems of photophysical and photochemical reactions.

Multidimensional spectroscopy is a powerful tool in exploring photo-induced dynamics and electron coupling processes in molecules. Here the authors demonstrate coherent two-dimensional electronic mass spectrometry on molecular beams and its application to photoionization studies of the NO₂ molecule.

Finite temperature dynamics of a Holstein polaron: The thermo-field dynamics approach

Combining the multiple Davydov D₂ Ansatz with the method of thermo-field dynamics, we study finite temperature dynamics of a Holstein polaron on a lattice. It has been demonstrated, using the hierarchy equations of motion method as a benchmark, that our approach provides an efficient, robust description of finite temperature dynamics of the Holstein polaron in the simultaneous presence of diagonal and off-diagonal exciton-phonon coupling. The method of thermo-field dynamics handles temperature effects in the Hilbert space with key numerical advantages over other treatments of finite-temperature dynamics based on quantum master equations in the Liouville space or wave function propagation with Monte Carlo importance sampling. While for weak to moderate diagonal coupling temperature increases inhibit polaron mobility, it is found that off-diagonal coupling induces phonon-assisted transport that dominates at high temperatures. Results on the mean square displacements show that band-like transport features dominate the diagonal coupling cases, and there exists a crossover from band-like to hopping transport with increasing temperature when including off-diagonal coupling. As a proof of concept, our theory provides a unified treatment of coherent and incoherent transport in molecular crystals and is applicable to any temperature.

Deconvolution of optical multidimensional coherent spectra

Optical coherent multidimensional spectroscopy is a powerful technique for unraveling complex and congested spectra by spreading them across multiple dimensions, removing the effects of inhomogeneity, and revealing underlying correlations. As the technique matures, the focus is shifting from understanding the technique itself to using it to probe the underlying dynamics in the system being studied. However, these dynamics can be difficult to discern because they are convolved with the nonlinear optical response of the system. Inspired by methods used to deblur images, we present a method for deconvolving the underlying dynamics from the optical response. To demonstrate the method, we extract the many-particle diffusion Green's functions for excitons in a semiconductor quantum well from two-dimensional coherent spectra.

Phase Cycling RT-TDDFT Simulation Protocol for Nonlinear XUV and X-ray Molecular Spectroscopy

Real-time time-dependent density functional theory (RT-TDDFT) provides a practical algorithm for propagating a many-electron system driven by external laser fields. The fields are included nonperturbatively in the propagation, and the molecular reduced single-electron density operator and various spectroscopic and diffraction signals can be computed directly, avoiding the expensive calculation of many-body states. Nonlinear optical signals contain contributions of multiple pathways. A phase cycling protocol is implemented in order to separate these pathways. Simulations of XUV four-wave mixing signals in the CO molecule are compared with ab initio sum-over-states calculations.

Frequency-domain coherent multidimensional spectroscopy when dephasing rivals pulsewidth: Disentangling material and instrument response

Ultrafast spectroscopy is often collected in the mixed frequency/time domain, where pulse durations are similar to system dephasing times. In these experiments, expectations derived from the familiar driven and impulsive limits are not valid. This work simulates the mixed-domain four-wave mixing response of a model system to develop expectations for this more complex field-matter interaction. We explore frequency and delay axes. We show that these line shapes are exquisitely sensitive to excitation pulse widths and delays. Near pulse overlap, the excitation pulses induce correlations that resemble signatures of dynamic inhomogeneity. We describe these line shapes using an intuitive picture that connects to familiar field-matter expressions. We develop strategies for distinguishing pulse-induced correlations from true system inhomogeneity. These simulations provide a foundation for interpretation of ultrafast experiments in the mixed domain.

Applications of the New Family of Coherent Multidimensional Spectroscopies for Analytical Chemistry

A new family of vibrational and electronic spectroscopies has emerged, comprising the coherent analogs of traditional analytical methods. These methods are also analogs of coherent multidimensional nuclear magnetic resonance (NMR) spectroscopy. This new family is based on creating the same quantum mechanical superposition states called multiple quantum coherences (MQCs). NMR MQCs are mixtures of nuclear spin states that retain their quantum mechanical phase information for milliseconds. The MQCs in this new family are mixtures of vibrational and electronic states that retain their phases for picoseconds or shorter times. Ultrafast, high-intensity coherent beams rapidly excite multiple states. The excited MQCs then emit bright beams while they retain their phases. Time-domain methods measure the frequencies of the MQCs by resolving their phase oscillations, whereas frequency-domain methods measure the resonance enhancements of the output beam while scanning the excitation frequencies. The resulting spectra provide multidimensional spectral signatures that increase the spectroscopic selectivity required for analyzing complex samples.

Quantum electron-vibrational dynamics at finite temperature: Thermo field dynamics approach

Quantum electron-vibrational dynamics in molecular systems at finite temperature is described using an approach based on the thermo field dynamics theory. This formulation treats temperature effects in the Hilbert space without introducing the Liouville space. A comparison with the theoretically equivalent density matrix formulation shows the key numerical advantages of the present approach. The solution of thermo field dynamics equations with a novel technique for the propagation of tensor trains (matrix product states) is discussed. Numerical applications to model spin-boson systems show that the present approach is a promising tool for the description of quantum dynamics of complex molecular systems at finite temperature.

Impact of Vibrational Coherence on the Quantum Yield at a Conical Intersection

We study the vibrationally coherent quantum dynamics of an electronic wave packet in the vicinity of a conical intersection within a three-state two-mode model. By transforming the coherent tuning and coupling modes into the bath, the underdamped dynamics of the resulting effective three-state model is solved efficiently by the numerically exact hierarchy equation of motion approach. The transient excited-state absorption and two-dimensional spectra reveal the impact of vibrational coherence on the relaxation pathways of the wave packet. We find that both the quantum yield and the isomerization rate are crucially influenced by the vibrational coherence of the wave packet. A less coherent wave packet can traverse the conical intersection more rapidly, while the resulting quantum yield is smaller. Finally, we show that repeated passages of the wave packet through the conical intersection can lead to measurable interference effects in the form of Stueckelberg oscillations.

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.