Towards an exact theory of linear absorbance and circular dichroism of pigment-protein complexes: Importance of non-secular contributions

Transient Chiral Dynamics in the Fenna-Matthews-Olson Complex Revealed by Two-Dimensional Circular Dichroism Spectroscopy

Chirality plays a pivotal role across scientific disciplines with profound implications spanning light-matter interactions, molecular recognition, and natural evolutionary processes. This study delves into the active influence of molecular chirality on exciton energy transfer within photosynthetic protein complexes, focusing on the Fenna-Matthews-Olson (FMO) complex. Employing two-dimensional circular dichroism (2DCD) spectroscopy, we investigate the transient chiral dynamics of excitons during energy transfer processes within the FMO complex. Our approach, incorporating pulse information into population dynamics based on the third-order response function, facilitates the calculation of 2DCD spectra and dynamics. This enables the extraction of chiral contributions to excitonic energy transfer and the examination of electronic wave functions. We demonstrate that 2DCD spectra offer excitation energies that are better resolved than those from conventional two-dimensional electronic spectroscopy. These findings deepen our understanding of exciton energy transfer mechanisms in natural photosynthesis, emphasizing the potential of 2DCD spectroscopy as a powerful tool for unraveling the chiral contribution to exciton dynamics.

First-principles simulation of excitation energy transfer and transient absorption spectroscopy in the CP29 light-harvesting complex

The dynamics of delocalized excitons in light-harvesting complexes (LHCs) can be investigated using different experimental techniques, and transient absorption (TA) spectroscopy is one of the most valuable methods for this purpose. A careful interpretation of TA spectra is essential for the clarification of excitation energy transfer (EET) processes occurring during light-harvesting. However, even in the simplest LHCs, a physical model is needed to interpret transient spectra as the number of EET processes occurring at the same time is very large to be disentangled from measurements alone. Physical EET models are commonly built by fittings of the microscopic exciton Hamiltonians and exciton-vibrational parameters, an approach that can lead to biases. Here, we present a first-principles strategy to simulate EET and transient absorption spectra in LHCs, combining molecular dynamics and accurate multiscale quantum chemical calculations to obtain an independent estimate of the excitonic structure of the complex. The microscopic parameters thus obtained are then used in EET simulations to obtain the population dynamics and the related spectroscopic signature. We apply this approach to the CP29 minor antenna complex of plants for which we follow the EET dynamics and transient spectra after excitation in the chlorophyll b region. Our calculations reproduce all the main features observed in the transient absorption spectra and provide independent insight on the excited-state dynamics of CP29. The approach presented here lays the groundwork for the accurate simulation of EET and unbiased interpretation of transient spectra in multichromophoric systems.

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.

Signatures of intramolecular vibrational and vibronic Q[Formula: see text]-Q[Formula: see text] coupling effects in absorption and CD spectra of chlorophyll dimers

An electron-vibrational coupling model that includes the vibronic (non-adiabatic) coupling between the Q[Formula: see text] and Q[Formula: see text] transitions of chlorophyll (Chl), created by Reimers and coworkers (Scientific Rep. 3, 2761, 2013) is extended here to chlorophyll dimers with interchlorophyll excitonic coupling. The model is applied to a Chl a dimer of the water-soluble chlorophyll binding protein (WSCP). As for isolated chlorophyll, the vibronic coupling is found to have a strong influence on the high-frequency vibrational sideband in the absorption spectrum, giving rise to a band splitting. In contrast, in the CD spectrum the interplay of vibronic coupling and static disorder leads to a strong suppression of the vibrational sideband in excellent agreement with the experimental data. The conservative nature of the CD spectrum in the low-energy region is found to be caused by a delicate balance of the intermonomer excitonic coupling between the purely electronic Q[Formula: see text] transition and the Q[Formula: see text] transition involving intramolecular vibrational excitations on one hand and the coupling to higher-energy electronic transitions on the other hand.

