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: a.olaya@ucl.ac.uk
bDepartment of Physics and Astronomy, VU University, 1081 HV Amsterdam, The Netherlands
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.
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.
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
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
(1) |
(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. 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.
(3) |
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. 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 kX,I = |VX,I|2SX,I where VX,I is the effective electronic coupling between the states and SX,I quantifies the spectral overlap:
(4) |
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
(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).
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 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 form
(6) |
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 (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.
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
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).
〈I1〉 = ≪|J|0≫ | (7) |
〈I2〉 = ≪|J − 2JJ|0≫ | (8) |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c7sc02983g |
This journal is © The Royal Society of Chemistry 2017 |