Controlling Band Alignment In Molecular Junctions: Utilizing Two

26 mins ago - We illustrate the effectiveness of this strategy by investigating the thermoelectric properties of the junctions with a model molecular ...
3 downloads 0 Views 2MB Size
Subscriber access provided by University of Massachusetts Amherst Libraries

C: Energy Conversion and Storage; Energy and Charge Transport

Controlling Band Alignment In Molecular Junctions: Utilizing Two-Dimensional Transition-Metal Dichalcogenides As Electrodes For Thermoelectric Devices Chengjun Jin, and Gemma C. Solomon J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b00464 • Publication Date (Web): 30 May 2018 Downloaded from http://pubs.acs.org on May 30, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Controlling Band Alignment In Molecular Junctions: Utilizing Two-dimensional Transition-Metal Dichalcogenides As Electrodes For Thermoelectric Devices Chengjun Jin and Gemma C. Solomon∗ Nano-Science Center and Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark E-mail: [email protected]

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract While chemical space is vast, there are many instances where chemical modifications make only insignificant changes to the current that can pass through a molecule. This insensitivity comes, in part, from the energy mismatch between molecular resonances and the Fermi level of the electrodes used. Here we present a strategy to overcome this problem by employing two-dimensional transition-metal dichalcogenides as electrodes. The work function of the electrodes can be tuned across the entire molecular energy range (from the highest occupied molecular orbital to the lowest unoccupied molecular orbital) at low bias with the appropriate choice of electrode material. We illustrate the effectiveness of this strategy by investigating the thermoelectric properties of the junctions with a model molecular system, as optimal thermoelectric performance requires a delicate balance between the electronic and heat transport properties. By using Van der Waals contacts between the binding groups and the electrodes, similar to the binding groups used in graphene junctions, we find that we can effectively suppress the phonon contribution to the heat transport while achieving high levels of electron transport and thermopower as the molecular level is tuned into near-resonance.

Introduction A long standing problem for the field of molecular electronics is the question of band-lineup or the alignment of the electrode Fermi energy and the molecular energy levels of interest. 1–7 This problem has captured the attention of chemists and physicists alike as it goes to the heart of understanding molecular electronic structure in the non-traditional environment of a single molecule bound between conducting electrodes and ensuring that the electronic functionality of interest is accessible in the experimentally relevant energy/bias range. One aspect of the problem is the treatment of the electronic structure, and considerable effort has been expended in developing more sophisticated electronic structure methods, such as GW 4,6 for transport applications. The other relates to the inherent mismatch between most

2

ACS Paragon Plus Environment

Page 2 of 24

Page 3 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

molecular resonances (where molecular conductance generally becomes more significant and chemical effects more pronounced) and the Fermi energy of the (usually gold) electrodes. While there has been some hope that gating of molecular energy levels might be able to address this issue by tuning the molecule into resonance, this solution has proven to be at best challenging and at worst totally impractical. 8,9 In this paper, we propose that alternative electrode materials will allow the band-lineup to be tuned across the experimentally relevant energy range at low bias. The range of electrode materials that have been employed to date is limited, with most effort focussing on metals (generally gold) 10,11 and some exploration of graphene. 12–16 Here we employ electrode materials that can effectively tune the band-lineup: two-dimensional transition-metal dichalcogenides (TMDs). We hypothesize that the combination of tunable band-lineup and a π-stacked molecular junction will offer significant scope for optimizing the properties of the system and test this with an application to molecular thermoelectrics as the delicate balance required between the electronic and heat transport properties makes this problem particularly challenging.

Molecular thermoelectrics - background A single molecule connected between two metallic electrodes has been proposed as a device for achieving a high figure of merit ZT. 17–20 ZT is defined as S 2 GT /κtotal , which comprises the temperature T, the thermopower S, the electronic conductance G and the thermal conductance κtotal = κel + κph , the sum of the electronic and phononic contributions. Strategies to find devices with ZT > 1, which are considered to be good thermoelectrics, 21 either suppress the thermal conductance κtotal or enhance the power factor S 2 G, or both. In molecular junctions, the typical thermopower S is around 10 µV/K, where S > 150 µV/K is considered for practical applications, 22,23 with an exceptionally high thermopower of 300 µV/K and high ZT of 4 predicted in the edge-over-edge zinc-porphyrin molecular wires. 24 In molecular thermoelectric devices, the phonon propagation needs to be suppressed due 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