Recent progress in atomistic modeling of light-harvesting complexes: a mini review

In this mini review, we focus on recent advances in the atomistic modeling of biological light-harvesting (LH) complexes. Because of their size and sophisticated electronic structures, multiscale methods are required to investigate the dynamical and spectroscopic properties of such complexes. The excitation energies, in this context also known as site energies, excitonic couplings, and spectral densities are key quantities which usually need to be extracted to be able to determine the exciton dynamics and spectroscopic properties. The recently developed multiscale approach based on the numerically efficient density functional tight-binding framework followed by excited state calculations has been shown to be superior to the scheme based on pure classical molecular dynamics simulations. The enhanced approach, which improves the description of the internal vibrational dynamics of the pigment molecules, yields spectral densities in good agreement with the experimental counterparts for various bacterial and plant LH systems. Here, we provide a brief overview of those results and described the theoretical foundation of the multiscale protocol.

Fingerprint and Universal Markovian Closure of Structured Bosonic Environments

We exploit the properties of chain mapping transformations of bosonic environments to identify a finite collection of modes able to capture the characteristic features, or fingerprint, of the environment. Moreover we show that the countable infinity of residual bath modes can be replaced by a universal Markovian closure, namely, a small collection of damped modes undergoing a Lindblad-type dynamics whose parametrization is independent of the spectral density under consideration. We show that the Markovian closure provides a quadratic speedup with respect to standard chain mapping techniques and makes the memory requirement independent of the simulation time, while preserving all the information on the fingerprint modes. We illustrate the application of the Markovian closure to the computation of linear spectra but also to nonlinear spectral response, a relevant experimentally accessible many body coherence witness for which efficient numerically exact calculations in realistic environments are currently lacking.

Accuracy of approximate methods for the calculation of absorption-type linear spectra with a complex system-bath coupling

The accuracy of approximate methods for calculating linear optical spectra depends on many variables. In this study, we fix most of these parameters to typical values found in photosynthetic light-harvesting complexes of plants and determine the accuracy of approximate spectra with respect to exact calculation as a function of the energy gap and interpigment coupling in a pigment dimer. We use a spectral density with the first eight intramolecular modes of chlorophyll a and include inhomogeneous disorder for the calculation of spectra. We compare the accuracy of absorption, linear dichroism, and circular dichroism spectra calculated using the Full Cumulant Expansion (FCE), coherent time-dependent Redfield (ctR), and time-independent Redfield and modified Redfield methods. As a reference we use spectra calculated with the Exact Stochastic Path Integral Evaluation method. We find the FCE method to be the most accurate for the calculation of all spectra. The ctR method performs well for the qualitative calculation of absorption and linear dichroism spectra when pigments are moderately coupled (∼15 cm-1), but ctR spectra may differ significantly from exact spectra when strong interpigment coupling (>100 cm-1) is present. The dependence of the quality of Redfield and modified Redfield spectra on molecular parameters is similar, and these methods almost always perform worse than ctR, especially when the interpigment coupling is strong or the excitonic energy gap is small (for a given coupling). The accuracy of approximate spectra is not affected by resonance with intramolecular modes for typical system-bath coupling and disorder values found in plant light-harvesting complexes.

Exact simulation of pigment-protein complexes unveils vibronic renormalization of electronic parameters in ultrafast spectroscopy

The primary steps of photosynthesis rely on the generation, transport, and trapping of excitons in pigment-protein complexes (PPCs). Generically, PPCs possess highly structured vibrational spectra, combining many discrete intra-pigment modes and a quasi-continuous of protein modes, with vibrational and electronic couplings of comparable strength. The intricacy of the resulting vibronic dynamics poses significant challenges in establishing a quantitative connection between spectroscopic data and underlying microscopic models. Here we show how to address this challenge using numerically exact simulation methods by considering two model systems, namely the water-soluble chlorophyll-binding protein of cauliflower and the special pair of bacterial reaction centers. We demonstrate that the inclusion of the full multi-mode vibronic dynamics in numerical calculations of linear spectra leads to systematic and quantitatively significant corrections to electronic parameter estimation. These multi-mode vibronic effects are shown to be relevant in the longstanding discussion regarding the origin of long-lived oscillations in multidimensional nonlinear spectra.

