Tight inner ring architecture and quantum motion of nuclei enable efficient energy transfer in bacterial light harvesting

Discrete Generalized Quantum Master Equations

Several derivative and integral approximations are explored for discretizing the Nakajima-Zwanzig generalized quantum master equation (NZ-QME or GQME) to obtain discrete quantum master equation (DQME) hierarchies and relations between discrete memory kernel and reduced density matrix (RDM) elements. It is shown that the simplest forward-difference approximation does not allow the reliable determination of the discrete kernel elements, even in the infinitesimal time-step limit, and that discrete kernels obtained in earlier work are flawed, although the procedure can be remedied. The various approximations give rise to DQMEs that differ in structure and in the RDM-kernel relationships. It is shown that the use of a more accurate discretization based on the midpoint derivative and midpoint integral approximations leads to a DQME that exhibits endpoint effects, which reflect the weaker impact of the bath on the RDM during the first time step and which parallel those encountered in the small matrix decomposition of the path integral (SMatPI) with a symmetric factorization of the short-time propagator. The features of the DQME hierarchies and RDM-kernel relations are illustrated through analytical examples involving a simple integrodifferential equation and a scalar GQME model, as well as numerical results for a two-level system (TLS) coupled to a harmonic bath.

Quantum algorithm for the simulation of non-Markovian quantum dynamics using Feynman-Vernon influence functional

In this work, we developed a quantum algorithm for the simulation of non-Markovian quantum dynamics based on the Feynman-Vernon's path integral formulation. The algorithm performs the full path sum and proves to be polynomial either in time or in space, compared to the same classical algorithm, which is exponential in time. In addition, the algorithm has no classical overhead and is equally applicable regardless of whether the temporal entanglement due to non-Markovianity is low or high, making it a unified framework for simulating non-Markovian dynamics in open quantum systems.

Investigation of quantum trajectories in photosynthetic light harvesting through a quantum stochastic approach

In natural photosynthesis systems, pigment-protein complexes harvest the photon from sunlight with near-unity quantum efficiency. These complexes show incredible properties that cannot be merely extrapolated from knowledge of their composition. Additionally, the environment perturbing the light-harvesting process significantly affects the mechanism of photosynthesis. This research investigates the photosystem II reaction center (PSII RC) from a new perspective which considers the restricted path of the exciton transfer, in the photosynthesis system, as a quantum trajectory picture with the quantum continuous measurement. In this work, the corridor path of exciton transfer dynamics satisfies the equation of motion, as the spin dynamics, which consists of precession, relaxation, and random force rapidly fluctuating spin splitting arising from the bath. Moreover, the width of the corridor is an important factor for restricting path dynamics resulting in the localization and decoherence phenomenon. Our method is to analyze exciton transfer dynamics through paths on the Bloch sphere, in order to investigate the propagating states in accordance with the weight functional which depends on the coupling parameter between the system and environment as the phonon bath. Our results show that the paths outside the width of the corridor have a considerably lower weight functional and decoherence functional than those inside the width. Therefore, the degrees of localization, the weight functional, and the decoherence functional are related. Furthermore, the simulation reveals three characteristics of exciton transfer: gradual transfer, no transfer, and rapid transfer, relying significantly on the coupling between the system and phonons.

Variational Quantum Algorithm for Non-Markovian Quantum Dynamics Using an Ensemble of Ehrenfest Trajectories

The simulation of non-Markovian quantum dynamics plays an important role in the understanding of charge and exciton dynamics in the condensed phase environment, yet such a simulation remains computationally expensive on classical computers. In this work, we develop a variational quantum algorithm that is capable of simulating non-Markovian quantum dynamics on quantum computers. The algorithm captures the non-Markovian effect by employing the Ehrenfest trajectories and Monte Carlo sampling of their thermal distribution. We test the algorithm with the spin-boson model on the quantum simulator, and the results match quantitatively with the exact ones. The algorithm naturally fits into the parallel computing platform of the NISQ devices and can be extended to anharmonic system-bath interactions and multistate systems.

Coherence in Chemistry: Foundations and Frontiers

Coherence refers to correlations in waves. Because matter has a wave-particle nature, it is unsurprising that coherence has deep connections with the most contemporary issues in chemistry research (e.g., energy harvesting, femtosecond spectroscopy, molecular qubits and more). But what does the word "coherence" really mean in the context of molecules and other quantum systems? We provide a review of key concepts, definitions, and methodologies, surrounding coherence phenomena in chemistry, and we describe how the terms "coherence" and "quantum coherence" refer to many different phenomena in chemistry. Moreover, we show how these notions are related to the concept of an interference pattern. Coherence phenomena are indeed complex, and ambiguous definitions may spawn confusion. By describing the many definitions and contexts for coherence in the molecular sciences, we aim to enhance understanding and communication in this broad and active area of chemistry.

Semiclassical dynamics in Wigner phase space II: Nonadiabatic hybrid Wigner dynamics