to the reduced dimensionality 25,26 or S 2 G needs to be enhanced due to the tunability of the energy level alignment at the molecule/metal interfaces. However it has been difficult to simultaneously control these two factors due to the correlation between the electronic and phononic properties. 27 π-stacked molecular junctions show significant phonon suppression as very few phonon modes can propagate across the non-bonded interface 28,29 and the ZT in such systems can be further improved if the electronic thermal conductance is reduced. 29 We note here that while the prior work on π-stacked molecular junctions for thermoelectrics is theoretical in nature, the geometries employed in that work and this present study are modeled on the structures hypothesized to exist in π-stacked junctions with graphene electrodes. 12,16

A strategy to tune band-lineup

Figure 1: An outline of the strategy employed in this paper. (a) The atomistic structure of the MoX2 (X= O, S, Se, Te) molecular junction. A six-carbon alkyne chain molecule is attached to two semi-infinite MoX2 electrodes via the Van der Waals forces. The chemical drawing of the molecule is shown below. Thermoelectric figure of merit can be improved by (b) a resonant shift in the electronic transmission function (from solid line to dashed line) to increase the power factor and (c) a suppression of the phononic transmission function (from solid line to dashed line). Figure 1 outlines the basic strategy proposed in this work to bring electronic features of interest into the low bias window. A model molecular system (a six-carbon alkyne chain) is connected to the two-dimensional TMDs electrodes via the Van der Waals interactions with anthracene linker groups, as shown in the Fig. 1 (a). The non-bonded contact provides 4

ACS Paragon Plus Environment

Page 4 of 24

Page 5 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

excellent electrical contact while at the same time being problematic for heat transport, making it possible to simultaneously enhance of the power factor S 2 G and the suppress of the thermal conductance κtotal . On one hand, the energy level alignment at the molecule/metal interfaces can be tuned by utilizing different metallic TMD materials with different work functions, as illustrated in Fig. 1 (b). On the other hand, due to the non-bonding nature at the molecule/metal interfaces, the heat transport is significantly blocked, as illustrated in Fig. 1 (c). Consequently, we find that molecular junctions utilizing certain metallic TMDs can allow ZT to be tuned across a significant range, in this case between almost 0 and 1.2 at room temperature, without any optimization of the molecular component (beyond the choice of a non-bonded system). More specifically, we have not fully optimized the system in terms of the values of S, G and κel in the energy window accessible with the TMDs.

Computational methods The junction structures are optimized until the forces are smaller than 0.02 eV/˚ A using density functional theory (DFT) with the Van der Waals exchange correlational functional 2 (vdW-DF2), 30 as implemented in the GPAW electronic structure code. 31 This gives the binding distances of ∼ 3.4 ˚ A for all the junctions, in agreement with the case of graphene. 32 The wave functions are sampled on a real space grid for optimizations. Double ζ polarized basis sets and PBE xc-functional are employed in the electronic transport calculations. The electronic transport calculations are performed with the Green’s function approach as implemented in the GPAW. 33 The electronic transmission function is calculated from the Landauer formula 34,35

Tel (E) = Tr[G(E)ΓL (E)G† (E)ΓR (E)],

5

ACS Paragon Plus Environment

(1)

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 24

where the retarded Green’s function is obtained from  −1 G(E) = ES − H − ΣL/R (E) .

(2)

Here S and H denote the overlap matrix and the Kohn-Sham (KS) Hamiltonian matrix in the atomic orbital basis, 36 respectively. ΓL/R = i[ΣL/R − Σ†L/R ] accounts for the level broadening due to the electrode coupling, where ΣL/R denotes the retarded lead self-energy of the left/right electrodes. The phononic transport calculations are performed in the Atomistix ToolKit version 2015.1, QuantumWise A/S. 37–39 Double ζ polarized basis sets and PBE xc-functional are employed in the phononic transport calculations. The phononic transmission function Tph (ω) can be calculated from Eq. 1 with the substitution H → K and S → ω 2 M, where K and M denote the force constant matrix and a diagonal matrix with the atomic masses. The force matrix K is calculated by the finite difference method as

KIµ,Jν =

FJν (QIµ ) − FJν (−QIµ ) , 2QIµ

(3)

which is constructed in a such way that each atom is displaced by QIµ = ±0.02 ˚ A in the direction µ = {x, y, z} to obtain the forces FJν (QIµ ) on the atom J 6= I in the direction ν. The intra-atomic elements are calculated by imposing momentum conservation KIµ,Iν = P − K6=I KIµ,Kν . Additionally, the phononic thermal conductance can be calculated from the Landauer-type formula 40

κph

h ¯ = 2πkB T 2

Z



dωω 2 Tph (ω)

0

e¯hω/kB T . (e¯hω/kB T − 1)2

(4)

The dimensionless figure of merit ZT = S 2 GT /(κel + κph ) is calculated to give a measure of the themoelectric properties of the molecular junctions, where T is the temperature.

6

