The origin of unidirectional charge separation in photosynthetic reaction centers: nonadiabatic quantum dynamics of exciton and charge in pigment-protein complexes

All-Atom Photoinduced Charge Transfer Dynamics in Condensed Phase via Multistate Nonlinear-Response Instantaneous Marcus Theory

Photoinduced charge transfer (CT) in the condensed phase is an essential component in solar energy conversion, but it is challenging to simulate such a process on the all-atom level. The traditional Marcus theory has been utilized for obtaining CT rate constants between pairs of electronic states but cannot account for the nonequilibrium effects due to the initial nuclear preparation. The recently proposed instantaneous Marcus theory (IMT) and its nonlinear-response formulation allow for incorporating the nonequilibrium nuclear relaxation to electronic transition between two states after the photoexcitation from the equilibrium ground state and provide the time-dependent rate coefficient. In this work, we extend the nonlinear-response IMT method for treating photoinduced CT among general multiple electronic states and demonstrate it in the organic photovoltaic carotenoid-porphyrin-fullerene triad dissolved in explicit tetrahydrofuran solvent. All-atom molecular dynamics simulations were employed to obtain the time correlation functions of energy gaps, which were used to generate the IMT-required time-dependent averages and variances of the relevant energy gaps. Our calculations show that the multistate IMT could capture the significant nonequilibrium effects due to the initial nuclear state preparation, and this is corroborated by the substantial differences between the population dynamics predicted by the multistate IMT and the Marcus theory, where the Marcus theory underestimates the population transfer. The population dynamics by multistate IMT is also shown to have a better agreement with the all-atom nonadiabatic mapping dynamics than the Marcus theory does. Because the multistate nonlinear-response IMT is straightforward and cost-effective in implementation and accounts for the nonequilibrium nuclear effects, we believe this method offers a practical strategy for studying charge transfer dynamics in complex condensed-phase systems.

Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers

Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair PD1-PD2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state PD1--PD2+ stabilizes energetically below the excitonic state. This underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

Observation of Coherent Symmetry-Breaking Vibration by Polarization-Dependent Femtosecond Spectroscopy

Understanding photoinduced chemical reactions beyond the Born-Oppenheimer paradigm requires a comprehensive examination of vibronic interactions. Although femtosecond studies have unveiled the influence of vibrational modes strongly coupled to ultrafast intramolecular reactions in the excited state, they often lack direct observations of how vibrations modulate electronic properties due to the rapid disappearance of reactants. To address this gap, our research investigates the dynamics of photoexcited molecules that do not react. Specifically, we focus on the coherent librational motion of molecular transition dipole moments, discovering that the coherent libration primarily originates from symmetry-breaking components in vibronically excited vibrational modes. Symmetry breaking motion can significantly impact the excited-state dynamics of highly symmetric molecules, potentially leading to nonadiabatic transitions. In essence, the data analysis framework introduced in this study can be harnessed to uncover potential reactivity in photoexcited molecules, further enhancing our understanding of the mechanisms governing these reactions.

Energy transfer from two luteins to chlorophylls in light-harvesting complex II study by using exciton models with phase correction

In light-harvesting complex II of plants, the two lutein pigments (LUT1 and LUT2) are always paired and an energy transfer pathway between them is believed to exist. However, it remains unclear whether this pathway is essential for the energy transfer between carotenoids and chlorophylls. In this work, we performed hybrid quantum mechanics/molecular mechanics simulations with Frenkel exciton models to investigate this energy transfer. The results show that the energy transfer pathways between the S2 state of LUT1 and CLAs are not affected by LUT2 S2. The energy transfer between LUT and chlorophyll-a (CLA) also follows a resonance mechanism. The two LUTs have different energy transfer pathways according to their energy gaps and coupling strengths with each CLA. The present work sheds light on the energy transfer pathways involved in the two LUTs.

Reaction Center Excitation in Photosystem II: From Multiscale Modeling to Functional Principles

