Activating MoS2 for pH-Universal Hydrogen Evolution Catalysis

Oct 25, 2017 - Department of Electrical and Computer Engineering, North Carolina State University, Raleigh, North Carolina 27695, United States. J. Am...
1 downloads 11 Views 1MB Size
Subscriber access provided by University of Florida | Smathers Libraries

Article

Activating MoS2 for pH-Universal Hydrogen Evolution Catalysis Guoqing Li, Du Zhang, Yifei Yu, Shengyang Huang, Weitao Yang, and Linyou Cao J. Am. Chem. Soc., Just Accepted Manuscript • DOI: 10.1021/jacs.7b07450 • Publication Date (Web): 25 Oct 2017 Downloaded from http://pubs.acs.org on October 25, 2017

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.

Journal of the American Chemical Society 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 22

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

Journal of the American Chemical Society

Activating MoS2 for pH-Universal Hydrogen Evolution Catalysis Guoqing Li1, Du Zhang2, Yifei Yu1, Shengyang Huang3, Weitao Yang2*, and Linyou Cao1,4,5* 1

Department of Materials Science and Engineering, North Carolina State University, Raleigh, NC 27695; 2Department of Chemistry, Duke University, Durham, NC 27708; 3 Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, 27695; 4 Department of Physics, North Carolina State University, Raleigh, NC 27695; 5Department of Electrical and Computer Engineering, North Carolina State University, Raleigh, NC 27695 ABSTRACT: MoS2 presents a promising catalyst for the hydrogen evolution reaction (HER) in water splitting, but its worse catalytic performance in neutral and alkaline media than in acidic environment may be problematic for practical application. This is because the other half reaction of water splitting, i.e., oxygen evolution reaction, often needs to be implemented in alkaline environment. Here we demonstrate a universal strategy that may be used to significantly improve the HER catalysis of MoS2 in all kinds of environment from acidic to alkaline, proton intercalation. Protons may be enabled to intercalate between monolayer MoS2 and underlying substrates or in the interlayer space of thicker MoS2 by two processes: electrochemically polarizing MoS2 at negative potentials (vs. RHE) in acidic media or immersing MoS2 into certain acid solutions like TFSI. The improvement in catalytic performance is due to the activity enhancement of the active sites in MoS2 by the intercalated protons, which might be related with the effect of the intercalated protons on electrical conductance and the adsorption energy of hydrogen atoms.

The

enhancement in catalytic activity by the intercalated proton is very stable even in neutral and alkaline electrolytes.

* To whom correspondence should be addressed.

ACS Paragon Plus Environment

1

Journal of the American Chemical Society

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 2 of 22

Email: [email protected], [email protected] INTRODUCTION The hydrogen evolution from water splitting is critical for the utilization of hydrogen energy, but its implementation has been delayed by the lack of low-cost high-performance catalysts that may be functional in all kinds of pH environment from acidic to alkaline.1-2 Earth abundant transition metal chalcogenide materials like molybdenum disulfide (MoS2) may provide a promising catalyst for the hydrogen evolution reaction (HER).1-2 However, despite a considerable amount of research effort,2-15 the catalytic performance of MoS2 remains unsatisfactory and is even worse in neutral and alkaline environment than in acidic media.16-17 Therefore, it is necessary to improve the catalytic performance of MoS2, in particular, in neutral or alkaline environment. This is important for practical application because the other half reaction of water splitting, i.e. oxygen evolution reaction, often requires to be implemented in neutral or alkaline electrolytes. Here we demonstrate a universal strategy to substantially improve the HER catalytic performance of MoS2 in all kinds of environment from acidic to alkaline, proton intercalation. The proton intercalation may be enabled by electrochemically polarizing MoS2 at negative potential (vs. RHE) in acidic electrolytes or immersing MoS2 into certain acid solution like TFSI. It may happen between monolayer MoS2 and underlying substrates or in the interlayer space of thicker MoS2. The improvement in catalytic performance results from the activity enhancement of the active sites in MoS2, including sulfur vacancies or edge sties18, by the intercalated protons. And this might be related with the effect of the intercalated protons on electrical conductance and the adsorption energy of hydrogen atoms. The enhancement effect of the intercalated protons for the catalytic activity is stable even in neutral and alkaline environment.

ACS Paragon Plus Environment

2

Page 3 of 22

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

Journal of the American Chemical Society

EXPERIMENTAL SECTION Synthesis and transfer of monolayer MoS2 films and flakes: MoS2 films were synthesized using a self-limiting chemical vapor deposition (CVD) process that we have recently developed19. Briefly molybdenum chloride (MoCl5) powder (99.99%, Sigma-Aldrich) was placed at the center of the furnace and sulfur powder (Sigma-Aldrich) were placed at the upstream entry of the furnace. Receiving substrates (sapphire) were placed in the downstream of the tube. Typical conditions include a temperature of 850 °C, a flow rate of 50 sccm, and a pressure around 2 Torr. Monolayer MoS2 flakes were grown using a different CVD20, in which MoO3 (99.99%, SigmaAldrich) instead of MoCl5 was used as the precursor. Typical growth was performed at 750 °C under a flow of Ar gas in rate of 100 sccm and ambient pressure. The transfer of the monolayers followed a surface-energy-assisted transfer approach that we have developed previously.21 Briefly, a layer of polystyrene (PS) was spin-coated on the as-grown monolayers. A water droplet was then dropped on the top of the polymer. Due to the different surface energies of the monolayer and the substrate, water molecules could penetrate under the monolayer, resulting the delamination of the PS-monolayer assembly. We could pick up the polymer/monolayer assembly with a tweezers and transferred it to glassy carbon substrates. Finally, PS was removed by rinsing with toluene several times. Characterizations: Raman and PL measurements were carried out by a Horiba xPlora system equipped with an excitation wavelength at 532 nm. XPS measurements were carried out at X-ray photoelectron spectroscope (SPECS System with PHOIBOS 150 analyzer using an Mg Kα X-ray source). SEM and XRD measurements were performed at scanning electron microscope (SEM, JEOL JSM-6400F) and Rigaku SmartLab X-ray diffractometer using Cu-Kα radiation, respectively. The electrochemical characterization was performed in 0.5 M H2SO4 using a CH

