# Origin of thermoelectric response fluctuations in single-molecule junctions

###### Abstract

The thermoelectric response of molecular junctions exhibits large fluctuations, as observed in recent experiments [e.g. Malen J. A. et al., Nano Lett. 10, 3406 (2009)]. These were attributed to fluctuations in the energy alignment between the highest occupied molecular orbital (HOMO) and the Fermi level at the electrodes. By analyzing these fluctuations assuming resonant transport through the HOMO level, we demonstrate that fluctuations in the HOMO level alone cannot account for the observed fluctuations in the thermopower, and that the thermo-voltage distributions obtained using the most common method, the Non-equilibrium Green’s function method, are qualitatively different than those observed experimentally. We argue that this inconsistency between the theory and experiment is due to the level broadening, which is inherently built-in to the method, and smears out any variations of the transmission on energy scales smaller than the level broadening. We show that although this smmearing only weakly affects the transmission, it has a large effect on the calculated thermopower. Using the theory of open quantum systems we account for both the magnitude of the variations and the qualitative form of the distributions, and show that they arise not only from variations in the HOMO-Fermi level offset, but also from variations of the local density of states at the contact point between the molecule and the electrode.

## I Introduction

Improving the thermoelectric (TE) energy conversion efficiency relies on increasing the thermoelectric response (Seebeck coefficient) and reducing the thermal conductance of the device Nitzan (2007); Dubi and Di Ventra (2011). Metal-single molecule-metal junctions (’molecular junctions’) seem promising in both aspects: their thermal conductance is small due to the mismatch between the vibrations of the electrodes and the molecule, and their TE response should be large due to the well-defined resonant structure of the electron transport (via the molecular HOMO or LUMO levels) Malen et al. (2010). This observation, along with the notion that TE response can shed light on transport mechanisms in molecular junctions Paulsson and Datta (2003); Koch et al. (2004); Segal (2005), have initiated in large number of experimental Reddy et al. (2007); Baheti et al. (2008); Malen et al. (2009a, b); Tan et al. (2011); Yee et al. (2011); Widawsky et al. (2012) and theoretical studies Murphy et al. (2008); Dubi and Di Ventra (2008); Liu and Chen (2009); Bergfield and Stafford (2009); Liu et al. (2009); Wang et al. (2010); Leijnse et al. (2010); Sergueev et al. (2011); Liu et al. (2011a, b); Quek et al. (2011); Stadler and Markussen (2011); Nikolic et al. (2012); Bürkle et al. (2012) on TE effects in molecular junctions.

Most notable are a series of impressive experiments in which TE conversion using a single-molecule junction was demonstrated Reddy et al. (2007); Baheti et al. (2008); Malen et al. (2009a, b, 2010). In these experiments, a single molecule (usually Benzene rings with various end groups) is trapped between an Au substrate and an Au scanning tunneling microscope (STM) tip, which are held at a constant temperature difference . A voltage bias (here called ’thermo-voltage’) is then applied between the tip and substrate to reach a state of zero current flowing through the molecule. The Seebeck coefficient is defined as (minus) the slope of as a function of in the linear response regime (i.e. ). In the experiments the value of is strongly fluctuating, and repeating the experiment many times results in a broad distribution of , as can be seen in Fig. 3 of Ref. Malen et al. (2009a) (and inset of Fig. 2(b) here). Using the typical (where the voltage distribution displays a maximum), the authors of Ref.Reddy et al. (2007); Malen et al. (2009a) find typical Seebeck coefficients of the order of V/K for single molecule junctions.

The observed variations in the thermoelectric response do not seem to be simple (Gaussian) noise, and the distributions of the thermo-voltage may have a well-defined double-peak structure Reddy et al. (2007); Malen et al. (2009a). The authors of Ref. Malen et al., 2009a concluded that the fluctuations (also reported in thermoelectric measurements of metal-fullerene-metal junctions Yee et al. (2011)) are caused by variations of the position of the HOMO level, which are translated into variations in the energy offset between the HOMO level and the Fermi energy of the contacts, , as the junction reconstructs. Assuming that the transport is dominated by a single resonant level and using the Landauer formalism (or Green’s function method) for thermoelectric transport (see, e.g. Refs. Di Ventra, 2008; Dubi and Di Ventra, 2011 ) they quantified the variations in and found them to be , of the order of the average offset itself.

