Equilibrium excited state and emission spectra of molecular aggregates from the hierarchical equations of motion approach

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.

Application of the imaginary time hierarchical equations of motion method to calculate real time correlation functions

We investigate the application of the imaginary time hierarchical equations of motion method to calculate real time quantum correlation functions. By starting from the path integral expression for the correlated system–bath equilibrium state, we first derive a new set of equations that decouple the imaginary time propagation and the calculation of auxiliary density operators. The new equations, thus, greatly simplify the calculation of the equilibrium correlated initial state that is subsequently used in the real time propagation to obtain the quantum correlation functions. It is also shown that a periodic decomposition of the bath imaginary time correlation function is no longer necessary in the new equations such that different decomposition schemes can be explored. The applicability of the new method is demonstrated in several numerical examples, including the spin-Boson model, the Holstein model, and the double-well model for proton transfer reaction.

Stark absorption and Stark fluorescence spectroscopies: Theory and simulations

Stark spectroscopy experiments are widely used to study the properties of molecular systems, particularly those containing charge-transfer (CT) states. However, due to the small transition dipole moments and large static dipole moments of the CT states, the standard interpretation of the Stark absorption and Stark fluorescence spectra in terms of the Liptay model may be inadequate. In this work, we provide a theoretical framework for calculations of Stark absorption and Stark fluorescence spectra and propose new methods of simulations that are based on the quantum-classical theory. In particular, we use the forward-backward trajectory solution and a variant of the Poisson bracket mapping equation, which have been recently adapted for the calculation of conventional (field-free) absorption and fluorescence spectra. For comparison, we also apply the recently proposed complex time-dependent Redfield theory, while exact results are obtained using the hierarchical equations of motion approach. We show that the quantum-classical methods produce accurate results for a wide range of systems, including those containing CT states. The CT states contribute significantly to the Stark spectra, and the standard Liptay formalism is shown to be inapplicable for the analysis of spectroscopic data in those cases. We demonstrate that states with large static dipole moments may cause a pronounced change in the total fluorescence yield of the system in the presence of an external electric field. This effect is correctly captured by the quantum-classical methods, which should therefore prove useful for further studies of Stark spectra of real molecular systems. As an example, we calculate the Stark spectra for the Fenna-Matthews-Olson complex of green sulfur bacteria.

Quantum-Classical Approach for Calculations of Absorption and Fluorescence: Principles and Applications

Absorption and fluorescence spectroscopy techniques provide a wealth of information on molecular systems. The simulations of such experiments remain challenging, however, despite the efforts put into developing the underlying theory. An attractive method of simulating the behavior of molecular systems is provided by the quantum-classical theory─it enables one to keep track of the state of the bath explicitly, which is needed for accurate calculations of fluorescence spectra. Unfortunately, until now there have been relatively few works that apply quantum-classical methods for modeling spectroscopic data. In this work, we seek to provide a framework for the calculations of absorption and fluorescence lineshapes of molecular systems using the methods based on the quantum-classical Liouville equation, namely, the forward-backward trajectory solution (FBTS) and the non-Hamiltonian variant of the Poisson bracket mapping equation (PBME-nH). We perform calculations on a molecular dimer and the photosynthetic Fenna-Matthews-Olson complex. We find that in the case of absorption, the FBTS outperforms PBME-nH, consistently yielding highly accurate results. We next demonstrate that for fluorescence calculations, the method of choice is a hybrid approach, which we call PBME-nH-Jeff, that utilizes the effective coupling theory [Gelzinis, A.; J. Chem. Phys. 2020, 152, 051103] to estimate the excited state equilibrium density operator. Thus, we find that FBTS and PBME-nH-Jeff are excellent candidates for simulating, respectively, absorption and fluorescence spectra of real molecular systems.

Strong Exciton-Vibrational Coupling in Molecular Assemblies. Dynamics Using the Polaron Transformation in HEOM Space

In Frenkel exciton dynamics of aggregated molecules, the polaron transformation (PT) technique leads to decoupling of diagonal elements in the subspace of excited electronic states from vibrations. In this article we describe for the first time how PT becomes applicable in the framework of the "Hierarchical Equations of Motion" (HEOM) approach for treatment of open quantum systems. We extend the concept of formulating operators in HEOM space by deriving hierarchical equations of PT which lead to a shift in the excited state potential energy surface to compensate its displacement. While the assumption of thermal equilibration of the vibrational oscillators, introduced by PT, results in a stationary state in a monomer, in a dimer under the same assumption nonequilibrium dynamics appears because of the interplay of the transfer process and vibrational equilibration. Both vertical transitions generating a vibrationally hot state and initially equilibrated vibrational oscillators evolve toward the same stationary asymptotic state associated with polaron formation. The effect of PT on the dynamics of this process depends on initial excitation and basis representation of the electronic system. The developed approach facilitates a generic formulation of quantum master equations involving perturbative treatment of polaron dynamics.

An extension of stochastic hierarchy equations of motion for the equilibrium correlation functions

A traditional stochastic hierarchy equations of motion method is extended into the correlated real-time and imaginary-time propagations, in this paper, for its applications in calculating the equilibrium correlation functions. The central idea is based on a combined employment of stochastic unravelling and hierarchical techniques for the temperature-dependent and temperature-free parts of the influence functional, respectively, in the path integral formalism of the open quantum systems coupled to a harmonic bath. The feasibility and validity of the proposed method are justified in the emission spectra of homodimer compared to those obtained through the deterministic hierarchy equations of motion. Besides, it is interesting to find that the complex noises generated from a small portion of real-time and imaginary-time cross terms can be safely dropped to produce the stable and accurate position and flux correlation functions in a broad parameter regime.