Multimode vibronic mixing in model photosynthetic systems revealed by numerically exact simulations is shown to strongly modify linear and non-linear optical responses and facilitate the persistence of coherent dynamics.

The atomistic modeling of light-harvesting complexes from the physical models to the computational protocol

The function of light-harvesting complexes is determined by a complex network of dynamic interactions among all the different components: the aggregate of pigments, the protein, and the surrounding environment. Complete and reliable predictions on these types of composite systems can be only achieved with an atomistic description. In the last few decades, there have been important advances in the atomistic modeling of light-harvesting complexes. These advances have involved both the completeness of the physical models and the accuracy and effectiveness of the computational protocols. In this Perspective, we present an overview of the main theoretical and computational breakthroughs attained so far in the field, with particular focus on the important role played by the protein and its dynamics. We then discuss the open problems in their accurate modeling that still need to be addressed. To illustrate an effective computational workflow for the modeling of light harvesting complexes, we take as an example the plant antenna complex CP29 and its H111N mutant.

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.

Time-dependent atomistic simulations of the CP29 light-harvesting complex

Light harvesting as the first step in photosynthesis is of prime importance for life on earth. For a theoretical description of photochemical processes during light harvesting, spectral densities are key quantities. They serve as input functions for modeling the excitation energy transfer dynamics and spectroscopic properties. Herein, a recently developed procedure is applied to determine the spectral densities of the pigments in the minor antenna complex CP29 of photosystem II, which has recently gained attention because of its active role in non-photochemical quenching processes in higher plants. To this end, the density functional-based tight binding (DFTB) method has been employed to enable simulation of the ground state dynamics in a quantum-mechanics/molecular mechanics (QM/MM) scheme for each chlorophyll pigment. Subsequently, the time-dependent extension of the long-range corrected DFTB approach has been used to obtain the excitation energy fluctuations along the ground-state trajectories also in a QM/MM setting. From these results, the spectral densities have been determined and compared for different force fields and to spectral densities from other light-harvesting complexes. In addition, time-dependent and time-independent excitonic Hamiltonians of the system have been constructed and applied to the determination of absorption spectra as well as exciton dynamics.

Normal mode analysis of spectral density of FMO trimers: Intra- and intermonomer energy transfer

The intermolecular contribution to the spectral density of the exciton-vibrational coupling of the homotrimeric Fenna-Matthews-Olson (FMO) light-harvesting protein of green sulfur bacteria P. aestuarii is analyzed by combining a normal mode analysis of the protein with the charge density coupling method for the calculation of local transition energies of the pigments. Correlations in site energy fluctuations across the whole FMO trimer are found at low vibrational frequencies. Including, additionally, the high-frequency intrapigment part of the spectral density, extracted from line-narrowing spectra, we study intra- and intermonomer exciton transfer. Whereas the intrapigment part of the spectral density is important for fast intramonomer exciton relaxation, the intermolecular contributions (due to pigment-environment coupling) determine the intermonomer exciton transfer. Neither the variations of the local Huang-Rhys factors nor the correlations in site energy fluctuations have a critical influence on energy transfer. At room temperature, the intermonomer transfer in the FMO protein occurs on a 10 ps time scale, whereas intramonomer exciton equilibration is roughly two orders of magnitude faster. At cryogenic temperatures, intermonomer transfer limits the lifetimes of the lowest exciton band. The lifetimes are found to increase between 20 ps in the center of this band up to 100 ps toward lower energies, which is in very good agreement with the estimates from hole burning data. Interestingly, exciton delocalization in the FMO monomers is found to slow down intermonomer energy transfer, at both physiological and cryogenic temperatures.

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.

