Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

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

Richard Stones a, Hoda Hossein-Nejad a, Rienk van Grondelle b and Alexandra Olaya-Castro *a
aDepartment of Physics and Astronomy, University College London, Gower Street, London, WC1E 6BT, UK. E-mail:
bDepartment of Physics and Astronomy, VU University, 1081 HV Amsterdam, The Netherlands

Received 7th July 2017 , Accepted 4th August 2017

First published on 4th August 2017

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.

Life on Earth is fueled by photosynthesis, the process by which plants, algae and certain bacteria convert solar energy into stable chemical energy.1 The initial electron transfer steps during solar energy conversion by these organisms take place in photosynthetic reaction centres (PRCs), sophisticated trans-membrane supramolecular pigment–protein complexes that exhibit a dual device-like functionality. Under illumination, a PRC complex effectively operates as Nature's solar cell1 where electronic excitations of chromophores are transformed into stable charge-separated states, with electron donor and electron acceptor separated by a few nanometres. This picosecond charge separation process occurs with near unit quantum efficiency implying that almost every quanta of energy absorbed results in charge separated across the PRC.2,3 The same chromophore–protein structure and energetic landscape promoting this quantum yield also favours a diode-like behaviour of all PRCs such that under an appropriate applied bias, electric current flows almost entirely in one direction.4 Their functional versatility, nanometre size, and near-unit quantum efficiency has motivated the exploration of PRCs as possible components of photovoltaic and photoelectrochemical cells5,6 as well as in biomolecular electronics.4,7,8

A step further in this field is the recent development of single-molecule techniques that allow measurement of the photocurrent through individual PRC complexes.9 Using cysteine group mutations, it has been possible to bind a photosystem I unit to a gold substrate and use a scanning probe gold-tip that acts as both an electrode and a localized light source to excite and measure the photocurrent of a fully functional PRC.9 Moreover, using a tapping mode atomic force microscope, it has been possible to confirm that electrons tunneling through a bacterial PRC, under an applied bias, follow the transfer pathway of the photochemical charge separation.4 These experiments open up a new platform to carry out further investigations on how the microscopic mechanisms underlying the function of PRCs affect the electric current output of a single PRC, as well as to reveal further details of such microscopic mechanisms. In particular, it is foreseeable that besides measuring the current–voltage features these techniques may allow characterization of current fluctuations and the associated counting statistics of electron transport in PRCs. In quantum transport setups10–12 it has been shown that such fluctuations can reveal intrinsic dynamical features of the quantum system through which electron transport occurs, including the influence of electron–phonon interactions13 and coherence.14 A theoretical study along the lines of counting statistics for light-harvesting complexes has been carried out15 but so far there has been no investigation of full counting statistics of charge transport in PRCs.

The photosystem II reaction centre (PSIIRC), present in higher plants, algae and cyanobacteria,1 is arguably the most important prototype to be considered as it is responsible for water splitting and production of all oxygen on Earth.16 Experimental evidence indicates that one of the distinguishing features of PSIIRC is the existence of at least two different charge separation pathways, one of which starts in the monomeric chlorophyll of the active D1 branch (ChlD1).17–19 This transfer pathway has been argued to be a deciding factor for the functional operation of PSIIRC as a water splitting complex20 and is therefore the focus of this paper. While the detailed quantum mechanical features underpinning charge separation in PRCs are still under scrutiny,20 there is a wealth of steady-state and time-resolved spectroscopy revealing the electronic state space and spectral density of fluctuations relevant for the formation of stable charge separated states.19–22 However, the implications of these microscopic mechanisms for the electric current output of a single PSIIRC are largely unknown. In fact, it is unclear whether these natural light-to-charge converting units are well suited for anthropogenic use: would photocells integrating the microscopic mechanisms of PRCs deliver optimal power? How do these mechanisms affect the statistics of electron transport? The answers to these questions hold the potential to provide valuable insight both on the biophysics of these systems and for the development of the next generation of bio-inspired energy technologies.

In this work we address these questions by envisioning a photocell device where a single PSIIRC using the ChlD1 pathway is placed between two electrodes and investigate how the microscopic mechanisms underlying the photochemical charge separation affect the electric current output and its fluctuations, under continuous incoherent illumination. By comparing the photocell operation under different spectral densities characterizing the interaction between electronic and vibrational degrees of freedom, we show that selective coupling to underdamped vibrations does not necessarily offer an advantage in terms of current and power output for this photocell but they lead to output currents with suppressed noise levels as quantified by a Fano factor less than one. A structured spectral density that includes coupling to well-resolved vibrations allows the current noise strength to probe the structure of the exciton manifold which transfers population to charge transfer states and leads to a sub-Poissonian statistics. This indicates that both the exciton manifold and the electron–vibration interactions in PSIIRC support a multi-step ordered electron transport process. Our work therefore puts forward a new possible functional role for the vibrational environment of PRCs, that of guaranteeing intrinsic current noise control for regulatory purposes. Our analysis also gives insight into the suitability of PSIIRCs and their microscopic principles for photovoltaic or nano-electronic applications.


A photocell based on the photosystem II reaction centre

The prototype complex we consider is the PSIIRC for which crystallography has provided the arrangement of the chromophores involved in primary charge transfer.23 As illustrated in Fig. 1(a), two central chlorophylls PD1 and PD2 (known as the special pair) are flanked by the accessory chlorophylls ChlD1 and ChlD2, and the pheophytins PheD1 and PheD2 to form two branches labeled by the proteins D1 and D2. Charge separation only occurs down the D1 branch.24,25 Nonlinear spectroscopy has revealed at least two different excited states, (PD1PD2ChlD1)* and (ChlD1PheD1)*, that give rise to two different pathways (denoted the PD1 and ChlD1 pathways respectively) for charge separation along the D1 branch.19,21,26 The likelihood of each depends on the specific protein configuration and the corresponding disorder of pigment excitation energies.19 Our analysis focuses on the ChlD1 pathway for which electron transfer is predominantly initiated from the state (ChlD1PheD1)*26 given that in this protein configuration the ChlD1 pigment is lowest in energy. For this pathway there is not yet evidence of coherent dynamics between excited states and charge transfer (CT) states.19,21,26 The ChlD1 pathway can thus be seen as the relaxation dynamics in the energy landscape illustrated in Fig. 1(b).
image file: c7sc02983g-f1.tif
Fig. 1 Photosystem II reaction centre photocell model. (a) Schematic diagram of a proposed experimental setup for the photocell unit. The isolated core chromophores of PSIIRC are positioned between a gold substrate and a gold coated scanning probe microscope tip which act as electrodes. A silicon field-effect transistor (SFET) placed near the drain electrode could be used to measure the current statistics. (b) Energy level diagram showing the electronic state space of the model. The red arrows represent rates connecting the ground and empty state to the excited state manifold. ΓL and ΓR connect the system to the source and drain leads respectively while γex represents a coupling to an optical field which excites the system from the ground state to the lowest energy exciton state. Green arrows represent Förster/Marcus rates for primary and secondary charge separation. The load between states α and β indicates the transition across which we calculate the output current of the photocell and its statistics.

