Spectra and Dynamics in the B800 Antenna: Comparing Hierarchical Equations, Redfield and Förster Theories

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.

Quantum chemical elucidation of a sevenfold symmetric bacterial antenna complex

The light-harvesting complex 2 (LH2) of purple bacteria is one of the most studied photosynthetic antenna complexes. Its symmetric structure and ring-like bacteriochlorophyll arrangement make it an ideal system for theoreticians and spectroscopists. LH2 complexes from most bacterial species are thought to have eightfold or ninefold symmetry, but recently a sevenfold symmetric LH2 structure from the bacterium Mch. purpuratum was solved by Cryo-Electron microscopy. This LH2 also possesses unique near-infrared absorption and circular dichroism (CD) spectral properties. Here we use an atomistic strategy to elucidate the spectral properties of Mch. purpuratum LH2 and understand the differences with the most commonly studied LH2 from Rbl. acidophilus. Our strategy exploits a combination of molecular dynamics simulations, multiscale polarizable quantum mechanics/molecular mechanics calculations, and lineshape simulations. Our calculations reveal that the spectral properties of LH2 complexes are tuned by site energies and exciton couplings, which in turn depend on the structural fluctuations of the bacteriochlorophylls. Our strategy proves effective in reproducing the absorption and CD spectra of the two LH2 complexes, and in uncovering the origin of their differences. This work proves that it is possible to obtain insight into the spectral tuning strategies of purple bacteria by quantitatively simulating the spectral properties of their antenna complexes.

The location of the low-energy states in Lhca1 favors excitation energy transfer to the core in the plant PSI-LHCI supercomplex

Lhca1 is one of the four pigment-protein complexes composing the outer antenna of plant Photosystem I-light-havesting I supercomplex (PSI-LHCI). It forms a functional dimer with Lhca4 but, differently from this complex, it does not contain 'red-forms,' i.e., pigments absorbing above 700 nm. Interestingly, the recent PSI-LHCI structures suggest that Lhca1 is the main point of delivering the energy harvested by the antenna to the core. To identify the excitation energy pathways in Lhca1, we developed a structure-based exciton model based on the simultaneous fit of the low-temperature absorption, linear dichroism, and fluorescence spectra of wild-type Lhca1 and two mutants, lacking chlorophylls contributing to the long-wavelength region of the absorption. The model enables us to define the locations of the lowest energy pigments in Lhca1 and estimate pathways and timescales of energy transfer within the complex and to the PSI core. We found that Lhca1 has a particular energy landscape with an unusual (compared to Lhca4, LHCII, and CP29) configuration of the low-energy states. Remarkably, these states are located near the core, facilitating direct energy transfer to it. Moreover, the low-energy states of Lhca1 are also coupled to the red-most state (red forms) of the neighboring Lhca4 antenna, providing a pathway for effective excitation energy transfer from Lhca4 to the core.

Stark absorption and Stark fluorescence spectroscopies: Theory and simulations

Stark spectroscopy experiments are widely used to study the properties of molecular systems, particularly those containing charge-transfer (CT) states. However, due to the small transition dipole moments and large static dipole moments of the CT states, the standard interpretation of the Stark absorption and Stark fluorescence spectra in terms of the Liptay model may be inadequate. In this work, we provide a theoretical framework for calculations of Stark absorption and Stark fluorescence spectra and propose new methods of simulations that are based on the quantum-classical theory. In particular, we use the forward-backward trajectory solution and a variant of the Poisson bracket mapping equation, which have been recently adapted for the calculation of conventional (field-free) absorption and fluorescence spectra. For comparison, we also apply the recently proposed complex time-dependent Redfield theory, while exact results are obtained using the hierarchical equations of motion approach. We show that the quantum-classical methods produce accurate results for a wide range of systems, including those containing CT states. The CT states contribute significantly to the Stark spectra, and the standard Liptay formalism is shown to be inapplicable for the analysis of spectroscopic data in those cases. We demonstrate that states with large static dipole moments may cause a pronounced change in the total fluorescence yield of the system in the presence of an external electric field. This effect is correctly captured by the quantum-classical methods, which should therefore prove useful for further studies of Stark spectra of real molecular systems. As an example, we calculate the Stark spectra for the Fenna-Matthews-Olson complex of green sulfur bacteria.

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.

Quantum-Classical Approach for Calculations of Absorption and Fluorescence: Principles and Applications

