# The effect of dephasing on thermoelectric efficiency of molecular junctions

###### Abstract

In this work we report the results of theoretical analysis of the effect of thermal environment on the thermoelectric efficiency of molecular junctions. The environment is represented by two thermal phonon baths associated with the electrodes which are kept at different temperatures. The analysis is carried out using the Buttiker model within the scattering matrix formalism to compute electron transmission through the system. This approach is further developed, so that the dephasing parameters are expressed in terms of relevant energies including the thermal energy, strengths of coupling between the molecular bridge and the electrodes and characteristic energies of electron-phonon interactions. It is shown that the latter significantly affect thermoelectric efficiency by destroying coherency of the electron transport through the considered system.

## I i. introduction

For the past three decades, significant efforts were applied to study thermal and thermoelectric transport in mesoscopic and nanoscale systems of various kinds including quantum dots and/or molecules attached to conducting electrodes. The latter serve as source and drain reservoirs for traveling charge carriers. Here, we concentrate on the analysis of thermoelectric properties of these systems. Below they are referred to as thermoelectric junctions. In part, the research interest appeared because thermoelectric junctions are expected to be useful in building up highly efficient energy conversion devices. Also, studies of thermoelectric properties of these systems can result in a deeper insight into the nature of general transport mechanisms and bring additional information on the electronic and vibrational excitation spectra of molecules 1 (); 2 (); 3 (); 4 (); 5 (); 6 (); 7 (); 8 (); 9 (); 10 (); 11 (); 12 (); 13 (); 14 (); 15 (). Presently, it is established that thermoelectric properties of quantum dots and/or molecules may be strongly affected by Coulomb interactions between charge carriers 16 (); 17 (); 18 (); 19 (); 20 (); 21 (); 22 (). Coulomb interactions lead to violation of the Wiedemann-Franz law in nanoscale thermoelectric junctions thus providing an enhancement of thermoelectric efficiency of these systems 21 (); 23 (). The thermal efficiency may be further increased due to the influence of quantum interference effects which may strongly affect electron transport characteristics 9 (); 24 (); 25 (); 26 (); 27 (); 28 (); 29 ().

In general studies of thermoelectric transport through molecules, quantum dots and similar systems one must imply that both atomic vibrations and charge carriers contribute to the energy transfer. Therefore, a unified description of the electrons and phonons dynamics is needed to thoroughly analyze thermoelectric properties of molecular junctions. The means for such analysis are provided by the nonequilibrium Green’s functions formalism (NEGF), as described in the review 30 (), and some other works (see e.g. 13 (); 31 (); 32 (); 33 (); 34 (); 35 ()). However, application of this formalism to realistic models simulating thermoelectric junctions is extremely difficult. Several simplified approaches were developed and used to analyze the specifics of heat transfer and other related phenomena in quantum dots and molecules taking into account the contribution of phonons and electron-phonon interactions 2 (); 13 (); 36 (); 37 (); 38 (). Nevertheless, these studies are not completed so far.

The phonons contributing to the charge and energy transfer may be subdivided in two classes: vibrational phonons associated with molecular vibrations and thermal phonons associated with random nuclear motions in the environment. In the present work we aim at theoretical analysis of the effect of thermal phonons on the thermoelectric characteristics of molecules and other similar systems. To carry on this analysis we combine NEGF with the approach first suggested by Buttiker to describe quantum transport through molecules 39 (). An important advantage of this approach is that it could be easily adapted to analyze various aspects of incoherent/inelastic transport through molecules (and some other mesoscopic systems) avoiding complicated and time-consuming methods based on more advanced formalisms.

## Ii ii. Main equations

For simplicity, in the following computations we simulate a molecule/quantum dot by a single level with the energy We assume that this single-level bridge is coupled to the pair of dissipative reservoirs, as shown in the Fig. 1. While on the bridge, an electron could be scattered into one of the reservoirs through the channels 3 and 4 (or 5 and 6) with a certain probability. In the reservoir, it undergoes inelastic scattering accompanied by phase breaking, and afterwards returns to the bridge with the same probability. In the present analysis, the reservoirs are treated as phonon baths representing thermal phonons associated with the left and right electrodes. Within the accepted model we imply that there is no phonon thermal conductance through the junction. This seems a reasonable assumption for experiments give low values of phonon thermal conductance in several thermoelectric molecular junctions 40 (); 41 (). This may be attributed to the fact that in many molecules the majority of vibrational transitions lie above the range determined by thermal energy when temperature takes on values of the order of or lower than the room temperature 14 (); 29 (). Within the Buttiker model, the electron transport through a thermoelectric junction is considered as combination of tunnelings through potential barriers separating the electrodes from the molecule/quantum dot and interaction with the reservoirs coupled to the bridge site.