Fluorescence line narrowing experiments27 have given evidence of a highly structured spectral density characterizing the interactions of an excited chromophore in PSIIRC with a wide range of vibrational motions as shown in Fig. 2. Although some of the sharp modes in this spectral density have frequencies that match energy gaps between exciton states,28 we will show that this feature is not entirely determinant for the mean current output of a PSIIRC-based photocell operating within the ChlD1 pathway. Together with the state-space aforementioned, this spectral density of fluctuations has provided a good fit for steady-state and transient spectroscopy of primary charge separation in PSIIRC.21,26,28

image file: c7sc02983g-f2.tif
Fig. 2 Photosystem II reaction centre spectral densities. The components of the spectral densities used in the PSIIRC photocell model. Mode parameters of the structured component JM(ω) are shown in Table S1. The Drude part JD(ω) is scaled relative to the high energy parts for clarity. The inset shows the spectral densities J1(ω) and J2(ω) which are used to approximate the full spectral density.

We therefore consider these features of the ChlD1 charge separation pathway to put forward two key questions: (i) do these naturally occurring electron–vibration interactions and their associated spectral density provide the best strategy for maximizing current output in a bio-inspired photocell? and (ii) how do the exciton manifold and the vibrational environment affect the statistics of electrons flowing through this PRC?

These questions are addressed by envisioning a theoretical photocell device in which a single PSIIRC unit operating with the ChlD1 pathway is placed between two electrodes, as illustrated in Fig. 1(a). We assume uni-directional electron flow29,30 through the photocell such that electrons are transferred from the source lead to the system and from the system to the drain lead with phenomenological rates ΓL and ΓR respectively. The sample is incoherently photo-excited at rate γex. Our model enforces the Coulomb blockade regime such that the probability of two electron occupancy of the photocell is negligible.31 Assuming low enough excitation rates to guarantee that only single excitation states are populated, the state space of the PSIIRC-based photocell, shown in Fig. 1(b), spans the following: the ground state |g〉, six exciton states |X1〉 to |X6〉, the initial CT state |ChlD1+PheD1〉 ≡ |I〉, the secondary CT state |PD1+PheD1〉 ≡ |α〉 and the positively charged state |PD1+PheD1〉 ≡ |β〉 which represents the ‘empty’ state of the system for counting statistics calculations. The six exciton states arise from diagonalization of the site part of Hamiltonian Hel (see Methods) which includes coherent electronic interactions among all the six core chromophores located in both the D1 and D2 branches of the PSIIRC.

To investigate the effects of the vibrational environment in the performance of the photocell, we aim to compare four cases corresponding to the four different spectral densities depicted in Fig. 2: (i) the case where the full structured spectral density J(ω) = JD(ω) + JM(ω) is considered, the cases where, besides the low-energy background JD(ω), we account for one and two well-resolved modes with spectral densities (ii) J1(ω) = JD(ω) + G1(ω) and (iii) J2(ω) = J1(ω) + G2(ω), and (iv) the case where only the smooth low energy component JD(ω) is included. The expressions for the different components are given by:21

image file: c7sc02983g-t1.tif(1)
image file: c7sc02983g-t2.tif(2)

J D(ω) is the Drude form of a spectral density describing an overdamped Brownian oscillator where λD and ΩD are the reorganisation energy and cut off frequency, respectively. Gi(ω) describes the spectral density of an underdamped mode coupled to an excited pigment, with λj, ωj and γj being the reorganisation energy, frequency and damping rate of mode j respectively. image file: c7sc02983g-t3.tif has been measured experimentally27 and includes 48 underdamped modes. For case (ii) described by J1(ω), we consider ω1 = 342 cm−1 and for case (iii) corresponding to J2(ω), we consider ω1 as well as ω2 = 742 cm−1. These two modes have been argued to be important for electron transfer along the alternative PD1 charge separation pathway28 and in our case their frequencies span all the energy gaps of the system. Hence, J2(ω) provides a good approximation to the full spectral density. All parameters for these spectral densities are detailed in ESI Note 1. By comparing these cases we attempt to address the question of how well ‘adapted’ these electron–vibrational interactions are for photovoltaics: if it were possible to decouple these well resolved nuclear motions from the charge separation process, would the resultant photocell exhibit a better current and power output? Besides its theoretical relevance, this comparison is experimentally motivated as such decoupling may be feasible via optical cavities.32,33

As we will discuss in the next section, the mean current and power of our photocell is determined by the steady-state population of the secondary CT state |α〉. A non-perturbative computation of the steady-steady state under the influence of the full spectral density J(ω) with its 48 sharp modes per electronic state is quite challenging and out of the scope of our computational capabilities. We can however compute non-perturbative dynamics and steady state of our photocell including coherent interactions among single-excitation states and under the influence of JD(ω), J1(ω) and J2(ω) using the hierarchical equations of motion34–36 (see Methods). For comparison we also investigate the dynamics and steady state predicted by a simpler Pauli master equation for state populations, with transfer rates as described in ref. 26 (see Methods). Theoretical justification of the validity of this approximate framework is presented in ESI Note 3. Fig. S2 and S4 show that the population dynamics and steady-state populations predicted by the accurate framework and the Pauli master equation agree qualitative and quantitatively. Furthermore it has been shown that this approximate scheme accurately reproduces transient and steady-state spectroscopy of the PSIIRC.26 The main difference we see is that, as expected, the accurate framework predicts some short-lived excitonic coherences (see Fig. S3) which we show to have negligible influence on the current and power delivered by our photocell.

