Coarse-Grained Approach to Simulate Signatures of Excitation Energy Transfer in Two-Dimensional Electronic Spectroscopy of Large Molecular Systems.

Two-dimensional electronic spectroscopy (2DES) has proven to be a highly effective technique in studying the properties of excited states and the process of excitation energy transfer in complex molecular assemblies, particularly in biological light-harvesting systems. However, the accurate simulation of 2DES for large systems still poses a challenge because of the heavy computational demands it entails. In an effort to overcome this limitation, we devised a coarse-grained 2DES method. This method encompasses the treatment of the entire system by dividing it into distinct weakly coupled segments, which are assumed to communicate predominantly through incoherent exciton transfer. We first demonstrate the efficiency of this method through simulation on a model dimer system, which demonstrates a marked improvement in calculation efficiency, with results that exhibit good concordance with reference spectra calculated with less approximate methods. Additionally, the application of this method to the light-harvesting antenna 2 (LH2) complex of purple bacteria showcases its advantages, accuracy, and limitations. Furthermore, simulating the anisotropy decay in LH2 induced by energy transfer and its comparison with experiments confirm that the method is capable of accurately describing dynamical processes in a biologically relevant system. This method presented lends itself to an extension that accounts for the effect of intrasegment relaxation processes on the 2DES spectra, which for computational efficiency are ignored in the implementation reported here. It is envisioned that the method will be employed in the future to accurately and efficiently calculate 2D spectra of more extensive systems, such as photosynthetic supercomplexes.

Ultrafast Energy Transfer in a Diatom Photosystem II Supercomplex.

The light-harvesting complexes (LHCs) of diatoms, specifically fucoxanthin-Chl a/c binding proteins (FCPs), exhibit structural and functional diversity, as highlighted by recent structural studies of photosystem II-FCP (PSII-FCPII) supercomplexes from different diatom species. The excitation dynamics of PSII-FCPII supercomplexes isolated from the diatom Thalassiosira pseudonana was explored using time-resolved fluorescence spectroscopy and two-dimensional electronic spectroscopy at room temperature and 77 K. Energy transfer between FCPII and PSII occurred remarkably fast (<5 ps), emphasizing the efficiency of FCPII as a light-harvesting antenna. The presence of long-wavelength chlorophylls may further help concentrate excitations in the core complex and increase the efficiency of light harvesting. Structure-based calculations reveal remarkably strong excitonic couplings between chlorophylls in the FCP antenna and between FCP and the PSII core antenna that are the basis for the rapid energy transfer.

Inhomogeneous energy transfer dynamics from iron-stress-induced protein A to photosystem I.

Cyanobacteria respond to iron limitation by producing the pigment-protein complex IsiA, forming rings associated with photosystem I (PSI). Initially considered a chlorophyll-storage protein, IsiA is known to act as an auxiliary light-harvesting antenna of PSI, increasing its absorption cross-section and reducing the need for iron-rich PSI core complexes. Spectroscopic studies have demonstrated efficient energy transfer from IsiA to PSI. Here we investigate the room-temperature excitation dynamics in isolated PSI-IsiA, PSI, IsiA monomer complexes and IsiA aggregates using two-dimensional electronic spectroscopy. Cross analyses of the data from these three samples allow us to resolve components of energy transfer between IsiA and PSI with lifetimes of 2-3 ps and around 20 ps. Structure-based Förster theory calculations predict a single major timescale of IsiA-PSI equilibration, that depends on multiple energy transfer routes between different IsiA subunits in the ring. Despite the experimentally observed lifetime heterogeneity, which is attributed to structural heterogeneity of the supercomplexes, IsiA is found to be a unique, highly efficient, membrane antenna complex in cyanobacteria.

Inter-subunit energy transfer processes in a minimal plant photosystem II supercomplex.

