Explicit correlated exciton-vibrational dynamics of the FMO complex

Systematic Coarse Graining of Environments for the Nonperturbative Simulation of Open Quantum Systems

Conducting precise electronic-vibrational dynamics simulations of molecular systems poses significant challenges when dealing with realistic environments composed of numerous vibrational modes. Here, we introduce a technique for the construction of effective phonon spectral densities that capture accurately open-system dynamics over a finite time interval of interest. When combined with existing nonperturbative simulation tools, our approach can reduce significantly the computational costs associated with many-body open-system dynamics.

Manipulation of photosynthetic energy transfer by vibrational strong coupling

Uncovering the mystery of efficient and directional energy transfer in photosynthetic organisms remains a critical challenge in quantum biology. Recent experimental evidence and quantum theory developments indicate the significance of quantum features of molecular vibrations in assisting photosynthetic energy transfer, which provides the possibility of manipulating the process by controlling molecular vibrations. Here, we propose and theoretically demonstrate efficient manipulation of photosynthetic energy transfer by using vibrational strong coupling between the vibrational state of a Fenna-Matthews-Olson (FMO) complex and the vacuum state of an optical cavity. Specifically, based on a full-quantum analytical model to describe the strong coupling effect between the optical cavity and molecular vibration, we realize efficient manipulation of energy transfer efficiency (from 58% to 92%) and energy transfer time (from 20 to 500 ps) in one branch of FMO complex by actively controlling the coupling strength and the quality factor of the optical cavity under both near-resonant and off-resonant conditions, respectively. Our work provides a practical scenario to manipulate photosynthetic energy transfer by externally interfering molecular vibrations via an optical cavity and a comprehensible conceptual framework for researching other similar systems.

Simulating Quantum Vibronic Dynamics at Finite Temperatures With Many Body Wave Functions at 0 K

For complex molecules, nuclear degrees of freedom can act as an environment for the electronic “system” variables, allowing the theory and concepts of open quantum systems to be applied. However, when molecular system-environment interactions are non-perturbative and non-Markovian, numerical simulations of the complete system-environment wave function become necessary. These many body dynamics can be very expensive to simulate, and extracting finite-temperature results—which require running and averaging over many such simulations—becomes especially challenging. Here, we present numerical simulations that exploit a recent theoretical result that allows dissipative environmental effects at finite temperature to be extracted efficiently from a single, zero-temperature wave function simulation. Using numerically exact time-dependent variational matrix product states, we verify that this approach can be applied to vibronic tunneling systems and provide insight into the practical problems lurking behind the elegance of the theory, such as the rapidly growing numerical demands that can appear for high temperatures over the length of computations.

Matrix Product State Simulations of Non-Equilibrium Steady States and Transient Heat Flows in the Two-Bath Spin-Boson Model at Finite Temperatures

Simulating the non-perturbative and non-Markovian dynamics of open quantum systems is a very challenging many body problem, due to the need to evolve both the system and its environments on an equal footing. Tensor network and matrix product states (MPS) have emerged as powerful tools for open system models, but the numerical resources required to treat finite-temperature environments grow extremely rapidly and limit their applications. In this study we use time-dependent variational evolution of MPS to explore the striking theory of Tamascelli et al. (Phys. Rev. Lett. 2019, 123, 090402.) that shows how finite-temperature open dynamics can be obtained from zero temperature, i.e., pure wave function, simulations. Using this approach, we produce a benchmark dataset for the dynamics of the Ohmic spin-boson model across a wide range of coupling strengths and temperatures, and also present a detailed analysis of the numerical costs of simulating non-equilibrium steady states, such as those emerging from the non-perturbative coupling of a qubit to baths at different temperatures. Despite ever-growing resource requirements, we find that converged non-perturbative results can be obtained, and we discuss a number of recent ideas and numerical techniques that should allow wide application of MPS to complex open quantum systems.

Exciton Delocalization in Indolenine Squaraine Aggregates Templated by DNA Holliday Junction Scaffolds