ACS Paragon Plus Environment

3

Journal of the American Chemical Society

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 4 of 22

Instrument electrochemical analyzer (Model CHI604D) with a saturated calomel reference electrode (SCE) using a one-compartment cell system. A Pt mesh (2.54cm x 2.54cm) or a graphite rod were used as counter electrode, and the results did not show any obvious difference (Fig. S1). Nitrogen gas was bubbled into the electrolyte throughout the experiment. The potential shift of the SCE is calibrated to be -0.262 V vs. RHE. Typical electrochemical characterizations of the monolayers were performed using linear sweeping from +0 V to -0.5V (vs. RHE) with a scan rate of 5 mV/s. The electrochemical characterization in neutral solution was performance in 1M PBS solution. The mild base solution (pH=12) was prepared by mixing 327.0 mL 0.05M Na2HPO4 (99%, Sigma-Aldrich) and 173.0 mL 0.1M sodium hydroxide (NaOH) (98.5%, SigmaAldrich). DFT Computation: All computations in this work were performed with the plane-wave software code VASP (Vienne ab initio simulation package)22using the Perdew-Burke-Ernzerhof (PBE) exchange-correlation energy functional23 with its corresponding pseudopotentials. The semiempirical D2 model 24 was applied to describe the Van der Waals interactions. Periodic boundary conditions were used in all three directions for our supercell simulation models. A vacuum of >10 Å was employed in the non-periodic directions. The H+-intercalated models were calculated with one electron removed from the corresponding neutral systems. An opposite background charge is always adopted which is uniformly distributed throughout the simulation box, which may make the absolute value of the total energy behave well under periodic boundary conditions. This background charge only shifts the Hamiltonian by a global constant, and whenever relative energies are calculated, this same constant always cancels out upon subtraction. A plane-wave cutoff of 500 eV was chosen. For the Mo edge model, the Gamma-centered 2×1×1 k-point grid was applied for the geometry optimizations, and the 6×1×1 k-point grid for the single point

ACS Paragon Plus Environment

4

Page 5 of 22

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

Journal of the American Chemical Society

energy calculations. For the S defect models, the Gamma-centered 2×2×1 k-point grid was applied for the geometry optimizations, and the 4×4×1 k-point grid for the single point energy calculations. The geometry optimizations were converged to 0.01 eV Å-1. The zero-point-energy (ZPE) correction was performed with the finite difference method. RESULTS AND DISCUSSION

Fig.1. Improved catalytic activity of monolayer MoS2 with electrochemical cycling in acidic electrolyte. Polarization curves of (a) typical monolayer MoS2 films, (b) monolayer MoS2 flakes in lateral size of around 2 µm, and (c) monolayer MoS2 film with < 2% sulfur vacancies with different cycles. The electrolyte is 0.5 M H2SO4. The number of cycles is given as shown. Insets of (a) and (b) are optical images of typical monolayer MoS2 films (in size of around 1 cm) and flakes (in size of around 2 µm). Corresponding Tafel slopes of these polarization curves are given in Fig. S2. The density of sulfur vacancies may be derived from XPS measurement as illustrated in Fig. S3. Fig. 1 shows the polarization curves collected from typical monolayer MoS2 films (Fig. 1a) and monolayer MoS2 flakes (Fig. 1b). The monolayer films and flakes were first grown on sapphire substrates using well-established chemical vapor deposition (CVD) processes19-20 and then transferred onto glassy carbon substrates for electrochemical characterizations.18, 21 The polarization curves were obtained using linear sweep voltammetry (LSV) from 0 V (vs RHE) to a voltage that can enable a current density of around 20 mA/cm2 as shown in Fig.1. Cyclic voltammetry (CV) was performed over similar voltages between each of the linear sweeps. Very interestingly, the catalytic performance of the film and the flakes both substantially increases with the cycling and eventually approaches to be stable. Usually the films request more cycles to

ACS Paragon Plus Environment

5

Journal of the American Chemical Society

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

Page 6 of 22

reach the stable performance than the flakes. We also find that more cycles are needed for the flakes with larger size than smaller flakes. Additionally, extending the cycling range down to more negative voltages may lower the number of the cycles needed to reach stable activity, but may not change the level of the stable activity (Fig. S4). Similar increase in the catalytic performance of MoS2 with electrochemical cycling has actually been reported previously,8, 15, 25 but no understanding for the underlying mechanism is available yet except vaguely calling it as an activation process.

Fig.2. p-doping to monolayer MoS2 by electrochemical cycling. (a) XPS, (b) Raman, and (c) photoluminescence (PL) spectra of monolayer MoS2 films before (blue curves) and after (red curves) the electrochemical cycling that may enable stable catalytic activities. (d) Schematic illustration for the major role of proton intercalation between the monolayer and underlying structures. The protons adsorbed on top of the monolayer may only have a minor role. The characterization results for monolayer MoS2 flakes before and after electrochemical cycling are given in Fig. S8.