Photosystem II (PSII) is an integral part of the photosynthesis machinery, in which several light-harvesting complexes rely on inter-complex excitonic energy transfer (EET) processes to channel energy to the reaction center. In this paper, we report on a direct observation of the inter-complex EET in a minimal PSII supercomplex from plants, containing the trimeric light-harvesting complex II (LHCII), the monomeric light-harvesting complex CP26, and the monomeric PSII core complex. Using two-dimensional (2D) electronic spectroscopy, we measure an inter-complex EET timescale of 50 picoseconds for excitations from the LHCII-CP26 peripheral antenna to the PSII core. The 2D electronic spectra also reveal that the transfer timescale is nearly constant over the pump spectrum of 600 to 700 nanometers. Structure-based calculations reveal the contribution of each antenna complex to the measured inter-complex EET time. These results provide a step in elucidating the full inter-complex energy transfer network of the PSII machinery.

Measuring the Ultrafast Correlation Dynamics of a Multilevel System Using the Center Line Slope Analysis in Two-Dimensional Electronic Spectroscopy.

In a two-dimensional (2D) optical spectrum of a multilevel system, there are diagonal peaks and off-diagonal cross-peaks that correlate the different levels. The time-dependent properties of these diagonal peaks and cross-peaks contain much information about the dynamics of the multilevel system. The time-dependent diagonal peakshape that depends on the spectral diffusion dynamics of the associated transition and characterized by the frequency-fluctuation correlation function (FFCF) is well studied. However, the time-dependent peakshape of a cross-peak that provides the correlation dynamics between different transitions is much less studied or understood. We derived the third-order nonlinear response functions that describe the cross-peaks in a 2D electronic spectrum of a multilevel system that arise from processes sharing a common ground state and/or from internal conversion and population transfer. We can use the center line slope (CLS) analysis to characterize the cross-peaks in conjunction with the diagonal peaks. This allows us to recover the frequency-fluctuation cross-correlation functions (FXCFs) between two transitions. The FXCF and its subsidiary quantities such as the initial correlation and the initial covariance between different transitions are important for studying the correlation effects between states in complex systems, such as energy-transfer processes. Furthermore, knowledge of how various molecular processes over different timescales affect simultaneously different transitions can also be obtained from the measured FXCF. We validated and tested our derived equations and analysis process by studying, as an example, the 2D electronic spectra of metal-free phthalocyanine in solution. We measured and analyzed the diagonal peaks of the Qx and Qy transitions and the cross-peaks between these two transitions of this multilevel electronic system and obtained the associated FFCFs and FXCFs. In this model system, we measured negative components of FXCF over the tens of picosecond timescale. This suggests that in phthalocyanine, the Qx and Qy transitions coupling with the solvent molecule motion are anticorrelated to each other.

An efficient time-domain implementation of the multichromophoric Förster resonant energy transfer method.

The excitation energy transfer (EET) process for photosynthetic antenna complexes consisting of subunits, each comprised of multiple chromophores, remains challenging to describe. The multichromophoric Förster resonance energy transfer theory is a popular method to describe the EET process in such systems. This paper presents a new time-domain method for calculating energy transfer based on the combination of multichromophoric Förster resonance energy transfer theory and the Numerical Integration of the Schrödinger Equation method. After validating the method on simple model systems, we apply it to the Light-Harvesting antenna 2 (LH2) complex, a light harvesting antenna found in purple bacteria. We use a simple model combining the overdamped Brownian oscillators to describe the dynamic disorder originating from the environmental fluctuations and the transition charge from the electrostatic potential coupling model to determine the interactions between chromophores. We demonstrate that with this model, both the calculated spectra and the EET rates between the two rings within the LH2 complex agree well with experimental results. We further find that the transfer between the strongly coupled rings of neighboring LH2 complexes can also be well described with our method. We conclude that our new method accurately describes the EET rate for biologically relevant multichromophoric systems, which are similar to the LH2 complex. Computationally, the new method is very tractable, especially for slow processes. We foresee that the method can be applied to efficiently calculate transfer in artificial systems as well and may pave the way for calculating multidimensional spectra of extensive multichromophoric systems in the future.

Measuring the Ultrafast Spectral Diffusion and Vibronic Coupling Dynamics in CdSe Colloidal Quantum Wells using Two-Dimensional Electronic Spectroscopy.

