recent developments in the fabrication and control of nanostructures and molecular devices have stimulated a growing interest in the study and research of heat conductance in such devices .
@xcite as efforts are being made to better utilize and control these devices , the understandings of the mechanisms for accumulation and dissipation of heat are of importance .
the literature on the topic considers different physical setups , differing in the process by which heat is mediated , the devices under investigation and the nature of the heat baths involved in the heat transfer process .  @xcite    we consider heat transfer between two baths held at different temperatures through a steady - state current of energy .
the linear response regime , which relies on the equilibrium properties of the system , arises in the limit of small gradient of temperatures .
a linear response study of a system similar to the one that is presented in this paper was done recently by entin - wohlman and others and the thermopower properties were calculated  @xcite . in the general case , however , the system is far from equilibrium , as the temperature gradient is finite .
the description of such systems is a challenging problem of great interest in current research , as many of the concepts and techniques used to describe equilibrium setups are inadequate .
exact solutions of systems far from equilibrium are particularly desirable as they may offer both a benchmark and an unbiased understanding of the underlying physics .
in this paper we present an asymptotically - exact calculation of the heat current through a molecular junction , under explicit nonequilibrium conditions manifested by a finite temperature gradient between the two baths to which the molecule is coupled  the electrons in a conduction band and a bosonic bath . in a typical molecular bridge ,
molecular orbitals are coupled simultaneously to the lead electrons and to the vibrational modes of the molecule , with the former degrees of freedom reduced to a single effective band in the absence of a bias voltage .
@xcite this coupling to the vibrational modes of the molecule is believed to have an essential role in heat transfer processes .
@xcite in the continuum limit , a minimal model for an unbiased molecular bridge therefore consists of a vibrational mode that is coupled by displacement to the conduction electrons at the origin , and is also coupled to a bosonic bath of vibrational modes , as described by the hamiltonian of eq .
( [ h ] ) .
the solution is asymptotically exact in the sense that it is governed by the smallness of the electron - phonon coupling @xmath0 with respect to the effective energy band of the conductance electrons .
the thermal properties of the bosonic hamiltonian of eq .
( [ h - bosonic ] ) , and similar hamiltonians describing heat current between harmonic baths , were studied before in different physical setups than presented here  @xcite .
the results we derive in this paper corroborate those studies , and adjust them to describe the electronic systems which we shall present .
the paper is organized as follows .
we begin in sec .
[ sec : model ] by presenting the model and the physical systems which it may represent , under suitable mappings . in sec .
[ sec : mapping ] we then introduce the nonequilibrium condition and map the model onto a form quadratic in bosons . this quadratic nature of the hamiltonian is then exploited , in sec .
[ sec : resonance_heat_current ] , in order to calculate exactly the heat current in the system . in sec .
[ sec : master_equation ] we then turn to a master equation approach , which is perturbative in nature , in order to calculate again the heat current .
this approach allows us to gain some useful physical understanding of the processes involved .
section  [ sec : off - resonance ] addresses the case where the localized fermionic level is taken off - resonance , which breaks the particle - hole symmetry of the problem and adds a linear term to our quadratic hamiltonian .
we then turn , in sec .
[ sec : conclusions ] to present our conclusions .
the model we consider consists of a single molecule coupled simultaneously to a conduction band and to a bosonic bath , which may represent , depending on the context , either the substrate phonons or any other continuum of bosonic degrees of freedom .
a finite temperature gradient is applied between the electronic and bosonic baths , generating a nonzero heat flow from the hotter bath to the colder one .
physically this model can describe any of a number of systems , three of which are portrayed schematically in fig .
[ fig : fig1 ] : a single molecule adsorbed on a metallic surface ( upper panel ) , a molecular bridge ( middle panel ) , and a single - molecule transistor embedded in one of the arms of an aharonov - bohm interferometer ( lower panel ) . under suitable conditions , to be specified below ,
each of these setups can be described by the generic continuum - limit hamiltonian  @xcite @xmath1 which is the focus of the present work . in the above hamiltonbian
, the one - dimensional fermionic field @xmath2 represents the conduction - electron degrees of freedom , @xmath3 denotes the molecular vibrational mode , and @xmath4 are the modes of the bosonic bath .
the fermionic field obeys canonical anticommutation relations @xmath5 subject to the regularization @xmath6 , where @xmath7 is a suitable short - distance cutoff . in the above model
we have restricted ourselves to the case where the coulomb repulsion between the lead and the localized fermionic level is negligible , and omitted that term .
the bosonic bath is characterized by the coupling function @xmath8 which is assumed to have the standard power - law form @xmath9 here , @xmath10 is a high - energy cutoff and @xmath11 is a dimensionless coupling constant parameterizing the coupling strength to the bosonic bath .
the power @xmath12 is of particular interest as it corresponds to an ohmic bath.@xcite the parameter @xmath0 describes the displacement coupling between the electrons and the vibrational mode , while @xmath13 represents ordinary potential scattering .    ) .
top to bottom : ( a ) a molecule adsorbed on a surface , with the bulk phonons serving as a bosonic thermal bath ; ( b ) a single - molecule transistor , where the phonons of the substrate are coupled to local vibrational mode ; ( c ) an aharonov - bohm interferometer with a molecular device embedded in one of its arms.,title="fig:",scaledwidth=44.0% ] ) .
top to bottom : ( a ) a molecule adsorbed on a surface , with the bulk phonons serving as a bosonic thermal bath ; ( b ) a single - molecule transistor , where the phonons of the substrate are coupled to local vibrational mode ; ( c ) an aharonov - bohm interferometer with a molecular device embedded in one of its arms.,title="fig:",scaledwidth=44.0% ] ) .
top to bottom : ( a ) a molecule adsorbed on a surface , with the bulk phonons serving as a bosonic thermal bath ; ( b ) a single - molecule transistor , where the phonons of the substrate are coupled to local vibrational mode ; ( c ) an aharonov - bohm interferometer with a molecular device embedded in one of its arms.,title="fig:",scaledwidth=44.0% ]    the generic model described by the hamiltonian of eq .
( [ h ] ) is found in discussions of different physical systems consisting of a single molecule which is coupled to large environments .
it has a long history that dates back to the 1970s , when it was proposed as a model for the electron - phonon coupling in mixed - valence compounds .
@xcite in the modern context of nanostructures it is expected to properly describe the physics of single - molecule devices away from coulomb - blockade valleys where a single unpaired spin resides on the molecule .
we shall now turn to present several of these setups .
the most direct realization of the model given by eq .
( [ h ] ) is that of a single molecule adsorbed on a metallic surface .
adsorbed molecules on surfaces have been intensely studied over the years , both theoretically and experimentally ( for a brief review see ref .
[ ] ) .    in such molecules , the vibrational modes of adsorbed molecules
relax by interaction with the surface conductance electrons as well as by coupling to the bulk phonons , and the interaction between the electrons and the vibrational mode plays an important role .
@xcite one of the common models used to describe such interactions  @xcite is the hamiltonian of eq .
( [ h ] ) , prior to the reduction of the relevant conductance electrons degrees of freedom to one - dimensional fields . adopting that description , the bosonic bath described by the operators @xmath14 and @xmath4 in our hamiltonian represent the bulk phonons .
the molecule is then brought into contact with an stm - tip or other electronic lead , and the conductance electrons couple to the electronic level on the molecule .
focusing on the electronic mode that couples to the molecule , we can map them onto the one - dimensional field that is represented by the electronic field operators @xmath15 and @xmath2 in the hamiltonian .
a detail of the process can be found in ref .  [ ] .
in such a system , generally one would take @xmath16 , as the conductions electrons are held in resonance with the electronic level on the molecule .      another system that can be described by the model of eq .
( [ h ] ) is that of a molecular bridge  a single molecule sandwiched between two electronic leads . in such single - molecule devices the electron - phonon interaction plays an important role , as the molecular orbitals are coupled simultaneously to the lead electrons and to the vibrational modes of the molecule itself .
the molecular bridge is typically placed on an insulating substrate which provides an additional phononic bath that couples to the molecular vibrational modes .
@xcite    focusing on a molecular junction held between leads with no bias voltage , one can choose a symmetric and anti - symmetric basis for the lead electrons , thus the molecule interacts with an effective single electronic band .
the molecule itself is modeled by a single spinless electronic level @xmath17 with energy @xmath13 which is coupled by displacement to a local vibrational mode @xmath3 with frequency @xmath18 . this vibrational mode
is further coupled by displacement to a bath of phonons @xmath4 . the level
is then coupled to the single band of electrons via a hopping matrix element @xmath19 .
the resulting hamiltonian is given by @xmath20 with @xmath21 and @xmath22 typically the hamiltonian is treated either in the weak - coupling limit using perturbation theory in @xmath0 , or using the lang - firsov transformation  @xcite and the polaronic approximation in the limit where @xmath19 is small . in the limit of a broad level close to resonance , however , this hamiltonian can be reduced  @xcite to the continuum limit hamiltonian of eq .
( [ h ] ) as described in detail in ref .
this mapping process give rise to an effective high - energy cutoff which is determined by the hybridization width of the level @xmath23 .
the effective bandwidth @xmath24 of the continuous field @xmath15 is related to the hybridization width by @xmath25 .
the final system considered is an aharonov - bohm interferometer with a molecular device embedded in one of its arms .
this setup is a more complicated variant of the molecular bridge described above , as there is an additional transmission channel between the two leads , and a magnetic flux @xmath26 threading the ring . lately
, this system has attracted a considerable interest , following the work of entin - wohlman and others  @xcite who computed the transport coefficients relating the linear - response electric and heat currents to the voltage bias and temperature gradient .    restricting our attention to the case where the two electronic leads are held in identical temperature and chemical potential , the hamiltonian describing such a device
is given  @xcite by @xmath27 , with @xmath28 here , @xmath29 describes the left ( l ) and right ( r ) lead electrons , @xmath30 describes the arm of the interferometer without the device , and @xmath31 describes the device itself and its coupling to the leads .
the general term @xmath32 may contain any local interactions at the dot that do not involve the lead electrons . in our case
, it will include the local vibrational mode , the electron - phonon coupling and the bosonic bath .
however , for the purpose of the mapping shown here , no further assumptions are needed on the form of @xmath32 beyond it not involving the lead electrons .
we begin by transferring to the symmetric and anti - symmetric basis , defining the new symmetric field @xmath33\ ] ] where @xmath34 .
the green function @xmath35 pertaining to this field at @xmath36 can now be calculated from the green functions of the left and right fields at the origin , given by matrix form as @xmath37^{-1 } & -te^{i\varphi } \\                                                  -te^{-i\varphi } &                                          \left[g_r(z)\right]^{-1}\end{array }          \right)^{-1}.\end{aligned}\ ] ] using this we write the dot s green function as @xmath38^{-1},\ ] ] where @xmath39 is the self - energy contribution of the @xmath32 .
close to resonance , the role of the symmetric electronic green function @xmath40 will be to renormalize @xmath13 and the level width .
therefore , our effective hamiltonian of eq .
( [ h ] ) generates the same correlation functions pertaining to the dot operators as the original hamiltonian , ensuring the validity of the calculations of the heat current . to this end
, we have to adjust the parameters of our effective hamiltonian such that the effective electronic band - width is @xmath41 with @xmath42 , where @xmath43 the renormalized level energy is effected by the flux threading the ring @xmath44 where @xmath45 .
as the flux @xmath26 contributes only to the effective @xmath13 , which will be shown to effect the heat current only at quartic orders or through the renormalization of the coupling coefficient @xmath0 , one concludes that the magnetic flux effects the heat current similarly .
moreover , its effect is an even function of the flux .
the nonequilibrium conditions that will give rise to the heat current will be manifested by assuming to hold the electronic and bosonic baths at different temperatures @xmath46 and @xmath47 , respectively .
this temperature gradient may be intentional and well controlled , or can be the by - product of some other dynamics that inevitably causes an imbalance between the two heat baths . in either case
we assume that all other energy relaxation mechanisms between the electronic and the bosonic bath are sufficiently inefficient such that all local relaxation can be regarded as mediated by the molecule . in the keldysh
spirit,@xcite we account for the temperature difference by taking the initial density operator in the distant past to have the form @xmath48 where @xmath49 and @xmath50 are the reciprocal temperatures and @xmath51 represent the two decoupled baths .
the system is then evolved in time according to the full hamiltonian @xmath52 until steady state is reached .
the electronic hamiltonian @xmath53 contains a natural high - energy cutoff or bandwidth @xmath54 , which , depending on the context , may represent either the actual conduction - electron bandwidth or the hybridization width of a certain molecular orbital .
our subsequent solution of the nonequilibrium state is confined to the weak - coupling regime , @xmath55 , which serves as a prerequisite for some of the realizations of hamiltonian of eq .
( [ h ] ) depicted in fig .  [ fig : fig1 ] . in this limit
one can apply abelian bosonization  @xcite to convert the hamiltonian and the initial density operator to a form quadratic in bosonic operators . specifically , the hamiltonian of eq .
( [ h ] ) takes the form @xmath56
\sum_{k > 0 } \xi_k ( a^{\dagger}_k + a_k ) \nonumber \\       & + ( b^{\dagger } + b )          \sum_n \lambda_n                 ( \gamma^{\dagger}_n + \gamma_n ) , \label{h - bosonic}\end{aligned}\ ] ] where the first two terms on the right - hand side correspond to the free hamiltonian terms @xmath53 and @xmath57 that appear in @xmath58 ; @xmath59 and @xmath60 with @xmath61 ( @xmath62 being the size of the system ) are canonical bosonic creation and annihilation operators corresponding to the fourier modes of the fermionic density ; @xmath63 equals @xmath64 ; and the coefficients @xmath65 , which have the dimension of one over length , are given by @xmath66 note that we have omitted in eq .
( [ h - bosonic ] ) the contribution of the @xmath67 mode of the fermionic density as it has no effect on our problem of interest . as for the coupling constants
@xmath68 and @xmath69 , these have the dimension of energy times length , and are given to linear order in @xmath0 and @xmath13 by  @xcite @xmath70 as we shall argue below , one must include higher orders in @xmath13 to account for its effect on the heat current .
our interest is in the steady - state heat current flowing between the two baths .
formally there are several heat - current operators one can define , e.g. , the heat current flowing into the bosonic bath or the heat current flowing out of the electronic bath , all of which must coincide in steady state . for convenience
we shall focus on the heat current flowing into the bosonic bath , whose corresponding operator @xmath71\ ] ] reads @xmath72 our goal is to evaluate the steady - state expectation value @xmath73 .
first we consider the case where @xmath74 , which corresponds for a molecular bridge to electronic resonance conditions . technically , this limit is somewhat easier to address as the bosonic hamiltonian is purely quadratic .
a nonzero @xmath69 introduces a term linear in bosonic operators to the hamiltonian of eq .
( [ h - bosonic ] ) , whose treatment requires some care .
we defer discussion of the off - resonance case to sec .
[ sec : off - resonance ] below .      to compute the heat current we begin by writing it in the form @xmath75 where @xmath76 in steady state all two - time correlation functions reduce to a function of the time difference only .
this allows one to convert to the energy domain by fourier transforming with respect to the time difference .
we apply this procedure to @xmath77 , and aiming at evaluating @xmath78 we resort to ordinary diagrammatic perturbation expansion to write an expression for @xmath79 . to this end , it is useful to define the free ( i.e. , @xmath80 ) green functions of the bath bosons in matrix form as @xmath81,\ ] ] for the retarded and advanced functions , and @xmath82,\ ] ] for the lesser and greater functions , where @xmath83 is the bose - einstein distribution function corresponding to the temperature @xmath84 .
similarly , we shall denote the fully dressed green functions of the molecular vibrational mode @xmath3 by @xmath85 and @xmath86 , for which the @xmath87 matrix notation @xmath88\ ] ] is used .
[ the same @xmath87 matrix notation applies to @xmath89 above . ]
having laid out the building blocks for the perturbation expansion , we rely on langreth theorem  @xcite to have the identity , at steady - state    @xmath90_{1 p }                    \left\ {                             \left [                                    g^a_{\gamma , n}(\epsilon )                             \right]_{22 }                           - \left [                                    g^a_{\gamma , n}(\epsilon )                             \right]_{11 }                    \right\ }                  + \left [                            g^{r}(\epsilon )                    \right]_{1 p }                    \left\ {                             \left [                                    g^>_{\gamma , n}(\epsilon )                             \right]_{22 }                           - \left [                                    g^>_{\gamma , n}(\epsilon )                             \right]_{11 }                    \right\ }           \right ] , \label{g^>_b , n}\ ] ]    which we shall now turn to evaluate .    focusing initially on the retarded and advanced green functions @xmath85 , these acquire the form @xmath91^{-1 } , \label{g^ra}\ ] ] where @xmath92 \nonumber \\     & + \sum_n \lambda_n^2              \left [                      \frac{1 }                           { \epsilon - \omega_n \pm i\eta }                    - \frac{1 }                           { \epsilon + \omega_n \pm i\eta }              \right ] , \label{s^ra}\end{aligned}\ ] ] are the corresponding self - energy functions .
there are two contributions to @xmath93 : one due to the coupling to the electronic bath [ the first term on the right - hand side of eq .
( [ s^ra ] ) ] , and another due to the coupling to the bosonic bath [ the second term on the right - hand side of eq .
( [ s^ra ] ) ] . denoting these two terms by @xmath94 and @xmath95 , respectively ,
the former can be expressed in a closed analytical form  @xcite in terms of the exponential integral function :  @xcite @xmath96 . \label{s^ra_e}\end{aligned}\ ] ] here , @xmath97 is the density of states per unit length , and @xmath98 equals @xmath99 . as for the second contribution @xmath95 , it can be conveniently expressed in terms of the coupling function @xmath100 , @xmath101 d\epsilon ' . \label{s^ra_b}\ ] ] for an ohmic bath with @xmath12 , eq .
( [ s^ra_b ] ) can be evaluated in closed analytic form to obtain @xmath102 .
\label{s^ra_b - ohmic}\ ] ]    proceeding to the lesser and greater green functions @xmath86 , these read @xmath103          \sigma^ { < , > } ( \epsilon )          g^{a}(\epsilon ) , \label{g^<>}\ ] ] with the lesser and greater self - energy functions @xmath104 .
\label{s^<>}\end{aligned}\ ] ] here , @xmath105 is the bose - einstein distribution function corresponding to the temperature @xmath106 . note that eq .
( [ s^ < > ] ) can be conveniently related to the two components of the retarded self - energy through @xmath107 which generalizes the standard equilibrium relation @xmath108 .    to evaluate @xmath109 of eq .
( [ g^>_b , n ] ) , it is useful to utilize the identities @xmath110_{1 p } & =       |g(\epsilon)|^2 2 \omega_0 ( \epsilon + \omega_0 )                         \sigma^{>}(\epsilon ) , \nonumber \\
\sum_{p = 1 , 2 }       \left [               g^{r}(\epsilon )       \right]_{1 p } & =       |g(\epsilon)|^2 ( \epsilon + \omega_0 )       \left [               \epsilon^2 - \omega_0 ^ 2               - 2\omega_0 \sigma^{a}(\epsilon )       \right ] , \end{aligned}\ ] ] which follow directly from eqs .
( [ g^ra ] ) and ( [ g^ < > ] ) .
here we have introduced the auxiliary function @xmath111 and made use of the fact that @xmath112^{\ast}$ ] . inserting these identities into eq .
( [ g^>_b , n ] ) one obtains    @xmath113             \left [                     \delta(\epsilon - \omega_n )                     + \delta(\epsilon + \omega_n )             \right ] , \ ] ]    which reduces by virtue of eq .  ( [ s^<>-via - s^r ] ) to @xmath114             \left [                     \delta(\epsilon - \omega_n )                     + \delta(\epsilon + \omega_n )             \right ] .\ ] ] plugging this result into eq .  ( [ j_q - via - g ] ) , employing the relations @xmath115 and @xmath116 = - [ n_e(\epsilon ) - n_b ( \epsilon ) ] $ ] , and exploiting the fact that @xmath117 is restricted to positive energies , we finally arrive at @xmath118 , \label{j_q - general}\ ] ]    where we have explicitly written out the function @xmath119 that appears in the integrand [ see eq .  ( [ g - def ] ) ] .
equation  ( [ j_q - general ] ) is the central result of this paper .
it provides an exact expression for the heat current corresponding to the bosonic hamiltonian of eq .
( [ h - bosonic ] ) , for a general coupling function @xmath117 . since @xmath120 $ ] with @xmath121 is positive definite for @xmath122 ( negative definite for @xmath123 ) , the heat current flows , as it physically should , from the hotter bath to the colder one .
below we analyze in detail the characteristics of @xmath124 in different temperature and coupling regimes .
the form of expression for the heat current bears a similarity to landauer formula .
other works studying thermal conductance through local contacts have arrived at a similar expressions or used a landauer - type expression as a starting point .
@xcite in the context of heat current between bosonic reservoirs this expression was derived in earlier works discussing hamiltonians similar to the one in eq .
( [ h - bosonic ] )  @xcite .
this type of expression survives even when accounting perturbatively for the interaction with vibrational modes .
@xcite      we begin with the low - temperature limit , @xmath125 ( throughout this paper we assume that @xmath126 ) . in this limit , the bose - einstein distribution functions @xmath127 and @xmath128 that enter the integrand of eq .
( [ j_q - general ] ) have decayed long before @xmath129 and @xmath119 have changed in any significant manner from their @xmath130 values .
thus , to a good approximation one can ( i ) set @xmath131 in the integrand of eq .
( [ j_q - general ] ) , and ( ii ) extend the upper integration limit to infinity . taking the coupling function @xmath117 to have the power - law form of eq .
( [ j - power - law ] ) this yields @xmath132 d \epsilon .\ ] ] the resulting integral can now be carried out analytically using @xmath133 where @xmath134 is the riemann zeta function .
@xcite this in turn gives @xmath135 with @xmath136 .
lastly , @xmath137 has the explicit expression @xmath138 } , \ ] ] allowing one to express the coefficient @xmath139 entirely in terms of the basic model parameters entering the bosonic hamiltonian of eq .
( [ h - bosonic ] ) .    as can be seen from eq .
( [ j_q - low - t ] ) , the low - temperature heat current shows a rather strong temperature dependence . in particular
, the linear - response heat conductance @xmath140 varies as @xmath141 , corresponding to @xmath142 for an ohmic bath .
this should be contrasted with the heat conductance of a generic noninteracting electronic tunnel junction , which varies linearly with @xmath143 at sufficiently low temperature .
@xcite for the heat conductance of a tunnel junction to display a superlinear temperature dependence of the form found here , its transmission coefficient must vanish in a power - law fashion at the fermi energy .
this power - law behavior of the heat current at low temperatures was also observed for systems in which the two bosonic reservoirs are connected by a system with few degrees of freedom such as in the spin - boson model  @xcite .
the similarity between these systems and the system under consideration here , at low - temperatures with respect to @xmath18 , stems from the fact that at this temperature regime only the lowest lying levels of the vibrational mode of the molecule are available for transferring energy between the baths .
as commented above , an ohmic bath with @xmath12 is of particular interest . focusing on this case and on the hierarchy @xmath144
, we devise below an analytical expression for the nonequilibrium heat current , encompassing the crossover from the low - temperature regime , @xmath145 , to the intermediate - temperature one , @xmath146 .
this expression is approximate as it employs a power series expansion of the self - energies , but yields accurate results within the regime where @xmath24 and @xmath10 are the largest energy scales in the system .    to this end , consider the function @xmath147 which enters the integrand of eq .
( [ j_q - general ] ) . since the bose - einstein distribution functions @xmath127 and @xmath128 decay on a scale far smaller than @xmath148 , it suffices to accurately represent @xmath147 for @xmath149 .
this allows one to expand the exact expressions for @xmath150 and @xmath151 [ see eqs .
( [ s^ra_e ] ) and ( [ s^ra_b - ohmic ] ) above ] in @xmath152 and @xmath153 to obtain @xmath154 , \\
\sigma^{r}_{b}(\epsilon ) & =       - 2 \pi \alpha \omega_c         \left [                i\pi y + 2                + { \cal o } \left (                                   y^2                           \right )         \right ] .\end{aligned}\ ] ] settling with linear orders in @xmath155 and @xmath156 , the function @xmath147 is well approximated for @xmath157 by @xmath158 where @xmath159 equals @xmath160 with the softened frequency @xmath161 ^ 2            } \label{omega}\ ] ] and the relaxation rate @xmath162 . \label{tau}\ ] ] the softened frequency @xmath163 characterizes the dressed excitations of the phonon , and thus serves as the energy scale which determines the cross over from the high - temperature to the low - temperature behavior . within this approximation for @xmath147 ,
the heat current for an ohmic bath becomes @xmath164 }                  { | ( \epsilon - z_{+ } )                     ( \epsilon - z_{- } ) |^2 } , \label{j_q - ohmic - int}\ ] ] where we have set in addition @xmath165 and extended the upper integration limit to infinity ( both approximations being well justified by the hierarchy @xmath166 ) .
the resulting integral in eq .
( [ j_q - ohmic - int ] ) can be performed in a closed analytic form in terms of the digamma function  @xcite @xmath167 .
skipping details of the algebra we quote here only the end result :    @xmath168                 + i\pi\tau \left (                                    t_{e } - t_{b }                            \right )        \right\ } .
\label{j_q - ohmic}\ ] ]    it is straightforward to confirm using the asymptotic expansion  @xcite @xmath169 that eq .
( [ j_q - ohmic ] ) properly reduces for @xmath170 to eq .
( [ j_q - low - t ] ) with @xmath12 , including the precise value of the prefactor @xmath139 .
more interesting is the limit @xmath171 , when the leading contribution to eq .
( [ j_q - ohmic ] ) crosses over to a linear dependence on the temperature difference @xmath172 .
thus , as the larger of the two temperatures exceeds the vibrational resonance energy , the heat current continues to increase linearly with @xmath173 .
this behavior markedly differs from that of a resonant electronic tunnel junction , whose heat current depends logarithmically on @xmath173 above the resonance energy .
the physical difference stems , as we show below , from the bosonic nature of the vibrational mode , which can be excited to exceedingly high energies by creating ever more phonons .
this should be contrasted with a resonant electronic level , which can only be empty or occupied .
having analyzed analytically certain limits , we now proceed to a complete numerical evaluation of the exact heat current of eq .
( [ j_q - general ] ) at arbitrary temperatures @xmath106 and @xmath84 . as the temperatures @xmath106 and
@xmath84 enter the expression for the heat current only through the term @xmath174 $ ] in the integrand , it is anti - symmetric under replacing them . as such
, we may restrict ourselves to calculations in which @xmath84 is held constant and the heat current is calculated for different values of @xmath106 , and generalize the results for opposite values by inverting the direction of the current .    in fig .
[ fig : j_q - low - t ] we have addressed the case where one bath ( the bosonic one ) is held at a constant low temperature @xmath175 , while scanning different values of @xmath106 higher than that temperature .
we have plotted the heat current in the case of a sub - ohmic bath ( with @xmath176 ) , an ohmic bath and a super - ohmic bath ( with @xmath177 ) .
the graph shows a clear crossover from a power - law behavior to a linear dependence as @xmath173 increases to values of the order of @xmath18 . in the inset of the graph we have plotted , on a log - log scale , the heat current for small values of @xmath173 for each type of bath , and a dashed line following the expected power - law behaviors at low temperatures given at eq .
( [ j_q - low - t ] ) .
the agreement between the calculated heat current and the expected one is excellent while @xmath178 , breaking at about @xmath179 .
figure  [ fig : j_q - high - t ] describes the heat current as well , but focuses on the case where the temperature of the bosonic bath is held at the high value of @xmath180 . in this case
the heat current about @xmath181 displays a linear dependence on @xmath173 , with different slope for every type of bosonic bath .    in fig .
[ fig : diff_conductance ] we have plotted the linear - response heat conductance @xmath182 , as defined in eq .
( [ conductance_defined ] ) , for different temperatures and types of baths .
as the temperatures rises , the conductance increases until it saturates at about @xmath183 .
the different exponents relating to the type of bosonic baths are evident in the low - temperature regime .     between the electronic and bosonic baths , as a function of the temperature difference between the baths . here
the temperature of the bosonic bath was held constant at a low value @xmath175 and the temperature of the electronic bath @xmath106 was changed .
the heat current was calculated for different types of bosonic baths , where we considered a sub - ohmic case ( @xmath176 , green ) , an ohmic case ( red ) and a super - ohmic case ( @xmath184 , blue ) .
inset : a log - log plot of the low - temperature regime .
the dashed lines are following the appropriate power - laws predicted in that regime by eq .
( [ j_q - low - t ] ) . here
@xmath124 is measured in units of the basic heat current @xmath185 , where @xmath186 , and @xmath187 are in the weak coupling regime .
the cutoffs used are @xmath188 . , scaledwidth=44.0% ]     and
we calculated the current for different values of @xmath106 .
the heat current was calculated for different types of bosonic baths , where we have considered a sub - ohmic case ( @xmath176 , green ) , an ohmic case ( red ) and a super - ohmic case ( @xmath184 , blue ) .
inset : a zoom over the regime about the point @xmath189 where the heat current displays a linear dependence on @xmath173 . here
@xmath124 is measured in units of the basic heat current @xmath185 , where @xmath186 , and @xmath187 are in the weak coupling regime .
the cutoffs used are @xmath188 .
, scaledwidth=44.0% ]     between the baths as a function of the temperature @xmath143 in which both bath are helds .
the conductance is calculated for different types of bosonic baths , where we have considered a sub - ohmic case ( @xmath176 , green ) , an ohmic case ( red ) and a super - ohmic case ( @xmath184 , blue ) .
here @xmath182 is measured in units of a basic heat conductance @xmath190 , where @xmath186 , and @xmath187 are in the weak coupling regime .
the cutoffs used are @xmath188 .
, scaledwidth=44.0% ]
although our solution for @xmath124 is formally exact within the bosonic hamiltonian of eq .
( [ h - bosonic ] ) , it is instructive to develop a more transparent physical picture that would , in particular , elucidate the source of distinction between the bosonic system under consideration and a conventional resonant tunnel junction . to this end , we devise below a master equation approach , applicable at weak coupling . this approach to address heat flow was previously applied by segal  @xcite , considering a different setup .
leijnse and others  @xcite have also used this method to calculate the thermopower properties in a similar setup .
the basic components of the theory are the probabilities @xmath191 to find the local boson @xmath3 at time @xmath19 in the state where @xmath192 .
these probabilities are connected at the golden - rule approximation by rate equations of the form @xmath193 ( terms with @xmath194 should be omitted for @xmath195 ) , with the transition rates  @xcite @xmath196 , \\
w_{n \to n - 1 } & = 2\pi n       \left\ {                f(\omega_0 ) [ 1 + n_e(\omega_0 ) ]       \right
. \nonumber \\     & \;\;\;\;\;\;\;\;\;\;\;\ ;       \left .
+ \lambda(\omega_0 ) [ 1 + n_b(\omega_0 ) ]       \right\ } .\end{aligned}\ ] ] here , @xmath197 equals @xmath198 .
the golden - rule approximation used here corresponds to lowest order perturbation theory in both the coupling to the phononic bath [ represented by the coupling function @xmath199 and the coupling to the bosonic bath [ represented by the coupling function @xmath117 ] .
it is useful to define at this point the reduced transition rates @xmath201 , \\
w_{\downarrow } & = 2\pi       \left\ {                f(\omega_0 ) [ 1 + n_e(\omega_0 ) ]                + \lambda(\omega_0 ) [ 1 + n_b(\omega_0 ) ]       \right\ } , \end{aligned}\ ] ] such that @xmath202 and @xmath203 . focusing on steady state ,
when @xmath204 , eq .
( [ master - equations ] ) can be recast as an infinite set of coupled linear equations , @xmath205 whose solution is @xmath206 with @xmath207 . here
, @xmath208 is a normalization factor which comes to ensure that @xmath209 .
thus , the probabilities @xmath210 obey a boltzmann - like form with an effective temperature @xmath211 defined by @xmath212    consider first the case of thermal equilibrium , when @xmath213 and @xmath214 . under these circumstances
@xmath215 hence @xmath211 equals @xmath143 irrespective of details of the two baths .
once a temperature gradient is applied between the two reservoirs , i.e. , @xmath216 , then @xmath211 falls in between @xmath217 and @xmath218 , as follows from the equality @xmath219 with @xmath220 }           { f(\omega_0 ) [ 1 + n_e(\omega_0 ) ]            + \lambda(\omega_0 ) [ 1 + n_b(\omega_0 ) ] } .
\label{q - def}\ ] ]    equation  ( [ t_eff - via - exponents ] ) can be solved analytically in the high - temperature limit , @xmath221 , where fermi s golden rule ( and thus our master equation approach ) is expected to apply .
specifically , expanding each of the exponents to linear order in @xmath18 one obtains @xmath222 this expression can further be simplified by noting that @xmath223 depends implicitly on @xmath106 and @xmath84 through the bose - einstein distribution functions @xmath224 and @xmath225 . plugging these relations into eq .
( [ q - def ] ) and inserting the resulting expression for @xmath223 into eq .
( [ t_eff - via - q ] ) , one finally arrives at @xmath226 where we have made repeated usage of the fact that @xmath221 .
it should be stressed that this result equally applies to all forms of the bosonic bath , be it ohmic , sub - ohmic or super - ohmic .
the steady - state solution to the probablities @xmath210 can be used in turn to calculate the heat current .
focusing again on the heat current flowing between the bosonic bath and the local phonon , the latter involves the transition rates @xmath227 to and from bosonic bath .
explicitly , the heat current takes the form @xmath228 where @xmath229.\end{aligned}\ ] ] using these rates the expression for the heat current gains the compact form @xmath230 ,      \label{j_q - master - eq}\ ] ] where @xmath231 is the average occupancy of the localized phonon . recalling that @xmath211 lies between the temperatures of the bosonic and electronic baths ,
it is clear that @xmath124 is positive ( negative ) when @xmath122 ( @xmath232 ) , which gives the correct direction of the heat flow .
the above result provides a transparent picture for the linear dependence of the heat current on the temperature gradient in the high - temperature regime . in this regime
, the expression for @xmath124 can be further approximated using eq .
( [ t_eff_high_t ] ) for the effective temperature , which results in @xmath233 hence , the linear dependence on the temperature gradient stems from the large occupancy of the localized phonon , which does not saturate with increasing temperature .
segal  @xcite , using the master equation approach , described a similar linear dependence of the heat flow on the temperature difference .
this should be contrasted with a spinless electronic resonant level , that can only be empty or singly occupied .
while eq .
( [ j_q - master - eq ] ) properly captures the physics of the high - temperature regime , it fails to produce the required power - law behavior in the low - temperature regime , @xmath234 .
indeed , replacing the bose - einstein functions with simple exponents and using eq .
( [ t_eff - via - exponents ] ) for the effective temperature , the heat current becomes @xmath235 where we have omitted the exponentially small bose - einstein distribution functions in the expression for @xmath223 [ see eq .
( [ q - def ] ) ] .
thus , the master equation approach predicts activated low - temperature behavior in place of the correct power - law form .
the fact that the master equation approach well describes the high - temperature regime but fails at low temperatures is by itself not surprising .
here we can exploit the exact solution to carefully examine the range of validity of the approach and the simple physical picture that it lends .
generally speaking , the quality of the master equation approach depends on two parameters : ( 1 ) the temperatures involved , and ( 2 ) the strength of the coupling constants .
naturally , the master equation approach loses quantitative accuracy as the coupling constants are increased since fermi s golden - rule corresponds to lowest - order perturbation theory . nevertheless ,
as we shall show , the qualitative picture can remain quite accurate even when the coupling constants are increased . to critically examine the master equation approach we shall focus on the three quantities : ( 1 ) the effective temperature @xmath211 defined by the average occupancy of the localized phonon [ see eq .
( [ t_eff_from_exact ] ) below ] ; ( 2 ) the boltzmann - like distribution of the localized phonon , and ( 3 ) the heat current @xmath124 .
as shown above , the master equation approach predicts an effective temperature @xmath211 which is manifest in the average phonon occupancy @xmath236 .
similarly , one can define an effective temperature from the exact phononic occupancy @xmath237 according to @xmath238 where @xmath163 is the softened phonon frequency , approximated by eq .
( [ omega ] ) .
usage of the softened frequency , rather than the bare one , comes to account for higher - order correction not included in the golden - rule approximation .
we emphasis that eq .
( [ t_eff_from_exact ] ) serves as an ad - hoc definition of an effective temperature , which does not , by itself , imply a boltzmann - like distribution
. we shall examine this latter point in the following subsection .
figure [ fig : t_eff ] compares the effective temperature extracted from the exact solution with the master equation result of eq .  (
[ t_eff - via - exponents ] ) , scanning different coupling strengths .
there is an excellent agreement at weak coupling which gradually deteriorates as the combined coupling to the two baths is increased .
the agreement is controlled by the decay time @xmath239 of eq .
( [ tau ] ) , which depends quadratically on the coupling @xmath0 to the electronic bath and linearly on the coupling @xmath11 to the bosonic one . for the strongest coupling considered ( right - hand side of the lower - right panel )
the deviation between the two curves is of the order of @xmath240 , which is still quite moderate .     and for different values of the coupling to the bosonic reservoir @xmath11 .
the red line represents the effective temperature calculated based on the master equation approach while the black line is the effective temperature extracted from the phononic occupancy as given in eq .
( [ t_eff_from_exact ] ) . in this graph @xmath241 ,
the bosonic bath was taken to be ohmic and the temperatures used are @xmath242 , @xmath243 , which lie well within the applicatbility temperature range of the master equation approximation . , scaledwidth=44.0% ]
the boltzmann - like distribution of the probablities @xmath210 predicts all the moments of the phononic occupancy @xmath244 .
the ratio between the moments can be expressed using the moments themselves , giving the ratio between the first two moments @xmath245 which does not directly depend on the effective temperature , and is true for every boltzmann - like distribution , regardless of the value of the temperature .
we shall use this ratio as a benchmark for examining how close the phononic distrubtion is to a boltzmann - like one , as the exact ratio is calculated using the phononic green functions at steady state @xmath246 +      \left|g^{<}_{bb}(t , t)\right|^2      \right\}. \nonumber\end{aligned}\ ] ] the calculated ratio deviates from the boltzmann - like one due to the presence of the term @xmath247 , which is identically zero for thermal distributions .    in fig .
[ fig : is_thermal ] we have plotted the exact ratio between the second and first moments of the phononic occupation and compared it to the expected ratio were the distribution was boltzmann - like . similarly to the comparison done for the effective temperature , the comparison was done as a function of the coupling strength . for small couplings
the exact ratio matches excellently the one expected from a boltzmann - like distribution , thus confirming the approximation .
as the coupling is increased , the exact ratio deviates from the botlzmann - like prediction , and the non - thermal nature of the distribution becomes prominent .
we feel it is important to stress , however , that the distribution is not thermal for any coupling strength .      as a third quantity to which we look in order to examine the validity of the master equation approximation
is the heat current between the two baths @xmath124 .
the exact expression for the heat current is given in its integral form in eq .
( [ j_q - general ] ) and the expression derived from the master equation approach given in eq .
( [ j_q - master - eq ] ) .
the comparison of the results of these two expressions is plotted in fig .
[ fig : j_q_master_eq ] where we have scanned different temperatures while holding the temperature difference between the baths constant . in accordance with our previous analysis of the results , the master equation reproduces quite well the exact result only at high temperatures , and deviating from it significantly as the temperatures are lowered .     as a function of the electron - phonon coupling @xmath0 and for different values of the coupling to the bosonic reservoir @xmath11 .
the red line represents the ratio as predicted by the master equation while the black line is the actual ratio calculated according to the exact solution . in this graph @xmath241 , the bosonic bath
was taken to be ohmic and the temperatures used are @xmath242 , @xmath243 , which lie well within the applicability temperature range of the master equation approximation . , scaledwidth=44.0% ]     and the approximate calculation relying on the master equation approach as given in eq .
( [ j_q - master - eq ] ) .
the temperature difference between the reservoirs is held fixed with @xmath248 and both temperatures are changed .
inset : the ratio of the approximated result from the exact one , showing a larger difference as the temperature is lowered . here
@xmath124 is measured in units of the basic heat current @xmath185 , where @xmath186 , and @xmath249 are in the weak coupling regime .
the cutoffs used are @xmath188 .
, scaledwidth=44.0% ]
so far only the resonant case where @xmath16 was considered . since @xmath13 can be controlled by using suitable gate voltages , studying its effects on the heat current is of particular interest . in this section
we address the case where @xmath13 is given some nonzero value , breaking the particle - hole symmetry of the model .
staying in accordance with the rest of our discussion , where we have assumed that the bandwidth @xmath24 is the largest energy scale of the system , we are interested in the regime where @xmath250 in addition to @xmath251
. we shall demonstrate that in this regime the leading contribution of @xmath13 to the heat current is either @xmath252 or @xmath253 .    before getting to the task of explicitly calculating the heat current for the off - resonant case , we first consider the symmetries of the hamiltonian of eq .
( [ h ] ) . if @xmath16 then the system holds a particle - hole symmetry , which breaks for nonzero values of @xmath13 . writing the hamiltonian as a function of @xmath13 and carrying out an particle - hole transformation @xmath254 one finds out that @xmath255 , while the expression for @xmath124 remains unchanged .
thus we conclude that @xmath124 is an even function of @xmath13 .
next we turn to the calculation done in sec .
[ sec : resonance_heat_current ] under resonance conditions and aim at adjusting it to the case where @xmath256 . since the @xmath13 term in the hamiltonian of eq .
( [ h - bosonic ] ) is linear in the operators @xmath257 and @xmath258 , one may express its effects on the heat current by appropriate corrections to the bare green functions pertaining to these operators
. we will then be able to use the functions , dressed by the @xmath13 term , to re - calculate the self - energies @xmath259 of eq .
( [ s^ra ] ) and @xmath260 of eq .
( [ s^ < > ] ) . to this end
, we consider a free hamiltonian to which we add an @xmath13 term @xmath261 where @xmath262 .
we will designate by @xmath263 the unperturbed green functions given , in energy space , by @xmath264 , \label{g0_ra_epsilon_d}\ ] ] @xmath265 , \label{g_les_grea_epsilon_d}\ ] ] and by @xmath266 the dressed functions with respect to the hamiltonian .
the hamiltonian can be diagonalized exactly by introducing new bosonic creation and annihilation operators @xmath267 , by which it takes the form @xmath268 this is a free hamiltonian with respect to the @xmath269 s .
the correlation functions between the original bosonic operators are related to the correlation functions of the @xmath269 s by the fact that @xmath270 we can thus adjust the lesser and greater green functions , in the energy domain , by adding the appropriate term proportional to @xmath271 , which reflects the fact that the extra term added by @xmath13 on the right - hand - side of eq .
( [ nonresonant_new_bosonic ] ) is independent of time @xmath272 we continue to note that the retarded and advanced green functions remain unchanged by this addition of non - zero @xmath13 , as the extra term is constant and drops out when the commutation relations are taken .
having arrived at the conclusion that the corrections due to the @xmath13 term exist only at @xmath273 , we point to the fact that the integrand in the expression for the heat current given in eq .
( [ j_q - general ] ) depends directly on @xmath274 .
therefore , all such terms do not contribute to the heat current , leaving it independent of the value of @xmath13 . as @xmath13 reflects the energy associated with the charging of the level
, one would expect that its value will directly affect the heat current .
we understand the independence of the latter on the value of @xmath13 as a result of the weak - coupling regime @xmath275 , and expect that as either of these values is increased , @xmath13 will play a role in determining the heat current . taking into account that the heat current is an even function of @xmath13 , we conclude our discussion in noting that the leading contribution in this regime is not lower than @xmath276 or @xmath277 .
in this paper we have presented an asymptotically exact calculation of the heat current between a bosonic bath and a fermionic bath , that is mediated by a single molecule .
the calculation is based on a mapping of continuous model given by the hamiltonian of eq .
( [ h ] ) onto a form quadratic in bosonic operators  @xcite .
this mapping is valid in the weak coupling regime , where @xmath278 , and @xmath24 is the effective electronic bandwidth .
this model may describe , under suitable mappings , several physically and experimentally relevant setups , the most relevant being a molecule adsorbed on a surface , a molecular junction and an aharonov - bohm interferometer with a molecular device embedded in one of its arms .
the exact calculation yields a landauer - type expression for the heat current , given in eq .
( [ j_q - general ] ) .
such an expression stands in accordance with previous works on thermal currents in confined nanostructures .
@xcite    at low temperatures , the heat current strongly depends on the nature of the bosonic bath . assuming that the bosonic bath has a power - law form and is characterized by the power @xmath279 , the low - temperature linear - response heat conductance varies as @xmath280 . at high temperatures , however , the heat current depends linearly on the temperature difference between the two baths , regardless of the nature of the power - law governing the bosonic bath .
the crossover between the low and high - temperatures regimes is at the scale of the softened vibrational mode frequency @xmath163 given in eq .
( [ omega ] ) .
the high - temperature behavior , which is markedly different than the transmission through a purely electronic system , is explained by the bosonic nature of vibrational mode , which can be excited to high energies by creating more phonons .
this is illustrated by a master equation analysis of the system , which is perturbative in nature but is justified in the high - temperature and weak coupling regime . in that regime
the master equation approach reproduces the heat current that was calculated exactly previously , and also offers an effective temperature that we assign to the local vibrational mode .
it should be stressed that even at that regime it does not have a thermal distribution , and the effective temperature is a useful illustrative approximation .    as our solution is exact only at weak electron - phonon coupling
, one would expect different features to appear in the heat current as the interaction strength is increased
. it would be interesting to compare our results with such an analysis , and see how our calculation persists into the strong - coupling regime .
we leave that to future work .
99 for recent reviews see , e.g. , _ introducing molecular electronics _ , edited by g.  cuniberti , g.  fagas , and k.  richter , lecture notes in physics vol .
680 ( springer , new york , 2005 ) ; m.  galperin , m.  a.  ratner , and a.  nitzan , j.  phys . : condens .
matter * 19 * , 103201 ( 2007 ) .
wang , j.  wang , and j.  t.  l , eur .
j.  b * 62 * , 381 - 404 ( 2008 ) .
d.  g.  cahill , w.  k.  ford , k.  e.  goodson , g.  d.  mahan , a.  majumdar , h.  j.  maris , r.  merlin and s.  r.  phillpot , j.  appl .  phys .  * 93 * , 793 ( 2003 ) .
d.  r.  ward , d.  a.  corley , j.  m.  tour and d.  natelson , nature nanotechnology * 6 * , 33 - 38 ( 2011 ) m.  galperin , a.  nitzan , and m.  a.  ratner , phys .
b * 75 * , 155312 ( 2007 ) .
o.  entin - wohlman , y.  imry and a.  aharony , phys .
b * 82 * , 115314 ( 2010 ) .
o.  entin - wohlman and a.  aharony , phys .  rev .
b * 85 * , 085401 ( 2012 ) .
a.  ozpineci and s.  ciraci , phys .
b * 63 * , 125415 ( 2001 ) .
d.  segal , a.  nitzan , and p.  hnggi , j.  chem .
* 119 * , 6840 ( 2003 ) .
l.  a.  wu and d.  segal , phys .
e * 83 * , 051114 ( 2011 ) . b.  sothmann , r.  sanchez , a.  n.  jordan and m.  bttiker , phys .
b. * 85 * , 205301 ( 2012 ) .
b.  sothmann and m.  bttiker , epl * 99 * , 27001 ( 2012 ) .
l.  i.  glazman and m.  raikh , jetp lett .
* 47 * , 452 ( 1988 ) .
g.  schulze , k.  j.  franke , a.  gagliardi , g.  romano , c.  s.  lin , a.  l.  rosa , t.  a.  niehaus , th .
frauenheim , a.  di carlo , a.  pecchia and j.  i.  pascual , phys .
* 100 * , 136801 ( 2008 ) ; g.  schulze , k.  j.  franke and j.  i.  pascual , new j.  phys .
* 10 * , 065005 ( 2008 ) ; k.  j.  franke and j.  i.  pascual , j.  phys . :  condens .
matter * 24 * , 394002 ( 2012 ) .
m.  tsutsui , m.  taniguchi , k.  yokota and t.  kawai , appl .
lett .  * 96 * , 103110 ( 2010 ) .
d.  segal , phys .
lett .  * 100 * , 105901 ( 2008 ) .
t.  ojanen and a.  p.  jauho , phys .
100 * , 155902 ( 2008 ) .
t.  ruokola and t.  ojanen , phys .
b. * 83 * , 045417 ( 2011 ) .
a.  dhar and d.  roy , j.  stat .
* 125 * 801 ( 2006 ) .
a.  dhar , adv .
phys .  * 57 * , 457 ( 2008 ) .
k.  saito and t.  kato , arxiv:0706.1234 we employ units in which @xmath281 .
a.  j.  legget , s.  chakravarty , a.  t.  dorsey , m.  p.  a.  fisher , a.  garg , and w.  zwerger , rev .
phys .  * 59 * , 1 ( 1987 ) .
d.  sherrington and s.  von molnr , solid st . comm . * 16 * , 1347 ( 1975 ) .
r.  r.  cavanagh , e.  j.  heilweil and j.  c.  stephenson , surface sciente * 299 * , 643 ( 1994 ) .
j.  w.  gadzuk , phys .
b * 24 * , 1651 ( 1981 ) ; e.  blaisten - barojas and j.  w.  gadzuk , j.  chem .  phys .
* 97 * , 862 ( 1992 ) .
y.  vinkler , a.  schiller , and n.  andrei , phys .
b * 85 * , 035411 ( 2012 ) .
h.  l.  edwards , q.  niu and a.  l.  de lozanne , appl .
lett .  * 63 * , 1815 ( 1993 ) ; h.  l.  edwards , q.  niu , g.  a.  georgakis and a.  l.  delozanne , phys .
b * 52 * , 5714 ( 1995 ) .
j.  r.  prance , c.  g.  smith , j.  p.  griffiths , s.  j.  chorley , d.  anderson , g.  a.  c.  jones , i.  farrer and d.  a.  ritchie , phys .
lett .  * 102 * , 146602 ( 2009 ) .
i.  g.  lang and yu .
a.  firsov , zh .  eksp .
* 43 * , 1843 ( 1962 ) [ sov .
jetp * 16 * , 1301 ( 1963 ) ] .
b.  dra and a.  halbritter , phys .
b. * 80 * , 155402 ( 2009 ) .
b.  dra , phys .
b. * 75 * , 245113 ( 2007 ) .
b.  dra and m.  gulcsi , phys .
b. * 78 * , 165111 ( 2008 ) .
w.  hofstetter , j.  knig and h.  schoeller , phys .
87 * , 156803 ( 2001 ) .
l.  v.  keldysh , sov .
jetp * 20 * , 1018 ( 1965 ) . f.  d.  m.  haldane , j.  phys .
c * 14 * , 2585 ( 1981 ) .
d.  c.  langreth , in _ linear and nonlinear electron transport in solids _
( plenum pres , new york , 1976 ) , vol .
17 of _ nato advanced study institute , series b : physics
_ , eds .
j.  t.   devreese and v.  e.  van doren .
see , e.g. , _ handbook of mathematical functions _ , eds . m.  abramowitz and i.  a.  stegun ( dover , new york , 1972 ) , chapter 5 .
see , e.g. , _ handbook of mathematical functions _
, eds . m.  abramowitz and i.  a.  stegun ( dover , new york , 1972 ) , chapter 23 .
f.  zhan , s.  denisov and p.  hnggi , phys .
b * 84 * , 195117 ( 2011 ) .
see , e.g. , _ handbook of mathematical functions _
, eds . m.  abramowitz and i.  a.  stegun ( dover , new york , 1972 ) , chapter 6 .
d.  segal , phys .
b * 73 * , 205415 ( 2006 ) .
m.  leijnse , m.  r.  wegewijs and k.  flensberg , phys .
b * 82 * , 045412 ( 2010 ) .