The Buttiker approach was applied to describe and analyze electron transport through molecules in several works (see e.g. Refs. 42 (); 43 (); 44 (). Following this approach, one can present particle fluxes outgoing from the system as linear combinations of incoming fluxes where the indexes label the channels for transport. For the adopted model

(1) |

In these equations, the coefficients are related to matrix elements of the scattering matrix namely: The scattering matrix expresses outgoing wave amplitudes in terms of incident ones To provide the charge conservation in the system, zero net current should flow in the channels linking the bridge site with the dephasing reservoirs, so we may write the following equations:

(2) |

To find the expressions for the matrix elements we first consider a subsystem including the left electrode with the associated dephasing reservoir and the bridge site. The expression for the matrix relating the wave amplitudes to the wave amplitudes has the form 39 ():

(3) |

where the phenomenological scattering probability corresponds to the reservoir associated with the left electrode. Also, an electron tunneling through a single potential barrier separating this electrode from the molecule/quantum dot is characterized by the transmission and reflection amplitudes and respectively). These are the matrix elements of a matrix:

(4) |

Combining Eqs. (3) and (4) one obtains the expression for the scattering matrix relating to :

(5) |

Here,

Now, we take into consideration the remaining elements of the original system. The matrix which relates to is 43 (); 44 ():

(6) |

where the scattering probability is associated with the right reservoir, and the transmission and reflection amplitudes characterize electron tunneling through the potential barrier between the molecule (bridge) and the right electrode. Using Eqs. (5) and (6) and excluding the wave amplitudes which correspond to the transport inside the system, we get the following expression for the scattering matrix:

(7) |

Here,

Solving the system of linear equations (1), (2) one obtains the following expression for the electron transmission which coincides with the corresponding result reported by D’Amato and Pastawski 42 ():

(8) |

Within the considered model,

(9) |

and is the matrix inversed with respect to matrix whose matrix elements are given by:

(10) |

In this expression, the following denotations are used:

(11) |

Assuming that both dephasing reservoirs are detached from the bridge the transport through the system becomes coherent and elastic. In this case, the electron transmission given by Eqs. (8)-(11) is reduced to a simple form:

(12) |

As known, the expression for the electron transmission in the case of coherent transport may be presented as follows:

(13) |

where are the retarded and advanced Green’s functions associated with the molecule/quantum dot bridging the electrodes, and self-energy terms describe the coupling between the electron of a certain spin orientation on the bridge and the corresponding electrode. For a symmetrically coupled system the expression for electron transmission may be reduced to the form:

(14) |

Provided that electron transport through the system is undisturbed by electron-phonon interactions, and disregarding spin-flip processes, the retarded Green’s function may be approximated as 45 ():

(15) |

Here, is the charging energy associated with Coulomb repulsion between the electrons on the molecular bridge/quantum dot and are one-particle occupation numbers:

(16) |

where is the lesser Green’s function for electrons on the bridge. Self-energy corrections appear in the expression for due to the coupling of the bridge to the electrodes. For example:

(17) |

In this expression, are single-electron energies on the electrode are coupling parameters characterizing the coupling of the electron states on the bridge to the electrodes and is an infinitesimal positive parameter. These self-energy terms are closely related to the previously introduced coupling strengths namely: The expressions (14)-(16) were repeatedly employed in studies of thermal transport through quantum dots (see e.g Ref. 26 (); 46 ()).

For a symmetrically coupled system, one may assume that the potential barriers separating the electrodes from the bridge are identical: Then the transmission amplitude could be easily expressed in terms of the corresponding Green’s functions:

(18) |

Within the Buttiker approach, the scattering probabilities are introduced as phenomenological parameters. However, these parameters may be given an explicit physical meaning by expressing them in terms of the relevant energies. In the considered system, dissipation and loss of coherency appear due to the interaction of charge carriers with thermal phonons associated with the electrodes and represented by the dephasing reservoirs. Therefore, as was suggested in an earlier work 47 (), one can approximate these parameters as follows:

(19) |

Here, represents the self-energy term originating from electron-phonon interactions occurring in the reservoir associated with the left/right electrode. Using NEGF and computing the relevant electron and phonon Green’s functions within the self-consistent Born approximation, one can arrive at a relatively simple expression for 30 ():

(20) |

In this expression, and are respectively the electron density of states associated with the bridge level and its steady state occupation, and is the phonon spectral function for the corresponding reservoir. We assume that the electrodes may be kept at different temperatures so we introduce phonon distribution functions where is the Boltzmann constant. Finally, the constant characterizes the coupling strength for electron interactions with the thermal phonons belonging to the bath

The particular form of the phonon spectral functions may be found basing on the molecular dynamic simulations. However, to qualitatively analyze the effect of dephasing on the thermoelectric transport, one may employ the approximation 48 ():

(21) |

where the parameter is related to the electron-phonon coupling strength, and characterizes the relaxation time for the thermal phonons.

The electron density of states includes self-energy corrections which appear due to electron-phonon interactions. Therefore, Eq. (20) is an integral equation for Substituting the approximation (17) into this equation, one may see that the major contribution to the integral over originates from the region where Omitting the terms in the arguments of all slowly varying terms in the integrand, we may reduce Eq. (20) to the form:

(22) |

where

(23) |

and is the Riemann’s function.

The suggested approach gives means to theoretically analyze the effect of thermal phonons on the thermoelectric properties of thermoelectric junctions. Using the obtained results given by Eqs. (7)-(21), one may compute electron transmission implying that the difference in the temperatures and can take on an arbitrary value. Therefore, these results may be employed to study thermoelectric properties of the considered systems beyond the linear regime. As known, nonlinear thermoelectric properties of molecular junctions and similar systems presently attract significant interest 33 (); 37 (); 46 (); 49 (); 50 (). However, in studies of thermoelectric characteristics of such systems beyond linear regime, one inevitably encounters a nontrivial task of introducing and defining the local temperature for the bridge which differs from temperatures associated with the electrodes. The definition of local temperature and related problems are thoroughly discussed in the recent review 49 ().

In the present work we avoid these difficulties by restricting further analysis with the linear temperature and bias regime. Also, we remark again that within the considered model the thermal conductivity associated with phonons is omitted for we do not include into consideration vibrational modes coupled to the bridge. Therefore, we may employ the following commonly used expressions for measurable thermoelectric characteristics:

(24) | ||||

(25) |

Here, and are electron electrical and thermal conductances, is the thermopower (Seebeck coefficient) and is the dimensionless thermoelectric figure of merit characterizing the efficiency of charge-driven cooling devices and/or heat-driven current generators. In deriving these expressions, it is assumed that and The integrals included in Eqs. (24), (25) are given by:

(26) |

where is the Fermi distribution function for the energy and the chemical potential characterizes the electrodes at zero bias. Coulomb interactions between electrons on the molecule/quantum dot may be accounted for by using appropriate expressions for the electron Green’s functions incorporated into the expression for the coherent transmission (13) as well as into the expression for the electron density of states

## Iii iii. Results and discussion

Specific thermoelectric properties of the considered systems depend on the relation of four relevant energies. These are the strength of coupling of the bridge to the electrodes the electron-phonon coupling strength the charging energy characterizing Coulomb interactions of electrons on the bridge, and the thermal energy It was established that the greater values of could be achieved in weakly coupled systems where the condition may be satisfied at reasonably low temperatures (see e.g. Ref. 21 ()), so in further analysis we assume that considered system complies with this condition.

Also, we assume that the considered quantum dot/molecule is symmetrically coupled to the electrodes and two thermal baths are identical Omitting for a while Coulomb interactions, one may derive a simple Lorentzian expression for the electron density of states:

(27) |

where Substituting this expression into Eq. (22), we may solve this equation and arrive at a reasonable asymptotic expression for

(28) |

where Using this result and the expression (19) for the scattering probabilities, we obtain:

(29) |

The parameter values vary between 0 and 1. When the bridge is detached from the reservoirs, and the electron transport is completely coherent and elastic. Within the opposite limit the transport is characterized by the overall phase randomization typical for inelastic sequential hopping. Within the adopted approach, the scattering probabilities depend on tunnel energy As well as electron transmission function they reach their maximum values and respectively) at This is shown in the left panels of the Fig. 2.

