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.

PathSum: A C++ and Fortran suite of fully quantum mechanical real-time path integral methods for (multi-)system + bath dynamics.

This paper reports the release of PathSum, a new software suite of state-of-the-art path integral methods for studying the dynamics of single or extended systems coupled to harmonic environments. The package includes two modules, suitable for system-bath problems and extended systems comprising many coupled system-bath units, and is offered in C++ and Fortran implementations. The system-bath module offers the recently developed small matrix path integral (SMatPI) and the well-established iterative quasi-adiabatic propagator path integral (i-QuAPI) method for iteration of the reduced density matrix of the system. In the SMatPI module, the dynamics within the entanglement interval can be computed using QuAPI, the blip sum, time evolving matrix product operators, or the quantum-classical path integral method. These methods have distinct convergence characteristics and their combination allows a user to access a variety of regimes. The extended system module provides the user with two algorithms of the modular path integral method, applicable to quantum spin chains or excitonic molecular aggregates. An overview of the methods and code structure is provided, along with guidance on method selection and representative examples.

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.

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

The efficient, directional transfer of absorbed solar energy between photosynthetic light-harvesting complexes continues to pose intriguing questions. In this work, we identify the pathways of energy flow between the B800 and B850 rings in the LH2 complex of Rhodopseudomonas molischianum using fully quantum mechanical path integral methods to simulate the excited-state dynamics of the 24 bacteriochlorophyll molecules and their coupling to 50 normal mode vibrations in each chromophore. While all pigments are identical, the tighter packing of the inner B850 ring is responsible for the thermodynamic stabilization of the inner ring. Molecular vibrations enable the 1-ps flow of energy to the B850 states, which would otherwise be kinetically inaccessible. A classical treatment of the vibrations leads to uniform equilibrium distribution of the excitation, with only 67% transferred to the inner ring. However, spontaneous fluctuations associated with the quantum motion of the nuclei increase the transfer efficiency to 90%.

Interference between Franck-Condon and Herzberg-Teller Terms in the Condensed-Phase Molecular Spectra of Metal-Based Tetrapyrrole Derivatives.

The commonly used Franck-Condon (FC) approximation is inadequate for explaining the electronic spectra of compounds that possess vibrations with substantial Herzberg-Teller (HT) couplings. Metal-based tetrapyrrole derivatives, which are ubiquitous natural pigments, often exhibit prominent HT activity. In this paper, we compare the condensed phase spectra of zinc-tetraphenylporphyrin (ZnTPP) and zinc-phthalocyanine (ZnPc), which exhibit vastly different spectral features in spite of sharing a common tetrapyrrole backbone. The absorption and emission spectra of ZnTPP are characterized by a lack of mirror symmetry and nontrivial temperature dependence. In contrast, mirror symmetry is restored, and the nontrivial temperature-dependent features disappear in ZnPc. We attribute these differences to FC-HT interference, which is less pronounced in ZnPc because of a larger FC component in the dipole moment that leads to FC-dominated transitions. A single minimalistic FC-HT vibronic model reproduces all the experimental spectral features of these molecules. These observations suggest that FC-HT interference is highly susceptible to chemical modification.

B800-to-B850 relaxation of excitation energy in bacterial light harvesting: All-state, all-mode path integral simulations.

We report fully quantum mechanical simulations of excitation energy transfer within the peripheral light harvesting complex (LH2) of Rhodopseudomonas molischianum at room temperature. The exciton-vibration Hamiltonian comprises the 16 singly excited bacteriochlorophyll states of the B850 (inner) ring and the 8 states of the B800 (outer) ring with all available electronic couplings. The electronic states of each chromophore couple to 50 intramolecular vibrational modes with spectroscopically determined Huang-Rhys factors and to a weakly dissipative bath that models the biomolecular environment. Simulations of the excitation energy transfer following photoexcitation of various electronic eigenstates are performed using the numerically exact small matrix decomposition of the quasiadiabatic propagator path integral. We find that the energy relaxation process in the 24-state system is highly nontrivial. When the photoexcited state comprises primarily B800 pigments, a rapid intra-band redistribution of the energy sharply transitions to a significantly slower relaxation component that transfers 90% of the excitation energy to the B850 ring. The mixed character B850* state lacks the slow component and equilibrates very rapidly, providing an alternative energy transfer channel. This (and also another partially mixed) state has an anomalously large equilibrium population, suggesting a shift to lower energy by virtue of exciton-vibration coupling. The spread of the vibrationally dressed states is smaller than that of the eigenstates of the bare electronic Hamiltonian. The total population of the B800 band is found to decay exponentially with a 1/e time of 0.5 ps, which is in good agreement with experimental results.