Absorption and fluorescence spectroscopy techniques provide a wealth of information on molecular systems. The simulations of such experiments remain challenging, however, despite the efforts put into developing the underlying theory. An attractive method of simulating the behavior of molecular systems is provided by the quantum-classical theory─it enables one to keep track of the state of the bath explicitly, which is needed for accurate calculations of fluorescence spectra. Unfortunately, until now there have been relatively few works that apply quantum-classical methods for modeling spectroscopic data. In this work, we seek to provide a framework for the calculations of absorption and fluorescence lineshapes of molecular systems using the methods based on the quantum-classical Liouville equation, namely, the forward-backward trajectory solution (FBTS) and the non-Hamiltonian variant of the Poisson bracket mapping equation (PBME-nH). We perform calculations on a molecular dimer and the photosynthetic Fenna-Matthews-Olson complex. We find that in the case of absorption, the FBTS outperforms PBME-nH, consistently yielding highly accurate results. We next demonstrate that for fluorescence calculations, the method of choice is a hybrid approach, which we call PBME-nH-Jeff, that utilizes the effective coupling theory [Gelzinis, A.; J. Chem. Phys. 2020, 152, 051103] to estimate the excited state equilibrium density operator. Thus, we find that FBTS and PBME-nH-Jeff are excellent candidates for simulating, respectively, absorption and fluorescence spectra of real molecular systems.

Artificial Photosynthesis: Is Computation Ready for the Challenge Ahead?

A tremendous effort is currently devoted to the generation of novel hybrid materials with enhanced electronic properties for the creation of artificial photosynthetic systems. This compelling and challenging problem is well-defined from an experimental point of view, as the design of such materials relies on combining organic materials or metals with biological systems like light harvesting and redox-active proteins. Such hybrid systems can be used, e.g., as bio-sensors, bio-fuel cells, biohybrid photoelectrochemical cells, and nanostructured photoelectronic devices. Despite these efforts, the main bottleneck is the formation of efficient interfaces between the biological and the organic/metal counterparts for efficient electron transfer (ET). It is within this aspect that computation can make the difference and improve the current understanding of the mechanisms underneath the interface formation and the charge transfer efficiency. Yet, the systems considered (i.e., light harvesting protein, self-assembly monolayer and surface assembly) are more and more complex, reaching (and often passing) the limit of current computation power. In this review, recent developments in computational methods for studying complex interfaces for artificial photosynthesis will be provided and selected cases discussed, to assess the inherent ability of computation to leave a mark in this field of research.

Absorption and Circular Dichroism Spectra of Molecular Aggregates With the Full Cumulant Expansion

The exciton Hamiltonian of multichromophoric aggregates can be probed by spectroscopic techniques such as linear absorption and circular dichroism. To compare calculated Hamiltonians to experiments, a lineshape theory is needed, which takes into account the coupling of the excitons with inter- and intramolecular vibrations. This coupling is normally introduced in a perturbative way through the cumulant expansion formalism and further approximated by assuming a Markovian exciton dynamics, for example with the modified Redfield theory. Here, we present the implementation of the full cumulant expansion (FCE) formalism (J. Chem. Phys.142, 2015, 09410625747060) to efficiently compute absorption and circular dichroism spectra of molecular aggregates beyond the Markov approximation, without restrictions on the form of exciton–phonon coupling. By employing the LH2 system of purple bacteria as a challenging test case, we compare the FCE lineshapes with the Markovian lineshapes obtained with the modified Redfield theory, showing that the latter presents a less satisfying agreement with experiments. The FCE approach instead accurately describes the lineshapes, especially in the vibronic sideband of the B800 peak. We envision that the FCE approach will become a valuable tool for accurately comparing model exciton Hamiltonians with optical spectroscopy experiments.

Coherent and incoherent theories for photosynthetic energy transfer

There is a remarkable characteristic of photosynthesis in nature, that is, the energy transfer efficiency is close to 100%. Recently, due to the rapid progress made in the experimental techniques, quantum coherent effects have been experimentally demonstrated. Traditionally, the incoherent theories are capable of calculating the energy transfer efficiency, e.g., (generalized) Förster theory and modified Redfield theory (MRT). However, in order to describe the quantum coherent effects in photosynthesis, one has to exploit coherent theories, such as hierarchical equation of motion (HEOM), quantum path integral, coherent modified Redfield theory (CMRT), small-polaron quantum master equation, and general Bloch-Redfield theory in addition to the Redfield theory. Here, we summarize the main points of the above approaches, which might be beneficial to the quantum simulation of quantum dynamics of exciton energy transfer (EET) in natural photosynthesis, and shed light on the design of artificial light-harvesting devices.

Quantum Chemical Modeling of the Photoinduced Activity of Multichromophoric Biosystems

