Few-Layer MoS2 p-Type Devices Enabled by Selective Doping Using

Jan 20, 2016 - Technical Physics Division, BARC, Mumbai 400085, India. ACS Nano , 2016, 10 (2), .... Covalent Nitrogen Doping and Compressive Strain i...
0 downloads 7 Views 4MB Size
Subscriber access provided by ORTA DOGU TEKNIK UNIVERSITESI KUTUPHANESI

Article 2

Few Layer MoS p-Type Devices Enabled by Selective Doping Using Low Energy Phosphorus Implantation Ankur Nipane, Debjani Karmakar, Naveen Kaushik, Shruti Karande, and Saurabh Lodha ACS Nano, Just Accepted Manuscript • DOI: 10.1021/acsnano.5b06529 • Publication Date (Web): 20 Jan 2016 Downloaded from http://pubs.acs.org on January 24, 2016

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 free 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 accessible to all readers and 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.

ACS Nano 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 27

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 Nano

Few Layer MoS2

p-Type

Devices Enabled by

Selective Doping Using Low Energy Phosphorus Implantation Ankur Nipane

†,

Debjani Karmakar

‡,

Naveen Kaushik

Saurabh Lodha

†,

Shruti Karande

†,

†∗

† Department of Electrical Engineering, IIT Bombay, Mumbai 400076, India ‡Technical Physics Division, BARC, Mumbai 400085, India E-mail: [email protected]

Abstract P

-type doping of MoS2 has proved to be a signicant bottleneck in the realization

of fundamental devices such as

p-n

junction diodes and p -type transistors due to its

intrinsic n -type behavior. We report a CMOS compatible, controllable and area selective phosphorus plasma immersion ion implantation (PIII) process for p -type doping of MoS2 . Physical characterization using SIMS, AFM, XRD and Raman techniques was used to identify process conditions with reduced lattice defects as well as low surface damage and etching, 4X lower than previous plasma based doping reports for MoS2 . A wide range of non-degenerate to degenerate p -type doping is demonstrated in MoS2 eld eect transistors exhibiting dominant hole transport. Nearly ideal and air stable, lateral homogeneous p-n junction diodes with a gate-tunable rectication ratio as high as 2×104 are demonstrated using area selective doping. Comparison of XPS data from unimplanted and implanted MoS2 layers shows a shift of 0.67 eV towards lower 1

ACS Paragon Plus Environment

and

ACS Nano

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

binding energies for Mo and S peaks indicating p -type doping. First-principles calculations using density functional theory techniques conrm p -type doping due to charge transfer originating from substitutional as well as physisorbed phosphorus in top few layers of MoS2 . Pre-existing sulphur vacancies are shown to enhance the doping level signicantly.

Keywords MoS2 , CMOS, plasma doping, hole transport, p-n junction Increased short channel eects in conventional bulk semiconductors at scaled dimensions has led to extensive research on 2D semiconductors such as molybdenum disulphide (MoS 2 ). Their atomistically thin dimensions enable excellent gate control. 1 A large on-current/ocurrent ratio (ION /IOF F >108 ) along with a low sub-threshold slope have been demonstrated through experiments 2 and simulations 3 for MoS2 eld eect transistors (FETs). Moreover, direct bandgap in monolayer form also enables high performance optoelectronics, exible electronics and energy harvesting applications. However, these applications require complementary devices with p as well as n -type regions, which have been particularly challenging to realize using MoS 2 because of its natural propensity for unipolar ( n -type) transport. Strong Fermi level pinning near the conduction band, intrinsic defects and ambient interactions result in strong n -type conduction in MoS 2 . 4 The use of traditional doping techniques for controlled chemical doping of MoS 2 , especially p -type, is challenging due to its atomistically thin dimensions. This has led to the exploration of dierent doping techniques. A few reports have demonstrated p -type doping of MoS 2 using charge transfer, surface functionalization, plasma treatment and substitutional CVD doping. 58 Amongst these, plasma treatment emerges as the most eective doping technique for few layer MoS 2 due to a wide range of doping that can be achieved with good control and selectivity. However reported plasma studies suer from large surface damage, high rate of etching and lack of precise control over dopant concentration and location. This makes it challenging to dope monolayer 2

ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27

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 Nano

or few layer (5-10) MoS 2 akes. In addition, the dopants are not compatible with complementary metal oxide semiconductor (CMOS) processes. For example, uorine plasmas 9 can easily remove (ash) resists used for selective doping in CMOS manufacturing. In this work, we demonstrate p -type doping of few layer MoS 2 with phosphorus (P) using a CMOS compatible plasma immersion ion implantation (PIII) process. Raman spectroscopy, atomic force microscopy (AFM) and x-ray diraction (XRD) characterization were rst used to identify PIII process conditions with low surface damage and etching. Back-gated transistors fabricated using varying implant energy and dose (plasma exposure time) show clear evidence of p -type conduction with transfer characteristics varying from ambipolar to complete p -type indicating non-degenerate to degenerate p -type doping. Gate tunable, air stable, lateral p-n junction diodes formed through area selective p -type doping of MoS 2 with excellent rectifying characteristics have also been demonstrated. Finally, x-ray photoelectron spectroscopy (XPS) characterization of MoS 2 akes with and without the PIII process as well as ab-initio simulations based on density functional theory (DFT) were used to validate and delineate the doping mechanism.

Results and Discussion This work employs an indigenously built PIII system designed for achieving shallow doping proles. 10 A schematic of the PIII system is shown in Figure 1a. PIII is a CMOS compatible process which uses an inductively coupled plasma (ICP) source to generate low energy high density plasma. The amplitude (0-2 keV), on-time and time period (frequency) of the pulsed DC bias applied to the sample holder can be varied to reduce etching and surface damage. Pulse biasing was preferred over continuous biasing as it creates less damage during the pulse on-time and neutralizes unwanted charges accumulated on the sample. CMOS compatible phosphorus was selected as the doping species due to its high anity for molybdenum in sulphur decient MoS 2 . 11 We believe that a large number of naturally occurring sulphur

3

ACS Paragon Plus Environment

