Delocalization Effects in Chlorophyll Fluorescence: Nonperturbative Line Shape Analysis of a Vibronically Coupled Dimer

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

Much can be learned about molecular aggregates by modeling their fluorescence-type spectra. In this study, we systematically describe the accuracy of various methods for simulating fluorescence-type linear spectra in a dimer system with a complex system-environment interaction, which serves as a model for various molecular aggregates, including most photosynthetic light-harvesting complexes (LHCs). We consider the approximate full cumulant expansion (FCE), complex time-dependent Redfield (ctR), time-independent Redfield, and modified Redfield methods and calculate their accuracy as a function of the site energy gap and coupling, excitonic energy gap, and dipole factor (i.e., type of spectrum). We find that the FCE method is the most accurate method for couplings smaller than 300 cm-1 at 300 K, but this method fails for very strong couplings or low temperatures due to inaccurate modeling of the equilibrium initial state. The ctR method performs well for the calculation of fluorescence and linear anisotropy spectra but poorer for circularly polarized fluorescence spectra or for all spectra when the coupling is strong (∼100cm-1). The Redfield and modified Redfield methods generally perform much more poorly than the ctR and FCE methods-especially for small excitonic energy gaps and strong couplings. We show that accurate modeling of the Stokes shift is crucial and present a version of the ctR method that treats both the Stokes shift and initial state correctly for the parameter ranges in plant LHCs. Apart from the application to LHCs, our results will be useful for the spectral characterization and design of organic molecular aggregates.

The interplay of excitonic delocalization and vibrational localization in optical lineshapes: A variational polaron approach

The dynamics of molecular excitonic systems are complicated by a competition between electronic coupling (which drives delocalization) and vibrational-electronic (vibronic) interactions (which tend to encourage electronic localization). A particular challenge of molecular systems is that they typically possess a large number of independent vibrations, with frequencies often spanning the entire spectrum of relevant electronic energy gaps. Recent spectroscopic observations and numerical simulations on a water-soluble chlorophyll-binding protein (WSCP) reveal a transition between two regimes of vibronic behavior, a Redfield-like regime in which low-frequency vibrations respond to a delocalized excitonic state, and a Förster-like regime where high-frequency vibrations act as incoherent excitations on individual pigments. Although numerical simulations can reproduce these effects, there is a need for a simple, systematic theory that accurately describes the smooth transition between these two regimes in experimental spectra. Here we address this challenge by generalizing the variational polaron transform approach of [Bloemsma et al., Chem. Phys. 481, 250 (2016)] to include arbitrary bath densities for systems with or without symmetry. We benchmark this theory against both numerical matrix-diagonalization methods and experimental 77 K fluorescence spectra for two WSCP variants, obtaining quite satisfactory agreement in both cases. We apply this theory to offer an explanation for the large loss in apparent electronic coupling in the WSCP Q57K mutant and to examine the likely impact of the interplay between excitonic delocalization and vibrational localization on vibrational sideband shapes and apparent coupling strengths in high-resolution optical spectra for chlorophyll-protein complexes such as WSCP.

Controlling Vibronic Coupling in Chlorophyll Proteins: The Effects of Excitonic Delocalization and Vibrational Localization

Vibrational-electronic (vibronic) coupling plays a critical role in excitation energy transfer in molecular aggregates and pigment-protein complexes (PPCs). But the interplay between excitonic delocalization and vibronic interactions is complex, often leaving even qualitative questions as to what conceptual framework (e.g., Redfield versus Förster theory) should be used to interpret experimental results. To shed light on this issue, we report here on the interplay between excitonic delocalization and vibronic coupling in site-directed mutants of the water-soluble chlorophyll protein (WSCP), as reflected in 77 K fluorescence spectra. Experimentally, we find that in PPCs where excitonic delocalization is disrupted (either by mutagenesis or heterodimer formation), the relative intensity of the vibrational sideband (VSB) in fluorescence spectra is suppressed by up to 37% compared to that of the native protein. Numerical simulations reveal that this effect results from the localization of high-frequency vibrations in the coupled system; while excitonic delocalization suppresses the purely electronic transition due to H-aggregate-like dipole-dipole interference, high-frequency vibrations are unaffected, leading to a relative enhancement of the VSB. By comparing VSB intensities of PPCs both in the presence and absence of excitonic delocalization, we extract a set of "local" Huang-Rhys (HR) factors for Chl a in WSCP. More generally, our results suggest a significant role for geometric effects in controlling energy-transfer rates (which depend sensitively on absorption/fluorescence line shapes) in molecular aggregates and PPCs.

