Twist-Angle-Dependent Optoelectronics in a Few ... - ACS Publications

Dec 18, 2018 - ABSTRACT: Lattice matching has been supposed to play an important role in the coupling between two materials in a vertical heterostruct...
0 downloads 0 Views 3MB Size
Subscriber access provided by University of South Dakota

Surfaces, Interfaces, and Applications

Twist-angle Dependent Optoelectronics in a fewlayer Transition-Metal Dichalcogenides Heterostructure Woosuk Choi, Imtisal Akhtar, Malik Abdul Rehman, Minwook Kim, Dongwoon Kang, Jongwan Jung, Yoon Myung, Jungcheol Kim, Hyeonsik Cheong, and Yongho Seo ACS Appl. Mater. Interfaces, Just Accepted Manuscript • DOI: 10.1021/acsami.8b15817 • Publication Date (Web): 18 Dec 2018 Downloaded from http://pubs.acs.org on December 19, 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 25 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

ACS Applied Materials & Interfaces

Twist-angle Dependent Optoelectronics in a few-layer Transition-Metal Dichalcogenides Heterostructure Woosuk Choi1†, Imtisal Akhtar1†, Malik Abdul Rehman1†, Minwook Kim1, Dongwoon Kang1, Jongwan Jung1, Yoon Myung1, Jungcheol Kim3, Hyeonsik Cheong3, Yongho Seo1,2* 1Department

of Nanotechnology and Advanced Material Engineering, Sejong University, Seoul, 05006, Korea 2 Graphene Research Institute and HMC, Sejong University, Seoul, 05006, Korea 3

Department of Physics, Sogang University, Seoul 04107, Korea

†These

authors contributed equally to this work author E-mail: [email protected], Tel.:+82-2-3408-3689 *Corresponding

Abstract Lattice matching has been supposed to play an important role in the coupling between two materials in vertical heterostructure (HS). To investigate this role, we fabricated a heterojunction device with a few layers of p-type WSe2 and n-type MoSe2 with different crystal orientation angles. The crystal orientations of WSe2 and MoSe2 were estimated by using highresolution X-ray diffraction. Heterojunction devices were fabricated with twist angles of 0˚, 15˚ and 30˚. The I-V curve of the sample with the twist angle of 0° under the dark condition showed a diode-like behavior. The strong coupling due to lattice matching caused a wellestablished p-n junction. In cases of 15° and 30° samples, the van der Waals gap was built due to lattice mismatching, which resulted in the formation of a potential barrier. However, when the LED light of 365 nm (3.4 eV) was illuminated, excited electrons and holes were possible to jump beyond the potential barrier and the current flowed well in both forward and reverse directions. The effects of the twist angle were analyzed by spectral responsivity and external quantum efficiency, where it was found that the untwisted HS exhibited higher sensitivity under IR illumination while the twisting effect was not noticeable under UV illumination. From the photoluminescence and Raman spectroscopy, it was confirmed that the twisted HS showed a weak coupling due to the lattice mismatch. Key words: transition-metal dichalcogenides; heterostructure; twist angle; optoelectronics; a few-layer TMD; p-n junction;

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

1. Introduction Recently, transition metal dichalcogenides (TMDs) having band gaps ranging from 1.1 to 1.9 eV have attracted considerable attention as promising two-dimensional (2D) layered materials, which can compensate for the lack of band gap in graphene. In particular, the TMDsbased p-n junction device stands as a key strategy for developing various kinds of applications, such as solar cells, photodetectors, light emitting diodes (LEDs) and photovoltaics.1-6 In order to improve the performance of those devices, diverse combinations of 2D hexagonal structure materials, such as graphene, hexagonal boron nitride (h-BN) and TMDs, have been suggested and studied. For instance, h-BN/graphene7, WSe2/graphene8, WSe2/MoS29, WS2/MoS210 and graphene/WSe2/MoS211 van der Waals (vdW) heterostructures (HSs) with vertical stacking, have been demonstrated. However, it has been difficult to match the lattices between two different 2D materials, in fabricating these HS devices. It was reported that weak vdW interaction between the 2D materials in HS fabricated by an artificial transfer induces an interlayer distance with mismatched lattices.12 A tunneling photocurrent depending on coupling strength was studied at the interface in MoS2/MoSe2 heterojunction devices by Luong, et al.13 In an effort to solve the lattice mismatch, many groups have fabricated HSs such as WS2-MoS214, WSe2/MoSe22, 15 and WTe2-MoTe216 with similar lattice parameters. Another approach has been attempted to match lattices by aligning the crystal orientation.17-18 In particular, the vertical HSs directly grown by chemical vapor deposition (CVD) have a strong coupling due to the perfect matching of the crystal directions between the respective materials, and exhibits higher performance than HSs fabricated by the transfer method.14 Heo et al. reported that the band structure in MoS2/WS2 monolayer stacks can be controlled from indirect- to direct-gap depending on the twist angle, and experimentally confirmed this via heteroepitaxial stacking growth.18 Unfortunately, the grown layers may have a lot of defects in general, the twist angle cannot be adjusted intentionally in the growth process,19-20 and the synthesis of the specific layer at specific location is not controlled. 21-22 To fabricate samples with specific twist angles, the crystal orientation of the single crystal must be determined. Transmission electron microscopy (TEM)18, 23-24 and second harmonic generation (SHG)25-26 were known as methods for determining the crystal orientation. TEM is capable of analyzing the crystal orientation and composition of various materials accurately. However, TEM analysis has a disadvantage as the sample cannot be inspected without

ACS Paragon Plus Environment

Page 2 of 25

Page 3 of 25 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

ACS Applied Materials & Interfaces