Oxygenic photosynthesis is the fundamental energy-converting process that utilizes sunlight to generate molecular oxygen and the organic compounds that sustain life. Protein-pigment complexes harvest light and transfer excitation energy to specialized pigment assemblies, reaction centers (RC), where electron transfer cascades are initiated. A molecular-level understanding of the primary events is indispensable for elucidating the principles of natural photosynthesis and enabling development of bioinspired technologies. The primary enzyme in oxygenic photosynthesis is Photosystem II (PSII), a membrane-embedded multisubunit complex, that catalyzes the light-driven oxidation of water. The RC of PSII consists of four chlorophyll a and two pheophytin a pigments symmetrically arranged along two core polypeptides; only one branch participates in electron transfer. Despite decades of research, fundamental questions remain, including the origin of this functional asymmetry, the nature of primary charge-transfer states and the identity of the initial electron donor, the origin of the capability of PSII to enact charge separation with far-red photons, i.e., beyond the "red limit" where individual chlorophylls absorb, and the role of protein conformational dynamics in modulating charge-separation pathways.In this Account, we highlight developments in quantum-chemistry based excited-state computations for multipigment assemblies and the refinement of protocols for computing protein-induced electrochromic shifts and charge-transfer excitations calibrated with modern local correlation coupled cluster methods. We emphasize the importance of multiscale atomistic quantum-mechanics/molecular-mechanics and large-scale molecular dynamics simulations, which enabled direct and accurate modeling of primary processes in RC excitation at the quantum mechanical level.Our findings show how differential protein electrostatics enable spectral tuning of RC pigments and generate functional asymmetry in PSII. A chlorophyll pigment on the active branch (ChlD1) has the lowest site energy in PSII and is the primary electron donor. The complete absence of low-lying charge-transfer states within the central pair of chlorophylls excludes a long-held assumption about the initial charge separation. Instead, we identify two primary charge separation pathways, both with the same pheophytin acceptor (PheoD1): a fast pathway with ChlD1 as the primary electron donor (short-range charge-separation) and a slow pathway with PD1PD2 as the initial donor (long-range charge separation). The low-energy spectrum is dominated by two states with significant charge-transfer character, ChlD1δ+PheoD1δ- and PD1δ+PheoD1δ-. The conformational dynamics of PSII allows these charge-transfer states to span wide energy ranges, pushing oxygenic photosynthesis beyond the "red limit". These results provide a quantum mechanical picture of the primary events in the RC of oxygenic photosynthesis, forming a solid basis for interpreting experimental observations and for extending photosynthesis research in new directions.

Triplet states in the reaction center of Photosystem II

In oxygenic photosynthesis sunlight is harvested and funneled as excitation energy into the reaction center (RC) of Photosystem II (PSII), the site of primary charge separation that initiates the photosynthetic electron transfer chain. The chlorophyll ChlD1 pigment of the RC is the primary electron donor, forming a charge-separated radical pair with the vicinal pheophytin PheoD1 (ChlD1+PheoD1-). To avert charge recombination, the electron is further transferred to plastoquinone QA, whereas the hole relaxes to a central pair of chlorophylls (PD1PD2), subsequently driving water oxidation. Spin-triplet states can form within the RC when forward electron transfer is inhibited or back reactions are favored. This can lead to formation of singlet dioxygen, with potential deleterious effects. Here we investigate the nature and properties of triplet states within the PSII RC using a multiscale quantum-mechanics/molecular-mechanics (QM/MM) approach. The low-energy spectrum of excited singlet and triplet states, of both local and charge-transfer nature, is compared using range-separated time-dependent density functional theory (TD-DFT). We further compute electron paramagnetic resonance properties (zero-field splitting parameters and hyperfine coupling constants) of relaxed triplet states and compare them with available experimental data. Moreover, the electrostatic modulation of excited state energetics and redox properties of RC pigments by the semiquinone QA- is described. The results provide a detailed electronic-level understanding of triplet states within the PSII RC and form a refined basis for discussing primary and secondary electron transfer, charge recombination pathways, and possible photoprotection mechanisms in PSII.

Understanding Primary Charge Separation in the Heliobacterial Reaction Center