Quantum biology revisited

Photosynthesis is a highly optimized process from which valuable lessons can be learned about the operating principles in nature. Its primary steps involve energy transport operating near theoretical quantum limits in efficiency. Recently, extensive research was motivated by the hypothesis that nature used quantum coherences to direct energy transfer. This body of work, a cornerstone for the field of quantum biology, rests on the interpretation of small-amplitude oscillations in two-dimensional electronic spectra of photosynthetic complexes. This Review discusses recent work reexamining these claims and demonstrates that interexciton coherences are too short lived to have any functional significance in photosynthetic energy transfer. Instead, the observed long-lived coherences originate from impulsively excited vibrations, generally observed in femtosecond spectroscopy. These efforts, collectively, lead to a more detailed understanding of the quantum aspects of dissipation. Nature, rather than trying to avoid dissipation, exploits it via engineering of exciton-bath interaction to create efficient energy flow.

FMOxFMO: Elucidating Excitonic Interactions in the Fenna-Matthews-Olson Complex with the Fragment Molecular Orbital Method

In order to study Förster resonance energy transfer (FRET), the fragment molecular orbital (FMO) method is extended to compute electronic couplings between local excitations via the excited state transition density model, enabling efficient calculations of nonlocal excitations in a large molecular system and overcoming the previous limitation of being able to compute only local excitations. The results of these simple but accurate models are validated against full quantum calculations without fragmentation. The developed method is applied to a very important photosynthetic pigment-protein complex, the Fenna-Matthews-Olson complex (FMOc), that is responsible for the energy transfer from a chlorosome to the reaction center in the green sulfur bacteria. Absorption and circular dichroism spectra of FMOc are simulated, and the role of the molecular environment on the excitations is revealed.

Multicolor Quantum Control for Suppressing Ground State Coherences in Two-Dimensional Electronic Spectroscopy

The measured multidimensional spectral response of different light harvesting complexes exhibits oscillatory features which suggest an underlying coherent energy transfer. However, making this inference rigorous is challenging due to the difficulty of isolating excited state coherences in highly congested spectra. In this work, we provide a coherent control scheme that suppresses ground state coherences, thus making rephasing spectra dominated by excited state coherences. We provide a benchmark for the scheme using a model dimeric system and numerically exact methods to analyze the spectral response. We argue that combining temporal and spectral control methods can facilitate a second generation of experiments that are tailored to extract desired information and thus significantly advance our understanding of complex open many-body structure and dynamics.

Macroorganisation and flexibility of thylakoid membranes

The light reactions of photosynthesis are hosted and regulated by the chloroplast thylakoid membrane (TM) — the central structural component of the photosynthetic apparatus of plants and algae. The two-dimensional and three-dimensional arrangement of the lipid–protein assemblies, aka macroorganisation, and its dynamic responses to the fluctuating physiological environment, aka flexibility, are the subject of this review. An emphasis is given on the information obtainable by spectroscopic approaches, especially circular dichroism (CD). We briefly summarise the current knowledge of the composition and three-dimensional architecture of the granal TMs in plants and the supramolecular organisation of Photosystem II and light-harvesting complex II therein. We next acquaint the non-specialist reader with the fundamentals of CD spectroscopy, recent advances such as anisotropic CD, and applications for studying the structure and macroorganisation of photosynthetic complexes and membranes. Special attention is given to the structural and functional flexibility of light-harvesting complex II in vitro as revealed by CD and fluorescence spectroscopy. We give an account of the dynamic changes in membrane macroorganisation associated with the light-adaptation of the photosynthetic apparatus and the regulation of the excitation energy flow by state transitions and non-photochemical quenching.

New homologues of Brassicaceae water-soluble chlorophyll proteins shed light on chlorophyll binding, spectral tuning, and molecular evolution