ACS Nano

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 1: (a) A schematic of the plasma immersion ion implantation (PIII) system used in this work, (b) SRIM simulation data showing implantation depth of P in MoS 2 for varying implant energies, (c) SIMS proles obtained from P-implanted MoS 2 showing shallow phosphorus doping with a sharp Gaussian nature compared to background sulphur, and, (d) cross-section schematics and energy band diagrams illustrating the impact of p -type doping in the top few layers of n -type MoS2 . vacancies in the exfoliated MoS 2 used in this study, further enhanced by plasma exposure, facilitate phosphorus doping. Stopping range of ions in matter (SRIM) simulations were used to estimate phosphorus implant depths in MoS 2 for varying implantation energies. 12 A projected range of 0-5 nm obtained from SRIM simulations for an implant energy range of 0-3 keV as shown in Figure 1b indicates shallow phosphorus doping contained within the top few MoS2 layers. Secondary ion mass spectroscopy (SIMS) data from implanted bulk MoS2 samples in Figure 1c shows a sharp Gaussian prole for phosphorus compared to the background sulphur. This indicates shallow surface doping and reinforces the SRIM

4

ACS Paragon Plus Environment

Page 4 of 27

Page 5 of 27

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 Nano

simulation results. The surface non-uniformity of the bulk crystal used for SIMS makes it dicult to measure the sputter depth accurately. However, a rst level approximation of sputter depth can be estimated using SRIM simulations. An extracted implantation depth of 3.5 nm for a 2 keV implant was used to determine the sputter depth for the SIMS data. P type doping of the top few MoS 2 layers due to either substitutional and/or physisorbed ++ species (P, PH3 , PH+ , etc.) is likely to deplete/invert the underlying n -type MoS2 2 , PH

layers through a charge transfer mechanism as shown through cross-section schematics and energy band diagrams in Figure 1d.

Figure 2: AFM micrographs and line scans of the same MoS 2 -SiO2 edge (a) before, and, (b) after a 2 keV PIII phosphorus implant indicating etching of the top 3-4 nm of the ake. AFM surface scans of the same MoS 2 ake (a) before, and, (b) after a 2 keV PIII phosphorus implant indicating a decrease of 0.52 nm in rms surface roughness. Implant process conditions resulting in minimal surface damage and etching were identied through AFM, Raman spectroscopy and XRD analysis of MoS 2 implanted with dierent energies. Figures 2a and 2b show etching of 3-4 nm from the top of the MoS 2 ake through comparative AFM line scans before and after implantation at 2 keV, the highest implant energy used in this work. This is nearly 4X lower when compared to previously employed plasma 5

ACS Paragon Plus Environment

ACS Nano

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 3: (a) Raman spectra, and, (b) FWHM values of Raman peaks ( E12g and A1g ) for unimplanted, 0 keV (unannealed), 1 keV (unannealed) and 2 keV (unannealed) MoS 2 samples. (c) Raman spectra, and, (d) FWHM values of Raman peaks ( E12g and A1g ) for unimplanted, 1 keV (unannealed) and 1 keV (annealed) MoS 2 samples. methods likely due to the pulsed DC bias (with variable on-time and frequency). 13 Figures 2c and 2d show rms surface roughness values measured using AFM on the same ake before and after a 2 keV implant. A value of 1.14 nm for the unimplanted sample is consistent with reported numbers. 14 2 keV implant decreases the roughness by 0.52 nm, again consistent with reports of decrease in Si surface roughness when etched in high density plasmas, 15 similar to the one used in our PIII chamber. This can be attributed to a higher ux of etch radicals seen by the crests as opposed to the valleys on the MoS 2 surface. Fluorine-based plasmas have been shown to increase the surface roughness of both Si 16 and MoS2 . 17 In this

6

ACS Paragon Plus Environment

Page 6 of 27

Page 7 of 27

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 Nano

context, the AFM results in this work are highly signicant as the rst report of reduced MoS2 surface roughness with plasma treatment. The reduced etching and improved surface roughness with PIII indicates that this method can also be extended to monolayer MoS 2 . The data for 2 keV implant suggests that 0 keV (no-bias) and 1 keV implants should result in lower MoS2 etch rates along with a smaller decrease in surface roughness. Nearly all the devices fabricated for electrical studies employed the 0 keV and 1 keV implant conditions for minimal surface etching and improved surface roughness. In addition, Raman spectra shown in Figure 3a indicate reduced implant-induced lattice distortions (defects) for lower energies. 18 Signicant broadening of the E12g peak indicates damage to the MoS 2 ake for the 2 keV implant process. However, a 7X reduction in the full width half maxima (FWHM) of the E12g peak indicates lower damage for 0 and 1 keV implants as shown in Figure 3b. Further, thermal annealing at 300 ◦ C after a 1 keV implant nearly recovers the full width half maxima (FWHM) for both E12g and A1g peaks to the values for a pristine, unimplanted ake as shown in Figure 3c and 3d. All p -type transistors and p-n junction diodes reported later in the manuscript were annealed at 300 ◦ C after the implant and before contact formation. Raman FWHM and AFM surface roughness data summarized in Table 1 shows the impact of varying energy (0, 1 and 2 keV) and annealing at 300 ◦ C on damage created in the MoS 2 samples by the PIII process. Low implant energy and thermal annealing result in a low damage process. Moreover, presence of the {002} crystalline peak for all three (0, 1 and 2 keV) implant conditions in the grazing and normal incidence angle XRD data of Figure S1 indicates that the plasma implant energy is too low to amorphize even the top few MoS 2 layers. Impact of the phosphorus implant on transport properties was evaluated by fabricating transistors using two dierent process sequences, (i) implant (doping) after, and, (ii) implant before source-drain contact patterning and metallization. Figures 4a and 4b show the impact of implantation on transfer characteristics of transistors where the implantation was carried out after source-drain contact formation. Thermal annealing for dopant activation

7

ACS Paragon Plus Environment

ACS Nano

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 27

Table 1: Comparison of Raman FWHM and AFM surface roughness values showing the impact of varying energy (0, 1 and 2 keV) and annealing at 300 ◦ C on damage created in MoS2 samples by the PIII process. Process Condition Unimplanted Implanted Implanted and annealed

