From Chemistry to Functionality:Trends for the Length Dependence of the Thermopower in Molecular Junctions.

From Chemistry to Functionality:
Trends for the Length Dependence of the Thermopower in Molecular Junctions.

Falco Hüser Nano Science Center and Department of Chemistry
University of Copenhagen, 2100 København Ø, Denmark
   Gemma C. Solomon gsolomon@nano.ku.dk Nano Science Center and Department of Chemistry
University of Copenhagen, 2100 København Ø, Denmark
July 19, 2019
Abstract

We present a systematic ab-initio study of the length dependence of the thermopower in molecular junctions. The systems under consideration are small saturated and conjugated molecular chains of varying length attached to gold electrodes via a number of different binding groups. Different scenarios are observed: linearly increasing and decreasing thermopower as function of the chain length as well as positive and negative values for the contact thermopower. Also deviation from the linear behaviour is found. The trends can be explained by details of the transmission, in particular the presence, position and shape of resonances from gateway states. We find that these gateway states do not only determine the contact thermopower, but can also have a large influence on the length-dependence itself. This demonstrates that simple models for electron transport do not apply in general and that chemical trends are hard to predict. Furthermore, we discuss the limits of our approach based on Density Functional Theory and compare to more sophisticated methods like self-energy corrections and the GW theory.

The perfect material for potential application in efficient thermoelectrical devices should be small in size, stable and, most of all, highly tunable. In these respects, organic molecules seem very promising candidates, since chemistry offers a vast number of possibilities for designing their electrical and thermal properties.Mahan and Sofo (1996); Humphrey and Linke (2005); Wang, Segalman, and Majumdar (2006); Venkataraman et al. (2006a, 2007); Kiguchi et al. (2010); Capozzi et al. (2014); Strange et al. (2015) In particular, simple carbon-based molecular wires bridging two gold electrodes via selected binding groups can provide systematic insight into fundamental electronic and thermal transport mechanisms. Ho Choi, Kim, and Frisbie (2008); Solomon et al. (2008, 2009); Hines et al. (2010); Lee et al. (2013) Ideally, mapping out chemical trends would be a new and easy way for predicting a molecule’s function. However, only little is known about the interplay between chemical structure and thermoelectric features like conductance or thermopower. The large number of degrees of freedom in molecular junctions, such as the junction geometry, conformation and the chemical structure, makes it particularly challenging to establish general trends.

Typical binding groups for gold electrodes are thiolsXu and Tao (2003) and amines.Venkataraman et al. (2006b) More recently, the thermopower of single molecules linked via isocyanide groupsTan et al. (2011) as well as through direct covalent \ceAu-\ceC bondsCheng et al. (2011) has been measured. Amines exhibit a well-defined weak electronic coupling due to bonding of the \ceN lone pair to an undercoordinated \ceAu atom, allowing for stable measurements. Strong coupling is given by thiol binding groups, where a single sulfur atom can bind to up to three gold atoms.Kristensen et al. (2008) However, uncertainties in the contact geometry lead to large variations in the measured and calculated conductances. A well-defined contact geometry is achieved through the direct -type \ceAu-\ceC bond to an undercoordinated gold atom. This also allows for direct coupling of the electrodes into the carbon backbone of the molecule.

In the linear response regime, the thermopower (also called Seebeck coefficient), is defined as

(1)

where is the thermally induced voltage at the steady state zero current due to an applied temperature difference, , across the junction.

It has been widely stated that a positive thermopower corresponds to hole transport or tunnelling through the highest occupied molecular orbital (HOMO) and a negative thermopower to electron transport or tunnelling through the lowest unoccupied molecular orbital (LUMO). The question, however, is how these orbitals of the molecule in the junction compare to the situation of the molecule in the gas phase. In the presence of metallic electrodes, the molecular orbitals (MOs) become renormalized and broadened.Neaton, Hybertsen, and Louie (2006) Charge transfer can lead to partial occupation numbers. Due to hybridization of localized MOs with metallic states at the interface with the electrodes, so-called gateway states can appear within the HOMO-LUMO gap.Widawsky et al. (2012); Kim et al. (2014); Khoo et al. (2015)

Experimentally, it has been found that the thermopower increases linearly with the length of a molecular chain, .Reddy et al. (2007); Malen et al. (2009) This is supported in theory for the off-resonant tunnelling regime, where the transmission is exponentially surpressed around the Fermi energyMcConnell (1961); Karlström, Strange, and Solomon (2014)

