Supporting information Fermi Level Pinning at Electrical Metal Contacts of Monolayer Molybdenum Dichalcogenides Changsik Kim 1,, Inyong Moon 1,, Daeyeong Lee 1, Min Sup Choi 1, Faisal Ahmed 1,2, Seunggeol Nam 3, Yeonchoo Cho 3, Hyeon-Jin Shin 3, Seongjun Park 3, and Won Jong Yoo 1, * 1 Samsung-SKKU Graphene Center (SSGC), SKKU Advanced Institute of Nano-Technology (SAINT), 2 School of Mechanical Engineering, Sungkyunkwan University, 2066, Seobu-ro, Jangan-gu, Suwon, Gyeonggi-do, 16419, Korea, 3 Device & System Research Center, Samsung Advanced Institute of Technology (SAIT), 130 Samsung-ro, Yeongtong-gu, Suwon, Gyeonggi-do 16676, Korea These authors contributed equally to this work. *Corresponding author, yoowj@skku.edu 1
Supporting Information 1. Pinning factor and charge neutrality level for hole SBH. According to the Schottky-Mott rule, p-type SBH (ϕ Bp ) decreases with increasing metal work function (ϕ m ). φ Bp = χ+ g φ m E (1) Here, χ is electron affinity. E g is band gap of semiconductor. The sum of electron affinity and band gap is valence band maximum. If p-type semiconductor wasn t pinned, the slope of p-type SBH vs. metal work function (S = dϕ Bp /dϕ m ) would be 1. In order to adopt the pinning factor and charge neutrality level (CNL, ϕ CNL ) for p-type SBH, we used following equations: φ S ( φ φ ) + ( χ + E g φ = Sφ + b (2) Bp = m CNL CNL ) m For S = 1, p-type SBH can be determined by the Schottky-Mott rule: ϕ Bp = E g + χ ϕ m. For S = 0, p-type SBH can be obtained from ϕ Bp = E g + χ ϕ CNL. Based on this equation, we calculated the pinning factor and the CNL of MoTe 2. 2
Supporting Information 2. Fabrication process Figure S1. (a) Monolayer MoS 2 or MoTe 2 exfoliated onto a 285 nm thick SiO 2 layer-covered Si wafer. (b) 2 µm wide channel defined by plasma etching of SF 6 /O 2. (c) Metal electrodes formed for the TLM measurement. (d) Back- to- back Schottky junction consisting of a reverse Schottky diode as source, a simple resistor as the channel, and a forward Schottky diode as the drain. 3
Supporting Information 3. Material details MoS 2 : Synthetic Molybdenum Disulfide obtained from 2D Semiconductors (http://www.2dsemiconductors.com/synthetic-molybdenum-disulfide-mos2/) 80 70 Bulk MoS 2 60 Intensity (a.u.) 50 40 30 20 10 384 408 0-10 360 380 400 420 440 Raman shift (cm -1 ) Figure S2. (a) Raman spectroscopy: The bulk MoS 2 flakes show two prominent Raman peaks at 384cm -1 (E 2g - in plane-) and 408cm -1 (A 1g out-of-plane) and the FWHM (full-widthat-half-maximum) is less than 4-5cm -1 showing highly crystalline nature. PL intensity (a.u.) 20 18 16 14 12 10 8 6 4 2 0 1.85-2 1.2 1.4 1.6 1.8 2.0 2.2 Photon energy (ev) Mono MoS2 Figure S2. (b) Photoluminescence (PL): For the monolayer of MoS 2, molybdenum disulfide possess direct band-gap is observed at 1.85eV. 4
MoTe 2 : MoTe 2 (Molybdenum Ditelluride) obtained from HQgraphene (http://hqgraphene.com/periodictableelements/mo.php) Figure S2. (c) Raman spectroscopy: The A 1g mode peak is at 171 cm -1 and E 1 2g mode peak is at 234 cm -1, for bulk MoTe 2 flakes. 5
Supporting Information 4. Extraction of Schottky barrier heights Derivation of 2D and 3D thermionic emission equations Thermionic emission: Carrier concentration: Fermi-Dirac distribution function: 3D density of states & thermionic emission: = 2 exp = =exp @ E > E C, E C E f >>kt = 4 2 h =2 2 h exp = exp exp, 2D density of states & thermionic emission: = 4 h = 4 h = 4 h exp = exp exp, = 8 Based on electron confinement of 2D structure, thermionic emission equation can be modified by the energy independent density of states. The current I defined by 2D thermionic emission equation employs the reduced power law from T 2 to T 3/2. The SBH for 2D transport contact is extracted from the slope of a linear fit to ln (I/T 3/2 ) as a function of 1/T. The SBH difference between 3D and 2D thermionic emission is ~ 0.01eV. h 6
Figure S3. (a) Back-to-back Schottky junction consisting of a reverse Schottky diode at the source, a simple resistor at the channel, and a forward Schottky diode at the drain. (b) Schematic plot of the thermionic emission properties, illustrating the dependence of SBH on gate voltage and temperature. Inset diagrams describe carrier transport mechanisms depending on temperature and VG. When applied VD is high, drain side tilted down. The end point of linearity is the true SBH at VFB. Under VFB, effective SBH increases linearly due to VG control. Note that unlike conventional silicon-based devices, in which conventional Schottky diodes are controlled by the Schottky barrier on only one side of the device, semiconductor devices fabricated using 2D materials cannot act as efficient Schottky diodes because both the source and drain metal electrodes create Schottky contacts, as illustrated in schematic diagram of figure S3a. Figure S3b plots the SBH and electron transport induced by applying a gate voltage (V G ) bias, assuming the presence of thermal energy. V G and the temperature affected carrier transport, which proceeded mainly through tunneling or thermionic emission mechanisms at 7
the interface between the metal and the 2D material. The thermal effects were highlighted by measuring the carrier transport in the high-temperature range 300 470 K, under which conditions thermionic emission occurs more actively than at low temperatures of 170 300 K. 8
Supporting Information 5-1. SBH and energy band values of MoS 2 Metal [ev] e-sbh [ev] #1 #2 #3 Average Ti 4.3 0.18 0.25 0.26 0.23 Cr 4.5 0.11 0.13 0.14 0.13 Au 5.2 0.35 0.29 0.32 0.32 Pd 5.6 0.37 0.23 0.29 0.30 Pinning factor 0.11 Table S1. SBH of MoS 2 CB [ev] 4.28 S 0.11 E G [ev] 1.88 b -0.271 VB [ev] 6.16 CNL [ev] 4.479 Table S2. Energy levels, pinning factor, y-intercept, and CNL of MoS 2 We fabricated 3 devices with different metal contacts and used average value for analyzing Fermi level pinning, as shown in table S1. 9
Supporting Information 5-2 SBH and energy band values of MoTe 2 Metal Low T range High T range h-sbh [ev] e-sbh [ev] h-sbh [ev] e-sbh [ev] Ti 4.3 0.27 0.03 0.15 0.19 Cr 4.5 0.23 - - 0.27 Au 5.2 0.14 0.06 0.47 0.12 Pd 5.6 0.2-0.12 0.05 Pinning factor -0.07 0.033 0.041-0.137 Table S3. SBH of MoTe 2 CB [ev] 3.71 S -0.07 E G [ev] 1.28 b 0.531 VB [ev] 4.99 CNL [ev] 4.772 Table S4. Energy levels, pinning factor, y-intercept, and CNL of MoTe 2 Unlike MoS 2 devices, MoTe 2 devices show ambipolar property with strong p-type characteristics and different temperature dependence on carrier type. The linearity of SBH is also harder to define than MoS 2. This is because thermionic currents of both carriers are overlapped by tunneling current of each carrier. Specially, in MoTe 2 /Au, no end point of linearity is observed at high temperature range for hole SBH due to V G limitation. In Table S4-2, we gathered hole and electron SBH for different temperature range. In order to describe hole Fermi level pinning, we used hole SBH at low temperature range. This indicates that temperature dependent transfer curves seem to have limitation to obtain the SBH for ambipolar materials. 10
Supporting Information 6. Detailed references Method Ref # Thickness Metal work function SBH Pinning factor Reference Transistor temperature dependence Diode Temperature Transistor Vd intercept bulk Sc 3.5 0.03 1 bulk Ti 4.2 0.05 bulk Ni 4.9 0.15 0.09 S. Das et al Nano Letters 13(1) 100-105 2013 bulk Pt 5.8 0.23 2 bulk Au 5.1 0.126 bulk Pd 5.4 0.4 0.91 N. Kaushik et al Applied Physics Letters 105(11) 113505 2014 3 3~7 Ni 5 0.26 3~7 Pd 5.6 0.33 0.12 3~7 Ni 5 0.18 3~7 Pd 5.6 0.2 0.03 S. Bhattacharjee et al IEEE Electron Device Letters 37(1) 119-122 2016 4 bulk Mo 4.5 0.14 G. Yoo et al IEEE Electron Device Letters 36(11) 1215-1218 2015 4layer Ti 4.3 0.095 0.11 5 4layer Pt 5.6 0.237 S. Lee et al Nano Letters 16(1) 276-281 2016 4layer Ti 4.3 0.03 6 13nm Co ferro 0.121 13nm Co ferro 0.027 A. Dankert et al ACS Nano 8(1) 476-482 2014 7 58nm Gr 0.022 J. Yoon et al Small 9(19) 3295-3300 2013 8 9 10 11 30nm bulk bulk mono Ti Cr Au Ti 4.3 4.5 5.1 4.33 0.05 0.087 0.3 0.3 30nm bulk bulk mono Ti Cr Gr/AU Ti 4.3 4.5 4.33 0.041 0.08-0.046 0.35 S. Kim et al Applied Physics A 117(2) 761-766 2014 K. Sano et al Japanese Journal of Applied Physics 55(3) 036501 2016 D. Qiu et al Scientific Reports 5 13743 2015 W. Liu et al ACS Nano 9(8) 7904-7912 2015 10~15 Ti 4.2 0.043 12 10~15 Ni 5 0.097 10~15 Au 5.2 0.12 0.21 N. Kaushik et al ACS applied materials & interfaces 8(1) 256-263 2016 10~15 Pd 5.5 0.37 13 Gr 4.5 0.23 Ti 0.4 J. Kwak et al Nano Letters 14(8) 4511-4516 2014 14 mono Ti 4.3 0.072 B. Liu et al ACS Nano 8(5) 5304-5314 2014 tri Py ferro -0.00572 15 mono Py ferro 0.08 W. Wang et al Scientific Reports 4 6928 2014 mono Py ferro 0.0027 CAFM Pt tip 16 50 Pt 5.3 0.3 F. Giannazzo et al Physical Review B 92(8) 081307 2015 Table S5. Experimentally obtained SBH for bulk and monolayer MoS 2 In this Table, the test method, thickness, device structure, and obtained SBH and pinning factor are shown in detail. Referring to this table, in this work we gathered bulk SBH by using the same method that measures drain current at flat band voltage for various temperatures. 11
Method Ref # Thickness Metal work function SBH Pinning factor Reference mono Mo 4.5 0.13 mono Ti 4.3 0.35 17 mono In 4.1 0.47 0.34 J. Kang Physical Review X 4(3) 031005 2014 mono Au 5.3 0.62 mono Pd 5.6 0.9 mono Sc 3.5-0.1 mono Mg 3.7 0.12 mono Al 4.27 0.23 mono Ti 4.33 0.41 mono Cr 4.5 0.45 18 mono Mo 4.6 0.47 mono Ru 4.71 0.43 0.29 Y. Guo et al ACS Applied Materials & Interfaces 7(46) 25709-25715 2015 mono Co 5 0.53 mono Pd 5.13 0.66 mono Ni 5.15 0.54 mono Pt 5.65 0.69 mono MoO3 6.6 0.84 mono Al 4.24 0.51 mono Ag 5.14 0.54 19 mono Au 5.76 0.88 mono Pd 5.83 0.85 0.27 C. Gong et al Nano Letters 14(4) 1717-1720 2014 mono Ir 5.89 0.86 mono Pt 6.12 1.01 mono Sc 3.593 0 mono Sc 3.593 0 mono Ti 4.427 0.187 mono Ti 4.427 0.216 mono Ag 4.489 0.212 0.33 mono Ni 5.222 0.633 mono Au 5.226 0.763 20 mono Pt 5.755 0.52 bi Sc 3.593 0 H. Zhong et al Sceintific Reports 621786 2016 bi Sc 3.593 0 bi Ti 4.427 0.096 bi Ti 4.427 0.161 bi Ag 4.489 0.138 0.27 bi Ni 5.222 0.612 bi Au 5.226 0.667 bi Pt 5.755 0.345 Calculation Table S6. Theoretical SBH of MoS 2 Calculation 18 mono Sc 3.5 0.89 mono Al 4.27 0.76 mono Cr 4.5 0.84 mono Ru 4.71 0.64 mono Co 5 0.66 mono Pd 5.13 0.72 mono Pt 5.65 0.49-0.17 Y. Guo et al ACS Applied Materials & Interfaces 7(46) 25709-25715 2015 Table S7. Theoretical SBH of MoTe 2 This table shows SBHs obtained from theoretical calculation. Interestingly, all the obtained pinning factors are around 0.3. From this table, SBH obtained only from monolayer 2D materials are used to compare with our data. In case of MoTe 2, we used hole SBH instead of electron SBH. 12
Supporting Information 7. The relation between SBH and R c Figure S4. The SBH and R c as function of V G. (a) MoS 2 /Ti contact, (b) MoS 2 /Cr contact, (c) MoS 2 /Au contact, (d) MoS2/Pd contact. The relation between R c and SBH is approximately linear for all the metals, as shown in figure S4. Only for the titanium, the trend is slightly deviating from the linear line, perhaps due to the formation of TiO 2 from the reaction of titanium and ambient oxygen. 13
Supporting Information 8. SBH and R c against V G before and after chemical doping Figure S5. The transfer curves and R c -V G, SBH-V G are plotted with before and after doping. (a), (c) are MoS 2 and (b), (d) are MoTe 2. Here, these chemical doping techniques, based on surface charge transfer dopants, can reduce R c and change carrier transport at the interface by affecting the area beneath metal contact. Their affected depth is around 1.5~3 nm which is deeper than the thickness of monolayer TMDCs. Figures S4 shows the R c against V G of monolayer MoS 2 and MoTe 2 FETs contacted with Cr. The R c and polarities are controlled well via chemical doping. The pristine MoS 2 (black color) was prepared first, and a 20 mm BV solution was then applied to induce degenerated n + doping (red color). A 20 mm AuCl 3 solution was used to induce a p-type doped state (blue color). We confirmed that a p-type transfer curve was clearly obtained, although the degree of doping was weaker than that obtained from BV due to the pristine 14
MoS 2 was n-type in nature. Figure S4b indicates that R c for the pristine MoS 2 was 10 2 10 10 kωµm across the entire V G range. After BV doping, R c decreased significantly, to 10 kωµm at V G = 40. After AuCl 3 doping R c was measured to be 56 kωµm at V G = 40. MoS 2 doped with BV and AuCl 3 provided low values of R c, but both R c are still high due to thickness dependence. Furthermore, the SBH decreased, as shown in Table S8. The SBH of a pristine MoS 2 -Cr FET was found to be 0.1 ev (SBH for electrons), and this value decreased to 0.01 ev (SBH for holes) after AuCl 3 doping. The properties of a MoTe 2 FET in contact with Cr are shown in figure S4e. After BV doping, pristine p-type curve changed to n-type with high on/off ratio. The R c decreased to 3 4 kωµm after AuCl 3 doping. This large decrease was attributed to a change in the SBH of MoTe 2, which decreased significantly from 0.15 to below 0 ev after AuCl 3 doping. The negative SBH indicated that no physical barrier formed at the interface between MoTe 2 and Cr, thereby providing a low R c of 3 4 kωµm across the entire V G range. V G MoS 2 MoTe 2 Pristine n-doping p-doping Pristine n-doping p-doping Contact resistance [kω μm ] -40 V 56 10 4 3 40 V 10 2 13 10 7 10 5 4 SBH [ev] @V FB 0.1 0.01 0.15 < 0 Table S8. A summary table listing the contact resistance (R c ) and SBH values of the devices. 15