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

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.

Simulating optical linear absorption for mesoscale molecular aggregates: An adaptive hierarchy of pure states approach

In this paper, we present dyadic adaptive HOPS (DadHOPS), a new method for calculating linear absorption spectra for large molecular aggregates. This method combines the adaptive HOPS (adHOPS) framework, which uses locality to improve computational scaling, with the dyadic HOPS method previously developed to calculate linear and nonlinear spectroscopic signals. To construct a local representation of dyadic HOPS, we introduce an initial state decomposition that reconstructs the linear absorption spectra from a sum over locally excited initial conditions. We demonstrate the sum over initial conditions can be efficiently Monte Carlo sampled and that the corresponding calculations achieve size-invariant [i.e., O(1)] scaling for sufficiently large aggregates while trivially incorporating static disorder in the Hamiltonian. We present calculations on the photosystem I core complex to explore the behavior of the initial state decomposition in complex molecular aggregates as well as proof-of-concept DadHOPS calculations on an artificial molecular aggregate inspired by perylene bis-imide to demonstrate the size-invariance of the method.

Tensor Network Path Integral Study of Dynamics in B850 LH2 Ring with Atomistically Derived Vibrations

The recently introduced multisite tensor network path integral (MS-TNPI) allows simulation of extended quantum systems coupled to dissipative media. We use MS-TNPI to simulate the exciton transport and the absorption spectrum of a B850 bacteriochlorophyll (BChl) ring. The MS-TNPI network is extended to account for the ring topology of the B850 system. Accurate molecular-dynamics-based description of the molecular vibrations and the protein scaffold is incorporated through the framework of Feynman-Vernon influence functional. To relate the present work with the excitonic picture, an exploration of the absorption spectrum is done by simulating it using approximate and topologically consistent transition dipole moment vectors. Comparison of these numerically exact MS-TNPI absorption spectra are shown with second-order cumulant approximations. The effect of temperature on both the exact and the approximate spectra is also explored.

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. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.