Nonperturbative non-Markovian quantum master equation: Validity and limitation to calculate nonlinear response functions

Simulation of the Pump-Probe Spectra and Excitation Energy Relaxation of the B850 Band of the LH2 Complex in Purple Bacteria

Ultrafast spectroscopic techniques have been vital in studying excitation energy transfer (EET) in photosynthetic light harvesting complexes. In this paper, we simulate the pump-probe spectra of the B850 band of the light harvesting complex 2 (LH2) of purple bacteria, by using the hierarchical equation of motion method and the optical response function approach. The ground state bleach, stimulated emission, and excited state absorption components of the pump-probe spectra are analyzed in detail. The laser pulse-induced population dynamics are also simulated to help understand the main features of the pump-probe spectra and the EET process. It is shown that the excitation energy relaxation is an ultrafast process with multiple time scales. The first 40 fs of the pump-probe spectra is dominated by the relaxation of the k = ±1 states to both the k = 0 and higher energy states. Dynamics on a longer time scale around 200 fs reflects the relaxation of higher energy states to the k = 0 state.

Role of Quantum Information in HEOM Trajectories

Open quantum systems often operate in the non-Markovian regime where a finite history of a trajectory is intrinsic to its evolution. The degree of non-Markovianity for a trajectory may be measured in terms of the amount of information flowing from the bath back into the system. In this study, we consider how information flows through the auxiliary density operators (ADOs) in the hierarchical equations of motion. We consider three cases for a range of baths, underdamped, intermediate, and overdamped. By understanding how information flows, we are able to determine the relative importance of different ADOs within the hierarchy. We show that ADOs sharing a common Matsubara axis behave similarly, while ADOs on different Matsubara axes behave differently. Using this knowledge, we are able to truncate hierarchies significantly, thus reducing the computation time, while obtaining qualitatively similar results. This is illustrated by comparing 2D electronic spectra for a molecule with an underdamped vibration subsumed into the bath spectral density.

Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions

Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system's complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.

Pathway selectivity in time-resolved spectroscopy using two-photon coincidence counting with quantum entangled photons

Ultrafast optical spectroscopy is a powerful technique for studying the dynamic processes of molecular systems in condensed phases. However, in molecular systems containing many dye molecules, the spectra can become crowded and difficult to interpret owing to the presence of multiple nonlinear optical contributions. In this work, we theoretically propose time-resolved spectroscopy based on the coincidence counting of two entangled photons generated via parametric down-conversion with a monochromatic laser. We demonstrate that the use of two-photon counting detection of entangled photon pairs enables the selective elimination of the excited-state absorption signal. This selective elimination cannot be realized with classical coherent light. We anticipate that the proposed spectroscopy will help simplify the spectral interpretation of complex molecular and material systems comprising multiple molecules.

Phonon Signatures in Photon Correlations

We show that the second-order, two-time correlation functions for phonons and photons emitted from a vibronic molecule in a thermal bath result in bunching and antibunching (a purely quantum effect), respectively. Signatures relating to phonon exchange with the environment are revealed in photon-photon correlations. We demonstrate that cross-correlation functions have a strong dependence on the order of detection giving insight into how phonon dynamics influences the emission of light. This work offers new opportunities to investigate quantum effects in condensed-phase molecular systems.

Quantum dynamics simulations of the 2D spectroscopy for exciton polaritons

We develop an accurate and numerically efficient non-adiabatic path-integral approach to simulate the non-linear spectroscopy of exciton-polariton systems. This approach is based on the partial linearized density matrix approach to model the exciton dynamics with explicit propagation of the phonon bath environment, combined with a stochastic Lindblad dynamics approach to model the cavity loss dynamics. Through simulating both linear and polariton two-dimensional electronic spectra, we systematically investigate how light-matter coupling strength and cavity loss rate influence the optical response signal. Our results confirm the polaron decoupling effect, which is the reduced exciton-phonon coupling among polariton states due to the strong light-matter interactions. We further demonstrate that the polariton coherence time can be significantly prolonged compared to the electronic coherence outside the cavity.

Structural Stability of Insulin Oligomers and Protein Association-Dissociation Processes: Free Energy Landscape and Universal Role of Water