Multichromophoric biosystems represent a broad family with very diverse members, ranging from light-harvesting pigment-protein complexes to nucleic acids. The former are designed to capture, harvest, efficiently transport, and transform energy from sunlight for photosynthesis, while the latter should dissipate the absorbed radiation as quickly as possible to prevent photodamages and corruption of the carried genetic information. Because of the unique electronic and structural characteristics, the modeling of their photoinduced activity is a real challenge. Numerous approaches have been devised building on the theoretical development achieved for single chromophores and on model Hamiltonians that capture the essential features of the system. Still, a question remains: is a general strategy for the accurate modeling of multichromophoric systems possible? By using a quantum chemical point of view, here we review the advancements developed so far highlighting differences and similarities with the single chromophore treatment. Finally, we outline the important limitations and challenges that still need to be tackled to reach a complete and accurate picture of their photoinduced properties and dynamics.

Successes & challenges in the atomistic modeling of light-harvesting and its photoregulation

Light-harvesting is a crucial step of photosynthesis. Its mechanisms and related energetics have been revealed by a combination of experimental investigations and theoretical modeling. The success of theoretical modeling is largely due to the application of atomistic descriptions combining quantum chemistry, classical models and molecular dynamics techniques. Besides the important achievements obtained so far, a complete and quantitative understanding of how the many different light-harvesting complexes exploit their structural specificity is still missing. Moreover, many questions remain unanswered regarding the mechanisms through which light-harvesting is regulated in response to variable light conditions. Here we show that, in both fields, a major role will be played once more by atomistic descriptions, possibly generalized to tackle the numerous time and space scales on which the regulation takes place: going from the ultrafast electronic excitation of the multichromophoric aggregate, through the subsequent conformational changes in the embedding protein, up to the interaction between proteins.

How Well Does the Hole-Burning Action Spectrum Represent the Site-Distribution Function of the Lowest-Energy State in Photosynthetic Pigment-Protein Complexes?

For the first time, we combined Monte Carlo and nonphotochemical hole burning (NPHB) master equation approaches to allow for ultrahigh-resolution (<0.005 cm^-1, smaller than the typical homogeneous line widths at 5 K) simulations of the NPHB spectra of dimers and trimers of interacting pigments. These simulations reveal significant differences between the zero-phonon hole (ZPH) action spectrum and the site-distribution function (SDF) of the lowest-energy state. The NPHB of the lowest-energy pigment, following the excitation energy transfer (EET) from the higher-energy pigments which are excited directly, results in the shifts of all excited states. These shifts affect the ZPH action spectra and EET times derived from the widths of the spectral holes burned in the donor-dominated regions. The effect is present for a broad variety of realistic antihole functions, and it is maximal at relatively low values of interpigment coupling (V ≤ 5 cm^-1) where the use of the Förster approximation is justified. These findings need to be considered in interpreting various optical spectra of photosynthetic pigment-protein complexes for which SDFs (describing the inhomogeneous broadening) are often obtained directly from the ZPH action spectra. Water-soluble chlorophyll-binding protein (WSCP) was considered as an example.

Dynamics of the mixed exciton and charge-transfer states in light-harvesting complex Lhca4: Hierarchical equation approach

We model the energy transfer dynamics in the Lhca4 peripheral antenna of photosystem I from higher plants. Equilibration between the bulk exciton levels of the antenna and the red-shifted charge-transfer (CT) states is described using the numerically inexpensive Redfield-Förster approach and exact hierarchical equation (HEOM) method. We propose a compartmentalization scheme allowing a quantitatively correct description of the dynamics with the Redfield-Förster theory, including the exciton-type relaxation within strongly coupled compartments and hopping-type migration between them. The Redfield-Förster method gives the kinetics close to the HEOM solution when treating the CT state as dynamically localized. We also demonstrate that the excited states strongly coupled with the CT should be considered as localized as well.

Exciton-vibrational resonance and dynamics of charge separation in the photosystem II reaction center