Type-II water-soluble chlorophyll (Chl) proteins (WSCPs) of Brassicaceae are promising models for understanding how protein sequence and structure affect Chl binding and spectral tuning in photosynthetic Chl-protein complexes. However, to date, their use has been limited by the small number of known WSCPs, which also limited understanding their physiological roles. To overcome these limitations, we performed a phylogenetic analysis to compile a more comprehensive and complete set of natural type-II WSCP homologues. The identified homologues were heterologously expressed in Escherichia coli, purified, tested for assembly with chlorophylls, and spectroscopically characterized. The analyses led to the discovery of previously unrecognized type-IIa and IIb subclass WSCPs, as well as of a new subclass that did not bind chlorophylls. Further analysis by ancestral sequence reconstruction yielded sequences of putative ancestors of the three subclasses, which were subsequently recombinantly expressed in E. coli, purified and characterized. Combining the phylogenetic and spectroscopic data with molecular structural information revealed distinct Chl-binding motifs, and identified residues critically impacting spectral tuning. The distinct Chl-binding properties of the WSCP archetypes suggest that the non-Chl-binding subclass evolved from a Chl-binding ancestor that most likely lost its Chl-binding capacity upon localization in the plant tissues with low Chl content. This dual evolutionary trajectory is consistent with WSCPs association with the Kunitz-type protease inhibitors superfamily, and indications of their inhibitory activity in response to various forms of stress in plants. These findings suggest new directions for exploring the physiological roles of WSCPs and the correlation, if any, between Chl-binding and protease inhibition functionality.

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.

Wavelength-Dependent Exciton-Vibrational Coupling in the Water-Soluble Chlorophyll Binding Protein Revealed by Multilevel Theory of Difference Fluorescence Line-Narrowing

One of the most powerful line-narrowing techniques used to unravel the homogeneous lineshapes of inhomogeneously broadened systems is difference fluorescence line-narrowing spectroscopy. When this spectroscopy was applied to multichromophoric systems so far, the spectra were analyzed by an effective two-level system approach, composed of the electronic ground state and the lowest exciton state. An effective Huang-Rhys factor was assigned for the coupling of this state to the vibrations. Here, we extend this approach by including a multilevel line shape theory, which takes into account the excitonic coupling between pigments and thereby the effect of the delocalization of the excited states explicitly. In this way, it becomes possible to extract the spectral density of the local exciton-vibrational coupling. The theory is applied to the recombinant water-soluble chlorophyll binding protein reconstituted with chlorophyll a or b and reveals a significant decrease of the Huang-Rhys factor of the local exciton-vibrational coupling with decreasing transition energy of the chlorophylls. This decrease could be due to the increase in steric interactions reducing the flexibility of the environment and red-shifting the site energy of the pigments.

Protein Vibration Effects on Primary Electron Transfer Dynamics in Rhodobacter sphaeroides Photosynthetic Reaction Center

Primary electron transfer (ET) in the chromophore subsystem in a bacterial reaction center (RC) is a unique process, and is coupled with the protein motion, which, like the ET, is caused by photoexcitation of these chromophores. ET is also coupled with dissipative processes, which are caused by interaction between chromophores and vibrations of its surrounding protein. We propose a new dynamics calculation method that accounts for both these effects of protein vibrations. Within this method, the photoinduced protein motion causes an addition of coherent component to the ET rate. We performed dynamics calculation using this method and parameters, which were determined from the ab initio wave functions of the chromophore subsystem and protein normal vibrational modes. We showed that it is this protein motion that causes oscillations in the time-dependencies of stimulated emission intensities and of absorption at 1020 nm. Moreover, the latter oscillations are related to the coherent component of the ET rate.

Construction of Multichromophoric Spectra from Monomer Data: Applications to Resonant Energy Transfer

We develop a model that establishes a quantitative link between the physical properties of molecular aggregates and their constituent building blocks. The relation is built on the coherent potential approximation, calibrated against exact results, and proven reliable for a wide range of parameters. It provides a practical method to compute spectra and transfer rates in multichromophoric systems from experimentally accessible monomer data. Applications to Förster energy transfer reveal optimal transfer rates as functions of both the system-bath coupling and intra-aggregate coherence.