We measure the ultrafast spectral diffusion, vibronic dynamics, and energy relaxation of a CdSe colloidal quantum wells (CQWs) system at room temperature using two-dimensional electronic spectroscopy (2DES). The energy relaxation of light-hole (LH) excitons and hot carriers to heavy-hole (HH) excitons is resolved with a time scale of ∼210 fs. We observe the equilibration dynamics between the spectroscopically accessible HH excitonic state and a dark state with a time scale of ∼160 fs. We use the center line slope analysis to quantify the spectral diffusion dynamics in HH excitons, which contains an apparent sub-200 fs decay together with oscillatory features resolved at 4 and 25 meV. These observations can be explained by the coupling to various lattice phonon modes. We further perform quantum calculations that can replicate and explain the observed dynamics. The 4 meV mode is observed to be in the near-critically damped regime and may be mediating the transition between the bright and dark HH excitons. These findings show that 2DES can provide a comprehensive and detailed characterization of the ultrafast spectral properties in CQWs and similar nanomaterials.

Ultrafast Excitation Energy Transfer Dynamics in the LHCII-CP29-CP24 Subdomain of Plant Photosystem II.

We measure the two-dimensional electronic spectra of the LHCII(M)-CP29-CP24 complex in photosystem II (PSII) and provide the first study of the ultrafast excitation energy transfer (EET) processes of an asymmetric and native light-harvesting assembly of the antenna of PSII. With comparisons to LHCII, we observe faster energy equilibrations in the intermediate levels of the LHCII(M)-CP29-CP24 complex at 662 and 670 nm. Notably, the putative "bottleneck" states in LHCII exhibit faster effective dynamics in the LHCII(M)-CP24-CP29 complex, with the average lifetime shortening from 2.5 ps in LHCII to 1.2 ps in the bigger assembly. The observations are supported by high-level structure-based calculations, and the accelerated dynamics can be attributed to the structural change of LHCII(M) in the bigger complex. This study shows that the biological functioning structures of the complexes are important to understand the overall EET dynamics of the PSII supercomplex.

Ultrafast excitation quenching by the oxidized photosystem II reaction center.

Photosystem II (PSII) is the pigment-protein complex driving the photoinduced oxidation of water and reduction of plastoquinone in all oxygenic photosynthetic organisms. Excitations in the antenna chlorophylls are photochemically trapped in the reaction center (RC) producing the chlorophyll-pheophytin radical ion pair P+ Pheo-. When electron donation from water is inhibited, the oxidized RC chlorophyll P+ acts as an excitation quencher, but knowledge on the kinetics of quenching is limited. Here, we used femtosecond transient absorption spectroscopy to compare the excitation dynamics of PSII with neutral and oxidized RC (P+). We find that equilibration in the core antenna has a major lifetime of about 300 fs, irrespective of the RC redox state. Two-dimensional electronic spectroscopy revealed additional slower energy equilibration occurring on timescales of 3-5 ps, concurrent with excitation trapping. The kinetics of PSII with open RC can be described well with previously proposed models according to which the radical pair P+ Pheo- is populated with a main lifetime of about 40 ps, which is primarily determined by energy transfer between the core antenna and the RC chlorophylls. Yet, in PSII with oxidized RC (P+), fast excitation quenching was observed with decay lifetimes as short as 3 ps and an average decay lifetime of about 90 ps, which is shorter than the excited-state lifetime of PSII with open RC. The underlying mechanism of this extremely fast quenching prompts further investigation.

Two-dimensional electronic spectroscopy of the Qx to Qy relaxation of chlorophylls a in photosystem II core complex.