We may conclude that no substantial amount of new active site is generated during the cycling process. Our previous studies have demonstrated that MoS2 have two kinds of active sites for the hydrogen evolution reaction, edge sites and sulfur vacancies. Sulfur vacancies are the major active sites of the film, while edge sites are the major active sites of the flakes.18 Should the observed dramatic improvement in catalytic activity be due to the generation of new active sites, we would expect a very substantial change in the number of edge sites or sulfur

ACS Paragon Plus Environment

6

Page 7 of 22

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

Journal of the American Chemical Society

vacancies, which could be realized by inducing cracks or lose of sulfur atoms. However, we did not observe any cracks in the film or flakes after the cycling from optical, SEM, and AFM measurements (Fig. S5-S6). We also did not find any substantial change in the stoichiometric ratio of Mo and S from XPS measurements. (Fig. 2a). These results indicate no substantial change in the number of edge sites and sulfur vacancies during the cycling. This conclusion is also supported by the very poor catalytic activity of the MoS2 films grown with little sulfur vacancies (Fig. 1c) and by the reasonable linear dependence of the exchange current density of the flakes on the length of the edge (Fig. S7), both of which indicate that the catalytic activity is eventually dictated by the as-grown active sites. Additionally, the Raman spectra of MoS2 show negligible change after the cycling, and this excludes out the possibility of phase change from 2H to 1T in the MoS2, which could substantially improve the catalytic activity and induce obvious changes in the Raman spectra.10-12, 26 Therefore, we may conclude that the improved catalytic performance during the cycling is not due to the generation of new active sites but due to catalytic activity enhancement of the as-grown active sites, which we refer as activation. To better understand the activation mechanism, we examine the PL, Raman, and XPS of monolayer MoS2 films and flakes before and after the cycling. All of the measurement results point toward that the cycling process may induce p-doping to MoS2. The PL of the MoS2 shows obviously stronger and blueshifts after the cycling (Fig. 2c). This may be ascribed to p-doping, which has previously been demonstrated able to enhance PL intensity and induce PL blueshift for monolayer MoS2.27 Monolayer MoS2 is intrinsically n-doped and p-doping may neutralize the free electrons to promote the emission of neutral excitons over trions (negatively charged

ACS Paragon Plus Environment

7

Journal of the American Chemical Society

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 22

excitons), the former showing stronger PL efficiency and shorter emission wavelengths than the latter (Fig. S9).27-29 Additionally, the A1g Raman peak blueshifts by ~ 0.5 cm-1, but no visible change may be found at the E12g peak after the cycling. This supports the notion of p-doping, as it is well known that the A1g peak of MoS2 tends to stiffen (blueshift) with p-doping due to strong susceptibility to the influence of electron-phonon coupling but the E12g peak may not be affected by doping.30 The p-doping is also supported by the XPS measurement, which shows a shift of around -0.14 eV in the S 2p and Mo 3d peaks after the cycling. This shift in binding energy is consistent with what was previously observed at p-doped MoS2.31-33

Fig.3. Negligible effect of the electrochemical cycling in neutral electrolytes. (a) Polarization curves, (b) Raman spectra, and (c) PL spectra of monolayer MoS2 films before and after 5000 cycles in 0.5 M solution of Na2SO4. The other cycling condition is identical to what used for the cycling of MoS2 in 0.5 M solution of H2SO4. Note that the polarization curves were collected in 0.5 M solution of H2SO4 while the cycling performed in a separate solution of 0.5 M Na2SO4. We may correlate the p-doping and catalytic activity enhancement of MoS2 to protons, more specifically, to the protons intercalated between the monolayer and underlying substrates. We performed electrochemical cycling at monolayer MoS2 in the solution of 0.5M Na2SO4 and found negligible improvement in the catalytic performance and trivial p-doping effect (Fig. 3).

ACS Paragon Plus Environment

8

Page 9 of 22

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

Journal of the American Chemical Society

This result indicates that protons play a dominant role in the p-doping effect and catalytic activity enhancement observed at the MoS2 cycled in the acidic electrolyte, in which the only difference from the 0.5M Na2SO4 solution lies in the cation. This is consistent with the result of our previous work. We previously demonstrated from both experimental and theoretical grounds that protons may p-dope monolayer MoS2.34 The previous study also indicates that the p-doping is mainly contributed by the protons intercalated between the monolayer and underlying substrates, while the protons adsorbed on the top of the monolayer may only have trivial doping effect.34 Therefore, the observed p-doping effect in the MoS2 strongly suggests that the intercalation of protons between the monolayer and underlying substrates as illustrated in Fig. 2d is enabled during the cycling. Our other experimental observation may provide more support or insight for the proton intercalation and its correlation with the improvement of catalytic activity. The catalytic activity and p-doping effect gradually increase with the cycling. This indicates that the protons adsorbed on top of the monolayer, which is expected to happen immediately once the monolayer is immersed into the acidic electrolyte, only have negligible contribution. Additionally, the micrometer-scale flakes request much less cycles to achieve stable performance than the centimeter-scale films (Fig.1a-b), and extending the cycling range down to more negative voltages may decrease the number of cycles needed to reach stable catalytic activity (Fig. S4). This indicates that the intercalation of protons, which is hydrated in the acidic electrolyte, between the monolayer and underlying substrates is not spontaneous and requests driving force like negative electrochemical potentials to overcome certain energy barrier. To further confirm the role of proton intercalation in the improvement of catalytic activity, we examine the HER catalysis of monolayer MoS2 (films and flakes) treated with TFSI solution (1,2-dichlorobenzene and 1,2-dichloroethane used as solvent). We immersed monolayer

ACS Paragon Plus Environment

9

Journal of the American Chemical Society

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 10 of 22

