Pathways and timescales of primary charge separation in the photosystem II reaction center as revealed by a simultaneous fit of time-resolved fluorescence and transient absorption

Unraveling quantum coherences mediating primary charge transfer processes in photosystem II reaction center

Photosystem II (PSII) reaction center (RC) is a unique complex that is capable of efficiently separating electronic charges across the membrane. The primary energy- and charge-transfer (CT) processes occur on comparable ultrafast timescales, which makes it extremely challenging to understand the fundamental mechanism responsible for the near-unity quantum efficiency of the transfer. Here, we elucidate the role of quantum coherences in the ultrafast energy and CT in the PSII RC by performing two-dimensional (2D) electronic spectroscopy at the cryogenic temperature of 20 kelvin, which captures the distinct underlying quantum coherences. Specifically, we uncover the electronic and vibrational coherences along with their lifetimes during the primary ultrafast processes of energy and CT. We construct an excitonic model that provides evidence for coherent energy and CT at low temperature in the 2D electronic spectra. The principles could provide valuable guidelines for creating artificial photosystems with exploitation of system-bath coupling and control of coherences to optimize the photon conversion efficiency to specific functions.

Short-Range Effects in the Special Pair of Photosystem II Reaction Centers: The Nonconservative Nature of Circular Dichroism

Photosystem II reaction centers extract electrons from water, providing the basis of oxygenic life on earth. Among the light-sensitive pigments of the reaction center, a central chlorophyll a dimer, known as the special pair, so far has escaped a complete theoretical characterization of its excited state properties. The close proximity of the special pair pigments gives rise to short-range effects that comprise a coupling between local and charge transfer (CT) excited states as well as other intermolecular quantum effects. Using a multiscale simulation and a diabatization technique, we show that the coupling to CT states is responsible for 45% of the excitonic coupling in the special pair. The other short-range effects cause a nonconservative nature of the circular dichroism spectrum of the reaction center by effectively rotating the electric transition dipole moments of the special pair pigments inverting and strongly enhancing their intrinsic rotational strength.

A quantum trajectory picture of single photon absorption and energy transport in photosystem II

We use quantum trajectory theory to study the dynamics of the first step in photosynthesis for a single photon interacting with photosystem II (PSII). By considering individual trajectories we are able to look beyond the ensemble average dynamics to compute the PSII system evolution conditioned upon individual photon counting measurements. Measurements of the transmitted photon beam strongly affects the system state, since detection of an outgoing photon confirms that the PSII must be in the electronic ground state, while a null measurement implies it is in an excited electronic state. We show that under ideal conditions, observing the null result transforms a state with a low excited state population to a state with nearly all population contained in the excited states. We study the PSII dynamics conditioned on such photon counting for both a pure excitonic model of PSII and a more realistic model with exciton-phonon coupling to a dissipative phononic environment. In the absence of such coupling, we show that the measured fluorescence rates show oscillations constituting a photon-counting witness of excitonic coherence. Excitonic coupling to the phonon environment has a strong effect on the observed rates of fluorescence, damping the oscillations. Addition of non-radiative decay and incoherent transitions to radical pair states in the reaction center to the phononic model allows extraction of a quantum efficiency of 92.5% from the long-time evolution, consistent with bulk experimental measurements.

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.

Towards the description of charge transfer states in solubilised LHCII using subsystem DFT

Light harvesting complex II (LHCII) in plants and green algae have been shown to adapt their absorption properties, depending on the concentration of sunlight, switching between a light harvesting and a non-harvesting or quenched state. In a recent work, combining classical molecular dynamics (MD) simulations with quantum chemical calculations (Liguori et al. in Sci Rep 5:15661, 2015) on LHCII, it was shown that the Chl611-Chl612 cluster of the terminal emitter domain can play an important role in modifying the spectral properties of the complex. In that work the importance of charge transfer (CT) effects was highlighted, in re-shaping the absorption intensity of the chlorophyll dimer. Here in this work, we investigate the combined effect of the local excited (LE) and CT states in shaping the energy landscape of the chlorophyll dimer. Using subsystem Density Functional Theory over the classical [Formula: see text]s MD trajectory we look explicitly into the excitation energies of the LE and the CT states of the dimer and their corresponding couplings. Upon doing so, we observe a drop in the excitation energies of the CT states, accompanied by an increase in the couplings between the LE/LE and the LE/CT states facilitated by a shorter interchromophoric distance upon equilibration. Both these changes in conjunction, effectively produces a red-shift of the low-lying mixed exciton/CT states of the supramolecular chromophore pair.

Characterization of excited states in time-dependent density functional theory using localized molecular orbitals

