Exciton dynamics in chromophore aggregates with correlated environment fluctuations

Simulating Chemistry on Bosonic Quantum Devices

Bosonic quantum devices offer a novel approach to realize quantum computations, where the quantum two-level system (qubit) is replaced with the quantum (an)harmonic oscillator (qumode) as the fundamental building block of the quantum simulator. The simulation of chemical structure and dynamics can then be achieved by representing or mapping the system Hamiltonians in terms of bosonic operators. In this Perspective, we review recent progress and future potential of using bosonic quantum devices for addressing a wide range of challenging chemical problems, including the calculation of molecular vibronic spectra, the simulation of gas-phase and solution-phase adiabatic and nonadiabatic chemical dynamics, the efficient solution of molecular graph theory problems, and the calculations of electronic structure.

Pursuing excitonic energy transfer with programmable DNA-based optical breadboards

DNA nanotechnology has now enabled the self-assembly of almost any prescribed 3-dimensional nanoscale structure in large numbers and with high fidelity. These structures are also amenable to site-specific modification with a variety of small molecules ranging from drugs to reporter dyes. Beyond obvious application in biotechnology, such DNA structures are being pursued as programmable nanoscale optical breadboards where multiple different/identical fluorophores can be positioned with sub-nanometer resolution in a manner designed to allow them to engage in multistep excitonic energy-transfer (ET) via Förster resonance energy transfer (FRET) or other related processes. Not only is the ability to create such complex optical structures unique, more importantly, the ability to rapidly redesign and prototype almost all structural and optical analogues in a massively parallel format allows for deep insight into the underlying photophysical processes. Dynamic DNA structures further provide the unparalleled capability to reconfigure a DNA scaffold on the fly in situ and thus switch between ET pathways within a given assembly, actively change its properties, and even repeatedly toggle between two states such as on/off. Here, we review progress in developing these composite materials for potential applications that include artificial light harvesting, smart sensors, nanoactuators, optical barcoding, bioprobes, cryptography, computing, charge conversion, and theranostics to even new forms of optical data storage. Along with an introduction into the DNA scaffolding itself, the diverse fluorophores utilized in these structures, their incorporation chemistry, and the photophysical processes they are designed to exploit, we highlight the evolution of DNA architectures implemented in the pursuit of increased transfer efficiency and the key lessons about ET learned from each iteration. We also focus on recent and growing efforts to exploit DNA as a scaffold for assembling molecular dye aggregates that host delocalized excitons as a test bed for creating excitonic circuits and accessing other quantum-like optical phenomena. We conclude with an outlook on what is still required to transition these materials from a research pursuit to application specific prototypes and beyond.

Mapping Molecular Hamiltonians into Hamiltonians of Modular cQED Processors

We introduce a general method based on the operators of the Dyson-Masleev transformation to map the Hamiltonian of an arbitrary model system into the Hamiltonian of a circuit Quantum Electrodynamics (cQED) processor. Furthermore, we introduce a modular approach to programming a cQED processor with components corresponding to the mapping Hamiltonian. The method is illustrated as applied to quantum dynamics simulations of the Fenna-Matthews-Olson (FMO) complex and the spin-boson model of charge transfer. Beyond applications to molecular Hamiltonians, the mapping provides a general approach to implement any unitary operator in terms of a sequence of unitary transformations corresponding to powers of creation and annihilation operators of a single bosonic mode in a cQED processor.

Determining interchromophore effects for energy transport in molecular networks using machine-learning algorithms

Nature uses chromophore networks, with highly optimized structural and energetic characteristics, to perform important chemical functions. Due to its modularity, predictable aggregation characteristics, and established synthetic protocols, structural DNA nanotechnology is a promising medium for arranging chromophore networks with analogous structural and energetic controls. However, this high level of control creates a greater need to know how to optimize the systems precisely. This study uses the system's modularity to produce variations of a coupled 14-Site chromophore network. It uses machine-learning algorithms and spectroscopy measurements to reveal the energy-transport roles of these Sites, paying particular attention to the cooperative and inhibitive effects they impose on each other for transport across the network. The physical significance of these patterns is contextualized, using molecular dynamics simulations and energy-transport modeling. This analysis yields insights about how energy transfers across the Donor-Relay and Relay-Acceptor interfaces, as well as the energy-transport pathways through the homogeneous Relay segment. Overall, this report establishes an approach that uses machine-learning methods to understand, in fine detail, the role that each Site plays in an optoelectronic molecular network.