Photocell current and power performance

We fix the rate ΓL at which electrons are injected and set γex to simulate excitation by concentrated solar radiation,37,38 ensuring detailed balance as specified in ESI Note 2. The generated steady-state current passing between the system and drain electrode is equivalent to the current flowing across a hypothetical load connecting the final states |α〉 and |β〉, which have an associated energy gap Eαβ = EαEβ. The voltage V across such a load quantifies the extractable energy from our photocell with final energy gap Eαβ and an expression for V can be derived following standard thermodynamic considerations of photocells39 and photochemical systems.40 Denoting the steady state populations of the secondary CT state ραα and the ‘empty’ state ρββ, the load voltage V can be expressed as:
image file: c7sc02983g-t4.tif(3)
where kB is the Boltzmann constant, T the temperature of the photocell and e the electric charge. The average current 〈I〉 and the power output P delivered by the photocell are given, respectively, by 〈I〉 = Rραα and P = 〈IV. By fixing all parameters except the rate ΓR at which electrons leave the system, we can then investigate the characteristic 〈I〉–V and PV curves which define the photovoltaic performance of a photocell.

Fig. 3 presents the characteristic curves for the four spectral densities depicted in Fig. 2. For JD(ω), J1(ω) and J2(ω), we present the results obtained both by the hybrid framework and by the approximate Pauli master equation, showing their remarkable agreement for the current and power predictions (consistent with the dynamics shown in Fig. S2 and S4). In all the cases, the limit of ΓR → 0 leads to 〈I〉 → 0 defining the maximum available voltage or open circuit regime eVoc which is proportional to the energy gap between the ground state and the state that is directly photo-excited38i.e. eVocE1Eg ≈ 1.8 eV. In the opposite limit, when ΓR → ∞, V → 0. In these two extremes the photocell delivers no power. For all spectral densities we observe that the current is constant at low voltages and drops off at a characteristic voltage V comparable to Eαβ when the spectral density is JD(ω). This characteristic voltage increases slightly as the spectral density includes more well-resolved modes.

image file: c7sc02983g-f3.tif
Fig. 3 Photocell mean current and power versus voltage for different spectral densities. (a) Current. (b) Power. Solid lines are calculated using HEOM hybrid model while dotted lines are calculated using the Pauli master equation. Calculations carried out at 300 K with excitation rate γex = 75 cm−1. See ESI Note 1 for all other parameters.

There are two remarkable features to highlight in Fig. 3. First, the constant current observed for small voltages and the maximum power are significantly lower for a photocell with a structured environment. We observe such a reduction even in the case of J1(ω) whose central frequency is quasi-resonant with several exciton energy gaps. This contrasts with the power enhancement predicted for a simple light-harvesting unit operating under coherent interactions between all states.41 Second, the current delivered by a photocell with J2(ω) is already quite close to that of J(ω), confirming that for the performance of the photocell, J2(ω) is a good approximation to the full spectral density.

The behaviour predicted for the current and power delivered by the PSIIRC photocell relies on two main features characterizing the ChlD1 charge separation pathway: the large reorganisation energies of the CT states and the concomitant weak coupling between the primary CT state and single-excitation states. These two facts lead to an incoherent population dynamics that is well described by second order perturbation theory in the electronic coupling and such that the rates of transfer between an exciton |X〉 and the primary CT state |I〉 are dominated by the convolution of their line shape functions.42,43 Specifically, the rate of transfer between these states is given by kX,I = |VX,I|2SX,I where VX,I is the effective electronic coupling between the states and SX,I quantifies the spectral overlap:

image file: c7sc02983g-t5.tif(4)
with ωX,I the energy gap between the states, λX(I) the corresponding reorganisation energies for each state and gX(I)(t) the associated line broadening functions. Full expressions for these functions can be found in ESI Note 2. Fig. S5 shows the overlap between the low-lying exciton |X1〉 and |I〉 for the three spectral densities. As the spectral density contains more peaks, the reorganization energy of both states increases. However, for all the cases the values for λI are about one order of magnitude larger than λX (see Table S5). The larger λI, the wider is the shift between donor and acceptor states and as such the overlap SX,I accounting for spectral resonances is reduced to yield a lower rate of transfer kX,I. Similar considerations apply for the transfer rate between the CT states |I〉 and |α〉 whose reorganisation energies satisfy λI < λα as discussed in ESI Note 4. In this scenario the condition of quasi-resonance between electronic gaps and well-resolved vibrations becomes irrelevant for population transfer. The reduced transfer rates from excitons to the intermediate CT state |I〉 can then be interpreted as a Zeno-like effect whereby the strongly coupled environment “measures” the population of the CT state at a very high rate thereby slowing transfer. The disparity between reorganisation energies for CT states and the donor exciton states ensures a downhill relaxation that in the biological context counteracts charge recombination.

Zero-frequency noise

Full performance characterization of a photocell cannot be limited to the steady-state current and power output. Steady-state current fluctuations quantified by second and higher order cumulants (〈Inc with n > 1) can give information on the microscopic mechanisms underlying correlations between elementary charge transfer events.10 While the theory of full counting statistics is well established for Markovian systems,29,31 for non-Markovian dynamics only a few comprehensive frameworks based on perturbative approaches have been put forward.44 To address this shortcoming, we have developed a non-perturbative formalism to compute the full counting statistics,45 which integrates the exact system dynamics provided by the hierarchical equations of motions with a recursive scheme that allows accurate computation of the current cumulants44 (see Methods).

We focus on the long-time limit or zero-frequency regime of the relative noise strength which is quantified by the second order Fano factor:46

image file: c7sc02983g-t6.tif(5)