Using two-dimensional electronic spectroscopy, we measured the Qx to Qy transfer dynamics of the chlorophyll a (Chl a) manifold in the photosystem II (PSII) monomeric core complex from Arabidopsis thaliana. A PSII monomeric core consists of 35 Chls a and no Chl b, thus allowing for a clear window to study Chl a Qx dynamics in a large pigment-protein complex. Initial excitation in the Qx band results in a transfer to the Qy band in less than 60 fs. Upon the ultrafast transfer, regardless of the excitation frequency within the Qx band, the quasi-transient absorption spectra are very similar. This observation indicates that Chl a's Qx to Qy transfer is not frequency selective. Using a simple model, we determined that this is not due to the lifetime broadening of the ultrafast transfer but predominantly due to a lack of correlation between the PSII core complex's Chl a Qx and Qy bands. We suggest the origin to be the intrinsic loss of correlation during the Qx to Qy internal conversion as observed in previous studies of molecular Chl a dissolved in solvents.

Two-Dimensional Electronic Spectroscopy of a Minimal Photosystem I Complex Reveals the Rate of Primary Charge Separation.

Photosystem I (PSI), found in all oxygenic photosynthetic organisms, uses solar energy to drive electron transport with nearly 100% quantum efficiency, thanks to fast energy transfer among antenna chlorophylls and charge separation in the reaction center. There is no complete consensus regarding the kinetics of the elementary steps involved in the overall trapping, especially the rate of primary charge separation. In this work, we employed two-dimensional coherent electronic spectroscopy to follow the dynamics of energy and electron transfer in a monomeric PSI complex from Synechocystis PCC 6803, containing only subunits A-E, K, and M, at 77 K. We also determined the structure of the complex to 4.3 Å resolution by cryoelectron microscopy with refinements to 2.5 Å. We applied structure-based modeling using a combined Redfield-Förster theory to compute the excitation dynamics. The absorptive 2D electronic spectra revealed fast excitonic/vibronic relaxation on time scales of 50-100 fs from the high-energy side of the absorption spectrum. Antenna excitations were funneled within 1 ps to a small pool of chlorophylls absorbing around 687 nm, thereafter decaying with 4-20 ps lifetimes, independently of excitation wavelength. Redfield-Förster energy transfer computations showed that the kinetics is limited by transfer from these red-shifted pigments. The rate of primary charge separation, upon direct excitation of the reaction center, was determined to be 1.2-1.5 ps-1. This result implies activationless electron transfer in PSI.

Electron Transfer Quenching of Rhodamine 6G by N-Methylpyrrole Is an Unproductive Process in the Photocatalytic Heterobiaryl Cross-Coupling Reaction.

In the heterobiaryl cross-coupling reaction between aryl halides (Ar-X) and N-methylpyrrole (N-MP) catalyzed by rhodamine 6G (Rh6G+) under irradiation with visible light, a highly active and long-lived (millisecond time range) rhodamine 6G radical (Rh6G•) is formed upon electron transfer from N,N-diisopropylethylamine (DIPEA) to Rh6G+. In this study, we utilized steady-state and time-resolved spectroscopy techniques to demonstrate the existence of another electron-transfer process occurring from the relatively electron-rich N-MP to photoexcited Rh6G+ that was neglected in the previous reports. In this case, the radical Rh6G• formed is short-lived and undergoes rapid recombination (nanosecond time-range), rendering it ineffective in reducing Ar-X to aryl radicals Ar• that can subsequently be trapped by N-MP. This is further demonstrated via two model reactions involving 4'-bromoacetophenone and 1,3,5-tribromobenzene with insignificant product yields after visible-light irradiation in the absence of DIPEA. The unproductive quenching of photoexcited Rh6G+ by N-MP leads to a lower concentration of photocatalyst available for competitive charge transfer with DIPEA and hence decreases the efficiency of the cross-coupling reaction.

Characterization of the ultrafast spectral diffusion and vibronic coherence of TIPS-pentacene using 2D electronic spectroscopy.