Localized molecular orbitals are often used for the analysis of chemical bonds, but they can also serve to efficiently and comprehensibly compute linear response properties. While conventional canonical molecular orbitals provide an adequate basis for the treatment of excited states, a chemically meaningful identification of the different excited-state processes is difficult within such a delocalized orbital basis. In this work, starting from an initial set of supermolecular canonical molecular orbitals, we provide a simple one-step top-down embedding procedure for generating a set of orbitals, which are localized in terms of the supermolecule but delocalized over each subsystem composing the supermolecule. Using an orbital partitioning scheme based on such sets of localized orbitals, we further present a procedure for the construction of local excitations and charge-transfer states within the linear response framework of time-dependent density functional theory (TDDFT). This procedure provides direct access to approximate diabatic excitation energies and, under the Tamm-Dancoff approximation, also their corresponding electronic couplings-quantities that are of primary importance in modeling energy transfer processes in complex biological systems. Our approach is compared with a recently developed diabatization procedure based on subsystem TDDFT using projection operators, which leads to a similar set of working equations. Although both of these methods differ in the general localization strategies adopted and the type of basis functions (Slaters vs Gaussians) employed, an overall decent agreement is obtained.

Structure-based model of fucoxanthin-chlorophyll protein complex: Calculations of chlorophyll electronic couplings

Diatoms are a group of marine algae that are responsible for a significant part of global oxygen production. Adapted to life in an aqueous environment dominated by the blue-green light, their major light-harvesting antennae-fucoxanthin-chlorophyll protein complexes (FCPs)-exhibit different pigment compositions than of plants. Despite extensive experimental studies, until recently the theoretical description of excitation energy dynamics in these complexes was limited by the lack of high-resolution structural data. In this work, we use the recently resolved crystallographic information of the FCP complex from Phaeodactylum tricornutum diatom [Wang et al., Science 363, 6427 (2019)] and quantum chemistry-based calculations to evaluate the chlorophyll transition dipole moments, atomic transition charges from electrostatic potential, and the inter-chlorophyll couplings in this complex. The obtained structure-based excitonic couplings form the foundation for any modeling of stationary or time-resolved spectroscopic data. We also calculate the inter-pigment Förster energy transfer rates and identify two quickly equilibrating chlorophyll clusters.

Searching for a Unique Exciton Model of Photosynthetic Pigment–Protein Complexes: Photosystem II Reaction Center Study by Differential Evolution

Studying the optical properties of photosynthetic pigment–protein complexes (PPCs) in the visible light range, both experimentally and theoretically, is one of the ways of gaining knowledge about the function of the photosynthetic machinery of living species. To simulate the PPC optical response, it is necessary to use semiclassical theories describing the effect of external fields–matter interaction, energy migration in molecular crystals, and electron–phonon coupling. In this paper, we report the results of photosystem II reaction center (PSIIRC) linear optical response simulations. Applying the multimode Brownian oscillator model and the theory of molecular excitons, we have demonstrated that the absorption, circular and linear dichroism, and steady-state fluorescence of PSIIRC can be accurately fitted with the help of differential evolution (DE), the multiparametric evolutionary optimization algorithm. To explore the effectiveness of DE, we used the simulated experimental data as the target functions instead of those actually measured. Only 2 of 10 DE strategies have shown the best performance of the optimization algorithm. With the best tuning parameters of DE/rand-to-best/1/exp strategy determined from the strategy tests, we found the exact solution for the PSIIRC exciton model and fitted the spectra with a reasonable convergence rate.

Phthalocyanine as a Bioinspired Model for Chlorophyll f-Containing Photosystem II Drives Photosynthesis into the Far-Red Region

The textbook explanation that P₆₈₀ pigments are the red limit to drive oxygenic photosynthesis must be reconsidered by the recent discovery that chlorophyll f (Chlf)-containing Photosystem II (PSII) absorbing at 727 nm can drive water oxidation. Two different families of unsymmetrically substituted Zn phthalocyanines (Pc) absorbing in the 700-800 nm spectral window and containing a fused imidazole-phenyl substituent or a fused imidazole-hydroxyphenyl group have been synthetized and characterized as a bioinspired model of the Chlf/Tyrosine_Z /Histidine₁₉₀ cofactors of PSII. Transient absorption studies in the presence of an electron acceptor and irradiating in the far-red region evidenced an intramolecular electron transfer process. Visible and FT-IR signatures indicate the formation of a hydrogen-bonded phenoxyl radical in ZnPc II-OH. This study sets the foundation for the utilization of a broader spectral window for multi-electronic catalytic processes with one of the most robust and efficient dyes.

Two-dimensional electronic Stark spectroscopy of a photosynthetic dimer