Correlated spectral fluctuations quantified by line shape analysis of fifth-order two-dimensional electronic spectra

Correlated spectral fluctuations were suggested to coordinate excitation transport inside natural light harvesting complexes. We demonstrate the capacities of 2D line shapes from fifth-order coherent electronic signals (R5-2D) to report on such fluctuations in molecular aggregates and present a stochastic approach to fluctuations in correlated site and bi-exciton binding energies in the optical dynamics of Frenkel excitons. The model is applied to R5-2D line shapes of a homodimer, and we show that the peak tilt dynamics are a measure for site energy disorder, inter-site correlation, and the strength of bi-exciton binding energy fluctuations.

Influence of Hydrophobicity on Excitonic Coupling in DNA-Templated Indolenine Squaraine Dye Aggregates

Control over the strength of excitonic coupling in molecular dye aggregates is a substantial factor for the development of technologies such as light harvesting, optoelectronics, and quantum computing. According to the molecular exciton model, the strength of excitonic coupling is inversely proportional to the distance between dyes. Covalent DNA templating was proved to be a versatile tool to control dye spacing on a subnanometer scale. To further expand our ability to control photophysical properties of excitons, here, we investigated the influence of dye hydrophobicity on the strength of excitonic coupling in squaraine aggregates covalently templated by DNA Holliday Junction (DNA HJ). Indolenine squaraines were chosen for their excellent spectral properties, stability, and diversity of chemical modifications. Six squaraines of varying hydrophobicity from highly hydrophobic to highly hydrophilic were assembled in two dimer configurations and a tetramer. In general, the examined squaraines demonstrated a propensity toward face-to-face aggregation behavior observed via steady-state absorption, fluorescence, and circular dichroism spectroscopies. Modeling based on the Kühn-Renger-May approach quantified the strength of excitonic coupling in the squaraine aggregates. The strength of excitonic coupling strongly correlated with squaraine hydrophobic region. Dimer aggregates of dichloroindolenine squaraine were found to exhibit the strongest coupling strength of 132 meV (1065 cm-1). In addition, we identified the sites for dye attachment in the DNA HJ that promote the closest spacing between the dyes in their dimers. The extracted aggregate geometries, and the role of electrostatic and steric effects in squaraine aggregation are also discussed. Taken together, these findings provide a deeper insight into how dye structures influence excitonic coupling in dye aggregates covalently templated via DNA, and guidance in design rules for exciton-based materials and devices.

Computational elucidations on the role of vibrations in energy transfer processes of photosynthetic complexes

Coupling between pigment excitations and nuclear movements in photosynthetic complexes is known to modulate the excitation energy transfer (EET) efficiencies. Toward providing microscopic information, researchers often apply simulation techniques and investigate how vibrations are involved in EET processes. Here, reports on such roles of nuclear movements are discussed from a theory perspective. While vibrations naturally present random thermal fluctuations that can affect energy transferring characteristics, they can also be intertwined with exciton structures and create more specific non-adiabatic energy transfer pathways. For reliable simulations, a bath model that accurately mimics a given molecular system is required. Methods for obtaining such a model in combination with quantum chemical electronic structure calculations and molecular dynamics trajectory simulations are discussed. Various quantum dynamics simulation tools that can handle pigment-to-pigment energy transfers together with their vibrational characters are also touched on. Behaviors of molecular vibrations often deviate from ideality, especially when all-atom details are included, which practically forces us to treat them classically. We conclude this perspective by considering some recent reports that suggest that classical descriptions of bath effects with all-atom details may still produce valuable information for analyzing sophisticated contributions by vibrations to EET processes.

Understanding Disorder, Vibronic Structure, and Delocalization in Electronically Coupled Dimers on DNA Duplexes

Structural DNA nanotechnology is a promising approach to create chromophore networks with modular structures and Hamiltonians to control the material's functions. The functional behaviors of these systems depend on the interactions of the chromophores' vibronic states, as well as interactions with their environment. To optimize their functions, it is necessary to characterize the chromophore network's structural and energetic properties, including the electronic delocalization in some cases. In this study, parameters of interest are deduced in DNA-scaffolded Cyanine 3 and Cyanine 5 dimers. The methods include steady-state optical measurements, physical modeling, and a genetic algorithm approach. The parameters include the chromophore network's vibronic Hamiltonian, molecular positions, transition dipole orientations, and environmentally induced energy broadening. Additionally, the study uses temperature-dependent optical measurements to characterize the spectral broadening further. These combined results reveal the quantum mechanical delocalization, which is important for functions like coherent energy transport and quantum information applications.

