Numerical method for nonlinear optical spectroscopies: Ultrafast ultrafast spectroscopy

A Quantum Algorithm from Response Theory: Digital Quantum Simulation of Two-Dimensional Electronic Spectroscopy

Multidimensional optical spectroscopies are powerful techniques to investigate energy transfer pathways in natural and artificial systems. Because of the high information content of the spectra, numerical simulations of the optical response are of primary importance to assist the interpretation of spectral features. However, the increasing complexity of the investigated systems and their quantum dynamics call for the development of novel simulation strategies. In this work, we consider using digital quantum computers. By combining quantum dynamical simulation and nonlinear response theory, we present a quantum algorithm for computing the optical response of molecular systems. The quantum advantage stems from the efficient quantum simulation of the dynamics governed by the molecular Hamiltonian, and it is demonstrated by explicitly considering exciton-vibrational coupling. The protocol is tested on a near-term quantum device, providing the digital quantum simulation of the linear and nonlinear response of simple molecular models.

Interpretations of High-Order Transient Absorption Spectroscopies

Transient absorption (TA) spectroscopy is an invaluable tool for determining the energetics and dynamics of excited states. When pump intensities are sufficiently high, TA spectra include both the generally desired third-order response and responses that are higher in order in the field amplitudes. Recent work demonstrated that pump-intensity-dependent TA measurements allow separating the orders of response, but the information content in those higher orders has not been described. We give a general framework for understanding high-order TA spectra. We extend to higher order the fundamental processes of standard TA: ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA). Each order introduces two new processes: SE and ESA from previously inaccessible highly excited states and negations of lower-order processes. We show the new spectral and dynamical information at each order and show how the relative signs of the signals in different orders can be used to identify which processes dominate.

Higher-Order Multidimensional and Pump-Probe Spectroscopies

Transient absorption and coherent two-dimensional spectroscopy are widely established methods for the investigation of ultrafast dynamics in quantum systems. Conventionally, they are interpreted in the framework of perturbation theory at the third order of interaction. Here, we discuss the potential of higher-(than-third-)order pump-probe and multidimensional spectroscopy to provide insight into excited multiparticle states and their dynamics. We focus on recent developments from our group. In particular, we demonstrate how phase cycling can be used in fluorescence-detected two-dimensional spectroscopy to isolate higher-order spectra that provide information about highly excited states such as the correlation of multiexciton states. We discuss coherently detected fifth-order 2D spectroscopy and its power to track exciton diffusion. Finally, we show how to extract higher-order signals even from ordinary pump-probe experiments, providing annihilation-free signals at high excitation densities and insight into multiexciton interactions.

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.

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.

Two-dimensional spectroscopy beyond the perturbative limit: The influence of finite pulses and detection modes

Ultra-fast and multi-dimensional spectroscopy gives a powerful looking glass into the dynamics of molecular systems. In particular, two-dimensional electronic spectroscopy (2DES) provides a probe of coherence and the flow of energy within quantum systems, which is not possible with more conventional techniques. While heterodyne-detected (HD) 2DES is increasingly common, more recently fluorescence-detected (FD) 2DES offers new opportunities, including single-molecule experiments. However, in both techniques, it can be difficult to unambiguously identify the pathways that dominate the signal. Therefore, the use of numerically modeling of 2DES is vitally important, which, in turn, requires approximating the pulsing scheme to some degree. Here, we employ non-perturbative time evolution to investigate the effects of finite pulse width and amplitude on 2DES signals. In doing so, we identify key differences in the response of HD and FD detection schemes, as well as the regions of parameter space where the signal is obscured by unwanted artifacts in either technique. Mapping out parameter space in this way provides a guide to choosing experimental conditions and also shows in which limits the usual theoretical approximations work well and in which limits more sophisticated approaches are required.

Automatic Feynman diagram generation for nonlinear optical spectroscopies and application to fifth-order spectroscopy with pulse overlaps

Perturbative nonlinear optical spectroscopies are powerful methods to understand the dynamics of excitonic and other condensed phase systems. Feynman diagrams have long provided the essential tool to understand and interpret experimental spectra and to organize the calculation of spectra for model systems. When optical pulses are strictly time ordered, only a small number of diagrams contribute, but in many experiments, pulse-overlap effects are important for interpreting results. When pulses overlap, the number of contributing diagrams can increase rapidly, especially with higher order spectroscopies, and human error is especially likely when attempting to write down all the diagrams. We present an automated Diagram Generator (DG) that generates all the Feynman diagrams needed to calculate any nth-order spectroscopic signal. We characterize all perturbative nonlinear spectroscopies by their associated phase-discrimination condition as well as the time intervals where pulse amplitudes are nonzero. Although the DG can be used to automate impulsive calculations, its greatest strength lies in automating finite pulse calculations where pulse overlaps are important. We consider third-order transient absorption spectroscopy and fifth-order exciton-exciton interaction 2D (EEI2D) spectroscopy, which are described by six or seven diagrams in the impulsive limit, respectively, but 16 or 240 diagrams, respectively, when pulses overlap. The DG allows users to automatically include all relevant diagrams at a relatively low computational cost, since the extra diagrams are only generated for the inter-pulse delays where they are relevant. For EEI2D spectroscopy, we show the important effects of including the overlap diagrams.

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.