Exciton delocalization plays a prominent role in the photophysics of molecular aggregates, ultimately governing their particular function or application. Deoxyribonucleic acid (DNA) is a compelling scaffold in which to template molecular aggregates and promote exciton delocalization. As individual dye molecules are the basis of exciton delocalization in molecular aggregates, their judicious selection is important. Motivated by their excellent photostability and spectral properties, here, we examine the ability of squaraine dyes to undergo exciton delocalization when aggregated via a DNA Holliday junction (HJ) template. A commercially available indolenine squaraine dye was chosen for the study given its strong structural resemblance to Cy5, a commercially available cyanine dye previously shown to undergo exciton delocalization in DNA HJs. Three types of DNA-dye aggregate configurations-transverse dimer, adjacent dimer, and tetramer-were investigated. Signatures of exciton delocalization were observed in all squaraine-DNA aggregates. Specifically, strong blue shift and Davydov splitting were observed in steady-state absorption spectroscopy and exciton-induced features were evident in circular dichroism (CD) spectroscopy. Strongly suppressed fluorescence emission provided additional, indirect evidence for exciton delocalization in the DNA-templated squaraine dye aggregates. To quantitatively evaluate and directly compare the excitonic Coulombic coupling responsible for exciton delocalization, the strength of excitonic hopping interactions between the dyes was obtained by simultaneously fitting the experimental steady-state absorption and CD spectra via a Holstein-like Hamiltonian, in which, following the theoretical approach of Kühn, Renger, and May, the dominant vibrational mode is explicitly considered. The excitonic hopping strength within indolenine squaraines was found to be comparable to that of the analogous Cy5 DNA-templated aggregate. The squaraine aggregates adopted primarily an H-type (dyes oriented parallel to each other) spatial arrangement. Extracted geometric details of the dye mutual orientation in the aggregates enabled a close comparison of aggregate configurations and the elucidation of the influence of dye angular relationship on excitonic hopping interactions in squaraine aggregates. These results encourage the application of squaraine-based aggregates in next-generation systems driven by molecular excitons.

Toward monitoring the dissipative vibrational energy flows in open quantum systems by mixed quantum-classical simulations

In open quantum system dynamics, rich information about the major energy relaxation channels and corresponding relaxation rates can be elucidated by monitoring the vibrational energy flow among individual bath modes. However, such calculations often become tremendously difficult as the complexity of the subsystem-bath coupling increases. In this paper, we attempt to make this task feasible by using a mixed quantum-classical method, the Poisson-bracket mapping equation with non-Hamiltonian modification (PBME-nH) [H. W. Kim and Y. M. Rhee, J. Chem. Phys. 140, 184106 (2014)]. For a quantum subsystem bilinearly coupled to harmonic bath modes, we derive an expression for the mode energy in terms of the classical positions and momenta of the nuclei, while keeping consistency with the energy of the quantum subsystem. The accuracy of the resulting expression is then benchmarked against a numerically exact method by using relatively simple models. Although our expression predicts a qualitatively correct dissipation rate for a range of situations, cases involving a strong vibronic resonance are quite challenging. This is attributed to the inherent lack of quantum back reaction in PBME-nH, which becomes significant when the subsystem strongly interacts with a small number of bath modes. A rigorous treatment of such an effect will be crucial for developing quantitative simulation methods that can handle generic subsystem-bath coupling.

Intramolecular vibrations enhance the quantum efficiency of excitonic energy transfer

We study the impact of underdamped intramolecular vibrational modes on the efficiency of the excitation energy transfer in a dimer in which each state is coupled to its own underdamped vibrational mode and, in addition, to a continuous background of environmental modes. For this, we use the numerically exact hierarchy equation of motion approach. We determine the quantum yield and the transfer time in dependence of the vibronic coupling strength, and in dependence of the damping of the incoherent background. Moreover, we tune the vibrational frequencies out of resonance with the excitonic energy gap. We show that the quantum yield is enhanced by up to 10% when the vibrational frequency of the donor is larger than at the acceptor. The vibronic energy eigenstates of the acceptor acquire then an increased density of states, which leads to a higher occupation probability of the acceptor in thermal equilibrium. We can conclude that an underdamped vibrational mode which is weakly coupled to the dimer fuels a faster transfer of excitation energy, illustrating that long-lived vibrations can, in principle, enhance energy transfer, without involving long-lived electronic coherence.

Numerical assessment for accuracy and GPU acceleration of TD-DMRG time evolution schemes