Role of Spatially Correlated Fluctuations in Photosynthetic Excitation Energy Transfer with an Equilibrium and a Nonequilibrium Initial Bath

The transfer of excitation energy in photosynthetic light-harvesting complexes has inspired growing interest for its scientific and engineering significance. Recent experimental findings have suggested that spatially correlated environmental fluctuations may account for the existence of long-lived quantum coherent energy transfer observed even at physiological temperature. In this paper, we investigate the effects of spatial correlations on the excitation energy transfer dynamics by including a nonequilibrium initial bath in a simulated donor-acceptor model. The initial bath state, which is assumed to be either equilibrium or nonequilibrium, is expanded in powers of coupling strength within the polaron formalism of a quantum master equation. The spatial correlations of bath fluctuations strongly influence the decay of coherence in the dynamics. The role of a nonequilibrium initial bath is also influenced by spatial correlations and becomes the most conspicuous for certain degrees of spatial correlations from which we propose a picture that the spatial correlations of bath fluctuations open up new energy transfer pathways, playing a role of protecting coherence. Besides, we apply the polaron master equation approach to study the dynamics in a two-site subsystem of the FMO complex and provide a practical example that shows the versatility of this approach.

Stochastically Realized Observables for Excitonic Molecular Aggregates

We show that a stochastic approach enables calculations of the optical properties of large 2-dimensional and nanotubular excitonic molecular aggregates. Previous studies of such systems relied on numerically diagonalizing the dense and disordered Frenkel Hamiltonian, which scales approximately as O(N3) for N dye molecules. Our approach scales much more efficiently as O(Nlog(N)), enabling quick study of systems with a million of coupled molecules on the micrometer size scale. We calculate several important experimental observables, including the optical absorption spectrum and density of states, and develop a stochastic formalism for the participation ratio. Quantitative agreement with traditional matrix diagonalization methods is demonstrated for both small- and intermediate-size systems. The stochastic methodology enables the study of the effects of spatial-correlation in site energies on the optical signatures of large 2D aggregates. Our results demonstrate that stochastic methods present a path forward for screening structural parameters and validating experiments and theoretical predictions in large excitonic aggregates.

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.

Discrimination of Diverse Coherences Allows Identification of Electronic Transitions of a Molecular Nanoring

The role of quantum coherence in photochemical functions of molecular systems such as photosynthetic complexes is a broadly debated topic. Coexistence and intermixing of electronic and vibrational coherences has been proposed to be responsible for the observed long-lived coherences and high energy transfer efficiency. However, clear experimental evidence of coherences with different origins operating at the same time has been elusive. In this work, multidimensional spectra obtained from a six-porphyrin nanoring system are analyzed in detail with support from theoretical modeling. We uncover a great diversity of separable electronic, vibrational, and mixed coherences and show their cooperation in shaping the spectroscopic response. The results permit direct assignment of electronic and vibronic states and characterization of the excitation dynamics. The clear disentanglement of coherences in molecules with extended π-conjugation opens up new avenues for exploring coherent phenomena and understanding their importance for the function of complex systems.

Light Harvesting for Organic Photovoltaics

The field of organic photovoltaics has developed rapidly over the last 2 decades, and small solar cells with power conversion efficiencies of 13% have been demonstrated. Light absorbed in the organic layers forms tightly bound excitons that are split into free electrons and holes using heterojunctions of electron donor and acceptor materials, which are then extracted at electrodes to give useful electrical power. This review gives a concise description of the fundamental processes in photovoltaic devices, with the main emphasis on the characterization of energy transfer and its role in dictating device architecture, including multilayer planar heterojunctions, and on the factors that impact free carrier generation from dissociated excitons. We briefly discuss harvesting of triplet excitons, which now attracts substantial interest when used in conjunction with singlet fission. Finally, we introduce the techniques used by researchers for characterization and engineering of bulk heterojunctions to realize large photocurrents, and examine the formed morphology in three prototypical blends.

Quantum Chemical Studies of Light Harvesting

The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.

Disorder-Induced Quantum Beats in Two-Dimensional Spectra of Excitonically Coupled Molecules