(2)

with parameters and , which determine the decay factor and the contact resistance, respectively. Then, the thermopower becomesPauly, Viljas, and Cuevas (2008)

(3)

Here, denotes the derivative with respect to energy. Eq. (3) implies that the length dependence of the thermopower (the slope) is given by , whereas the contact thermopower, , only depends on .

Recent tight-binding calculations support the common assumption that is related to the molecular backbone and position of the MO energies, since it is independent of the contact coupling.Karlström, Strange, and Solomon (2014) , on the other hand, is directly associated with states on the binding groups and has a strong dependence on the coupling strength. Depending on the position of the Fermi level and the binding group energy, four different scenarios for the length dependence of the thermopower can in principal be found: positive or negative contact thermopower with increasing or decreasing values with length. A case, where the thermopower changes sign with increasing molecular length has been measured recently for the first time for oxidized oligothiophenes.Dell et al. (2015)

I Method

The transmission, , is calculated within the Landauer-Büttiker formalism for coherent transportMeir and Wingreen (1992) using standard nonequilibrium Green’s Functions techniques based on Density Functional Theory (NEGF-DFT).Taylor, Guo, and Wang (2001); Xue, Datta, and Ratner (2002); Brandbyge et al. (2002); Thygesen, Bollinger, and Jacobsen (2003) Within this framework, the thermopower is given as

(4)

where is the elementary charge, the absolute temperature and the Fermi distribution. In the limit of small temperature gradients, , Eq. (4) reduces toButcher (1990); Paulsson and Datta (2003); Segal (2005); Malen et al. (2010)

(5)

where is the Boltzman constant. All thermopower values presented in this work have been calculated from Eq. (4) at room temperature (). Eq. (5) illustrates that the thermopower is proportional to the negative logarithmic derivative of the transmission at the Fermi energy of the electrodes, .

All calculations have been carried out with the electronic structure code GPAWEnkovaara et al. (2010) using the generalized gradient PBE exchange-correlation functional.Perdew, Burke, and Ernzerhof (1996) A double-zeta polarized basis set of diffuse basis functions with a confinement-energy shift of was employed.Strange et al. (2011) The junctions were modelled in supercells with 7-8 layers of gold containing 16 atoms each arranged in the fcc structure. A -point sampling was applied in the directions of the electrode surface and the transmission was averaged over all points.Thygesen and Jacobsen (2005); Strange et al. (2011) The coupling to the bulk electrodes is obtained using standard techniques.Strange et al. (2008, 2011)

The junction geometries are sketched in Fig. 1: Nitrogen and carbon atoms were chosen to bind to a single gold adatom at the tip of a small pyramid, whereas sulfur atoms were placed in a fcc-hollow site of a flat gold surface. For each binding group, alkanes (saturated carbon chains), alkenes (conjugated carbon chains with alternating single and double bonds) and alkynes (conjugated carbon chains with alternating single and triple bonds) with 2, 4, 6 and 8 carbon units were modelled. Geometry optimization was performed for all molecules (including the tip atoms in the pyramid configuration) on ideal gold surfaces.

Figure 1: Left side: Chemical structure of (i) alkane, (ii) alkene and (iii) alkyne chains. \ceX = \ceS, \ceNH2, \ceNC, - (direct coupling). is the number of \ce(CH2)2, \ce(CH)2 and \ce(C)2 units, respectively. Right side: Junction geometries for (a) thiol-, (b) isocyanide- and (c) amine-end groups and (d) direct coupling.

Ii Results

ii.1 Alkane chains

Figure 2: Calculated thermopower for alkane chains with amine-, isocyanide- and sulfur-binding groups as well as direct coupling in the junction geometries of Fig. 1. Straight lines are linear fits to , where is the number of \ce(CH2)2 units. All fitted values are listed in Table 1. The first data point is not considered in the fitting for direct coupling. The dotted horizontal line serves as guide for the eye at .

The thermopower for the alkane chains is shown in Fig. 2. It increases roughly linearly with length for all binding groups with a rather small prefactor of around per methyl unit. The onset, however, is very different for the four binding groups considered here: \ceNH2 gives a small and \ceNC a large negative value, respectively, whereas the onset is close to for sulfur and slightly positive for direct coupling. This results in different scenarios: Positive and increasing thermopower for thiol binding groups as well as direct coupling, and negative and increasing thermopower for isocyanide binding groups. For the amine linkers, the thermopower changes sign with increasing length, going from negative values for short chains to positive values for longer chains ( 4 methyl units).

