Excitation Energy Transfer Efficiency in Fluctuating Environments: Role of Quantum Coherence in the Presence of Memory Effects

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.

Memory effects in the efficiency control of energy transfer under incoherent light excitation in noisy environments

Fluctuations in energy gap and coupling constants between chromophores can play an important role in absorption and energy transfer across a collection of two-level systems. In photosynthesis, light-induced quantum coherence can affect the efficiency of energy transfer to the designated "trap" state. Theoretically, the interplay between fluctuations and coherence has been studied often, employing either a Markovian or a perturbative approximation. In this study, we depart from these approaches to incorporate memory effects by using Kubo's quantum stochastic Liouville equation. We introduce the effects of decay of the created excitation (to the ground state) on the desired propagation and trapping that provides a direction of flow of the excitation. In the presence of light-induced pumping, we establish a relation between the efficiency, the mean survival time, and the correlation decay time of the bath-induced fluctuations. A decrease in the steady-state coherence during the transition from the non-Markovian regime to the Markovian limit results in a decrease in efficiency. As in the well-known Haken-Strobl model, the ratio of the square of fluctuation strength to the rate plays a critical role in determining the mechanism of energy transfer and in shaping the characteristics of the efficiency profile. We recover a connection between the transfer flux and the imaginary part of coherences in both equilibrium and excited bath states, in both correlated and uncorrelated bath models. We uncover a non-monotonic dependence of efficiency on site energy heterogeneity for both correlated and uncorrelated bath models.

The influence of a Hamiltonian vibration vs a bath vibration on the 2D electronic spectra of a homodimer

We elucidate the influence of the system–bath boundary placement within an open quantum system, with emphasis on the two-dimensional electronic spectra, through the application of the hierarchical equations of motion formalism for an exciton system. We apply two different models, the Hamiltonian vibration model (HVM) and bath vibration model (BVM), to a monomer and a homodimer. In the HVM, we specifically include the vibronic states in the Hamiltonian capturing vibronic quenching, whereas in the BVM, all vibrational details are contained within the bath and described by an underdamped spectral density. The resultant spectra are analyzed in terms of energetic peak position and thermodynamic broadening precision in order to evaluate the efficacy of the two models. The HVM produces 2D spectra with accurate peak positional information, while the BVM is well suited to modeling dynamic peak broadening. For the monomer, both models produce equivalent spectra in the limit where additional damping associated with the underdamped vibration in the BVM approaches zero. This is supported by analytical results. However, for the homodimer, the BVM spectra are redshifted with respect to the HVM due to an absence of vibronic quenching in the BVM. The computational efficiency of the two models is also discussed in order to inform us of the most appropriate use of each method.