TIPS-pentacene is a small-molecule organic semiconductor that is widely used in optoelectronic devices. It has been studied intensely owing to its ability to undergo singlet fission. In this study, we aim to develop further understanding of the coupling between the electronic and nuclear degrees of freedom of TIPS-pentacene (TIPS-Pn). We measured and analyzed the 2D electronic spectra of TIPS-Pn in solutions. Using center line slope (CLS) analysis, we characterized the frequency-fluctuation correlation function of the 0-0 vibronic transition. Strong oscillations in the CLS values were observed for up to 5 ps with a frequency of 264 cm-1, which are attributable to a large vibronic coupling with the TIPS-Pn ring-breathing vibrational mode. In addition, detailed analysis of the CLS values allowed us to retrieve two spectral diffusion lifetimes, which are attributed to the inertial and diffusive dynamics of solvent molecules. Amplitude beating analysis also uncovered couplings with another vibrational mode at 1173 cm-1. The experimental results can be described using the displaced harmonic oscillator model. By comparing the CLS values of the simulated data with the experimental CLS values, we estimated a Huang-Rhys factor of 0.1 for the ring-breathing vibrational mode. The results demonstrated how CLS analysis can be a useful method for characterizing the strength of vibronic coupling.

An Exciton Dynamics Model of Bryopsis corticulans Light-Harvesting Complex II.

Bryopsis corticulans is a marine green macroalga adapted to the intertidal environment. It possesses siphonaxanthin-binding light-harvesting complexes of photosystem II (LHCII) with spectroscopic properties markedly different from the LHCII in plants. By applying a phenomenological fitting procedure to the two-dimensional electronic spectra of the LHCII from B. corticulans measured at 77 K, we can extract information about the excitonic states and energy-transfer processes. The fitting method results in well-converged parameters, including excitonic energy levels with their respective transition dipole moments, spectral widths, energy-transfer rates, and coupling properties. The 2D spectra simulated from the fitted parameters concur very well with the experimental data, showing the robustness of the fitting method. An excitonic energy-transfer scheme can be constructed from the fitting parameters. It shows the rapid energy transfer from chlorophylls (Chls) b to a at subpicosecond time scales and a long-lived state in the Chl b region at around 659 nm. Three weakly connected terminal states are resolved at 671, 675, and 677 nm. The lowest state is higher in energy than that in plant LHCII, which is probably because of the fewer number of Chls a in a B. corticulans LHCII monomer. Modeling based on existing Hamiltonians for the plant LHCII structure with two Chls a switched to Chls b suggests several possible Chl a-b replacements in comparison with those of plant LHCII. The adaptive changes result in a slower energy equilibration in the complex, revealed by the longer relaxation times of several exciton states compared to those of plant LHCII. The strength of our phenomenological fitting method for obtaining excitonic energy levels and energy-transfer network is put to the test in systems such as B. corticulans LHCII, where prior knowledge on exact assignment and spatial locations of pigments are lacking.

Observing the Fluctuation Dynamics of Dative Bonds Using Two-Dimensional Electronic Spectroscopy.

We perform two-dimensional electronic spectroscopy on chlorophyll (Chl) a and b molecules in aprotic solvents of different Lewis basicity. By analyzing the ultrafast spectral diffusion dynamics of the Qy transition, we show that a certain timescale of the spectral diffusion dynamics is affected by the solvents' Lewis basicity. Control experiments with Chlorin-e6-a Chl molecule analog-and ab initio time-dependent density functional theory calculations confirm that we are directly probing the fluctuation dynamics of the dative bond between the solvent's lone pair and the Mg2+ center in Chls that is responsible for the Lewis basicity. The observation is indicative of dative bond length and angular fluctuations with timescales ranging between ∼30 and 150 ps and the dative bond-strength-dependent perturbation on the Qy transition frequency of Chls.

Hierarchical Equations of Motion Simulation of Temperature-Dependent Two-Dimensional Electronic Spectroscopy of the Chlorophyll a Manifold in LHCII.

We simulated two-dimensional electronic spectra (2DES) of the chlorophyll a manifold of light-harvesting complex II (LHCII) at various temperatures (77, 110, 150, 190, 230, 273, and 293 K) using the hierarchical equations of the motion-phase matching approach. We confirm the main excitation energy transfer pathways assignments within the chlorophyll a manifold of LHCII measured in a recent work (J. Phys. Chem. B 2019, 123, 6765-6775). The calculated transfer rates are also in general agreement with the measured rates. We also provided theoretical confirmation for the experimental assignments, as uphill and downhill energy transfer processes, of 2D spectral features that were reported in recent experimental reports. These temperature-dependent features were also ascertained to follow the detailed-balance principle.

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.