The homodimeric reaction center of heliobacteria retains features of the ancestral reaction center and can thus provide insights into the evolution of photosynthesis. Primary charge separation is expected to proceed in a two-step mechanism along either of the two reaction center branches. We reveal the first charge-separation step from first-principles calculations based on time-dependent density functional theory with an optimally tuned range-separated hybrid and ab initio Born-Oppenheimer molecular dynamics: the electron is most likely localized on the electron transfer cofactor 3 (EC3, OH-chlorophyll a), and the hole on the adjacent EC2. Including substantial parts of the surrounding protein environment into the calculations shows that a distinct structural mechanism is decisive for the relative energetic positioning of the electronic excitations: specific charged amino acids in the vicinity of EC3 lower the energy of charge-transfer excitations and thus facilitate efficient charge separation. These results are discussed considering recent experimental insights.

Alternative Fast and Slow Primary Charge-Separation Pathways in Photosystem II

Photosystem-II (PSII) is a multi-subunit protein complex that harvests sunlight to perform oxygenic photosynthesis. Initial light-activated charge separation takes place at a reaction centre consisting of four chlorophylls and two pheophytins. Understanding the processes following light excitation remains elusive due to spectral congestion, the ultrafast nature, and multi-component behaviour of the charge-separation process. Here, using advanced computational multiscale approaches which take into account the large-scale configurational flexibility of the system, we identify two possible primary pathways to radical-pair formation that differ by three orders of magnitude in their kinetics. The fast (short-range) pathway is dominant, but the existence of an alternative slow (long-range) charge-separation pathway hints at the evolution of redundancy that may serve other purposes, adaptive or protective, related to formation of the unique oxidative species that drives water oxidation in PSII.

Quantum mechanical analysis of excitation energy transfer couplings in photosystem II

We evaluated excitation energy transfer (EET) coupling (J) between all pairs of chlorophylls (Chls) and pheophytins (Pheos) in the protein environment of photosystem II based on the time-dependent density functional theory with a quantum mechanical/molecular mechanics approach. In the reaction center, the EET coupling between Chls PD1 and PD2 is weaker (|J(PD1/PD2)| = 79 cm-1), irrespective of a short edge-to-edge distance of 3.6 Å (Mg-to-Mg distance of 8.1 Å), than the couplings between PD1 and the accessory ChlD1 (|J(PD1/ChlD2)| = 104 cm-1) and between PD2 and ChlD2 (|J(PD2/ChlD1)| = 101 cm-1), suggesting that PD1 and PD2 are two monomeric Chls rather than a "special pair". There exist strongly coupled Chl pairs (|J| > ∼100 cm-1) in the CP47 and CP43 core antennas, which may be candidates for the red-shifted Chls observed in spectroscopic studies. In CP47 and CP43, Chls ligated to CP47-His26 and CP43-His56, which are located in the middle layer of the thylakoid membrane, play a role in the "hub" that mediates the EET from the lumenal to stromal layers. In the stromal layer, Chls ligated to CP47-His466, CP43-His441, and CP43-His444 mediate the EET from CP47 to ChlD2/PheoD2 and from CP43 to ChlD1/PheoD1 in the reaction center. Thus, the excitation energy from both CP47 and CP43 can always be utilized for the charge-separation reaction in the reaction center.

Electron-Transfer Route in the Early Oxidation States of the Mn4CaO5 Cluster in Photosystem II

The electron transfer from the oxygen-evolving Mn4CaO5 cluster to the electron acceptor D1-Tyr161 (TyrZ) is a prerequisite for water oxidation and O2 evolution. Here, we analyzed the electronic coupling in the rate-limiting electron-transfer transitions using a combined quantum mechanical/molecular mechanical/polarizable continuum model approach. In the S0 to S1 transition, the electronic coupling between the electron-donor Mn3(III) and TyrZ is small (2 meV). In contrast, the electronic coupling between the dangling Mn4(III) and TyrZ is significantly large (172 meV), which suggests that the electron transfer proceeds from Mn3(III) to TyrZ via Mn4(III). In the S1 to S2 transition, the electronic coupling between Mn4(III) and TyrZ is also larger (124 meV) than that between Mn1(III) and TyrZ (1 meV), which favors the formation of the open-cubane S2 conformation with Mn4(IV) over the formation of the closed-cubane S2 conformation with Mn1(IV). In the S0 to S1 and S1 to S2 transitions, the Mn4 d-orbital and the TyrZ π-orbital are hybridized via D1-Asp170, which suggests that D1-Asp170 commonly provides a dominant electron-transfer route.

