Simulation of the two-dimensional electronic spectra of the Fenna-Matthews-Olson complex using the hierarchical equations of motion method

Transient Chiral Dynamics in the Fenna-Matthews-Olson Complex Revealed by Two-Dimensional Circular Dichroism Spectroscopy

Chirality plays a pivotal role across scientific disciplines with profound implications spanning light-matter interactions, molecular recognition, and natural evolutionary processes. This study delves into the active influence of molecular chirality on exciton energy transfer within photosynthetic protein complexes, focusing on the Fenna-Matthews-Olson (FMO) complex. Employing two-dimensional circular dichroism (2DCD) spectroscopy, we investigate the transient chiral dynamics of excitons during energy transfer processes within the FMO complex. Our approach, incorporating pulse information into population dynamics based on the third-order response function, facilitates the calculation of 2DCD spectra and dynamics. This enables the extraction of chiral contributions to excitonic energy transfer and the examination of electronic wave functions. We demonstrate that 2DCD spectra offer excitation energies that are better resolved than those from conventional two-dimensional electronic spectroscopy. These findings deepen our understanding of exciton energy transfer mechanisms in natural photosynthesis, emphasizing the potential of 2DCD spectroscopy as a powerful tool for unraveling the chiral contribution to exciton dynamics.

Efficient formulation of multitime generalized quantum master equations: Taming the cost of simulating 2D spectra

Modern 4-wave mixing spectroscopies are expensive to obtain experimentally and computationally. In certain cases, the unfavorable scaling of quantum dynamics problems can be improved using a generalized quantum master equation (GQME) approach. However, the inclusion of multiple (light-matter) interactions complicates the equation of motion and leads to seemingly unavoidable cubic scaling in time. In this paper, we present a formulation that greatly simplifies and reduces the computational cost of previous work that extended the GQME framework to treat arbitrary numbers of quantum measurements. Specifically, we remove the time derivatives of quantum correlation functions from the modified Mori-Nakajima-Zwanzig framework by switching to a discrete-convolution implementation inspired by the transfer tensor approach. We then demonstrate the method's capabilities by simulating 2D electronic spectra for the excitation-energy-transfer dimer model. In our method, the resolution of data can be arbitrarily coarsened, especially along the t2 axis, which mirrors how the data are obtained experimentally. Even in a modest case, this demands O(103) fewer data points. We are further able to decompose the spectra into one-, two-, and three-time correlations, showing how and when the system enters a Markovian regime where further measurements are unnecessary to predict future spectra and the scaling becomes quadratic. This offers the ability to generate long-time spectra using only short-time data, enabling access to timescales previously beyond the reach of standard methodologies.

Transient chiral dynamics revealed by two-dimensional circular dichroism spectroscopy

Chirality has been considered as one of the key factors in the evolution of life in nature. It is important to uncover how chiral potentials of molecular systems play vital role in fundamental photochemical processes. Here, we investigate the role of chirality in photoinduced energy transfer in a model dimeric system, where the monomers are excitonically coupled. To observe transient chiral dynamics and energy transfer, we employ circularly polarized laser pulses in two-dimensional electronic spectroscopy to construct the two-dimensional circular dichroism (2DCD) spectral maps. Tracking time-resolved peak magnitudes in 2DCD spectra allows one to identify chirality induced population dynamics. The dynamics of energy transfer is revealed by the time-resolved kinetics of cross peaks. However, the differential signal of 2DCD spectra shows the magnitude of cross peaks is dramatically reduced at initial waiting time, which indicates the weak chiral interactions between two monomers. The downhill energy transfer is resolved by presenting a strong magnitude of cross peak in 2DCD spectra after long waiting time. The chiral contribution towards coherent and incoherent energy-transfer pathways in the model dimer system is further examined via control of excitonic couplings between two monomers. Applications are made to study the energy-transfer process in the Fenna-Matthews-Olson complex. Our work uncovers the potential of 2DCD spectroscopy to resolve the chiral-induced interactions and population transfers in excitonically coupled systems.

Coherent excitation energy transfer in model photosynthetic reaction center: Effects of non-Markovian quantum environment