Implant Energy (keV) 0 1 2 1

FWHM E12g 5.6 6.5 11 71.7 6.5

(cm −1 ) A11g 6.9 7.5 10.7 12.2 7.3

RMS surface roughness (nm) 1.14 0.62 -

was intentionally avoided in these devices as this could also modify the source-drain contact resistance. Unimplanted MoS 2 exhibits strong electron conduction ( n -type) with a high

ION /IOF F ratio (∼106 ) likely due to intrinsic sulphur vacancies, a low electron injection barrier due to Fermi level pinning near the conduction band and strong electrostatic coupling with gate voltage. However we observe signicant reduction in current along with a positive threshold voltage ( VT ) shift towards higher gate voltage after p -type doping of a few layers at the top as shown in Figure 4c. Figure 4d illustrates a cross-section schematic of the nal device. A decrease in current with the P implant can be attributed to reduction in the number of n -type layers due to p -type counter doping as well as etching of the top few layers. Also,

p -type layers at the top deplete the n -type MoS2 underneath of charge carriers and hence a larger gate voltage (positive VT shift) is now required to turn on (accumulate) the bottom

n -MoS2 transistor channel. These eects are also depicted through energy band diagrams plotted in Figure 4e for A-A' and B-B' cuts of Figure 4d. Although this data suggests p -type doping of the top few MoS 2 layers, signicant p -type transport was not observed since the contacts are still on top of n -MoS2 and strongly pinned near the conduction band. 19 Transistors where blanket implantation of the MoS 2 akes was done before sourcedrain contact formation were also fabricated and characterized. More ecient hole injection through tunneling from the contacts into p -MoS2 enables p -type transport in this case. The transistors were annealed at 300 ◦ C after the implant and before contact formation. As shown in Figure 5a, high implant dose (time=35 sec) at 1 keV results in weak gate control 8

ACS Paragon Plus Environment

Page 9 of 27

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 Nano