MoS2 into the TFSI solution for minutes34-35 and then took it out for linear sweep voltammetry in acidic electrolytes. No electrochemical cycling is involved between the different linear sweeps. According to our previous studies, the immersion into TFSI solution may lead to the intercalation of protons between monolayer MoS2 and underlying substrates.34 The spontaneous protonintercalation solution could be due to the smaller size of non-hydrated protons in the TFSI solution , which may make the intercalation easier. We confirm the successful intercalation of protons after the TFSI treatment with PL, Raman, and XPS measurements, which shows similar enhancement and/or shifts as observed at the cycling in acidic media ( Fig. 4a-c). We have also confirmed no obvious change in the composition, structure, and morphology at the monolayers during the TFSI treatment,34 excluding the possibility of generation of new active sites or phase change. The TFSI treatment may enable enhancement in the catalytic activity similar to what observed at the electrochemical cycling. The catalytic performance increases with the number of TSFI treatment and eventually turns to be stable (Fig. 4d). The stable catalytic activity caused by the TFSI treatment is very similar to what caused by the cycling (Fig. 4d). We may ascribe the activity improvement by the TFSI treatment to protons as well, because no major activity improvement may be observed at the MoS2 treated by Li-TFSI solutions (Fig. S10), which have the same species as the TFSI except the cations. The gradual increase of catalytic activity with the number of treatment also suggests that the protons adsorbed at the surface of the monolayer, which is expected to happen immediately, do not play any major role. All of these indicate similar activation mechanisms for both TFSI-treated and electrochemically cycled MoS2, i.e. intercalation of protons. It is worthwhile to point out that the PL of the TFSI-treated monolayer shows negligible change after the catalytic characterization. This indicates that the intercalated protons are not reduced to hydrogen atoms during the electrocatalytic reaction, which would

ACS Paragon Plus Environment

10

Page 11 of 22

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

Journal of the American Chemical Society

otherwise lead to much weaker PL intensity as hydrogen atoms are known able to n-dope the monolayer.36-37

Fig. 4. Effect of TFSI treatment on monolayer MoS2. (a) Polarization curves of monolayer MoS2 films with different times of TFSI treatment. Each TFSI treatment is done at 100 ºC for 10mins. The dashed line is the polarization curve collected from a similar film experienced extensive electrochemical cycling. (b) PL spectra of monolayer MoS2 films before (black curve) and after (yellow curve) TFSI treatment. The red dotted curve is the PL spectrum collected from the TFSI-treated MoS2 film experienced multiple linear sweeps. (c) Raman and (d) XPS spectra of monolayer MoS2 films before (black curve) and after (yellow curve) TFSI treatment. The dashed lines in (c) indicate the positions of the Raman peaks before the TFSI treatment. Similar results for the TFSI-treated monolayer MoS2 flakes are given in Fig. S11. Except the space between monolayer and underlying substrates, protons may also be enabled to intercalate in the interlayer space of thicker MoS2 (Fig. 5a). To illustrate this notion, we examine the catalytic performance of bilayer MoS2 films. Similar to monolayer MoS2, the catalytic activity of bilayer MoS2 substantially increases with the electrochemical cycling in acidic electrolyte and eventually turns to be stable (Fig. 5b). We believe that this improvement in catalytic activity is collectively contributed by the protons intercalated underneath the bilayer and in the interlayer space. The proton intercalation in the interlayer space is evidenced by the decrease of the frequency difference between the A1g and E12g peaks after the cycling (Fig. 5c). The Raman frequency difference is related with the coupling of neighboring layers, and the decrease in the frequency difference indicates a decrease of interlayer coupling in the bilayer,38-39 which suggests the intercalation of protons into the interlayer space. The PL of the bilayer

ACS Paragon Plus Environment

11

Journal of the American Chemical Society

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 12 of 22

blueshifts and shows increase in efficiency, consistent with the expectation for reduced interlayer coupling (Fig. 5d).40-41 To further confirm the interlayer intercalation of protons, we examine the effect of TFSI treatment on the catalytic activity of bilayer films. Our previous studies have demonstrated that TFSI treatment may lead to proton intercalation between MoS2 and underlying substrates, but less efficient in enabling the proton intercalation into the interlayer space of multilayer 2D materials due to the stronger interaction between 2D layers.34 Therefore, the TFSI treatment is expected to mainly induce proton intercalation between bilayer and the substrates. Indeed, while both TFSI treatment and electrochemical cycling lead to similar catalytic activity improvement at monolayer MoS2 (Fig. 4a), the electrochemical cycling may enable greater activity improvement at bilayer MoS2 than the TFSI treatment (Fig. 5b). The TFSI treatment also gives rise to smaller changes in the Raman and PL of bilayer MoS2 than the electrochemical cycling (Fig. 5c-d). By electrochemically cycling the TFSI-treated bilayer in acidic electrolytes may further improve the catalytic activity. All these results indicate that the electrochemical cycling may enable proton intercalation into the interlayer space of bilayer MoS2 (Fig. 5a).

Fig. 5. Proton intercalation at bilayer MoS2. (a) Schematic illustration for the different proton intercalation configuration at the bilayer MoS2 experiencing electrochemically cycling and TFSI treatment. (b) Polarization curves of bilayer MoS2 as a function of electrochemical cycling (solid curves) and TFSI treatment (dashed curves). Typical initial polarization curve of bilayer MoS2 with neither electrochemical cycling nor TSFI treatment is also plotted (black curve). (c) Raman

ACS Paragon Plus Environment

12

Page 13 of 22

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

Journal of the American Chemical Society