Stark spectroscopy, which measures changes in the linear absorption of a sample in the presence of an external DC electric field, is a powerful experimental tool for probing the existence of charge-transfer (CT) states in photosynthetic systems. CT states often have small transition dipole moments, making them insensitive to other spectroscopic methods, but are particularly sensitive to Stark spectroscopy due to their large permanent dipole moment. In a previous study, we demonstrated a new experimental method, two-dimensional electronic Stark spectroscopy (2DESS), which combines two-dimensional electronic spectroscopy (2DES) and Stark spectroscopy. In order to understand how the presence of CT states manifest in 2DESS, here, we perform computational modeling and calculations of 2DESS as well as 2DES and Stark spectra, studying a photosynthetic dimer inspired by the photosystem II reaction center. We identify specific cases where qualitatively different sets of system parameters produce similar Stark and 2DES spectra but significantly different 2DESS spectra, showing the potential for 2DESS to aid in identifying CT states and their coupling to excitonic states.

New Perspectives on Photosystem II Reaction Centres

We apply the differential optical spectroscopy techniques of circular polarisation of luminescence (CPL) and magnetic CPL (MCPL) to the study of isolated reaction centres (RCs) of photosystem II (PS II). The data and subsequent analysis provide insights into aspects of the RC chromophore site energies, exciton couplings, and heterogeneities. CPL measurements are able to identify weak luminescence associated with the unbound chlorophyll-a (Chl-a) present in the sample. The overall sign and magnitude of the CPL observed relates well to the circular dichroism (CD) of the sample. Both CD and CPL are reasonably consistent with modelling of the RC exciton structure. The MCPL observed for the free Chl-a luminescence component in the RC samples is also easily understandable, but the MCPL seen near 680nm at 1.8K is anomalous, appearing to have a narrow, strongly negative component. A negative sign is inconsistent with MCPL of (exciton coupled) Qy states of either Chl-a or pheophytin-a (Pheo-a). We propose that this anomaly may arise as a result of the luminescence from a transient excited state species created following photo-induced charge separation within the RC. A comparison of CD spectra and modelling of RC preparations having a different number of pigments suggests that the non-conservative nature of the CD spectra observed is associated with the ‘special pair’ pigments PD1 and PD2.

Dissipation-Assisted Matrix Product Factorization

Charge and energy transfer in biological and synthetic organic materials are strongly influenced by the coupling of electronic states to a highly structured dissipative environment. Nonperturbative simulations of these systems require a substantial computational effort, and current methods can only be applied to large systems if environmental structures are severely coarse grained. Time evolution methods based on tensor networks are fundamentally limited by the times that can be reached due to the buildup of entanglement in time, which quickly increases the size of the tensor representation, i.e., the bond dimension. In this Letter, we introduce a dissipation-assisted matrix product factorization (DAMPF) method that combines a tensor network representation of the vibronic state within a pseudomode description of the environment where a continuous bosonic environment is mapped into a few harmonic oscillators under Lindblad damping. This framework is particularly suitable for a tensor network representation, since damping suppresses the entanglement growth among oscillators and significantly reduces the bond dimension required to achieve a desired accuracy. We show that dissipation removes the "time-wall" limitation of existing methods, enabling the long-time simulation of large vibronic systems consisting of 10-50 sites coupled to 100-1000 underdamped modes in total and for a wide range of parameter regimes. For these reasons, we believe that our formalism will facilitate the investigation of spatially extended systems with applications to quantum biology, organic photovoltaics, and quantum thermodynamics.

A theoretical study on the dynamics of light harvesting in the dimeric photosystem II core complex: regulation and robustness of energy transfer pathways

Here we present our theoretical investigations into the light reaction in the dimeric photosystem II (PSII) core complex. An effective model for excitation energy transfer (EET) and primary charge separation (CS) in the PSII core complex was developed, with model parameters constructed based on molecular dynamics (MD) simulation data. Compared to experimental results, we demonstrated that this model faithfully reproduces the absorption spectra of the RC and core light-harvesting complexes (CP43 and CP47) as well as the full EET dynamics among the chromophores in the PSII core complex. We then applied master equation simulations and network analysis to investigate detailed EET plus CS dynamics in the system, allowing us to identify key EET pathways and produce a coarse-grained cluster model for the light reaction in the dimeric PSII core complex. We show that non-equilibrium energy transfer channels play important roles in the efficient light harvesting process and that multiple EET pathways exist between subunits of PSII to ensure the robustness of light harvesting in the system. Furthermore, we revealed that inter-monomer energy transfer dominated by the coupling between the two CLA625 molecules enables efficient energy exchange between two CP47s in the dimeric PSII core complex, which leads to significant energy pooling in the CP47 domain during the light reaction. Our study provides a blueprint for the design of light harvesting in the PSII core and show that a structure-based approach using molecular dynamics simulations and quantum chemistry calculations can be effectively utilized to elucidate the dynamics of light harvesting in complex photosynthetic systems.