This ratio between the zero-frequency second order current cumulant and the mean current quantifies the deviation of the underlying statistical process from a Poissonian distribution. A Fano factor of 1 indicates a Poissonian process without correlations among charge transfer events. Deviations from 1 are interpreted as either super-Poissonian (F(1) > 1) or sub-Poissonian (F(2) < 1), regimes that can be associated with highly fluctuating or more stable currents, respectively.

Fig. 4 reports both non-perturbative and approximate results for F(2)versus V at room temperature for JD(ω) (Fig. 4(c)) and J2(ω) (Fig. 4(b)). For simplicity, we have omitted the results for J1(ω) as they simply follow the same trend. The qualitative and quantitative agreement between these curves indicate that zero-frequency noise properties of the photocell for these spectral densities is dominated by a population dynamics that is well-captured by the approximate Pauli framework. This is also consistent with the fact that exciton coherences arising from the interaction with a slowly relaxing bath or well-resolved vibrational motions decay much faster than the time scale on which the steady state is reached as discussed in the ESI. These arguments then justify the use of the approximate scheme to obtain insights into the behaviour of F(2) for the full spectral density J(ω) as shown in Fig. 4(b).

image file: c7sc02983g-f4.tif
Fig. 4 Fano factor versus voltage. (a) Fano factor for the PSIIRC photocell with the structured spectral density J(ω) for a modified (slow) (dotted line) and the measured (solid line) secondary charge transfer rate. (b) Fano factor for the PSIIRC photocell with spectral density J2(ω) containing two underdamped vibrations. (c) Fano factor for the PSIIRC photocell with a smooth low-energy vibrational environment JD(ω). Calculations carried out at 300 K with excitation rate γex = 75 cm−1. See ESI Note 1 for all other parameters.

In mesoscopic and quantum systems the zero-frequency Fano factor has proven to be very sensitive to the structure of the state space47 as well as to system-environment interactions.13 This is precisely what is indicated by the results for the Fano factor of the current through our photocell device shown in Fig. 4(a) and (b), which show that the noise profiles for J2(ω) and J(ω) each have a single minimum but are not symmetric with respect to the voltage at which this minimum occurs. In both cases we see that for V → 0 we have F(2) < 1 while in the opposite limit of VEα,β we obtain F(2) → 1. This indicates that at small voltages the electron transport is slightly correlated in all cases. As the spectral density exhibits more structure, the noise levels are overall lower. For instance, for the full spectral density we have F(2) (V = 0) = 0.90 while for J2(ω), we have F(2) (V = 0) = 0.95. Similarly, the minimum of F(2) reaches lower values as the spectral density acquires more structure i.e. F(2) = 0.55 for J(ω) which is less than the values obtained for J2(ω) i.e. F(2) = 0.6 and for JD(ω) i.e. F(2) = 0.7. In all the cases the sub-Poissonian behaviour is a manifestation of the Coulomb blockade regime where the presence of an electron in the system prevents another one entering until the system is empty. However, the more “ordered” transport observed for spectral densities with more well-resolved spectral features relies on the rapid population transport among excitons induced by such structured vibrational environments. The rates of transport to CT states, while reduced, are still comparable to the transfer among excitons such that the statistics of transitions from states |α〉 to |β〉 samples the manifold of exciton states donating population to the primary CT state. This hypothesis is confirmed by our analysis of the energy scale determining the voltage at which F(2) is minimum in each case of Fig. 4.

The analytic form of the Fano factor for our photocell model is too cumbersome to give any insight into the conditions determining the minima in Fig. 4. In order to rationalise the minimum in each curve it is useful to consider the case of a single resonant level (SRL) in the infinite bias limit.29 The dynamics of this system in the basis {|occupied〉, |empty〉} is governed by a Liouvillian with matrix elements [script L]11 = −[script L]21 = ΓR and [script L]22 = −[script L]12 = ΓL such that the Fano factor as a function of the voltage of a load across the occupied and empty states (eqn (3)) takes the form

image file: c7sc02983g-t7.tif(6)
where E0 is the energy gap between the occupied and empty states. It is simple to show that this expression is equivalent to writing the Fano factor in terms of ΓL and ΓR as image file: c7sc02983g-t8.tif (cf.eqn (45) in ref. 29 and see ESI Note 6). Fig. S2 shows that F(2)(V) for the SRL exhibits a single minimum, just as in our PSIIRC photocell. The minimum occurs when Vmin = E0 which is equivalent to the condition of ΓL/ΓR = ρoccupied/ρempty = 1 as shown in ESI Note 6. In this case, however, the function is symmetric about Vmin approaching 1 at both large and small voltages and indicating that electron transfer events in these extremes are uncorrelated. Based on this, we can say that near the Vmin the noise in our PSIIRC is approximately equal to that of an effective SRL with occupied level |α*〉, empty level |β〉 and renormalized energy gap image file: c7sc02983g-t9.tif that determines Vmin. Denoting as Ejk the energy gap between states |j〉 and |k〉 of our photocell, we notice that for the case of JD(ω) the Fano factor has a minimum for VminEI,β ≈ 1.50 eV while for the full spectral density VminEX6,β ≈ 1.56 eV and for the two mode spectral density J2(ω) we have, as expected, a value in between. This indicates that in the photocell with the full spectral density, the minimum noise samples the largest energy gap between |β〉 and the exciton manifold while for JD(ω) the noise only witnesses the energy gap up to the intermediate CT state. This is consistent with the fact that the rate of transfer among excitons in the D1 branch are larger for J(ω) than for the other two spectral densities.

As mentioned above, the most important feature of Fig. 4(a) and (b) is the non-symmetric profile F(2) with respect to Vmin. We now show that the rate limiting this asymmetric behaviour is the secondary charge transfer rate. To do this we consider the situation where the rate of secondary CT transfer kI,α is set by hand to a very low value compared to relaxation rates within the exciton manifold and the rates between the excitons and primary CT state for the full spectral density. In this case, the population of |α〉 is so slow that all the internal transfers from the exciton manifold to the primary CT can be described as a single step process i.e. there is no sampling of the exciton manifold and the Fano factor tends to 1 for small and larger voltages as shown by the dotted line in Fig. 4(a). This means that the system behaves as a SRL for all values of V with renormalized image file: c7sc02983g-t10.tif (cf.Fig. 4(a) with Fig. S6).