In this paper we argue that besides the variations in the offset energy between the HOMO level and the Fermi energy, there is an additional contribution to the fluctuations in the thermopower, which is the strong variations of the LDOS (as a function of energy) at the point of contact between the molecule and the electrodes. We use a toy model to verify this effect. The main results are (I) Within the Landauer or Green’s function formalism for thermoelectric transport, any variations in the transmission function which are smaller than the self-energy (or level broadening) are suppressed. The result is a smooth transmission function and a smooth thermopower (as a function of energy).(II) The thermopower is extremely sensitive to small fluctuations of in the transmission function. As a result, the suppression of fluctuations in the transmission function, inherent to the Green’s function method, strongly affects on the calculated thermopower, and results in an over-estimation for the fluctuations in the HOMO level-Fermi level energy offset fluctuations. (III) The thermo-voltage histogram obtained using the Green’s function method is qualitatively different from the histograms obtained experimentally, and (IV) all of the above draw-backs of the Green’s function formalism can be overcome by using the open quantum system approach to calculating the thermo-electric response in molecular junctions.

The paper is organized as follows. In Sec. II we derive the estimation for the fluctuations in the HOMO-Fermi level energy offset from the experimental values of thermopower for typical molecular junctions, using a simple description of the molecular junction as a resonant level and using the Landauer formula. In Sec. III we present the open quantum systems approach, and using a simple model for the molecular junction we show that this approach for calculation of thermopower yields results which are in qualitative agreement with the experimental results. In Sec. IV we show that, for the same model, the Green’s function approach gives very different results, due to the inherent suppression of fluctuations built in to the method (due to the presence of a self-energy of level broadening). We demonstrate that the thermopower is very sensitive to fluctuations, and thus the suppression of (small) fluctuations in the transmission strongly affects the calculation of thermo-power. Sec. V is devoted to summary and discussion.

## Ii Analysis of thermopower fluctuations using the Landauer formula for a resonant level

As a starting point, we consider the data presented in Fig. 3 of Ref. Malen et al., 2009a (also inset of Fig. 2 here), where the thermo-voltage histograms of a 2’,5’-dimethyl-4,4”-tribenzenedithiol (DMTBDT) molecular junction are presented. We observe several notable features: (i) Noise is present in the measurement: even at there is a finite thermo-voltage distribution, with an average signal offset of a few tens of microvolts, (ii) the typical form of the thermo-voltage distribution seems Gaussian, (iii) a double-peak structure of the distribution emerges, indicating two dominant values of the Seebeck coefficient, V/K and V/K.

In order to estimate the variation in the molecular HOMO required to obtain these values, we use the standard Landauer formalism . Within this formalism, the conductance is given by , and the Seebeck coefficient can be approximated as

(1) |

where is the electron energy, is the electron transmission function, is the proton charge, the Boltzmann constant and is the average temperature.

We proceed by assuming that the transport through molecule is characterized by a transmission through a resonant level. To justify this, as well as the toy model considered in the next section, we briefly review the literature on theoretical calculations of thermopower in molecular junctions. The most common theoretical tool to calculate the transmission (and from it the thermopower) is the combination of non-equilibrium Green’s function (NEGF) and density-functional theory (DFT) Di Ventra (2008), where the transmission function is evaluated using the Green’s function as obtained from DFT. Within this approach, the transmission function is calculated as a function of gate voltage, and from it the thermopower is calculated, also as a function of gate voltage (see, e.g. Refs. Liu and Chen, 2009; Stadler and Markussen, 2011; Tan et al., 2011; Quek et al., 2011; Balachandran et al., 2012). While in transport experiments gating of a molecule can be achieved Lee et al. (2003); Park et al. (2000); Li et al. (2006), such a dependence of the thermopower on the gate voltage was never studied experimentally, and experiments always carry a statistical nature as discussed above. However, to our knowledge the statistical nature of thermopower in molecular junctions has never been addressed (although it has been addressed in the context of transmission Reuter et al. (2012) or thermopower of atomic wires Pauly et al. (2011)). Nevertheless, to justify our model we observe that in the theoretical calculations, the transmission function seems always to resemble a resonant level, at least close to the HOMO or LUMO levels. Liu and Chen (2009); Stadler and Markussen (2011); Tan et al. (2011); Quek et al. (2011); Balachandran et al. (2012). Thus, considering a resonant level model is a common description even of realistic molecules with electron interactions.