Quantum superposition of molecular electronic states is very fragile because of thermal energy fluctuations and the static conformational disorder induced by the intimate surrounding of constituent molecules of the system. However, the nature of the long-lived quantum beats, observed in time-resolved spectra of molecular aggregates at physiological conditions, is still being debated. We present our study of the conditions when long-lived electronic quantum coherences originating from recently proposed inhomogeneous broadening mechanism are enhanced and reflected in the two-dimensional electronic spectra of the excitonically coupled molecular dimer. We show that depending on the amount of inhomogeneous broadening, the excitonically coupled molecular system can establish long-lived electronic coherences, caused by a disordered subensemble, for which the dephasing due to static energy disorder becomes significantly reduced. On the basis of these considerations, we present explanations for why the electronic or vibrational coherences were or were not observed in a range of recent experiments.

Role of coherent vibrations in energy transfer and conversion in photosynthetic pigment-protein complexes

Oscillatory features of two-dimensional spectra of photosynthetic pigment-protein complexes during few picoseconds after electronic excitations of chlorophylls in various pigment-proteins were recently related to the coherent nuclear vibrations. It has been also speculated that the vibrations may assist the excitonic energy transfer and charge separation, hence contributing to energy transport and energy conversion efficiency. Here, we consider three theoretical approaches usually used for characterization of the excitation dynamics and charge separation, namely Redfield, Förster, and Marcus model descriptions, regarding this question. We show that two out of the three mechanisms require explicit resonances of excitonic splittings and the nuclear vibration frequencies. However, the third one related to the electron transfer is in principle off resonant.

Efficient tool to calculate two-dimensional optical spectra for photoactive molecular complexes

We combine the coherent modified Redfield theory (CMRT) with the equation of motion-phase matching approach (PMA) to calculate two-dimensional photon-echo spectra for photoactive molecular complexes with an intermediate strength of the coupling to their environment. Both techniques are highly efficient, yet they involve approximations at different levels. By explicitly comparing with the numerically exact quasiadiabatic path integral approach, we show for the Fenna-Matthews-Olson complex that the CMRT describes the decay rates in the population dynamics well, but final stationary populations and the oscillation frequencies differ slightly. In addition, we use the combined CMRT+PMA to calculate two-dimensional photon-echo spectra for a simple dimer model. We find excellent agreement with the exact path integral calculations at short waiting times where the dynamics is still coherent. For long waiting times, differences occur due to different final stationary states, specifically for strong system-bath coupling. For weak to intermediate system-bath couplings, which is most important for natural photosynthetic complexes, the combined CMRT+PMA gives reasonable results with acceptable computational efforts.

Vibrational beatings conceal evidence of electronic coherence in the FMO light-harvesting complex

In biological light harvesting, solar energy is captured by photosynthetic antennae for subsequent storage into chemical bonds. The remarkable efficiency reached in transferring the energy between the collection and storage events recently has been attributed to long-lived electronic coherence present in such antennae systems. We present numerical simulations indicating that the spectroscopic transients that supported this hypothesis are not induced by electronic coherence but instead are caused by vibrational (nuclear) motion in the electronic ground state potential. Besides emphasizing the significance of such nuclear modes, our findings stimulate a reconsideration of the role of electronic coherence in promoting energy transfer in natural photosynthesis. Furthermore, they require us to rethink how energy transfer efficiency is reflected in spectral signals.

Theoretical study on excited states of bacteriochlorophyll a in solutions with density functional assessment

The excited-state properties of bacteriochlorophyll (BChl) a in triethylamine, 1-propanol, and methanol are investigated with the time-dependent density functional theory by using the quantum mechanical and molecular mechanical reweighting free energy self-consistant field method. It is found that no prevalent density functionals can reproduce the experimental excited-state properties, i.e., the absorption and reorganization energies, of BChl a in the solutions. The parameter μ in the range-separated hybrid functional is therefore optimized to reproduce the differences of the absorption energies in the solutions. We examine the origin of the differences of the absorption energies in the solutions and find that sensitive balance between contributions of structural changes and solute-solvent interactions determines the differences. The accurate description of the excitation with the density functional with the adjusted parameter is therefore essential to the understanding of the excited-state properties of BChl a in proteins and also the mechanism of the photosynthetic systems.

Attenuation of excitation decay rate due to collective effect