Electron Transfer Route between Quinones in Type-II Reaction Centers

In photosynthetic reaction centers from purple bacteria (PbRCs) and photosystem II (PSII), the photoinduced charge separation is terminated by an electron transfer between the primary (QA) and secondary (QB) quinones. Here, we investigate the electron transfer route, calculating the superexchange coupling (HQA-QB) for electron transfer from QA to QB in the protein environment. HQA-QB is significantly larger in PbRC than in PSII. In superexchange electron tunneling, the electron transfer via unoccupied molecular orbitals of the nonheme Fe complex (QA → Fe → QB) is pronounced in PbRC, whereas the electron transfer via occupied molecular orbitals (Fe → QB followed by QA → Fe) is pronounced in PSII. The significantly large HQA-QB is caused by a water molecule that donates the H-bond to the ligand Glu-M234 in PbRC. The corresponding water molecule is absent in PSII due to the existence of D1-Tyr246. HQA-QB increases in response to the Ser-L223···QB H-bond formation caused by an extension of the H-bond network, which facilitates charge delocalization over the QB site. This explains the observed discrepancy in the QA-to-QB electron transfer between PbRC and PSII, despite their structural similarity.

Quasiparticle Self-Consistent GW-Bethe-Salpeter Equation Calculations for Large Chromophoric Systems

The GW-Bethe–Salpeter equation (BSE) method is promising for calculating the low-lying excitonic states of molecular systems. However, so far it has only been applied to rather small molecules and in the commonly implemented diagonal approximations to the electronic self-energy, it depends on a mean-field starting point. We describe here an implementation of the self-consistent and starting-point-independent quasiparticle self-consistent (qsGW)-BSE approach, which is suitable for calculations on large molecules. We herein show that eigenvalue-only self-consistency can lead to an unfaithful description of some excitonic states for chlorophyll dimers while the qsGW-BSE vertical excitation energies (VEEs) are in excellent agreement with spectroscopic experiments for chlorophyll monomers and dimers measured in the gas phase. Furthermore, VEEs from time-dependent density functional theory calculations tend to disagree with experimental values and using different range-separated hybrid (RSH) kernels does change the VEEs by up to 0.5 eV. We use the new qsGW-BSE implementation to calculate the lowest excitation energies of the six chromophores of the photosystem II (PSII) reaction center (RC) with nearly 2000 correlated electrons. Using more than 11,000 (6000) basis functions, the calculation could be completed in less than 5 (2) days on a single modern compute node. In agreement with previous TD-DFT calculations using RSH kernels on models that also do not include environmental effects, our qsGW-BSE calculations only yield states with local characters in the low-energy spectrum of the hexameric complex. Earlier works with RSH kernels have demonstrated that the protein environment facilitates the experimentally observed interchromophoric charge transfer. Therefore, future research will need to combine correlation effects beyond TD-DFT with an explicit treatment of environmental electrostatics.

Excimer evolution hampers symmetry-broken charge-separated states

Achieving long-lived symmetry-broken charge-separated states in chromophoric assemblies is quintessential for enhanced performance of artificial photosynthetic mimics. However, the occurrence of energy trap states hinders exciton and charge transport across photovoltaic devices, diminishing power conversion efficiency. Herein, we demonstrate unprecedented excimer formation in the relaxed excited-state geometry of bichromophoric systems impeding the lifetime of symmetry-broken charge-separated states. Core-annulated perylenediimide dimers (SC-SPDI2 and SC-NPDI2) prefer a near-orthogonal arrangement in the ground state and a π-stacked foldamer structure in the excited state. The prospect of an excimer-like state in the foldameric arrangement of SC-SPDI2 and SC-NPDI2 has been rationalized by fragment-based excited state analysis and temperature-dependent photoluminescence measurements. Effective electronic coupling matrix elements in the Franck-Condon geometry of SC-SPDI2 and SC-NPDI2 facilitate solvation-assisted ultrafast symmetry-breaking charge-separation (SB-CS) in a high dielectric environment, in contrast to unrelaxed excimer formation (Ex*) in a low dielectric environment. Subsequently, the SB-CS state dissociates into an undesired relaxed excimer state (Ex) due to configuration mixing of a Frenkel exciton (FE) and charge-separated state in the foldamer structure, downgrading the efficacy of the charge-separated state. The decay rate constant of the FE to SB-CS (k FE→SB-CS) in polar solvents is 8-17 fold faster than that of direct Ex* formation (k FE→Ex*) in non-polar solvent (k FE→SB-CS≫k FE→Ex*), characterized by femtosecond transient absorption (fsTA) spectroscopy. The present investigation establishes the impact of detrimental excimer formation on the persistence of the SB-CS state in chromophoric dimers and offers the requisite of conformational rigidity as one of the potential design principles for developing advanced molecular photovoltaics.