Association and dissociation of proteins are important biochemical events. In this Feature Article, we analyze the available studies of these processes for insulin oligomers in aqueous solution. We focus on the solvation of the insulin monomer in water, stability and dissociation of its dimer, and structural integrity of the hexamer. The intricate role of water in solvation of the dimer- and hexamer-forming surfaces, in long-range interactions between the monomers and the stability of the oligomers, is discussed. Ten water molecules inside the central cavity stabilize the structure of the insulin hexamer. We discuss how different order parameters can be used to understand the dissociation of the insulin dimer. The calculation of the rate using a recently computed multidimensional free energy provides considerable insight into the interplay between protein and water dynamics.

Ultrafast Charge Transfer and Relaxation at a Donor-Acceptor Interface

The efficiency of charge separation in organic photovoltaic materials is crucially determined by the underlying dynamics of the charge transfer (CT) excitons and their dissociation into free electrons and holes. To unravel the main principles of the underlying mechanism on a molecular level, we construct a toy model of electronically coupled donors interacting with a manifold of CT exciton states. In particular, we set up a ladder of CT site energies to model the exciton dissociation. To mimic the complexity of the exciton dynamics at the donor-acceptor interface, the electronic CT manifold is designed to include two vibrational modes that are vibronically coupled to the excitons. We examine the impact of the electronic and vibrational coherences and the structure of the vibronic manifold on the transfer efficiency and charge recombination. Optimal configurations of the vibronic CT manifold are revealed. In particular, the rate of charge recombination can be minimized when the transient dynamics are carefully explored. Such a toy model can be used as a guide for the design of organic materials for efficient photovoltaic devices.

A quantum Langevin equation approach for two-dimensional electronic spectra of coupled vibrational and electronic dynamics

We present an efficient method to simulate two-dimensional (2D) electronic spectra of condensed-phase systems with an emphasis on treating quantum nuclear wave packet dynamics explicitly. To this end, we combine a quantum Langevin equation (QLE) approach for dissipation and a perturbative scheme to calculate three-pulse photon-echo polarizations based on wave packet dynamics under the influence of external fields. The proposed dynamical approach provides a consistent description of nuclear quantum dynamics, pulse-overlap effects, and vibrational relaxation, enabling simulations of 2D electronic spectra with explicit and non-perturbative treatment of coupled electronic-nuclear dynamics. We apply the method to simulate 2D electronic spectra of a displaced-oscillator model in the condensed phase and discuss the spectral and temporal evolutions of 2D signals. Our results show that the proposed QLE approach is capable of describing vibrational relaxation, decoherence, and vibrational coherence transfer, as well as their manifestations in spectroscopic signals. Furthermore, vibrational quantum beats specific for excited-state vs ground-state nuclear wave packet dynamics can also be identified. We anticipate that this method will provide a useful tool to conduct theoretical studies of 2D spectroscopy for strong vibronically coupled systems and to elucidate intricate vibronic couplings in complex molecular systems.

Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM)

An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.

Exact description of excitonic dynamics in molecular aggregates weakly driven by light

We present a rigorous theoretical description of excitonic dynamics in molecular light-harvesting aggregates photoexcited by weak-intensity radiation of arbitrary properties. While the interaction with light is included up to the second order, the treatment of the excitation-environment coupling is exact and results in an exact expression for the reduced excitonic density matrix that is manifestly related to the spectroscopic picture of the photoexcitation process. This expression takes fully into account the environmental reorganization processes triggered by the two interactions with light. This is particularly important for slow environments and/or strong excitation-environment coupling. Within the exponential decomposition scheme, we demonstrate how our result can be recast as the hierarchy of equations of motion (HEOM) that explicitly and consistently includes the photoexcitation step. We analytically describe the environmental reorganization dynamics triggered by a delta-like excitation of a single chromophore and demonstrate how our HEOM, in appropriate limits, reduces to the Redfield equations comprising a pulsed photoexcitation and the nonequilibrium Förster theory. We also discuss the relation of our formalism to the combined Born-Markov-HEOM approaches in the case of excitation by thermal light.

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.

Linear and nonlinear spectroscopy from quantum master equations