On the performance of a photosystem II reaction centre-based photocell

The photosystem II reaction centre is the photosynthetic complex responsible for oxygen production on Earth. Its water splitting function is particularly favoured by the formation of a stable charge separated state via a pathway that starts at an accessory chlorophyll. Here we envision a photovoltaic device that places one of these complexes between electrodes and investigate how the mean current and its fluctuations depend on the microscopic interactions underlying charge separation in the pathway considered. Our results indicate that coupling to well resolved vibrational modes does not necessarily offer an advantage in terms of power output but can lead to photo-currents with suppressed noise levels characterizing a multi-step ordered transport process. Besides giving insight into the suitability of these complexes for molecular-scale photovoltaics, our work suggests a new possible biological function for the vibrational environment of photosynthetic reaction centres, namely, to reduce the intrinsic current noise for regulatory processes.

Primary Charge Separation in the Photosystem II Reaction Center Revealed by a Global Analysis of the Two-dimensional Electronic Spectra

The transfer of electronic charge in the reaction center of Photosystem II is one of the key building blocks of the conversion of sunlight energy into chemical energy within the cascade of the photosynthetic reactions. Since the charge transfer dynamics is mixed with the energy transfer dynamics, an effective tool for the direct resolution of charge separation in the reaction center is still missing. Here, we use experimental two-dimensional optical photon echo spectroscopy in combination with the theoretical calculation to resolve its signature. A global fitting analysis allows us to clearly and directly identify a decay pathway associated to the primary charge separation. In particular, it can be distinguished from regular energy transfer and occurs on a time scale of 1.5 ps under ambient conditions. This technique provides a general tool to identify charge separation signatures from the energy transport in two-dimensional optical spectroscopy.

Quantum design of photosynthesis for bio-inspired solar-energy conversion

Photosynthesis is the natural process that converts solar photons into energy-rich products that are needed to drive the biochemistry of life. Two ultrafast processes form the basis of photosynthesis: excitation energy transfer and charge separation. Under optimal conditions, every photon that is absorbed is used by the photosynthetic organism. Fundamental quantum mechanics phenomena, including delocalization, underlie the speed, efficiency and directionality of the charge-separation process. At least four design principles are active in natural photosynthesis, and these can be applied practically to stimulate the development of bio-inspired, human-made energy conversion systems.

Electrostatic Asymmetry in the Reaction Center of Photosystem II

The exciton Hamiltonian of the chlorophyll (Chl) and pheophytin (Pheo) pigments in the reaction center (RC) of photosystem II is computed based on recent crystal structures by using the Poisson-Boltzmann/quantum-chemical method. Computed site energies largely confirm a previous model inferred from fits of optical spectra, in which Chl_D1 has the lowest site energy, while that of Pheo_D1 is higher than that of Pheo_D2. The latter assignment has been challenged recently under reference to mutagenesis experiments. We argue that these data are not in contradiction to our results. We conclude that Chl_D1 is the primary electron donor in both isolated RCs and intact core complexes at least at cryogenic temperatures. The main source of asymmetry in site energies is the charge distribution in the protein. Because many small contributions from various structural elements have to be taken into account, it can be assumed that this asymmetry was established in evolution by global optimization of the RC protein.

Efficiency of energy funneling in the photosystem II supercomplex of higher plants

The investigation of energy transfer properties in photosynthetic multi-protein networks gives insight into their underlying design principles. Here, we discuss the excitonic energy transfer mechanisms of the photosystem II (PS-II) C2S2M2 supercomplex, which is the largest isolated functional unit of the photosynthetic apparatus of higher plants. Despite the lack of a definite energy gradient in C2S2M2, we show that the energy transfer is directed by relaxation to low energy states. C2S2M2 is not organized to form pathways with strict energetically downhill transfer, which has direct consequences for the transfer efficiency, transfer pathways and transfer limiting steps. The exciton dynamics is sensitive to small changes in the energetic layout which, for instance, are induced by the reorganization of vibrational coordinates. In order to incorporate the reorganization process in our numerical simulations, we go beyond rate equations and use the hierarchically coupled equation of motion approach (HEOM). While transfer from the peripheral antenna to the proteins in proximity to the reaction center occurs on a faster time scale, the final step of the energy transfer to the RC core is rather slow, and thus the limiting step in the transfer chain. Our findings suggest that the structure of the PS-II supercomplex guarantees photoprotection rather than optimized efficiency.