Figure 3: Transmission curves for alkane chains with (a) thiol, (b) isocyanide and (c) amine binding groups and 2, 4, 6 and 8 methyl units. The dotted line indicates the Fermi level, which is set to .

The reasons for these different behaviours can be seen from the corresponding transmission curves presented in Fig. 3: For thiol end groups, a broad gateway state resonance originating from hybridization of states located on the sulfur atoms with \ceAu d-states appears within the HOMO-LUMO gap about below the Fermi level resulting in a negative slope of and thus a positive thermopower. With increasing chain length, the gap becomes deeper without changing the shape or position of the gateway state resonance. Since the center of the gap falls off faster than the resonance peak height, the slope of the transmission curve at becomes steeper and the thermopower larger.

This situation is similar for the chains with direct coupling (see Supporting Information), where the gateway resonance is much broader. For the shortest chain, however, the hybridization is so strong that the transmission becomes very flat with values close to 1 around the Fermi energy. In this case, the transport mechanism is resonant tunnelling and Eqs. (2) and (3) do not apply.Cheng et al. (2011) Therefore, the first data point in Fig. 2 is off the straight line.

For the isocyanides on the other hand, the transmission is dominated by a relatively sharp gateway state resonance slightly above the Fermi energy giving positive slopes and large negative Seebeck coefficients. As for the previous molecules, the peak height does not drop off as fast as the center of the gap for longer chains. Additionally, the center of this peak moves slightly up in energy. This results in an increase of the slope of the transmission curve and thus the thermopower.

In the case of amine-end groups, the Fermi level is located around the middle of the HOMO-LUMO gap, where the transmission curve is very flat. This is why the obtained values for S are small and even subtle variations in the position of a resonance imply a change of sign of the thermopower.

alkane-\ceNH2
alkane-\ceNC
alkane-\ceS
alkane-direct
alkene-\ceNH2
alkene-\ceNC
alkene-\ceS
alkene-direct
alkyne-\ceNH2
alkyne-\ceNC
alkyne-\ceS
alkyne-direct
Table 1: Fitted contact thermopower, , and length dependence, .

ii.2 Alkene chains

Figure 4: Same as Fig. 2 for alkene chains. For the amine-end groups, no linear fitting could be performed.

The results for alkene chains are plotted in Fig. 4. Again, the calculated values for the thermopower for the amine-end groups are very small, but we cannot see a linear dependence with length. In the corresponding transmission curves (see SI), a sharp very resonance appears at around above the Fermi energy, which arises from the LUMO of the molecule in the gas phase. This orbital is localized on the carbon backbone of the molecule and on the hydrogen atoms of the amine groups that stick out in the direction of the -system of the chain, thereby forming a hyperconjugated state. There is no contribution on the nitrogen atoms. Therefore, this state does not couple to the gold electrodes when the molecule is put into the junction in the given contact geometry. This situation makes it impossible to distinguish between backbone and end group effects on the transmission properties. A possible workaround would be to rotate the hydrogen atoms out of the -system of the backbone. However, we could not find an energetically stable conformation or junction geometry for that.

Figure 5: Transmission for alkene-\ceNC chains. Inset: Zoom in around the Fermi energy.

A new situation occurs for the isocyanide end-groups: the thermopower has a large negative onset and decreases even more with length. As shown in Fig. 5, there is a sharp resonance from the LUMO slightly above reaching a transmission value of . It does neither decrease nor shift in energy for longer chains. Therefore, as the transmission gap deepens, the slope of the transmission becomes steeper and the thermopower decreases.

For the thiol-end groups, the thermopower is again positive and increasing as for the alkane chains, but with larger values for the onset and slope of . A small negative contact thermopower is found for the directly coupled chains, although the thermopower itself is positive for all molecules and increases with length. As opposed to the alkanes, we find a large spread of values for with both negative and positive values.

ii.3 Alkyne chains

Figure 6: Same as Fig. 2 for alkyne chains. For the amine-end groups, no linear fitting could be performed.

The results for alkyne chains are presented in Fig. 6 and are similar to the alkenes. There are four different scenarios for the length dependence of the thermopower: positive and increasing (\ceS), negative and decreasing (\ceNC), negative and increasing (direct) and changing sign from positive to negative (\ceNH2).