Excitation energy transfer (EET) and electron transfer (ET) are crucially involved in photosynthetic processes. In reality, the photosynthetic reaction center constitutes an open quantum system of EET and ET, which manifests an interplay of pigments, solar light and phonon baths. So far theoretical studies have been mainly based on master equation approaches in the Markovian condition. The non-Markovian environmental effect, which may play a crucial role, has not been sufficiently considered. In this work, we propose a mixed dynamic approach to investigate this open system. The influence of phonon bath is treated via the exact dissipaton equation of motion (DEOM) while that of photon bath is via the Lindblad master equation. Specifically, we explore the effect of non-Markovian quantum phonon bath on the coherent transfer dynamics and its manipulation on the current--voltage behavior. Distinguished from the results of completely Markovian Lindblad equation and those adopting classical environment description, the mixed DEOM--Lindblad simulations exhibittransfer coherence up to a few hundreds femtosecondsand the related environmental manipulation effect on current.These non-Markovian quantum coherent effects may be extended tomore complex and realistic systems and be helpful to thedesign of organic photovoltaic devices.

Theoretical Study of the Spectral Differences of the Fenna-Matthews-Olson Protein from Different Species and Their Mutants

The structural basis for the spectral differences between the Fenna-Matthews-Olson (FMO) proteins from Chlorobaculum tepidum (C. tepidum) and Prosthecochloris aestuarii 2K (P. aestuarii) is yet to be fully understood. Mutation-induced perturbation to the exciton structure and the optical spectra of the complex provide a suitable means to investigate the critical role played by the protein scaffold. In this work, we have performed quantum-mechanics/molecular-mechanics calculations over the molecular dynamics simulation trajectories with the polarized protein-specific charge scheme for both wild-type FMOs and two mutants. Our result reveals that a single-point mutation in the vicinity of BChl 6, namely, W183F of C. tepidum, significantly affects the absorption spectrum, resulting in a switch of the absorption spectrum from type 2, for which the 806 nm band is more pronounced than the 815 nm band, to type 1, for which the 815 nm band is pronounced. Our observations agree with the single-point mutation experiments reported by Saer et al. (Biochim. Biophys. Acta, Bioenerg. 2017, 1858, 288-296) and Khmelnitskiy et al. (J. Phys. Chem. Lett. 2018, 9, 3378-3386). In contrast, the absorption spectrum of the P. aestuarii experiences the opposite transition (from type 1 to type 2) upon the same mutation. Furthermore, by comparing the contributions of individual pigments to the spectra in the wild type and its mutant, we find that a single-point mutation near BChl 6 not only induces changes in excitation energy of BChl 6 per se but also affects the excitonic structures of the neighboring BChls 5 and 7 through strong interpigment electronic couplings, resulting in a significant change in the absorption spectra.

Correlated vibration-solvent effects on the non-Condon exciton spectroscopy

Excitation energy transfer is crucially involved in a variety of systems. During the process, the non-Condon vibronic coupling and the surrounding solvent interaction may synergetically play important roles. In this work, we study the correlated vibration-solvent influences on the non-Condon exciton spectroscopy. Statistical analysis is elaborated for the overall vibration-plus-solvent environmental effects. Analytic solutions are derived for the linear absorption of monomer systems. General simulations are accurately carried out via the dissipaton-equation-of-motion approach. The resulted spectra in either the linear absorption or strong field regime clearly demonstrate the coherence enhancement due to the synergetic vibration-solvent correlation.

Efficient propagation of the hierarchical equations of motion using the Tucker and hierarchical Tucker tensors

We develop new methods to efficiently propagate the hierarchical equations of motion (HEOM) by using the Tucker and hierarchical Tucker (HT) tensors to represent the reduced density operator and auxiliary density operators. We first show that by employing the split operator method, the specific structure of the HEOM allows a simple propagation scheme using the Tucker tensor. When the number of effective modes in the HEOM increases and the Tucker representation becomes intractable, the split operator method is extended to the binary tree structure of the HT representation. It is found that to update the binary tree nodes related to a specific effective mode, we only need to propagate a short matrix product state constructed from these nodes. Numerical results show that by further employing the mode combination technique commonly used in the multi-configuration time-dependent Hartree approaches, the binary tree representation can be applied to study excitation energy transfer dynamics in a fairly large system including over 10⁴ effective modes. The new methods may thus provide a promising tool in simulating quantum dynamics in condensed phases.

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.

Effect of disorder and polarization sequences on two-dimensional spectra of light-harvesting complexes

Two-dimensional electronic spectra (2DES) provide unique ways to track the energy transfer dynamics in light-harvesting complexes. The interpretation of the peaks and structures found in experimentally recorded 2DES is often not straightforward, since several processes are imaged simultaneously. The choice of specific pulse polarization sequences helps to disentangle the sometimes convoluted spectra, but brings along other disturbances. We show by detailed theoretical calculations how 2DES of the Fenna-Matthews-Olson complex are affected by rotational and conformational disorder of the chromophores.