We study a series of N oscillators, each coupled to its nearest neighbors, and linearly to a phonon field through the oscillator's number operator. We show that the Hamiltonian of a pair of adjacent oscillators, or a dimer, within the series of oscillators can be transformed into a form in which they are collectively coupled to the phonon field as a composite unit. In the weak coupling and rotating-wave approximation, the system behaves effectively as the trilinear boson model in the one excitation subspace of the dimer subsystem. The reduced dynamics of the one excitation subspace of the dimer subsystem coupled weakly to a phonon bath is similar to that of a two-level system, with a metastable state against the vacuum. The decay constant of the subsystem is proportional to the dephasing rate of the individual oscillator in a phonon bath, attenuated by a factor that depends on site asymmetry, intersite coupling, and the resonance frequency between the transformed oscillator modes, or excitons. As a result of the collective effect, the excitation relaxation lifetime is prolonged over the dephasing lifetime of an individual oscillator coupled to the same bath.

Enhancement of vibronic and ground-state vibrational coherences in 2D spectra of photosynthetic complexes

A vibronic-exciton model is applied to investigate the recently proposed mechanism of enhancement of coherent oscillations due to mixing of electronic and nuclear degrees of freedom. We study a dimer system to elucidate the role of resonance coupling, site energies, vibrational frequency and energy disorder in the enhancement of vibronic-exciton and ground-state vibrational coherences, and to identify regimes where this enhancement is significant. For a heterodimer representing two coupled bachteriochloropylls of the FMO complex, long-lived vibronic coherences are found to be generated only when the frequency of the mode is in the vicinity of the electronic energy difference. Although the vibronic-exciton coherences exhibit a larger initial amplitude compared to the ground-state vibrational coherences, we conclude that, due to the dephasing of the former, both type of coherences have a similar magnitude at longer population time.

Designs for molecular circuits that use electronic coherence

The mounting evidence of recent years regarding long-lived coherent dynamics of electronic excitations in several light-harvesting antenna proteins suggests the possibility of realizing and exploiting light-initiated quantum dynamics in synthetic molecular devices based on electronic energy transfer. Inspired by the field of molecular logic, we focus this discussion on the prospect of using quantum coherence to control the direction of energy flow in a molecular circuit. As a prototype system we consider a circuit consisting of three chromophores that deliver energy to two trap chromophores. Our aim is to control to which trap the energy is more likely to be delivered. This is achieved by switching one of the circuit chromophores ON and OFF from the system, such that the direction of energy flow substantially changes from the ON and OFF states of the circuit. We find that quantum coherence can allow a significant ability to direct energy transfer in the circuit. However, when realistic levels of noise are added, quantum coherence only slightly improves the ability to direct electronic energy in comparison to a classical hopping mechanism.

Spatial and Electronic Correlations in the PE545 Light-Harvesting Complex

The recent discovery of long-lasting quantum coherence effects in photosynthetic pigment-protein complexes has challenged our view of the role that protein motions play in light-harvesting processes. Several groups have suggested that correlated fluctuations involving the pigments site energies and couplings could be at the origin of such unexpected behavior. Here we combine molecular dynamics simulations with quantum mechanics/molecular mechanics calculations to analyze the degree of correlated fluctuations in the PE545 complex of Rhodomonas sp. strain CS24. We find that correlations between the motions of the chromophores, which are significantly assisted by the water solvent, do not translate into appreciable site energy correlations but do lead to significant cross-correlations of energies and couplings. Such behavior, not observed in a recent study on the Fenna-Mathews-Olson complex, seems to provide phycobiliproteins with an additional fundamental mechanism to control quantum coherence and light-harvesting efficiency compared with chlorophyll-containing complexes.

Normal mode analysis of the spectral density of the Fenna-Matthews-Olson light-harvesting protein: how the protein dissipates the excess energy of excitons

We report a method for the structure-based calculation of the spectral density of the pigment-protein coupling in light-harvesting complexes that combines normal-mode analysis with the charge density coupling (CDC) and transition charge from electrostatic potential (TrEsp) methods for the computation of site energies and excitonic couplings, respectively. The method is applied to the Fenna-Matthews-Olson (FMO) protein in order to investigate the influence of the different parts of the spectral density as well as correlations among these contributions on the energy transfer dynamics and on the temperature-dependent decay of coherences. The fluctuations and correlations in excitonic couplings as well as the correlations between coupling and site energy fluctuations are found to be 1 order of magnitude smaller in amplitude than the site energy fluctuations. Despite considerable amplitudes of that part of the spectral density which contains correlations in site energy fluctuations, the effect of these correlations on the exciton population dynamics and dephasing of coherences is negligible. The inhomogeneous charge distribution of the protein, which causes variations in local pigment-protein coupling constants of the normal modes, is responsible for this effect. It is seen thereby that the same building principle that is used by nature to create an excitation energy funnel in the FMO protein also allows for efficient dissipation of the excitons' excess energy.