Small Matrix Quantum-Classical Path Integral.

The quantum-classical path integral (QCPI) is a rigorous formulation of nonadiabatic dynamics, where the dynamical interaction between a quantum system and its environment is captured consistently through classical trajectories driven by forces along quantum paths of the system. In this Letter, we develop a small matrix decomposition (SMatQCPI) that eliminates the tensor storage requirements of the iterative QCPI algorithm. In the case of a system coupled to a harmonic bath, SMatQCPI provides fully quantum mechanical propagation, which also reduces the computational cost to that of a single QCPI step. Further, the SMatQCPI matrices only need to account for quantum contributions to decoherence, allowing high efficiency in challenging regimes of incoherent dynamics. Overall, this new composite algorithm combines the best features of two powerful path integral formulations and offers a versatile tool for simulating condensed phase quantum dynamics.

Franck-Condon and Herzberg-Teller Signatures in Molecular Absorption and Emission Spectra.

Some molecules of chemical and biological significance possess vibrations with significant Herzberg-Teller (HT) couplings, which render the Franck-Condon (FC) approximation inadequate and cause the breakdown of the well-known mirror-image symmetry between linear absorption and emission spectra. Using a model two-state system with displaced harmonic potential surfaces, we show analytically that the FC-HT interference gives rise to asymmetric intensity modification, which has the same sign for all transitions on one side of the 0-0 absorption line and the opposite sign in the equivalent fluorescence transitions, while the trend is exactly reversed for all transitions on the other side the 0-0 line. We examine the dependence of the absorption-emission asymmetry on the mode frequency, Huang-Rhys factor, and dipole moment parameters to show the recovery of symmetry with particular combinations of parameters and a crossover from fluorescence to absorption dominance. We illustrate the analytical predictions through numerically exact calculations in models of one and two discrete vibrational modes and in the presence of a harmonic dissipative bath. In addition to homogeneous broadening effects, we identify large asymmetric shifts of absorption and emission band maxima, which can produce the illusion of unequal frequencies in the ground and excited potential surfaces as well as a nontrivial modulation of spectral asymmetry by temperature, which results from the enhancement of transitions on one side of the 0-0 line. These findings will aid the interpretation of experimental spectra in HT-active molecular systems.

Synthetic Control of Exciton Dynamics in Bioinspired Cofacial Porphyrin Dimers.

Understanding how the complex interplay among excitonic interactions, vibronic couplings, and reorganization energy determines coherence-enabled transport mechanisms is a grand challenge with both foundational implications and potential payoffs for energy science. We use a combined experimental and theoretical approach to show how a modest change in structure may be used to modify the exciton delocalization, tune electronic and vibrational coherences, and alter the mechanism of exciton transfer in covalently linked cofacial Zn-porphyrin dimers (meso-beta linked ABm-ß and meso-meso linked AAm-m). While both ABm-ß and AAm-m feature zinc porphyrins linked by a 1,2-phenylene bridge, differences in the interporphyrin connectivity set the lateral shift between macrocycles, reducing electronic coupling in ABm-ß and resulting in a localized exciton. Pump-probe experiments show that the exciton dynamics is faster by almost an order of magnitude in the strongly coupled AAm-m dimer, and two-dimensional electronic spectroscopy (2DES) identifies a vibronic coherence that is absent in ABm-ß. Theoretical studies indicate how the interchromophore interactions in these structures, and their system-bath couplings, influence excitonic delocalization and vibronic coherence-enabled rapid exciton transport dynamics. Real-time path integral calculations reproduce the exciton transfer kinetics observed experimentally and find that the linking-modulated exciton delocalization strongly enhances the contribution of vibronic coherences to the exciton transfer mechanism, and that this coherence accelerates the exciton transfer dynamics. These benchmark molecular design, 2DES, and theoretical studies provide a foundation for directed explorations of nonclassical effects on exciton dynamics in multiporphyrin assemblies.

