electron transfer ( et ) is a very important process in biology , chemistry and physics @xcite .
the most well known et theory is the one of marcus @xcite .
of special interest is the et in configurations where a bridge ( b ) between donor ( d ) and acceptor ( a ) mediates the transfer . on this kind of et we will focus in this paper .
the primary step of et in bacterial photosynthetic reaction centers is of this type @xcite and a lot of work in this direction was done after the structure of the protein - pigment complex of the photosynthetic reaction centers of purple bacteria was clarified in 1984 @xcite .
many artificial systems especially self - organized porphyrin complexes have been developed to model this bacterial photosynthetic reaction center @xcite .
bridge - mediated et reactions can occur via different mechanisms @xcite : incoherent sequential transfer in which the bridge level is populated or coherent superexchange @xcite in which the mediating bridge level is not populated but nevertheless necessary for the transfer . changing
a building block of the complex @xcite or changing the environment @xcite can modify which mechanism is mainly at work .
actually , there is still a discussion in literature whether sequential transfer and superexchange are limiting cases of one process @xcite or whether they are two processes which can coexist @xcite . to clarify
which mechanism is present in an artificial system one can systematically vary the energetics of the complex . in experiments
this is done by substituting parts of the complexes @xcite or by changing the polarity of the solvent @xcite .
also the geometry and size of the bridge block can be varied , and in this way the length of the subsystem through which the electron has to be transfered @xcite can be changed .
superexchange occurs due to coherent mixing of the three or more states of the system @xcite .
the et rate in this channel depends algebraically on the differences between the energy levels @xcite and decreases exponentially with increasing bridge length @xcite .
when incoherent effects such as dephasing dominate the transfer is mainly sequential @xcite , i.  e. , the levels are occupied mainly in sequential order @xcite .
the dependence on the differences between the energy levels is exponential @xcite .
an increase of the bridge length induces only a small reduction in the et rate @xcite .
this is why sequential transfer is the desired process in molecular wires @xcite .    in the superexchange case
the dynamics is mainly hamiltonian and can be described on the basis of the schrdinger equation .
the physically important results can be obtained by perturbation theory @xcite and , most successfully , by the semiclassical marcus theory @xcite .
the complete system dynamics can be directly extracted by numerical diagonalization of the hamiltonian @xcite . in case of sequential transfer
the influence of an environment has to be taken into account .
there are quite a few different ways how to include an environment modeled by a heat bath .
the simplest phenomenological descriptions are based on the einstein coefficients or on the imaginary terms in the hamiltonian @xcite , as well as on the fokker - planck or langevin equations @xcite .
the most accurate but also numerically most expensive way is the path integral method @xcite .
this has been applied to bridge - mediated et especially in the case of bacterial photosynthesis @xcite .
bridge - mediated et has also been investigated using redfield theory @xcite , by propagating a density matrix ( dm ) in liouville space @xcite and other methods ( e.  g.   @xcite ) .
the purpose of the present investigation is to present a simple , analytically solvable model based on the dm formalism @xcite and apply it to a porphyrin - quinone complex which is taken as a model system for the bacterial photosynthetic reaction center .
the master equation which governs the dm evolution as well as the appropriate relaxation coefficients can be derived from such basic informations as system - environment coupling strength and spectral density of the environment @xcite . in the present model relaxation
is introduced in a way similar to redfield theory but in site representation not in eigenstate representation .
a discussion of advantages and disadvantages these representations has been given elsewhere @xcite .
the equations for the dm are the same as in the generalized stochastic liouville equation ( gsle ) model @xcite for exciton transfer which is an extension of the haken - strobl - reineker ( hsr ) model @xcite to a model with a quantum bath . here
we give an analytic solution to these equations .
the present equations for the dm obtained are also similar to those of ref .
@xcite where relaxation is introduced in a phenomenological fashion but only a steady - state solution is found in contrast to the model introduced here .
in addition the present model is applied to a concrete system .
a comparison of the et time with the bath correlation time allows us to regard three time intervals of system dynamics : the interval of memory effects , the dynamical interval , and the kinetic , long - time interval @xcite . in the framework of dm theory
one can describe the et dynamics in all three time intervals . however , often it is enough to find the solution in the kinetic interval for the explanation of experiments within the time resolution of most experimental setups , as has been done in ref .
the master equation is analytically solvable only for simple models , for example @xcite .
most investigations are based on the numerical solution of this equation @xcite . here
we perform numerical as well as approximate analytical calculations for a simple model .
since the solution can be easily obtained , the influence of all parameters on the et can be examined .
the paper is organized as follows . in the next section
we introduce the model of a supermolecule which we use to describe et processes .
the properties of an isolated supermolecule are modeled in subsection [ chem - isolate ] , as well as the static influence of the environment .
the dynamical influence of bath fluctuations is discussed and modeled by a heat bath of harmonic oscillators in section [ chem - sb ] .
the reduced dm equation of motion ( rdmem ) describing the excited state dynamics is presented in subsection [ chem - rdm ] . in subsection
[ chem - scaling ] the system parameter dependence on the solvent dielectric constant is discussed for different models of solute - solvent interaction . in subsection [ chem - parameters ]
system parameters are determined .
the methods and results of the numerical and analytical solutions of the rdmem are presented in section [ chem - results ] .
the dependencies of the et rate and final acceptor population on the system parameters are given for the numerical and analytical solutions in subsection [ chem - super ] .
the analysis of the physical processes in the system is also performed there . in subsection [ chem - solvents ]
we discuss the dependence of the et rate on the solvent dielectric constant for different models of solute - solvent interaction and compare the calculated et rates with the experimentally measured ones .
the advantages and disadvantages of the presented method in comparison with the gsle model @xcite and the method of davis et  al .
@xcite are analyzed in subsection [ chem - davis ] . in the conclusions
the achievements and possible extensions of this work are discussed .
the photoinduced et in supermolecules consisting of three sequentially connected molecular blocks , i.  e. , donor , bridge , and acceptor , ( @xmath4 ) is analyzed .
the donor is not able to transfer its charge directly to acceptor because of their spatial separation .
donor and acceptor can exchange their charges only through b. in the present investigation the supermolecule consists of free - base tetraphenylporphyrin ( @xmath1 ) as donor , zinc substituted tetraphenylporphyrin ( @xmath2 ) as bridge , and p - benzoquinone as acceptor @xcite . in each of those molecular blocks we consider only two molecular orbitals ( @xmath5 ) , the highest occupied molecular orbital ( homo ) and the lowest unoccupied molecular orbital ( lumo ) @xcite .
each of these orbitals can be occupied by an electron or not , denoted by @xmath6 or @xmath7 , respectively .
this model allows us to describe four states of each molecular block , the neutral ground state @xmath8 , the neutral excited state @xmath9 , the positively charged ionic state @xmath10 , and the negatively charged ionic state @xmath11 .
@xmath12 , @xmath13 , and @xmath14 describe the creation , annihilation , and number of electrons in orbital @xmath15 , respectively , while @xmath16 gives the number of electrons in a molecular block .
the number of particles in the whole supermolecule is conserved , @xmath17 .
each of the electronic states has its own vibrational substructure . as a rule for the porphyrin - containing systems
the time of vibrational relaxation is two orders of magnitude faster than the characteristic et time @xcite .
because of this we assume that only the ground vibrational states play a role and we do not include the vibrational substructure . a comparison of the models with and without vibrational substructure has been given elsewhere @xcite .
below we consider the evolution of single charge - transfer exciton states in the system .
for the full description of the system one also should include photon modes to describe for example the fluorescence from the lumo to the homo in each molecular block transferring an excitation to the electro - magnetic field .
but the rates of fluorescence and recombination are small in comparison to other processes for porphyrin - type systems @xcite .
when fluorescence does not have to be taken into account , all states except @xmath18 ( @xmath19 ) , @xmath20 ( @xmath21 ) , and @xmath22 ( @xmath23 ) remain essentially unoccupied , while those three take part in the intermolecular transport process . in this case
the number of states coincides with the number of sites in the system and we label the states @xmath24 instead of @xmath25 .    for the description of the et and other dynamical processes in the system placed in a dissipative environment we introduce the hamiltonian @xmath26 where @xmath27 describes the supermolecule , @xmath28 the dissipative bath , and @xmath29 their interaction .
we are interested in the kinetic limit of the excited state dynamics here . for this limit
we assume that the relaxation of the solvent takes only a very short time compared to the system times of interest .
the effect of the solvent is twofold . on one hand the system dynamics is perturbed by the solvent state fluctuations , independent of the system states .
@xmath29 shall only reflect the dynamical influence of the fluctuations leading to dissipative processes as discussed in the next subsection . on the other hand
the system states are shifted in energy @xcite , @xmath30 due to the static influence of the solvent which is determined by the relaxed value of the solvent polarization and in general also includes the non - electrostatic contributions such as van - der - waals attraction , short - range repulsion , and hydrogen bonding @xcite . in eq .
( [ 2-hami ] ) the energy of free and noninteracting blocks @xmath31 , is given by the energies @xmath32 of orbitals @xmath33 in the independent electron approximation @xcite .
the @xmath32 are chosen to reproduce the ground - state  excited - state transitions e.
g.  @xmath34 , which change only a little for different solvents @xcite and are assumed to be constants here . to determine @xmath32 one starts from fully ionized double bonds in each molecular block @xcite ,
calculates the one - particle states and fills these orbitals with two electrons each starting from the lowest energy . by exciting , removing , adding the last electron to the model system one obtains the energy of the excited , oxidized , reduced molecular block in the independent particle approximation .    the inter - block hopping term @xmath35\ ] ] in eq .
( [ 2-hami ] ) includes the hopping operators @xmath36 and the coherent couplings @xmath37 .
we assume @xmath38 because there is no direct connection between donor and acceptor .
the scaling of @xmath37 for different solvents is discussed in subsection [ chem - scaling ] .    the electrostatic interaction @xmath39 scales like energies of a system of charges in a single or multiple cavity surrounded by a medium with static dielectric constant @xmath40 according to the classical reaction field theory @xcite .
here we consider two models of scaling . in the first model
each molecular block is in an individual cavity in the dielectric . for this case
the electrostatic energy reads @xmath41 @xmath42 takes the electron interaction into account while bringing an additional charge onto the block @xmath43 and thus describes the energy to create an isolated ion .
this term depends on the characteristic radius @xmath44 of the molecular block .
the interaction between the ions @xmath45 depends on the distance between the molecular blocks @xmath46 .
both distances @xmath44 and @xmath46 are also used in marcus theory @xcite .
the term @xmath47 reflects the interaction of charges inside the supermolecule which is weakend by the reaction field according to the born formula @xcite @xmath48    in the second model , considering the supermolecule as one object placed in a single cavity of constant radius one has to use the onsager term @xcite .
this term is state selective , it gives a contribution only for the states with nonzero dipole moment , i.e. , charge separation .
defining the static dipole moment operator as @xmath49 we obtain @xmath50 , with onsager scaling @xmath51      one can express the dynamic part of the system - bath interaction as @xmath52 here @xmath53 denotes the field of the electrostatic displacement at point @xmath54 induced by the system transition dipole moment @xmath55 @xcite . the field of the environmental polarization is denoted as @xmath56 , where @xmath57 is the @xmath58th dipole of the environment and @xmath59 its position .
only fluctuations of the environment polarization @xmath60 influence the system dynamics .
averaged over the angular dependence the interaction reads @xcite @xmath61    the dynamical influence of the solvent is described with a thermal bath model .
the deviation @xmath62 of @xmath63 from its mean value is determined by temperature induced fluctuations . for unpolar solvents described by a set of harmonic oscillators the diagonalization of their interaction yields a bath of harmonic oscillators with different frequencies @xmath64 and effective masses @xmath65 . in the case of a polar solvent the dipoles are interacting rotators as , e.g. used to describe magnetic phenomena @xcite .
the elementary excitation of each frequency can again be characterized by an appropriate harmonic oscillator .
so we use generalized coordinates of solvent harmonic oscillator modes @xmath66 for polar as well as unpolar solvents .
the occupation of the @xmath67th state of the @xmath68th oscillator is defined by the equilibrium dm @xmath69 } \delta_{ij}$ ] .
all mutual orientations and distances of solvent molecules have equal probability .
an average over all spatial configurations is performed .
the interaction hamiltonian  ( [ averaged_value_of_interaction ] ) is written in a form which is bilinear in system and bath operators : @xmath70 \left [ \sum_\lambda k_\lambda ( \hat a^+_\lambda +      \hat a_\lambda ) \right ] s_{\rm sb } \label{phase}\end{aligned}\ ] ] @xmath71 denotes the interaction intensity between the bath mode @xmath72 of frequency @xmath73 and the quantum transition between the lumos of molecules @xmath43 and @xmath74 with frequency @xmath75 .
the scaling function @xmath76 reflects the properties of the solvent .
explicit expressions for the solvent influence are still under discussion in the literature @xcite .
the interaction of the system with the bath of harmonic oscillators describes the irradiative energy transfer from the system to the solvent as modeled by eq .
( [ phase ] ) . for the description of the dynamics we use the reduced dm which can be obtained from the full dm @xmath77 by tracing over the environmental degrees of freedom @xmath78 @xcite with the evolution operator technique @xcite , restricting ourselves to the second order cumulant expansion @xcite . here
we apply the markov approximation , i.e. ,  we restrict ourselves to the limit of long times .
furthermore , we replace the discrete set of bath modes with a continuous one .
to do so one has to introduce the spectral density of bath modes @xmath79 .
finally one obtains the following master equation @xmath80\right)_{\kappa \lambda } + 2\delta_{\kappa \lambda }    \sum\limits_\mu { } \left\ { \gamma_{\mu{}\kappa } \left [ n(\omega_{\mu{}\kappa } ) + 1      \right ] + \gamma_{\kappa \mu { } } n(\omega_{\kappa \mu { } } ) \right\ } \sigma_{\mu{}\mu { } }    \nonumber \\ & & - \sum\limits_\mu { } \left\ { \gamma_{\mu{}\kappa } \left [        n(\omega_{\mu{}\kappa } ) + 1 \right ] + \gamma_{\kappa \mu { } } n(\omega _ { \kappa        \mu { } } ) + \gamma_{\mu{}\lambda } \left [ n(\omega_{\mu{}\lambda } ) + 1 \right ] +      \gamma _ { \lambda \mu { } } n(\omega _ { \lambda \mu { } } ) \right\ } \sigma_{\kappa      \lambda } \nonumber \\ & & + \left\ { \gamma _ { \lambda \kappa } \left [        2n(\omega _ { \lambda \kappa } ) + 1 \right ] + \gamma_{\kappa \lambda }      \left [ 2n(\omega_{\kappa \lambda } ) + 1 \right ] \right\ } \sigma _ { \lambda      \kappa } , \label{rwa - short}\end{aligned}\ ] ] where @xmath81^{-1}$ ] denotes the bose - einstein distribution .
the damping constant @xmath82 reflects the coupling of the transition @xmath83 to a bath mode of the same frequency .
it depends on the density of bath modes @xmath84 at the transition frequency @xmath85 and on the transition dipole moments @xmath86 .
a rdmem of similar structure was used for the description of exciton transfer in the haken , strobl , and reineker ( hsr ) model @xcite and the generalized stochastic liouville equation ( gsle ) model @xcite .
the hsr method originating from the stochastic bath model , is valid only in the high temperature limit @xcite .
the gsle method @xcite appeals to the quantum bath model with system - bath coupling of the form @xmath87 , which modulates the system transition frequency . in ref .
@xcite the equations for exciton motion are derived using the projection operator technique . taking the different system - bath coupling we have derived the rdmem which coincides with gsle @xcite .
both gsle and our rdmem are able to describe correctly finite temperatures .
below we neglect the last term of eq .
( [ rwa - short ] ) corresponding to the @xmath88 term in the hsr and gsle models because the rotating wave approximation ( rwa ) is applied .    for the sake of convenience of analytical and numerical calculations
we replace @xmath89 and the population of the corresponding bath mode @xmath90 with the dissipative transitions @xmath91 and the corresponding dephasings @xmath92 with this one can express the rdmem ( [ rwa - short ] ) in the form @xmath93 the parameters controlling the transitions between the selected states are discussed in subsection  [ chem - parameters ] .
the relaxation coefficients eq .
( [ g = kv ] ) include the second power of the scaling function @xmath76 because one constructs the relaxation term of eq .
( [ rwa - short ] ) with the second power of the interaction hamiltonian .
the physical meaning of @xmath94 is similar to the interaction of the system dipole with a surrounding media .
that is why it is reasonable to use the onsager expression  ( [ onsager - scaling ] ) for @xmath76 . in the work of mataga , kaifu , and koizumi  @xcite the interaction energy between the system dipole and the media scales in leading order as @xmath95,\ ] ] where @xmath96 denotes the optical dielectric constant .
from a recent paper of georgievskii , hsu , and marcus @xcite we extract @xmath97 for the multiple cavities model assuming @xmath98 . in terms of a scaling function
it can be expressed as @xmath99    as we have already argued in @xcite the coherent coupling @xmath37 between two electronic states scales with @xmath40 and @xmath100 too , because a coherent transition in the system is accompanied by a transition of the environment state which is larger for solvents with larger polarity .
as discussed above we neglect the vibrational substructure of each electronic state because the vibrational relaxation is about two orders of magnitude faster than the characteristic et time .
but in contrast to the model with vibrational substructure the present model does not involve any reaction coordinate . to reproduce the results of the more elaborate model with vibrational substructure one has to scale the electronic couplings @xmath37 with the franck - condon overlap elements @xmath101 between the vibrational ground states of each pair of electronic surfaces @xmath102 where @xmath103 is the coupling of electronic states of the isolated molecule .
for the calculation of the franck - condon factors one has to introduce the leading ( mean ) environment oscillator frequency @xmath104 . here
@xmath105 is used which is similar to the frequency of the c - c stretching mode . with this scaling one
implicitly introduces a reaction coordinate into the model .
the dynamics of the system is controlled by the following parameters : energies of system states @xmath106 , coherent couplings @xmath37 , and damping constants @xmath89 .    on the basis of the spectral data @xcite and taking reference energy @xmath107 we determine @xmath108  ev ( in @xmath109 ) .
we take the energy of the state with et to @xmath3 from reference : @xmath110 @xcite .
further rempel et al .
@xcite estimate the coupling of initially excited and charged bridge states @xmath111 and the coupling of the two states with charge separation @xmath112 .
the values of the couplings are essentially lower than the energy differences between the relevant system states @xmath113 this is the reason to remain in site representation instead of eigenstate representation @xcite .
the damping constants are found with help of the analytical solution derived at the end of the next section to be @xmath114 .
the typical radius of the porphyrin ring is about @xmath115 @xcite , while the distance @xmath46 between the blocks of @xmath0 reaches @xmath116 @xcite , @xmath117 , @xmath118 .
the main parameter which controls et in a triad is the energy of the state @xmath119 .
this state has a big dipole moment because of its charge separation and is therefore strongly influenced by the solvent .
because of the special importance of this value we calculate it for the different solvents as a matrix element of the system hamiltonian ( 2 ) .
the calculated values of the energies of the @xmath120 state for some solutions are shown in table  [ tab:2 ] .
the time evolution of the et in the supermolecule is described by solving numerically and analytically eqs .
( [ tosolve1])-([tosolve2 ] ) with the initial condition of donor excitation with a @xmath121 pulse of appropriate frequency , i.e. , the donor population is set to one .    for the numerical simulation
we express the system of eqs .
( [ tosolve1])-([tosolve2 ] ) in the form @xmath122 where @xmath123 is a vector of dimension @xmath124 for the model with @xmath125 system states and the super - operator @xmath126 is a matrix of dimension @xmath127 .
we find an exponential growth of the acceptor population @xmath128,\ ] ] where for the solvent mthf @xmath129 and @xmath130 .
the population @xmath131 which corresponds to charge localization on the bridge does not exceed @xmath132 .
this means that in this case the superexchange mechanism dominates over the sequential transfer mechanism . besides it ensures the validity of characterizing the system dynamics with @xmath133 and @xmath134 dt \right\}^{-1}.\end{aligned}\ ] ]    the alternative analytical approach is performed in the kinetic limit @xmath135 in laplace space the inequality ( [ kinetic - limit ] ) reads @xmath136 , where @xmath137 denotes the laplace variable .
it is equivalent to replacing the factor @xmath138 in the laplace transform of eqs .
( [ tosolve1])-([tosolve2 ] ) with @xmath139 .
this trick allows to substitute the expressions  ( [ tosolve2 ] ) for non - diagonal elements of the dm into eq .
( [ tosolve1 ] ) . after this elimination
we describe the coherent transitions to which the non - diagonal elements contribute by redefinition of the diagonal rdmem  ( [ tosolve1 ] ) @xmath140 the transition coefficients @xmath141 contain dissipative and coherent contributions @xmath142^{-1}.\ ] ] now it is assumed that the bridge is not populated .
this allows us to find the acceptor population in the form of eq .
( [ single - exponent ] ) , where @xmath143^{-1}.\end{aligned}\ ] ] the value of @xmath144 can be found comparing the experimentally determined et rate and eq .
( [ rate ] ) . to calculate
@xmath145 would require a microscopic model . to avoid a microscopic consideration
we simply take the same @xmath89 for all transitions between excited states .
the value of et for @xmath0 in mthf is found by rempel et al .
@xcite to be @xmath146 .
if the bridge state has a rather high energy one can neglect thermally activated processes .
@xmath147 is negligibly small with respect to @xmath148 . in this case
our result  ( [ rate ] ) reads @xmath149 with the relation @xmath150 and the experimental @xmath151 one obtains @xmath152 .
the fit of the numerical solution of eqs .
( [ tosolve1])-([tosolve2 ] ) to the experimental @xmath153 in mthf gives the same value .
so the damping constants are fixed for a specific solvent and for other solvents they are calculated with the scaling functions . with this method
the et was found to occur with dominance of the superexchange mechanism with rates @xmath154 for cyclo and @xmath155 for @xmath156 .
to discuss how the transfer mechanism depends on the change of parameters we calculate the system dynamics varying one parameter at a time .
the dependencies of @xmath153 and @xmath133 on @xmath148 , @xmath147 and @xmath157 , @xmath158 are shown in figs . 2 and 3 .
the change of each parameter influences the transfer in a different way .
in particular , @xmath153 depends quadratically on @xmath148 from @xmath159 to @xmath160 in fig .  2 .
below it saturates at the lower bound @xmath161 .
this corresponds to a crossover of the et mechanism from superexchange to sequential transfer .
but , due to the big energy difference between donor and bridge states the sequential transfer efficiency is extremely low .
this is displayed by @xmath162 . in the region @xmath163 both mechanisms contribute to @xmath153 .
the decrease of @xmath133 in this region corresponds to coherent back transfer .
the et rate depends on @xmath147 in a similar way . at rather high values of @xmath148 , @xmath164 the relation  ( [ inequality_w > v ] )
is no more valid
. for this regime one has to use eigenstate instead of site representation because the wavefunctions are no more localized @xcite .
the variation @xmath157 , @xmath158 near the experimental values shows similar behavior of @xmath165 and @xmath166 ( see fig .
we independently vary @xmath157 and @xmath158 .
both , @xmath165 and @xmath166 increase linear until the saturation value @xmath167 at @xmath168 is reached .
there is qualitative agreement between the numerical and analytical values . in eq .
( [ fit ] ) infinite time is approximated by @xmath169 and so one can not obtain et rates lower than this limit .    the physical meaning of the et rate dependence on @xmath170 seems to be transparent . at small values of @xmath170 a part of the population coherently
oscillates back and forth between the states .
the increase of the dephasing @xmath171 quenches the coherence and makes the transfer irreversible .
so transfer becomes faster up to a maximal value .
for the whole range of @xmath170 , depopulations @xmath172 , @xmath173 and thermally activated transitions @xmath174 , @xmath175 always remain smaller than the coherent couplings , therefore they do not play an essential role .
next , the similarity of the dependencies on @xmath157 and @xmath158 will be discussed on the basis of eq .
( [ rate ] ) . in the limit
@xmath176 thermally activated processes with @xmath177 vanish and so @xmath178 , while depopulations with @xmath179 remain constant @xmath180 .
the condition @xmath181 allows us to neglect @xmath182 in comparison with @xmath183 . with these simplifications eq .
( [ rate ] ) becomes @xmath184 i.e.  symmetric with respect to @xmath157 and @xmath158 .    to the largest extent
the mechanism of transfer depends on the bridge energy @xmath185 as presented in fig .
4 . in different regions
one observes different types of dynamics . for large bridge energy @xmath186
the numerical and analytical solutions do not differ from each other .
the transfer occurs with the superexchange mechanism .
the et rate reaches a maximal value of @xmath187 for low bridge energies .
while the bridge energy approaches the donor energy the sequential transfer starts to contribute to the et process .
the traditional scheme of sequential transfer is obtained when donor , bridge , and acceptor levels are arranged in a cascade . in this region
the analytical solution need not coincide with the numerical solution because the used approximations are no more valid .
for equal bridge and acceptor energies @xmath153 displays a small resonance peak in fig .
when the bridge energy is lower than the acceptor energy the population gets trapped at the bridge .
the finite @xmath153 for @xmath188 does not mean et because @xmath189 .
for the dynamic time interval @xmath190 a part of the population tunnels force and back to the acceptor with @xmath153 .
the analytical solution  ( [ rate ] ) gives a constant rate for the regime @xmath191 , while the numerical solution of eqs .
( [ tosolve1])-([tosolve2 ] ) is instable , because such coherent oscillations of population can not be described by eq .
( [ single - exponent ] ) and @xmath153 can not be fitted with eq .
( [ fit ] ) . in fig .
4 the regime @xmath191 occurs for small @xmath192 while @xmath193 is kept constant and for large @xmath193 while @xmath192 remains constant .
the energy dependence of the final population has a transparent physical meaning for the whole range of energy .
a large bridge energy ensures the transition of the whole population to the acceptor . in the intermediate case , when the bridge has the same energy as the acceptor , final population spreads itself over these two states @xmath194 . lowering the bridge
even more the whole population remains on the bridge as the lowest state of the system .
the dependence of the et rate on the acceptor energy @xmath195 in fig .
4 remains constant while the acceptor energy lies below the bridge energy .
increase of @xmath193 up to @xmath196 gives the maximal @xmath197 .
when @xmath193 increases further the acceptor becomes the highest level in the system and therefore the population can not remain on it .      for the application of the results to various solvents and comparison with experiment one should use the scaling for energy , coherent couplings , and damping constants as discussed above .
the combinations of the energy scaling in subsection [ chem - isolate ] and damping constants scalings in subsection  [ chem - scaling ] are represented in fig .
an increase in @xmath40 from @xmath198 to @xmath199 leads to an increase of the et rate , no matter which scaling is used .
further increase of @xmath40 induces saturation for the onsager - mataga scaling and even a small decrease within the applied approximations an increase in the solvent polarizability and , consequently , of its dielectric constant lowers the bridge and acceptor energies and increases the system - bath interaction and , consequently , the relaxation coefficients .
it induces a smooth rise of the et rate for the onsager - mataga scaling . on the other hand large @xmath40 leads to essentially different polarisational states of the environment for the supermolecule states with different dipole moment
this reduces the coherent couplings , see eq .
( [ v_scale ] ) , leading for the born - marcus scaling to a small decrease of @xmath153 for large @xmath40 .
the et rate with this scaling comes closer to the experimental value @xmath200 .
this gives a hint that the model of individual cavities for each molecular block is closer to reality than the model with a single cavity for the whole supermolecule .
below we consider born scaling eq .
( [ born - scaling ] ) for the system energies and marcus scaling eq .
( [ marcus - scaling ] ) for the damping constant to compare the calculated et rates with the measured ones .
for the solvents cyclo , mthf , and @xmath156 one obtains the relative bridge energies @xmath201 , @xmath202 , and @xmath203 , respectively .
the calculated et rate coincides with the experimental value @xcite for @xmath0 in cyclo , see table  [ tab:2 ] . for @xmath156
the numerical et rate is approximately 30% faster than the experimental value .
it has to be noted that a value for the damping rates can be chosen such that the calculated curve almost passes through all three experimental error bars .
on the other hand an error in the present calculation could be due to ( i ) absence of vibrational substructure of the electronic states in the present model ; ( ii ) incorrect dependence of system states energies on the solvent properties ; ( iii ) opening of additional transfer channels not mentioned in the scheme shown in fig .
each of these possibilities needs some comments .    _
( i ) : the incorporation of the vibrational substructure will result in a complication of the model @xcite giving a more complicated et rate dependence on the energy of the electronic states and dielectric constant .
it should yield the maximal et rate for nonequal energies of electronic states , namely for the activationless case when the energy difference equals the reorganization energy . for a comparison of the models with and without vibrational substructure
see ref .
@xcite .
_ ad _ ( ii ) : effects as the solvation shell @xcite do need a molecular dynamics simulation .
the total influence of the solvent is , probably , reflected in an energy shift between the spectroscopically observable states @xmath204 and @xmath205 @xcite .    _ ad _ ( iii ) : a solvent with large @xmath40 can bring high - lying system states closer to the ones included in fig .  1 .
e.g. ,  because of its larger dipole moment , @xmath206 is strongly influenced by the solvent .      as discussed above the rdmem
are very similar to those of the gsle model .
this is an extension of the hsr theory in which a classical bath is used and for which analytical solutions are available @xcite .
we are not aware of any analytical solution of the gsle model as presented here .
also this model has not been applied to similar et processes .
the numerical steady - state method used by davis et al .
@xcite is an attractive one due to its simplicity , but unlike our method it is not able to give information about the time evolution of the system .
we use a similar approach derived within a redfield - like theory .
but we consider dephasing and depopulation between each pair of levels .
in contrast davis et al .
incorporate relaxation phenomenologically only to selected levels , dephasing @xmath207 occurs between excited levels , while depopulation @xmath208 takes place only for the sink from acceptor to the ground state . the advantage of the approach of davis et al .
is the possibility to investigate the et rate dependence for the bridge consisting of more than one molecular block .
this was not the goal of the present work but it can be extended into this direction .
we are interested in the et in a concrete molecular complex with realistic parameter values and realistic possibilities to modify those parameter .
our results as well as the results of davis show that et can occur as coherent ( with the superexchange mechanism ) or dissipative process ( with the sequential transfer mechanism ) .
we have performed a study of the et in the supermolecular complex @xmath209 within the dm formalism .
the determined analytical and numerical et rates are in reasonable correspondence with the experimental data .
the superexchange mechanism of et dominates over the sequential transfer .
we have investigated the stability of the model varying one parameter at a time .
the qualitative character of the transfer is stable with respect to a local change of system parameter .
the crossover between the two transfer mechanisms can be induced by lowering the bridge energy .
the relation of the theory presented here to other theoretical approaches to et has been discussed .
the calculations performed in the framework of the present formalism can be extended in the following directions : ( i ) considerations beyond the kinetic limit .
the vibrational substructure has to be included into the model as well as solvent dynamics and , probably , non - markovian rdmem .
( ii ) enlargement of the number of molecular blocks in the complex .
( iii ) initial excitation of states with rather high energy should open additional transfer channels .
j.  jortner and m.  bixon ( eds . ) , adv .
* 106&107 * ( 1999 ) .
d.  devault , _ quantum mechanical tunneling in biological systems _
( cambridge university press , cambridge , 1993 ) . v.  balzani and f.  scandola , _ supramolecular photochemistry _ ( ellis horwood , chichester , 1991 )
. m.  newton , chem . rev . * 91 * , 767 ( 1991 ) .
p.  f. barbara , t.  j. meyer , and m.  a.  ratner , j. phys .
chem . * 100 * , 13148 ( 1996 ) .
r.  a.  marcus , j.  chem .  phys . * 24 * , 966 ( 1956 ) ; rev .  mod
. phys .  * 65 * , 599 ( 1993 ) ; r.  a.  marcus and n.  sutin , biochim .
biophys .
acta * 811 * , 265 ( 1985 ) .
m.  bixon , j.  jortner , and m.  e.  michel - beyerle , biochim .
biophys .
acta * 1056 * , 301 ( 1991 ) ; chem .  phys . ,
389 ( 1995 ) .
j. deisenhofer , o. epp , k. miki , r. huber and h. michel , j. mol .
biol .  * 180 * , 385 ( 1984 ) .
m.  r. wasielewski , chem . rev . * 92 * , 345 ( 1992 ) ; m.  r. wasielewski , d.  g.  johnson , w.  a.  svec , k.  m.  kersey , d.  e.  cragg , and d.  w.  minsek , in : _ photochemical energy conversion _ , eds . j.  norris and d.  meisel ( elsevier , amsterdam , 1989 ) p.  135 ; m.  r. wasielewski , m.  p.  niemczyk , d.  g.  johnson , w.  a.  svec , and d.  w.  minsek , _ tetrahedron _ * 45 * , 4785 ( 1989 ) .
d.  g.  johnson , m.  p.  niemczyk , d.  w.  minsek , g.  p.  wiederrecht , w.  a.  svec , g.  l.  gaines iii , and m.  r.  wasielevski , j. am .
soc . * 115 * , 5692 ( 1993 ) .
w.  davis , m.  wasielewski ,  m. ratner , v.  mujica , and a. nitzan , j. phys . chem . * 101 * , 6158 ( 1997 ) . s.  s. scourtis and s.  mukamel , chem . phys . * 197 * , 367 ( 1995 ) .
m.  schreiber , c.  fuchs , and r.  scholz , j. lumin .
* 76&77 * , 482 ( 1998 ) . h.  m. mcconnel , j. chem . phys . * 35 * , 508 ( 1961 ) .
r.  kosloff and m.  a. ratner , isr . j. chem . *
30 * , 45 ( 1990 ) .
u.  rempel , b.  von  maltzan , and c.  von  borczyskowski , chem .
* 245 * , 253 ( 1995 ) .
e.  i.  zenkevich , v.  n.  knyukshto , a.  m.  shulga , v.  a.  kuzmitsky , v.  i.  gael , e.  g.  levinson , and a.  f.  mironov , j. lumin . * 75 * , 229 ( 1997 ) ; e.  i.  zenkevich , a.  m.  shulga , s.  m.  bachilo , u.  rempel , j.  von  richthofen , and c.  von  borczyskowski , j. lumin . *
76&77 * , 354 ( 1998 ) ; a.  chernook , u.  rempel , c.  von  borczyskowski , a.  m.  shulga , and e.  i.  zenkevich , chem .
. lett . * 254 * , 229 ( 1996 ) ; j. phys .
chem . * 100 * , 1918 ( 1996 ) ; ber . bunsenges .
phys . chem . *
100 * , 2065 ( 1996 ) .
a.  i. burstein and y.  georgievskii , j. chem . phys . *
100 * , 7319 ( 1994 ) .
h.  sumi and t.  kakitani , chem .
lett .   * 252 * , 85(1996 ) ; h.  sumi , j. electroan .
* 438 * , 11 ( 1997 ) .
b.  paulson , k.  pramod , p.  eaton , g.  l.  closs , and j.  r.  miller , j. chem .
phys . * 97 * , 13042 ( 1993 ) .
n.  s. hush , coord .
chem . rev . * 64 * , 135 ( 1985 ) .
r.  langen i.  chang , j.  p.  germanas , j.  h.  richards , j.  r.  winkler , and h.  b.  gray , science * 268 * , 1733 ( 1995 ) ; c.  j. murphy , m.  r.  arkin , y.  jenkins , n.  d.  ghatlia , s.  h.  bossman , n.  j.  turro , and j.  k.  barton , science * 262 * , 1025 ( 1993 ) .
v.  mujica , m.  kemp , and m.  a. ratner , j. chem .
* 101 * , 6849 ( 1994 ) ; * 101 * , 6856 ( 1994 ) ; m.  kemp , v.  mujica , and m.  a. ratner , j. chem
. phys . * 101 * , 5172 ( 1994 ) ; v.  mujica , m.  kemp , a.  roitberg , and m.  a.  ratner , j. chem .
phys . * 104 * , 7296 ( 1996 ) ; m.  p. samanta , w.  tian , s. datta , j.  i.  henderson , and c.  p.  kubiak , phys . rev .
b * 53 * , r7626 ( 1996 ) .
m.  a. ratner , j. phys .
* 94 * , 4877 ( 1990 ) ; r.  j. miller and j.  v. beitz , j. chem . phys . * 74 * , 6749 ( 1981 ) . j.  w. evenson and m.  karplus , j. chem
96 * , 5272 ( 1992 ) ; science , * 262 * , 1247 ( 1993 ) .
w.  t. pollard , a.  k. felts , and r.  a. friesner , adv .
93 * , 77 ( 1996 ) ; a.  k. felts , w.  t. pollard , and r.  a.  friesner , j. phys . chem . * 99 * , 2929 ( 1995 ) .
w. b. davis , w. a. svec , m. a. ratner , and m. r. wasielewski , nature * 396 * , 60 ( 1998 ) .
r.  cave and m.  d. newton , chem .
* 249 * , 15 ( 1996 ) .
m. bixon and j. jortner , j. chem . phys . *
107 * , 5154 ( 1997 ) .
u.  weiss , _ quantum dissipative systems _ ( world scientific , singapore , 1999 ) , 2nd ed .
r.  loudon , _ the quantum theory of light _ ( clarendon , oxford , 1973 ) ; j.  r.  klauder and e.  c.  g.  sudarshan , _ fundamentals of quantum optics _
( benjamin , new york , 1968 ) r.  egger , c.  h.  mak , and u.  weiss , phys .
e * 50 * , r655 ( 1994 ) ; c.  h.  mak and r.  egger , chem .
lett .  * 238 * , 149 ( 1995 ) ; n.  makri , e.  sim , e.  makarov , and m.  topaler , proc .
usa * 93 * , 3926 ( 1996 ) ; e.  sim and n.  makri , j.  phys .  chem .
b * 101 * , 5446 ( 1997 ) .
a.  k.  felts , w.  t.  pollard , and r.  a.  friesner , j. phys .
chem .  * 99 * , 2929 ( 1995 ) . s. h. lin , m. hayashi , s. suzuki , x. gu , w. xiao , and m. sugawara , chem .
phys .  * 197 * , 435 ( 1995 ) .
o. khn , v. rupasov , and s. mukamel , j. chem .
, 5821 ( 1996 ) .
k. blum , _ density matrix theory and applications _ ( plenum , new york , 1996 ) , 2nd ed .
a.  g. redfield , _
ibm j. res .
_ * 1 * , 19 ( 1957 ) ; adv .
. reson . * 1 * , 1 ( 1965 ) . v.  may and m.  schreiber , phys .
a * 45 * , 2868 ( 1992 ) .
o.  khn , v.  may , and m.  schreiber , j. chem . phys . * 101 * , 10404 ( 1994 ) . j.  m. jean , r.  a. friesner , and g.  r. fleming , j. chem
. phys . * 96 * , 5827 ( 1992 ) ; w.  t. pollard and r.  a. friesner , j. chem . phys . * 100 * , 5054 ( 1994 ) ; j.  m. jean and g.  r. fleming , j. chem . phys .
* 103 * , 2092 ( 1995 ) ; j.  m. jean , j. chem .
phys . * 104 * , 5638 ( 1996 ) .
a.  nitzan , chem .
* 41 * , 163 ( 1979 ) ; r.  kosloff and s.  a. rice , j. chem .
phys . * 72 * , 4591 ( 1980 ) .
r.  kosloff , m.  a. ratner , and w.  b. davis , j. chem .
phys . * 106 * , 7036 ( 1997 ) .
m.  schreiber and d.  kilin , in : _ proc .
excitonic processes in condensed matter _ , ed
. m. schreiber ( dresden university press , dresden , 1996 ) p. 331
, http://xxx.lanl.gov/abs/cond-mat/0001061 .
w.  b.  davis , m.  r.  wasielewski , r.  kosloff , and m.  a.  ratner , j. phys .
a * 102 * 9360 ( 1998 ) .
v. kenkre and p. reineker , _ exciton dynamics in molecular crystals and aggregates _
( springer , berlin , 1982 ) .
h. haken , p. reineker , z. phys .
* 249 * , 253 ( 1972 ) .
d.  a.  weitz , s.  garoff , j.  i.  gersten , and a.  nitzan , j. chem .
phys . * 78 * , 5324 ( 1983 ) .
d.  kilin and m.  schreiber , j. lumin .
* 76&77 * 433
( 1998 ) , http://xxx.lanl.gov/abs/cond-mat/0001052 .
m.  gouterman , j. mol . spectrosc .
* 6 * , 138 ( 1961 ) .
m. schreiber , d. kilin , and u. kleinekathfer , j. lumin .
* 83&84 * , 235 ( 1999 ) , htp://xxx.lanl.gov / abs / quant - ph/9912068 .
d.  d.  fraser , j.  r.  bolton , j. phys .
chem . * 98 * , 1626 ( 1994 ) .
v.  m.  agranovich and m.  d.  galanin , _ electronic excitation energy transfer in condensed matter _
( north - holland , amsterdam , 1982 ) ; l.  d.  landau , e.  m.  lifshitz , _ lehrbuch der theoretischen physik , band viii _ , ( akademie - verlag , berlin 1985 ) .
y.  georgievskii , c .- p .
hsu , r.  a.  marcus , j. chem .
phys . * 110 * , 5307 ( 1999 ) .
o.  christiansen and k.  mikkelsen , j.  chem . phys . * 110 * , 8348 ( 1999 ) .
r.  p.  feynman , _ the feynman lectures on physics vol .
iii _ , ( addison - wesley , reading , 1963 ) . c.  j.  f.  bttcher , _ theory of electric polarization _
( elsevier , amsterdam , 1973 ) vol . 1 .
m.  karelson , g.  h.  f.  diercksen , in : _ problem solving in computational molecular science : molecules in different environments _ , eds .
s.  wilson and g.  h.  f.  diercksen ( kluwer , dordrecht , 1997 ) p.  195 .
k.  yosida , _ theory of magnetism _
( springer , berlin , 1996 ) .
s.  v.  tyablikov , _ methods of the quantum theory of magnetism _ ( plenum , new york , 1967 ) .
a.  apanasevich , s.  ya .
kilin , a.  p.  nizovtsev , and n.  s.  onishchenko , j. opt .
b * 3 * , 587 ( 1986 )
. b. yoon , j.m .
deutch , and j.h .
freed , j. chem .
phys . * 62 * , 4687 ( 1975 ) ; a.  isihara , _ statistical physics _ , ( academic press , new york , 1971 ) .
.energy of the charged bridge state in different solvents and corresponding et rates ( for calculations born scaling  ( [ born - scaling ] ) of energy and marcus scaling  ( [ marcus - scaling ] ) of dissipation are used ) .
mthf denotes 2-methyl - tetrahydrofuran , and cyclo denotes cyclohexane . [ cols="^,^,^,^",options="header " , ]