Spectral tuning of light-harvesting complex II in the siphonous alga Bryopsis corticulans and its effect on energy transfer dynamics.

Light-harvesting complex II (LHCII) from the marine green macroalga Bryopsis corticulans is spectroscopically characterized to understand the structural and functional changes resulting from adaptation to intertidal environment. LHCII is homologous to its counterpart in land plants but has a different carotenoid and chlorophyll (Chl) composition. This is reflected in the steady-state absorption, fluorescence, linear dichroism, circular dichroism and anisotropic circular dichroism spectra. Time-resolved fluorescence and two-dimensional electronic spectroscopy were used to investigate the consequences of this adaptive change in the pigment composition on the excited-state dynamics. The complex contains additional Chl b spectral forms - absorbing at around 650 nm and 658 nm - and lacks the red-most Chl a forms compared with higher-plant LHCII. Similar to plant LHCII, energy transfer between Chls occurs on timescales from under hundred fs (mainly from Chl b to Chl a) to several picoseconds (mainly between Chl a pools). However, the presence of long-lived, weakly coupled Chl b and Chl a states leads to slower exciton equilibration in LHCII from B. corticulans. The finding demonstrates a trade-off between the enhanced absorption of blue-green light and the excitation migration time. However, the adaptive change does not result in a significant drop in the overall photochemical efficiency of Photosystem II. These results show that LHCII is a robust adaptable system whose spectral properties can be tuned to the environment for optimal light harvesting.

Insights into the mechanisms and dynamics of energy transfer in plant light-harvesting complexes from two-dimensional electronic spectroscopy.

During the past two decades, two-dimensional electronic spectroscopy (2DES) and related techniques have emerged as a potent experimental toolset to study the ultrafast elementary steps of photosynthesis. Apart from the highly engaging albeit controversial analysis of the role of quantum coherences in the photosynthetic processes, 2DES has been applied to resolve the dynamics and pathways of energy and electron transport in various light-harvesting antenna systems and reaction centres, providing unsurpassed level of detail. In this paper we discuss the main technical approaches and their applicability for solving specific problems in photosynthesis. We then recount applications of 2DES to study the exciton dynamics in plant and photosynthetic light-harvesting complexes, especially light-harvesting complex II (LHCII) and the fucoxanthin-chlorophyll proteins of diatoms, with emphasis on the types of unique information about such systems that 2DES is capable to deliver. This article is part of a Special Issue entitled Light harvesting, edited by Dr. Roberta Croce.

Revealing the excitation energy transfer network of Light-Harvesting Complex II by a phenomenological analysis of two-dimensional electronic spectra at 77 K.

Energy equilibration in light-harvesting antenna systems normally occurs before energy is transferred to a reaction center. The equilibration mechanism is a characteristic of the excitation energy transfer (EET) network of the antenna. Characterizing this network is crucial in understanding the first step of photosynthesis. We present our phenomenology-based analysis procedure and results in obtaining the excitonic energy levels, spectral linewidths, and transfer-rate matrix of Light-Harvesting Complex II directly from its 2D electronic spectra recorded at 77 K with waiting times between 100 fs to 100 ps. Due to the restriction of the models and complexity of the system, a unique EET network cannot be constructed. Nevertheless, a recurring pattern of energy transfer with very similar overall time scales between spectral components (excitons) is consistently obtained. The models identify a "bottleneck" state in the 664-668 nm region although with a relatively shorter lifetime (∼4-6 ps) of this state compared to previous studies. The model also determines three terminal exciton states at 675, 677-678, and 680-681 nm that are weakly coupled to each other. The excitation energy equilibration between the three termini is found to be independent of the initial excitation conditions, which is a crucial design for the light-harvesting complexes to ensure the energy flow under different light conditions and avoid excitation trapping. We proposed two EET schemes with tentative pigment assignments based on the interpretation of the modeling results together with previous structure-based calculations and spectroscopic observables.