We investigate the accuracy of the second-order time-convolutionless (TCL2) quantum master equation for the calculation of linear and nonlinear spectroscopies of multichromophore systems. We show that even for systems with non-adiabatic coupling, the TCL2 master equation predicts linear absorption spectra that are accurate over an extremely broad range of parameters and well beyond what would be expected based on the perturbative nature of the approach; non-equilibrium population dynamics calculated with TCL2 for identical parameters are significantly less accurate. For third-order (two-dimensional) spectroscopy, the importance of population dynamics and the violation of the so-called quantum regression theorem degrade the accuracy of TCL2 dynamics. To correct these failures, we combine the TCL2 approach with a classical ensemble sampling of slow microscopic bath degrees of freedom, leading to an efficient hybrid quantum-classical scheme that displays excellent accuracy over a wide range of parameters. In the spectroscopic setting, the success of such a hybrid scheme can be understood through its separate treatment of homogeneous and inhomogeneous broadening. Importantly, the presented approach has the computational scaling of TCL2, with the modest addition of an embarrassingly parallel prefactor associated with ensemble sampling. The presented approach can be understood as a generalized inhomogeneous cumulant expansion technique, capable of treating multilevel systems with non-adiabatic dynamics.

Linear and third- and fifth-order nonlinear spectroscopies of a charge transfer system coupled to an underdamped vibration

We study hole, electron, and exciton transports in a charge transfer system in the presence of underdamped vibrational motion. We analyze the signature of these processes in the linear and third-, and fifth-order nonlinear electronic spectra. Calculations are performed with a numerically exact hierarchical equations of motion method for an underdamped Brownian oscillator spectral density. We find that combining electron, hole, and exciton transfers can lead to non-trivial spectra with more structure than with excitonic coupling alone. Traces taken during the waiting time of a two-dimensional (2D) spectrum are dominated by vibrational motion and do not reflect the electron, hole, and exciton dynamics directly. We find that the fifth-order nonlinear response is particularly sensitive to the charge transfer process. While third-order 2D spectroscopy detects the correlation between two coherences, fifth-order 2D spectroscopy (2D population spectroscopy) is here designed to detect correlations between the excited states during two different time periods.

Quantum Mechanical Wave Packet Dynamics at a Conical Intersection with Strong Vibrational Dissipation

We derive a reduced model for the nonadiabatic quantum dynamics of an electronic wave packet moving through a conical intersection in the presence of strong vibrational damping. Starting from the dissipative two-state two-model model, we transform the tuning and the coupling mode to the bath. The resulting quantum two-state model with two highly structured environments is solved numerically exactly in the regime of strong vibrational damping. We find negative cross peaks in the ultrafast optical 2D spectra as clear signatures of the conical intersection. They arise from secondary excitations of the wave packet after having passed through the photophysical energy funnel. This feature is in agreement with recent transient absorption measurements of rhodopsin.

Simulating signatures of two-dimensional electronic spectra of the Fenna-Matthews-Olson complex: By using a numerical path integral

A framework for simulating electronic spectra from photon-echo experiments is constructed by using a numerical path integral technique. This method is non-Markovian and nonperturbative and, more importantly, is not limited by a fixed form of the spectral density functions of the environment. Next, a two-dimensional (2D) third-order electronic spectrum of a dimer system is simulated. The spectrum is in agreement with the experimental and theoretical results previously reported [for example, M. Khalil, N. Demirdöven, and A. Tokmakoff, Phys. Rev. Lett. 90, 047401 (2003)]. Finally, a 2D third-order electronic spectrum of the Fenna-Matthews-Olson (FMO) complex is simulated by using the Debye, Ohmic, and Adolphs and Renger spectral density functions. It is shown that this method can clearly produce the spectral signatures of the FMO complex by using only the Adolphs and Renger spectral density function. Plots of the evolution of the diagonal and cross-peaks show that they are oscillating with the population time.

On the accuracy of coherent modified Redfield theory in simulating excitation energy transfer dynamics