The time dependent density matrix renormalization group (TD-DMRG) has become one of the cutting edge methods of quantum dynamics for complex systems. In this paper, we comparatively study the accuracy of three time evolution schemes in the TD-DMRG, the global propagation and compression method with the Runge-Kutta algorithm (P&C-RK), the time dependent variational principle based methods with the matrix unfolding algorithm (TDVP-MU), and with the projector-splitting algorithm (TDVP-PS), by performing benchmarks on the exciton dynamics of the Fenna-Matthews-Olson complex. We show that TDVP-MU and TDVP-PS yield the same result when the time step size is converged and they are more accurate than P&C-RK4, while TDVP-PS tolerates a larger time step size than TDVP-MU. We further adopt the graphical processing units to accelerate the heavy tensor contractions in the TD-DMRG, and it is able to speed up the TDVP-MU and TDVP-PS schemes by up to 73 times.

High-dimensional exciton-vibrational wave-packet dynamics in the FMO complex. influence of site-specific spectral densities

The correlated exciton-vibrational dynamics of the Fenna-Matthews-Olson (FMO) complex is studied using Multi-layer Multi-configuration Time-dependent Hartree (ML-MCTDH) wavepacket propagation. Exciton populations and coherences are shown to be sensitive to the details of the spectral density.

Optimizing co-operative multi-environment dynamics in a dark-state-enhanced photosynthetic heat engine

We analyze the role of coherent, non-perturbative system-bath interactions in a photosynthetic heat engine. Using the reaction-coordinate formalism to describe the vibrational phonon-environment in the engine, we analyze the efficiency around an optimal parameter regime predicted in earlier studies. We show that, in the limit of high-temperature photon irradiation, the phonon-assisted population transfer between bright and dark states is suppressed due to dephasing from the photon environment, even in the Markov limit where we expect the influence of each bath to have an independent and additive effect on the dynamics. Manipulating the phonon bath properties via its spectral density enables us to identify both optimal low- and high-frequency regimes where the suppression can be removed. This suppression of transfer and its removal suggests that it is important to consider carefully the non-perturbative and cooperative effects of system-bath environments in designing artificial photosynthetic systems and also that manipulating inter-environmental interactions could provide a new multidimensional "lever" by which photocells and other types of quantum devices can be optimized.

Finite-temperature time-dependent variation with multiple Davydov states

The Dirac-Frenkel time-dependent variational approach with Davydov Ansätze is a sophisticated, yet efficient technique to obtain an accurate solution to many-body Schrödinger equations for energy and charge transfer dynamics in molecular aggregates and light-harvesting complexes. We extend this variational approach to finite temperature dynamics of the spin-boson model by adopting a Monte Carlo importance sampling method. In order to demonstrate the applicability of this approach, we compare calculated real-time quantum dynamics of the spin-boson model with that from numerically exact iterative quasiadiabatic propagator path integral (QUAPI) technique. The comparison shows that our variational approach with the single Davydov Ansätze is in excellent agreement with the QUAPI method at high temperatures, while the two differ at low temperatures. Accuracy in dynamics calculations employing a multitude of Davydov trial states is found to improve substantially over the single Davydov Ansatz, especially at low temperatures. At a moderate computational cost, our variational approach with the multiple Davydov Ansatz is shown to provide accurate spin-boson dynamics over a wide range of temperatures and bath spectral densities.

Machine learning for quantum dynamics: deep learning of excitation energy transfer properties

Understanding the relationship between the structure of light-harvesting systems and their excitation energy transfer properties is of fundamental importance in many applications including the development of next generation photovoltaics. Natural light harvesting in photosynthesis shows remarkable excitation energy transfer properties, which suggests that pigment-protein complexes could serve as blueprints for the design of nature inspired devices. Mechanistic insights into energy transport dynamics can be gained by leveraging numerically involved propagation schemes such as the hierarchical equations of motion (HEOM). Solving these equations, however, is computationally costly due to the adverse scaling with the number of pigments. Therefore virtual high-throughput screening, which has become a powerful tool in material discovery, is less readily applicable for the search of novel excitonic devices. We propose the use of artificial neural networks to bypass the computational limitations of established techniques for exploring the structure-dynamics relation in excitonic systems. Once trained, our neural networks reduce computational costs by several orders of magnitudes. Our predicted transfer times and transfer efficiencies exhibit similar or even higher accuracies than frequently used approximate methods such as secular Redfield theory.

