Coherent excitation transfer via the dark-state channel in a bionic system

Light absorption and energy transfer were studied in a bionic system with donors and an acceptor. In the optimal case of uniform couplings, this seemingly complicated system was reduced to a three-level Λ-type system. With this observation, we showed that the efficiency of energy transfer through a dark-state channel, which is free of the spontaneous decay of the donors, was dramatically improved. After the overall average transfer time was evaluated, it was revealed that the time required for transfer through this channel was not optimal. To find the optimal parameters for the present system, we defined a new quantity, output power, which characterizes the average mean output of the whole transfer process. We estimated the optimal parameters of the system to achieve the maximum output power. The splitting behaviour of the maximum power may be used to explain the phenomenon of photosynthetic systems primarily absorbing two colours of light. A molecular ring structure used for photosynthesis in plants could influence the design of future solar panel technologies. H. Dong and co-workers at the Chinese Academy of Sciences in Beijing have demonstrated a theoretical light-harvesting model that mimics photosynthesis in plants. A central acceptor surrounded by a ring of donor molecules transforms light into energy through the excitation of electrons. The researchers modelled donors coupled to photons of a particular energy and then examined the energy transfer process between the acceptor and donors of different molecular states. According to their results, the presence of a ‘dark state’ — which occurs when a molecule can no longer absorb or emit photons — significantly improves the overall efficiency of energy transfer in these artificial photosystems.


INTRODUCTION
The natural process of light harvesting in plants, algae and bacteria may provide the key to tackling the global problem of the lack of longterm, viable energy sources. 1,2 An understanding of the high-efficiency energy-transfer mechanism of the light-harvesting process would aid in the design of a new generation of clean solar energy sources. Recently, the long-time coherent properties of excitations in lightharvesting systems have been observed in experiments. [3][4][5][6] This coherence can be preserved in these structures for long periods of time, even at room temperature. It seems that the high efficiency of the light conversion is related to these coherent properties, which are quantum in nature. Therefore, the physical mechanism of photosynthesis assisted by quantum coherence has generated much experimental [3][4][5][6] and theoretical [7][8][9][10][11][12][13][14][15][16][17][18][19][20][21] research.
An inherent mechanism for high-efficiency energy transfer is, to our best knowledge, still unknown. One possible reason is the optimisation of the spatial structures in light-harvesting systems. X-ray analysis 22,23 has revealed some common elements shared by different light-harvesting complexes in nature, one of which is a ring structure with a reaction centre. 24,25 It is well known that natural selection preserves the features most advantageous for a biological system in its current environment. Therefore, it is believed that this ring structure is beneficial during the light-harvesting process. The mechanism behind this optimal structure may account for the high efficiency of the natural light-harvesting process. Thus, an investigation of the mechanisms present in similar systems would aid in the design of artificial light-harvesting systems with self-assembling molecular arrays and in the development of quantum dot arrays for future top-down semiconductor fabrication techniques.
In this paper, we will study a generic model, which is similar to the type I light-harvesting complex (LH I), a centralized acceptor surrounded by coupled donors arranged in a ring. In the natural lightharvesting process, the light-capture process and the excitation transfer occur simultaneously. To mimic the natural process, we include the light-capture process in the model by coupling the donors with photons in a single mode. By showing that the present model with homogeneous coupling can be reduced to the well-studied L system, we find that the overall transfer efficiency is insensitive to the decay of the donors when the eigenfrequency of the acceptor is resonant with the light. This discovery implies that the dark-state effect suppresses the excitation population of the noisy donors so that the energy transfer efficiency is dramatically improved. Otherwise, the loss of the donor excitations would largely decrease the efficiency of transfer to the acceptor. Figure 1a, our model concerning light capture and coherent excitation transfer is similar to the structure of the LH I. The one-dimensional circular array consisting of N donors is also analogous to the ring structures, such as the B800 or the B850 rings, in the light-harvesting complex II (LH II). The acceptor is placed at the centre of the ring. Similar to the pigment molecules in the natural LH I ring, each donor can be modelled as a two-level system (TLS) e j T and g j T with energy level spacing e i (i51, …, N). In practice, a TLS is a proper approximation for a single-excitation case only. In most natural conditions, the LHs always have only one excitation, 14 thus our discussion concerns only the zero-and one-excitation subspaces. Because the rate of hopping between the non-adjacent sites would be small, we consider only the adjacent hopping of excitations of strength g on the ring. Visible light constitutes the majority of the energy in the solar spectrum and has a wavelength of approximately 500 nm. In the cells of photosynthetic bacteria, the B800 and the B850 rings usually have radii of approximately 4-6 nm. 24 Thus, BChl molecules, the units of the B800 and the B850 rings, are coupled uniformly to incident light of frequency v, which is described by the creation (annihilation) operator b 1 (b).