Hierarchical equations of motion method based on Fano spectrum decomposition for low temperature environments

The hierarchical equations of motion (HEOM) method has become one of the most popular methods for the studies of the open quantum system. However, its applicability to systems at ultra-low temperatures is largely restrained by the enormous computational cost, which is caused by the numerous exponential functions required to accurately characterize the non-Markovian memory of the reservoir environment. To overcome this problem, a Fano spectrum decomposition (FSD) scheme has been proposed recently [Cui et al., J. Chem. Phys. 151, 024110 (2019)], which expands the reservoir correlation functions using polynomial-exponential functions and hence greatly reduces the size of the memory basis set. In this work, we explicitly establish the FSD-based HEOM formalisms for both bosonic and fermionic environments. The accuracy and efficiency of the FSD-based HEOM are exemplified by the calculated low-temperature dissipative dynamics of a spin-boson model and the dynamic and static properties of a single-orbital Anderson impurity model in the Kondo regime. The encouraging numerical results highlight the practicality and usefulness of the FSD-based HEOM method for general open systems at ultra-low temperatures.

Does electronic coherence enhance anticorrelated pigment vibrations under realistic conditions?

The light-harvesting efficiency of a photoactive molecular complex is largely determined by the properties of its electronic quantum states. Those, in turn, are influenced by molecular vibrational states of the nuclear degrees of freedom. Here, we reexamine two recently formulated concepts that a coherent vibronic coupling between molecular states would either extend the electronic coherence lifetime or enhance the amplitude of the anticorrelated vibrational mode at longer times. For this, we study a vibronically coupled dimer and calculate the nonlinear two-dimensional (2D) electronic spectra that directly reveal electronic coherence. The time scale of electronic coherence is initially extracted by measuring the antidiagonal bandwidth of the central peak in the 2D spectrum at zero waiting time. Based on the residual analysis, we identify small-amplitude long-lived oscillations in the cross-peaks, which, however, are solely due to groundstate vibrational coherence, regardless of having resonant or off-resonant conditions. Our studies neither show an enhancement of the electronic quantum coherence nor an enhancement of the anticorrelated vibrational mode by the vibronic coupling under ambient conditions.

Two-dimensional electronic vibrational spectroscopy and ultrafast excitonic and vibronic photosynthetic energy transfer

2-Dimensional electronic vibrational spectroscopy presents a novel experimental and theoretical approach to study energy transfer.

On the exact truncation tier of fermionic hierarchical equations of motion

The hierarchical equations of motion (HEOM) theory is in principle exact for describing the dissipative dynamics of quantum systems linearly coupled to Gaussian environments. In practice, the hierarchy needs to be truncated at a finite tier. We demonstrate that, for general systems described by the fermionic HEOM, the (n+L̃)th-tier truncation with L̃=2N_σN_ν yields the exact density operators up to the nth tier. Here, N_σ = 2 for fermionic systems and N_ν is the system degrees of freedom. For noninteracting systems, L̃ is further reduced by half. Such an exact termination pattern originates from the Pauli exclusion principle for fermions, and it holds true regardless of the system-environment coupling strength, the number of coupling reservoirs, or the specific scheme employed to unravel the environment memory contents. The relatively small L̃ emphasizes the nonperturbative nature of the HEOM theory. We also propose a simplified HEOM approach to further reduce the memory cost for practical calculations.

Two-dimensional photon echoes reveal non-Markovian energy transfer in an excitonic dimer

We show that strong non-Markovian effects can be revealed by the steady-state two-dimensional (2D) photon echo spectra at asymptotic waiting times. For this, we use a simple dimer toy model that is strongly coupled to a harmonic bath with parameters typical for photoactive biomolecules. We calculate the 2D photon echo spectra employing both the numerically exact hierarchy equation of motion and the quasiadiabatic path integral approach and compare these results with approximate results from a time-nonlocal quantum master equation approach. While the latter correctly reproduces the exact population dynamics at long times, it fails at the same time to correctly describe the 2D photon echo spectra at long waiting times. The differences show that non-Markovian effects are much more important for the steady-state 2D photon echoes than for the equilibrium populations. Thus, accurate theoretical descriptions of the energy transfer dynamics in biomolecular complexes have to be based on numerically exact simulations of the environmental fluctuations when nonlinear response functions are analyzed.

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.

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.

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.