Construction of Multichromophoric Spectra from Monomer Data: Applications to Resonant Energy Transfer

We develop a model that establishes a quantitative link between the physical properties of molecular aggregates and their constituent building blocks. The relation is built on the coherent potential approximation, calibrated against exact results, and proven reliable for a wide range of parameters. It provides a practical method to compute spectra and transfer rates in multichromophoric systems from experimentally accessible monomer data. Applications to Förster energy transfer reveal optimal transfer rates as functions of both the system-bath coupling and intra-aggregate coherence.

Effect of Pulse Shaping on Observing Coherent Energy Transfer in Single Light-Harvesting Complexes

Recent experimental and theoretical studies have revealed that quantum coherence plays an important role in the excitation energy transfer in photosynthetic light-harvesting (LH) complexes. Inspired by the recent single-molecule two-color double-pump experiment, we theoretically investigate the effect of pulse shaping on observing coherent energy transfer in the single bacterial LH2 complex. It is found that quantum coherent energy transfer can be observed when the time delay and phase difference between the two laser pulses are controlled independently. However, when the two-color pulses are generated using the pulse-shaping method, how the laser pulses are prepared is crucial to the observation of quantum coherent energy transfer in single photosynthetic complexes.

Simulated two-dimensional electronic spectroscopy of the eight-bacteriochlorophyll FMO complex

The Fenna-Matthews-Olson (FMO) protein-pigment complex acts as a molecular wire conducting energy between the outer antenna system and the reaction center; it is an important photosynthetic system to study the transfer of excitonic energy. Recent crystallographic studies report the existence of an additional (eighth) bacteriochlorophyll a (BChl a) in some of the FMO monomers. To understand the functionality of this eighth BChl, we simulated the two-dimensional electronic spectra of both the 7-site (apo form) and the 8-site (holo form) variant of the FMO complex from green sulfur bacteria, Prosthecochloris aestuarii. By comparing the spectrum, it was found that the eighth BChl can affect two different excitonic energy transfer pathways: (1) it is directly involved in the first apo form pathway (6 → 3 → 1) by passing the excitonic energy to exciton 6; and (2) it facilitates an increase in the excitonic wave function overlap between excitons 4 and 5 in the second pathway (7 → 4,5 → 2 → 1) and thus increases the possible downward sampling routes across the BChls.

Förster resonance energy transfer, absorption and emission spectra in multichromophoric systems. I. Full cumulant expansions and system-bath entanglement

We study the Förster resonant energy transfer rate, absorption and emission spectra in multichromophoric systems. The multichromophoric Förster theory (MCFT) is determined from an overlap integral of generalized matrices related to the donor's emission and acceptor's absorption spectra, which are obtained via a full 2nd-order cumulant expansion technique developed in this work. We calculate the spectra and MCFT rate for both localized and delocalized systems, and calibrate the analytical 2nd-order cumulant expansion with the exact stochastic path integral method. We present three essential findings: (i) The role of the initial entanglement between the donor and its bath is found to be crucial in both the emission spectrum and the MCFT rate. (ii) The absorption spectra obtained by the cumulant expansion method are nearly identical to the exact spectra for both localized and delocalized systems, even when the system-bath coupling is far from the perturbative regime. (iii) For the emission spectra, the cumulant expansion can give reliable results for localized systems, but fail to provide reliable spectra of the high-lying excited states of a delocalized system, when the system-bath coupling is large and the thermal energy is small. This paper also provides a simple golden-rule derivation of the MCFT, reviews existing methods, and motivates further developments in the subsequent papers.

Förster resonance energy transfer, absorption and emission spectra in multichromophoric systems. III. Exact stochastic path integral evaluation

A numerically exact path integral treatment of the absorption and emission spectra of open quantum systems is presented that requires only the straightforward solution of a stochastic differential equation. The approach converges rapidly enabling the calculation of spectra of large excitonic systems across the complete range of system parameters and for arbitrary bath spectral densities. With the numerically exact absorption and emission operators, one can also immediately compute energy transfer rates using the multi-chromophoric Förster resonant energy transfer formalism. Benchmark calculations on the emission spectra of two level systems are presented demonstrating the efficacy of the stochastic approach. This is followed by calculations of the energy transfer rates between two weakly coupled dimer systems as a function of temperature and system-bath coupling strength. It is shown that the recently developed hybrid cumulant expansion (see Paper II) is the only perturbative method capable of generating uniformly reliable energy transfer rates and emission spectra across a broad range of system parameters.

Förster resonance energy transfer, absorption and emission spectra in multichromophoric systems. II. Hybrid cumulant expansion

We develop a hybrid cumulant expansion method to account for the system-bath entanglement in the emission spectrum in the multi-chromophoric Förster transfer rate. In traditional perturbative treatments, the emission spectrum is usually expanded with respect to the system-bath coupling term in both real and imaginary time. This perturbative treatment gives a reliable absorption spectrum, where the bath is Gaussian and only the real-time expansion is involved. For the emission spectrum, the initial state is an entangled state of the system plus bath. Traditional perturbative methods are problematic when the excitations are delocalized and the energy gap is larger than the thermal energy, since the second-order expansion cannot predict the displacement of the bath. In the present method, the real-time dynamics is carried out by using the 2nd-order cumulant expansion method, while the displacement of the bath is treated more accurately by utilizing the exact reduced density matrix of the system. In a sense, the hybrid cumulant expansion is based on a generalized version of linear response theory with entangled initial states.