and (d) PL spectra of bilayer MoS2 with no treatment (black), electrochemically cycled 3000 times (blue), and treated with TFSI solution (brown). The interlayer intercalation of protons and its effect on the catalytic activity may also be observed at even thicker MoS2 (Fig. 6). We synthesized MoS2 sheets in thickness of tens of nanometers and lateral size of micrometers on glassy carbon substrates following a process we developed previously,8 and then examined its catalytic performance as a function of electrochemical cycling. Similar to monolayer and bilayer MoS2, the catalytic activity of the MoS2 nanosheets substantially increases with the electrochemical cycling in acidic electrolytes, and eventually tends to be stable (Fig. 6a). And smaller nanosheets usually request less cycles to reach stable performance. We have confirmed no observerable change in the crystalline structures and composition of the nanosheets (Fig. 6b). However, XRD measurement shows that the (002) peak slightly shift toward smaller angle and shows obvious decrease in the intensity (Fig. 6c). Quantitative analysis indicates that the interlayer spacing increases by around 0.01 nm on average. Based on what we have discussed, it is reasonable to ascribe this increase in interlayer spacing to the intercalation of protons.

Fig.6. Proton intercalation at thick MoS2 nanosheets. (a) Polarization curves of MoS2 nanosheets as a function of electrochemical cycling in acidic electrolyte. The number of cycles is given as shown. Inset, SEM image for a typical MoS2 nanosheets (b) Raman spectra, and (c) XRD spectra of the MoS2 nanosheets before and after the electrochemical cycling.

ACS Paragon Plus Environment

13

Journal of the American Chemical Society

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

Page 14 of 22

Fig. 7. Activation of HER catalysis at MoS2 in neutral and alkaline environment. (a)Polarization curves of monolayer MoS2 films without activation (black), with activation (red), and cycled for 2000 times after the activation (blue) at neutral solution pH = 7. (b)Polarization curves of monolayer MoS2 films without activation (black), with activation (red), and cycled for 2000 times after the activation (blue) at alkaline solution pH = 12. The activation process is performing cyclic voltammetry at the MoS2 film in acidic electrolyte till stable catalytic activity.

We believe that the enhancement effect of the intercalated protons is an important reason for the worse catalytic performance of MoS2 in neutral and alkaline environment than in acidic media as reported previously.16-17 The electrochemical cycling in neutral or alkaline electrolyte may not give rise to any substantial increase in the catalytic performance of MoS2, as there is no enough hydrogen ions available for intercalation (Fig. S12). However, the catalytic activity of MoS2 in neutral and alkaline environment may also be substantially improved by the proton intercalation. We activate MoS2 by performing electrochemical cycling in acidic electrolytes first, and then transfer the monolayer in neutral or alkaline electrolytes for evaluation of catalytic activity. The activated MoS2 monolayer shows much higher catalytic performance than the counterpart without being activated (Fig. 7). The activated MoS2 also shows remarkable stability in neutral and alkaline solutions with no obvious decrease in catalytic activity over > 2000 cycles.

ACS Paragon Plus Environment

14

Page 15 of 22

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

Journal of the American Chemical Society

The stable enhancement effect results from the persistent presence of intercalated protons even in alkaline electrolyte, as evidenced by the negligible change in the Raman and PL of the MoS2 during the HER (Fig. S13).

There are two reasons for the intercalated protons not being

neutralized. First, it may be difficult for hydroxyl ions to intercalate into the confined space underneath monolayers or between the layers of thicker MoS2 as the size of hydroxyl ions is relatively big. Our experimental result suggests that the efficiency of the intercalation is inversely proportional to the size of the ions. For instance, the application of negative potential is necessary to drive the intercalation of hydrated protons in acidic electrolytes, but not necessary for the non-hydrated protons in TFSI, because hydrated protons may have larger sizes than nonhydrated ones. Additionally, we previously demonstrated that the cations with larger sizes like Na+ or K+ cannot intercalate underneath monolayer MoS2.34 Second, the intercalated protons might be chemically adsorbed onto the active sites, such as sulfur vacancies and edge sites, and this could prevent the intercalated protons to come out from the confined space. The intercalated protons may affect the catalytic activity through multiple effects. One effect lies in the charge transfer between the proton and MoS2, as indicated by the p-doping effect. The charge transfer may affect the catalytic activity by changing the electronic structure of the active site.42-44 Previous studies have demonstrated that the effect of charge transfer may be examined from the adsorption energy of hydrogen atoms.42-44 We perform first principle calculation for the adsorption free energy of hydrogen atoms, which is widely used as a figure of

ACS Paragon Plus Environment

15

Journal of the American Chemical Society

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 16 of 22

merit for the HER catalytic activity,4, 45 at the sulfur vacancies and edge sites of MoS2 with the presence of protons. The calculation indicates that the adsorption energy turns to be closer to 0 with the presence of H+ nearby as shown in Table 1 (see Fig. S13 and S14 for the structure used in the calculations). This indicates better HER catalytic activities at both sulfur vacancies and edge sites with the presence of protons, consistent with what we observed in experiments. It is worthwhile to point that the charge transfer between proton and MoS2 may induce p-doping effect, but this does not mean that any cations with p-doping effect would certainly be able to enhance the catalytic activity. For instance, our previous studies have demonstrated that Li-TFSI treatment may provide p-doping to monolayer MoS2,34 but we do not observe any improvement in the catalytic activity of MoS2 after Li-TFSI treatment (Fig. S10). Besides from the charge transfer effect, the intercalated protons might also affect the catalytic activity by increasing the electric conductance of MoS2 in both the vertical direction and horizontal directions. This is evidenced by the significant increase in exchange current density after the proton intercalation (Fig. S2 and S7), as recent studies have demonstrated that the exchange current density of MoS2 is related with the electrical conductance.15 More studies are necessary to further confirm and understand the effect of intercalated protons on the electrical conductance of MoS2. Table 1. Adsorption free energies of H atom Site