We present an approximate semiclassical (SC) framework for mixed quantized dynamics in Wigner phase space in a two-part series. In the first article, we introduced the Adiabatic Hybrid Wigner Dynamics (AHWD) method that allows for a few important "system" degrees of freedom to be quantized using high-level double Herman-Kluk SC theory while describing the rest (the "bath") using classical-limit linearized SC theory. In this second article, we extend our hybrid Wigner dynamics to nonadiabatic processes. The resulting Nonadiabatic Hybrid Wigner Dynamics (NHWD) has two variants that differ in the choice of degrees of freedom to be quantized. Specifically, we introduce NHWD(E) where only the electronic state variables are quantized and the NHWD(V) where both electronic state variables and a handful of strongly coupled nuclear modes are quantized. We show that while NHWD(E) proves accurate for a wide range of scattering models and spin-boson models, systems where a few nuclear modes are strongly coupled to electronic states require NHWD(V) to accurately capture the long-time dynamics. Taken together, we show that AHWD and NHWD represent a new framework for SC simulations of high-dimensional systems with significant quantum effects.

Improved memory truncation scheme for quasi-adiabatic propagator path integral via influence functional renormalization

Accurately simulating non-Markovian quantum dynamics in system-bath coupled problems remains challenging. In this work, we present a novel memory truncation scheme for the iterative quasi-adiabatic propagator path integral (iQuAPI) method to improve accuracy. Conventional memory truncation in iQuAPI discards all influence functional beyond a certain time interval, which is not effective for problems with a long memory time. Our proposed scheme selectively retains the most significant parts of the influence functional using the density matrix renormalization group algorithm. We validate the effectiveness of our scheme through simulations of the spin-boson model across various parameter sets, demonstrating faster convergence and improved accuracy compared to the conventional scheme. Our findings suggest that the new memory truncation scheme significantly advances the capabilities of iQuAPI for problems with a long memory time.

ACE: A general-purpose non-Markovian open quantum systems simulation toolkit based on process tensors

We describe a general-purpose computational toolkit for simulating open quantum systems, which provides numerically exact solutions for composites of zero-dimensional quantum systems that may be strongly coupled to multiple, quite general non-Markovian environments. It is based on process tensor matrix product operators (PT-MPOs), which efficiently encapsulate environment influences. The code features implementations of several PT-MPO algorithms, in particular Automated Compression of Environments for general environments comprised of independent modes as well as schemes for generalized spin boson models. The latter includes a divide-and-conquer scheme for periodic PT-MPOs, which enable million time step simulations for realistic models. PT-MPOs can be precalculated and reused for efficiently probing different time-dependent system Hamiltonians. They can also be stacked together and combined to provide numerically complete solutions of small networks of open quantum systems. The code is written in C++ and is fully controllable by configuration files, for which we have developed a versatile and compact human-readable format.

An ensemble variational quantum algorithm for non-Markovian quantum dynamics

Many physical and chemical processes in a condensed phase environment exhibit non-Markovian quantum dynamics. As such simulations are challenging on classical computers, we developed a variational quantum algorithm that is capable of simulating non-Markovian dynamics on noisy intermediate-scale quantum (NISQ) devices. We used a quantum system linearly coupled to its harmonic bath as the model Hamiltonian. The non-Markovianity is captured by introducing auxiliary variables from the bath trajectories. With Monte Carlo sampling of the bath degrees of freedom, finite temperature dynamics is produced. We validated the algorithm on a simulator and demonstrated its performance on an IBM quantum device. The framework developed naturally adapts to any anharmonic bath with non-linear coupling to the system, and is also well suited for simulating spin chain dynamics in a dissipative environment.

Examining Sound, Light, and Vibrations as Tools to Manage Microbes and Support Holobionts, Ecosystems, and Technologies

The vast array of interconnected microorganisms across Earth's ecosystems and within holobionts has been called the "Internet of Microbes." Bacteria and archaea are masters of energy and information collection, storage, transformation, and dissemination using both "wired" and wireless (at a distance) functions. Specific tools affecting microbial energy and information functions offer effective strategies for managing microbial populations within, between, and beyond holobionts. This narrative review focuses on microbial management using a subset of physical modifiers of microbes: sound and light (as well as related vibrations). These are examined as follows: (1) as tools for managing microbial populations, (2) as tools to support new technologies, (3) as tools for healing humans and other holobionts, and (4) as potential safety dangers for microbial populations and their holobionts. Given microbial sensitivity to sound, light, and vibrations, it is critical that we assign a higher priority to the effects of these physical factors on microbial populations and microbe-laden holobionts. We conclude that specific sound, light, and/or vibrational conditions are significant therapeutic tools that can help support useful microbial populations and help to address the ongoing challenges of holobiont disease. We also caution that inappropriate sound, light, and/or vibration exposure can represent significant hazards that require greater recognition.

Quantum information scrambling and chemical reactions