Towards an ab initio description of the optical spectra of light-harvesting antennae: application to the CP29 complex of photosystem II

Light-harvesting pigment-protein complexes (PPC) represent the fundamental units through which the photosynthetic organisms absorb sunlight and funnel the energy to the reaction centre for carrying out the primary energy conversion reactions of photosynthesis. Here we apply a multiscale computational strategy to a specific PPC present in the photosystem II of plants and algae (CP29) to investigate in what detail should the environment effects due to protein and membrane/solvent be included for an accurate description of optical spectra. We find that a refinement of the crystal structure is needed before any meaningful quantum chemical calculations of pigment transition energies can be performed. For this purpose we apply classical molecular dynamics simulations of the PPC within its natural environment and we perform ab initio computations of the exciton Hamiltonian of the complex, including the environment either implicitly by the polarizable continuum model (PCM) or explicitly using the polarizable QM/MM methodology (MMPol). However, PCM essentially leads to an unspecific redshift of all transition energies, and MMPol is able to reveal site-specific changes in the optical properties of the pigments. Based on the latter and the excitonic couplings obtained within a polarizable QM/MM methodology, optical spectra are calculated, which are in good qualitative agreement with experimental data. A weakness of the approach is however found in the overestimation of the fluctuations of the excitonic parameters of the pigments along the MD trajectory. An explanation for such a finding in terms of the limits of the force fields commonly used for protein cofactors is presented and discussed.

Role of coherent vibrations in energy transfer and conversion in photosynthetic pigment-protein complexes

Oscillatory features of two-dimensional spectra of photosynthetic pigment-protein complexes during few picoseconds after electronic excitations of chlorophylls in various pigment-proteins were recently related to the coherent nuclear vibrations. It has been also speculated that the vibrations may assist the excitonic energy transfer and charge separation, hence contributing to energy transport and energy conversion efficiency. Here, we consider three theoretical approaches usually used for characterization of the excitation dynamics and charge separation, namely Redfield, Förster, and Marcus model descriptions, regarding this question. We show that two out of the three mechanisms require explicit resonances of excitonic splittings and the nuclear vibration frequencies. However, the third one related to the electron transfer is in principle off resonant.

Optimal Energy Transfer in Light-Harvesting Systems

Photosynthesis is one of the most essential biological processes in which specialized pigment-protein complexes absorb solar photons, and with a remarkably high efficiency, guide the photo-induced excitation energy toward the reaction center to subsequently trigger its conversion to chemical energy. In this work, we review the principles of optimal energy transfer in various natural and artificial light harvesting systems. We begin by presenting the guiding principles for optimizing the energy transfer efficiency in systems connected to dissipative environments, with particular attention paid to the potential role of quantum coherence in light harvesting systems. We will comment briefly on photo-protective mechanisms in natural systems that ensure optimal functionality under varying ambient conditions. For completeness, we will also present an overview of the charge separation and electron transfer pathways in reaction centers. Finally, recent theoretical and experimental progress on excitation energy transfer, charge separation, and charge transport in artificial light harvesting systems is delineated, with organic solar cells taken as prime examples.

A comparison between plant photosystem I and photosystem II architecture and functioning

Oxygenic photosynthesis is indispensable both for the development and maintenance of life on earth by converting light energy into chemical energy and by producing molecular oxygen and consuming carbon dioxide. This latter process has been responsible for reducing the CO2 from its very high levels in the primitive atmosphere to the present low levels and thus reducing global temperatures to levels conducive to the development of life. Photosystem I and photosystem II are the two multi-protein complexes that contain the pigments necessary to harvest photons and use light energy to catalyse the primary photosynthetic endergonic reactions producing high energy compounds. Both photosystems are highly organised membrane supercomplexes composed of a core complex, containing the reaction centre where electron transport is initiated, and of a peripheral antenna system, which is important for light harvesting and photosynthetic activity regulation. If on the one hand both the chemical reactions catalysed by the two photosystems and their detailed structure are different, on the other hand they share many similarities. In this review we discuss and compare various aspects of the organisation, functioning and regulation of plant photosystems by comparing them for similarities and differences as obtained by structural, biochemical and spectroscopic investigations.

Explicit calculation of the excited electronic states of the photosystem II reaction centre

The excited states of the photosystem II reaction centre cofactors have been calculated as a single “supermolecule”. Charge transfer states are shown to be dependent on electrostatic environment.

Dynamic protein conformations preferentially drive energy transfer along the active chain of the photosystem II reaction centre