All further transmission curves are given in the Supporting Information. The fitted values for the contact thermopower and the slope are listed in Table 1.

Iii Discussion

iii.1 Contact thermopower

By comparing all of our results, we do see a trend for the contact thermopower, , for the different binding groups: For thiols, it is positive for all systems studied here. This is in agreement with experimental findings for alkanedithiolsMalen et al. (2009); Guo, Zhou, and Tao (2013) and phenyldithiols.Reddy et al. (2007); Tan, Sadat, and Reddy (2010); Tan et al. (2011); Balachandran et al. (2013) Large negative values are obtained in all cases for isocyanide binding groups. This observation matches with measurements for \ceNC-terminated phenyls.Tan et al. (2011)

The situation is not entirely clear for the directly coupled chains: While all calculated values for the thermopower itself were positive, we find both positive and negative contact thermopowers in the linear fitting. This is a result from the appearance of very broad gateway resonance close to the Fermi energy, as mentioned in the previous section for alkane chains. For short lengths, resonant rather than off-resonant tunneling occurs.

The contact thermopower is expected to be around for amine-terminated chains, which would agrees with measurements on phenyldiaminesMalen et al. (2009) but no linear length dependence could be found for alkenes and alkynes. The calculated values for the thermopower itself, however, are all very close to .

iii.2 Length dependence

The slope, does not follow a simple rule of thumb, since neither its value nor even its sign can be directly predicted from the properties of the backbone alone. In particular, we find a large variation for the conjugated molecules. However, this lies within the predictions of the tunneling model as described in Ref. 26: exhibits a strong energy dependence and its derivative, , changes sign around the Fermi level. Since end groups can strongly influence the position of MO energies with respect to (alteration of the band lineup or level alignment) Xue and Ratner (2004); Zotti et al. (2010); Balachandran et al. (2013) they also effect . The length dependence of the thermopower is thus not only a property of the backbone, as previously assumed, but rather of the backbone + binding groups. In the presence of gateway states, it even depends on the whole junction.

Early measurements reported positive and decreasing thermopower with length for alkanedithiolsMalen et al. (2009), whereas positive and increasing values were measured more recentlyGuo, Zhou, and Tao (2013), the latter being consistent with our calculations. We note that we have also calculated the thermopower for \ceS-terminated chains in the top-geometry (as for the other binding groups). While the transmission curves look very different (see SI), we still find a positive and increasing thermopower. However, our findings suggest that the length dependence can be very sensitive to details in the coupling, which might explain the differences in the two experiments.

In the experiments of Ref. 48, a positive and increasing thermopower was found for directly \ceAu-\ceC coupled alkane chains. The dependence was nonlinear and a saturation of with molecular length was observed for a length of 10 methyl units. However, this data set contained only three points. Thus, we don’t see a contradiction to our results.

iii.3 Methodology

It is well known that DFT tends to overestimate conductances by up to several orders of magnitude and that higher levels of theory, such as the GW method, are required to bring calculated values in good agreement with experimental results.Strange and Thygesen (2011); Strange et al. (2011); Jin et al. (2013); Markussen, Jin, and Thygesen (2013) Even though the thermopower is most likely a more robust quantity, since it is found to be less sensitive to errors in method,Ke et al. (2009) DFT cannot be expected to give quantitatively accurate results. As discussed here, the thermopower depends mainly on the presence, position and shape of resonances from gateway states. DFT usually places molecular resonances too close to the Fermi level, leading to large errors in transport calculations. This is a consequence of the self-interaction error inherent in most functionals, which results in an inaccurate molecule-lead charge transfer.Ke, Baranger, and Yang (2007) Another problem is the underestimation of the HOMO-LUMO gap with DFT, resulting in far too low values for the conductance. Nevertheless, it is possible to predict qualitative trends. In particular, correct signs of the thermopower and the linear prefactor, , have been found.Ke et al. (2009); Quek et al. (2011); Widawsky et al. (2012)

Good results for the conduction and thermopower could also been obtained within the DFT +  approach, which corrects for static correlation and image charge effects.Quek et al. (2009, 2011); Widawsky et al. (2013); Kim et al. (2014) This method shifts all occupied (unoccupied) molecular orbitals down (up) in energy by a constant value. Thus, it moves the molecular resonances further away from the Fermi level, which results in a deepening of the transmission gap. It is, however, not clear how gateway states are affected. In any case, DFT +  fails for systems where charge transfer from the electrodes leads to partial occupation of otherwise unoccupied molecular orbitals.Khoo et al. (2015) This turned out to be the case for many of the systems in our study. A detailed comparison is given in the Supporting Information.