For a given resonant molecular level with energy and contact Fermi energy , the typical form of the transmission function is a resonant Lorentzian, , where is the energy offset and is the contact-induced effective level broadening. The resulting Seebeck coefficient is of the form

(2) |

Noting the double-peak structure of the voltage distribution, we thus consider that there is a typical shift in the HOMO level . Taking the conductance value to be , one needs to solve simultaneously These equations yield (at ambient temperature) eV, eV and eV, such that , similar to the value obtained in Ref. Malen et al., 2009a. This means that the molecular level alignment with the contact Fermi level is changing by electron-volts at every reconstruction of the junction. This large variation (of the order of itself) was attributed to variations in the junction contact geometry and intermolecular interactions. However, a recent study of transition voltage spectroscopy in molecular junctions (see Ref. Beebe et al., 2006 for description of the method) demonstrated that the variation of are of the order of eV, Guo et al. (2011) much smaller than the variation required to generate the value above. This discrepancy suggests that there is an additional factor contributing to the variations in thermopower.

The discrepancy describes above hints that there is a flaw in the analysis of the thermopower fluctuations using the Landauer formula with the resonant level model. As we will show in the following sections, the flaw is that the Landauer formula inherently suppresses any variations in the transmission function (as a function of molecular level energy ) which are on an energy scale smaller than the level broadening (or electron self-energy), and since the thermopower is very sensitive to variations in the transmission, this suppression strongly affects the thermopower. Thus, fluctuations in sample a thermopower function which is "too smooth", resulting in the over-estimate of .

## Iii Open quantum system approach to thermopower

In order to identify the origin of variations in the thermopower, we study a toy model of a molecular junction. We use the method of open quantum systems (OQS) (which was described in detail in previous publications Pershin et al. (2008); Dubi and Di Ventra (2008)) to show that the thermopower exhibits variations as a function of which are on the same energy scale as the variations in the LDOS at the molecule-lead point of contact, and show that the thermo-voltage histogram resulting in fluctuations in is qualitatively similar to the experimental ones. In the next section we will compare the result ontained from the OQS method to those obtained by using NEGF method.

Our toy model consists of two finite electrodes with a molecular bridge between them (upper panel of Fig. 1), attached to reservoirs with different temperatures and for the left and right edges respectively. The molecule is described by a simple chain with four atomic orbitals. The OQS method enables one to calculate the electronic density in this non-equilibrium situation (a finite temperature difference between the electrodes). From the electronic density, one can calculate (via the Poisson equation) the voltage difference between the electrodes as a function of the temperature difference . Repeating the calculation for different values of we then find a curve , and the linear slope is the thermopower (or Seebeck coefficient) .

Essentially, OQS theory is a mapping of the system Hamiltonian (which in general includes all electronic degrees of freedom, electron interactions, junction geometry etc.) onto a master equation for the single-particle density matrix of the Lindblad form Pershin et al. (2008),

(3) |

We consider a non-interacting tight-binding Hamiltonian of the form , where

(4) |

are the tight-binding Hamiltonians of the left and right leads respectively ( is the hopping integral, which is taken as the unit energy, eV from here on), and

(5) |

is the Hamiltonian for the wire, which includes the usual hopping integral and an energy which can be tuned (for instance using a gate electrode experimentally). The energy is measured with respect to the electrodes’ Fermi energy which is set as the zero energy ().

The coupling between the electrodes and the molecular wire is described by

(6) |

describes the coupling between the left (right) lead to the wire, with being the creation operator for an electron at the point of contact between the left (right) lead and the wire, and () destroys an electron at the left-most (right-most) sites of the wire (we take here ). The external environment(s) are accounted for in Eq. 3 by the second term, which has the form

(7) |

where are the Lindblad -operators which encode the properties of the environment, i.e. its temperatures and position (i.e. left or right electrode). An appropriate form for the -operators is Pershin et al. (2008)