One longstanding puzzle concerning photosystem II, a core component of photosynthesis, is that only one of the two symmetric branches in its reaction centre is active in electron transfer. To investigate the effect of the photosystem II environment on the preferential selection of the energy transfer pathway (a prerequisite for electron transfer), we have constructed an exciton model via extensive molecular dynamics simulations and quantum mechanics/molecular mechanics calculations based on a recent X-ray structure. Our results suggest that it is essential to take into account an ensemble of protein conformations to accurately compute the site energies. We identify the cofactor CLA606 of active chain as the most probable site for the energy excitation. We further pinpoint a number of charged protein residues that collectively lower the CLA606 site energy. Our work provides insights into the understanding of molecular mechanisms of the core machinery of the green-plant photosynthesis.

Cofactor-mediated energy and electron transfer in photosystem II occurs preferentially through one branch of the reaction centre, despite there being a symmetric path available. Here, the authors use computational methods to determine the influence of protein conformation on this selectivity.

From localized excited States to excitons: changing of conceptions of primary photosynthetic processes in the twentieth century

A short description of two theories of the primary photosynthetic processes is given. Generally accepted in 1950s-1990s, the localized excited states theory has been changed to the modern exciton theory. Appearance of the new experimental data and the light-harvesting complex crystal structure are reasons why the exciton theory has become important. The bulk of data for the old theory and outstanding experiments that have been the driving force for a new theory are discussed in detail.

Simulations of the two-dimensional electronic spectroscopy of the photosystem II reaction center

We report simulations of the two-dimensional electronic spectroscopy of the Q(y) band of the D1-D2-Cyt b559 photosystem II reaction center at 77 K. We base the simulations on an existing Hamiltonian that was derived by simultaneous fitting to a wide range of linear spectroscopic measurements and described within modified Redfield theory. The model obtains reasonable agreement with most aspects of the two-dimensional spectra, including the overall peak shapes and excited state absorption features. It does not reproduce the rapid equilibration from high energy to low energy excitonic states evident by a strong cross-peak below the diagonal. We explore modifications to the model to incorporate new structural data and improve agreement with the two-dimensional spectra. We find that strengthening the system-bath coupling and lowering the degree of disorder significantly improves agreement with the cross-peak feature, while lessening agreement with the relative diagonal/antidiagonal width of the 2D spectra. We conclude that two-dimensional electronic spectroscopy provides a sensitive test of excitonic models of the photosystem II reaction center and discuss avenues for further refinement of such models.

Two-dimensional electronic spectroscopy using incoherent light: theoretical analysis

Electronic energy transfer in photosynthesis occurs over a range of time scales and under a variety of intermolecular coupling conditions. Recent work has shown that electronic coupling between chromophores can lead to coherent oscillations in two-dimensional electronic spectroscopy measurements of pigment-protein complexes measured with femtosecond laser pulses. A persistent issue in the field is to reconcile the results of measurements performed using femtosecond laser pulses with physiological illumination conditions. Noisy-light spectroscopy can begin to address this question. In this work we present the theoretical analysis of incoherent two-dimensional electronic spectroscopy, I((4)) 2D ES. Simulations reveal diagonal peaks, cross peaks, and coherent oscillations similar to those observed in femtosecond two-dimensional electronic spectroscopy experiments. The results also expose fundamental differences between the femtosecond-pulse and noisy-light techniques; the differences lead to new challenges and new opportunities.

Mixed exciton-charge-transfer states in photosystem II: Stark spectroscopy on site-directed mutants

We investigated the electronic structure of the photosystem II reaction center (PSII RC) in relation to the light-induced charge separation process using Stark spectroscopy on a series of site-directed PSII RC mutants from the cyanobacterium Synechocystis sp. PCC 6803. The site-directed mutations modify the protein environment of the cofactors involved in charge separation (P(D1), P(D2), Chl(D1), and Phe(D1)). The results demonstrate that at least two different exciton states are mixed with charge-transfer (CT) states, yielding exciton states with CT character: (P(D2)(δ)(+)P(D1)(δ)(-)Chl(D1)) (673 nm) and (Chl(D1)(δ)(+)Phe(D1)(δ)(-)) (681 nm) (where the subscript indicates the wavelength of the electronic transition). Moreover, the CT state P(D2)(+)P(D1)(-) acquires excited-state character due to its mixing with an exciton state, producing (P(D2)(+)P(D1)(-))(δ) (684 nm). We conclude that the states that initiate charge separation are mixed exciton-CT states, and that the degree of mixing between exciton and CT states determines the efficiency of charge separation. In addition, the results reveal that the pigment-protein interactions fine-tune the energy of the exciton and CT states, and hence the mixing between these states. This mixing ultimately controls the selection and efficiency of a specific charge separation pathway, and highlights the capacity of the protein environment to control the functionality of the PSII RC complex.