ACS Paragon Plus Environment

Page 7 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Defining an integral 41 2 Lm (µ) = h



  ∂f (E, µ, T ) dETel (E)(E − µ) − , ∂E −∞

Z

m

(5)

where f (E, µ, T ) is the Fermi-Dirac distribution function, allows all the relevant quantities to be calculated. The electronic conductance, the thermopower and the electronic thermal conductance can be calculated in terms of L as

G = e2 L0 (EF ),

(6)

L1 (EF ) , eT L0 (EF )

(7)

1 L2 (L2 − 1 ). T L0

(8)

S=

κel =

We note that the quantities in Equ. 6, 7 and 8 are defined within the linear response limit, and are thus applicable when ∆T /T is small, where ∆T = TL − TR  T and T is the averaged temperature of the left and right electrodes. The non-linear effects are expected to be of minor importance, since the experiments are typically done with T ≈ 300 K and |∆T | ≤ 30 K.

Results and discussion The choice of TMD Two-dimensional transition-metal dichalcogenides (TMDs) exist in either metallic or semiconducting phases, depending on the experimental details. 43 While the semiconducting phases have been widely studied, the metallic phases are of interest in our study because TMDs can then function as metallic electrodes in our molecular junctions. A recent com7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

XY-averaged VH (eV)

(a) Φ=5.03 eV

Z (Å)

(b)

−1

Molecule

−2

TMDs

−3

E - Evac (eV)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 24

-EA

MoTe2

−4

MoSe2

−5

MoS2

−6 MoO2

−7 −8

-IP

−9

Energy levels

Figure 2: (a) The XY-averaged Hartree potential along the Z direction for the MoS2 slab. The black dashed line indicates the Fermi level. The work function is indicated by the red dotted line. The atomistic structure of MoS2 is embedded to show its position in the slab calculation. (b) The electron affinity (EA) and ionization potential (IP) of the gas-phase molecule and the work function of 216 of free-standing metallic TMDs. The electron affinity (EA) and ionization potential (IP) are calculated from the ∆SCF method. 42 The work functions of MoX2 (X = O, S, Se, Te) are color-coded.

8

ACS Paragon Plus Environment

Page 9 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

putational screening study on the two-dimensional transition-metal dichalcogenides (TMDs) shows that many TMD phases are metallic and thermodynamically stable. 44 Moreover it is experimentally feasible to fabricate these defect-free metallic TMDs. 43,45,46 The metallic phase has shown many interesting applications such as the low-contact resistance, 43 the high hydrogen evolution reaction activity. 47 To the best of our knowledge, we are the first to calculate the work functions of these metallic TMDs using a DFT slab model, which has been very successfully applied to calculate the work functions of conventional metals. 48 The work function of a free-standing TMD is calculated with a slab model, as shown in Fig. 2 (a) for the case of MoS2 . The work function is obtained from the difference between the Fermi level and the Hartree potential in the vacuum. Fig. 2 (b) shows the work functions of 216 metallic free-standing TMDs, 49 where the work functions of MoX2 (X = O, S, Se, Te) are highlighted. This distribution has a mean of 5.4 eV and a standard deviation of 1.1 eV, ranging from 3.6 eV to 8.1 eV. Additionally the electron affinity (EA) and ionization potential (IP) of the gas-phase molecule, calculated from the ∆SCF method, 42 are also shown in Fig. 2 (b). From the simple picture where the molecular levels and the work function of TMD are rigidly combined, we expect to have different energy level alignments at the molecule/metal interfaces because of the wide range of TMD work functions. In the molecular system of interest here, the lowest unoccupied molecular orbital (LUMO) resonance can be probed by using TMDs with small work functions. Consequently it is possible to tune the power factor S 2 G as illustrated in Fig. 1 (b). We note that the realistic energy level alignment at the molecule/metal interface is more complicated than this simple picture because of the molecule-metal interactions at the interface. 50 Therefore the energy level alignment in the metal/molecule/metal junctions are investigated using DFT as the next step.

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Electron Transport The ability to tune the energy level alignment between the metal Fermi level and the molecular levels, gives the possibility of changing the zero-bias electronic conductance in a controlled way. Here, four metallic TMDs with different work functions (4.2 eV to 6.6 eV) are considered, namely MoX2 (X= O, S, Se, Te) as shown in Fig. 2 (b), to test the hypothesis that controlled energy level alignment can be used to optimize thermoelectric properties. Because of the weak Van der Waals interactions between the molecule and the TMD, the energy level alignment at the metal/molecule interfaces are approximately the superposition of the energy levels of the gas-phase molecule and the energy levels of the TMD. 51 This is reflected in the fact that the electronic transport gaps among the molecular junctions are similar to the Kohn-Sham highest occupied molecular orbital and lowest unoccupied molecular orbital (KS HOMO-LUMO) gap of the gas-phase molecule, which is ∼ 1.4 eV indicated by the vertical lines in Fig. 3 (a). While it is clear that while the transmission resonances are, to a large extent, simply shifted with the changing work function, the transmission line shapes also change significantly. The broadenings of the molecular energy levels (responsible for the line shape) are controlled by the strength of the interaction between the molecule and the electrode, and the electrode density of states (DOS). The the DOS of the MoX2 around the Fermi energy is dominated by Mo 4d states and X p states, with qualitatively similar features shifted in a similar fashion to the work function shift. Importantly the DOS of MoX2 is not flat around the Fermi energy, so the work function shift also induces a change in the broadening of the resonances and the variation seen in the transmission line shapes. From the point-of-view of a one-level Lorentzian model for the molecular transport, the position of the LUMO is more important for controlling the magnitude of the electronic conductance than the energy level broadening, as reflected in the electronic transmission shown in the Fig. 3 (a), where the position of the LUMO is indicated by the dashed vertical line and the Fermi level is at 0 eV. From MoO2 to MoTe2 , the Fermi level moves towards 10