destroying the device. The SHG method is adequate to estimate the crystal orientation for single layer materials. However, it is known that the resolution of the crystal orientation angle is low and second harmonic power for the multilayer is much lower than that for a single layer, which makes the analysis difficult.27-29 On the contrary, the high-resolution X-ray diffraction (HRXRD) method had advantages of simple measurement, high angle resolution (0.1˚ resolution), and non-destructive inspection. In this paper, we used HR-XRD to find the single crystal direction of p-type WSe2 and n-type MoSe2 and the fabricated few-layer HS devices with twist angles of 0˚, 15˚ and 30˚. By analyzing the fabricated devices, we investigated the effect of the twist angle and how it differs from HS composed of monolayers. Through this analysis, we have provided evidence of how important the lattice matching is in vertically coupled HS devices.

2. Experimental Section 2.1 Sample fabrication mechanism A special method was developed to fabricate HS devices with specific twist-angles. Starting from a single crystal flake (~1 mm, purchased from HQ Graphene), micro-mechanical exfoliation using Scotch tape was performed to obtain a few-layer WSe2 or MoSe2, instead of chemical intercalation.30-33 Scotch tape with exfoliated flakes was mounted on the HR-XRD to determine the crystal orientation. In the repeated exfoliation process, the Scotch tape was moved in translational motions only, without twisting to keep the single crystalline orientation for all flakes. After determining the crystal orientation, the Scotch tape was attached to a SiO2 (300 nm)/Si substrate to transfer a few-layers of WSe2 (or MoSe2) on it. The thin flakes of WSe2 or MoSe2 on the Si substrate were found under an optical microscope, and an atomic force microscope (AFM) was used to select smooth and clean flakes. Thereafter, we picked up the few-layer MoSe2 on a polypropylene carbonate (PPC) film via a so-called ‘pick-up’ method.34 The few-layer WSe2 on the other substrate was picked up in the same way after adjusting the twist angles (0˚, 15˚ and 30˚) on the same PPC with MoSe2. (For details, see Supporting data Fig. S1) The twist angle for a representative sample was confirmed by second harmonic generation, as shown in the Supporting Information Fig. S9. Afterwards, they were transferred onto the other Si substrate and it was placed in a chloroform solution for 20 minutes to remove the PPC film and rinsed with isopropanol. For electrical

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

characterization of an HS device, a polymethyl methacrylate was coated as a resist, and e-beam lithography was performed to draw electrode patterns. Finally, e-beam evaporation was used to deposit Pd(20 nm)/Au(40 nm) on WSe2 and Ti(20 nm)/Au(40 nm) on MoSe2. In order to reduce the contact resistance avoiding the Schottky barrier, Pd was deposited on the p-type WSe2 while Ti was deposited on the n-type MoSe2.35-36 Surface oxidation or contamination was not found in the fabrication, and absence of 1T phase was confirmed by x-ray photoelectron spectroscopy, as shown in Supporting Information Fig. S8.