Intramolecular Vibrations in Excitation Energy Transfer: Insights from Real-Time Path Integral Calculations.

Excitation energy transfer (EET) is fundamental to many processes in chemical and biological systems and carries significant implications for the design of materials suitable for efficient solar energy harvest and transport. This review discusses the role of intramolecular vibrations on the dynamics of EET in nonbonded molecular aggregates of bacteriochlorophyll, a perylene bisimide, and a model system, based on insights obtained from fully quantum mechanical real-time path integral results for a Frenkel exciton Hamiltonian that includes all vibrational modes of each molecular unit at finite temperature. Generic trends, as well as features specific to the vibrational characteristics of the molecules, are identified. Weak exciton-vibration (EV) interaction leads to compact, near-Gaussian densities on each electronic state, whose peak follows primarily a classical trajectory on a torus, while noncompact densities and nonlinear peak evolution are observed with strong EV coupling. Interaction with many intramolecular modes and increasing aggregate size smear, shift, and damp these dynamical features.

Time Evolution of Bath Properties in Spin-Boson Dynamics.

The dynamical behaviors of a two-level system (TLS) coupled to a harmonic dissipative bath has been studied extensively using a variety of analytical and numerical methods. The focus of the vast majority of these studies has been on the properties of the TLS, averaged with respect to the bath degrees of freedom. In this work, we use real-time path integral methods to probe the behavior of select bath degrees of freedom during the dynamics of a symmetric two-level system (TLS) coupled to a dissipative bath by calculating system-bath densities (SBD) and coordinate expectation values. Overall, the SBD motion on each diabatic state is simpler than the motion of the total density. In the weak coupling regime, which characterizes the parameters of oscillators that comprise such a bath, the SBD on each TLS state remains primarily compact and Gaussian-like, such that its peak is well characterized by the mode expectation value. In the absence of a dissipative environment, nonadiabatic density depletion leads to spikes in coordinate expectation values. The evolution of the SBD peak trajectory for two discrete modes exhibits Lissajous patterns with frequency-dependent shapes that strongly resemble classical trajectory motion on a torus. These patterns become more complex when the coupling of the mode to the TLS is increased outside of this regime, leading to persistent small amplitude oscillations in the TLS populations characterized by a very slow decay and SBD trajectories that exhibit behaviors reminiscent of chaotic classical systems. Indirect coupling to a dissipative bath has a stabilizing effect on the dynamics, eliminating spikes, synchronizing the SBD motion on the two diabatic states and regularizing the SBD trajectory to simple rectangular Lissajous-like shapes with a slowly shrinking boundary, regardless of the mode frequencies.

Electronic-vibrational density evolution in a perylene bisimide dimer: mechanistic insights into excitation energy transfer.

The process of excitation energy transfer (EET) in molecular aggregates is etched with the signatures of a multitude of electronic and vibrational time scales that often are extremely difficult to resolve. The effect of the motion associated with one molecular vibration on that of another is fundamental to the dynamics of EET. In this paper we present simple theoretical ideas along with fully quantum mechanical calculations to develop a comprehensive mechanistic picture of EET in terms of the time evolution of electronic-vibrational densities (EVD) in a perylene bisimide (PBI) dimer, where 28 intramolecular normal modes couple to the ground and excited electronic states of each molecule. The EVD motion exhibits a plethora of dynamical features, which impart physical justification for the composite effects observed in the EET dynamics. Weakly coupled vibrations lead to classical-like motion of the EVD center on each electronic state, while highly nontrivial EVD characteristics develop under moderate or strong exciton-vibration interaction, leading to the formation of split or crescent-shaped densities, as well as density retention that slows down energy transfer and creates new peaks in the electronic populations. Pronounced correlation effects are observed in two-mode projections of the EVD, as a consequence of indirect vibrational coupling between uncoupled normal modes induced by the electronic coupling. Such indirect coupling depends on the strength of exciton-vibration interactions as well as the frequency mismatch between the two modes and leaves nontrivial signatures in the electronic population dynamics. The collective effects of many vibrational modes cause a partial smearing of these features through dephasing.