ACS Paragon Plus Environment

Page 10 of 24

Page 11 of 24

(a)

10 0 10 -2

MoO2

10 0

T el

10 -2

MoS2

10 0 10 -2

MoSe2

10 0 10 -2 10 -4 −1.5

(b)

MoTe 2 −1.0

−0.5

0.0

0.5

1.0

1.5

5.5

6.0

6.5

7.0

E - E f (eV)

10 0 10 -1

G (G 0 )

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

10 -2 10 -3 10 -4 4.0

4.5

5.0

Work function Φ (eV)

Figure 3: (a) The electronic transmission functions Tel (E) of MoX2 (X = O, S, Se, Te) junctions. The vertical lines indicate the corresponding HOMO and LUMO positions and the Fermi energy is set at 0 eV. (b) The zero-bias electronic conductance as a function of the work function of MoX2 (X = O, S, Se, Te) junctions.

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the LUMO as a result of the decreasing work function of the TMDs. According to Eq. 6, the zero-bias electronic conductance can be enhanced when the work function of the TMDs decreases as shown in Fig. 3 (b). The electronic conductance at the Fermi energy varies from 10−4 G0 to 10−1 G0 across the series.

Heat transport The weak Van der Waals interactions between the molecule and the electrode are expected to block the heat transport at the molecule/electrode interfaces as they did for the case of graphene electrodes, 29 and indeed this is found to be the case. The phononic transmission functions are shown in Fig. 4 (a). It is clear that the phononic transmission is suppressed across the whole energy range and is most significant below 100 meV, which is similar to the case of graphene as electrodes. 29 The thermal conductance varies from 14.2 pW/K to 38.2 pW/K at room temperature as shown in Fig. 4 (b). Of all the systems considered, MoO2 exhibits the largest thermal conductance as there is reduced mass miss-match between O and C atoms (compared with either Mo or other choices for X). The Debye frequencies of bulk MoX2 (X = O, S, Se, Te) calculated from DFT are ∼ 90 meV, 53 meV, 42 meV and 38 meV, respectively, in agreement with the cut-off observed in the phonon transmission functions. In the the larger X systems (MoS2 , MoSe2 and MoTe2 ), we effectively have two mechanisms of phonon transport suppression in play: the effect of non-bonded contact and energy mismatch between the phonon modes in the electrode and the molecule.The suppression of the phonon thermal conductance to this extent means that the total thermal conductance is now dominated by the electronic part. This leads to a significant modulation of ZT, as shown later.

Thermoelectric properties We now discuss the thermoelectric properties in terms of the thermopower S, the power factor S 2 G, the thermal conductance κtotal = κe + κph and the figure of merit ZT. Fig. 5 shows 12

ACS Paragon Plus Environment

Page 12 of 24

Page 13 of 24

(a) 0.8 0.4

MoO 2

T ph (ω)

0.0 0.8 0.4

MoS 2

0.0 0.8 0.4

MoSe 2

0.0 0.8 0.4 0.0

(b) κph (pW/K)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

MoTe 2 0

50

100

150

Phonon energy ħω (meV)

200

50 40 30 20 10 0

0

200

400

600

Temperature (K)

800

1000

Figure 4: (a) The phononic transmission functions Tph (ω) of MoX2 (X = O, S, Se, Te) junctions. (b) the corresponding thermal conductance κph as a function of the temperature.

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

these quantities for four junctions as a function of energy. Because typical thermoelectric measurements are performed at low bias voltages, the thermoelectric properties around the