Origin of non-conservative circular dichroism of the CP29 antenna complex of photosystem II

Structure-based explanation of the non-conservative circular dichroism spectrum of the photosynthetic antenna complex CP29 of photosystem II.

Fine Tuning of Chlorophyll Spectra by Protein-Induced Ring Deformation

The ability to tune the light-absorption properties of chlorophylls by their protein environment is the key to the robustness and high efficiency of photosynthetic light-harvesting proteins. Unfortunately, the intricacy of the natural complexes makes it very difficult to identify and isolate specific protein-pigment interactions that underlie the spectral-tuning mechanisms. Herein we identify and demonstrate the tuning mechanism of chlorophyll spectra in type II water-soluble chlorophyll binding proteins from Brassicaceae (WSCPs). By comparing the molecular structures of two natural WSCPs we correlate a shift in the chlorophyll red absorption band with deformation of its tetrapyrrole macrocycle that is induced by changing the position of a nearby tryptophan residue. We show by a set of reciprocal point mutations that this change accounts for up to 2/3 of the observed spectral shift between the two natural variants.

Hole-Burning Spectroscopy on Excitonically Coupled Pigments in Proteins: Theory Meets Experiment

A theory for the calculation of resonant and nonresonant hole-burning (HB) spectra of pigment–protein complexes is presented and applied to the water-soluble chlorophyll-binding protein (WSCP) from cauliflower. The theory is based on a non-Markovian line shape theory (Renger and MarcusJ. Chem. Phys.2002, 116, 9997) and includes exciton delocalization, vibrational sidebands, and lifetime broadening. An earlier approach by Reppert (J. Phys. Chem. Lett.2011, 2, 2716) is found to describe nonresonant HB spectra only. Here we present a theory that can be used for a quantitative description of HB data for both nonresonant and resonant burning conditions. We find that it is important to take into account the excess energy of the excitation in the HB process. Whereas excitation of the zero-phonon transition of the lowest exciton state, that is, resonant burning allows the protein to access only its conformational substates in the neighborhood of the preburn state, any higher excitation gives the protein full access to all conformations present in the original inhomogeneous ensemble. Application of the theory to recombinant WSCP from cauliflower, reconstituted with chlorophyll a or chlorophyll b, gives excellent agreement with experimental data by Pieper et al. (J. Phys. Chem. B2011, 115, 405321417356) and allows us to obtain an upper bound of the lifetime of the upper exciton state directly from the HB experiments in agreement with lifetimes measured recently in time domain 2D experiments by Alster et al. (J. Phys. Chem. B2014, 118, 352424627983).

Förster resonance energy transfer, absorption and emission spectra in multichromophoric systems. I. Full cumulant expansions and system-bath entanglement

We study the Förster resonant energy transfer rate, absorption and emission spectra in multichromophoric systems. The multichromophoric Förster theory (MCFT) is determined from an overlap integral of generalized matrices related to the donor's emission and acceptor's absorption spectra, which are obtained via a full 2nd-order cumulant expansion technique developed in this work. We calculate the spectra and MCFT rate for both localized and delocalized systems, and calibrate the analytical 2nd-order cumulant expansion with the exact stochastic path integral method. We present three essential findings: (i) The role of the initial entanglement between the donor and its bath is found to be crucial in both the emission spectrum and the MCFT rate. (ii) The absorption spectra obtained by the cumulant expansion method are nearly identical to the exact spectra for both localized and delocalized systems, even when the system-bath coupling is far from the perturbative regime. (iii) For the emission spectra, the cumulant expansion can give reliable results for localized systems, but fail to provide reliable spectra of the high-lying excited states of a delocalized system, when the system-bath coupling is large and the thermal energy is small. This paper also provides a simple golden-rule derivation of the MCFT, reviews existing methods, and motivates further developments in the subsequent papers.