Probing Electron-Transfer Times in Photosynthetic Reaction Centers by Hole-Burning Spectroscopy

A brief discussion is presented of transient hole-burned (HB) spectra (and the information that they provide) obtained for isolated reaction centers (RCs) from wild-type (WT) Rhodobacter sphaeroides, RCs containing zinc-bacteriochlorophylls (Zn-BChls), and RCs of Photosystem II (PSII) from spinach and Chlamydomonas reinhardtii . The shape of the spectral density and the strength of electron-phonon coupling in bacterial RCs are discussed. We focus, however, on heterogeneity of isolated PS II RCs from spinach and, in particular, Chlamydomonas reinhardtii , site energies of active (electron acceptor) and inactive pheophytins, the nature of the primary electron donor(s), and the possibility of multiple charge-separation (CS) pathways in the isolated PSII RC. We conclude with comments on current efforts in HB spectroscopy in the area of photosynthesis and future directions in HB spectroscopy.

Primary electron donor(s) in isolated reaction center of photosystem II from Chlamydomonas reinhardtii

Isolated reaction centers (RCs) from wild-type Chlamydomonas (C.) reinhardtii of Photosystem II (PSII), at different levels of intactness, were studied to provide more insight into the nature of the charge-separation (CS) pathway(s). We argue that previously studied D1/D2/Cytb559 complexes (referred to as RC680), with ChlD1 serving as the primary electron donor, contain destabilized D1 and D2 polypeptides and, as a result, do not provide a representative model system for the intact RC within the PSII core. The shapes of nonresonant transient hole-burned (HB) spectra obtained for more intact RCs (referred to as RC684) are very similar to P(+)QA(-) - PQA absorbance difference and triplet minus singlet spectra measured in PSII core complexes from Synechocystis PCC 6803 [Schlodder et al. Philos. Trans. R. Soc. London, Ser. B2008, 363, 1197]. We show that in the RC684 complexes, both PD1 and ChlD1 may serve as primary electron donors, leading to two different charge separation pathways. Resonant HB spectra cannot distinguish the CS times corresponding to different paths, but it is likely that the zero-phonon holes (ZPHs) observed in the 680-685 nm region (corresponding to CS times of ∼1.4-4.4 ps) reveal the ChlD1 pathway; conversely, the observation of charge-transfer (CT) state(s) in RC684 (in the 686-695 nm range) and the absence of ZPHs at λB > 685 nm likely stem from the PD1 pathway, for which CS could be faster than 1 ps. This is consistent with the finding of Krausz et al. [Photochem. Photobiol. Sci.2005, 4, 744] that CS in intact PSII core complexes can be initiated at low temperatures with fairly long-wavelength excitation. The lack of a clear shift of HB spectra as a function of excitation wavelength within the red-tail of the absorption (i.e., 686-695 nm) and the absence of ZPHs suggest that the lowest-energy CT state is largely homogeneously broadened. On the other hand, in usually studied destabilized RCs, that is, RC680, for which CT states have never been experimentally observed, ChlD1 is the most likely electron donor.

Nonlinear optical absorption of photosynthetic pigment molecules in leaves

A mathematical formulation of the relationship between optical absorption coefficient of photosynthetic pigment molecules and light intensity was developed. It showed that physical parameters of photosynthetic pigment molecule (i.e., light absorption cross-section of photosynthetic pigment molecule, its average lifetime in the excited state, total photosynthetic pigment molecules, the statistical weight, or degeneracy of energy level of photosynthetic pigment molecules in the ground state and in the excited state) influenced on both the light absorption coefficient and effective light absorption cross-section of photosynthetic pigment molecules. Moreover, it also showed that both the light absorption coefficient and effective light absorption cross-section of photosynthetic pigment molecules were not constant, they decreased nonlinearly with light intensity increasing. The occupation numbers of photosynthetic pigment molecules in the excited states increased nonlinearly with light intensity increasing.

Modeling of multi-exciton transient absorption spectra of protochlorophyllide aggregates in aqueous solution