Without H+

With H+

Mo edge

0.100

0.044

S vacancy

-0.090

-0.052

CONCLUSION We have demonstrated that the intercalation of protons presents a universal strategy to significantly improve the HER catalytic performance of MoS2 for all kinds of environment from

ACS Paragon Plus Environment

16

Page 17 of 22

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

Journal of the American Chemical Society

acidic to alkaline. This strategy may be generally applied to MoS2 either with atomically thin dimension such as monolayer or bilayer or with nanostructured features such as nanosheets. Protons can be enabled to intercalate between MoS2 and substrates or into the interlayer space of MoS2 by two approaches, electrochemically polarizing MoS2 at negative potentials (vs. RHE) in acidic media or immersing MoS2 in certain acid solutions like TFSI. This improvement in catalytic performance result from the catalytic activity enhancement of the active sites in MoS2 by the intercalated protons. The intercalated protons may affect the catalytic activity by enabling more favorable adsorption energy of hydrogen atoms and enhancing the electrical conductance. More studies are necessary to obtain better understanding. This result represents significant advance in fundamental understanding for the HER catalysis, and points toward a new way for the development of cost-effective high performance catalysts for HER at all kinds of pH environments. While this study focuses on MoS2, protons, and hydrogen evolution, the notion of using intercalation of species to tune the catalytic performance may generally be applied to other ions or molecules, other layered materials, and other catalytic reactions.

ACS Paragon Plus Environment

17

Journal of the American Chemical Society

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 18 of 22

ACKNOWLEDGEMENTS This work was supported by CCDM, an EFRC funded by U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), under award #DE-SC0012575. Y.Z. acknowledges the support by DOE, Office of Science, BES, Materials Sciences and Engineering Division, under contract #DE-SC0012704. The authors acknowledge the use of the Analytical Instrumentation Facility (AIF) at North Carolina State University, which is supported by the State of North Carolina and the National Science Foundation. COMPETING FINANCIAL INTERESTS The authors declare no competing financial interests. Supporting Information The Supporting Information is available free of charge on the ACS Publications website at DOI: xxxxx. 1. 2. 3. 4. 5.

Effect of counter electrodes on the catalytic measurement Tafel plots of the catalytic characterization. Catalyic performance of MoS2 flakes with different sizes. Effect of electrochemical cycling at smaller voltage ranges. Structural, compositional, Raman, and PL characterization of MoS2 after electrochemical cycling or TFSI treatment. 6. Effect of Li-TFSI treatment. 7. Effect of electrochemical cycling in neutral and alkaline solutions. 8. Geometrical feature of the structures used in first principle calculation.

ACS Paragon Plus Environment

18

Page 19 of 22

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

Journal of the American Chemical Society

REFERENCE 1. Laursen, A. B.; Kegnaes, S.; Dahl, S.; Chorkendorff, I., Molybdenum sulfides-efficient and viable materials for electro - and photoelectrocatalytic hydrogen evolution. Energ Environ Sci 2012, 5 (2), 5577-5591. 2. Merki, D.; Hu, X. L., Recent developments of molybdenum and tungsten sulfides as hydrogen evolution catalysts. Energ Environ Sci 2011, 4, 3878-3888. 3. Gao, M. R.; Chan, M. K. Y.; Sun, Y. G., Edge-terminated molybdenum disulfide with a 9.4-angstrom interlayer spacing for electrochemical hydrogen production. Nat Commun 2015, 6. 4. Jaramillo, T. F.; Jorgensen, K. P.; Bonde, J.; Nielsen, J. H.; Horch, S.; Chorkendorff, I., Identification of active edge sites for electrochemical H-2 evolution from MoS2 nanocatalysts. Science 2007, 317 (5834), 100-102. 5. Karunadasa, H. I.; Montalvo, E.; Sun, Y. J.; Majda, M.; Long, J. R.; Chang, C. J., A Molecular MoS2 Edge Site Mimic for Catalytic Hydrogen Generation. Science 2012, 335 (6069), 698-702. 6. Kibsgaard, J.; Chen, Z. B.; Reinecke, B. N.; Jaramillo, T. F., Engineering the surface structure of MoS2 to preferentially expose active edge sites for electrocatalysis. Nat Mater 2012, 11 (11), 963-969. 7. Kong, D.; Wang, H.; Cha, J. J.; Pasta, M.; Koski, K. J.; Yao, J.; Cui, Y., Synthesis of MoS2 and MoSe2 Films with Vertically Aligned Layers. Nano Letters 2013, 13 (3), 1341-1347. 8. Li, Y.; Yu, Y.; Huang, Y.; Nielsen, R. A.; Goddard, W. A.; Li, Y.; Cao, L., Engineering the Composition and Crystallinity of Molybdenum Sulfide for High-Performance Electrocatalytic Hydrogen Evolution. ACS Catalysis 2015, 5 (1), 448-455. 9. Li, Y. G.; Wang, H. L.; Xie, L. M.; Liang, Y. Y.; Hong, G. S.; Dai, H. J., MoS2 Nanoparticles Grown on Graphene: An Advanced Catalyst for the Hydrogen Evolution Reaction. Journal of the American Chemical Society 2011, 133 (19), 7296-7299. 10. Lukowski, M. A.; Daniel, A. S.; Meng, F.; Forticaux, A.; Li, L.; Jin, S., Enhanced Hydrogen Evolution Catalysis from Chemically Exfoliated Metallic MoS2 Nanosheets. J. Am. Chem. Soc. 2013, 135, 10274–10277. 11. Voiry, D.; Yamaguchi, H.; Li, J.; Silva, R.; Alves, D. C. B.; Fujita, T.; Chen, M.; Asefa, T.; Shenoy, V. B.; Eda, G.; Chhowalla, M., Enhanced catalytic activity in strained chemically exfoliated WS2 nanosheets for hydrogen evolution. Nat. Mater. 2013, 12, doi:10.1038/nmat3700. 12. Wang, H. T.; Lu, Z. Y.; Xu, S. C.; Kong, D. S.; Cha, J. J.; Zheng, G. Y.; Hsu, P. C.; Yan, K.; Bradshaw, D.; Prinz, F. B.; Cui, Y., Electrochemical tuning of vertically aligned MoS2