A Simple Expression for the Screening of Excitonic Couplings between Chlorophylls as Inferred for Photosystem I Trimers

The Coulomb coupling between transition densities of the pigments in photosynthetic pigment-protein complexes, termed excitonic coupling, is a key factor for the description of optical spectra and energy transfer. A challenging question is the quantification of the screening of the excitonic coupling by the optical polarizability of the environment. We use the equivalence between the sophisticated quantum chemical polarizable continuum (PCM) model and the simple electrostatic Poisson-TrEsp approach to analyze the distance and orientation dependence of the dielectric screening between chlorophylls in photosystem I trimers. On the basis of these calculations we find that the vacuum couplings Vmn(0) and the couplings in the dielectric medium Vmn=fmnVmn(0) are related by the empirical screening factor fmn=0.60+39.6θ(|κmn|-1.17)exp(-0.56Rmn/Å), where κmn is the usual orientational factor of the dipole-dipole coupling between the pigments, Rmn is the center-to-center distance, and the Heaviside-function θ(|κmn|-1.17) ensures that the exponential distance dependence only contributes for in-line type dipole geometries. We are confident that the present expression can be applied also to other pigment-protein complexes with chlorophyll or related pigments of similar shape. The variance between the Poisson-TrEsp and the approximate coupling values is found to decrease by a factor of 8 and 3-4 using the present expression, instead of an exponential distance dependent or constant screening factor, respectively, assumed previously in the literature.

Deciphering between enhanced light emission and absorption in multi-mode porphyrin cavity polariton samples

It remains unclear how the collective strong coupling of cavity-confined photons to the electronic transitions of molecular chromophore leverages the distinct properties of the polaritonic constituents for future technologies. In this study, we design, fabricate, and characterize multiple types of Fabry-Pérot (FP) mirco-resonators containing copper(II) tetraphenyl porphyrin (CuTPP) to show how cavity polariton formation affects radiative relaxation processes in the presence of substantial non-Condon vibronic coupling between two of this molecule's excited electronic states. Unlike the prototypical enhancement of Q state radiative relaxation of CuTPP in a FP resonator incapable of forming polaritons, we find the light emission processes in multimode cavity polariton samples become enhanced for cavity-exciton energy differences near those of vibrations known to mediate non-Condon vibronic coupling. We propose the value of this detuning is consistent with radiative relaxation of Herzberg-Teller polaritons into collective molecular states coupled to the cavity photon coherently. We contrast the feature stemming from light emission from the HT polariton state with those that occur due to polariton-enhanced light absorption. Our results demonstrate the landscape of molecular and photonic interactions enabled by cavity polariton formation using complex chromophores and how researchers can design resonators to leverage these interactions to characterize and control polaritonic properties.

PigmentHunter: A point-and-click application for automated chlorophyll-protein simulations

Chlorophyll proteins (CPs) are the workhorses of biological photosynthesis, working together to absorb solar energy, transfer it to chemically active reaction centers, and control the charge-separation process that drives its storage as chemical energy. Yet predicting CP optical and electronic properties remains a serious challenge, driven by the computational difficulty of treating large, electronically coupled molecular pigments embedded in a dynamically structured protein environment. To address this challenge, we introduce here an analysis tool called PigmentHunter, which automates the process of preparing CP structures for molecular dynamics (MD), running short MD simulations on the nanoHUB.org science gateway, and then using electrostatic and steric analysis routines to predict optical absorption, fluorescence, and circular dichroism spectra within a Frenkel exciton model. Inter-pigment couplings are evaluated using point-dipole or transition-charge coupling models, while site energies can be estimated using both electrostatic and ring-deformation approaches. The package is built in a Jupyter Notebook environment, with a point-and-click interface that can be used either to manually prepare individual structures or to batch-process many structures at once. We illustrate PigmentHunter's capabilities with example simulations on spectral line shapes in the light harvesting 2 complex, site energies in the Fenna-Matthews-Olson protein, and ring deformation in photosystems I and II.