The initial charge separation step in oxygenic photosynthesis

Photosystem II is crucial for life on Earth as it provides oxygen as a result of photoinduced electron transfer and water splitting reactions. The excited state dynamics of the photosystem II-reaction center (PSII-RC) has been a matter of vivid debate because the absorption spectra of the embedded chromophores significantly overlap and hence it is extremely difficult to distinguish transients. Here, we report the two-dimensional electronic-vibrational spectroscopic study of the PSII-RC. The simultaneous resolution along both the visible excitation and infrared detection axis is crucial in allowing for the character of the excitonic states and interplay between them to be clearly distinguished. In particular, this work demonstrates that the mixed exciton-charge transfer state, previously proposed to be responsible for the far-red light operation of photosynthesis, is characterized by the Chl_D1⁺Phe radical pair and can be directly prepared upon photoexcitation. Further, we find that the initial electron acceptor in the PSII-RC is Phe, rather than P_D1, regardless of excitation wavelength.

The photosystem II reaction center (PSII-RC) is a model system to understand the initial steps of photosynthesis, but its excited state dynamics is difficult to disentangle with most spectroscopic methods. Here the authors perform a two-dimensional electronic-vibrational spectroscopic study of PSII-RC, providing detailed insight into such dynamics and into the mechanism of charge separation.

The Electronic Origin of Far-Red-Light-Driven Oxygenic Photosynthesis

Photosystem‐II uses sunlight to trigger charge separation and catalyze water oxidation. Intrinsic properties of chlorophyll a pigments define a natural “red limit” of photosynthesis at ≈680 nm. Nevertheless, charge separation can be triggered with far‐red photons up to 800 nm, without altering the nature of light‐harvesting pigments. Here we identify the electronic origin of this remarkable phenomenon using quantum chemical and multiscale simulations on a native Photosystem‐II model. We find that the reaction center is preorganized for charge separation in the far‐red region by specific chlorophyll–pheophytin pairs, potentially bypassing the light‐harvesting apparatus. Charge transfer can occur along two distinct pathways with one and the same pheophytin acceptor (Pheo_D1). The identity of the donor chlorophyll (Chl_D1 or P_D1) is wavelength‐dependent and conformational dynamics broaden the sampling of the far‐red region by the two charge‐transfer states. The two pathways rationalize spectroscopic observations and underpin designed extensions of the photosynthetically active radiation limit.

Energetic Basis and Design of Enzyme Function Demonstrated Using GFP, an Excited-State Enzyme

The past decades have witnessed an explosion of de novo protein designs with a remarkable range of scaffolds. It remains challenging, however, to design catalytic functions that are competitive with naturally occurring counterparts as well as biomimetic or nonbiological catalysts. Although directed evolution often offers efficient solutions, the fitness landscape remains opaque. Green fluorescent protein (GFP), which has revolutionized biological imaging and assays, is one of the most redesigned proteins. While not an enzyme in the conventional sense, GFPs feature competing excited-state decay pathways with the same steric and electrostatic origins as conventional ground-state catalysts, and they exert exquisite control over multiple reaction outcomes through the same principles. Thus, GFP is an "excited-state enzyme". Herein we show that rationally designed mutants and hybrids that contain environmental mutations and substituted chromophores provide the basis for a quantitative model and prediction that describes the influence of sterics and electrostatics on excited-state catalysis of GFPs. As both perturbations can selectively bias photoisomerization pathways, GFPs with fluorescence quantum yields (FQYs) and photoswitching characteristics tailored for specific applications could be predicted and then demonstrated. The underlying energetic landscape, readily accessible via spectroscopy for GFPs, offers an important missing link in the design of protein function that is generalizable to catalyst design.