To conclude, our results at room temperature indicate that the PSIIRC-based photocell with the structured vibrational environment delivers less power than a photocell with an unstructured environment, yet this is accompanied by a suppression of current fluctuations. This noise reduction is supported by system–environment interactions that achieve rapid population transfer between excitons while still ensuring a downhill transfer to CT states.


Present single-molecule technologies demonstrating the ability to manipulate and measure the photocurrent of single PRCs have motivated us to investigate how the microscopic physical mechanisms underlying the function of these complexes may affect their performance as components in a photovoltaic cell. To approach this question we have brought together biological and physical perspectives by considering a PSIIRC-based photocell in a protein configuration that is argued to favour water-splitting20 and in which charge separation is initiated at the accessory ChlD1.17–19 This configuration is so far characterized by the lack of coherent delocalization between excitons and CT states.19

Our results show that a structured environment assisting electron transfer in our model PSIIRC-based photocell device acts to reduce the current and power output in comparison to a situation where electron–vibration interactions are described by a simple smooth low-energy background function. This reduction is a manifestation of a Zeno-like effect whereby the weak donor–acceptor electronic coupling concomitant with the stronger coupling of CT states to well-resolved high energy modes lead to slower transfer rates to CT states. These observations suggest that while PSIIRC complexes operating in the ChlD1 pathway may favour water oxidation under in vivo conditions, they may not necessarily be well suited for maximizing current output in single-molecule photovoltaics. Notwithstanding, the predicted reduction in the average photocurrent upon inclusion of coupling to well-resolved high-energy modes is not detrimental for the biological operation of PRCs. In the biological scenario it is more important to inhibit charge recombination and to ensure the captured energy is not wasted. In the ChlD1 transfer pathway, stronger coupling of CT states to these well-resolved vibrational motions ensures downhill relaxation thereby helping to prevent charge recombination.

For anthropogenic purposes of obtaining the largest current out of these units regardless of its fluctuations, the best strategy then may be to decouple specific vibrational motions from electronic states. Indeed, modifications of the electron–nuclei interactions can be achieved by strong coupling of pigment–protein complex to a confined optical cavity mode.32,33 In particular, ref. 33 has shown that the energy exchange of electronic transitions with a strongly coupled optical mode could help suppressing reorganisation energy of the nuclei thereby increasing the rate of electron transfer reactions. Alternatively, one can select PSIIRCs operating in the PD1PD2 pathway where coherent delocalization across exciton and CT states has been probed19,28,48,49 and which will lead to an enhancement as predicted in ref. 41.

While no advantage is obtained in terms of mean current and power output, strong coupling to well resolved vibrational modes results in a reduction of current fluctuations of our PSIIRC-based photocell. This lower noise strength obeying a sub-Poissonian statistics and the associated ordered electron transport is promoted by the exciton manifold and signals out the multi-step nature of the transport process. Preliminary calculations (not shown) for a photocell with delocalized states across excitons and CT states indicate that such a noise reduction maybe a general feature.

From the electronic-device viewpoint, reducing any kind of noise is always a desirable feature to guarantee device resolution; this includes intrinsic noise due to the inherently probabilistic nature of the process. Hence the device functionality of this noise reduction appears to be straightforward: to improve precision in the current delivered. More interesting is to discuss the possible advantages of such noise reduction in the biological context. It is well known that noise and its control is crucial across all scales in biology.50,51 For instance, it has been discussed that biochemical processes that are inherently stochastic include mechanisms to control intrinsic noise and, in particular, to reduce it for regulatory processes.51 Indeed, electron transfer events in photosynthetic reaction centres belong to a larger family of stochastic transport processes in biology, some of which have already been predicted to exhibit mechanisms suppressing fluctuations below the Poisson level.52 Moreover, increased complexity in biological networks has been linked to intrinsic noise reduction.53 We therefore argue that, for biological function of PRCs, the coupling to well-resolved vibrations and the predicted noise reduction could indeed have a regulatory function. In these systems, the final stable CT state |α〉 donates an electron to quinone B and once this reduction happens, the PRC is unable to handle an excitation during a finite time. Having single electrons delivered at regular (ordered) time intervals (with narrower fluctuations of waiting times) as opposed to randomly (Poisson-like process) could avoid wasting excitations during such overly long blocking periods.

The experimental implementation of our proposal assumes PSIIRC units that have been modified to have no quinones as has been done for the protein samples used in ref. 28. This will ensure that attachments of electrodes to individual PSIIRC units are at the level of the electron donor and electron acceptor pigments. We envision metal–protein junctions and a scanning tip microscopy setup as those that have been realised for photosystem I (PSI) units9 and for reaction centre-enriched purple bacterial membranes.4 Genetic manipulation of both the oxidizing and reducing sites could allow covalent attachment of the protein to the electrodes across which the photo-current could be measured.9 To measure the elementary transfer events from the electron acceptor site to the drain electrode we envision a device, such as a silicon field-effect transistor,54,55 capable of detecting single charges and feasible to be integrated in the scanning tip setup at room temperature as shown in Fig. 1(a).

The main limitation of isolated natural photosynthetic proteins for realistic, long-lived photovoltaic applications lies in the photodamage they experience. In PSIIRC this occurs in the D1 protein resulting in a lifetime as short as tens of minutes.56,57 Emerging organic alternatives, such as the synthesis of man-made protein maquettes58 has opened the possibility of building nanometric units accurately mimicking the structure of photosynthetic systems yet displaying enhanced photostability. Merging this area with photovoltaics may unleash an unforeseen remarkable development.

From the theoretical view point a few remarks must be made. As specified in the Methods section, we assume a simplified phenomenological description of the coupling between the electronic system and the leads. This approximation has been used in similar systems with arbitrary electron–phonon couplings59,60 and has been shown to give relevant physical insight. However, this treatment neglects the possibility of the leads coupling to localized vibronic states rather than to bare electronic degrees of freedom. The full extent of this effect in our system remains to be investigated. The inclusion of incoherent transfer of population to CT states via a Lindblad dissipator with rates given by Generalized Förster maybe be questioned. While this form guarantees physical consistency it negates the possibility of any coherent transfer even if very weak. A full treatment of the problem would include coherent coupling to the primary and secondary CT states within the HEOM. There are also additional questions about the thermodynamic consistency of calculating the power output across with the phenomenologically modeled load as raised in ref. 61. It will therefore be important to investigate entropy production62 for our model photocell to assess its consistency with the second law.