Quantum phase synchronization via exciton-vibrational energy dissipation sustains long-lived coherence in photosynthetic antennas

The lifetime of electronic coherences found in photosynthetic antennas is known to be too short to match the energy transfer time, rendering the coherent energy transfer mechanism inactive. Exciton-vibrational coherence time in excitonic dimers which consist of two chromophores coupled by excitation transfer interaction, can however be much longer. Uncovering the mechanism for sustained coherences in a noisy biological environment is challenging, requiring the use of simpler model systems as proxies. Here, via two-dimensional electronic spectroscopy experiments, we present compelling evidence for longer exciton-vibrational coherence time in the allophycocyanin trimer, containing excitonic dimers, compared to isolated pigments. This is attributed to the quantum phase synchronization of the resonant vibrational collective modes of the dimer, where the anti-symmetric modes, coupled to excitonic states with fast dephasing, are dissipated. The decoupled symmetric counterparts are subject to slower energy dissipation. The resonant modes have a predicted nearly 50% reduction in the vibrational amplitudes, and almost zero amplitude in the corresponding dynamical Stokes shift spectrum compared to the isolated pigments. Our findings provide insights into the mechanisms for protecting coherences against the noisy environment.

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.

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.

Bioexcitons by Design: How Do We Get There?

Biological pigment-protein complexes (PPCs) exhibit a remarkable ability to tune the optical properties of biological excitons (bioexcitons) through specific pigment-protein interactions. While such fine-tuning allows natural systems (e.g., photosynthetic proteins) to carry out their native functions with near-optimal performance, native function itself is often suboptimal for applications such as biofuel production or quantum technology development. This perspective offers a look at near-term prospects for the rational reoptimization of PPC bioexcitons for new functions using site-directed mutagenesis. The primary focus is on the "structure-spectrum" challenge of understanding the relationships between structural features and spectroscopic properties. While recent examples demonstrate that site-directed mutagenesis can be used to tune nearly all key bioexciton parameters (e.g., site energies, interpigment couplings, and electronic-vibrational interactions), critical challenges remain before we achieve truly rational design of bioexciton properties.

Light Emission from Vibronic Polaritons in Coupled Metalloporphyrin-Multimode Cavity Systems

In this study, we explore how one can use cavity polariton formation and a non-Condon vibronic coupling mechanism to form a type of hybrid light-matter state we denote as Herzberg-Teller (HT) vibronic polaritons. We use simple models to define the basic characteristics of these hybrid light-matter excitations including their dispersive energies. Experimentally, we find evidence of HT polaritons in the light emission spectra from copper(II) tetraphenylporphyrin (CuTPP) molecules strongly coupled to both single and multimode Fabry-Perot resonator structures. For specific resonator designs, we find evidence of significant enhancement of light emission from a short-lived sing-doublet state of CuTPP, which couples to a higher energy singlet state via a non-Condon vibronic mechanism. The results of a two-state model support the conclusion that this enhancement and the temperature-dependent dispersion of the light emission peak energy stem from radiative relaxation into cavity photon states dressed by collective vibrations of the molecules participating in polariton formation. These results show how researchers can leverage the complex interplay of electronic and nuclear degrees of freedom in light absorbing molecules to form a vaster array of coherent light-matter states and potentially transform platforms in optoelectronic and photocatalytic technologies.

Towards a quantitative description of excitonic couplings in photosynthetic pigment-protein complexes: quantum chemistry driven multiscale approaches