The ultimate regularity of quantum mechanics creates a tension with the assumption of classical chaos used in many of our pictures of chemical reaction dynamics. Out-of-time-order correlators (OTOCs) provide a quantum analog to the Lyapunov exponents that characterize classical chaotic motion. Maldacena, Shenker, and Stanford have suggested a fundamental quantum bound for the rate of information scrambling, which resembles a limit suggested by Herzfeld for chemical reaction rates. Here, we use OTOCs to study model reactions based on a double-well reaction coordinate coupled to anharmonic oscillators or to a continuum oscillator bath. Upon cooling, as one enters the tunneling regime where the reaction rate does not strongly depend on temperature, the quantum Lyapunov exponent can approach the scrambling bound and the effective reaction rate obtained from a population correlation function can approach the Herzfeld limit on reaction rates: Tunneling increases scrambling by expanding the state space available to the system. The coupling of a dissipative continuum bath to the reaction coordinate reduces the scrambling rate obtained from the early-time OTOC, thus making the scrambling bound harder to reach, in the same way that friction is known to lower the temperature at which thermally activated barrier crossing goes over to the low-temperature activationless tunneling regime. Thus, chemical reactions entering the tunneling regime can be information scramblers as powerful as the black holes to which the quantum Lyapunov exponent bound has usually been applied.

Kink Sum for Long-Memory Small Matrix Path Integral Dynamics

The small matrix decomposition of the real-time path integral (SMatPI) allows for numerically exact and efficient propagation of the reduced density matrix (RDM) for system-bath Hamiltonians. Its high efficiency lies in the small size of the SMatPI matrices employed in the iterative algorithm, whose size is equal to that of the full RDM. By avoiding the storage and multiplication of large tensors, the SMatPI algorithm is applicable in multistate systems under a variety of conditions. The main computational effort is the evaluation of path sums within the entangled memory length to construct the SMatPI matrices. A number of methods are available for this task, each with its own favorable parameter regime, but calculations with strong system-bath coupling and long memory at low temperatures remain out of reach. The present paper evaluates the path sums by binning the paths (in forward time only) based on their amplitudes, which depend on the number and type of kinks they contain. The algorithm is very efficient, leading to a dramatic acceleration of path sums and significantly extending the accessible memory length in the most challenging regimes.

Exact Non-Markovian Quantum Dynamics on the NISQ Device Using Kraus Operators

The theory of open quantum systems has many applications ranging from simulating quantum dynamics in condensed phases to better understanding quantum-enabled technologies. At the center of theoretical chemistry are the developments of methodologies and computational tools for simulating charge and excitation energy transfer in solutions, biomolecules, and molecular aggregates. As a variety of these processes display non-Markovian behavior, classical computer simulation can be challenging due to exponential scaling with existing methods. With quantum computers holding the promise of efficient quantum simulations, in this paper, we present a new quantum algorithm based on Kraus operators that capture the exact non-Markovian effect at a finite temperature. The implementation of the Kraus operators on the quantum machine uses a combination of singular value decomposition (SVD) and optimal Walsh operators that result in shallow circuits. We demonstrate the feasibility of the algorithm by simulating the spin-boson dynamics and the exciton transfer in the Fenna-Matthews-Olson (FMO) complex. The NISQ results show very good agreement with the exact ones.

Impact of Spatial Inhomogeneity on Excitation Energy Transport in the Fenna-Matthews-Olson Complex

The dynamics of the excitation energy transfer (EET) in photosynthetic complexes is an interesting question both from the perspective of fundamental understanding and the research in artificial photosynthesis. Over the past decade, very accurate spectral densities have been developed to capture spatial inhomogeneities in the Fenna-Matthews-Olson (FMO) complex. However, challenges persist in numerically simulating these systems, both in terms of parameterizing them and following their dynamics over long periods of time because of long non-Markovian memories. We investigate the dynamics of FMO with the exact treatment of various theoretical spectral densities using the new tensor network path integral-based methods, which are uniquely capable of addressing the difficulty of long memory length and incoherent Förster theory. It is also important to be able to analyze the pathway of EET flow, which can be difficult to identify given the non-trivial structure of connections between bacteriochlorophyll molecules in FMO. We use the recently introduced ideas of relating coherence to population derivatives to analyze the transport process and reveal some new routes of transport. The combination of exact and approximate methods sheds light on the role of coherences in affecting the fine details of the transport and promises to be a powerful toolbox for future exploration of other open systems with quantum transport.

Coherence Maps and Flow of Excitation Energy in the Bacterial Light Harvesting Complex 2

We present and analyze coherence maps [ J. Phys. Chem. B 2022, 126, 9361-9375] to investigate the quantum coherences that are created, sustained, and damped by vibrational modes during the transfer of excitation energy from the B800 (outer) to the B850 (inner) ring of the light harvesting complex 2 (LH2) of purple bacteria with a variety of initial conditions. The reduced density matrix of the 24-pigment complex, where the ground and excited electronic states of each bacteriochlorophyll are explicitly coupled to 50 intramolecular vibrations at room temperature, is obtained from fully quantum-mechanical small matrix path integral (SMatPI) calculations. The coherence maps show a very rapid localization within the outer ring, accompanied by the formation of inter-ring quantum superpositions that evolve to a partial quantum delocalization at equilibrium, and quantify in state-to-state detail the flow of energy within the complex.