The dynamics of charge separation in the photosystem II reaction center (PSII-RC) in the presence of intramolecular vibrations with their frequency matching the energy gap between the exciton state acting as the primary electron donor and the first charge-transfer (CT) state are investigated. A reduced PSII-RC 4-state model explicitly including a CT state is analyzed within Redfield relaxation theory in the multidimensional exciton-vibrational (vibronic) basis. This model is used to study coherent energy/electron transfers and their spectral signatures obtained by two-dimensional electronic spectroscopy (2DES). Modeling of the time-resolved 2D frequency maps obtained by wavelet analysis reveals the origins of the coherences which produce the observed oscillating features in 2DES and allows comparing the lifetimes of the coherences. The results suggest faster excitonic decoherence as compared with longer-lived vibronic oscillations. The emerging picture of the dynamics unravels the role of resonant vibrations in sustaining the effective energy conversion in the PSII-RC. We demonstrate that the mixing of the exciton and CT states promoted by a resonant vibrational quantum allows faster penetration of excitation energy into the CT with subsequent dynamic localization at the bottom of the CT potential induced by the remaining non-resonant nuclear modes. The degree of vibration-assisted mixing and, correspondingly, the rate of primary charge separation, increases significantly in the case of electron-vibrational resonance. The observed features illustrate the principles of quantum design of the photosynthetic unit. These principles are connected with the phenomenon of coherent mixing within vibronic eigenstates, increasing the effectiveness of charge separation not only upon coherent and impulsive laser excitation utilized in the 2DES experiment, but also under natural conditions under non-coherent non-impulsive solar light illumination.

Robust light harvesting by a noisy antenna

Photosynthetic light harvesting can be very efficient in solar energy conversion while taking place in a highly disordered and noisy physiological environment. This efficiency is achieved by the ultrafast speed of the primary photosynthetic processes, which is enabled by a delicate interplay of quantum effects, thermodynamics and environmental noise. The primary processes take place in light-harvesting antennas built from pigments bound to a fluctuating protein scaffold. Here, we employ ultrafast single-molecule spectroscopy to follow fluctuations of the femtosecond energy transfer times in individual LH2 antenna complexes of purple bacteria. By combining single molecule results with ensemble spectroscopy through a unified theoretical description of both, we show how the protein fluctuations alter the excitation energy transfer dynamics. We find that from the thirteen orders of magnitude of possible timescales from picoseconds to minutes, the relevant fluctuations occur predominantly on a biological timescale of seconds, i.e. in the domain of slow protein motion. The measured spectra and dynamics can be explained by the protein modulating pigment excitation energies only. Moreover, we find that the small spread of pigment mean energies allows for excitation delocalization between the coupled pigments to survive. These unique features provide fast energy transport even in the presence of disorder. We conclude that this is the mechanism that enables LH2 to operate as a robust light-harvester, in spite of its intrinsically noisy biological environment.

The Role of Resonant Vibrations in Electronic Energy Transfer

Nuclear vibrations play a prominent role in the spectroscopy and dynamics of electronic systems. As recent experimental and theoretical studies suggest, this may be even more so when vibrational frequencies are resonant with transitions between the electronic states. Herein, a vibronic multilevel Redfield model is reported for excitonically coupled electronic two-level systems with a few explicitly included vibrational modes and interacting with a phonon bath. With numerical simulations the effects of the quantized vibrations on the dynamics of energy transfer and coherence in a model dimer are illustrated. The resonance between the vibrational frequency and energy gap between the sites leads to a large delocalization of vibronic states, which then results in faster energy transfer and longer-lived mixed coherences.

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.

Mixing of exciton and charge-transfer states in light-harvesting complex Lhca4

Structure-based modeling of spectra of the wild-type Lhca4 and NH mutant enables us to build the exciton model of the complex that includes a charge-transfer state mixed with the excited-state manifold.

On the accuracy of coherent modified Redfield theory in simulating excitation energy transfer dynamics

In this study, we investigate the accuracy of a recently developed coherent modified Redfield theory (CMRT) in simulating excitation energy transfer (EET) dynamics. The CMRT is a secular non-Markovian quantum master equation that is derived by extending the modified Redfield theory to treat coherence dynamics in molecular excitonic systems. Herein, we systematically survey the applicability of the CMRT in a large EET parameter space through the comparisons of the CMRT EET dynamics in a dimer system with the numerically exact results. The results confirm that the CMRT exhibits a broad applicable range and allow us to locate the specific parameter regimes where CMRT fails to provide adequate results. Moreover, we propose an accuracy criterion based on the magnitude of second-order perturbation to characterize the applicability of CMRT and show that the criterion summarizes all the benchmark results and the physics described by CMRT. Finally, we employ the accuracy criterion to quantitatively compare the performance of CMRT to that of a small polaron quantum master equation approach. The comparison demonstrates the complementary nature of these two methods, and as a result, the combination of the two methods provides accurate simulations of EET dynamics for the full parameter space investigated in this study. Our results not only delicately evaluate the applicability of the CMRT but also reveal new physical insights for factors controlling the dynamics of EET that should be useful for developing more accurate and efficient methods for simulations of EET dynamics in molecular aggregate systems.