Finally, our work has focused on the zero-frequency noise showing that, in this case, it is dominated by the population dynamics as confirmed by our comparison between the hybrid and the approximate frameworks. An extension of our study to investigate finite-frequency noise63 could therefore be a suitable alternative to obtain signatures of quantum coherence. More generally, current statistics measurements also potentially offer a non-invasive, single system level probe of charge transfer phenomena in a wide range of biological64,65 and chemical systems66,67 ranging from charge transfer along molecular wires made from DNA strands64 to general donor–bridge–acceptor systems66,67 or to unveil vibrational mechanisms for odour receptors.65


Dynamical evolution of electronic excitations

We consider an exciton dynamics described by image file: c7sc02983g-t11.tif where i corresponds to the basis of single-excitation states of the six core chromophores i.e. {|PD1〉, |PD2〉, |ChlD1〉, |ChlD2〉, |PheD1〉, |PheD2〉} Tij is the electronic coupling between states and ei are onsite energies given in Table S1. The six eigenstates of Hel are denoted as |X1〉 to |X6〉 with corresponding eigenenergies EX1 to EX6 in ascending order. The electronic operators |i〉〈i| couple linearly with coupling gi to identical baths of harmonic oscillators image file: c7sc02983g-t12.tif. The strength of the system-bath interaction is quantified by the spectral density that will be of the form J(ω), JD(ω), J1(ω) or J2(ω).

To describe the full PSIIRC photocell dynamics under the operation conditions illustrated in Fig. 1 we consider two frameworks. (1) A hybrid framework that accounts for a non-perturbative approach to the exciton dynamics using the hierarchical equations of motion34–36 in combination with incoherent transfer rates68 to and from all other states. The non-perturbative expansion of the exciton dynamics is used to account accurately for the effects of JD(ω), J1(ω) or J2(ω). Incoherent rates connecting the exciton states with the rest of states in the photocell are defined using a Lindblad dissipator coupled to each auxiliary density matrix in the expansion.68 ESI Note 2 presents further details of the hierarchical expansion of exciton dynamics under this scheme. Converged dynamics are obtained by terminating the hierarchical expansion at level N = 8 for JD(ω) and level N = 5 for J1(ω) and J2(ω). Only the K = 0 Matsubara term was explicitly accounted for, though a Markovian truncation term for Matsubara frequencies was included to capture some finite temperature effects.35

(2) A Pauli master equation for electronic state populations is also considered, similarly to the approach followed on ref. 26. The Pauli rate equations have the form |〉〉 = M|P〉〉, where |P〉〉 is a vector of state populations in the basis {|g〉, |X1〉, …, |X6〉, |I〉, |α〉, |β〉} and M is a stochastic matrix containing the rates for transfer between these electronic states. Modified Redfield theory as presented in ESI Note 2 is used to compute population transfer among exciton states.43,69 In both frameworks (1) and (2) we assume weak and incoherent coupling from excited states to the primary CT state as well as weak coupling between charge transfer states. The transfer from exciton states |Xn〉 to the intermediate CT state |I〉 are given by Generalised Förster theory, and Förster-like rates are used to describe transfer between the CT states |I〉 and |α〉.70 Other incoherent rates are described in ESI Note 2. Theoretical validity of this framework is discussed in ESI Note 3 along with a systematic comparison of the predictions of frameworks (1) and (2).

Theory of full counting statistics

We envisage our photocell positioned between source and drain leads which supply or remove electrons from the system respectively. The leads are taken as weakly coupled fermionic reservoirs and their influence is described by phenomenological Lindblad-type dissipators,30 as specified in ESI Note 5. With weak coupling to the leads, the theory of full counting statistics29,31,63 provides a framework to investigate the cumulants of the current passing through the system. This framework is applied to both the hybrid, non-perturbative approach, and the approximate Pauli model we use to describe the dynamics. For both models a time-local master equation [small sigma, Greek, dot above](t) = [scr M, script letter M]σ(t) can be constructed, where the state of the system σ(t) is propagated through time by an operator [scr M, script letter M]. This dynamical equation is augmented by a counting field χ, used to single out the incoherent transition across which the electron statistics is counted. For our photocell this is the transition from state |α〉 to state |β〉 where an electron is transferred to the drain lead. This leads to the time propagator [scr M, script letter M](χ) = [scr M, script letter M]0 + eiχ[scr M, script letter M]J where [scr M, script letter M]0 describes the time evolution of the system between counting events and [scr M, script letter M]J is the jump matrix describing hopping events between the system and the drain lead. Further details on the calculation of non-perturbative electron counting statistics using the hierarchical equations of motion45 are given in ESI Note 5. The zero-frequency cumulants are encoded in the probability distribution of the number of electrons that hop into the drain lead in some long time period.29 A recursive scheme is then followed which generates zero-frequency current cumulants up to any order44 and expresses them in terms of the jump matrix [scr M, script letter M]J and the pseudo-inverse [scr R, script letter R] of the time propagator. The mean 〈I1〉 and noise 〈I2〉 are given by
I1〉 = ≪[0 with combining tilde]|[scr M, script letter M]J|0≫(7)
I2〉 = ≪[0 with combining tilde]|[scr M, script letter M]J − 2[scr M, script letter M]J[scr R, script letter R][scr M, script letter M]J|0≫(8)
where ≪[0 with combining tilde]| and |0≫ are the left and right steady state eigenvectors of the time propagator.

Author contributions statement

A. O.-C. designed the research, R. S. and H. H.-N. carried out the calculations. RS, H. H.-N., RvG and A. O.-C. analysed the results and wrote the manuscript.

Conflicts of interest

The authors declare no competing financial interests.