Much better results are expected for the GW method. This, however, is computationally extremely challenging and we were not able to converge results properly with respect to the basis set and the size of the interacting region for many of the structures.

Iv Conclusions

We have calculated the thermopower of small molecular chains attached to gold electrodes in well-defined and consistent binding geometries using Density Functional Theory. By varying the degree of conjugation and binding group, we have investigated the question to what extent their chemical properties determine the length dependence of the thermopower and if predictions can be made. A simple assumption that the thermopower depends linearly on the length and that this dependence is given by the nature of the molecular backbone does not hold. Instead, we find that in many cases resonances from gateway states dominate the electron transmission close to the Fermi level of the electrodes and thus the transport properties of the molecule. It is therefore the binding group, that not only dictates the contact thermopower, but also influences the length dependence to a large degree. A situation, where the gateway states are completely decoupled from the backbone in order to seperate their effects would be desirable, but could not be realized with reasonable junction geometries for many of the short molecules. The importance of gateway resonances is diminished for longer chains. However, we expect their influence to be neglibile only for much greater lengths than studied here. Overall, we can identify different scenarios for the length dependence. Amongst others, some where the thermopower changes sign with length.

Although DFT is not expected to yield quantitatively accurate results, we believe that it gives correct trends. Including static self-energy corrections in the so-called DFT +  approach did not prove to be feasible in this broad study, since its validity is restricted to cases where interactions between molecule and electrode surface are sufficiently low. Unfortunately, the self-consistent GW method did not turn out to be generally applicable either, as many of the calculations were practically impossible to converge.

In summary, we have demonstrated that, although it is very difficult to separate end group from backbone effects, the thermopower can be tuned chemically, even for the simplest systems. Other chemical modifications may lead to further variations in the values for the thermopower. However, the influence from the binding groups cannot be neglected for short molecules.

V Acknowledgements

This project has received fundings from the European Research Council under the European Union’s (EU) Seventh Framework Program (FP7/2007-2013)/ERC Grant Agreement No. 258806.

Vi Supplementary Information

Transmission curves obtained with DFT and DFT +  for all 60 molecules (including an alternative binding geometry for thiols).