EFFECTIVE HAMILTONIAN As illustrated in
To simulate the major function of this natural process, we adopt the uniform-site-energy case in the present model, namely, e i 5e. Then, the model Hamiltonian reads where H D~X j e i A z z ð h:c:Þ are the corresponding Hamiltonian of a donor and its interaction with the acceptor, respectively. And h.c. denotes the hermitian conjugate term here and thereafter. Here, e i z~e j T i Sgj and A zẽ j T A Sgj represent the creation of an excitation on the jth donor and on the acceptor, respectively. The physical implementation of the present structure could be a quantum dot array, where each dot has a size of approximately 2-10 nm and the interparticle spacing is approximately 10 nm. 26 The recent experimental synthesis of a 12porphyrin ring 27 also creates the possibility of designing a structure resembling a natural light-harvesting element.
In this paper, we primarily consider the effective transfer process and explore the advantage of the spatial configuration of the ring type. The collective excitation of the donor ring is described by the Fourier transformation whereẽ e kẽ e z k À Á is the annihilation (creation) operator of the collective mode with a definite momentum k. Here, the summation is over all the discrete momenta k n 52p(n21)/N, with n51, …, N. Indeed, in the large-N limit, we can show that ½ẽ e k ,ẽ e z k' ?d kk' , thus demonstrating that the collective excitations behave as bosons. 28 In terms of the boson-like operatorsẽ e kẽ e z k À Á , the donor Hamiltonian is seemingly diagonalized as which represents an energy band with N sub-energy levels. We note that the interaction term is rewritten as H DA~ffi ffiffiffi N p jẽ e 0 A z zh:c: ð Þ , which shows that only the zero-mode excitation with k50 is coupled to the acceptor, while the others described byẽ e kẽ e z k À Á (k=0) are decoupled with the acceptor. Therefore, the above Fourier transformation separates the total Hamiltonian into two decoupled parts: where v 0 5e12g2ik and v A~eA {iC. We must note that the dissipation rates k and C have been phenomenologically introduced to describe the loss of excitations for the donors and the acceptor, respectively. We also emphasize that only the single-excitation subspace is considered in the physical situation discussed above. Before the following discussion, we offer three remarks on the implications of the above effective Hamiltonian: 1. The coupling is strengthened by a factor ffiffiffiffi N p , which improves both the speed of light capture by the donors and the rate of excitation transfer to the accepter from the donor ring. 2. The zero-mode excitation described byẽ e z 0 has an energy E k~0ẽ z2g. Acting on the ground state 0 of the single localized excitations 1 j T~g 1 ,:::,g j{1 e j g jz1 :::g N T in the jth donor (j51, 2, …, N). It is similar to the single-magnon state in a spin-wave system. It has been proven numerically that the collective excited initial state 1 k~0 j T would result in the maximum efficiency with respect to any other mode. 8 3. In the single-excitation subspace, the photon-donor-acceptor system could be described as a three-level L-type system, as illustrated in Figure 1b. The boxed area is the corresponding singleexcitation subspace. A single photon can excite a donor, while the excitation of a donor can be transferred to the acceptor. Interestingly, the incident light only couples to the zero mode of the donor ring. Because all other modes are decoupled from the capture process and do not contribute to the light-harvesting process, the excitation energy is transferred to the acceptor only through the zero-mode channel. Therefore, the transfer efficiency would improve dramatically, if natural light-harvesting systems were optimized to contain a zero mode.

NUMERICAL CALCULATION
In the single-excitation case, the evolution of the excitation is constrained in the subspace spanned by 1 b j T: 1,0,0 j T, 1 D j T: 0,1,0 j T and 1 A j T: 0,0,1 j T, which represent the direct product state a,b,c j T: a j T6 b j T6 c j T of the photon, the excitation of the donor ring in the zero mode and the acceptor, respectively. Let w(t) j T be a single-excitation wave function with amplitudes u(t), v(t) and w(t) corresponding to the above basis vectors. The Schrödinger equation is Figure 1 (a) The donor-acceptor system with dissipation from both the donors (rate k) and the acceptor (rate C). In the natural light-harvesting process, the rate C represents the speed at which the excitation of the acceptor is utilized by outside agents in a chemical reaction; (b) the reduced energy spectrum in the zero-and the one-excitation subspaces with its structural optimisation.
Coherent excitation transfer via dark-state channel H Dong et al 2 The energy transfer is usually understood as the transportation of excitations from the donors to the acceptor. 8,17,18 The overall transfer efficiency 8 is given by an integral In practice, the rate of excitations hopping from a donor to the acceptor and the rate of excitation transfer to an outside agent should both be larger than the rate of excitation dissipation for excited donors, i.e., ffiffiffiffi N p j?k and C . k. As is typically done, we will set the parameters 24 to be e A /j512, J/j50.1, g/j50.3, e/j511.4 and C/j50.3 with j&10 ps {1 . In most cases, we set the total donor number to be N58. We illustrate the dependence of the efficiency on the detuning D5v2e A and the dissipation parameter k/C in Figure 2a.
There exists a large peak at D50 for a given k/C, as explicitly shown in Figure 2b. For a system comprising only a donor and an acceptor, where the donor is initially excited, there is an upper bound on the overall transfer efficiency g ' max~C = Czk ð Þ. For instance, this upper bound could have a value of g ' max k=C~0:5 ð Þ2=3. In the present discussion, where light capture is included, the efficiency actually surpasses the upper bound. We also note that the efficiency increases as N becomes larger and that there is no optimal N in the present model due to our neglect drawback of the donor number, which is discussed elsewhere.
It is also significant that the peak value of k/C remains nearly unchanged for the situation of resonance D50, as illustrated in Figure 2c. This observation reflects the fact that the efficiency is insensitive to the decay of the donor ring. To understand this phenomenon, we present the amplitudes of the photon mode, the donor zero mode and the acceptor in Figure 3a for the case of resonance D50. At resonance, the excitation amplitude of the donor is highly suppressed. The suppression of the amplitude eventually reduces the energy dissipation and increases the efficiency beyond the bound g ' max . We also note that the amplitude of the donor is strongly suppressed at resonance (D50), with a comparison with the non-resonance case (D52) illustrated in Figure 3b. This observation also explains the appearance of a peak at resonance. We explain the reason for the observed phenomenon here and leave the physical mechanism for the next section.
In the transfer process, another important quantity is the average transfer time, which is defined as To effectively utilize the energy, the excitation should be transferred with a high efficiency and also within a short time scale. In Figure 4a, we graph the average transfer time as a function of the dissipation rate of the donor k and the detuning D. It can be readily seen in Figure 4a that the optimal frequencies for quick transfer are not at resonance (v5e12g) but split into two, a phenomenon that is known as Rabi splitting in quantum optics. The two peaks can be determined by exactly diagonalizing the photon-assisted donor-acceptor system in the single-excitation subspace as As illustrated in Figure 4b, the transfer time reaches its minimum optimal point when v5v 6 , while the efficiency is not at its maximum. Under these circumstances, finding an optimal parameter for both the short transfer time and the high efficiency is not possible. We will resolve this dilemma by defining a more practical quantity in the following discussion.

THEORETICAL ANALYSIS
The above discovery that the transfer efficiency is insensitive to the dissipation from the donors at resonance can be explained according to the dark state, which has been widely investigated in quantum optics. 29 To this end, we write the dark state in the single-excitation subspace for C50, where h t ð Þã . It has been demonstrated that a perfect transfer of the population between the two low-lying energy levels can be achieved by adiabatically tuning the Rabi frequencies J and j. [30][31][32] In this process, the excitation of the upper energy level is suppressed to avoid dissipation of the excitations.
This dark-state-based mechanism persists in the present artificial system with a bionic structure. For this system, the evolution wave function is where x i is the corresponding eigenvalue of the matrix M, whose expressions are not explicitly written down, as they are too lengthy. The detailed calculation is shown in the supplementary information.
If we choose the parameters as we did previously, e A /j512, g/j50.3, k/j50.1, C/j50.3, e/j511.4 and J/j50.1, then the eigenstates at resonance v5e A are written as which is almost 1 under the practical condition J =j. In the capture and the transfer process, the population of excited donors is suppressed to avoid dissipation, as illustrated in Figure 3b. In the dark state, the system decays very slowly at the rate of c D &C J =j ð Þ 2 , while it dissipates quickly in the bright state ( E 2 j T, E 3 j T) at the approximate rate of c B & kzC ð Þ=2. Thus, for the large time scale t B .1/c B , the majority of the transfer is carried out by the dark state. We graph the ratio of the amplitudes of the photonic state and the acceptor in Figure 5. The asymptotic value of this ratio for long times is approximately j= , which is the value for the dark state in the initial state. This asymptotic behaviour is clearly a confirmation of our present explanation about the intrinsic physical mechanism of the excitation transfer process.

OUTPUT POWER
In the previous discussions, we have proven that the efficiency (proportional to [11(J/j) 2 ] 21/2 ) can be improved via the dark-state mechanism. However, we have also shown that the decay rate of photons at resonance is suppressed by a factor (J/j) 2 . We have reached the dilemma that both the efficiency and the average transfer time cannot be optimized simultaneously.
In fact, such a dilemma has been reached in investigations of heat engines; 33 the Carnot heat engine converts heat into work with maximum efficiency over an infinitely long time. In practice, one would be concerned more with the output power, which characterizes the output energy within a unit time. In the present case, we introduce a similar quantity P5g/t, called the mean transfer power, which characterizes the ability of the excitation energy to be transferred. In Figure 6, we demonstrate the mean power P as a function of the detuning between the incident light frequency and the acceptor excitation energy for different donor dissipation rates k5C, 1.5C and 2C. The optimal frequencies are v59.1j and v514.7j, which are different from the optimal values of both the efficiency and the average transfer time. Actually, it has been observed that the photosynthetic systems chosen absorb sunlight mainly from two domains of the solar spectrum. For example, green plants primarily use red-and blue-coloured photons. The present model can be utilized to justify this observation. However, for a practical system, the estimation of the exact frequencies is beyond the scope of the present model, because it concerns very complicated biological environments.

CONCLUSIONS
We have studied the light capture and excitation transfer process in a generic model consisting of donors and an acceptor. We proved that the present model and the well-studied three-level L-type system are equivalent. The existence of a dark-state channel allows the energy to be transferred efficiently, with the dissipation of donor excitations effectively supressed. We also found two optimal incident frequencies for the maximum output power and discussed their possible connection to the favouritism of two colours in the natural light-harvesting process. In this study, we dealt with the dissipation of an excitation by phenomenologically introducing an imaginary part to the Hamiltonian of a donor. However, in reality, this dissipation is always connected to the vibrational degrees of freedom, which may be accounted for by considering the dimerized structures of the LH I and the LH II, as discussed in Ref. 19. Also, it is worth investigating the effect of the vibrational freedom in a designed artificial system.  Figure 5 The ratio of the amplitudes of the donor mode and the acceptor: | w(t) / v(t) | . This ratio tends to j/(j 2 1J 2 ) 1/2 as time goes to infinity.