As follows from Eq. (29), the character of the electron transport is determined by the value of the dimensionless parameter Transport remains nearly coherent when On the contrary, the strong dephasing/dissipation occurs when takes on values significantly greater than 1. To find a suitable estimate for one needs to approximate the Riemann’s function included into expression (23). The approximation depends of the relation between the energies and As discussed in an earlier work 43 (), the effect of the thermal bath on the electron transport is significantly more pronounced when the lifetime of thermal excitations is sufficiently long Under this condition, one may apply the estimation Correspondingly, This shows that the maximum value of the scattering probabilities is determined with two parameters, namely, and We remark that in the absence of electron-phonon interactions regardless of the energy value, and In general, increases when the temperature rises, and it takes on greater values when the electron-phonon interactions are getting stronger as illustrated in the Fig. 2. The enhancement of is accompanied by the decrease of maximum value of electron transmission These results have an obvious physical sense because in the considered situation the phase randomization is inseparable from inelastic scattering of electrons by thermal phonons hindering electron transport through the system.

It was first shown by Sofo and Mahan 51 () and then confirmed in several later works (see e.g. Ref. 21 ()) that the figure of merit diverges when approaches zero provided that the system is characterized by zero phonon contribution to the thermal conductance, and the effects of Coulomb interactions between electrons on the bridge are disregarded. The results obtained in the present work agree with this conclusion. In the right bottom panel of the Fig. 3, the divergence of within the limit in the absence of the electron-phonon interactions is clearly illustrated. The junction figure of merit is limited due to the effect of thermal phonons associated with the electrodes. The stronger these interactions are, the lower in magnitude maximum values of both and thermopower become. Comparing the present results with those reported in Ref. 21 () one may presume that the thermal phonons take the part of phonon thermal conductance (which equals zero for the considered system) in limiting the maximum value of and removing the divergence. Unlike the thermopower remains finite at small values of even when the electron-phonon interactions are disregarded, as illustrated in the Fig. 3 (see right top panel). This means that the divergence of originates from extremely strong violation of the Wiedemann-Franz law resulting in the divergence of the Lorentz ratio.