Origin of vibrational features in the excitation energy transfer dynamics of perylene bisimide J-aggregates.

We investigate the role of intramolecular normal mode vibrations in the excitation energy transfer (EET) dynamics of perylene bisimide J-aggregates composed of 2 or 25 units using numerically exact methods. The calculations employ a Frenkel exciton Hamiltonian where the ground and excited electronic states of each molecular unit are coupled to 28 intramolecular normal mode vibrations at various temperatures. The electronic populations exhibit strong damping effects, a lengthening of the EET time scale, and complex dynamical patterns, which depend on aggregate length, temperature, as well as electronic and vibrational initial conditions and which are not additive. The early evolution is dominated by high-frequency vibrational modes, but all modes are responsible for the observed dynamics after the initial 25 fs. Overall, we observe significant changes in the electronic populations upon varying the temperature between 0 and 600 K. With a Franck-Condon (FC) initial excitation, a strongly coupled vibrational mode introduces new peaks to the dimer populations, which show very weak temperature sensitivity. The first of these peaks is also seen in the long aggregate, but subsequent recurrences appear strongly quenched and merged. These structures are drastically altered if a non-FC initial condition is assumed. Additional insights are obtained from the diagonal elements of the dimer electronic-vibrational reduced density matrix. We find that the vibronic peaks result from depletion of the crossing region during the early coherent evolution of the vibrational density away from the crossing point, which allows the premature back-transfer of excitation to the initially excited unit.

Real-Time Path Integral Simulation of Exciton-Vibration Dynamics in Light-Harvesting Bacteriochlorophyll Aggregates.

The mechanism of excitation energy transfer in photoexcited bacteriochlorophyll (BChl) aggregates poses intriguing questions, which have important implications for the observed efficiency of photosynthesis. We investigate this process through fully quantum mechanical calculations of exciton-vibration dynamics in chains and rings of BChl a molecules, with parameters characterizing the B850 ring of the LH2 complex of photosynthetic bacteria. The calculations are performed using the modular path integral methodology, which allows the exact treatment of 50 intramolecular vibrations on each pigment using parameters obtained from spectroscopic Huang-Rhys factors with computational effort that scales linearly with aggregate length. Our results indicate that the interplay between electronic and vibrational time scales leads to the rapid suppression but not the overdamping of electronic coherence, which facilitates the spreading of excitation energy throughout the aggregate.

Modular path integral for finite-temperature dynamics of extended systems with intramolecular vibrations.

The modular decomposition of the path integral is a linear-scaling, numerically exact algorithm for calculating dynamical properties of extended systems composed of multilevel units with local couplings. In a recent article, we generalized the method to wavefunction propagation in aggregates characterized by non-diagonal couplings between adjacent units. Here, we extend the method to the calculation of reduced density matrices in aggregates where each unit includes an arbitrary number of coupled harmonic bath modes, which may describe intramolecular normal mode vibrations, at finite temperature. The effects of harmonic modes are included through influence functional factors, which involve analytical expressions that we derive. Representative applications to spin arrays described by the Heisenberg Hamiltonian with dissipative interactions and to J-aggregates of perylene bisimide, where all coupled normal modes are treated explicitly, are presented.

Modular path integral for discrete systems with non-diagonal couplings.

The modular decomposition of the path integral, which leads to linear scaling with the system length, is extended to Hamiltonians with intermonomer couplings that are not diagonalizable in any single-particle basis. An optimal factorization of the time evolution operator is identified, which minimizes the number of path integral variables while ensuring high accuracy and preservation of detailed balance. The modular path integral decomposition is described, along with a highly efficient tensor factorization of the path linking process. The algorithm is illustrated with applications to a model of coupled spins and a Frenkel exciton chain.