ACS Paragon Plus Environment

19

Journal of the American Chemical Society

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 20 of 22

nanofilms and its application in improving hydrogen evolution reaction. P Natl Acad Sci USA 2013, 110 (49), 19701-19706. 13. Xie, J.; Zhang, J.; Li, S.; Grote, F.; Zhang, X.; Zhang, H.; Wang, R.; Lei, Y.; Pan, B.; Xie, Y., Controllable Disorder Engineering in Oxygen-Incorporated MoS2 Ultrathin Nanosheets for Efficient Hydrogen Evolution. Journal of the American Chemical Society 2013, 135 (47), 17881-17888. 14. Ye, G. L.; Gong, Y. J.; Lin, J. H.; Li, B.; He, Y. M.; Pantelides, S. T.; Zhou, W.; Vajtai, R.; Ajayan, P. M., Defects Engineered Monolayer MoS2 for Improved Hydrogen Evolution Reaction. Nano Letters 2016, 16 (2), 1097-1103. 15. Yu, Y.; Huang, S.-Y.; Li, Y.; Steinmann, S. N.; Yang, W.; Cao, L., Layer-Dependent Electrocatalysis of MoS2 for Hydrogen Evolution. Nano Letters 2014, 14 (2), 553-558. 16. Pua, Z.; Liua, Q.; Asirib, A. M.; Luo, Y.; Sun, X.; He, Y., 3D macroporous MoS2thin film: in situ hydrothermal preparation andapplication as a highly active hydrogen evolution electrocatalyst at all pH values. Electrochimica Acta 2015, 168, 133-138. 17. Xing, Z.; Yang, X.; Asiri, A. M.; Sun, X., Three-Dimensional Structures of MoS2@Ni Core/Shell Nanosheets Array toward Synergetic Electrocatalytic Water Splitting. ACS Appl. Mater. Interfaces 2016, 8 14521. 18. Guoqing Li; Du Zhang; Qiao Qiao; Yifei Yu; David Peterson; Abdullah Zafar; Raj Kumar; Stefano Curtarolo; Frank Hunte; Steve Shannon; Yimei Zhu; Weitao Yang; Cao, a. L., All The Catalytic Active Sites of MoS2 for Hydrogen Evolution. J. Am. Chem. Soc. 2016, 138, 16632−16638. 19. Yu, Y. F.; Li, C.; Liu, Y.; Su, L. Q.; Zhang, Y.; Cao, L. Y., Controlled Scalable Synthesis of Uniform, High-Quality Monolayer and Few-layer MoS2 Films. Sci Rep-Uk 2013, 3. 20. Yu, Y.; Hu, S.; Su, L.; Huang, L.; Liu, Y.; Jin, Z.; Purezky, A. A.; Geohegan, D. B.; Kim, K. W.; Zhang, Y., Equally Efficient Interlayer Exciton Relaxation and Improved Absorption in Epitaxial and Nonepitaxial MoS2/WS2 Heterostructures. Nano Letters 2014, 15 (1), 486491. 21. Gurarslan, A.; Yu, Y.; Su, L.; Yu, Y.; Suarez, F.; Yao, S.; Zhu, Y.; Ozturk, M.; Zhang, Y.; Cao, L., Surface-Energy-Assisted Perfect Transfer of Centimeter-Scale Monolayer and FewLayer MoS2 Films onto Arbitrary Substrates. ACS nano 2014, 8 (11), 11522-11528. 22. Kresse, G.; Furthmuller, J., Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys Rev B 1996, 54 (16), 11169-11186. 23. Perdew, J. P.; Burke, K.; Ernzerhof, M., Generalized gradient approximation made simple. Phys Rev Lett 1996, 77 (18), 3865-3868. 24. Grimme, S., Semiempirical GGA-type density functional constructed with a longrange dispersion correction. J Comput Chem 2006, 27 (15), 1787-1799.

ACS Paragon Plus Environment

20

Page 21 of 22

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

Journal of the American Chemical Society

25.

Tsai, C.; Li, H.; Park, S.; Park, J.; Han, H. S.; Norskov, J. K.; Zheng, X. L.; Abild-Pedersen,

