Nonperturbative quantum simulation of time-resolved nonlinear spectra: Methodology and application to electron transfer reactions in the condensed phase

Smolyak Scheme for solving the Schrödinger equation: Application to Malonaldehyde in Full Dimensionality

In 1963 Smolyak introduced an approach to overcome the exponential scaling with respect to the number of variables of the direct product size [S. A. Smolyak Soviet Mathematics Doklady, 4, 240 (1963)]. The main idea is to replace a single large direct product by a sum of selected small direct products. It was first used in quantum dynamics in 2009 by Avila and Carrington [G. Avila and T. Carrington, J. Chem. Phys., 131, 174103 (2009)]. Since then, several calculations have been published by Avila and Carrington and by other groups. In the present study, and to push the limit to larger and more complex systems, this scheme is combined with the use of an on-the-fly calculation of the kinetic energy operator and a Block-Davidson procedure to obtain eigenstates in our home-made Fortran codes, ElVibRot and Tnum-Tana. This was applied to compute the tunneling splitting of malonaldehyde in full dimensionality (21D) using the potential of Mizukami et al. [W. Mizukami, S. Habershon, and D.P. Tew, J. Chem. Phys. 141, 1443-10 (2014)]. Our tunneling splitting calculations, 21.7±0.3 cm-1 and 2.9±0.1 cm-1 , show excellent agreement with the experimental values, 21.6 cm-1 and 2.9 cm-1 for the normal isotopologue and the mono-deuterated one, respectively.

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.

Efficient quantum dynamics simulations of complex molecular systems: A unified treatment of dynamic and static disorder

We present a unified and highly numerically efficient formalism for the simulation of quantum dynamics of complex molecular systems, which takes into account both temperature effects and static disorder. The methodology is based on the thermo-field dynamics formalism, and Gaussian static disorder is included into simulations via auxiliary bosonic operators. This approach, combined with the tensor-train/matrix-product state representation of the thermalized stochastic wave function, is applied to study the effect of dynamic and static disorders in charge-transfer processes in model organic semiconductor chains employing the Su-Schrieffer-Heeger (Holstein-Peierls) model Hamiltonian.

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.

Non-adiabatic quantum dynamics without potential energy surfaces based on second-quantized electrons: Application within the framework of the MCTDH method

A first principles quantum formalism to describe the non-adiabatic dynamics of electrons and nuclei based on a second quantization representation (SQR) of the electronic motion combined with the usual representation of the nuclear coordinates is introduced. This procedure circumvents the introduction of potential energy surfaces and non-adiabatic couplings, providing an alternative to the Born-Oppenheimer approximation. An important feature of the molecular Hamiltonian in the mixed first quantized representation for the nuclei and the SQR representation for the electrons is that all degrees of freedom, nuclear positions and electronic occupations, are distinguishable. This makes the approach compatible with various tensor decomposition Ansätze for the propagation of the nuclear-electronic wavefunction. Here, we describe the application of this formalism within the multi-configuration time-dependent Hartree framework and its multilayer generalization, corresponding to Tucker and hierarchical Tucker tensor decompositions of the wavefunction, respectively. The approach is applied to the calculation of the photodissociation cross section of the HeH⁺ molecule under extreme ultraviolet irradiation, which features non-adiabatic effects and quantum interferences between the two possible fragmentation channels, He + H⁺ and He⁺ + H. These calculations are compared with the usual description based on ab initio potential energy surfaces and non-adiabatic coupling matrix elements, which fully agree. The proof-of-principle calculations serve to illustrate the advantages and drawbacks of this formalism, which are discussed in detail, as well as possible ways to overcome them. We close with an outlook of possible application domains where the formalism might outperform the usual approach, for example, in situations that combine a strong static correlation of the electrons with non-adiabatic electronic-nuclear effects.

Effects of high pulse intensity and chirp in two-dimensional electronic spectroscopy of an atomic vapor

The effects of high pulse intensity and chirp on two-dimensional electronic spectroscopy signals are experimentally investigated in the highly non-perturbative regime using atomic rubidium vapor as clean model system. Data analysis is performed based on higher-order Feynman diagrams and non-perturbative numerical simulations of the system response. It is shown that higher-order contributions may lead to a fundamental change of the static appearance and beating-maps of the 2D spectra and that chirped pulses enhance or suppress distinct higher-order pathways. We further give an estimate of the threshold intensity beyond which the high-intensity effects become visible for the system under consideration.