Protochlorophyllide (Pchlide) is a natural porphyrin, a precursor of chlorophyll, synthesized by plants for its photosynthetic apparatus. The pigment spontaneously forms aggregates when dissolved in neat water solution. We present here calculations of the transient absorption spectra and its comprising components (ground-state bleach, stimulated emission, and excited-state absorption) for a strongly excitonically coupled linear chain of four Pchlide chromophores, using exciton theory with phenomenological Gaussian line shapes and without energetic disorder. A refined multiexciton model that includes static disorder is applied to fit the experimental power-dependent transient absorption spectra of aqueous protochlorophyllide and the kinetics for delay times up to 20 ps after photoexcitation. We show that population up to the 4-exciton manifold is sufficient to explain the pronounced saturation of the bleaching and the shape changes in the instantaneous, t = 0.2 ps transient spectra when the pulse energy is increased from 10 to 430 nJ per pulse. The decay of the multiexciton manifold is relatively slow and is preceded by a spectroscopically distinct process. We suggest that the exciton states in the Pchlide aggregates are mixed with charge-transfer states (CTS) and that the population and repopulation of the CTS coupled to the exciton states explains the relatively slow decay of the multiexciton manifold. The relevance of our results to the optical properties and dynamics of natural photosynthetic complexes and the possible physical origin of CTS formation are discussed.

Energy-transfer and charge-separation pathways in the reaction center of photosystem II revealed by coherent two-dimensional optical spectroscopy

The excited state dynamics and relaxation of electrons and holes in the photosynthetic reaction center of photosystem II are simulated using a two-band tight-binding model. The dissipative exciton and charge carrier motions are calculated using a transport theory, which includes a strong coupling to a harmonic bath with experimentally determined spectral density, and reduces to the Redfield, the Förster, and the Marcus expressions in the proper parameter regimes. The simulated third order two-dimensional signals, generated in the directions -k(1)+k(2)+k(3), k(1)-k(2)+k(3), and k(1)+k(2)-k(3), clearly reveal the exciton migration and the charge-separation processes.

Multiple charge-separation pathways in photosystem II: modeling of transient absorption kinetics

We explain the transient absorption kinetics (E. Romero, I. H. M. van Stokkum, V. I. Novoderezhkin, J. P. Dekker, R. van Grondelle, Biochemistry 2010, 49, 4300) measured for isolated reaction centers of photosystem II at 77 K upon excitation of the primary donor band (680 nm). The excited-state dynamics is modeled on the basis of the exciton states of 6 cofactors coupled to 4 charge-transfer (CT) states. One CT state (corresponding to charge separation within the special pair) is supposed to be strongly coupled with the excited states, whereas the other radical pairs are supposed to be localized. Relaxation within the strongly coupled manifold and transfer to localized CT's are described by the modified Redfield and generalized Förster theories, respectively. A simultaneous and quantitative fit of the 680, 545, and 460 nm kinetics (corresponding to respectively the Q(y) transitions of the red-most cofactors, Q(x) transition of pheophytin, and pheophytin anion absorption) enables us to define the pathways and time scales of primary electron transfer. A consistent modeling of the data is only possible with a Scheme where charge separation occurs from both the accessory chlorophyll and from the special pair, giving rise to fast and slow components of the pheophytin anion formation, respectively.

Quantum oscillatory exciton migration in photosynthetic reaction centers

The harvesting of solar energy and its conversion to chemical energy is essential for all forms of life. The primary photon absorption, transport, and charge separation events, which trigger a chain of chemical reactions, take place in membrane-bound photosynthetic complexes. Whether quantum effects, stemming from entanglement of chromophores, persist in the energy transport at room temperature, despite the rapid decoherence effects caused by environment fluctuations, is under current active debate. If confirmed, these may explain the high efficiency of light harvesting and open up numerous applications to quantum computing and information processing. We present simulations of the photosynthetic reaction center of photosystem II that clearly establish oscillatory energy transport at room temperature originating from interference of quantum pathways. These signatures of quantum transport may be observed by two dimensional coherent optical spectroscopy.

Physical origins and models of energy transfer in photosynthetic light-harvesting

We perform a quantitative comparison of different energy transfer theories, i.e. modified Redfield, standard and generalized Förster theories, as well as combined Redfield-Förster approach. Physical limitations of these approaches are illustrated and critical values of the key parameters indicating their validity are found. We model at a quantitative level the spectra and dynamics in two photosynthetic antenna complexes: in phycoerythrin 545 from cryptophyte algae and in trimeric LHCII complex from higher plants. These two examples show how the structural organization determines a directed energy transfer and how equilibration within antenna subunits and migration between subunits are superimposed.

Theory of excitation energy transfer: from structure to function

This mini-review summarizes our current theoretical knowledge about excitation energy transfer in pigment-protein complexes. The challenge for theory lies in the complexity of these systems and in the fact that the pigment-pigment and the pigment-protein interactions are of equal magnitude. The first part of this review contains an introduction to the theory of light harvesting and to structure- based calculations of the parameters of the theory. The second part provides a discussion of the standard Förster and Redfield theories of excitation energy transfer, which are valid in the limit of weak and strong pigment-pigment coupling, respectively. Afterward, we provide a description of recent extensions of the standard theories and discuss challenging problems to be solved in the future.