References

  • Mahan and Sofo (1996) G. D. Mahan and J. O. Sofo, “The best thermoelectric.” Proc. Natl. Acad. Sci. U.S.A. 93, 7436–7439 (1996).
  • Humphrey and Linke (2005) T. E. Humphrey and H. Linke, “Reversible thermoelectric nanomaterials.” Phys. Rev. Lett. 94, 096601 (2005).
  • Wang, Segalman, and Majumdar (2006) R. Y. Wang, R. A. Segalman,  and A. Majumdar, “Room temperature thermal conductance of alkanedithiol self-assembled monolayers.” Appl. Phys. Lett. 89, 173113 (2006).
  • Venkataraman et al. (2006a) L. Venkataraman, J. E. Klare, C. Nuckolls, M. S. Hybertsen,  and M. L. Steigerwald, “Dependence of single-molecule junction conductance on molecular conformation.” Nature 442, 904–907 (2006a).
  • Venkataraman et al. (2007) L. Venkataraman, Y. S. Park, A. C. Whalley, C. Nuckolls, M. S. Hybertsen,  and M. L. Steigerwald, “Electronics and chemistry: Varying single-molecule junction conductance using chemical substituents.” Nano Lett. 7, 502–506 (2007).
  • Kiguchi et al. (2010) M. Kiguchi, H. Nakamura, Y. Takahashi, T. Takahashi,  and T. Ohto, ‘‘Effect of anchoring group position on formation and conductance of a single disubstituted benzene molecule bridging au electrodes: Change of conductive molecular orbital and electron pathway.” J. Phys. Chem. C 114, 22254–22261 (2010).
  • Capozzi et al. (2014) B. Capozzi, Q. Chen, P. Darancet, M. Kotiuga, M. Buzzeo, J. B. Neaton, C. Nuckolls,  and L. Venkataraman, “Tunable charge transport in single-molecule junctions via electrolytic gating.” Nano Lett. 14, 1400–1404 (2014).
  • Strange et al. (2015) M. Strange, J. S. Seldenthuis, C. J. O. Verzijl, J. M. Thijssen,  and G. C. Solomon, “Interference enhanced thermoelectricity in quinoid type structures.” J. Chem. Phys. 142, 084703 (2015).
  • Ho Choi, Kim, and Frisbie (2008) S. Ho Choi, B. Kim,  and C. D. Frisbie, “Electrical resistance of long conjugated molecular wires.” Science 320, 1482–1486 (2008).
  • Solomon et al. (2008) G. C. Solomon, D. Q. Andrews, R. H. Goldsmith, T. Hansen, M. R. Wasielewski, R. P. Van Duyne,  and M. A. Ratner, “Quantum interference in acyclic systems: Conductance of cross-conjugated molecules.” J. Am. Chem. Soc. 130, 17301–17308 (2008).
  • Solomon et al. (2009) G. C. Solomon, D. Q. Andrews, R. P. Van Duyne,  and M. A. Ratner, “Electron transport through conjugated molecules: When the system only tells part of the story.” ChemPhysChem 10, 257–64 (2009).
  • Hines et al. (2010) T. Hines, I. Diez-Perez, J. Hihath, H. Liu, Z.-S. Wang, J. Zhao, G. Zhou, K. Müllen,  and N. Tao, “Transition from tunneling to hopping in single molecular junctions by measuring length and temperature dependence.” J. Am. Chem. Soc. 132, 11658–11664 (2010).
  • Lee et al. (2013) W. Lee, K. Kim, W. Jeong, L. A. Zotti, F. Pauly, J. C. Cuevas,  and P. Reddy, “Heat dissipation in atomic-scale junctions.” Nature 498, 209–212 (2013).
  • Xu and Tao (2003) B. Xu and N. J. Tao, “Measurement of single-molecule resistance by repeated formation of molecular junctions.” Science 301, 1221–1223 (2003).
  • Venkataraman et al. (2006b) L. Venkataraman, J. E. Klare, I. W. Tam, C. Nuckolls, M. S. Hybertsen,  and M. L. Steigerwald, “Single-molecule circuits with well-defined molecular conductance,” Nano Lett. 6, 458–462 (2006b).
  • Tan et al. (2011) A. Tan, J. Balachandran, S. Sadat, V. Gavini, B. D. Dunietz, S.-Y. Jang,  and P. Reddy, ‘‘Effect of length and contact chemistry on the electronic structure and thermoelectric properties of molecular junctions.” J. Am. Chem. Soc. 133, 8838–8841 (2011).
  • Cheng et al. (2011) Z.-L. Cheng, R. Skouta, H. Vazquez, J. R. Widawsky, S. Schneebeli, W. Chen, M. S. Hybertsen, R. Breslow,  and L. Venkataraman, “In situ formation of highly conducting covalent au-c contacts for single-molecule junctions.” Nat. Nanotechnol. 6, 353–357 (2011).
  • Kristensen et al. (2008) I. S. Kristensen, D. J. Mowbray, K. S. Thygesen,  and K. W. Jacobsen, “Comparative study of anchoring groups for molecular electronics: Structure and conductance of au-s-au and au-nh-au junctions,” J. Phys.: Condens. Matter 20, 374101 (2008).
  • Neaton, Hybertsen, and Louie (2006) J. B. Neaton, M. S. Hybertsen,  and S. G. Louie, “Renormalization of molecular electronic levels at metal-molecule interfaces.” Phys. Rev. Lett. 97, 216405 (2006).
  • Widawsky et al. (2012) J. R. Widawsky, P. Darancet, J. B. Neaton,  and L. Venkataraman, “Simultaneous determination of conductance and thermopower of single molecule junctions.” Nano Lett. 12, 354–358 (2012).
  • Kim et al. (2014) T. Kim, P. Darancet, J. R. Widawsky, M. Kotiuga, S. Y. Quek, J. B. Neaton,  and L. Venkataraman, ‘‘Determination of energy level alignment and coupling strength in 4,4’-bipyridine single-molecule junctions.” Nano Lett. 14, 794–798 (2014).
  • Khoo et al. (2015) K. H. Khoo, Y. Chen, S. Li,  and S. Y. Quek, “Length dependence of electron transport through molecular wires – a first principles perspective.” Phys. Chem. Chem. Phys. 17, 77–96 (2015).
  • Reddy et al. (2007) P. Reddy, S.-Y. Jang, R. A. Segalman,  and A. Majumdar, ‘‘Thermoelectricity in molecular junctions.” Science 315, 1568–1571 (2007).
  • Malen et al. (2009) J. A. Malen, P. Doak, K. Baheti, T. D. Tilley, R. A. Segalman,  and A. Majumdar, “Identifying the length dependence of orbital alignment and contact coupling in molecular heterojunctions.” Nano Lett. 9, 1164–1169 (2009).
  • McConnell (1961) H. M. McConnell, “Intramolecular charge transfer in aromatic free radicals,” J. Chem. Phys. 35, 508–515 (1961).
  • Karlström, Strange, and Solomon (2014) O. Karlström, M. Strange,  and G. C. Solomon, “Understanding the length dependence of molecular junction thermopower.” J. Chem. Phys. 140, 044315 (2014).
  • Pauly, Viljas, and Cuevas (2008) F. Pauly, J. K. Viljas,  and J. C. Cuevas, “Length-dependent conductance and thermopower in single-molecule junctions of dithiolated oligophenylene derivatives: A density functional study.” Phys. Rev. B 78, 035315 (2008).
  • Dell et al. (2015) E. J. Dell, B. Capozzi, J. Xia, L. Venkataraman,  and L. M. Campos, “Molecular length dictates the nature of charge carriers in single-molecule junctions of oxidized oligothiophenes.” Nat. Chem. 7, 209–214 (2015).
  • Meir and Wingreen (1992) Y. Meir and N. S. Wingreen, “Landauer formula for the current through an interacting electron region.” Phys. Rev. Lett. 68, 2512–2515 (1992).
  • Taylor, Guo, and Wang (2001) J. Taylor, H. Guo,  and J. Wang, “Ab initio modeling of quantum transport properties of molecular electronic devices.” Phys. Rev. B 63, 245407 (2001).
  • Xue, Datta, and Ratner (2002) Y. Xue, S. Datta,  and M. A. Ratner, “First-principles based matrix green’s function approach to molecular electronic devices: General formalism.” Chem. Phys. 281, 151–170 (2002).
  • Brandbyge et al. (2002) M. Brandbyge, J.-L. Mozos, P. Ordejón, J. Taylor,  and K. Stokbro, “Density-functional method for nonequilibrium electron transport.” Phys. Rev. B 65, 165401 (2002).
  • Thygesen, Bollinger, and Jacobsen (2003) K. S. Thygesen, M. V. Bollinger,  and K. W. Jacobsen, “Conductance calculations with a wavelet basis set.” Phys. Rev. B 67, 115404 (2003).
  • Butcher (1990) P. N. Butcher, ‘‘Thermal and electrical transport formalism for electronic microstructures with many terminals.” J. Phys.: Condens. Matter 2, 4869–4878 (1990).
  • Paulsson and Datta (2003) M. Paulsson and S. Datta, “Thermoelectric effect in molecular electronics.” Phys. Rev. B 67, 241403 (2003).
  • Segal (2005) D. Segal, “Thermoelectric effect in molecular junctions: A tool for revealing transport mechanisms.” Phys. Rev. B 72, 165426 (2005).
  • Malen et al. (2010) J. A. Malen, S. K. Yee, A. Majumdar,  and R. A. Segalman, “Fundamentals of energy transport, energy conversion, and thermal properties in organic-inorganic heterojunctions.” Chem. Phys. Lett. 491, 109–122 (2010).
  • Enkovaara et al. (2010) J. Enkovaara, C. Rostgaard, J. J. Mortensen, J. Chen, M. Dułak, L. Ferrighi, J. Gavnholt, C. Glinsvad, V. Haikola, H. A. Hansen, H. H. Kristoffersen, M. Kuisma, A. H. Larsen, L. Lehtovaara, M. Ljungberg, O. Lopez-Acevedo, P. G. Moses, J. Ojanen, T. Olsen, V. Petzold, N. A. Romero, J. Stausholm-Møller, M. Strange, G. A. Tritsaris, M. Vanin, M. Walter, B. Hammer, H. Häkkinen, G. K. H. Madsen, R. M. Nieminen, J. K. Nørskov, M. Puska, T. T. Rantala, J. Schiøtz, K. S. Thygesen,  and K. W. Jacobsen, “Electronic structure calculations with gpaw: A real-space implementation of the projector augmented-wave method.” J. Phys.: Condens. Matter 22, 253202 (2010).
  • Perdew, Burke, and Ernzerhof (1996) J. P. Perdew, K. Burke,  and M. Ernzerhof, “Generalized gradient approximation made simple,” Phys. Rev. Lett. 77, 3865–3868 (1996).
  • Strange et al. (2011) M. Strange, C. Rostgaard, H. Häkkinen,  and K. S. Thygesen, “Self-consistent gw calculations of electronic transport in thiol- and amine-linked molecular junctions.” Phys. Rev. B 83, 115108 (2011).
  • Thygesen and Jacobsen (2005) K. S. Thygesen and K. W. Jacobsen, “Interference and k-point sampling in the supercell approach to phase-coherent transport.” Phys. Rev. B 72, 033401 (2005).
  • Strange et al. (2008) M. Strange, I. S. Kristensen, K. S. Thygesen,  and K. W. Jacobsen, “Benchmark density functional theory calculations for nanoscale conductance.” J. Chem. Phys. 128, 114714 (2008).
  • Guo, Zhou, and Tao (2013) S. Guo, G. Zhou,  and N. Tao, “Single molecule conductance, thermopower, and transition voltage.” Nano Lett. 13, 4326–4332 (2013).
  • Tan, Sadat, and Reddy (2010) A. Tan, S. Sadat,  and P. Reddy, “Measurement of thermopower and current-voltage characteristics of molecular junctions to identify orbital alignment.” Appl. Phys. Lett. 96, 013110 (2010).
  • Balachandran et al. (2013) J. Balachandran, P. Reddy, B. D. Dunietz,  and V. Gavini, “End-group influence on frontier molecular orbital reorganization and thermoelectric properties of molecular junctions.” J. Phys. Chem. Lett. 4, 3825–3833 (2013).
  • Xue and Ratner (2004) Y. Xue and M. A. Ratner, “End group effect on electrical transport through individual molecules: A microscopic study.” Phys. Rev. B 69, 085403 (2004).
  • Zotti et al. (2010) L. A. Zotti, T. Kirchner, J.-C. Cuevas, F. Pauly, T. Huhn, E. Scheer,  and A. Erbe, “Revealing the role of anchoring groups in the electrical conduction through single-molecule junctions.” Small 6, 1529–1535 (2010).
  • Widawsky et al. (2013) J. R. Widawsky, W. Chen, H. Vázquez, T. Kim, R. Breslow, M. S. Hybertsen,  and L. Venkataraman, “Length-dependent thermopower of highly conducting au-c bonded single molecule junctions.” Nano Lett. 13, 2889–2894 (2013).
  • Strange and Thygesen (2011) M. Strange and K. S. Thygesen, ‘‘Towards quantitative accuracy in first-principles transport calculations: The gw method applied to alkane/gold junctions.” Beilstein J. Nanotechnol. 2, 746–754 (2011).
  • Jin et al. (2013) C. Jin, M. Strange, T. Markussen, G. C. Solomon,  and K. S. Thygesen, “Energy level alignment and quantum conductance of functionalized metal-molecule junctions: Density functional theory versus gw calculations.” J. Chem. Phys. 139, 4307 (2013).
  • Markussen, Jin, and Thygesen (2013) T. Markussen, C. Jin,  and K. S. Thygesen, “Quantitatively accurate calculations of conductance and thermopower of molecular junctions.” Phys. Status Solidi B 250, 2394–2402 (2013).
  • Ke et al. (2009) S.-H. Ke, W. Yang, S. Curtarolo,  and H. U. Baranger, “Thermopower of molecular junctions: An ab initio study.” Nano Lett. 9, 1011–1014 (2009).
  • Ke, Baranger, and Yang (2007) S.-H. Ke, H. U. Baranger,  and W. Yang, ‘‘Role of the exchange-correlation potential in ab initio electron transport calculations.” J. Chem. Phys. 126, 201102 (2007).
  • Quek et al. (2011) S. Y. Quek, H. J. Choi, S. G. Louie,  and J. B. Neaton, “Thermopower of amine-gold-linked aromatic molecular junctions from first principles.” ACS Nano 5, 551–557 (2011).
  • Quek et al. (2009) S. Y. Quek, H. J. Choi, S. G. Louie,  and J. B. Neaton, ‘‘Length dependence of conductance in aromatic single-molecule junctions.” Nano Lett. 9, 3949–3953 (2009).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
316066
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description