200 100 0 −100 −200

MoO 2

MoS 2

MoSe 2

MoTe 2

15 10

κtotal

(nW/K)

b

S 2 G (k 2 /h)

S

(µV/K)

Fermi energy will be the focus here.

5 0 5 4 3 2 1 0

1.5

ZT

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 24

1.0 0.5 0.0

−1.0−0.5 0.0 0.5 1.0−1.0−0.5 0.0 0.5 1.0−1.0−0.5 0.0 0.5 1.0−1.0−0.5 0.0 0.5 1.0

E - E (eV) f

E - E (eV)

E - E (eV)

f

f

E - E (eV) f

Figure 5: The thermopower S, the power factor S 2 G, the thermal conductance κtotal and the thermoelectric figure of merit ZT of MoX2 (X = O, S, Se, Te) junctions. The thermopower can be calculated from the slope of the transmission function at the Fermi energy, see Eq. 7. The four junctions have a negative thermopower, as expected because there is the LUMO-mediated transport in all cases. MoS2 has the largest thermopower while MoTe2 has thermopower of nearly 0 µV/K as a result of zero slope in the electronic transmission function of MoTe2 at the Fermi energy (the system is on resonance). From MoO2 to MoTe2 , the thermopower increases to the maximum (MoS2 ) and is gradually decreased to 0 µV/K (MoTe2 ). The magnitude of the thermopower plays a critical rule in determining the power factor. It results in a nearly zero power factor for MoTe2 although its electronic conductance is close to 1 G0 . The combination of large S 2 and large G leads 14

ACS Paragon Plus Environment

Page 15 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

to the power factor of up to 9.3 k2b /h in the case of MoS2 . The total thermal conductance is, in general, very low ranging from 0.2 to 4.1 nW/K. Because the phononic part is suppressed in these systems to the order of pW/K, the main contribution comes from the electronic part on the order of nW/K. The large electronic thermal conductance can be understood from the Wiedemann-Franz law, 52 recently shown to hold down to metallic nanowires 53 where the κel is proportional to the G. Additionally, the operation of the Wiedemann-Franz law is also reflected in the fact that the κtotal is positively correlated with the G from MoO2 to MoTe2 . By taking the ratio between S 2 G and κtotal to calculate ZT, we can see that the MoS2 shows the most promising high ZT up to 1.2 at the Fermi energy. The nearly zero ZT in the MoTe2 case is due to the nearly zero S, while the extremely small κtotal in the MoO2 case does not leads to a large ZT due to the small S 2 G.

Conclusions In conclusion, the tunability of energy level alignment at the molecule/metal interfaces is achieved by tailoring the work function of the electrode. As a result, a large power factor S 2 G is achieved without further system optimization. On the other hand, the thermal conductance is significantly suppressed due to the weak Van der Waals interactions between the molecule and the TMD materials. We show that with the appropriate choice of electrodes, ZT for a given molecule can be tuned significantly, in this case from almost zero to 1.2, with the possibility of achieving an even higher ZT with optimization of the molecular bridge. This opens the possibility of accessing whichever energy offers the most promising combination of S, G and κel for a given molecule across a wide energy range. While there is undoubtedly a possibility of further optimization of the molecular S and G, the limiting factor for the molecules considered here appears to be the magnitude of κel . A dramatic reduction in κel without a corresponding reduction in G requires a violation of the Wiedemann-Franz law. This is outside the scope of this work but remains a very interesting

15

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

direction for future investigation in molecular thermoelectrics. While we investigated thermoelectric properties here as they are particularly sensitive to band lineup, the principles we have outlined are general. Although gold electrodes have been the popular choice for many molecular electronics investigations, molecular features of interest do not always occur in the energy range accessible from the gold Fermi energy. The series of electrodes studied here show us just how sensitive thermoelectric properties can be to band-lineup, for example in the stark difference between ZT calculated for MoSe2 and MoTe2 . We should stress here however, that the limitations in the accuracy of any DFT calculation for these systems mean that we cannot make definitive predictions, for example discriminating between the choice of MoSe2 or MoTe2 for an experiment. While the many-body perturbation methods such as GW could result in more accurate energy level alignments at the weaklycoupled molecule/metal interface, 4,6 this remains a challenge for large molecules in junctions. Nonetheless, we demonstrate the feasibility of tuning thermoelectric properties in this study and anticipate that DFT is able to qualitatively capture the trends. The binding groups that have been developed to allow Van der Waals contact between molecules and graphene electrodes can be extended to TMDs and this extension allows molecular properties to be tuned into (near) resonance. Many different electronic effects, be they due to conformational switching, chemical substituents or different binding geometries, are all expected to be more significant if the systems can be taken out of the far-off-resonant regime. The chemical identity of the electrode will certainly have consequences for the details of the transmission, seen here in the transmission line shapes, but for these systems, the first order effect was simply to shift the energy levels relative to the Fermi energy. While the use of TMD electrodes will certainly present challenges, as did metal electrodes and graphene electrodes, these systems offer the possibility of a better energetic match between the Fermi energy and molecular energy levels hinting at increased chemical control of transport properties.