F., Electrochemical generation of sulfur vacancies in the basal plane of MoS2 for hydrogen evolution. Nat Commun 2017, 8. 26. Voiry, D.; Salehi, M.; Silva, R.; Fujita, T.; Chen, M.; Asefa, T.; Shenoy, V. B.; Eda, G.; Chhowalla, M., Conducting MoS2 Nanosheets as Catalysts for Hydrogen Evolution Reaction. Nano Letters 2013, 13 (12), 6222-6227. 27. Yifei Yu; Yiling Yu; Chao Xu; Yong-Qing Cai; Liqin Su; Yong Zhang; Yong-Wei Zhang; Kenan Gundogdu; Cao, L., Engineering Substrate Interactions for High Luminescence Efficiency of Transition-Metal Dichalcogenide Monolayers. Adv. Func. Mater. 2016, DOI: 10.1002/adfm.201600418. 28. Mak, K. F.; He, K. L.; Lee, C.; Lee, G. H.; Hone, J.; Heinz, T. F.; Shan, J., Tightly bound trions in monolayer MoS2. Nat Mater 2013, 12 (3), 207-211. 29. Ross, J. S.; Wu, S. F.; Yu, H. Y.; Ghimire, N. J.; Jones, A. M.; Aivazian, G.; Yan, J. Q.; Mandrus, D. G.; Xiao, D.; Yao, W.; Xu, X. D., Electrical control of neutral and charged excitons in a monolayer semiconductor. Nat Commun 2013, 4. 30. Chakraborty, B.; Bera, A.; Muthu, D. V. S.; Bhowmick, S.; Waghmare, U. V.; Sood, A. K., Symmetry-dependent phonon renormalization in monolayer MoS2 transistor. Phys Rev B 2012, 85, 161403(R). 31. Mouri, S.; Miyauchi, Y.; Matsuda, K., Tunable photoluminescence of monolayer MoS2 via chemical doping. Nano Letters 2013, 13 (12), 5944-5948. 32. Nipane, A.; Karmakar, D.; Kaushik, N.; Karande, S.; Lodha, S., Few-Layer MoS2 p-Type Devices Enabled by Selective Doping Using Low Energy Phosphorus Implantation. Acs Nano 2016, 10 (2), 2128-2137. 33. Yang, L. M.; Majumdar, K.; Liu, H.; Du, Y. C.; Wu, H.; Hatzistergos, M.; Hung, P. Y.; Tieckelmann, R.; Tsai, W.; Hobbs, C.; Ye, P. D., Chloride Molecular Doping Technique on 2D Materials: WS2 and MoS2. Nano Letters 2014, 14 (11), 6275-6280. 34. Yu, Y.; Huang, L.; Li, G.; Barrette, A.; Cai, Y.-Q.; Yu, Y.; Gundogdu, K.; Zhang, Y.-W.; Cao, L., Enhancing Multifunctionalities of Transition Metal Dichalcogenide Monolayers via Cation Intercalation. ACS Nano 2017, 11, 9390-9396. 35. Amani, M.; Lien, D.-H.; Kiriya, D.; Xiao, J.; Azcatl, A.; Noh, J.; Madhvapathy, S. R.; Addou, R.; Santosh, K.; Dubey, M., Near-unity photoluminescence quantum yield in MoS2. Science 2015, 350 (6264), 1065-1068. 36. Shi, H.; Pan, H.; Zhang, Y.-W.; Yakobson, B. I., Strong ferromagnetism in hydrogenated monolayer MoS2tuned by strain. Phys. Rev. B 2013, 88, 205305.

ACS Paragon Plus Environment

21

Journal of the American Chemical Society

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

37.

Page 22 of 22

Cai, Y.; Zhaoqiang Bai; Pan, H.; Feng, Y. P.; Yakobson, B. I.; Zhang, Y.-W., Constructing

metallic nanoroads on a MoS2monolayer via hydrogenation. Nanoscale 2014, 6 (1691– 1697). 38. Lee, C.; Yan, H.; Brus, L. E.; Heinz, T. F.; Hone, J.; Ryu, S., Anomalous Lattice Vibrations of Single- and Few-Layer MoS2. Acs Nano 2010, 4 (5), 2695-2700. 39. Rice, C.; Young, R.; Zan, R.; Bangert, U.; Wolverson, D.; Georgiou, T.; Jalil, R.; Novoselov, K., Raman-scattering measurements and first-principles calculations of straininduced phonon shifts in monolayer MoS 2. Phys Rev B 2013, 87 (8), 081307. 40. Mak, K. F.; Lee, C.; Hone, J.; Shan, J.; Heinz, T. F., Atomically Thin MoS2: A New DirectGap Semiconductor. Phys Rev Lett 2010, 105 (13). 41. Splendiani, A.; Sun, L.; Zhang, Y. B.; Li, T. S.; Kim, J.; Chim, C. Y.; Galli, G.; Wang, F., Emerging Photoluminescence in Monolayer MoS2. Nano Letters 2010, 10 (4), 1271-1275. 42. Chen, W.; Santos, E. J. G.; Zhu, W. G.; Kaxiras, E.; Zhang, Z. Y., Tuning the Electronic and Chemical Properties of Monolayer MoS2 Adsorbed on Transition Metal Substrates. Nano Letters 2013, 13 (2), 509-514. 43. Tsai, C.; Abild-Pedersen, F.; Norskov, J. K., Tuning the MoS2 Edge-Site Activity for Hydrogen Evolution via Support Interactions. Nano Letters 2014, 14 (3), 1381-1387. 44. Zhong, H. X.; Quhe, R. G.; Wang, Y. Y.; Ni, Z. Y.; Ye, M.; Song, Z. G.; Pan, Y. Y.; Yang, J. B.; Yang, L.; Lei, M.; Shi, J. J.; Lu, J., Interfacial Properties of Monolayer and Bilayer MoS2 Contacts with Metals: Beyond the Energy Band Calculations. Sci Rep-Uk 2016, 6. 45. Hinnemann, B.; Moses, P. G.; Bonde, J.; Jorgensen, K. P.; Nielsen, J. H.; Horch, S.; Chorkendorff, I.; Norskov, J. K., Biornimetic hydrogen evolution: MoS2 nanoparticles as catalyst for hydrogen evolution. Journal of the American Chemical Society 2005, 127 (15), 5308-5309.

ACS Paragon Plus Environment

22