The authors would like to thank Clive Emary, Elisabet Romero, Vladimir Novoredezkin and Jeroem Elzerman for helpful discussions. Financial support from the Engineering and Physical Sciences Research Council (EPSRC UK) Grant EP/G005222/1 and from the EU FP7 Project PAPETS (GA 323901) is gratefully acknowledged.


  1. R. Blankenship, Molecular Mechanisms of Photosynthesis, Wiley-Blackwell, 2001 Search PubMed.
  2. R. Blankenship, D. Tiede, J. Barber, G. Brudvig, G. Fleming, M. Ghirardi, M. Gunner, W. Junge, D. Kramer, A. Melis, T. Moore, C. Moser, D. Nocera, A. Nozik, D. Ort, W. Parson, R. Prince and R. Sayre, Science, 2011, 332, 805–809 CrossRef CAS PubMed.
  3. E. Wientjes, H. van Amerongen and R. Croce, Nat. Commun., 2013, 117, 11200–11208 CAS.
  4. M. Kamran, V. Friebe, J. Delgado, T. Aartsma, R. Frese and M. Jones, Nat. Commun., 2015, 6, 653 Search PubMed.
  5. O. Yehezkeli, R. Tel-Vered, D. Michaeli, I. Willner and R. Nechushtai, Photosynth. Res., 2014, 120, 71–85 CrossRef CAS PubMed.
  6. M. Grätzal, Nature, 2001, 414, 338–344 CrossRef PubMed.
  7. B. Reiss, D. Hanson and M. Firestone, Biotechnol. Prog., 2007, 23, 985–989 CrossRef CAS PubMed.
  8. T. Mikayama, T. Miyashita, K. Iida, Y. Suemori and M. Nango, Mol. Cryst. Liq. Cryst., 2006, 445, 291–296 CAS.
  9. D. Gerster, J. Reichert, H. Bi, J. Barth, S. Kaniber, A. Holleitner, I. Visoly-Fisher, S. Sergani and I. Carmeli, Nat. Nanotechnol., 2012, 7, 673–676 CrossRef CAS PubMed.
  10. Y. Blanter and M. Büttiker, Phys. Rep., 2000, 336, 1–166 CrossRef CAS.
  11. G. Kießlich, E. Schöll, T. Brandes, F. Hohls and R. Haug, Phys. Rev. Lett., 2007, 99, 206602 CrossRef PubMed.
  12. N. Ubbelohde, C. Fricke, C. Flindt, F. Hohls and R. Haug, Nat. Commun., 2012, 3, 612 CrossRef PubMed.
  13. F. Haupt, F. Cavaliere, R. Fazio and M. Sassetti, Phys. Rev. B: Condens. Matter Mater. Phys., 2006, 74, 205328 CrossRef.
  14. W. Belzig, Phys. Rev. B: Condens. Matter Mater. Phys., 2005, 71, 1–4 CrossRef.
  15. H. Hossein-Nejad, A. Olaya-Castro, F. Fassioli and G. Scholes, New J. Phys., 2013, 15, 083056 CrossRef.
  16. A. W. Rutherford and P. Faller, Philos. Trans. R. Soc. London, Ser. B, 2003, 358, 245–253 CrossRef CAS PubMed.
  17. M. Groot, N. Pawlowicz, L. van Wilderen, J. Breton, I. van Stokkum and R. van Grondelle, Proc. Natl. Acad. Sci. U. S. A., 2005, 102, 13087–13092 CrossRef CAS PubMed.
  18. A. Holzwarth, M. Müller, M. Reus, M. Nowaczyk, J. Sander and M. Rögner, Proc. Natl. Acad. Sci. U. S. A., 2006, 103, 6895–6900 CrossRef CAS PubMed.
  19. E. Romero, I. van Stokkum, V. Novoderezhkin and J. Dekker, Biochem., 2010, 49, 4300–4307 CrossRef CAS PubMed.
  20. T. Renger and E. Schlodder, J. Photochem. Photobiol., B, 2011, 104, 126–141 CrossRef CAS PubMed.
  21. V. Novoderezhkin, E. Andrizhiyevskaya, J. Dekker and R. van Grondelle, Biophys. J., 2005, 89, 1464–1481 CrossRef CAS PubMed.
  22. V. Novoderezhkin, A. Yakovlev, R. van Grondelle and V. Shuvalov, J. Phys. Chem. B, 2004, 108, 7445–7457 CrossRef CAS.
  23. Y. Umena, K. Kawakami, J. Shen and N. Kamiya, Nature, 2011, 473, 55–60 CrossRef CAS PubMed.
  24. B. Diner and F. Rappaport, Annu. Rev. Plant Biol., 2002, 53, 551–580 CrossRef CAS PubMed.
  25. M. Steffen, K. Lao and S. Boxer, Science, 1994, 264, 810–816 CAS.
  26. V. Novoderezhkin, E. Romero, J. Dekker and R. van Grondelle, ChemPhysChem, 2011, 12, 681–688 CrossRef CAS PubMed.
  27. P. Peterman, H. van Amerongen, R. van Grondelle and J. Dekker, Proc. Natl. Acad. Sci. U. S. A., 1998, 95, 6128–6133 CrossRef.
  28. E. Romero, R. Augulis, V. Novoderezhkin, M. Ferretti, J. Thieme, D. Zigmantas and R. van Grondelle, Nat. Phys., 2014, 10, 676–682 CrossRef CAS PubMed.
  29. D. Marcos, C. Emary, T. Brandes and R. Aguado, New J. Phys., 2010, 12, 123009 CrossRef.
  30. U. Harbola, M. Esposito and S. Mukamel, Phys. Rev. B: Condens. Matter Mater. Phys., 2006, 74, 235309 CrossRef.
  31. D. Bagrets and Y. Nazarov, Phys. Rev. B: Condens. Matter Mater. Phys., 2003, 67, 085316 CrossRef.
  32. D. Coles, Y. Yang, Y. Wang, R. Grant, R. Taylor, S. Saikin, A. Aspuru-Guzik, D. Lidzey, J. Kuo-Hsiang Tang and J. Smith, Nat. Commun., 2014, 5, 5561 CrossRef CAS PubMed.
  33. F. Herrera and F. Spano, Phys. Rev. Lett., 2016, 116, 238301 CrossRef PubMed.
  34. A. Ishizaki and G. Fleming, J. Chem. Phys., 2009, 130, 234111 CrossRef PubMed.
  35. Q. Shi, L. Chen, G. Nan, R. Xu and Y. Yan, J. Chem. Phys., 2009, 130, 084105 CrossRef PubMed.
  36. Y. Tanimura, J. Chem. Phys., 2012, 137, 22A550 CrossRef PubMed.
  37. K. Dorfman, D. Voronine, S. Mukamel and M. Scully, Proc. Natl. Acad. Sci. U. S. A., 2013, 110, 2746–2751 CrossRef CAS PubMed.
  38. C. Creatore, M. Parker, S. Emmott and A. Chin, Phys. Rev. Lett., 2013, 111, 253601 CrossRef CAS PubMed.
  39. W. Shockley and H. Queisser, J. Appl. Phys., 1961, 32, 510 CrossRef CAS.
  40. R. Ross and M. Calvin, Biophys. J., 1967, 7, 595–614 CrossRef CAS PubMed.
  41. N. Killoran, S. Huelga and M. Plenio, J. Chem. Phys., 2015, 143, 155102 CrossRef CAS PubMed.
  42. T. Förster, Modern Quantum Chemistry, Part III. B. Action of Light and Organic Crystals, Academic Press, New York, 1965, pp. 93–137 Search PubMed.
  43. M. Yang and G. Fleming, Chem. Phys., 2002, 275, 355–372 CrossRef CAS.
  44. C. Flindt, T. Novotný, A. Braggio and A.-P. Jauho, Phys. Rev. B: Condens. Matter Mater. Phys., 2010, 82, 155407 CrossRef.
  45. R. Stones and A. Olaya-Castro, arXiv:1705.02320, 2017.
  46. U. Fano, Phys. Rev., 1947, 72, 26–29 CrossRef CAS.
  47. J. Egues, S. Hershfield and J. Wilkins, Phys. Rev. B: Condens. Matter Mater. Phys., 1994, 49, 13517–13527 CrossRef CAS.
  48. F. Fuller, J. Pan, A. Gelzinis, V. Butkus, S. Seckin Senlik, D. Wilcox, C. Yocum, D. Valkunas, L. Abramavicius and J. Ogilvie, Nat. Chem., 2014, 6, 706–711 CAS.
  49. E. Romero, V. Novoderezhkin and R. van Grondelle, Nature, 2017, 543, 355–365 CrossRef CAS PubMed.
  50. L. Tsimring, Rep. Prog. Phys., 2014, 77, 026601 CrossRef PubMed.
  51. J. Raser and E. O'Shea, Science, 2005, 309, 2010–2013 CrossRef CAS PubMed.
  52. R. Brunetti, F. Affinito, C. Jacoboni, E. Piccinini and M. Rudan, J. Comput. Electron., 2007, 6, 391–394 CrossRef CAS.
  53. L. Cardelli, A. Csikász-Nagy, N. Dalchau, M. Tribastone and M. Tschaikowski, Sci. Rep., 2016, 6, 20214 CrossRef CAS PubMed.
  54. K. Nishiguchi and A. Fujiwara, Nanotechnology, 2009, 20, 175201 CrossRef CAS PubMed.
  55. K. Nishiguchi, Y. Ono and A. Fujiwara, Appl. Phys. Lett., 2011, 98, 193502 CrossRef.
  56. K. Brinkert, F. Le Formal, X. Li, J. Durrant, A. Rutherford and A. Fantuzzia, Biochim. Biophys. Acta, 2016, 1857, 1497–1505 CrossRef CAS PubMed.
  57. M. Edelman and A. Mattoo, Photosynth. Res., 2008, 98, 609–620 CrossRef CAS PubMed.
  58. B. Lichtenstein, C. Bialas, J. Cerda, B. Fry, P. Dutton and C. Moser, Angew. Chem., Int. Ed., 2015, 54, 13626–13629 CrossRef CAS PubMed.
  59. A. Braggio, C. Flindt and T. Novotný, J. Stat. Mech., 2009, 2009, P01048 CrossRef.
  60. D. Santamore, N. Lambert and F. Nori, Phosphorus Res. Bull., 2013, 87, 075422 Search PubMed.
  61. D. Gelbwaser-Klimovsky and A. Aspuru-Guzik, Chem. Sci., 2017, 8, 1008 RSC.
  62. M. Esposito, K. Lindenberg and C. Van den Broeck, New J. Phys., 2010, 12, 013013 CrossRef.
  63. C. Emary, D. Marcos, R. Aguado and T. Brandes, Phys. Rev. B: Condens. Matter Mater. Phys., 2007, 76, 161404 CrossRef.
  64. L. Xiang, J. Palma, C. Bruot, V. Mujica, M. Ratner and N. Tao, Nat. Chem., 2015, 7, 221–226 CrossRef CAS PubMed.
  65. M. Franco, L. Turin, A. Mershin and E. Skoulakis, Proc. Natl. Acad. Sci. U. S. A., 2011, 108, 3797–3802 CrossRef CAS PubMed.
  66. M. Delor, I. Sazanovich, M. Torie and J. Weinstein, Acc. Chem. Res., 2015, 48, 1131–1139 CrossRef CAS PubMed.
  67. A. Bakulin, R. Lovrincic, X. Yu, O. Selig, H. Bakker, Y. Rezus, P. Nayak, A. Fonari, V. Coropceanu, J. Bredas and D. Cahen, Nat. Commun., 2015, 6, 7880 CrossRef CAS PubMed.
  68. C. Kreisbeck, T. Kramer, M. Rodríguez and B. Hein, J. Chem. Theory Comput., 2011, 7, 2166–2174 CrossRef CAS PubMed.
  69. W. M. Zhang, T. Meier, V. Chernyak and S. Mukamel, J. Chem. Phys., 1998, 108, 7763–7774 CrossRef CAS.
  70. E. O'Reilly, PhD thesis, University College London, 2014.


Electronic supplementary information (ESI) available. See DOI: 10.1039/c7sc02983g

This journal is © The Royal Society of Chemistry 2017