A multilayer multiconfiguration time-dependent Hartree simulation of the reaction-coordinate spin-boson model employing an interaction picture

The multilayer multiconfiguration time-dependent Hartree method is applied in an interaction picture to simulate dynamics of the spin-boson model in the reaction-coordinate representation. The use of the interaction picture allows a more effective description of correlation effects, especially when the coupling strength between the reaction coordinate and the bath is very strong. Examples show that in most physical regimes the efficiency is improved significantly, in some cases up to several orders of magnitude. This opens up new avenues for studying quantum dynamical problems.

Multilayer Multiconfiguration Time-Dependent Hartree Theory

Multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) theory is a rigorous and powerful method to simulate quantum dynamics in complex many-body systems. This approach extends the original MCTDH theory of Meyer, Manthe, and Cederbaum to include dynamically contracted layers in a recursive way, within which the equations of motion are determined from the Dirac-Frenkel variational principle. This paper presents the general derivation of the theory and analyzes the important features that make the ML-MCTDH method numerically efficient. Furthermore, we discuss the generalization of the theory to treat many-body identical particles (fermions or bosons) as well as calculating energy eigenstates via the improved relaxation method.

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.

Quantum dynamics with sparse grids: a combination of Smolyak scheme and cubature. Application to methanol in full dimensionality

Quantum dynamical approaches based on product-grids are limited to the studies of molecular systems with few degrees of freedom, typically less than ten. Recently, Avila et al. [G. Avila, T. Carrington, J. Chem. Phys., 131 (2009) 174103] have introduced the Smolyak scheme [S.A. Smolyak, Sov. Math. Dokl., 4 (1963) 240], which considerably reduces the size of the grids. This approach has pushed back the present calculation limits on the vibrational spectra of polyatomic molecules. In the present study, we have developed an extension of the standard Smolyak scheme in which this scheme is combined with multidimensional grids, such as cubatures, to obtain new sparse grids. This scheme has been applied to the study of the torsional energy levels of methanol in full dimensionality (12D).

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.

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

"Four-wave-mixing" is the generic name for a family of nonlinear electronic and vibrational spectroscopies. These techniques are widely used to explore dissipation, dephasing, solvation, and interstate coupling mechanisms in various material systems. Four-wave-mixing spectroscopy needs a firm theoretical support, because it delivers information on material systems indirectly, through certain transients, which are measured as functions of carrier frequencies, durations, and relative time delays of the laser pulses. The observed transients are uniquely determined by the three-pulse-induced third-order polarization. There exist two conceptually different approaches to the calculation of the nonlinear polarization. In the standard perturbative approach to nonlinear spectroscopy, the third-order polarization is expressed in terms of the nonlinear response functions. As the material systems become more complex, the evaluation of the response functions becomes cumbersome and the calculation of the signals necessitates a number of approximations. Herein, we review alternative methods for the calculation of four-wave-mixing signals, in which the relevant laser pulses are incorporated into the system Hamiltonian and the driven system dynamics is simulated numerically exactly. The emphasis is on the recently developed equation-of-motion phase-matching approach (EOM-PMA), which allows us to calculate the three-pulse-induced third-order polarization in any phase-matching direction by performing three (with the rotating wave approximation) or seven (without the rotating wave approximation) independent propagations of the density matrix. The EOM-PMA is limited to weak laser fields (its domain of validity is equivalent to the approach based on the third-order response functions) but allows for arbitrary pulse durations and automatically accounts for pulse-overlap effects. As an illustration, we apply the EOM-PMA to the calculation of optical three-pulse photon-echo two-dimensional (2D) signals for a generic model system, which represents a characteristic photophysical dynamics of large molecules or chromophores in condensed phases. The EOM-PMA is easy to implement and can straightforwardly be incorporated into any computational scheme, which provides the time-dependent density matrix or wave function of the material system of interest. In particular, EOM-PMA-based computer codes can efficiently be implemented on parallel computers. The EOM-PMA facilitates considerably the computation of four-wave-mixing signals and 2D spectra, in both vibrational and electronic spectroscopy. The EOM-PMA can be extended to higher order optical responses, e.g., heterodyned 3D IR, transient 2D IR, and other six-wave-mixing techniques.