molecular junctions are electronic devices that consist of an organic quantum dot that is contacted by two macroscopic leads .
modern nanotechnology allows for the reliable fabrication of systems where the dots are single aromatic rings , molecular wires@xcite , c@xmath0 fullerenes@xcite or carbon nanotubes.@xcite they are a promising candidates in the search for further miniaturization of electronic and thermoelectric devices.@xcite    transport through such systems is determined by the discrete levels of the dot , whose position relative to the fermi energy can be tuned , e.g. , with the help of a third ( gate ) electrode .
@xcite the level broadening depends on the dot - lead coupling strength , which can be manipulated via the lead distance or through the choice of different metal - molecule anchoring groups .
@xcite in addition to these two external parameters , electron - phonon ( ep ) interaction influences transport through molecular junctions : when it is occupied by charge carriers , the molecule may undergo structural changes or vibrations that correspond to the excitation of local optical phonons of considerable energy .
they show up as vibrational signatures in the current - voltage characteristics of the device  @xcite .
a temperature difference between the leads induces a current of charge carriers across the junction .
this thermoelectric effect is measured by recording , for constant temperature difference , the voltage bias necessary to cancel this current .
@xcite the quotient of the temperature difference and the thermovoltage , the so - called thermopower , can be used to probe the systems fermi energy@xcite and the vibronic structure of the molecule s state.@xcite
most of the experimental findings are well understood within a linear response formulation of the thermopower.@xcite however , the applied temperature differences can be tens of degrees kelvin , i.e. larger than the dot - lead coupling energy .
some features of the measured voltage histograms , such as side peaks and temperature dependent widths,@xcite are not accounted for in the linear theory .
that is why recently , the discussion of the thermopower has been extended to the nonlinear regime.@xcite then the question arises , how the thermoelectric coefficients can be generalized .
our approach in the present work is motivated by the experimental situation : for a given temperature bias we determine the thermovoltage numerically by minimizing the thermally induced charge current .
we base our calculations on the anderson - holstein model . here , the organic molecule is represented by a single energy level and a local optical mode , which is linearly coupled to the electron on this level .
the quantum dot is connected to two macroscopic leads , while the local mode is coupled to a phonon bath .
the current between the leads is given by the interacting dot density of states.@xcite based on such models , different methods have been applied , such as the numerical@xcite and functional renormalization group@xcite , rate equations,@xcite master equations@xcite and nonequilibrium keldysh green functions@xcite to describe transport through the dot for small - to - large dot - lead coupling and weak - to - strong ep interaction . in the antiadiabatic , strong ep coupling regime ,
typically a lang - firsov transformation@xcite is applied , based on the exact solution of the isolated dot .
it predicts the formation of a local polaron , which reduces the effective dot - lead coupling exponentially . this could be beneficial for the thermoelectric response of the system .
@xcite for practical applications , however , a moderate level broadening is needed to ensure usable power output .
moreover , long electron residence times and strong ep interaction may lead to the accumulation of energy at the dot and , consequently , to its degradation.@xcite because of this , the regime of comparable electronic and phononic time scales and intermediate ep interaction becomes interesting .    to account for the polaronic character of the dot - state away from the strong ep coupling , antiadiabatic limit
, we use an approach that is based on a variational form of the lang - firsov transformation.@xcite for the polaron problem , the variational lang - firsov approach has been proven to give reliable results in the whole electron - phonon coupling and phonon frequency regime , even in the most physically difficult polaron crossover region.@xcite here , the polaron variational parameter is determined by minimizing the relevant thermodynamical potential .    in previous work ,  @xcite we considered the steady - state current response of the quantum dot to a finite voltage difference between the leads . in doing so
, we assumed that the temperatures of the leads and of the phonon bath are all equal .
the calculations in these papers were based on the kadanoff - baym formalism ,  @xcite which relies on the relation between the real - time response functions and the non - equilibrium green functions of the complex time variables .
the present work will focus on the thermoelectric effects induced by a finite temperature difference between the leads .
section  [ secmodel ] introduces the model and the variational ansatz we employ to describe a vibrating quantum dot that is coupled to a phonon heat bath and two macroscopic leads at different temperatures .
because for such a setup the temperature is not constant throughout the system , the reasonings of our earlier approach  @xcite have to be modified .
that is why in sec .
[ secresponse ] we generalize the kadanoff - baym method such that the steady - state equation of the response functions , as well as their formal solution deduced in ref .  , are applicable to the present case . in sec .
[ secself ] we then derive an approximation for the polaronic self - energy , and sec .  [ seccurrent ] provides the relation between the polaronic and electronic spectral functions , the latter of which enters the current formula . in sec .
[ secvar ] we derive the thermodynamic potential which will be used to determine the variational parameter numerically .
section  [ secnumerics ] presents our numerical results .
the main conclusions and future prospects can be found in sec .
[ secsummary ] .
is coupled to two macroscopic leads , each in its own thermal equilibrium with different temperatures @xmath1 and chemical potentials @xmath2 .
the dot electrons interact with an optical phonon mode of energy @xmath3 that is coupled to a phonon bath at temperature @xmath4 . ]
we consider a molecular quantum dot in the three - terminal configuration depicted in fig .
[ fig1 ] , which will be modelled by the following hamiltonian : @xmath5 here , the quantum dot is represented by a single energy level @xmath6 with the fermionic operators @xmath7 . it interacts with a local optical phonon mode @xmath8 of energy @xmath3 , where the ep coupling strength is given by the so - called polaron binding energy , @xmath9 by the last term in eq .
( [ equhamiltonianstart ] ) , the dot is coupled to left ( @xmath10 ) and right ( @xmath11 ) leads , each of which contains @xmath12 free electrons with the energies @xmath13 and the corresponding fermionic operators @xmath14 , respectively . in equilibrium ,
the dot - lead system is characterized by a common chemical potential @xmath15 and a global temperature @xmath16 .    in the nonequilibrium situation ,
a voltage difference between the leads is described by adding to ( [ equhamiltonian ] ) the term @xmath17 with the voltage bias @xmath18 where @xmath19 is the electron charge .
in addition , we consider a temperature difference between the macroscopic leads , whereby each lead is supposed to stay in its own thermal equilibrium with the temperatures @xmath20 and @xmath21 , respectively .
we also suppose that the local oscillator is coupled to its own heat bath , which has the temperature @xmath22 . in accordance with entin - wohlman  @xcite
we assume that the coupling to this heat bath ( indicated as @xmath23 in fig .
[ fig1 ] ) far exceeds the ep - coupling strength , so that the phonon population at the dot is given by the bose - einstein distribution @xmath24 in the following calculations , we suppose that @xmath25 and @xmath26 in general .
moreover , all energies will be taken with respect to the equilibrium chemical potential , i.e. , @xmath27 .
we raise the voltage bias symmetrically around @xmath15 , i.e. , @xmath28 and @xmath29 . to account for the polaron localization at the dot
, we apply to the hamiltonian ( [ equhamiltonianstart ] ) an incomplete lang - firsov transformation @xmath30 ,  @xcite introducing the variational parameter @xmath31 $ ] : @xmath32 the transformed hamiltonian then reads @xmath33 where the renormalization of the dot level and the interaction coefficients depends on the parameter @xmath34 : @xmath35 now @xmath36 and @xmath37 are the operators of dressed electrons ( in analogy to polarons ) and the shifted local oscillator .
the original electron and oscillator operators read @xmath38 and @xmath39 , respectively .      in our previous works ,
@xcite we considered the electron current response of a molecular quantum dot that was initially in equilibrium with the leads as well as with the phononic bath at a common , fixed temperature .
our calculations were based on the kadanoff - baym theory .
@xcite at first , we summarize the aspects of this non - equilibrium response theory that are essential for the present setting .
the relevant response functions are represented by the real - time green functions of the dot - operators @xmath7 and the lead operators @xmath14 , e.g. @xmath40 the functions @xmath41 of the `` mixed '' operators @xmath14 and @xmath7 are defined in an analogous way . in eqs.([equdefresponse])([equdefresponsegtr ] ) the time dependence of @xmath42 is determined by @xmath43 .
the symbol @xmath44 means the standard time - ordering operator so that the function @xmath45 is equal to the functions @xmath46 and @xmath47 for @xmath48 and @xmath49 , respectively .
the statistical average @xmath50 corresponds to the equilibrium state at the temperature @xmath16 before the disturbance was turned on .
going to the interaction representation , the heisenberg operators @xmath51 are expressed as @xmath52 the real - time response functions may be deduced using the equations of motion for the nonequilibrium green functions of the complex time variables @xmath53 , @xmath54 $ ] , defined as @xmath55 the time dependence of the operators in eq .
( [ equdefg ] ) is determined by @xmath56 only , while the external disturbance is explicit in the time - ordered exponential operator @xmath57 .
the operation @xmath58 orders the operators by the imaginary parts of the times @xmath59 and @xmath60 , so that @xmath61 for @xmath62 and @xmath63 for @xmath64 . to find the relation between the functions @xmath65 and @xmath66 , the function @xmath67 is considered for @xmath64 : @xmath68 where @xmath69 the continuation of @xmath70 and @xmath71 to analytic functions of the time variables leads to the identification of @xmath70 with @xmath71 in the limit @xmath72
consequently , the connection of the analytic functions @xmath73 and @xmath74 is given by @xmath75 with similar relations for the functions @xmath76 and @xmath77 .
it is evident that the derivation of eq .
( [ equlimitlessgtr ] ) outlined above does not refer to some special properties of the statistical ensemble , nor to the physical meaning of @xmath78 . in this way , it is possible to generalize the definition ( [ equdefg ] ) for @xmath67 assuming the mean value @xmath79 to be unspecified and the complex - time variable to be defined in the interval @xmath80 $ ] , where the time - ordering parameter @xmath81 has no specific physical meaning .
we assume that before the disturbance ( [ equhint ] ) was turned on , the system was in a steady state with the temperatures @xmath20 , @xmath21 of the left and right leads , and @xmath4 of the phonon bath .
the function @xmath67 defined in this way does not have the properties of the temperature ( @xmath78 ) dependent green function in eq .
( [ equdefg ] ) , it rather represents a functional of the ordered operators which we use to determine the real - time response funcitons @xmath45 .    to do this , we define , in analogy to the self energy , the function @xmath82 by the equation @xmath83\bullet & g_{dd }   ( \bar t , t_2;u , t_0)\nonumber \\ & = \delta(t_1-t_2)\;,\label{equdyson}\end{aligned}\ ] ] with the inverse zeroth - order function @xmath84 in eq .
( [ equdyson ] ) the matrix multiplication `` @xmath85 '' is defined by the integration @xmath86 containing the time - ordering parameter @xmath81 . the @xmath87-function of complex arguments is understood with respect to this integration .
the inverse function to @xmath67 is given as @xmath88 .
the deduction of the steady - state equations for the real - time functions @xmath89 and @xmath90 by the limiting procedure @xmath72 is analogous to that given in ref .  .
defining the fourier transformations according to kadanoff - baym,@xcite @xmath91 and similarly for @xmath92 ) , the solution of the steady - state equations may be written as follows : @xmath93 with the non - equilibrium polaronic spectral function @xmath94      according to the preceeding section , the concrete form of the steady - state solution for the special choice of interactions is determined by the functions @xmath95 . to find an explicit expression for @xmath95 , we start with the equations of motion for @xmath67 and @xmath96 , which are given by the commutators of the operators @xmath7 and @xmath14 with @xmath56 . as a purely formal device , we add to @xmath97 in eq .
( [ equhint ] ) the interaction with fictitious external fields @xmath98 .
the equations of motion of @xmath67 and @xmath96 are then expressed by means of the functional derivatives of @xmath99 with respect to these fields .
the resulting equations for @xmath95 are solved iteratively .
we then let @xmath100 and perform the limit @xmath101 . in the following calculations we will use the self - energy function after the first iteration step:@xcite @xmath102 \big \}\ ; , \label{equsigma1lessgtr}\end{aligned}\ ] ] where we have defined @xmath103 the bose - function @xmath104 in ( [ equsigma1lessgtr ] )
contains the phonon - bath temperature @xmath4 according to eq .
( [ equbosebath ] ) , while the zeroth - order functions of the leads depend on the different temperatures @xmath21 and @xmath20 : @xmath105\;,\end{aligned}\ ] ] with the lead fermi - functions @xmath106 the fourier - transformation of eq .
( [ equsigma1lessgtr ] ) leads to @xmath107 \big \ } \ ; , \nonumber\end{aligned}\ ] ] with the renormalized dot - lead coupling function @xmath108 which depends on the lead density of states : @xmath109 from eq .
( [ equsigmalessfourier1 ] ) , @xmath110 results by interchanging @xmath111 and @xmath112 .
note that eq .
( [ equsigmalessfourier1 ] ) is similar to our result for the self - energy in ref .  , but now @xmath113 is replaced by the individual lead fermi - functions @xmath114 that were defined in ( [ equfermileads ] ) . if we insert the approximation ( [ equsigmalessfourier1 ] ) into eq .
( [ equdefgamma ] ) , we find @xmath115 \gamma_{a}^{(1)}&(\omega;u)= i_{0}(\kappa )   \widetilde\gamma^{(0)}_a(\omega)\label{equgamma1a}\\ & + \sum_{s\ge 1 } i_{s}(\kappa)2\sinh(s\theta)\nonumber\\ & \times \big [   \widetilde\gamma^{(0)}_a(\omega - s\omega_0)\big(n_b(s\omega_0)+1-f_a(\omega+u_a - s\omega_0)\big ) \nonumber\\ & + \widetilde\gamma^{(0)}_a(\omega+s\omega_0)\big ( n_b(s\omega_0)+f_a(\omega+u_a+s\omega_0 ) \big)\big ] \;.\nonumber\end{aligned}\ ] ] from the functions @xmath116 and @xmath117 , the first order spectral function @xmath118 , the distribution function @xmath119 and the response functions @xmath120 follow according to eqs .
( [ equafrac ] ) , ( [ equbarf ] ) and ( [ equgless ] ) , respectively .
the operator of the particle current from lead @xmath121 to the dot reads @xmath122\;.\label{equcurrentop}\end{aligned}\ ] ] its mean value @xmath123 is given by the real - time response functions @xmath124 of the electron operators , which are defined in analogy to eqs .
( [ equdefresponseless ] ) and ( [ equdefresponsegtr ] ) , e.g. @xmath125 we determine @xmath124 based on the equation of motion of the corresponding function @xmath126 of complex - time variables : @xmath127 the current @xmath128 results as @xmath129 where @xmath130 and @xmath131 are the electronic dot spectral function and response function . in the steady state regime
our approximation conserves the particle current : @xmath132 . in the following , identical
leads are assumed and we work in the wide band limit where @xmath133 .
then the total particle current through the dot , @xmath134 , is given by @xmath135\;.\label{equcurrent } \end{aligned}\ ] ] based on eq .
( [ equcurrent ] ) , we determine the electron current through the quantum dot numerically .
we note that in the present work we are only interested in the current response to finite voltage biases and temperature differences between the leads .
the generation of an electronic response merely by the interaction with a hot phonon bath would require different , energy - dependent densities of states in the leads , as was shown in ref .  .
in our calculations however , we consider identical lead densities of states . that is why @xmath136 for all @xmath4 , as long as @xmath137 and @xmath138 .    to find a relation between the polaronic and electronic functions , we decouple the fermionic and bosonic degrees of freedom in the electronic response function as is customary,@xcite @xmath139 as stated in sec .
[ secmodel ] the local oscillator is supposed to be strongly coupled to the heat bath .
accordingly we neglect the influence of the ep interaction on the dynamics of the phonon subsystem and evaluate the bosonic correlation function in ( [ equdefgelec ] ) using eq .
( [ equbosebath ] ) .
the first order electronic response functions then read @xmath140g_{dd}^{(1)\lessgtr}(\omega\pm s\omega_0;u ) \nonumber\\ & + n_b(s\omega_0)g_{dd}^{(1)\lessgtr}(\omega\mp s\omega_0;u )   \big ) \big\}\;.\nonumber\end{aligned}\ ] ] the corresponding electronic spectral function @xmath141 follows according to the steady - state equation ( [ equdefa ] ) as @xmath142 & \quad\mathrm{e}^{-\widetilde g^2\coth\theta}\big\ { i_{0}(\kappa)a^{(1)}(\omega;u ) + \sum_{s\ge 1 } i_{s}(\kappa)2\sinh(s\theta ) \nonumber \\ & \quad\times\big ( [ n_b(s\omega_0)+\bar f^{(1)}(\omega+s\omega_0;u)]a^{(1)}(\omega+ s\omega_0;u ) \nonumber \\ & \quad+ [ n_b(s\omega_0)+1-\bar f^{(1)}(\omega - s\omega_0;u)]a^{(1)}(\omega- s\omega_0;u)\big ) \big\}\;.\nonumber\end{aligned}\ ] ]      to determine the optimal variational parameter @xmath143 , we have to minimize the relevant thermodynamic potential .
this poses a problem , since for finite @xmath144 , the effective temperature determining the statistics of the dot electron is not known .
it will be given by the constitution of the steady state due to the coupling of the dot with the surroudings .
in the present paper , we suppose @xmath34 to be mainly determined by the ep interaction terms contained in eq .
( [ equhamiltonianstart ] ) , since the variational parameter was introduced to characterize the polaron - like quasilocalization of the dot electron .
therefore we assume the equilibrium thermodynamic potential of the system , before the temperature- and voltage differences were turned on , to be a reasonable first approximation for the variational function . the use of the well - known `` @xmath145 trick '' to take the interaction terms from @xmath56 into account results in@xcite @xmath146 the function @xmath147 represents the first order equilibrium spectral function if the interaction coefficients in eq .
( [ equhamiltonian ] ) are multiplied by the factor @xmath145 .
if we write @xmath147 in terms of @xmath148 , we can carry out the @xmath145-integration and find @xmath149 via @xmath150 and @xmath151 , the thermodynamic potential @xmath152 is a function of @xmath34 . for given parameters @xmath153 , @xmath154 , @xmath3 , @xmath6 and @xmath16 , the optimal variational parameter @xmath143 will be determined in equilibrium , i.e. , for @xmath155 and @xmath138 , by minimizing @xmath152 .
we then keep @xmath143 fixed and calculate the self energy , the spectral function and the resulting particle current for finite voltages or temperature differences between the leads .
depending on the type of molecular junction , the energies of the characteristic vibrational modes range from the order of @xmath156mev in small molecules@xcite to several @xmath157mev in @xmath158 molecules@xcite . in the following calculations ,
the corresponding model parameter @xmath3 will be used as the unit of energy , i.e. , we keep @xmath159 fixed and set @xmath160 , @xmath161 and @xmath162 .
we assume identical leads and work in the wide band approximation , so that the dot - lead coupling is determined by a single parameter : @xmath133 .
furthermore , we consider low equilibrium temperatures , @xmath163 . then , according to eq .
( [ equtildegamma ] ) , the renormalized dot - lead coupling parameter is approximately given as @xmath164 usually , the ratio of the bare dot - lead coupling parameter @xmath154 to the phonon energy @xmath3 is used to determine whether the system is in the adiabatic ( @xmath165 ) or anitadiabatic ( @xmath166 ) regime .
recently , eidelstein _
et al._@xcite argued that for strong ep interaction , i.e. , @xmath167 , the antiadiabatic regime can be extended to @xmath168 as long the exponential renormalization of the dot - lead coupling parameter in ( [ equtildegamma2 ] ) is so strong that @xmath169 . in this case , the `` mobility '' of passing charge carriers is reduced so far that the local oscillator is still fast enough to adjust to them individually . because then the physics is essentially those of the antiadiabatic regime , it can be called the `` extended antiadiabatic regime''.@xcite only when @xmath170 approaches @xmath3 the system crosses over to the adiabatic regime .    to investigate the latter situation , in the present work we consider comparable electronic and phononic timescales @xmath171 and moderate ep coupling @xmath172 .
since in our approach , the polaronic renormalization in eq .
( [ equtildegamma2 ] ) also depends on the value of the variational parameter @xmath31 $ ] , we can interpolate between the extended adiabatic regime and the aforementioned crossover regime .      , @xmath173 , and @xmath174 .
( a ) electronic equilibrium spectral function for @xmath175 , @xmath176 , @xmath177 , @xmath178 and @xmath179 .
the variational calculation ( black lines ) is compared to the small polaron limit @xmath180 ( red dashed lines ) and to the result for @xmath181 ( blue dot - dashed lines ) . in the latter two cases we set @xmath6 in such a way that @xmath150 is the same as in the variational calculation .
( b ) same as ( a ) , but for @xmath182 , @xmath183 , @xmath184 , @xmath185 and @xmath186 .
( c ) and ( d ) for the same parameters as in ( a ) and ( b ) , thermally induced current @xmath187 as a function of the temperature difference @xmath144 between the leads .
, title="fig : " ] + , @xmath173 , and @xmath174 .
( a ) electronic equilibrium spectral function for @xmath175 , @xmath176 , @xmath177 , @xmath178 and @xmath179 .
the variational calculation ( black lines ) is compared to the small polaron limit @xmath180 ( red dashed lines ) and to the result for @xmath181 ( blue dot - dashed lines ) . in the latter two cases we set @xmath6 in such a way that @xmath150 is the same as in the variational calculation .
( b ) same as ( a ) , but for @xmath182 , @xmath183 , @xmath184 , @xmath185 and @xmath186 .
( c ) and ( d ) for the same parameters as in ( a ) and ( b ) , thermally induced current @xmath187 as a function of the temperature difference @xmath144 between the leads .
, title="fig : " ]    to understand the thermoelectric response of the dot we first investigate the dots electronic spectral function in equilibrium .
thereby we keep the parameters @xmath159 , @xmath173 , @xmath138 , @xmath155 , and @xmath174 fixed and consider two parameter sets for @xmath154 and @xmath6 .
the first set is @xmath175 and @xmath176 , which means that the electronic and phononic subsystems react on a similar timescale and the dot level acts as a tunneling barrier between the leads . for these parameters
our variational calculation yields an optimal @xmath177 . as a consequence , the effective dot level given in eq .
( [ equdefeta ] ) is lowered to near resonance , @xmath178 .
moreover , the effective dot - lead coupling is reduced by about half , @xmath188 , and we are in the aforementioned crossover regime.@xcite    for the second parameter set we , reduced the dot - lead coupling to @xmath182 ( meaning the phononic subsystem is the faster one ) and raised the bare dot level to @xmath183 .
this results in @xmath184 , @xmath186 , and @xmath185 , so that the renormalized dot level is still far from resonance .
the electronic spectral functions for the first and second parameter set are presented as the black curves in figs .
[ fig2](a ) and [ fig2](b ) , respectively .
we compare them to two limiting cases : the red dashed curves represent the small polaron limit , where we set @xmath180 instead of the optimal @xmath143 .
then @xmath30 in eq .  ( [ equlft ] ) corresponds to the complete lang - firsov transformation . in this case
, we have @xmath189 and @xmath190 , respectively , and the system is described as being in the extended antiadiabatic regime .
the blue dot - dashed curves follow from setting @xmath181 , which represents a rigid quantum dot without ep interaction . in both limiting cases we set @xmath6 in such a way , that @xmath150 is the same as in the respective variational calculation .
in general , for @xmath181 , the electronic spectral function features a single broad band centered at @xmath191 . from setting @xmath192 in eqs .
( [ equgamma1 ] ) and ( [ equafrac ] ) , we see that it has a lorentzian shape and its width is given by the bare dot - lead coupling parameter @xmath154 . for finite ep interaction and @xmath180 ,
the polaronic character of the dot state is signalled by the appearance of several narrow side bands , which are given by the terms with @xmath193 in eq .
( [ equelectronic ] ) .
these side bands represent the emission and absorption of optical phonons by the incident electrons .
their width is determined mainly by the small parameter @xmath170 and their maxima are located multiple integers of @xmath3 away from the central ( zero - phonon ) peak .
note however , that the phonon peak corresponding to @xmath194 is suppressed in figs .
[ fig2](a ) and [ fig2](b ) , and only a small shoulder at @xmath195 remains .
this is a consequence of pauli blocking and can be understood from eqs .
( [ equsigmalessfourier1 ] ) and ( [ equgamma1a ] ) ; for small @xmath144 , we have @xmath196 and the fermi functions of both leads are nearly step - like . then the third and fourth lines in eqs .
( [ equsigmalessfourier1 ] ) and ( [ equgamma1a ] ) do not contribute in the region @xmath197 $ ] . from eq .
( [ equbarf ] ) we see that @xmath198 in this region and , consequently , the third and fourth line in eq .
( [ equelectronic ] ) do not contribute to the spectral function .
this `` floating '' of the phonon bands has been discussed for linear electric transport through molecular quantum dots.@xcite as we will see , it is also the reason why the low temperature thermoelectric response is determined by the shape and position of the zero - phonon peak alone .
the variational dot state , due to the moderate renormalization , features broad , overlapping side bands .
still , we find a considerable shift of spectral weight to higher energies and a suppression of the zero - phonon peak .
the resulting reduction of the low - energy tunneling rate is a main consequence of the polaron formation and is known as the franck - condon blockade.@xcite      in the following we calculate the particle current through the junction that results from a finite temperature difference @xmath144 between the leads .
let us again consider the two sets of parameters used in sec .
[ secstate ] .
keeping the respective variational parameters @xmath143 fixed we increase the temperature difference @xmath144 , whereby @xmath199 and @xmath25 .
this induces a net particle current through the junction , which is depicted as a function of @xmath144 in figs .
[ fig2](c ) and [ fig2](d ) . for the first parameter set , we find that the ep interaction enhances the thermally induced current for small temperature differences
. however , it reduces the current throughout the temperature range for the second parameter set .
this can be understood from eq .
( [ equcurrent ] ) : the net current through the dot depends on the electronic nonequilibrium spectral function and the difference of the fermi distribution functions of the left and right lead . for zero voltage bias and small temperature differences , the term @xmath200 in eq .
( [ equcurrent ] ) differs from zero only in the narrow region @xmath201 $ ] and changes sign at @xmath202 .
physically , this means that the temperature difference between the leads induces a flow of hot charge carriers from the left to the right lead , since @xmath200 is positive for @xmath203 .
this current , however , is compensated by a counterflow of cold carriers at @xmath204 , where @xmath200 is negative .
since @xmath205 in our calculations , dot states above the fermi energy have more spectral weight than the states below the fermi energy .
the net particle current is positive , i.e. , it goes from the left to the right lead .
the magnitude of the current is determined by the relative weight of the two flows , i.e. , by the slope of the dot spectral function around the fermi energy .
that is why the effect of the polaron formation on the thermally induced current crucially depends on the specific values of @xmath150 and @xmath170 .
for the first parameter set , the dot level @xmath178 lies near the fermi surface .
in this situation , the polaronic renormalization of the effective dot - lead coupling increases the slope of @xmath206 at @xmath202 , as can be seen in fig .
[ fig2](a ) .
consequently , in fig .
[ fig2](c ) the low temperature thermoelectric current grows with respect to the noninteracting case . for the second parameter set ( corresponding to the non - resonant situation ) the effect
is reversed : the slope of @xmath206 decreases and the thermoelectric response of the vibrating molecule is smaller than that of the rigid quantum dot .    for both parameter sets we find that , as @xmath207 , the maximum thermocurrent is largest in the case of zero ep interaction .
this can also be understood from eq .
( [ equcurrent ] ) . for large @xmath144 , the region @xmath201 $ ] grows and @xmath200 is finite far from the fermi surface . now the thermoelectric current is given not by the slope of @xmath206 at @xmath202 but by the ratio of the integrated spectral weight below and above the fermi surface . because @xmath4 is low , the phonon bands in @xmath208 are weighted according to a poisson distribution with the parameter @xmath209 . as @xmath34 grows , a considerable portion of the total spectral weight is shifted to the phonon bands at @xmath210 . even for large temperature differences ,
their contribution to the thermocurrent will be exponentially small , and the maximum current decreases with respect to the rigid quantum dot .
( b ) and ( d ) , i.e. , @xmath159 , @xmath173 , @xmath174 , @xmath182 , @xmath183 , @xmath184 , @xmath185 and @xmath186 .
( a ) electronic spectral functions for fixed @xmath211 and @xmath138 .
we compare the variational calculation ( black line ) to the small polaron limit @xmath180 ( red dashed line ) and to the result for @xmath181 ( blue dot - dashed line ) .
( b ) current @xmath187 as a function of the voltage difference @xmath212 between the leads for fixed @xmath211 .
( c ) thermovoltage @xmath213 as a function of @xmath144.,title="fig : " ] + ( b ) and
( d ) , i.e. , @xmath159 , @xmath173 , @xmath174 , @xmath182 , @xmath183 , @xmath184 , @xmath185 and @xmath186 . ( a ) electronic spectral functions for fixed @xmath211 and @xmath138 .
we compare the variational calculation ( black line ) to the small polaron limit @xmath180 ( red dashed line ) and to the result for @xmath181 ( blue dot - dashed line ) .
( b ) current @xmath187 as a function of the voltage difference @xmath212 between the leads for fixed @xmath211 .
( c ) thermovoltage @xmath213 as a function of @xmath144.,title="fig : " ] + ( b ) and ( d ) , i.e. , @xmath159 , @xmath173 , @xmath174 , @xmath182 , @xmath183 , @xmath184 , @xmath185 and @xmath186 . ( a ) electronic spectral functions for fixed @xmath211 and @xmath138 .
we compare the variational calculation ( black line ) to the small polaron limit @xmath180 ( red dashed line ) and to the result for @xmath181 ( blue dot - dashed line ) .
( b ) current @xmath187 as a function of the voltage difference @xmath212 between the leads for fixed @xmath211 .
( c ) thermovoltage @xmath213 as a function of @xmath144.,title="fig : " ]    we have seen how , in general , a temperature difference @xmath144 between the leads will induce a particle current through the dot . in a typical experiment ,
this thermoelectric effect is measured by applying a voltage difference @xmath212 in such a way that the thermally induced current is compensated : @xmath214 for finite @xmath144 and @xmath212 .
the determination of the so - called thermovoltage , @xmath213 , and its dependence on the ep coupling will be the subject of the following numerical calculations .
thereby , we use the second parameter set from sec .
[ secstate ] , i.e. @xmath159 , @xmath173 , @xmath174 , @xmath183 and @xmath182 , resulting in @xmath184 , @xmath186 and @xmath185 .    to understand the mechanism , we first consider a finite but fixed temperature difference @xmath211 .
for @xmath138 , the corresponding electronic spectral function is shown as the black curve in fig .
[ fig3](a ) , where it is again compared to the noninteracting case and the small polaron limit .
we calculate the total current @xmath187 as a function of the voltage bias @xmath212 , whereby @xmath215 and @xmath216 .
the results are depicted in fig .
[ fig3](b ) .    for negative voltages ,
the current signal is in qualitative agreement with the results of previous works.@xcite in a nutshell , the formation of a polaron like state at the dot reduces the zero - bias conductance due to the franck - condon blockade , but for growing voltage we find steps in the current signal whenever @xmath217 , with @xmath218 $ ] . here
the chemical potential of the right lead crosses the phonon side bands in @xmath208 and resonant transport of electrons takes place via the emission and subsequent absorption of an equal number of phonons .
the right lead is the cold one and has a step - like fermi function . because of this
, the width of the current steps is mainly determined by the width of the phonon bands in the spectral function , i.e. , by @xmath170 . in the variational calculation
the bands overlap considerably ( see fig .
[ fig3](a ) ) and the resonant steps in the current - voltage signal are smeared out .    for positive voltage bias
we find no current steps even in the small polaron limit .
now it is the chemical potential of the hot left lead that crosses the renormalized dot level and the phonon bands .
the width of the resonant tunneling steps in @xmath219 does not depend on @xmath170 alone but also on the width of the soft fermi - surface .
since @xmath211 is of the order of @xmath3 , the current steps are smeared out and therefore no longer discernible .    based on the current - voltage signal
we now determine the thermovoltage @xmath213 as a function of @xmath144 numerically and show the result in fig .
[ fig3](c ) .
we find that for small temperature differences , the absolute value of @xmath213 grows strongest for @xmath180 .
this is a consequence of the reduced electrical conductance of the quantum dot in the small polaron picture ; for a given @xmath144 , relatively large voltages are necessary to compensate the corresponding thermocurrent .
for large temperature differences , we find a maximum value @xmath220 for all three data sets given in fig .
[ fig3](c ) .
this voltage corresponds to the postition of the zero - phonon step in @xmath219 in fig .
[ fig3](b ) . here , the systems conductance grows considerably and any thermally induced current can easily be compensated by a slight growth in @xmath212 .      , @xmath173 , @xmath174 , @xmath182 , as in fig .
( a ) thermovoltage @xmath213 as a function of the bare dot level @xmath6 for @xmath221 .
we compare the variational calculation ( black line ) to the small polaron limit @xmath180 ( red dashed line ) and the case with zero ep coupling ( blue dash - dot line ) .
( b ) optimal variational parameter @xmath143 as a function of the bare dot level .
( c ) same as ( a ) but with @xmath211 .
( d ) electronic spectral functions for @xmath222 , @xmath223 and several temperature differences.,title="fig : " ] + , @xmath173 , @xmath174 , @xmath182 , as in fig .  [ fig3 ] . ( a ) thermovoltage @xmath213 as a function of the bare dot level @xmath6 for @xmath221 .
we compare the variational calculation ( black line ) to the small polaron limit @xmath180 ( red dashed line ) and the case with zero ep coupling ( blue dash - dot line ) .
( b ) optimal variational parameter @xmath143 as a function of the bare dot level .
( c ) same as ( a ) but with @xmath211 .
( d ) electronic spectral functions for @xmath222 , @xmath223 and several temperature differences.,title="fig : " ] + , @xmath173 , @xmath174 , @xmath182 , as in fig .  [ fig3 ] . ( a )
thermovoltage @xmath213 as a function of the bare dot level @xmath6 for @xmath221 .
we compare the variational calculation ( black line ) to the small polaron limit @xmath180 ( red dashed line ) and the case with zero ep coupling ( blue dash - dot line ) .
( b ) optimal variational parameter @xmath143 as a function of the bare dot level .
( c ) same as ( a ) but with @xmath211 .
( d ) electronic spectral functions for @xmath222 , @xmath223 and several temperature differences.,title="fig : " ] + , @xmath173 , @xmath174 , @xmath182 , as in fig .  [ fig3 ] . ( a ) thermovoltage @xmath213 as a function of the bare dot level @xmath6 for @xmath221 .
we compare the variational calculation ( black line ) to the small polaron limit @xmath180 ( red dashed line ) and the case with zero ep coupling ( blue dash - dot line ) .
( b ) optimal variational parameter @xmath143 as a function of the bare dot level .
( c ) same as ( a ) but with @xmath211 .
( d ) electronic spectral functions for @xmath222 , @xmath223 and several temperature differences.,title="fig : " ]    in the previous sections , we chose the bare dot level @xmath6 in such a way , that the renormalized level @xmath150 was equal in the scenarios with and without ep interaction . in this way , we concentrated on the influence of the renormalized dot - lead coupling on the thermovoltage .
however , in an experimental situation , the dot level may be manipulated , e.g. , by applying a gate voltage through a third electrode , in order to optimize the thermoelectric effect .
to investigate this situation , we now consider @xmath6 as our free parameter .
for each @xmath6 we temporarily set @xmath138 and @xmath155 to determine the optimal parameter @xmath143 .
we then keep @xmath143 fixed , set @xmath144 to a finite value and calculate the thermovoltage @xmath213 .
the remaining parameters are the same as in sec .
[ seccompvol ] , i.e. , @xmath159 , @xmath173 , @xmath174 , and @xmath182 .
first , we consider a small temperature difference @xmath221 .
our result for the thermovoltage as a function of the bare dot level is presented in fig .
[ fig4](a ) . in general
, @xmath224 features two resonances of opposite sign and goes to zero when the renormalized dot level crosses the fermi surface , i.e. , when @xmath225 . at this point
the spectral function is symmetrical around @xmath202 , and the thermoelectric flow and counterflow between the leads cancel exactly . as @xmath226 grows , the net thermocurrent and , consequently , the thermovoltage increases . in accordance with our results in sec .
[ seccompvol ] , the maximum thermovoltage is largest in the small polaron picture due to the strong renormalization of the dot - lead coupling . when @xmath226 grows further and the zero phonon peak shifts away from the fermi surface , the dot density of states near @xmath202 decreases , and the thermoelectric effect vanishes again .    in fig .
[ fig4](a ) , the small polaron result is shifted from the @xmath181 curve by the value of the polaron binding energy @xmath153 , which can be understood from setting @xmath180 in eq .
( [ equdefeta ] ) . in both cases ,
the thermovoltage signal runs linearly through zero .
we note that for growing @xmath16 ( not shown here ) its slope reduces , while the position of the positive ( negative ) resonance shifts to lower ( higher ) @xmath6 . in this
regard , our result resembles the sawtooth - like thermopower signal that was predicted in ref .   and experimentally measured in ref .  .
there , the periodicity of the thermopower oscillations was determined by the difference in the ground - state energies for different numbers of electrons on the dot . in our model
however , we only account for a single dot electron .
that is why we only observe a single `` tooth '' , i.e. , only two thermovoltage resonances in fig .
[ fig4](a ) .
moreover , we find no side peaks in the thermoelectric signal that could be attributed to the phonon side bands in @xmath208 .
this is a consequence of the floating effect discussed in sec .
[ secstate ] ; since the thermovoltages in fig
.  [ fig4](a ) are small , the expression @xmath227 in eq .
( [ equcurrent ] ) differs from zero only within the region @xmath228 $ ] . in this region ,
the phonon side bands in @xmath206 are suppressed and do not contribute to the thermoelectric transport .
that is why the resonance signal in fig .
[ fig4](a ) is given by the varying position of the zero - phonon peak only .
consequently , in fig .
[ fig4](a ) the variational calculation with its moderate renormalization of the dot - lead coupling predicts a weaker resonance signal than the small polaron picture , which is also shifted from the result for the rigid dot by less than @xmath153 .
we find that in contrast to the two limiting cases , the strength of the resonance now depends on the sign of @xmath213 .
this can be understood from fig .
[ fig4](b ) , where we show the optimal variational parameter as a function of the dot level . for @xmath205
, we have @xmath229 like in the previous sections . when the renormalized dot level crosses the fermi surface , @xmath143 approaches @xmath230 and a small transient polaron forms at the dot . then the thermoelectric effect increases .
figure  [ fig4](c ) shows the thermovoltage signals after having raised the temperature difference to @xmath211 .
as expected , for all three cases , the heights and widths of the resonances grow with respect to fig .
[ fig4](a ) .
more importantly , for finite ep coupling the curves now feature side bands at a distance of @xmath3 from their maximum resonances .
this becomes possible because for large temperature differences the floating condition for the phonon bands is relaxed : since @xmath231 , the fermi surface of the hot left lead softens considerably .
now the third and fourth lines in eq .
( [ equelectronic ] ) contribute to the spectral function even for @xmath228 $ ] .
for example , fig .
[ fig4](d ) shows the electronic spectral function in the variational calculation for @xmath222 and several temperature differences . since @xmath232
, we have @xmath233 . for growing @xmath144 , a small
peak appears near the fermi - surface that is related to the polaronic state with one phonon at @xmath234 .
when , for varying @xmath6 , this peak crosses the fermi surface , the side peak in the thermovoltage signal in fig .
[ fig4](c ) appears .
note , however , that the floating condition still holds for the cold right lead , as can be seen from the shoulder appearing near @xmath235 in fig .
[ fig4](d ) .
as the thermovoltage grows , the window @xmath228 $ ] closes and the shoulder shifts towards @xmath202 .
in this work , we investigated the steady - state thermoelectric transport through a vibrating molecular quantum dot in the crossover regime far from the antiadiabatic limit . within a kadanoff - baym formalism that is generalized to account for different lead temperatures ,
the nonequilibrium dot self - energy was calculated to second order in the dot - lead interaction coefficient . in order to account for the polaronic character of the dot state
, we applied a variational small - polaron transformation and determined the degree of transformation by minimizing the relevant thermodynamic potential .
in essence , we calculated the current induced by a finite temperature difference between the leads . for small temperature differences , the influence of the electron - phonon ( ep ) interaction strongly depends on the specific system parameters ; for a near - resonant dot level , the interacting quantum dot acts as a more efficient energy filter and the thermocurrent increases . in the tunneling regime , however , we found a reduction of the thermocurrent due to the decreasing density of states near the fermi surface .    in order to relate our results to experiment , we determined the thermovoltage required to compensate the thermally induced particle current at a given temperature difference .
we found that in principle the franck - condon blockade boosts the thermovoltage through the reduction of the systems electrical conductance .
for intermediate ep coupling , the small polaron picture overestimates this effect .
finally , we determined the thermovoltage as a function of the dots energy level . because of pauli blocking , we found no phonon features in the low - temperature thermoelectric signal .
however , our variational calculation predicts an asymmetrical line shape due to the formation of a small polaron as the dot level drops beneath the equilibrium fermi energy . for large temperature differences of the order of the phonon energy , pauli blocking is relaxed and
the thermovoltage signal features multiple resonances that can be attributed to resonant transport through vibrational dot states .
the present study should be considered a first step in applying our variational ansatz to the thermoelectric transport through molecular junctions .
although it captures the essential physics , it must be extended in several directions .
most importantly , we have yet to consider energy transport through the junction
. the efficiency of energy deposition by the excitation of local phonons will strongly depend on the effective ep interaction , i.e. , on the polaronic character of the dot state . to investigate the subsequent heating or cooling of the molecule
, the effective temperature of the dot electrons has to be determined numerically , for which several methods have been proposed.@xcite this should have some influence on the variational parameter @xmath34 as well . in our previous works ,
a voltage - dependent variational parameter was responsible for the junctions negative differential conductance .
therefore , it would be highly desirable to carry this nonlinear behavior over to the discussion of thermoelectric transport .    another worthwhile extension of our work concerns coulomb interaction effects .
recently , andergassen _
et al._@xcite argued for the negative@xmath236 anderson model ( which neglects the coupling to the phonon degrees of freedom however ) that the resulting charge kondo effect leads to a large enhancement of the linear response thermopower due to the highly asymmetric dot spectral function .
this effect was shown to be tunable applying a gate voltage .
it would be interesting to reexamine this problem for a model with additional ep interaction , using our variational scheme , since  as we have demonstrated in sec .
[ seclevel]the effective ep interaction ( described by the variational parameter ) is strongly influenced by the gate voltage . in this connection
we like to stress that for a combined holstein - hubbard quantum dot model it has been shown that strong ep coupling may result in a net attractive coulomb interaction.@xcite then , depending on the energy of the dot level , a variational ansatz might be able to interpolate between the positive@xmath236 holstein - hubbard dot model and the effective anisotropic kondo model regime investigated in ref .  .
46ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]  + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * , ( ) @noop * * ,   ( ) @noop * * , ( ) @noop * * ,   ( ) @noop * * , ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop   @noop * * , ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) ,   @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop _ _  ( , ,  ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( )