(8) |

where are the Fermi distributions of the left and right leads (with the corresponding temperature), the chemical potential, and are the overlap integrals between the and states on the left (L) and right (R) edges of the electrodes. This form for the -operators guarantees that at equilibrium (i.e. ) the diagonal elements of the density matrix are described by a Fermi function.

We note that the Lindblad equation assumes a Markov approximation for the environment, which is reasonable since we are considering a system at room temperature (where quantum memory effect of the environment should not be important), and since we are interested in the steady state (and not the dynamics). The numerical calculation is performed in the following way: (i) From the tight-binding Hamiltonian the single-particle states and energies are calculated, and the -operators are constructed according to Eq. (8) for different electrode temperatures and . (ii) The -operators are inserted to Eq. (3) and Eq. (7), which are then solved in the steady-state (i.e. for ). (iii) From the diagonal elements steady-state solution for and the wave functions, the electron density is calculated. (iv) By solving the Poisson equation, the voltage at the center of the electrodes is calculated. The voltage difference between the electrodes is the thermo-voltage, as it is induced by the temperature difference. the slope of the thermo-voltage as a function of the temperature difference is the thermopower.

In Fig. 1 we plot the the Seebeck coefficient (solid blue line) and the electron occupation of the molecular chain (dashed orange line) of a model molecular junction as a function of a gate potential applied to the junction, i.e. as a function of the molecular energy level with respect to the Fermi level (which is set as the zero energy). The junction is composed of a series of four atomic orbitals connected to square two dimensional electrodes (upper panel). Numerical parameters are: electrode size , tight-binding hopping integral eV (corresponding to a band-width of eV), electrode-molecule coupling is eV, Temperature is room temperature, and the electrodes are kept at half filling.

The first thing to be noted is that whenever a molecular orbital crosses the Fermi energy, the molecular occupation changes by one, and the thermopower changes sign. This is in accord with the known results from the Landauer formula for thermopower, and is due to the change from electron-dominated to hole-dominated transport every time the Fermi level is crossed (which also corresponds to a transmission resonance). The second important feature is that, in contrast to the result expected from the regular Landauer formula Paulsson and Datta (2003); Dubi and Di Ventra (2011), the Seebeck coefficient exhibits strong variations with gate voltage.

Before we proceed to discussing the origin of these variations, it is useful to compare the properties of the Seebeck coefficient as obtained from the OQS method to those known from the Landauer formalism. In Fig. 1(b) the Seebeck coefficient as a function of gate voltage (same as in Fig. 1(a)) is plotted for different values of the coupling between the electrode and the molecule, and meV. As expected from the Landauer formalism Quek et al. (2011), there is little effect to the coupling on the magnitude of the Seebeck coefficient. In Fig. 1(c) the Seebeck coefficient is plotted as a function of temperature for various molecular chain lengths. An inhomogeneous temperature dependence is found, originating from a crossover from coherent to incoherent transport, again in agreement with results obtained from the Landauer formalism (e.g. Segal (2005)). These results demonstrate that the main physical features of the thermopower which are present in the Landauer formalism, also appear within the framework of the OQS theory.

We now turn to calculating the distribution of the thermo-voltage across the molecular junction. To obtain the distribution, we calculate the temperature-difference induced voltage , taking the gate voltage of the molecule (i.e. the HOMO-Fermi energy offset) to be a random variable, normally distributed around eV with a width eV. While these values are rather tdifferent then experimental values (probably is bigger in experiments) we point that we are aiming at qualitative similarity to experiment, to point the origin of the variations, rather than to analyze realistic junctions. To mimic the experiment, we also apply a small variation (two degrees Kelvin) to the electrode temperatures (note that in the experiments, a finite thermo-voltage distibution was observed even at ,see inset of Fig. 2(b), indicating the existence of a small temperature difference, probably due to Johnson noise Malen (). This does not affect our results or conclusions).

In Fig. 2(a) the distribution of thermo-voltage is plotted for temperature differences K. The distributions qualitatively resemble the experimental distributions, exhibiting a broad double-peak structure. To understand the origin of the distribution shapes, in the right panel of Fig. 2(a) the Seebeck coefficient is again plotted (solid line). Due to the strong sensitivity of on , even a relatively small variation in (of the order of eV) can include several maxima of , giving rise to the different peaks in the distributions.