16

ACS Paragon Plus Environment

Page 16 of 24

Page 17 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Acknowledgements This work was supported by the funding from the Villum Foundation, Carlsberg Foundation and the Danish Council for Independent Research — Natural Sciences.

References (1) Xue, Y.; Datta, S.; Ratner, M. A. Charge transfer and ”band lineup” in molecular electronic devices: A chemical and numerical interpretation. The Journal of Chemical Physics 2001, 115, 4292–4299. (2) Perrin, M. L.; Verzijl, C. J. O.; Martin, C. A.; Shaikh, A. J.; Eelkema, R.; van Esch, J. H.; van Ruitenbeek, J. M.; Thijssen, J. M.; van der Zant, H. S. J.; Duli´c, D. Large tunable image-charge effects in single-molecule junctions. Nature Nanotechnology 2013, 8, 282 EP –. (3) Heimel, G.; Romaner, L.; Br´edas, J.-L.; Zojer, E. Interface Energetics and Level Alignment at Covalent Metal-Molecule Junctions: π-Conjugated Thiols on Gold. Phys. Rev. Lett. 2006, 96, 196806. (4) Neaton, J. B.; Hybertsen, M. S.; Louie, S. G. Renormalization of Molecular Electronic Levels at Metal-Molecule Interfaces. Phys. Rev. Lett. 2006, 97, 216405. (5) DellAngela, M.; Kladnik, G.; Cossaro, A.; Verdini, A.; Kamenetska, M.; Tamblyn, I.; Quek, S. Y.; Neaton, J. B.; Cvetko, D.; Morgante, A. et al. Relating Energy Level Alignment and Amine-Linked Single Molecule Junction Conductance. Nano Letters 2010, 10, 2470–2474. (6) Strange, M.; Rostgaard, C.; H¨akkinen, H.; Thygesen, K. S. Self-consistent GW calculations of electronic transport in thiol- and amine-linked molecular junctions. Phys. Rev. B 2011, 83, 115108. 17

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(7) Jin, C.; Thygesen, K. S. Dynamical image-charge effect in molecular tunnel junctions: Beyond energy level alignment. Phys. Rev. B 2014, 89, 041102. (8) Song, H.; Kim, Y.; Jang, Y. H.; Jeong, H.; Reed, M. A.; Lee, T. Observation of molecular orbital gating. Nature 2009, 462, 1039 EP –. (9) Prins, F.; Barreiro, A.; Ruitenberg, J. W.; Seldenthuis, J. S.; Aliaga-Alcalde, N.; Vandersypen, L. M. K.; van der Zant, H. S. J. Room-Temperature Gating of Molecular Junctions Using Few-Layer Graphene Nanogap Electrodes. Nano Letters 2011, 11, 4607–4611. (10) Venkataraman, L.; Klare, J. E.; Tam, I. W.; Nuckolls, C.; Hybertsen, M. S.; Steigerwald, M. L. Single-Molecule Circuits with Well-Defined Molecular Conductance. Nano Letters 2006, 6, 458–462. (11) Kim, T.; Vzquez, H.; Hybertsen, M. S.; Venkataraman, L. Conductance of Molecular Junctions Formed with Silver Electrodes. Nano Letters 2013, 13, 3358–3364. (12) Prins, F.; Barreiro, A.; Ruitenberg, J. W.; Seldenthuis, J. S.; Aliaga-Alcalde, N.; Vandersypen, L. M. K.; van der Zant, H. S. J. Room-Temperature Gating of Molecular Junctions Using Few-Layer Graphene Nanogap Electrodes. Nano Letters 2011, 11, 4607–4611. (13) Wang, G.; Kim, Y.; Choe, M.; Kim, T.-W.; Lee, T. Graphene-Based Molecular Devices: A New Approach for Molecular Electronic Junctions with a Multilayer Graphene Electrode (Adv. Mater. 6/2011). Advanced Materials 2011, 23, 683–683. (14) Jia, C.; Ma, B.; Xin, N.; Guo, X. Carbon Electrode–Molecule Junctions: A Reliable Platform for Molecular Electronics. Accounts of Chemical Research 2015, 48, 2565– 2575.

18

ACS Paragon Plus Environment

Page 18 of 24