Figure 4: (a) Before, and, (b) after implant IDS vs VGS transfer characteristics of MoS 2 transistors with the implant carried out after source-drain metallization. (c) Comparative analysis of transfer characteristics before and after the implant for a xed VDS , (d) crosssection schematic of the nal device, and, (e) energy band diagrams along A-A' and B-B' cuts in the schematic of (d). at negative VGS and high current values indicating degenerate p -type doping. 20 Similar observation of weak gate modulation has been reported by Javey et al. 21 for degenerate n -type doping of MoS2 and by Zhao et al. 22 for degenerate p -type doping of WSe 2 . A quantitative estimate of the degenerate hole sheet doping concentration is given later in this manuscript. Non-degenerate doping obtained by lowering the dose (time=25 sec) at 0 keV energy results in ambipolar characteristics (Figure 5a) with hole conduction dominating for negative gate voltages and electron conduction for positive gate voltages. The change from electron to hole transport at negative VGS conrms p -type doping. Figure 5b shows the p -type transfer characteristics for the non-degenerately doped transistor with varying VDS and Figure 5c shows a cross-section schematic of the device with cuts along A-A' and B-B'. The ambipolar transfer characteristics can be explained through energy band diagrams for A-A' and B-B' (Figure

9

ACS Paragon Plus Environment

ACS Nano

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 5: IDS vs VGS transfer characteristics at a xed VDS of 3 V for (a) degenerately and non-degenerately p -doped MoS2 transistors with the implant and anneal carried out before source-drain metallization, and, (b) non-degenerately doped MoS 2 transistor with varying VDS . (c) Cross-section schematic of the nal device, and, (d) energy band diagrams along A-A' and B-B' cuts in the schematic of (c). 5d). For negative VGS , free hole concentration in the top p -MoS2 layers (A-A') exceeds free electron concentration in the bottom n -MoS2 layers resulting in dominant p -type transport. Similar arguments explain n -type transport for positive VGS . The overlap of electron and hole conduction along with electrostatic screening due to the intermediate n -MoS2 layers leads to a relatively low ION /IOF F ratio (∼102 ) for the hole transport regime at negative

VGS . Electrostatic screening can also be understood in terms of a dierence in capacitive coupling between the channel and gate for unimplanted and implanted transistors as explained in Figure S2 (supplementary information). Electron transport through the bottom

n -MoS2 layers for positive VGS is similar to the case of implanting after contact formation where reduction in number of n -channel layers, counter doping and charge depletion due to

10

ACS Paragon Plus Environment

Page 10 of 27

Page 11 of 27

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 Nano

p -type MoS2 layers on top leads to a reduction in current level with the implant. Quantitative analysis of transfer characteristics for the degenerately and non-degenerately doped p -type transistors in Figure 5a yields source-drain resistance ( RSD ) corrected 23 peak eld eect hole mobilities of 8.4 cm 2 /V-s and 137.7 cm2 /V-s respectively. These correspond to sheet hole concentrations of 2.4 ×1012 cm−2 and 1.3×1010 cm−2 respectively at VGS = 0 V. Similar combinations of mobility and sheet concentrations have been reported previously for non-degenerate 24 as well as degenerate 25 doping. The decrease in mobility with increasing doping is consistent with previous reports on Nb-doped p -MoS2 . 26 Assuming a range of 7 (entire ake)- 0.7 (monolayer) nm for thickness of the p -type doped MoS2 layers, we get 1.9×1016 - 1.9×1017 cm−3 and 3.4×1019 - 3.4×1020 cm−3 as the ranges for hole doping concentrations in the non-degenerately and degenerately doped transistors respectively. The classication of these sheet concentration values as non-degenerate and degenerate is consistent with the theoretical threshold of 7.7 ×1017 cm−3 for degenerate hole doping in MoS 2 . 27 Calculations for arriving at this theoretical value along with the detailed methodology used to extract mobilities and sheet concentrations are shown in the supplementary information. The demonstrated doping concentration range of 10 17 -1020 cm−3 (considering monolayer MoS 2 ) corresponding to hole sheet densities of 10 10 -1012 cm−2 is the widest reported till date using plasma treatment.

Figure 6: (a) 3D schematic of a lateral p-n junction diode at the time of the phosphorus implant, and, (b) SEM image of a fully processed p-n junction diode sample. Lateral p-n junction diodes were also fabricated through area selective PIII doping of 11

ACS Paragon Plus Environment

ACS Nano

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 7: (a) I-V characteristics of lateral p-n junction diode for varying bottom gate bias exhibit high rectication ratio (upto 2 ×104 ) and good ideality ( ∼1.2), (b) diode rectication ratio as a function of applied gate bias, (c) cross-section schematic of the nal p-n junction device, and, (d) energy band diagrams along A-A' and B-B' cuts in the schematic of (c). EBL patterned high quality multilayer akes ( ∼5-10 nm) exfoliated on Si-SiO 2 substrates as shown in the 3D schematic (Figure 6a), followed by resist stripping and annealing in nitrogen (N2 ) ambient at 300 ◦ C for 15 min. Finally, contacts were formed using EBL patterning, metallization and lifto. Figure 6b shows a scanning electron microscopy (SEM) image of a sample diode. Lateral homojunction p-n diodes can support higher levels of current density due to the lack of band-osets and higher in-plane mobilities when compared to vertical heterojunction diodes. 28 I -V characteristics in Figure 7a show a rectication ratio (RR) as high as 2×104 (Vbias = ±3 V) and an ideality factor of 1.2 at VBG = -80 V. The diode rectication ratio can be modulated using VBG as shown in Figure 7b. Figure 7c shows a cross-section schematic of the nal device while Figure 7d shows the energy band diagrams along the A-A' and B-B' cuts for the top p-n junction diode and the bottom n MoS2 transistor respectively. Variable diode rectication can be explained by comparing the 12

ACS Paragon Plus Environment

Page 12 of 27

Page 13 of 27

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 Nano

current components of the top p-n junction and the bottom n -type transistor at dierent gate and diode voltages. With increasing positive VBG the transistor channel turns on with high electron density and adds to the reverse bias leakage current of the diode. At the same time the high electron density reduces series resistance of the quasi-neutral n -MoS2 region resulting in higher forward bias on-current. As VBG is made zero and then negative, the bottom transistor turns o exponentially below VT resulting in an exponential decrease in diode reverse leakage that leads to an exponential increase in rectication ratio as seen in Figures 7a and 7b. This is also accompanied by a smaller decrease in on-current of the diode due to an increase in series resistance because of lower electron density. This suggests that true electrical performance of the p-n junction is observed for VBG = -80 V when the bottom transistor is completely switched o. Hence 2 ×104 refers to the absolute rectication ratio of the diode, which is also the highest achieved with plasma based doping of MoS 2 till date. To rule out the possibility of rectication due to dierent Schottky barrier heights at the contacts, output characteristics of unimplanted ( n -type) and doped (p -type, implantation before source-drain contact formation) transistors were measured as shown in Figure S3. Ti/Au contacts in both cases are ohmic in nature with symmetric I-V characteristics over ±1V of applied VDS . They are either p -type (hole transport modulated with negative

VBG ) or n -type (electron transport modulated with positive VBG ) depending on whether the channel is doped p -type or unimplanted n -type, respectively. The I-V characteristics indicate ecient electron (hole) transport/injection across the source-drain contact barriers for the unimplanted (implanted) sample for positive as well as negative VDS . Hence, rectication observed in the diode sample (Figure 7a) can only be attributed to the formation of a stable p-n junction. Impact of thermal annealing on diode and transistor characteristics can be seen in Figure S4 (a, b). A 9X increase in diode rectication ratio and decrease in ideality factor indicates signicant repair of implant-induced defects. For back-gated transistors, 13% increase in on-current and 3.5X increase in hole sheet concentration along with reduced gate control indicates increase in dopant activation and defect reduction due to thermal an-

13

ACS Paragon Plus Environment

ACS Nano

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

nealing. However, a 3X reduction in peak eld eect mobility was also observed which can be attributed to increased scattering due to higher sheet concentration. The implanted diodes show good air stability with minimal change in I-V characteristics over 7 days of ambient exposure as shown in Figure S4c. PIII process oers the ability to control dopant prole and concentration by varying energy and dose conditions. Transistors that received a xed implantation dose at dierent energies after source-drain contact formation were fabricated to demonstrate this capability of PIII. Figure 8a shows the IDS -VGS data indicating a monotonic decrease in transistor

ION with increasing implantation energy and a signicant increase in VT for the 1 and 2 keV higher energy implants. Figure 8b shows the average normalized change in ION and

VT obtained from measurements on multiple devices for each energy. Larger drop in current level and increase in VT with increasing implant energy can be explained by the reduction in number of n -type layers and their electron charge density due to the deeper p -type counter doping (see Figure 4 and accompanying text). Figure 8c shows the impact of varying implant dose (time) on electrical performance of p-n diode samples. Increasing the dose at constant implant energy of 0 keV increases (4X) the rectication ratio due to the higher doping concentration resulting in higher diusion on-current and reduction in series resistance. The high degree of controllability oered by the PIII process makes it useful for a large variety of device applications. X-ray photoelectron spectroscopy (XPS) characterization was done on MoS 2 samples before and after implantation to conrm and understand the doping mechanism. Figure 9a conrms the presence of phosphorus after implantation. Additionally, the P 2p 3/2 peak occurs at a binding energy (BE) of 133.4 eV which matches the peak position for P 2p 3/2 in molybdenum phosphide (MoP). 29 This indicates that the implanted phosphorus interacts chemically with MoS 2 and is present at a substitutional/interstitial site. BE shifts in XPS spectra for core level electrons have been used previously to identify the doping type in MoS2 . 30,31 Figures 9b and 9c compare the BEs for Mo 3d 5/2 , Mo 3d3/2 , S 2p3/2 and S 2p1/2

14

ACS Paragon Plus Environment

Page 14 of 27

Page 15 of 27

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 Nano

Figure 8: (a)IDS -VGS data of transistors implanted with varying energies but xed phosphorus dose, (b) normalized change in transistor ION and shift in VT versus implant energy, and, (c) I-V characteristics of p-n diodes fabricated using dierent implant doses at xed energy of 0 keV. Longer implant time (40 s) corresponds to higher dose.

Figure 9: (a) XPS spectrum for MoS 2 implanted with phosphorus showing presence of the P 2p3/2 peak. XPS spectra for (b) molybdenum (Mo 3d 5/2 and Mo 3d3/2 ), and (c) sulphur (S 2p3/2 and S 2p1/2 ) peaks before and after implantation showing a downward shift of 0.67 eV in binding energy for all four peaks. peaks before and after implantation on the same ake to avoid ambiguities arising from ake-to-ake variations. All four peaks shift by ∼0.67 eV to lower BEs with the implant indicating a shift in the Fermi level towards the valence band. This is equivalent to doping the MoS2 layers p -type. Hence, the XPS data conrms that phosphorus acts as a p -type dopant in MoS2 primarily occupying substitutional/interstitial sites. First principle calculations using DFT were also used to gain additional insights into the doping mechanism. Phosphorus has been theoretically predicted to result in p -type doping of MoS2 as a substitutional dopant. 32 However, the use of a PH 3 plasma in our experiments

15

ACS Paragon Plus Environment

ACS Nano

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 dope MoS2 p -type suggests three possibilities, viz. (1) P substitutional doping either by replacing S or occupying a S-vacancy (SV), (2) P-adatom adsorbed on the MoS 2 surface, and/or, (3) PH3 , PH2 + and PH++ ions adsorbed on the MoS 2 surface. Presence of S or Mo-vacancies (MoV) will also inuence the doping nature of all these three cases. Since MoS2 has a higher propensity for SV, the simulation study focuses more on analyzing these doping cases in the presence of SVs. All these three cases have been constructed for a freestanding MoS2 monolayer (1H-MoS 2 ) and the ionic positions and lattice parameters of the full assembly have been geometrically optimized following the conjugate gradient algorithm until the Hellmann-Feynman forces acting on each atom are less than 0.01 eV/Å. 33 Although the P-substituted, P-adatom and PH 3 -adsorbed cases result in a lowering of the total energy compared to 1H-MoS 2 , adsorption of PH2 + and PH++ ions results in an unstable system with much higher total energy. In addition, the relaxed system with these ions, without any selective dynamics, results in only the P-adatom adsorbed on MoS 2 surface with the H-ions getting stripped o from P and placed closer to S. Therefore, for the third case, we have only dealt with PH3 -adsorption. A comparison of the total energy for all these cases indicates that P-substitution and P-adatom are energetically more favourable than PH 3 -adsorption. Presence of a nearest neighbour (NN) SV favours the P-substituted, P-adatom and PH 3 adsorbed cases by 981, 937 and 35 meV respectively than the corresponding without-vacancy systems. Presence of MoV favors formation of only the rst two cases by 655 and 562 meV respectively. Comparison of ionic binding energies (IBE) for these three cases with pristine 1H-MoS2 , as presented in Table S1, indicates a considerable reduction in the IBE of Mo and S for all three cases, which gets enhanced further in the presence of vacancies. This reduction in IBE is consistent with the reduction in BEs observed in XPS results. A detailed discussion on IBE is presented in the supplementary information. In addition to this, impact of interlayer separation on p -type doping of MoS 2 using phosphorus is also discussed in supplementary information. Here, P atom and PH 3 molecule placed in the interlayer space of bilayer MoS2 show signicant doping in presence of SV as shown in Figure S5.

16

ACS Paragon Plus Environment

Page 16 of 27

Page 17 of 27

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 Nano

Figure 10: OPDOS and APDOS plots for the three cases: (a) P-substituted, (b) P-adatom, and, (c) PH3 -adsorbed. For each panel Figure (i) represents the unimplanted 1H-MoS 2 DOS for comparison. Figures (ii) and (iii) depict the DOS without and with SV respectively. For all these three cases P partial DOS (green) has been multiplied by ten times to improve its visibility over the same scale. States close to the Fermi level have been highlighted to indicate possible changes resulting from P incorporation in the system. We have investigated the orbital (OPDOS) and atom-projected (APDOS) density of states for all three cases mentioned above. Figure 10a depicts OP and APDOS plots for the P-substituted case. We have substituted one single S in 1H-MoS 2 by a P atom corresponding to a doping percentage of 3.125 %. In Figure 10a-(i), total DOS and OPDOS for Mo- 4d and S-3p are plotted for 1H-MoS 2 , showing a bandgap of ∼1.8 eV with strongly hybridized Mo-

4d and S-3p orbital levels constituting the top of the valence band (VBT) and bottom of the conduction band (CBB). For P-substituted case, Figure 10a-(ii) in the same panel depicts a shift of the Fermi level ( EF ) towards the valence band indicating p -type doping. This system shows a dierence in the Mo-S (2.42 Å) and Mo-P (2.39 Å) bond-lengths after geometrical optimization. The impact of shorter bond-length can be clearly observed from the overlapping charge spheres around Mo and P in the charge-density plot of Figure 11a. There is no such overlap for Mo-S bonds in the same isosurface scale. This dierence in

17

ACS Paragon Plus Environment

ACS Nano

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 11: Charge-density plots for (a) P-substituted without SV, (b) P-adatom without SV, and, (3) PH3 -adsorbed without SV. bond-length increases in presence of a NN SV. The Mo-P bond length reduces to 2.37 Å which leads to more charge sphere overlap. Thus, although the presence of SV initially leads to n -type doping within the system, vacancy-induced structural distortion leads to signicant charge transfer from Mo to P. To obtain a visible eect on OPDOS of Figure 10a(iii), we have kept the concentration of P-dopant (6.25 %) as double the concentration of SV (3.125 %). Due to charge transfer from Mo- 4d to P-3p, Mo-4d levels (marked red in gure) become partially unoccupied compared to Figure 10a-(i). These unoccupied levels are highly hybridized with P- 3p levels (green) indicating the charge transfer direction. Presence of SV results in an additional shift in EF towards the valence band indicating increased p -type doping, as can be seen from Figure 10a-(iii). The levels are highlighted near EF in all the gures to indicate the extent of p -type doping. DOS plots for adsorption cases (2) and (3) are also presented in Figures 10b and 10c respectively. For P-adatom, the charge-density plot (Figure 11b) implies that a signicant amount of the charge transfer happens from Mo to P via NN S atoms due to van der Waals-type bonding of P over S-layer of MoS 2 . Figure 10b-(ii) indicates that Mo and S occupied levels shift to a more bound state because of this charge transfer. This particular case shows no p -type doping without the presence of a vacancy. In presence of SV, after relaxation, the P-adatom is seen to settle at the SV position but with a shorter bond-length ( ∼2.37 Å) with the Mo-ion and thereby charge transfer from Mo to P 18

ACS Paragon Plus Environment

Page 18 of 27

Page 19 of 27

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 Nano

leads to partially unoccupied Mo- 4d levels hybridized with P- 3p near EF . Qualitatively, this case becomes exactly similar to a P-substituted system, as can be inferred from the close resemblance between the DOS gures presented in Figures 10a-(ii) and 10b-(iii). The extent of p -type doping appears to be slightly larger for the latter case because of shorter Mo-P bond length compared to the P-substituted case. From an experimental perspective this also implies that presence of SV promotes conversion of any adsorbed P into a substituted one. Adsorption of PH3 , however, is not energetically favourable for p -type doping. PH3 does not lead to signicant charge transfer from 1H-MoS 2 because of its attened nature of self-charge distribution. As can be seen from Figure 10c-(ii), adsorbed PH 3 only leads to some additional gap-states mostly resulting from Mo- 4d and P-3p states. In presence of SV, similar to the P-adatom case, PH 3 molecule also tries to occupy the SV, albeit remaining at a slightly higher level than the S-plane, leading to little charge transfer from the nearest neighbour Mo resulting in only a small shift of EF towards the valence band, as seen in Figure 10c-(iii). Comparison of Figures 11b and 11c shows that the amount of charge transfer to P-adatom is far larger than for the PH 3 -adsorbed case. This is consistent with the lesser IBE of P for case (3) (Table S1). Thus, phosphorus implantation using a low energy PIII process leads to a mixed state of all the above three possible congurations with more favourable chances for P-substitution, resulting in p -type doping of 1H-MoS 2 .

Conclusion In conclusion, we have demonstrated p -type doping of MoS 2 using a low damage phosphorus PIII process. High level of control and area selectivity enabled by PIII is shown to result in a wide range of doping concentration and potential use in various device applications. Dominant p -type transport in eld eect transistors as well as air stable, lateral p-n junction diodes with a high rectication ratio (2 ×104 ) demonstrate electrical proof of selective p type doping of MoS 2 . XPS analysis conrms the p -nature of phosphorus doping which is

19

ACS Paragon Plus Environment

ACS Nano

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

further reinforced by DFT simulation data for substitutional and adsorbed P in MoS 2 layers, especially in the case of pre-existing S vacancies. Although data for monolayer MoS 2 has not been shown, the low damage PIII process can be easily adapted to fabricate monolayer devices. These ndings open up the roadmap for development of complementary ( p and n ) and p-n junction based MoS 2 devices as well as the use of PIII doping technique for other emerging 2D materials.

Methods Transistor fabrication started with reactive ion etching of thermally grown SiO 2 on a degenerately doped Si substrate. MoS 2 was mechanically exfoliated from molybdenite crsytals (SPI supplies) using the scotch tape method. Large area, multilayer akes were identied using optical, Raman and atomic force microscopy. Transistors were fabricated on the selected akes using bilayer resist (PMMA and EL-9) with electron beam lithography (EBL, Raith 150-Two) for source-drain contact patterning. In this work, gold (Au) has been used as the primary contact metal with a thin Ti adhesion layer, since it has been shown to form ohmic contacts on MoS 2 . Samples were then transferred to the PIII chamber for the phosphorus implant through PH 3 /He plasma exposure either before or after source-drain contact patterning and metallization according to the sample type. The samples were kept at an appropriate distance from the plasma source to minimize surface damage. The indigenously built PIII system uses a pulsed DC bias source with a bias range of 0 to 5 kV. Chamber pressure was kept constant at 1.1E-1 mbar with RF power set at 1000 W for all the implant runs. A process ow similar to the one used for transistor fabrication was followed for the diode samples. Multilayer akes exfoliated on SiO 2 /Si substrates were rst patterned for area selective plasma exposure using EBL for forming the p-n junction. PIII chamber was used to dope MoS2 p -type using phosphorus through a PH 3 /He plasma implant. The resist

20

ACS Paragon Plus Environment

Page 20 of 27

Page 21 of 27

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 Nano

was then stripped by soaking the samples in acetone for 2 hours. After resist stripping, the samples were annealed in N 2 ambient at 300 ◦ C for 15 min using a 4 inch quartz furnace. Finally, contacts to the p and n regions of the MoS2 ake were fabricated using EBL, metallization and lifto. Room temperature electrical transport studies were conducted using a probe station connected to an Agilent BX-51 parameter analyzer. Photoemission XPS studies were done using a PHI Versaprobe system equipped with an Al k α source. SIMS doping proles were acquired using TOF-SIMS (PHI Trift V Nano) equipped with LMIG and Cesium ion guns. XRD analysis was performed using RIGAKU Smartlab where both normal and GIA modes were used. Experimental observation of p -type doping of MoS 2 through a phosphorus implant was veried and analyzed with the help of spin-polarized density functional calculations. We have used projector augmented wave (PAW) potentials with Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional within generalized gradient approximation (GGA), as implemented in Vienna Ab-initio Simulation Package (VASP). 34 The plane wave cut-o energy was 500 eV and Monkhorst-Pack (MP) meshes of 17 x 17 x 1 and 5 x 5 x 3 were used for complete geometry optimization and electronic structure calculations of substitutional and adsorbed cases respectively. Van der Waals interactions were incorporated using DFT-D2 method of Grimme as in VASP. We have used a 4 x 4 x 1 MoS 2 monolayer to simulate dierent cases resulting from the various possible experimental scenarios.

Acknowledgement The authors acknowledge the Confederation of Indian Industry (CII) and the Department of Science and Technology (DST), Government of India for funding this work. Debjani Karmakar acknowledges the support of BARC-ANUPAM supercomputing facility.

21

ACS Paragon Plus Environment

ACS Nano

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

Supporting Information Available Supporting information gives details regarding (i) impact of phosphorus implant on the crystallinity of MoS2 using XRD characterization, (ii) capacitor models to explain gate voltage dependence of the transfer characteristics of unimplanted ( n -type) and implanted ( p -type) MoS2 transistors, (iii) extraction methodology for hole eld eect mobilities and sheet concentrations from IDS -VGS data and calculation of theoretical limit for degenerate hole doping in MoS2 , (iv) the output characteristics of unimplanted and implanted MoS 2 transistors showing ohmic nature of contacts to n and p -type MoS2 respectively, (iv) the impact of thermal annealing and air exposure on electrical characteristics of diodes and transistors, (v) DFT analysis of ionic binding energies for various experimental scenarios of P incorporation in the MoS2 system and their correlation with observed XPS binding energy shifts, and, (vi) DFT analysis of the impact of interlayer separation on p -type doping of MoS 2 . This material is available free of charge via the Internet at http://pubs.acs.org/ .

References 1. Mak, K. F.; Lee, C.; Hone, J.; Shan, J.; Heinz, T. F. Atomically Thin MoS2 : A New Direct-Gap Semiconductor. Phys. Rev. Lett. 2010, 105, 136805136808. 2. Radisavljevic, B.; Radenovic, A.; Brivio, J.; Giacometti, V.; A.Kis, Single-Layer MoS 2 Transistors. Nat. Nanotechnol. 2011, 6, 147150. 3. Yoon, Y.; Ganapathi, K.; Salahuddin, S. How Good Can Monolayer MoS 2 Transistors Be? Nano Lett. 2011, 11, 37683773, PMID: 21790188. 4. Kaushik, N.; Nipane, A.; Basheer, F.; Dubey, S.; Grover, S.; Deshmukh, M. M.; Lodha, S. Schottky Barrier Heights for Au and Pd Contacts to MoS 2 . Appl. Phys. Lett. 2014, 105, 113505113508.

22

ACS Paragon Plus Environment

Page 22 of 27

Page 23 of 27

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 Nano

5. Chen, M.; Nam, H.; Wi, S.; Priessnitz, G.; Gunawan, I. M.; Liang, X. Multibit Data Storage States Formed in Plasma-Treated MoS 2 Transistors. ACS Nano 2014, 8, 4023 4032, PMID: 24680193. 6. Mouri, S.; Miyauchi, Y.; Matsuda, K. Tunable Photoluminescence of Monolayer MoS 2 via Chemical Doping. Nano Lett. 2013, 13, 59445948, PMID: 24215567. 7. Najmaei, S.; Zou, X.; Er, D.; Li, J.; Jin, Z.; Gao, W.; Zhang, Q.; Park, S.; Ge, L.; Lei, S.; Kono, J.; Shenoy, V. B.; Yakobson, B. I.; George, A.; Ajayan, P. M.; Lou, J. Tailoring the Physical Properties of Molybdenum Disulde Monolayers by Control of Interfacial Chemistry. Nano Lett. 2014, 14, 13541361, PMID: 24517325. 8. Zhang, K.; Feng, S.; Wang, J.; Azcatl, A.; Lu, N.; Addou, R.; Wang, N.; Zhou, C.; Lerach, J.; Bojan, V.; Kim, M. J.; Chen, L.-Q.; Wallace, R. M.; Terrones, M.; Zhu, J.; Robinson, J. A. Manganese Doping of Monolayer MoS 2 : The Substrate is Critical. Nano

Lett. 2015, 15, 65866591, PMID: 26349430. 9. Nam, H.; Wi, S.; Rokni, H.; Chen, M.; Priessnitz, G.; Lu, W.; Liang, X. MoS 2 Transistors Fabricated via Plasma-Assisted Nanoprinting of Few-Layer MoS 2 Flakes into Large-Area Arrays. ACS Nano 2013, 7, 58705881, PMID: 23790007. 10. Ray, A.; Nori, R.; Bhatt, P.; Lodha, S.; Pinto, R.; Rao, V. R.; Jomard, F.; NeumannSpallart, M. Optimization of a Plasma Immersion Ion Implantation Process for Shallow Junctions in Silicon. J. Vac. Sci. Technol. A 2014, 32, 061302061307. 11. Carvalho, A.; Neto, A. H. C. Donor and Acceptor Levels in Semiconducting TransitionMetal Dichalcogenides. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 89, 081406 081410. 12. Ziegler, J. F. SRIM-2003. Nucl. Instrum. Methods Phys. Res., Sect. B 2004, 219 - 220, 1027  1036.

23

ACS Paragon Plus Environment

ACS Nano

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

13. Chen, M.; Nam, H.; Wi, S.; Ji, L.; Ren, X.; Bian, L.; Lu, S.; Liang, X. Stable Few-Layer MoS2 Rectifying Diodes Formed by Plasma-Assisted Doping. Appl. Phys. Lett. 2013,

103, 142110142113. 14. Du, Y.; Liu, H.; Neal, A.; Si, M.; Ye, P. Molecular Doping of Multilayer MoS 2 FieldEect Transistors: Reduction in Sheet and Contact Resistances. IEEE Electron Device

Lett. 2013, 34, 13281330. 15. Martin, M.; Cunge, G. Surface Roughness Generated by Plasma Etching Processes on Silicon. J. Vac. Sci. Technol., B: Microelectron. Nanometer Struct.Process., Meas.,

Phenom. 2008, 26, 12811288. 16. Kokkoris, G.; Constantoudis, V.; Angelikopoulos, P.; Boulousis, G.; Gogolides, E. Dual Nanoscale Roughness on Plasma-Etched Si Surfaces: Role of Etch Inhibitors. Phys. Rev.

B: Condens. Matter Mater. Phys. 2007, 76, 193405193408. 17. Chen, M.; Wi, S.; Nam, H.; Priessnitz, G.; Liang, X. Eects of MoS 2 Thickness and Air Humidity on Transport Characteristics of Plasma-Doped MoS 2 Field-Eect Transistors.

J. Vac. Sci. Technol., B: Microelectron. Nanometer Struct.Process., Meas., Phenom.

2014, 32, 152105152110. 18. Kang, N.; Paudel, H. P.; Leuenberger, M. N.; Tetard, L.; Khondaker, S. I. Photoluminescence Quenching in Single-Layer MoS 2 via Oxygen Plasma Treatment. J. Phys. Chem.

C 2014, 118, 2125821263. 19. Gong, C.; Colombo, L.; Wallace, R. M.; Cho, K. The Unusual Mechanism of Partial Fermi Level Pinning at Metal-MoS 2 Interfaces. Nano Lett. 2014, 14, 17141720, PMID: 24660782. 20. Kiriya, D.; Tosun, M.; Zhao, P.; Kang, J. S.; Javey, A. Air-Stable Surface Charge Transfer Doping of MoS2 by Benzyl Viologen. J. Am. Chem. Soc. 2014, 136, 78537856, PMID: 24836497. 24

ACS Paragon Plus Environment

Page 24 of 27

Page 25 of 27

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 Nano

21. Fang, H.; Tosun, M.; Seol, G.; Chang, T. C.; Takei, K.; Guo, J.; Javey, A. Degenerate nDoping of Few-Layer Transition Metal Dichalcogenides by Potassium. Nano Lett. 2013,

13, 19911995, PMID: 23570647. 22. Zhao, P.; Kiriya, D.; Azcatl, A.; Zhang, C.; Tosun, M.; Liu, Y.-S.; Hettick, M.; Kang, J. S.; McDonnell, S.; KC, S.; Guo, J.; Cho, K.; Wallace, R. M.; Javey, A. Air Stable p-Doping of WSe 2 by Covalent Functionalization. ACS Nano 2014, 8, 1080810814, PMID: 25229426. 23. Kaushik, N.; Karmakar, D.; Nipane, A.; Karande, S.; Lodha, S. Interfacial n-Doping Using an Ultrathin TiO 2 Layer for Contact Resistance Reduction in MoS 2 . ACS Appl.

Mater. Interfaces 2016, 8, 256263, PMID: 26649572. 24. Park, H.-Y.; Lim, M.-H.; Jeon, J.; Yoo, G.; Kang, D.-H.; Jang, S. K.; Jeon, M. H.; Lee, Y.; Cho, J. H.; Yeom, G. Y.; Jung, W.-S.; Lee, J.; Park, S.; Lee, S.; Park, J.-H. Wide-Range Controllable n-Doping of Molybdenum Disulde (MoS 2 ) through Thermal and Optical Activation. ACS Nano 2015, 9, 23682376, PMID: 25692499. 25. Fang, H.; Chuang, S.; Chang, T. C.; Takei, K.; Takahashi, T.; Javey, A. HighPerformance Single Layered WSe 2 p-FETs with Chemically Doped Contacts. Nano Lett.

2012, 12, 37883792, PMID: 22697053. 26. Laskar, M. R.; Nath, D. N.; Ma, L.; Lee, E. W.; Lee, C. H.; Kent, T.; Yang, Z.; Mishra, R.; Roldan, M. A.; Idrobo, J.-C.; Pantelides, S. T.; Pennycook, S. J.; Myers, R. C.; Wu, Y.; Rajan, S. P-Type Doping of MoS 2 Thin Films Using Nb. Appl. Phys. Lett. 2014, 104, 092104092107. 27. Suh, J.; Park, T.-E.; Lin, D.-Y.; Fu, D.; Park, J.; Jung, H. J.; Chen, Y.; Ko, C.; Jang, C.; Sun, Y.; Sinclair, R.; Chang, J.; Tongay, S.; Wu, J. Doping against the Native Propensity of MoS2 : Degenerate Hole Doping by Cation Substitution. Nano Lett. 2014, 14, 6976 6982, PMID: 25420217. 25

ACS Paragon Plus Environment

ACS Nano

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

28. Lee, C.-H.; Lee, G.-H.; der Zande, A. M. V.; 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, 676681. 29. Xiao, P.; Sk, M. A.; Thia, L.; Ge, X.; Lim, R. J.; Wang, J.-Y.; Lim, K. H.; Wang, X. Molybdenum Phosphide as an Ecient Electrocatalyst for the Hydrogen Evolution Reaction. Energy Environ. Sci. 2014, 7, 26242629. 30. Yang, L.; Majumdar, K.; Liu, H.; Du, Y.; Wu, H.; Hatzistergos, M.; Hung, P. Y.; Tieckelmann, R.; Tsai, W.; Hobbs, C.; Ye, P. D. Chloride Molecular Doping Technique on 2D Materials : WS 2 and MoS2 . Nano Lett. 2014, 14, 62756280, PMID: 25310177. 31. Lin, J. D.; Han, C.; Wang, F.; Wang, R.; Xiang, D.; Qin, S.; Zhang, X.-A.; Wang, L.; Zhang, H.; Wee, A. T. S.; Chen, W. Electron-Doping-Enhanced Trion Formation in Monolayer Molybdenum Disulde Functionalized with Cesium Carbonate. ACS Nano

2014, 8, 53235329, PMID: 24785254. 32. Dolui, K.; Rungger, I.; Das, P. C.; Sanvito, S. Possible Doping Strategies for MoS 2 Monolayers : An Ab Initio Study. Phys. Rev. B: Condens. Matter Mater. Phys. 2013,

88, 075420075428. 33. Kresse, G.; Furthmüller, J. Ecient Iterative Schemes for Ab Initio Total-Energy Calculations using a Plane-Wave Basis Set. Phys. Rev. B: Condens. Matter Mater. Phys.

1996, 54, 1116911186. 34. Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B:

Condens. Matter Mater. Phys. 1993, 47, 558561.

26

ACS Paragon Plus Environment

Page 26 of 27

Page 27 of 27

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 Nano

Graphical TOC Entry

Implant = 1 keV Energy

27

ACS Paragon Plus Environment