Long-Range Electron Tunneling from the Primary to Secondary Quinones in Photosystem II Enhanced by Hydrogen Bonds with a Nonheme Fe Complex

The mechanisms governing the long-range electron tunneling from the primary (Q_A) to secondary (Q_B) quinones in photosystem II are clarified by analyzing superexchange pathways through a nonheme Fe complex, using a quantum mechanics/molecular mechanics/polarizable continuum model approach. The electron tunneling rate is evaluated using the Marcus-Levich-Jortner theory considering electronic coupling, energy difference, and Franck-Condon factor. The superexchange Q_A → Q_B electron tunneling is enhanced by hybridized σ/σ* orbitals of histidines (D2-His214 and D1-His215) via penetration of the wave function into hydrogen bonds with both Q_A and Q_B. Despite a large energy gap to the intermediate states, the contributions of the histidine σ/σ* orbitals to the superexchange coupling are larger than those of π/π* orbitals. Fe²⁺ is not an essential component for the Q_A → Q_B electron tunneling because hybridized histidine molecular orbitals can be coupled with both Q_A and Q_B simultaneously in the absence of Fe d orbitals.

Electrostatic profiling of photosynthetic pigments: implications for directed spectral tuning

Photosynthetic pigment-protein complexes harvest solar energy with a high quantum efficiency. Protein scaffolds are known to tune the spectral properties of embedded pigments principally through structured electrostatic environments. Although the physical nature of electrostatic tuning is straightforward, the precise spatial principles of electrostatic preorganization remain poorly explored for different protein matrices and incompletely characterized with respect to the intrinsic properties of different photosynthetic pigments. In this work, we study the electronic structure features associated with the lowest excited state of a series of eight naturally occurring (bacterio)chlorophylls and pheophytins to describe the precise topological differences in electrostatic potentials and hence determine intrinsic differences in the expected mode and impact of electrostatic tuning. The difference electrostatic potentials between the ground and first excited states are used as fingerprints. Both the spatial profile and the propensity for spectral tuning are found to be unique for each pigment, indicating spatially and directionally distinct modes of electrostatic tuning. The results define a specific partitioning of the protein matrix around each pigment as an aid to identify regions with a maximal impact on spectral tuning and have direct implications for dimensionality reduction in protein design and engineering. Thus, a quantum mechanical basis is provided for understanding, predicting, and ultimately designing sequence-modified or pigment-exchanged biological systems, as suggested for selected examples of pigment-reconstituted proteins.

Excitation Energy Transfer between bodipy Dyes in a Symmetric Molecular Excitonic Seesaw

We examine the redistribution of energy between electronic and vibrational degrees of freedom that takes place between a π-conjugated oligomer, a phenylene-butadiynylene, and two identical boron-dipyrromethene (bodipy) end-caps using femtosecond transient absorption spectroscopy, single-molecule spectroscopy, and nonadiabatic excited-state molecular dynamics (NEXMD) modeling techniques. The molecular structure represents an excitonic seesaw in that the excitation energy on the oligomer backbone can migrate to either one end-cap or the other, but not to both. The NEXMD simulations closely reproduce the characteristic time scale for redistribution of electronic and vibrational energy of 2.2 ps and uncover the vibrational modes contributing to the intramolecular relaxation. The calculations indicate that the dihedral angle between the bodipy dye and the oligomer change upon excitation of the oligomer. Single-molecule experiments reveal a difference in photoluminescence lifetime of the bodipy dyes depending on whether they are excited by direct absorption or by redistribution of energy from the backbone. This difference in lifetime may be attributed to the difference in dihedral angle. The simulations also suggest that a strong coupling can occur between the two end-caps, giving rise to a reversible shuttling of excitation energy between them. Strong coupling should lead to a pronounced loss in polarization memory of the fluorescence since the oligomer backbone tends to be slightly distorted and the two bodipy transition dipoles have different orientations. A sensitive single-molecule technique is presented to test for such coupling. However, although redistribution of electronic and vibrational energy between the end-caps can occur, it appears to be unidirectional and irreversible, suggesting that an additional localization mechanism is at play which is, as yet, not fully accounted for in the simulations.