To understand the origin of the sensitivity of , in the right panel of Fig. 2(a) we plot the local density of states at the point of contact between the electrode and the molecular chain (dashed line). One can see that the variations of and the LDOS vary on the same energy scale. Note that these variations are a surface effect, and are not due to the electrode level spacing (the average level spacing is ). We thus conclude that it is the variations in the LDOS (which are on an energy scale much smaller than the energy difference between the molecular levels) that give rise to the sensitivity of (we note that the total DOS of the electrodes is a much smoother function with no observed variations on these energy scales). We point that the LDOS variations are not a finite-size effect, but rather are a surface effect. Calculating the LDOS to systems as large as (where the level spacing is eV), we found that local variations in the LDOS persist with roughly the same energy scale of eV (although their magnitude somewhat changes).

On the other hand, when considering the thermopower as obtained by the Landauer formula, the local variations in the LDOS are not taken into account, or rather they are smeared by the electron self energy which is reflected through the level broadening in Eq. (2). The resulting is a much smoother function. As a first demonstration of this effect, here we use the formula for transport and thermopower through a resonant level, Eq. 1. The positions of the levels and the level broadenings are obtained by fitting the data of Fig. 1(a). In the right panel of Fig. 2(b) is plotted using Eq. (2) and the parameters of the junction (i.e. resonances and widths fitted to Fig. 1(a)). For comparison, obtained from OQS theory is plotted as a bright line in the background. The resulting distributions are plotted in Fig. 2(b) for different . The distributions obtained using Eq. 2 are quite different from the experimental results (inset of Fig. 2(b)), both in terms of shape and in the lack of the double-peak structure.

## Iv Green’s function analysis of the toy model

Let us now examine the same model within the Green’s function formalism. In this formalism, and for our non-interacting toy model, the transmission function is given by Meir and Wingreen (1992); Di Ventra (2008)

(9) |

where are the retarded and advanced Green’s functions, and represent the level broadening due to the electrodes, typically a few hundreds meVs Bergfield and Stafford (2009); Bürkle et al. (2012). In the basis of atomic orbitals (as the Hamiltonian is written) are diagonal matrices, with if is in the left or right edges of the electrodes, and zero otherwise. We take eV in the numerical example below. is the self-energy, which in the non-interacting case is only due to the electrodes, and hence .

Once the transmission is calculated, the thermopower can be calculated directly using Eq. (1). In Fig. 3(a) the transmission function and the its logarithmic derivative (proportional to the thermopower ) of the toy model are plotted (as a function of the gate energy). For comparison, the local density of states at the point of contact between the wire and the electrode is plotted in the upper inset (the dashed line in the inset is the DOS when broadened by eV). As seen, the oscillations in the LDOS are completely smeared in the transmission (plotted as a function of the gate energy). For comparison, the local density of states is plotted in the upper inset (the dashed line in the inset is the DOS when broadened by eV). As seen, the oscillations in the LDOS are completely smeared in the transmission function and hence in the thermopower. The lower inset shows the transmission function on a log scale, verifying that the variations are really smeared out. To demonstrate that this smearing has a large effect on the thermopower, we consider an artificial transmission function with eV. This transmission exhibits oscillations on an energy scale . However, they are chosen in such a way that they are hardly visible if the transmission is observed in its full scale, as seen from Fig. 3(b), where is plotted. In the inset the transmission is plotted on a log scale, and only there the oscillations are visible. On the other hand, the thermopower exhibits very strong oscillations of considerable size even if the oscillations are hardly observed in the transmission functions. This is due to the extreme sensitivity of to local variations of , reflecting its derivative structure.

## V Summary and discussion

In summary, we have analyzed the variations in the thermoelectric response of a metal-single molecule-metal junctions, based on detailed examination of experimental results. The experimental results, namely the width of the thermo-voltage variations and the shape of the thermo-voltage distributions, cannot be accounted for by only assuming variations in the misalignment between the molecular HOMO and the electrode Fermi energy. Using the theory of open quantum systems we qualitatively reproduce the experimental results, and show that they may originate in a combination of the level misalignment variations and variations of the local density of states at the point of contact between the molecule and the electrode, specifically the STM tip in typical molecular junction experiments.