In this study, we investigate the accuracy of a recently developed coherent modified Redfield theory (CMRT) in simulating excitation energy transfer (EET) dynamics. The CMRT is a secular non-Markovian quantum master equation that is derived by extending the modified Redfield theory to treat coherence dynamics in molecular excitonic systems. Herein, we systematically survey the applicability of the CMRT in a large EET parameter space through the comparisons of the CMRT EET dynamics in a dimer system with the numerically exact results. The results confirm that the CMRT exhibits a broad applicable range and allow us to locate the specific parameter regimes where CMRT fails to provide adequate results. Moreover, we propose an accuracy criterion based on the magnitude of second-order perturbation to characterize the applicability of CMRT and show that the criterion summarizes all the benchmark results and the physics described by CMRT. Finally, we employ the accuracy criterion to quantitatively compare the performance of CMRT to that of a small polaron quantum master equation approach. The comparison demonstrates the complementary nature of these two methods, and as a result, the combination of the two methods provides accurate simulations of EET dynamics for the full parameter space investigated in this study. Our results not only delicately evaluate the applicability of the CMRT but also reveal new physical insights for factors controlling the dynamics of EET that should be useful for developing more accurate and efficient methods for simulations of EET dynamics in molecular aggregate systems.

Reduced quantum dynamics with arbitrary bath spectral densities: hierarchical equations of motion based on several different bath decomposition schemes

We investigated applications of the hierarchical equation of motion (HEOM) method to perform high order perturbation calculations of reduced quantum dynamics for a harmonic bath with arbitrary spectral densities. Three different schemes are used to decompose the bath spectral density into analytical forms that are suitable to the HEOM treatment: (1) The multiple Lorentzian mode model that can be obtained by numerically fitting the model spectral density. (2) The combined Debye and oscillatory Debye modes model that can be constructed by fitting the corresponding classical bath correlation function. (3) A new method that uses undamped harmonic oscillator modes explicitly in the HEOM formalism. Methods to extract system-bath correlations were investigated for the above bath decomposition schemes. We also show that HEOM in the undamped harmonic oscillator modes can give detailed information on the partial Wigner transform of the total density operator. Theoretical analysis and numerical simulations of the spin-Boson dynamics and the absorption line shape of molecular dimers show that the HEOM formalism for high order perturbations can serve as an important tool in studying the quantum dissipative dynamics in the intermediate coupling regime.

Simulation of two-dimensional electronic spectra of phycoerythrin 545 at ambient temperature

By using a hierarchical equations-of-motion approach, we reproduce the two-dimensional electronic spectra of phycoerythrin 545 from Rhodomonas CS24 at ambient temperature (294 K). The simulated spectra are in agreement with the experimental results reported in Wong et al. (Nat. Chem. 2012, 4, 396). The evolutions of cross peaks for rephasing spectra and diagonal peaks for nonrephasing spectra have also been plotted. The peaks oscillate with the population times, with frequencies, phases, and amplitudes of the oscillating curves also being qualitatively consistent with the experimental results.

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.

Dynamical steady States in driven quantum systems

We derive dynamical equations for a driven, dissipative quantum system in which the environment-induced relaxation rate is comparable to the Rabi frequency, avoiding assumptions on the frequency dependence of the environmental coupling. When the environmental coupling varies significantly on the scale of the Rabi frequency, secular or rotating wave approximations break down. We avoid these approximations, yielding dynamical steady states which account for the interaction between driven quantum dots and their phonon environment. The theory, which is motivated by recent experimental observations, qualitatively and quantitatively describes the transition from asymmetric unsaturated resonances at weak driving to population inversion at strong driving.

Non-Markovianity: initial correlations and nonlinear optical measurements

By extending the response function approach developed in nonlinear optics, we analytically derive an expression for the non-Markovianity in the time evolution of a system in contact with a quantum mechanical bath, and find a close connection with the directly observable nonlinear optical response. The result indicates that memory in the bath-induced fluctuations rather than in the dissipation causes non-Markovianity. Initial correlations between states of the system and the bath are shown to be essential for a correct understanding of the non-Markovianity. These correlations are included in our treatment through a preparation function.

Optical properties, excitation energy and primary charge transfer in photosystem II: theory meets experiment

In this review we discuss structure-function relationships of the core complex of photosystem II, as uncovered from analysis of optical spectra of the complex and its subunits. Based on descriptions of optical difference spectra including site directed mutagenesis we propose a revision of the multimer model of the symmetrically arranged reaction center pigments, described by an asymmetric exciton Hamiltonian. Evidence is provided for the location of the triplet state, the identity of the primary electron donor, the localization of the cation and the secondary electron transfer pathway in the reaction center. We also discuss the stationary and time-dependent optical properties of the CP43 and CP47 subunits and the excitation energy transfer and trapping-by-charge-transfer kinetics in the core complex.