In analyzing temperature dependencies of thermoelectric characteristics of systems consisting of a molecule/quantum dot linking two electrodes, it was established that usually the figure of merit is a nonmonotonous function of temperature (see e.g. Refs. 9 (); 19 (); 20 (); 21 (); 25 (); 46 (). Also, the present results show that at low temperatures increases as the temperature enhances and it reaches a maximum value at certain temperature As further rises, the figure of merit decreases approaching zero when the temperature significantly exceeds This is illustrated in the Fig. 4. The value of the optimal temperature as well as the corresponding value of is determined by the relation between the coupling energies and Assuming that is fixed, one observes that takes on greater values and the optimal temperature becomes higher as increases. On the contrary, enhancement of at fixed leads to a significant decrease of and shifts to a lower value. Thus the electron interactions with the thermal baths suppress values. Molecular vibrations may affect thermoelectric efficiency of the considered nanoscale systems in a similar way, as discussed in several works (see e.g. Refs. 13 (); 21 (); 30 ()).

The character of temperature dependence of displayed in the Fig. 4 indicates that while the transition from coherent and elastic tunneling to the dissipative transport significantly reduces values, the general character of temperature dependence of the figure of merit remains unchanged. At low temperatures, erosion of the sharp step in the Fermi distribution functions for the electrodes occurring at creates better opportunities for the electron tunneling through the system. However, when the temperature exceeds a certain value, the same process starts to hinder electron transport. Also, at sufficiently strong electron-phonon interactions, the peak value of the electron transmission decreases bringing further reduction of the thermoelectric efficiency.

Although considerable efforts are applied to reach understanding of combined effects of electron-electron and electron-phonon interactions on the thermoelectric transport, this subject is not fully investigated so far. Now, we reconsider the above results taking into account previously disregarded Coulomb interactions between electrons on the bridge of a thermoelectric junction. Then the expression (15) for the electron Green’s function may be employed to compute the scattering probabilities and, ultimately, the electron transmission and measurable characteristics of thermoelectric transport. In further analysis we assume that the linker (molecule/quantum dot) is weakly coupled to the electrodes so that the charging energy significantly exceeds the coupling parameter The results of these computations are displayed in the Figs. 5,6.

As shown in the Fig. 5, the dependence of of undergoes significant changes as the temperature increases. At low temperatures, exhibits two pairs of peaks of unequal height situated near and respectively. At higher temperatures two peaks making a pair cling together so that each pair is transformed to a single peak. At sufficiently high temperatures, these peaks become nearly equal in height, and their tops are shifted farther away from each other. The curves displayed in the left panel of the Fig. 5 are plotted assuming that temperature is noticeably lower than the temperature providing the maximum value of as shown in the Fig. 4. We cannot explicitly compare the results represented in these figures because the curves shown in the Fig. 4 are plotted disregarding electron-electron interactions. However, we may conjecture that further increase of temperature accompanied by intensification of scattering processes will bring furthermost rise of peaks as well as it happens in the case when one neglects electron-electron interactions. Also, we may expect that the increase in the peaks heights would be replaced by their reduction as the temperature would exceed a certain value. An explicit effect of electron-phonon interactions on the figure of merit is shown in the right panel of the Fig. 5. Again, one may observe the suppression of originating from these interactions.

Further illustration of the influence of thermal phonons on the figure of merit is presented in the Fig. 6. The curves shown in the left panel of this figure are plotted at a moderately high temperature when the adjacent peaks are already merged, so that exhibit only two maxima. Omitting electron-phonon interactions, one observes a significant difference in the peaks heights. This difference originates from the characteristic features of electron density of states on the bridge level which are manifested in the characteristics of coherent electron transport. As the electrons interaction with thermal phonons strengthens, the peaks heights become leveled. At fixed temperature, maximum value of is determined by the relation between the charging energy and coupling strengths and We remark that the suppression of due to electron-phonon interactions may be replaced by its promotion which occurs due to a combined effect of electron-electron and electron-phonon interactions 33 (). However, this increase of thermoelectric efficiency is expected to appear when electron-phonon interactions and Coulomb repulsion between electrons are comparable in strength. These conditions are different from those considered in the present work. Disregarding for a while the effect of phonons, one observes that takes on greater values within the limits of low and high values of charging energy, and it drops at intermediate values of This behavior was previously described and explained within the sequential hopping approximation for the electron transmission through a thermoelectric junction 21 (). Electron-phonon interactions promote washing out of these features.

## Iv iv. Conclusion

In conclusion, we remark that thorough studies of thermoelectric properties of nanoscale systems taking into account both electron and phonon transport as well as diverse effects arising due to electron-electron and electron-phonon interactions are not completed so far. In several earlier works this theoretical research was carried out employing single-particle scattering approach pioneered by Landauer in the context of charge transport in mesoscopic and nanoscale systems. These ideas were generalized to phonon transport through nanoscale junctions 30 (); 32 (); 42 (); 52 (); 53 (). Within this approach, transport characteristics of a considered nanoscale system are expressed in terms of electron and phonon transmission functions. The latter were computed employing several methods, including some based on scattering matrices formalism 42 (); 54 (). Later, these methods were mostly abandoned in favor of more advanced formalisms such as NEGF and/or various modifications of quantum rate equations. However, potential usefulness of the approaches based on scattering theory is not exhausted so far.

These approaches have an advantage of being computationally simple and less time and effort consuming than advanced formalisms. At the same time, their shortcomings could be largely removed by incorporating some NEGF based results into a computational scheme. In the present work we suggest such approach, and we employ it to theoretically analyze some effects of electron-phonon interactions on the efficiency of nanoscale thermoelectric junctions. Presently, various manifestations of electron-phonon interactions in thermoelectric transport characteristics of nanoscale molecular junctions are already explored, and the research is still going on. However, the research efforts were and still are mostly concentrated on the effects arising from vibrational modes on the molecules linking the electrodes. Less attention was paid to the influence of thermal phonons associated with random nuclear motions in the ambience. Here, we focus on the analysis of thermal phonons on the electron transport. We show that direct interaction of electrons with thermal phonons assuming that these phonons are assembled in two baths associated with the electrodes. may significantly affect thermoelectric efficiency of molecular junctions and similar nanoscale systems.

Specifically, we show that electron-phonon interactions assist the increase of the scattering probabilities thus destroying the coherence of electron transport and promoting energy dissipation. When the electron-phonon coupling becomes sufficiently strong, this brings a significant suppression of both thermopower and thermoelectric figure of merit thus worsening thermoelectric efficiency of a considered system. This effect is illustrated in the Figs. 3,4. We remark that and thermal energy appear as cofactors in the expression for the scattering probability (see Eq. (29)), so they affect it in a similar way. When either or increases, this results in strengthening of dephasing in the electron transport. However, entire effects of these two parameters on the thermoelectric properties of considered systems are unidentical. While the strengthening of electron-phonon interactions always leads to reduction of the rise of temperature can promote the figure of merit increase provided that temperature does not exceed a certain value. This may be explained by the fact that besides affecting the intensity of scattering, the temperature influences distributions of electrons in the electrodes, and it may either assist or hinder their transport through the system. Also, we analyzed the combined effect of electron-electron and electron-phonon interactions on thermoelectric properties. Obtained results agree with those reported in the earlier works 21 (); 33 (). In particular, it was confirmed that exhibits a minimum at a certain value of the charging energy which becomes less distinct at stronger values of

The suggested computational scheme may be generalized to include vibrational modes. For this purpose, one may add to the adopted model simulating a thermoelectric junction an extra reservoir representing vibrons. Also, one may mimic the molecular bridge by a set of energy levels thus opening the way to studies of interference effects. For a realistic molecular junction, relevant energies may be computed using either density functional theory or other method of electronic structure calculations. Finally, the proposed scheme may appear helpful in studies of thermoelectric transport beyond linear regime in temperature. On these grounds, we believe that presented method and results could help to reach better understanding of some important aspects of thermal transport in molecular junctions and similar nanoscale systems.

Acknowledgments: This work was supported by NSF-DMR-PREM 0934195 and NSF-EPS-1010094. The author thank G. M. Zimbovsky for help with the manuscript.

## References

- (1) M. Paulsson and S. Datta, Phys. Rev. B 67, 241403 (2003).
- (2) J. Koch, F. von Oppen, Y. Oreg, and E. Sela, Phys. Rev. B 70, 195107 (2004).
- (3) X. Zheng, W. Zheng, Y. Wei, Z. Zeng, and J. Wang, J. Chem. Phys. 121, 8537 (2004).
- (4) D. Segal, Phys. Rev. B 72, 165426 (2005).
- (5) P. Reddy, S.-Y. Jang, R. A. Segalman, A. Majumdar, Science 315, 1568 (2007).
- (6) K. Baheti, J. A. Malen. P. Doak, P. Reddy, S.-Y. Jang,T. D. Tilley, A. Majumdar, and R. A. Segalman, Nano Lett. 8 715 (2008).
- (7) F. Pauli, J. K. Viljas, and J. C. Cuevas, Phys. Rev. B 78, 035315 (2008). M. Galperin, A. Nitzan, and M. A. Ratner, Mol. Phys. 106, 394 (2008).
- (8) M. Galperin, A. Nitzan, and M. A. Ratner, Mol. Phys. 106, 397 (2008).
- (9) C. M Finch, V. M Garcia-Suarez, and C. J Lambert, Phys. Rev. B 79, 033405 (2009).
- (10) Y. S. Liu, Y. R. Chen, Y. C. Chen, ACS Nano, 3, 3497 (2009).
- (11) J. P. Bergfield, M. A. Solis, and C. A. Stafford, ACS Nano 4, 5314 (2010).
- (12) E. Pop, Nano Research 3, 146 (2010).
- (13) S. Y. Quek, H. J. Choi, S. G. Louie, and J. B. Neaton, ACS Nano 5, 551 (2010).
- (14) D. Nozaki, H. Sevincli, W. Li, R. Gutierrez, and G. Cuniberti, Phys. Rev. B 81, 235406 (2010).
- (15) M. Burkle, L. A. Zotti, J. K Viljas, D. Vonlanthen, A. Mishchenko, T. Wandlowski, M. Mayor, G. Schon, and F. Pauly, Phys. Rev. B 86, 115304 (2012).
- (16) B. Kubala and J. Konig, Phys. Rev. B 73,195316 (2006).
- (17) X. Zianni, Phys. Rev. B 75, 045344 (2007).
- (18) J. P. Heremans, C. M. Thrush, and D. T. Morelli, Phys. Rev. B 70, 115334 (2004).
- (19) M Krawiec and K. L. Wysokinski, Phys. Rev. B 73, 075307 (2006).
- (20) R. Swirkowicz, M. Wierzbicki, and J. Barnas, Phys. Rev. B 80, 195403 (2009).
- (21) P. Murphy, S. Mukerjee, and J. Moore, Phys. Rev. B 78, 161406(R) (2008).
- (22) J. Ren, J. X. Zhu, J. E. Gubernatis, C. Wang, B. Li, Phys. Rev. B 85, 155443 (2012).
- (23) B. Kubala, J. Konig, and J. Pecola, Phys. Rev. Lett. 100, 066801 (2008).
- (24) M. Yoshida, L. N. Olivera, Physica B 404, 3312 (2009).
- (25) M. Wierzbicki and R. Swirkowicz, Phys. Rev. B 84, 075410 (2011).
- (26) P. Trocha and J. Barnas, Phys. Rev. B 85, 085408 (2012).
- (27) R. Hartle, M. Butzin, O. Rubio-Pons, M. Thoss, Phys. Rev. Lett. 107, 046802 (2011).
- (28) P. Trocha and J. Barnas, Phys. Rev. B 76, 165432 (2007).
- (29) O. Karlstrom, H. Linke, G. Karlstrom, and A. Wacker, Phys. Rev. B 84, 113415 (2011).
- (30) M. Galperin, M. A. Ratner, and A. Nitzan, J. Phys.: Condens. Matter 19, 103201 (2007).
- (31) N. Sergueev, D. Roubtsov, and H. Guo, Phys. Rev. Lett. 95, 146803 (2005).
- (32) J. T. Lu and J.-S. Wang, Phys. Rev. B 76, 165418 (2007).
- (33) O. Entin-Wohlman, Y. Imry, and A. Aharony, Phys. Rev. B 82 115314 (2010).
- (34) D. M.-T. Kuo and Y.-C.Chang, Phys. Rev. B 81, 205321 (2010).
- (35) N. Sergueev, S. Shin, M. Kaviany, and B. Dunietz, Phys. Rev. B 83, 195415 (2011).
- (36) T. Markussen, A.-P. Jauho, and M. Brandbyge, Phys. Rev. B 79, 035415 (2009).
- (37) M. Leijnse, M. R. Wegewijs, and K. Flensberg, Phys. Rev. B 82, 045412 (2010).
- (38) Y. Dubi and M. Di Ventra, Nano Lett. 9, 97 (2009).
- (39) M. Buttiker, Phys. Rev. B 33, 3020 (1986).
- (40) Z. Wang, J. A. Carter, A. Lagutchev, Y. K. Koh, and N. H. Seong, Science 317, 787 (2007).
- (41) N. Mingo, Phys. Rev. B 74, 125402 (2006).
- (42) J. L. D’Amato and H. M. Pastawski, Phys. Rev. B 41, 7411 (1990).
- (43) X.-Q. Li and Y.-J. Yan, J. Chem. Phys. 115, 4169 (2001).
- (44) N. A. Zimbovskaya and M. M. Kuklja, J. Chem. Phys. 131, 114703 (2009).
- (45) Y. Meir, N. S. Wingreen, and P. A. Lee, Phys. Rev. Lett. 66, 3048 (1991).
- (46) M. Wierzbicki and R. Swirkowicz, Phys. Rev. B 82, 165334 (2010).
- (47) N. A. Zimbovskaya, J. Chem. Phys. 123, 114708 (2005).
- (48) C. D. Mahan, Many-Particle Physics (Plenum, New York, 2000).
- (49) Y. Dubi and M. Di Ventra, Rev. Mod. Phys. 83, 131 (2011).
- (50) Y.-C. Chang and D.-M. T. Kuo, Phys. Rev. B 77, 245412 (2008).
- (51) G. D. Sofo and J. O. Mahan, PNAS 93, 7436 (1996).
- (52) L. G. C. Rego, G. Kirczenow, Phys. Rev. Lett. 81, 232 (1998).
- (53) D. Segal, A. Nitzan, and P. Hanggi, J. Chem. Phys. 119, 6840 (2003).
- (54) M. Di Ventra, N. D. Lang, Phys. Rev. B 65, 045402 (2001).