To put it differently, we found that in order to explain the thermopower variations in molecular junctions the electronic transmission function cannot have a simple Lorentzian form, but rather should have a more complicated form which includes variations on a scale smaller than the HOMO-LUMO gap. Such variations may be induced by the variations in the LDOS at the tip-molecule point of contact, originating from the STM tip structure, impurities or trapped states, etc. The use of the Green’ function formula smears out all variations on energy scale smaller than the level broadening (typically a few hundred meVs). This smearing has a very small effect on the transmission and hence the conductance of the molecular junction, but have a strong effect on the thermopower. Thus, caution needs to be taken when using the NEGF method for calculation of themopower. We note a recent paper Evans and Voorhis (2009) arguing that the NEGF method has limitations even when calculating the conductance of a molecular junction (although the argument is different then presented in this paper).

The effects of STM tip structure on local transport measurements have been extensively discussed in the STM literature Tersoff (1990); Pelz (1991); Hofer et al. (2003); Passoni et al. (2009); Kwapinski and Jalochowski (2010), and were recognized to have an important contribution to the overall charge transport, via the tip LDOS. Our results imply that a similar importance of STM tip structure appears in thermo-electric measurements. Consequently, a direct extraction of relevant parameters such as the HOMO-Fermi level offset from thermo-electric measurements requires taking the electronic properties of the metal-molecule interface into account.

The author wishes to thank J. Malen for valuable comments on the manuscript.

## References