A structure-based quantitative calculation of excitonic couplings between photosynthetic pigments has to describe the dynamical polarization of the protein/solvent environment of the pigments, giving rise to reaction field and screening effects. Here, this challenging problem is approached by combining the fragment molecular orbital (FMO) method with the polarizable continuum model (PCM). The method is applied to compute excitonic couplings between chlorophyll a (Chl a) pigments of the water-soluble chlorophyll-binding protein (WSCP). By calibrating the vacuum dipole strength of the 0-0 transition of the Chl a chromophores according to experimental data, an excellent agreement between calculated and experimental linear absorption and circular dichroism spectra of WSCP is obtained. The effect of the mutual polarization of the pigment ground states is calculated to be very small. The simple Poisson-Transition-charge-from-Electrostatic-potential (Poisson-TrEsp) method is found to accurately describe the screening part of the excitonic coupling, obtained with FMO/PCM. Taking into account that the reaction field effects of the latter method can be described by a scalar constant leads to an improvement of Poisson-TrEsp that is expected to provide the basis for simple and realistic calculations of optical spectra and energy transfer in photosynthetic light-harvesting complexes. In addition, we present an expression for the estimation of Huang-Rhys factors of high-frequency pigment vibrations from experimental fluorescence line-narrowing spectra that takes into account the redistribution of oscillator strength by the interpigment excitonic coupling. Application to WSCP results in corrected Huang-Rhys factors that are less than one third of the original values obtained by the standard electronic two-state analysis that neglects the above redistribution. These factors are important for the estimation of the dipole strength of the 0-0 transition of the chromophores and for the development of calculation schemes for the spectral density of the exciton-vibrational coupling.

Accurate prediction of mutation-induced frequency shifts in chlorophyll proteins with a simple electrostatic model

Photosynthetic pigment-protein complexes control local chlorophyll (Chl) transition frequencies through a variety of electrostatic and steric forces. Site-directed mutations can modify this local spectroscopic tuning, providing critical insight into native photosynthetic functions and offering the tantalizing prospect of creating rationally designed Chl proteins with customized optical properties. Unfortunately, at present, no proven methods exist for reliably predicting mutation-induced frequency shifts in advance, limiting the method's utility for quantitative applications. Here, we address this challenge by constructing a series of point mutants in the water-soluble chlorophyll protein of Lepidium virginicum and using them to test the reliability of a simple computational protocol for mutation-induced site energy shifts. The protocol uses molecular dynamics to prepare mutant protein structures and the charge density coupling model of Adolphs et al. [Photosynth. Res. 95, 197-209 (2008)] for site energy prediction; a graphical interface that implements the protocol automatically is published online at http://nanohub.org/tools/pigmenthunter. With the exception of a single outlier (presumably due to unexpected structural changes), we find that the calculated frequency shifts match the experiment remarkably well, with an average error of 1.6 nm over a 9 nm spread in wavelengths. We anticipate that the accuracy of the method can be improved in the future with more advanced sampling of mutant protein structures.

Exciton Origin of Color-Tuning in Ca2+-Binding Photosynthetic Bacteria

Flexible color adaptation to available ecological niches is vital for the photosynthetic organisms to thrive. Hence, most purple bacteria living in the shade of green plants and algae apply bacteriochlorophyll a pigments to harvest near infra-red light around 850–875 nm. Exceptions are some Ca²⁺-containing species fit to utilize much redder quanta. The physical basis of such anomalous absorbance shift equivalent to ~5.5 kT at ambient temperature remains unsettled so far. Here, by applying several sophisticated spectroscopic techniques, we show that the Ca²⁺ ions bound to the structure of LH1 core light-harvesting pigment–protein complex significantly increase the couplings between the bacteriochlorophyll pigments. We thus establish the Ca-facilitated enhancement of exciton couplings as the main mechanism of the record spectral red-shift. The changes in specific interactions such as pigment–protein hydrogen bonding, although present, turned out to be secondary in this regard. Apart from solving the two-decade-old conundrum, these results complement the list of physical principles applicable for efficient spectral tuning of photo-sensitive molecular nano-systems, native or synthetic.