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.

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 cell[1] 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 cells[5,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 setups[10-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 interactions[13] and coherence.[14] A theoretical study along the lines of counting statistics for light-harvesting complexes has been carried out[15] 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 complex[20] 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.

Results

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).

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 experiments[27] 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]

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 flow[29,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]

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. G(ω) describes the spectral density of an underdamped mode coupled to an excited pigment, with λ, ω and γ being the reorganisation energy, frequency and damping rate of mode j respectively. has been measured experimentally[27] 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 pathway[28] 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 motion[34-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 photocells[39] 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: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 = eΓ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 P–V 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-excited[38] i.e. eVoc ≈ E1 – Eg ≈ 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.

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 k = |V|2S where V is the effective electronic coupling between the states and S quantifies the spectral overlap:with ω the energy gap between the states, λ the corresponding reorganisation energies for each state and g(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 λ are about one order of magnitude larger than λ (see Table S5†). The larger λ, the wider is the shift between donor and acceptor states and as such the overlap S accounting for spectral resonances is reduced to yield a lower rate of transfer k. Similar considerations apply for the transfer rate between the CT states |I and |α whose reorganisation energies satisfy λ < λ 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 (Ic 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 cumulants[44] (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]

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).

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 space[47] 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 V ≫ E 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 ℒ11 = –ℒ21 = ΓR and ℒ22 = –ℒ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 formwhere 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 (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 that determines Vmin. Denoting as E 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 Vmin ≈ E ≈ 1.50 eV while for the full spectral density Vmin ≈ E ≈ 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 k 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 (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.

Discussion

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-splitting[20] 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 probed[19,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) units[9] 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 maquettes[58] 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 couplings[59,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 production[62] 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 noise[63] 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 biological[64,65] and chemical systems[66,67] ranging from charge transfer along molecular wires made from DNA strands[64] to general donor–bridge–acceptor systems[66,67] or to unveil vibrational mechanisms for odour receptors.[65]

Methods

Dynamical evolution of electronic excitations

We consider an exciton dynamics described by where i corresponds to the basis of single-excitation states of the six core chromophores i.e. {|PD1, |PD2, |ChlD1, |ChlD2, |PheD1, |PheD2} T is the electronic coupling between states and e are onsite energies given in Table S1.† The six eigenstates of Hel are denoted as |X1 to |X6 with corresponding eigenenergies E to E in ascending order. The electronic operators |ii| couple linearly with coupling g to identical baths of harmonic oscillators . 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 motion[34-36] in combination with incoherent transfer rates[68] 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 |X 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 statistics[29,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 σ̇(t) = ℳσ(t) can be constructed, where the state of the system σ(t) is propagated through time by an operator ℳ. 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 ℳ(χ) = ℳ0 + eiℳJ where ℳ0 describes the time evolution of the system between counting events and ℳ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 motion[45] 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 order[44] and expresses them in terms of the jump matrix ℳJ and the pseudo-inverse ℛ of the time propagator. The mean I1 and noise I2 are given by where ≪0| 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.

Click here for additional data file.