Origin of long-lived coherences in light-harvesting complexes

A vibronic exciton model is applied to explain the long-lived oscillatory features in the two-dimensional (2D) electronic spectra of the Fenna-Matthews-Olson (FMO) complex. Using experimentally determined parameters and uncorrelated site energy fluctuations, the model predicts oscillations with dephasing times of 1.3 ps at 77 K, which is in a good agreement with the experimental results. These long-lived oscillations originate from the coherent superposition of vibronic exciton states with dominant contributions from vibrational excitations on the same pigment. The oscillations obtain a large amplitude due to excitonic intensity borrowing, which gives transitions with strong vibronic character a significant intensity despite the small Huang-Rhys factor. Purely electronic coherences are found to decay on a 200 fs time scale.

Weak exciton scattering in molecular nanotubes revealed by double-quantum two-dimensional electronic spectroscopy

The two-exciton manifold of a double-wall cylindrical molecular aggregate is studied using a coherent third order optical technique. Experiments reveal the anharmonic character of the exciton bands. Atomistic simulations of the exciton-exciton scattering show that the excitons can be treated as weakly coupled hard-core bosons. The weak coupling stems from the extended exciton delocalization made possible by the nanotube geometry.

Atomistic study of the long-lived quantum coherences in the Fenna-Matthews-Olson complex

A remarkable amount of theoretical research has been carried out to elucidate the physical origins of the recently observed long-lived quantum coherence in the electronic energy transfer process in biological photosynthetic systems. Although successful in many respects, several widely used descriptions only include an effective treatment of the protein-chromophore interactions. In this work, by combining an all-atom molecular dynamics simulation, time-dependent density functional theory, and open quantum system approaches, we successfully simulate the dynamics of the electronic energy transfer of the Fenna-Matthews-Olson pigment-protein complex. The resulting characteristic beating of populations and quantum coherences is in good agreement with the experimental results and the hierarchy equation of motion approach. The experimental absorption, linear, and circular dichroism spectra and dephasing rates are recovered at two different temperatures. In addition, we provide an extension of our method to include zero-point fluctuations of the vibrational environment. This work thus presents, to our knowledge, one of the first steps to explain the role of excitonic quantum coherence in photosynthetic light-harvesting complexes based on their atomistic and molecular description.

Exact quantum master equation for a molecular aggregate coupled to a harmonic bath

We consider a molecular aggregate consisting of N identical monomers. Each monomer comprises two electronic levels and a single harmonic mode. The monomers interact with each other via dipole-dipole forces. The monomer vibrational modes are bilinearly coupled to a bath of harmonic oscillators. This is a prototypical model for the description of coherent exciton transport, from quantum dots to photosynthetic antennae. We derive an exact quantum master equation for such systems. Computationally, the master equation may be useful for the testing of various approximations employed in theories of quantum transport. Physically, it offers a plausible explanation of the origins of long-lived coherent optical responses of molecular aggregates in dissipative environments.

From atomistic modeling to excitation transfer and two-dimensional spectra of the FMO light-harvesting complex

The experimental observation of long-lived quantum coherences in the Fenna-Matthews-Olson (FMO) light-harvesting complex at low temperatures has challenged general intuition in the field of complex molecular systems and provoked considerable theoretical effort in search of explanations. Here we report on room-temperature calculations of the excited-state dynamics in FMO using a combination of molecular dynamics simulations and electronic structure calculations. Thus we obtain trajectories for the Hamiltonian of this system which contains time-dependent vertical excitation energies of the individual bacteriochlorophyll molecules and their mutual electronic couplings. The distribution of energies and couplings is analyzed together with possible spatial correlations. It is found that in contrast to frequent assumptions the site energy distribution is non-Gaussian. In a subsequent step, averaged wave packet dynamics is used to determine the exciton dynamics in the system. Finally, with the time-dependent Hamiltonian, linear and two-dimensional spectra are determined. The thus-obtained linear absorption line shape agrees well with experimental observation and is largely determined by the non-Gaussian site energy distribution. The two-dimensional spectra are in line with what one would expect by extrapolation of the experimental observations at lower temperatures and indicate almost total loss of long-lived coherences.