Page 19 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(15) Jia, C.; Migliore, A.; Xin, N.; Huang, S.; Wang, J.; Yang, Q.; Wang, S.; Chen, H.; Wang, D.; Feng, B. et al. Covalently bonded single-molecule junctions with stable and reversible photoswitched conductivity. Science 2016, 352, 1443–1445. (16) Gehring, P.; Harzheim, A.; Spi`ece, J.; Sheng, Y.; Rogers, G.; Evangeli, C.; Mishra, A.; Robinson, B. J.; Porfyrakis, K.; Warner, J. H. et al. Field-Effect Control of Graphene– Fullerene Thermoelectric Nanodevices. Nano Letters 2017, 17, 7055–7061. (17) Reddy, P.; Jang, S.-Y.; Segalman, R. A.; Majumdar, A. Thermoelectricity in Molecular Junctions. Science 2007, 315, 1568–1571. (18) Paulsson, M.; Datta, S. Thermoelectric effect in molecular electronics. Phys. Rev. B 2003, 67, 241403. (19) Sadeghi, H.; Sangtarash, S.; Lambert, C. J. Enhanced Thermoelectric Efficiency of Porous Silicene Nanoribbons. Scientific Reports 2015, 5, 9514. (20) Rincon-Garcia, L.; Ismael, A. K.; Evangeli, C.; Grace, I.; Rubio-Bollinger, G.; Porfyrakis, K.; Agrait, N.; Lambert, C. J. Molecular design and control of fullerenebased bi-thermoelectric materials. Nat Mater 2016, 15, 289–293. (21) Vining, C. B. An inconvenient truth about thermoelectrics. Nat Mater 2009, 8, 83–85. (22) Rincon-Garcia, L.; Evangeli, C.; Rubio-Bollinger, G.; Agrait, N. Thermopower measurements in molecular junctions. Chem. Soc. Rev. 2016, 45, 4285–4306. (23) Cui, L.; Miao, R.; Jiang, C.; Meyhofer, E.; Reddy, P. Perspective: Thermal and thermoelectric transport in molecular junctions. The Journal of Chemical Physics 2017, 146, 092201. (24) Noori, M.; Sadeghi, H.; Lambert, C. J. High-performance thermoelectricity in edgeover-edge zinc-porphyrin molecular wires. Nanoscale 2017, 9, 5299–5304.

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(25) Minnich, A. J.; Dresselhaus, M. S.; Ren, Z. F.; Chen, G. Bulk nanostructured thermoelectric materials: current research and future prospects. Energy Environ. Sci. 2009, 2, 466–479. (26) Famili, M.; Grace, I.; Sadeghi, H.; Lambert, C. J. Suppression of Phonon Transport in Molecular Christmas Trees. ChemPhysChem 2017, n/a–n/a. (27) Murphy, P.; Mukerjee, S.; Moore, J. Optimal thermoelectric figure of merit of a molecular junction. Phys. Rev. B 2008, 78, 161406. (28) Kirˇsanskas, G.; Li, Q.; Flensberg, K.; Solomon, G. C.; Leijnse, M. Designing -stacked molecular structures to control heat transport through molecular junctions. Applied Physics Letters 2014, 105, 233102. (29) Li, Q.; Strange, M.; Duchemin, I.; Donadio, D.; Solomon, G. C. A Strategy to Suppress Phonon Transport in Molecular Junctions Using -Stacked Systems. The Journal of Physical Chemistry C 2017, 121, 7175–7182. (30) Lee, K.; Murray, E. D.; Kong, L.; Lundqvist, B. I.; Langreth, D. C. Higher-accuracy van der Waals density functional. Phys. Rev. B 2010, 82, 081101. (31) Enkovaara, J.; Rostgaard, C.; Mortensen, J. J.; Chen, J.; Duak, M.; Ferrighi, L.; Gavnholt, J.; Glinsvad, C.; Haikola, V.; Hansen, H. A. et al. Electronic structure calculations with GPAW: a real-space implementation of the projector augmented-wave method. Journal of Physics: Condensed Matter 2010, 22, 253202. (32) Bjrk, J.; Hanke, F.; Palma, C.-A.; Samori, P.; Cecchini, M.; Persson, M. Adsorption of Aromatic and Anti-Aromatic Systems on Graphene through Stacking. The Journal of Physical Chemistry Letters 2010, 1, 3407–3412. (33) Strange, M.; Kristensen, I. S.; Thygesen, K. S.; Jacobsen, K. W. Benchmark den-

20

ACS Paragon Plus Environment

Page 20 of 24

Page 21 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

sity functional theory calculations for nanoscale conductance. The Journal of Chemical Physics 2008, 128, 114714. (34) Meir, Y.; Wingreen, N. S. Landauer formula for the current through an interacting electron region. Phys. Rev. Lett. 1992, 68, 2512–2515. (35) Thygesen, K. S. Electron transport through an interacting region: The case of a nonorthogonal basis set. Phys. Rev. B 2006, 73, 035309. (36) Larsen, A. H.; Vanin, M.; Mortensen, J. J.; Thygesen, K. S.; Jacobsen, K. W. Localized atomic basis set in the projector augmented wave method. Phys. Rev. B 2009, 80, 195112. (37) Brandbyge, M.; Mozos, J.-L.; Ordej´on, P.; Taylor, J.; Stokbro, K. Density-functional method for nonequilibrium electron transport. Phys. Rev. B 2002, 65, 165401. (38) Soler, J. M.; Artacho, E.; Gale, J. D.; Garca, A.; Junquera, J.; Ordejn, P.; SnchezPortal, D. The SIESTA method for ab initio order- N materials simulation. Journal of Physics: Condensed Matter 2002, 14, 2745. (39) Markussen, T.; Jauho, A.-P.; Brandbyge, M. Surface-Decorated Silicon Nanowires: A Route to High-ZT Thermoelectrics. Phys. Rev. Lett. 2009, 103, 055502. (40) Dubi, Y.; Di Ventra, M. Colloquium : Heat flow and thermoelectricity in atomic and molecular junctions. Rev. Mod. Phys. 2011, 83, 131–155. (41) Esfarjani, K.; Zebarjadi, M.; Kawazoe, Y. Thermoelectric properties of a nanocontact made of two-capped single-wall carbon nanotubes calculated within the tight-binding approximation. Phys. Rev. B 2006, 73, 085406. (42) Rostgaard, C.; Jacobsen, K. W.; Thygesen, K. S. Fully self-consistent GW calculations for molecules. Phys. Rev. B 2010, 81, 085103.

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(43) Kappera, R.; Voiry, D.; Yalcin, S. E.; Branch, B.; Gupta, G.; Mohite, A. D.; Chhowalla, M. Phase-engineered low-resistance contacts for ultrathin MoS2 transistors. Nat Mater 2014, 13, 1128–1134. (44) Rasmussen, F. A.; Thygesen, K. S. Computational 2D Materials Database: Electronic Structure of Transition-Metal Dichalcogenides and Oxides. The Journal of Physical Chemistry C 2015, 119, 13169–13183. (45) Voiry, D.; Goswami, A.; Kappera, R.; de Carvalho Castro e, S.; Kaplan, D.; Fujita, T.; Chen, M.; Asefa, T.; Chhowalla, M. Covalent functionalization of monolayered transition metal dichalcogenides by phase engineering. Nat Chem 2015, 7, 45–49. (46) Voiry, D.; Mohite, A.; Chhowalla, M. Phase engineering of transition metal dichalcogenides. Chem. Soc. Rev. 2015, 44, 2702–2712. (47) Li, H.; Tsai, C.; Koh, A. L.; Cai, L.; Contryman, A. W.; Fragapane, A. H.; Zhao, J.; Han, H. S.; Manoharan, H. C.; Abild-Pedersen, F. et al. Activating and optimizing MoS2 basal planes for hydrogen evolution through the formation of strained sulphur vacancies. Nat Mater 2016, 15, 48–53. (48) Skriver, H. L.; Rosengaard, N. M. Surface energy and work function of elemental metals. Phys. Rev. B 1992, 46, 7157–7168. (49) Computational Materials Repository, https://cmr.fysik.dtu.dk/c2dm/c2dm.html, accessed 23 May 2018. (50) Cahen, D.; Kahn, A.; Umbach, E. Energetics of molecular interfaces. Materials Today 2005, 8, 32 – 41. (51) Jin, C.; Rasmussen, F. A.; Thygesen, K. S. Tuning the Schottky Barrier at the Graphene/MoS2 Interface by Electron Doping: Density Functional Theory and ManyBody Calculations. The Journal of Physical Chemistry C 2015, 119, 19928–19933. 22

ACS Paragon Plus Environment

Page 22 of 24

Page 23 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(52) Franz, R.; Wiedemann, G. Ueber die Wrme-Leitungsfhigkeit der Metalle. Annalen der Physik 1853, 165, 497–531. (53) Cui, L.; Jeong, W.; Hur, S.; Matt, M.; Kl¨ockner, J. C.; Pauly, F.; Nielaba, P.; Cuevas, J. C.; Meyhofer, E.; Reddy, P. Quantized thermal transport in single-atom junctions. Science 2017, 355, 1192–1195.

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Graphic

24

ACS Paragon Plus Environment

Page 24 of 24