Effect of Time-Delayed Feedback on the Interaction of a Dimer System with its Environment

In this work, we report modeling of non-Markovian open quantum systems, consisting of an excitonic dimer that displays memory effect due to time delayed interaction with its environment. We, indeed investigate the effect of these time delays on quantum coherence and excitation dynamical behavior in the time domain generally considered for photosynthetic experiments (few hundred femtoseconds). In particular, we show that the coherence is maintained for periods proportional to time delays. Additionally, if delay is taken into account, coupling to the environment can be tuned to lower values, unlike in previous studies. This kind of intriguing effect can, therefore, when generalized to complete systems, permit more control on the experimental parameters, which may lead to more accurate description of the photosynthetic energy transfer functioning and subsequent applications in artificial photovoltaic research.

Simulation of Quantum Dynamics of Excitonic Systems at Finite Temperature: an efficient method based on Thermo Field Dynamics

Quantum electron-vibrational dynamics in molecular systems at finite temperature is described using an approach based on Thermo Field Dynamics theory. This formulation treats temperature effects in the Hilbert space without introducing the Liouville space. The solution of Thermo Field Dynamics equations with a novel technique for the propagation of Tensor Trains (Matrix Product States) is implemented and discussed. The methodology is applied to the study of the exciton dynamics in the Fenna-Mathews-Olsen complex using a realistic structured spectral density to model the electron-phonon interaction. The results of the simulations highlight the effect of specific vibrational modes on the exciton dynamics and energy transfer process, as well as call for careful modeling of electron-phonon couplings.

Multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach to the correlated exciton-vibrational dynamics in the FMO complex

The coupled quantum dynamics of excitonic and vibrational degrees of freedom is investigated for high-dimensional models of the Fenna-Matthews-Olson complex. This includes a seven- and an eight-site model with 518 and 592 harmonic vibrational modes, respectively. The coupling between local electronic transitions and vibrations is described within the Huang-Rhys model using parameters that are obtained by discretization of an experimental spectral density. Different pathways of excitation energy flow are analyzed in terms of the reduced one-exciton density matrix, focussing on the role of vibrational and vibronic excitations. Distinct features due to both competing time scales of vibrational and exciton motion and vibronically assisted transfer are observed. The question of the effect of initial state preparation is addressed by comparing the case of an instantaneous Franck-Condon excitation at a single site with that of a laser field excitation.

Influence of weak vibrational-electronic couplings on 2D electronic spectra and inter-site coherence in weakly coupled photosynthetic complexes

Coherence oscillations measured in two-dimensional (2D) electronic spectra of pigment-protein complexes may have electronic, vibrational, or mixed-character vibronic origins, which depend on the degree of electronic-vibrational mixing. Oscillations from intrapigment vibrations can obscure the inter-site coherence lifetime of interest in elucidating the mechanisms of energy transfer in photosynthetic light-harvesting. Huang-Rhys factors (S) for low-frequency vibrations in Chlorophyll and Bacteriochlorophyll are quite small (S ≤ 0.05), so it is often assumed that these vibrations influence neither 2D spectra nor inter-site coherence dynamics. In this work, we explore the influence of S within this range on the oscillatory signatures in simulated 2D spectra of a pigment heterodimer. To visualize the inter-site coherence dynamics underlying the 2D spectra, we introduce a formalism which we call the "site-probe response." By comparing the calculated 2D spectra with the site-probe response, we show that an on-resonance vibration with Huang-Rhys factor as small as S = 0.005 and the most strongly coupled off-resonance vibrations (S = 0.05) give rise to long-lived, purely vibrational coherences at 77 K. We moreover calculate the correlation between optical pump interactions and subsequent entanglement between sites, as measured by the concurrence. At 77 K, greater long-lived inter-site coherence and entanglement appear with increasing S. This dependence all but vanishes at physiological temperature, as environmentally induced fluctuations destroy the vibronic mixing.

Quantum effects in ultrafast electron transfers within cryptochromes

Cryptochromes and photolyases are flavoproteins that may undergo ultrafast charge separation upon electronic excitation of their flavin cofactors.