2.2 Device measurements (HR-XRD, AFM, Raman, PL, I-V curve) High-resolution X-ray diffraction (PANalytical X'pert PRO) was used to confirm the crystal orientation of WSe2 and MoSe2. While 2θ was fixed at 37.8 ˚ corresponding to (103) plane, two other angles, Ψ and Φ, (defined in the Supporting data Fig. S2) were scanned independently. The scan size of Ψ was 10˚ to 80˚ and the step size was 3˚, and the scan size of Φ was 0˚ to 360˚ and step size was 1˚. Ψ was fixed first, Φ was rotated from 0˚ to 360˚, and Ψ was moved to the next step. For an accurate thickness measurement of the samples, commercial AFMs (XE-100 by Park systems Inc., n-Tracer by Nano Focus Inc.) were used in contact mode with a scan rate of 0.3 Hz. Photoluminescence (PL) was measured using a micro Raman system (Nanofinder 30) with an argon laser (532 nm, 0.3 mW) under ambient conditions. The Raman spectrum was measured using another micro Raman system (Renishaw, InVia systems) with a laser line of 514.5 nm (argon ion) under ambient conditions. For comparison, PL and Raman were measured at three different regions: WSe2 only, MoSe2 only, and the overlapped area in the heterostructure devices. In order to investigate electrical and optoelectronic properties, I-V curves and transfer curves were measured under different LED lights (365 nm at 36 mW/cm2, 530 nm at 30.5 mW/cm2, and 850 nm at 102.39 mW/cm2) illumination. All electrical measurements were performed at room temperature under vacuum (~ 10 mTorr) conditions.

3. Results and discussion 3.1 HR-XRD for crystal alignment To investigate the effect of the crystal twist angle, the crystal directions of the few-layers WSe2 and MoSe2 were determined before fabricating HS with the twist angles. In both cases

ACS Paragon Plus Environment

Page 4 of 25

Page 5 of 25 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

ACS Applied Materials & Interfaces

of WSe2 and MoSe2, the angle 2θ = 37.8˚ corresponds to crystal plane (103) with a strong peak,37 and the incident beam and detector were fixed at this angle. The crystal flakes on the pieces of Scotch tape were measured, and the pole figures as functions of (Ψ, Φ) are shown in Figure 1(a). Six peaks appeared because WSe2 and MoSe2 form a hexagonal lattice structure. At a fixed Ψ=58˚ value with the highest peaks in the data, the intensity as a function of Φ is shown in Figure 1(b). The peaks in the line profile represent the crystal direction with the spacing of 60˚ between the peaks. As these results indicate sharp peaks, it was confirmed that the crystal orientation of flakes were aligned well. From these line profiles, the tilt angles of WSe2 and MoSe2 were estimated as 26.3˚ and 7.5˚, respectively. Based on these tilt angles, the vertical p-n junctions of WSe2 and MoSe2 with twist angles of 0˚, 15˚, and 30˚ were fabricated as illustrated in Fig. 1(c). The expected lattice structures of WSe2 and MoSe2 with twist angles 0˚, 15˚ and 30˚ as shown in Fig. 1(d).

Figure 1 High-resolution x-ray diffraction data of WSe2 and MoSe2 and schematics of WSe2/MoSe2 heterostructure p-n junction. (a) HR-XRD results of WSe2 and MoSe2 are plotted in the pole figure. (b) The line profile of the Φ scan is plotted at fixed Ψ with the highest intensity. (c) Schematic drawing of the WSe2-MoSe2 device. WSe2-MoSe2 heterostructure p-n junction on SiO2 substrate was illuminated by LED lights (365 nm, 560 nm, 850 nm) in a vacuum at room temperature. (d) Illustration of crystal orientation with different twist angles (0˚, 15˚ and 30˚).

As a preliminary test, field-effect transistors (FET) with WSe2 and MoSe2, separately were fabricated, and transfer curves and the optoelectronic characteristics were measured to check the carrier types. (Figure S3-4) These devices were made of a few layers of WSe2 and MoSe2,

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

which were confirmed by the Raman spectra. The I-V characteristics of WSe2 showed that the dark current was low. (Figure S3(c)). When the LED light was illuminated on the WSe2 device, the current increases remarkably as the photon energy was increased with short wavelength. The transfer curves depending on the gate voltage revealed that the current raised dramatically as the gate voltage became negative as shown in Fig. S3(d). This confirmed that WSe2 was a p-type semiconductor. On the other hand, when the gate voltage was less than -40 V, the current was reduced as the shorter wavelength LED was illuminated. This implies that the holes generated by the field effect were recombined with the electrons excited by high energy photon illumination. This unexpected behavior has been observed by other groups, which was explained by the dissociation of adsorbed molecules on the TMD surface.38 The on/off ratio depending on the gate voltage in our sample was 106, which is similar to the monolayer WSe2 based FET as reported previously.35 In the case of MoSe2 devices, Ti was selected as the contact metal that lowered the Schottky barrier, and it was confirmed that the current amounted to 74 nA, when 1 V was applied in a dark condition. (See Fig. S4) While the current was slightly increased from the dark state at 530 nm and 850 nm illuminations, it was increased noticeably only under the 365 nm light illumination (3.4 eV). (Fig. S4(c)) Because the few layer MoSe2 had an indirect band gap, electrons and holes were considerably generated under the high-energy illumination. The current was increased when the gate voltage was positive in the transfer curve, confirming that MoSe2 was an n-type semiconductor (Fig. S4(d)). Contrary to WSe2, the larger current under the 365 nm illumination was due to the increased number of electrons generated by the higher photon energy.

3.2 Microscopic image and I-V curve Figure 2(a) shows an optical micrograph of a typical heterostructure device, where MoSe2 is on the right, and WSe2 on the top side. In the middle of WSe2 and MoSe2, the light blue area was the vertically overlapped heterostructure, and the yellow areas were the electrodes at both ends. Figures 2(b-c) are images of the surface morphology of the heterostructure obtained by AFM. The thicknesses of WSe2 and MoSe2 were measured around 3 nm from the line profile shown in Fig. 2(d). All the HS devices were made of WSe2 and MoSe2 with this thickness corresponding to 4 layers, as the thickness of a single layer is 0.7 nm.

ACS Paragon Plus Environment

Page 6 of 25

Page 7 of 25 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

ACS Applied Materials & Interfaces

Figure 2 Microscopic analysis of the WSe2/MoSe2 heterostructure p-n junction. (a) The optical micrograph of the WSe2/MoSe2 device. (b) 3-dimensional and (c) plane images of the AFM surface morphology of HS. (d) Line profiles of WSe2 and MoSe2 were plotted in red and blue lines, respectively

Figures 3(a, b) show the I-V curves of representative samples with different twist angles. The I-V curves in the dark state of the 0° sample showed a diode-like behavior of the p-n junction, and when the voltage was 1 V, the current was 0.62 nA as shown in Fig. 3(a), which was larger than the p-n junctions of the monolayer.11, 39 As a result, it was confirmed that the p-type WSe2 and the n-type MoSe2 were coupled to form a p-n junction with a built-in potential. (i) The band diagrams of the p-n junction without bias are shown in Fig. 3(c). A Kelvin probe measurement revealed that the work function for multilayer MoSe2 is lower than that for multilayer WSe2. (Supporting Information Fig. S7) As the Fermi level of p-type WSe2 is lower than that of n-type MoSe2, the energy bands for WSe2 is shifted upward due to the band bending. (ii) When the forward bias was applied, the potential barrier decreased, and the electrons of MoSe2 in the conduction band were easily transferred to WSe2. Also, the holes of WSe2 in the

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

valence band can be easily transferred to MoSe2 under the dark state. (iii) On the contrary, when a reveres bias was applied, a large potential barrier was formed, and the current was reduced. Under the LED illumination, the current of the 0° sample was increased in the forward direction more than that in the dark state as the electron and hole pairs are generated by the photon energy. In Fig 3(a), (ii) the larger photocurrents were induced under shorter wavelength (365 and 560 nm) illuminations. As the shorter wavelength photons have higher energy, the more electrons in the valence band can be excited to the conduction band. (iii) On the other hand, the potential barrier became higher in the reverse bias, and almost no current flowed in the dark condition. When it was illuminated by high-energy photons, the limited numbers of electrons and holes were generated near the junction, and a low current flowed through the potential barrier.

Figure 3 I-V curves of heterostructures with different twist angles ((a) 0˚, (b) 15˚) under LED illumination with different wavelengths: 365 nm, 560 nm, and 850 nm. Schematic energy band diagrams are shown for the (c) p-n junction HS with a staggered band alignment and (d) HS with a vdW gap at different bias voltages: (i) zero bias, (ii) forward bias, and (iii) reverse bias. Electrons (-) and holes (+) in the dark state were indicated as white symbols. Red symbols and dotted lines represent the electron and hole generation in the 850 nm LED illumination, while green in 530 nm and blue in 365 nm.

The I-V curves of the HSs with twist angles 15˚ and 30˚ shown in Figure 3(b) and S5

ACS Paragon Plus Environment

Page 8 of 25

Page 9 of 25 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

ACS Applied Materials & Interfaces

exhibited a completely different behavior, compared to the 0° sample. In the case of two twisted samples under the dark condition, the photocurrent hardly flowed in the order of ~10 pA. (Fig. S5(c)) It was reported that the interlayer distance at the interface of HS depends on the twist angle.17 While the vdW HS has a high-symmetry-stacking with the strongest interlayer coupling strength for twist angle of 0° or 60°, HSs with the other angles have weak coupling strength due to mismatched atomic alignment. The coupling strength also can be changed by thermal annealing.13 We assume that the weak coupling expands the vdW gap with a potential barrier, and the p-n junction may not be formed normally. 12, 40 Thus, a potential barrier at the interface blocked the carriers under the dark condition, and the conductance was decreased as shown in Fig. 3(b). However, the photocurrent appeared noticeably only under UV (365 nm) illumination. Due to its high energy, the electrons and holes generated by a 365 nm photon can go beyond the potential barrier formed by the vdW gap as shown in Fig. 3(d).40 Therefore, the I-V curves under the UV illumination showed linear behavior according to the band diagram without band bending. In Figure S5(a), the current of the 30˚ sample was even lower than the 15˚ sample, due to the higher potential barrier being induced as the interlayer distance was increased.17

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

Figure 4 A comparison of (a-c) spectral responsivity and external quantum efficiency with different the twist angle under the LED illuminations ((a, d, g) 365 nm, (b, e, h) 560 nm, (c, f, i) 365 nm). The responsivity and EQE data in (d-f) were obtained at a fixed bias of 1 V, while EQE data in (g-i) were obtained at 0 V.

Spectral responsivity (R) and external quantum efficiency (EQE) were estimated to compare the twist angle effects with different wavelength illuminations. Spectral responsivity is defined as the photocurrent ratio to the intensity of the illuminated light. It can be expressed by:

R=

𝐼𝑝ℎ 𝑃0 × 𝑆

,

where Iph is the difference between the illuminated (Iillum) and dark (Idark) current i.e. (Iph = Iillum – Idark), P0 is the illumination power density, and S is the effective illuminated area. The spectral responsivity and external quantum efficiency vary depending on the wavelength of light, the intensity of light, bias voltage, and gate voltage (Vg). The intensity of light was fixed as described in the experimental section. R was calculated under bias 1 V and Vg = 0 V. Figure

ACS Paragon Plus Environment

Page 10 of 25

Page 11 of 25 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

ACS Applied Materials & Interfaces

4(a) shows the responsivity in IR (850 nm) illumination depending on the twist angle. Four samples for each twist angle were prepared and measured individually, and the median values are plotted together. In particular, it had lower responsivity at 15˚ and 30˚, which implied that the junction was not well formed due to the lattice mismatching. According to Nayak et al. the lattice parameter and interlayer distance of the HS lattice vary depending on the twist angle in the range from 0˚ to 60˚ with the maximum distance at 30˚.17 In Heo et al.’s paper, the band structure changes from the direct bandgap to the indirect band gap as the interlayer distance increases.18 It was plausible that the twisted samples have lower spectral responsivities, but under UV (365 nm) illumination there was no consistent tendency except the experimental error, as the photons had high energy above the potential barrier. (Fig. 4(c)) The large experimental error was attributed to the delicacy in fabrication of the HS device. Responsivities of photovoltaic devices based on monolayer TMDs reported by others were in large variations from 0.1 mA/W to 880 A/W, in deed.41-43 EQE can be estimated from the responsivity simply by: EQE =

𝑅 × ℎ𝑐 𝑒λ

,

where h is the Planck’s constant, c the speed of light, e the charge of the electron, and λ the wavelength of the light illuminated. EQE under IR illumination was around a few percent, as the light has a low-energy, but the untwisted device (0˚) exhibited a remarkably higher value compared with the other twisted devices as shown in Fig. 4(d). EQE value of the untwisted device under visible illumination (530 nm) was around 100% as shown in Fig. 4(e), which is much higher than that reported by Lee et al.11 Indeed, EQE can be greater than 100%, especially under a bias, because the applied voltage provides carriers, and excitation of photons increases the conductivity of the junction. Also, the remarkable increment of EQE from IR to visible illumination is in good agreement with the paper.11 The data of EQE measured at zero bias voltage are shown in Fig. 4(g-i), where the angle dependency trend is similar to the EQE at 1 V bias, though the absolute values are much lower. Under UV illumination, the twist angle dependence of EQE was weakened due to the same reason for the responsivity as shown in Fig. 4(f).

3.3 Photoluminescence and Raman Spectroscopy To investigate the excitons in the heterojunctions with different twist angles between

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

MoSe2 and WSe2, we measured the PL at room temperature. In the Fig 5(a), two peaks (peak A at 1.5 eV and peak I at 1.7 eV) in the WSe2 spectrum indicate that the number of layers was around four.29, 44 The peak A representing the direct band gap was weakened compared with the monolayer, and peak I corresponding the indirect band gap was lowered compared with bulk state. These results were consistent with the number of layers obtained from AFM data. MoSe2 spectrum showed a peak near 1.6 eV, corresponding to the K-K direct band gap. This PL peak around 1.5–1.6 eV exhibits its maximum intensity for a single-layer MoSe2, and the few-layer MoSe2 show gradual decrement with the number of layers.45 In the PL spectrum of the MoSe2/WSe2 heterostructures, the peak for MoSe2 was clear, but WSe2 peaks were suppressed. The reason why the MoSe2 peak was much stronger than the WSe2 peaks is supposed to be due to the geometry having MoSe2 on the top surface.

Figure 5 Photoluminescence and Raman spectra of MoSe2/WSe2 vertical heterostructures. (a) PL

ACS Paragon Plus Environment

Page 12 of 25

Page 13 of 25 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

ACS Applied Materials & Interfaces

spectra of individual few layers WSe2 (blue), MoSe2 (red) and HS (black). (b) PL spectra of HSs with different twist angles (0˚, 15˚ and 30˚). (c) Raman spectra of individual few layers WSe2 (blue), MoSe2 (red), and HS (black). (d) Raman spectra of HSs with different twist angles: 0˚(orange), 15˚(green) and 30˚(purple).

Figure 5 (b) shows the PL spectra for HSs depending on the twist angle. The weak peak near 1.68 eV observed in all angles is supposed to be red-shifted from peak I of WSe2, due to interlayer radiative recombination.23 Another red-shifted PL signal (1.51 eV) from peak A of the WSe2 was clear particularly in the 30˚ sample. In the 30˚ sample of which the interlayer distance was considered to be large, WSe2 and MoSe2 were less influenced by each other, and peak A of WSe2 was strong. A similar result for twist angle 29˚ was reported by Nayak et al.17 On the other hand, an interlayer exciton peak at 1.3 eV for monolayer based WSe2/MoSe2 HS was reported by others,17, 46-48 which was not found in our sample. The peak at 1.3 eV has a low intensity even in monolayer based HS, and it is only noticeable for the samples with the twist angles of 0˚ and 60˚.17 The spatially separated interlayer excitons were not expected in HS composed of a few layers, based on our experimental results. Raman spectroscopy was adopted to investigate the phonon oscillation of the lattice of the WSe2/MoSe2 heterostructures. The Raman shift in high frequency range from 200 cm-1 to 420 cm-1 was measured including the in-plane and out-of-plane modes of WSe2 and MoSe2. The Raman spectra of the WSe2, MoSe2, and MoSe2/WSe2 heterojunction are shown in Fig. 5(c). The intensities and positions of E12g (and A1g, 250.6 cm-1) and 2LA (261.1 cm-1) modes in the Raman spectrum of WSe2 correspond to four-layers, roughly.29, 49 This result was similar to the estimated number of layers from the thickness measured by AFM. The Raman spectra of MoSe2 including A1g mode (243.6 cm-1) and E12g mode (287.1 cm-1) also show a similar number of layers when referred to other papers.45 The Raman spectra of the WSe2/MoSe2 heterojunction was similar to that of MoSe2 with the biggest peak of A1g mode of MoSe2 because MoSe2 was in the upper position. In addition, the E12g and 2LA mode peaks of WSe2 at the bottom of the heterostructures were found. The phonon mode modification at HS may not significantly influence the Raman data, considering the thickness of the multilayer devices. In Figure 5(d), the Raman spectra of the heterostructures were compared according to the twist angle. Due to different intensities of the Raman spectrum for each sample, the average values of the five samples were calculated for each twist angle. It can be seen that the position and intensity of the A1g and E12g modes of MoSe2 in HS were almost the same as those in MoSe2 only. However, the intensity of the E12g and 2LA modes of WSe2 were changed. In

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

particular, both of the peaks were increased in the 30˚ sample. In order to analyze it statistically, the positions and intensities of the peaks in HSs were compared as shown in the Figure S6 with 5 samples for each twist angle. The positions of Raman modes of WSe2 and MoSe2 were similar to those in a monolayer based HS as shown Figure S6(b-c).17 As the variation of the Raman shifts of each mode for different samples were within 2 to 3 cm-1, it can be seen that no significant effect of the twist angle was observed.29, 49 However, when the intensities of the peaks were analyzed in detail, E12g and 2LA modes of WSe2 were strong, particularly in the 30˚ sample. (Figure S6(a)) A similar behavior confirming our result was found in the monolayer heterojunction, where the WSe2 E2g mode peak at 29˚ twist angle was especially high.17 As a result, it was confirmed that both Raman and PL signals of WSe2 were strong at the 30˚ sample, and these results can be attributed to the large vdW gap at the twist angle 30˚.

4. Conclusions We observed that a large-scale photocurrent was generated at the HSs composed of fewlayers TMDs as compared with that of monolayers. When the twist angle was 0˚, it was found that the p-n junction was well formed because the interlayer distance was short. In cases of 15˚ and 30˚, a potential barrier was induced due to the vdW gap, and electrons and holes excited by low energy photon illumination could not tunnel the barrier. However, when the UV was illuminated, the highly excited electrons and holes could jump beyond the potential barrier, and the bipolar current flowed in both forward and reverse biases. It was difficult to distinguish the different Raman and PL peak shifts depending on the twist angle, but it was confirmed that the intensities of WSe2 peaks below MoSe2 layer became clear at the 30˚ sample. This implies that the twisted orientation causes decoupling of the heterojunction. The importance of crystal orientation matching in such vertical heterojunction devices has been recognized, and lattice matching technique can be applied in optoelectronic devices to strengthen junction coupling. For the optoelectronic device application, scalability or high yield is a critical issue, and the chemical intercalation32-33 or CVD deposition50 can be applied instead of mechanical exfoliation for that purpose.

Supporting Information Sample fabrication process of WSe2-MoSe2 heterostructure p-n junction; HR-XRD data of

ACS Paragon Plus Environment

Page 14 of 25

Page 15 of 25 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

ACS Applied Materials & Interfaces

disarranged crystalline flakes; Optoelectronic properties of FETs, I-V curve of heterostructure with different twist angles; Raman signatures of MoSe2/WSe2 vertical heterostructures; second harmonic generation analysis: X-ray photoelectron spectroscopy result; Kelvin probe force microscopy images;

ACKNOWLEDGMENT This research was supported by the Priority Research Centers Program (2010-0020207) and Basic Science Research Program (2017R1A2B4002379) through the National Research Foundation of Korea (NRF) funded by the Ministry of Education. This work was also supported by the World Class 300 Project (S2561932) of the MOTIE, MSS(Korea).

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

Reference (1) Jariwala, D.; Davoyan, A. R.; Wong, J.; Atwater, H. A. Van Der Waals Materials for Atomically-Thin Photovoltaics: Promise and Outlook. Acs Photonics 2017, 4 (12), 2962-2970. (2) Flory, N.; Jain, A.; Bharadwaj, P.; Parzefall, M.; Taniguchi, T.; Watanabe, K.; Novotny, L. A Wse2/Mose2 Heterostructure Photovoltaic Device. Applied Physics Letters 2015, 107 (12), 123106. (3) Withers, F.; Del Pozo-Zamudio, O.; Mishchenko, A.; Rooney, A. P.; Gholinia, A.; Watanabe, K.; Taniguchi, T.; Haigh, S. J.; Geim, A. K.; Tartakovskii, A. I.; Novoselov, K. S. Light-Emitting Diodes by Band-Structure Engineering in Van Der Waals Heterostructures. Nat Mater 2015, 14 (3), 301-306. (4) Li, H. M.; Lee, D.; Qu, D.; Liu, X.; Ryu, J.; Seabaugh, A.; Yoo, W. J. Ultimate Thin Vertical P-N Junction Composed of Two-Dimensional Layered Molybdenum Disulfide. Nat Commun 2015, 6, 6564. (5) Zhang, W. J.; Wang, Q. X.; Chen, Y.; Wang, Z.; Wee, A. T. S. Van Der Waals Stacked 2d Layered Materials for Optoelectronics. 2d Mater 2016, 3 (2), 022001. (6) Yu, W. J.; Vu, Q. A.; Oh, H.; Nam, H. G.; Zhou, H. L.; Cha, S.; Kim, J. Y.; Carvalho, A.; Jeong, M.; Choi, H.; Neto, A. H. C.; Lee, Y. H.; Duan, X. F. Unusually Efficient Photocurrent Extraction in Monolayer Van Der Waals Heterostructure by Tunnelling through Discretized Barriers. Nat Commun 2016, 7, 13278 (7) Dean, C. R.; Young, A. F.; Meric, I.; Lee, C.; Wang, L.; Sorgenfrei, S.; Watanabe, K.; Taniguchi, T.; Kim, P.; Shepard, K. L.; Hone, J. Boron Nitride Substrates for High-Quality Graphene Electronics. Nat

Nanotechnol 2010, 5 (10), 722-726. (8) Kim, K.; Larentis, S.; Fallahazad, B.; Lee, K.; Xue, J.; Dillen, D. C.; Corbet, C. M.; Tutuc, E. Band Alignment in Wse2-Graphene Heterostructures. Acs Nano 2015, 9 (4), 4527-32. (9) Furchi, M. M.; Pospischil, A.; Libisch, F.; Burgdorfer, J.; Mueller, T. Photovoltaic Effect in an Electrically Tunable Van Der Waals Heterojunction. Nano Lett 2014, 14 (8), 4785-4791. (10) Hill, H. M.; Rigosi, A. F.; Rim, K. T.; Flynn, G. W.; Heinz, T. F. Band Alignment in Mos2/Ws2 Transition Metal Dichalcogenide Heterostructures Probed by Scanning Tunneling Microscopy and Spectroscopy. Nano Lett 2016, 16 (8), 4831-4837. (11) Lee, C. H.; Lee, G. H.; van der Zande, A. M.; Chen, W.; Li, Y.; Han, M.; Cui, X.; Arefe, G.; Nuckolls, C.; Heinz, T. F.; Guo, J.; Hone, J.; Kim, P. Atomically Thin P-N Junctions with Van Der Waals Heterointerfaces. Nat Nanotechnol 2014, 9 (9), 676-81. (12) Allain, A.; Kang, J. H.; Banerjee, K.; Kis, A. Electrical Contacts to Two-Dimensional Semiconductors.

Nat Mater 2015, 14 (12), 1195-1205. (13) Luong, D. H.; Lee, H. S.; Neupane, G. P.; Roy, S.; Ghimire, G.; Lee, J. H.; Vu, Q. A.; Lee, Y. H. Tunneling Photocurrent Assisted by Interlayer Excitons in Staggered Van Der Waals Hetero-Bilayers.

Adv Mater 2017, 29 (33), 1701512. (14) Gong, Y. J.; Lin, J. H.; Wang, X. L.; Shi, G.; Lei, S. D.; Lin, Z.; Zou, X. L.; Ye, G. L.; Vajtai, R.; Yakobson, B. I.; Terrones, H.; Terrones, M.; Tay, B. K.; Lou, J.; Pantelides, S. T.; Liu, Z.; Zhou, W.; Ajayan, P. M. Vertical and in-Plane Heterostructures from Ws2/Mos2 Monolayers. Nat Mater 2014, 13 (12), 1135-

ACS Paragon Plus Environment

Page 16 of 25

Page 17 of 25 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

ACS Applied Materials & Interfaces

1142. (15) Gong, Y. J.; Lei, S. D.; Ye, G. L.; Li, B.; He, Y. M.; Keyshar, K.; Zhang, X.; Wang, Q. Z.; Lou, J.; Liu, Z.; Vajtai, R.; Zhou, W.; Ajayan, P. M. Two-Step Growth of Two-Dimensional Wse2/Mose2 Heterostructures. Nano Lett 2015, 15 (9), 6135-6141. (16) Zheng, Z.; Wang, X.; Mi, W. Strain and Electric-Field Tunable Valley States in 2d Van Der Waals Mote2/Wte2 Heterostructures. J Phys Condens Matter 2016, 28 (50), 505003. (17) Nayak, P. K.; Horbatenko, Y.; Ahn, S.; Kim, G.; Lee, J. U.; Ma, K. Y.; Jang, A. R.; Lim, H.; Kim, D.; Ryu, S.; Cheong, H.; Park, N.; Shin, H. S. Probing Evolution of Twist-Angle-Dependent Interlayer Excitons in Mose2/Wse2 Van Der Waals Heterostructures (Vol 11, Pg 4041, 2017). Acs Nano 2017,

11 (9), 9566-9566. (18) Heo, H.; Sung, J. H.; Cha, S.; Jang, B. G.; Kim, J. Y.; Jin, G.; Lee, D.; Ahn, J. H.; Lee, M. J.; Shim, J. H.; Choi, H.; Jo, M. H. Interlayer Orientation-Dependent Light Absorption and Emission in Monolayer Semiconductor Stacks. Nat Commun 2015, 6, 7372. (19) Surrente, A.; Dumcenco, D.; Yang, Z.; Kuc, A.; Jing, Y.; Heine, T.; Kung, Y. C.; Maude, D. K.; Kis, A.; Plochocka, P. Defect Healing and Charge Transfer-Mediated Valley Polarization in Mos2/Mose2/Mos2 Trilayer Van Der Waals Heterostructures. Nano Lett 2017, 17 (7), 4130-4136. (20) Najmaei, S.; Yuan, J. T.; Zhang, J.; Ajayan, P.; Lou, J. Synthesis and Defect Investigation of TwoDimensional Molybdenum Disulfide Atomic Layers. Accounts Chem Res 2015, 48 (1), 31-40. (21) Chen, L.; Liu, B. L.; Abbas, A. N.; Ma, Y. Q.; Fang, X.; Liu, Y. H.; Zhou, C. W. Screw-DislocationDriven Growth of Two-Dimensional Few-Layer and Pyramid-Like Wse2 by Sulfur-Assisted Chemical Vapor Deposition. Acs Nano 2014, 8 (11), 11543-11551. (22) Ling, X.; Lee, Y. H.; Lin, Y. X.; Fang, W. J.; Yu, L. L.; Dresselhaus, M. S.; Kong, J. Role of the Seeding Promoter in Mos2 Growth by Chemical Vapor Deposition. Nano Lett 2014, 14 (2), 464-472. (23) Fang, H.; Battaglia, C.; Carraro, C.; Nemsak, S.; Ozdol, B.; Kang, J. S.; Bechtel, H. A.; Desai, S. B.; Kronast, F.; Unal, A. A.; Conti, G.; Conlon, C.; Palsson, G. K.; Martin, M. C.; Minor, A. M.; Fadley, C. S.; Yablonovitch, E.; Maboudian, R.; Javey, A. Strong Interlayer Coupling in Van Der Waals Heterostructures Built from Single-Layer Chalcogenides. P Natl Acad Sci USA 2014, 111 (17), 61986202. (24) Li, Y. T.; Huang, L.; Li, B.; Wang, X. T.; Zhou, Z. Q.; Li, J. B.; Wei, Z. M. Co-Nucleus 1d/2d Heterostructures with Bi2s3 Nanowire and Mos2 Monolayer: One-Step Growth and Defect-Induced Formation Mechanism. Acs Nano 2016, 10 (9), 8938-8946. (25) Chen, H.; Corboliou, V.; Solntsev, A. S.; Choi, D. Y.; Vincenti, M. A.; de Ceglia, D.; de Angelis, C.; Lu, Y.; Neshev, D. N. Enhanced Second-Harmonic Generation from Two-Dimensional Mose2 on a Silicon Waveguide. Light Sci Appl 2017, 6 (10), e17060. (26) Hu, L.; Wei, D.; Huang, X. Second Harmonic Generation Property of Monolayer Tmdcs and Its Potential Application in Producing Terahertz Radiation. J Chem Phys 2017, 147 (24), 244701. (27) Kumar, N.; Najmaei, S.; Cui, Q. N.; Ceballos, F.; Ajayan, P. M.; Lou, J.; Zhao, H. Second Harmonic Microscopy of Monolayer Mos2. Phys Rev B 2013, 87 (16), 161403. (28) Kim, D. H.; Lim, D. Optical Second-Harmonic Generation in Few-Layer Mose2. J Korean Phys

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

Soc 2015, 66 (5), 816-820. (29) Zeng, H.; Liu, G. B.; Dai, J.; Yan, Y.; Zhu, B.; He, R.; Xie, L.; Xu, S.; Chen, X.; Yao, W.; Cui, X. Optical Signature of Symmetry Variations and Spin-Valley Coupling in Atomically Thin Tungsten Dichalcogenides. Sci Rep 2013, 3, 1608. (30) Mayorga-Martinez, C. C.; Ambrosi, A.; Eng, A. Y. S.; Sofer, Z.; Pumera, M. Transition Metal Dichalcogenides (Mos2, Mose2, Ws2 and Wse2) Exfoliation Technique Has Strong Influence Upon Their Capacitance. Electrochem Commun 2015, 56, 24-28. (31) Tan, S. M.; Sofer, Z.; Luxa, J.; Pumera, M. Aromatic-Exfoliated Transition Metal Dichalcogenides: Implications for Inherent Electrochemistry and Hydrogen Evolution. Acs Catal 2016, 6 (7), 45944607. (32) Luxa, J.; Vosecky, P.; Mazanek, V.; Sedmidubsky, D.; Pumera, M.; Sofer, Z. Cation-Controlled Electrocatalytical Activity of Transition-Metal Disulfides. Acs Catal 2018, 8 (4), 2774-2781. (33) Luxa, J.; Vosecky, P.; Mazinek, V.; Sedmidubsky, D.; Pumera, M.; Lazar, P.; Sofer, Z. Layered Transition-Metal Ditellurides in Electrocatalytic Applications-Contrasting Properties. Acs Catal 2017,

7 (9), 5706-5716. (34) Wang, L.; Meric, I.; Huang, P. Y.; Gao, Q.; Gao, Y.; Tran, H.; Taniguchi, T.; Watanabe, K.; Campos, L. M.; Muller, D. A.; Guo, J.; Kim, P.; Hone, J.; Shepard, K. L.; Dean, C. R. One-Dimensional Electrical Contact to a Two-Dimensional Material. Science 2013, 342 (6158), 614-617. (35) Fang, H.; Chuang, S.; Chang, T. C.; Takei, K.; Takahashi, T.; Javey, A. High-Performance Single Layered Wse2 P-Fets with Chemically Doped Contacts. Nano Lett 2012, 12 (7), 3788-3792. (36) Pudasaini, P. R.; Oyedele, A.; Zhang, C.; Stanford, M. G.; Cross, N.; Wong, A. T.; Hoffman, A. N.; Xiao, K.; Duscher, G.; Mandrus, D. G.; Ward, T. Z.; Rack, P. D. High-Performance Multilayer Wse2 Field-Effect Transistors with Carrier Type Control. Nano Res 2018, 11 (2), 722-730. (37) Wang, J.; Ma, C.; Mu, X.; Cai, W.; Liu, L.; Zhou, X.; Hu, W.; Hu, Y. Construction of Multifunctional Mose2 Hybrid Towards the Simultaneous Improvements in Fire Safety and Mechanical Property of Polymer. J Hazard Mater 2018, 352, 36-46. (38) Son, S. B.; Kim, Y.; Kim, A.; Cho, B.; Hong, W. K. Ultraviolet Wavelength-Dependent Optoelectronic Properties in Two-Dimensional Nbse2-Wse2 Van Der Waals Heterojunction-Based Field-Effect Transistors. Acs Appl Mater Inter 2017, 9 (47), 41537-41545. (39) Tsai, M. L.; Li, M. Y.; Retamal, J. R. D.; Lam, K. T.; Lin, Y. C.; Suenaga, K.; Chen, L. J.; Liang, G.; Li, L. J.; He, J. H. Single Atomically Sharp Lateral Monolayer P-N Heterojunction Solar Cells with Extraordinarily High Power Conversion Efficiency. Adv Mater 2017, 29 (32), 1701168. (40) Yan, X.; Liu, C. S.; Li, C.; Bao, W. Z.; Ding, S. J.; Zhang, D. W.; Zhou, P. Tunable Snse2/Wse2 Heterostructure Tunneling Field Effect Transistor. Small 2017, 13 (34), 1701478. (41) Lopez-Sanchez, O.; Lembke, D.; Kayci, M.; Radenovic, A.; Kis, A. Ultrasensitive Photodetectors Based on Monolayer Mos2. Nat Nanotechnol 2013, 8 (7), 497-501. (42) Baugher, B. W. H.; Churchill, H. O. H.; Yang, Y. F.; Jarillo-Herrero, P. Optoelectronic Devices Based on Electrically Tunable P-N Diodes in a Monolayer Dichalcogenide. Nat Nanotechnol 2014, 9 (4), 262-267.

ACS Paragon Plus Environment

Page 18 of 25

Page 19 of 25 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

ACS Applied Materials & Interfaces

(43) Mak, K. F.; Shan, J. Photonics and Optoelectronics of 2d Semiconductor Transition Metal Dichalcogenides. Nat Photonics 2016, 10 (4), 216-226. (44) Xu, K.; Wang, Z.; Du, X.; Safdar, M.; Jiang, C.; He, J. Atomic-Layer Triangular Wse2 Sheets: Synthesis and Layer-Dependent Photoluminescence Property. Nanotechnology 2013, 24 (46), 465705. (45) Tongay, S.; Zhou, J.; Ataca, C.; Lo, K.; Matthews, T. S.; Li, J. B.; Grossman, J. C.; Wu, J. Q. Thermally Driven Crossover from Indirect toward Direct Bandgap in 2d Semiconductors: Mose2 Versus Mos2.

Nano Lett 2012, 12 (11), 5576-5580. (46) Lin, Y. C.; Ghosh, R. K.; Addou, R.; Lu, N.; Eichfeld, S. M.; Zhu, H.; Li, M. Y.; Peng, X.; Kim, M. J.; Li, L. J.; Wallace, R. M.; Datta, S.; Robinson, J. A. Atomically Thin Resonant Tunnel Diodes Built from Synthetic Van Der Waals Heterostructures. Nat Commun 2015, 6, 7311. (47) Wilson, N. R.; Nguyen, P. V.; Seyler, K.; Rivera, P.; Marsden, A. J.; Laker, Z. P.; Constantinescu, G. C.; Kandyba, V.; Barinov, A.; Hine, N. D.; Xu, X.; Cobden, D. H. Determination of Band Offsets, Hybridization, and Exciton Binding in 2d Semiconductor Heterostructures. Sci Adv 2017, 3 (2), e1601832. (48) Rivera, P.; Schaibley, J. R.; Jones, A. M.; Ross, J. S.; Wu, S.; Aivazian, G.; Klement, P.; Seyler, K.; Clark, G.; Ghimire, N. J.; Yan, J.; Mandrus, D. G.; Yao, W.; Xu, X. Observation of Long-Lived Interlayer Excitons in Monolayer Mose2-Wse2 Heterostructures. Nat Commun 2015, 6, 6242. (49) Zhao, W. J.; Ghorannevis, Z.; Amara, K. K.; Pang, J. R.; Toh, M.; Zhang, X.; Kloc, C.; Tan, P. H.; Eda, G. Lattice Dynamics in Mono- and Few-Layer Sheets of Ws2 and Wse2. Nanoscale 2013, 5 (20), 9677-9683. (50) Shehzad, M. A.; Hussain, S.; Lee, J.; Jung, J.; Lee, N.; Kim, G.; Seo, Y. Study of Grains and Boundaries of Molybdenum Diselenide and Tungsten Diselenide Using Liquid Crystal. Nano Lett 2017, 17 (3), 1474-1481.

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

Table of Contents graphic

ACS Paragon Plus Environment

Page 20 of 25

Page 21 of 25 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

ACS Applied Materials & Interfaces

Figure 1 329x153mm (150 x 150 DPI)

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

Figure 2 256x191mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 22 of 25

Page 23 of 25 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

ACS Applied Materials & Interfaces

Figure 3 335x185mm (150 x 150 DPI)

ACS Paragon Plus Environment

ACS Applied Materials & Interfaces 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

Figure 4 346x280mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 24 of 25

Page 25 of 25 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

ACS Applied Materials & Interfaces

Figure 5 230x196mm (150 x 150 DPI)

ACS Paragon Plus Environment