- Nitzan (2007) A. Nitzan, Science 317, 759 (2007).
- Dubi and Di Ventra (2011) Y. Dubi and M. Di Ventra, Reviews of Modern Physics 83, 131 (2011).
- Malen et al. (2010) J. Malen, S. Yee, A. Majumdar, and R. Segalman, Chemical Physics Letters 491, 109 (2010).
- Paulsson and Datta (2003) M. Paulsson and S. Datta, Phys. Rev. B 67, 241403 (2003).
- Koch et al. (2004) J. Koch, F. von Oppen, Y. Oreg, and E. Sela, Phys. Rev. B 70, 195107 (2004).
- Segal (2005) D. Segal, Phys. Rev. B 72, 165426 (2005).
- Reddy et al. (2007) P. Reddy, S. Jang, R. Segalman, and A. Majumdar, Science 315, 1568 (2007).
- Baheti et al. (2008) K. Baheti, J. Malen, P. Doak, P. Reddy, S. Jang, T. Tilley, A. Majumdar, and R. Segalman, Nano letters 8, 715 (2008).
- Malen et al. (2009a) J. Malen, P. Doak, K. Baheti, T. Tilley, A. Majumdar, and R. Segalman, Nano letters 9, 3406 (2009a).
- Malen et al. (2009b) J. Malen, P. Doak, K. Baheti, T. Tilley, R. Segalman, and A. Majumdar, Nano letters 9, 1164 (2009b).
- Tan et al. (2011) A. Tan, J. Balachandran, S. Sadat, V. Gavini, B. Dunietz, S. Jang, and P. Reddy, Journal of the American Chemical Society (2011).
- Yee et al. (2011) S. Yee, J. Malen, A. Majumdar, and R. Segalman, Nano letters (2011).
- Widawsky et al. (2012) J. R. Widawsky, P. Darancet, J. B. Neaton, and L. Venkataraman, Nano Letters 12, 354 (2012).
- Murphy et al. (2008) P. Murphy, S. Mukerjee, and J. Moore, Phys. Rev. B 78, 161406 (2008).
- Dubi and Di Ventra (2008) Y. Dubi and M. Di Ventra, Nano Letters 9, 97 (2008).
- Liu and Chen (2009) Y.-S. Liu and Y.-C. Chen, Phys. Rev. B 79, 193101 (2009).
- Bergfield and Stafford (2009) J. P. Bergfield and C. A. Stafford, Nano Letters 9, 3072 (2009).
- Liu et al. (2009) Y.-S. Liu, Y.-R. Chen, and Y.-C. Chen, ACS Nano 3, 3497 (2009).
- Wang et al. (2010) R.-Q. Wang, L. Sheng, R. Shen, B. Wang, and D. Y. Xing, Phys. Rev. Lett. 105, 057202 (2010).
- Leijnse et al. (2010) M. Leijnse, M. R. Wegewijs, and K. Flensberg, Phys. Rev. B 82, 045412 (2010).
- Sergueev et al. (2011) N. Sergueev, S. Shin, M. Kaviany, and B. Dunietz, Phys. Rev. B 83, 195415 (2011).
- Liu et al. (2011a) Y.-S. Liu, B. C. Hsu, and Y.-C. Chen, The Journal of Physical Chemistry C 115, 6111 (2011a).
- Liu et al. (2011b) Y.-S. Liu, H.-T. Yao, and Y.-C. Chen, The Journal of Physical Chemistry C 115, 14988 (2011b).
- Quek et al. (2011) S. Y. Quek, H. J. Choi, S. G. Louie, and J. B. Neaton, ACS Nano 5, 551 (2011).
- Stadler and Markussen (2011) R. Stadler and T. Markussen, The Journal of Chemical Physics 135, 154109 (2011).
- Nikolic et al. (2012) B. Nikolic, K. Saha, T. Markussen, and K. Thygesen, Journal of Computational Electronics 11, 78 (2012).
- Bürkle et al. (2012) M. Bürkle, L. A. Zotti, J. K. Viljas, D. Vonlanthen, A. Mishchenko, T. Wandlowski, M. Mayor, G. Schön, and F. Pauly, Phys. Rev. B 86, 115304 (2012).
- Di Ventra (2008) M. Di Ventra, Electrical Transport in Nanoscale Systems (Cambridge University Press, 2008).
- Balachandran et al. (2012) J. Balachandran, P. Reddy, B. D. Dunietz, and V. Gavini, The Journal of Physical Chemistry Letters 3, 1962 (2012).
- Lee et al. (2003) J.-O. Lee, G. Lientschnig, F. Wiertz, M. Struijk, R. A. J. Janssen, R. Egberink, D. N. Reinhoudt, P. Hadley, and C. Dekker, Nano Lett. 3, 113 (2003).
- Park et al. (2000) H. Park, J. Park, A. K. L. Lim, E. H. Anderson, A. P. Alivisatos, and P. L. McEuen, Nature 407, 57 (2000).
- Li et al. (2006) X. Li, B. Xu, X. Xiao, X. Yang, L. Zang, and N. Tao, Faraday Discuss. 131, 111 (2006).
- Reuter et al. (2012) M. G. Reuter, M. C. Hersam, T. Seideman, and M. A. Ratner, Nano Letters 12, 2243 (2012).
- Pauly et al. (2011) F. Pauly, J. K. Viljas, M. Bürkle, M. Dreher, P. Nielaba, and J. C. Cuevas, Phys. Rev. B 84, 195420 (2011).
- Beebe et al. (2006) J. M. Beebe, B. Kim, J. W. Gadzuk, C. Daniel Frisbie, and J. G. Kushmerick, Phys. Rev. Lett. 97, 026801 (2006).
- Guo et al. (2011) S. Guo, J. Hihath, I. Diez-Perez, and N. Tao, Journal of the American Chemical Society 133, 19189 (2011).
- Pershin et al. (2008) Y. Pershin, Y. Dubi, and M. Di Ventra, Physical Review B 78, 054302 (2008).
- (38) J. Malen, , private communications.
- Meir and Wingreen (1992) Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 (1992).
- Evans and Voorhis (2009) J. S. Evans and T. V. Voorhis, Nano Letters 9, 2671 (2009).
- Tersoff (1990) J. Tersoff, Physical Review B 41, 1235 (1990).
- Pelz (1991) J. Pelz, Physical Review B 43, 6746 (1991).
- Hofer et al. (2003) W. Hofer, A. Foster, and A. Shluger, Reviews of Modern Physics 75, 1287 (2003).
- Passoni et al. (2009) M. Passoni, F. Donati, A. Bassi, C. Casari, and C. Bottani, Physical Review B 79, 045404 (2009).
- Kwapinski and Jalochowski (2010) T. Kwapinski and M. Jalochowski, Surface Science 604, 1752 (2010).