Effect of Intercalated Metals on the Electrocatalytic ... - ACS Publications

Nov 9, 2017 - Department of Chemistry, Temple University, 1901 North 13th Street, ... Institute for Computational Molecular Science, Temple University...
0 downloads 0 Views 7MB Size
This article is made available for a limited time sponsored by ACS under the ACS Free to Read License, which permits copying and redistribution of the article for non-commercial scholarly purposes.

Effect of Intercalated Metals on the Electrocatalytic Activity of 1T-MoS2 for the Hydrogen Evolution Reaction Nuwan H. Attanayake,†,§ Akila C. Thenuwara,†,§ Abhirup Patra,‡,§ Yaroslav V. Aulin,†,§ Thi M. Tran,† Himanshu Chakraborty,†,§,∥ Eric Borguet,†,§ Michael L. Klein,†,§,∥ John P. Perdew,‡,§ and Daniel R. Strongin*,†,§ †

Department of Chemistry, Temple University, 1901 North 13th Street, Philadelphia, Pennsylvania 19122, United States Department of Physics, Temple University 1925 North 12th Street, Philadelphia, Pennsylvania 19122, United States § Center for Computational Design of Functional Layered Materials (CCDM), Temple University, 1925 North 12th Street, Philadelphia, Pennsylvania 19122, United States ∥ Institute for Computational Molecular Science, Temple University, 1925 North 12th Street, Philadelphia, Pennsylvania 19122, United States

ACS Energy Lett. 2018.3:7-13. Downloaded from pubs.acs.org by 5.189.204.75 on 10/15/18. For personal use only.



S Supporting Information *

ABSTRACT: We show that intercalation of cations (Na+, Ca2+, Ni2+, and Co2+) into the interlayer region of 1T-MoS2 is an effective strategy to lower the overpotential for the hydrogen evolution reaction (HER). In acidic media the onset potential for 1TMoS2 with intercalated ions is lowered by ∼60 mV relative to that for pristine 1T-MoS2 (onset of ∼180 mV). Density functional theory (DFT) calculations show a lowering in the Gibbs free energy for H-adsorption (ΔGH) on these intercalated structures relative to intercalant-free 1T-MoS2. The DFT calculations suggest that Na+ intercalation results in a ΔGH close to zero. Consistent with calculation, experiments show that the intercalation of Na+ ions into the interlayer region of 1T-MoS2 results in the lowest overpotential for the HER.

T

Precious metals such as Pt are still the most efficient HER catalysts, but they are not ideal for mass scale H2 production via water electrolysis because of their scarcity and high cost. Toward the goal of economic feasibility there is a significant research effort designed to find efficient and earth abundant electrocatalytic materials that can catalyze the HER at the lowest possible overpotential. Two-dimensional (2D) transition-metal dichalcogenides (TMDs) have shown significant promise for the HER in acidic media.7−9 Within this material class, WS2 and MoS2 are the most studied.10−12 This research effort has been inspired in part by the excellent catalytic activity of nano- and micro-sized MoS2 structures in petroleum refining industries.13 TMDs are layered structures with van der Waals forces binding the sheets together. The activity of the TMDs vary depending on structure and morphology. For example, bulk MoS2 was thought to be inactive toward HER until nano structures were realized and found to be excellent HER catalysts.8,14,15 Three polymorphs of MoS2 have been identified. The polymorphic semiconducting 2H type with a trigonal

he release of CO2 into the atmosphere through anthropogenic activities, such as the combustion of fossil fuels, has severe consequences for both flora and fauna. A proposed solution to this growing problem would be to use green energy from renewable energy sources such as wind and solar.1−3 Capture and storage of solar energy as a renewable energy source would be one of the most important steps toward turning away from the use of fossil fuels. Storing solar energy in batteries, however, is not currently a plausible solution because of their low energy densities. Hence, there is significant interest in developing methods to store solar energy in the form of chemical fuels. In this context, splitting water into hydrogen and oxygen using solar energy would be one such method to harness solar energy and store it in chemical bonds (i.e., hydrogen).4−6 Solar-induced water splitting has generally been investigated in two ways: (1) the direct photochemical splitting of water over suitable semiconductor catalyst and (2) the use of photovoltaics to harness solar energy and to then split water with an electrolyzer.1 One crucial step in this latter effort is the hydrogen evolution reaction (HER) step: +



2H + 2e → H 2

Received: September 12, 2017 Accepted: November 9, 2017 Published: November 9, 2017

(1) © 2017 American Chemical Society

7

DOI: 10.1021/acsenergylett.7b00865 ACS Energy Lett. 2018, 3, 7−13

Letter

http://pubs.acs.org/journal/aelccp

Letter

ACS Energy Letters prismatic coordination to six sulfur atoms was shown to have active sites on the edges of under-coordinated Mo atoms. The basal planes of these materials were found in general to be inactive toward the HER,14,16 while edge sites and sulfur anion vacancies provided sites for increased HER activity. 17 Furthermore, it was shown that introduction of lattice strain in the layered sheets increased the HER activity.18 Also, of note is that the interplanar conductivity of the TMD is 2200 times less than the conductivity along the basal plane.10 The paucity of active sites and low out-of-plane conductivity are key issues that limit the electrocatalytic activity of the bulk 2H phase of MoS2 (i.e., 2H-MoS2) for the HER. Toward the goal of making MoS2 a more efficient HER catalyst, it has been shown that the 2H-MoS2 phase can be transformed into the polymorphic metallic 1T type MoS2 (via the intercalation of Li ions19−23) with octahedral coordination of S around the Mo. Lithium intercalation followed by exfoliation yields stable nanosheets (NS) of 1T-MoS2, with lattice distortions.24−26 These distorted metallic 1T-MoS2 NS show good electrochemical activity relative to the 2H phase on both the edge sites and the basal plane.26,27 The exfoliated NS bear a surface negative charge that enables a stable colloidal solution in water, although this charge has been proposed to reduce their electrocatalytic efficiency.7,24 A hypothesis tested in this research is that the electrocatalytic activity of the layered 1T-MoS2 NS for the HER can be enhanced by controlling its electronic and structural properties via the intercalation of metal cations. Prior studies have already shown examples where the electrocatalytic behavior for water splitting of particular layered materials is sensitive to the nature of metal/cations in the interlayer region.28−30 To test this hypothesis for 1T-MoS2 we determined the electrocatalytic HER activity of this layered material individually intercalated with Co2+, Ni2+, Ca2+, and Na+. Density functional theory (DFT) calculations were employed to gain insight into the Gibbs free energy of adsorption for hydrogen (ΔGH) on each of the intercalated samples. The value of ΔGH has been widely used as a descriptor for evaluating catalysts for HER.14,31−33 In particular, an improved catalytic activity is observed for ΔGH values close to zero. To prepare 1T-MoS2 samples we used a well-established method that yields lithium intercalated bulk 2H-MoS2 by using strong reducing agents such as LiBH4 and n-butyl lithium.19 The intercalated material, LixMoS2, was then exfoliated in water to produce metallic monolayer NS of MoS2.34 Prior studies have shown that these NS can be restacked by adding cations to a solution containing NS.35 To prepare intercalated samples for this study, solutions containing either Na+, Ca2+, Ni2+, or Co2+ ions were added separately to a solution containing 2 mg/mL 1T-MoS2 NS. This addition led to a precipitated phase. The solution was centrifuged, and the precipitate was washed three times with deionized (DI) water (see the Supporting Information). Samples individually prepared with Na2+, Ca2+, Ni2+, and Co2+ will be referred to hereafter as Na/1T-MoS2, Ca/1TMoS2, Ni/1T-MoS2, and Co/1T-MoS2, respectively. X-ray diffraction (XRD) was used to determine the interlayer spacing of the various samples. Figure 1 exhibits diffractograms for bulk 2H-MoS2 and 1T-MoS2 NS that were collected for analysis upon drying and 1T-MoS2 NS that immediately precipitated from solution after individual exposure to Co2+, Ni2+, Ca2+, and Na+. Inspection of the data shows that bulk 2HMoS2 exhibits a (002) reflection associated with an interlayer

Figure 1. Comparison of (001) and (002) XRD reflections of bulk 2H-MoS2, 1T-MoS2, and metal cation intercalated 1T-MoS2.

spacing of 6.36 Å.36,37 Exfoliated MoS2 sheets that stacked upon drying exhibit an interlayer spacing of 12.87 Å, due to a water double layer that persists in the interlayer region.37 The comparatively weak Bragg reflections associated with 1T-MoS2 NS (see Figures 1 and S2 for complete diffractograms) after exfoliation and drying suggest that this material likely has poor long-range order (relative to the other samples). An analysis of the XRD data does show that the interlayer distance is a function of the intercalated metal cation. On the basis of an analysis of the (001) reflection, the presence of Na+ results in the smallest interlayer distance (12.33 Å), while Ni2+ results in the greatest interlayer distance (12.73 Å). To confirm that the intercalated samples maintained the 1T phase, Raman spectroscopy was carried out (Figure 2). The

Figure 2. Raman spectra of 2H-MoS2 (black),1T-MoS2 (green), and metal ion intercalated 1T-MoS2 samples.

distinct phonon modes at 154 (J1), 219 (J2) and 327 cm−1 (J3), which are characteristic of metallic 1T-MoS2, are attributed to the formation of a super lattice structure.35 The presence of the 1T super lattice structure was also observed in selective area electron diffraction (SAED) patterns of both 1T-MoS2 and intercalated 1T-MoS2 samples (Figure S3).27,38,39 The broadening of 404 (A1g) and 380 cm−1 (E12g) modes and the suppressed intensity of E12g are also typical for the 1T phase. These results suggest that the metallic nature of the 1T NS is 8

DOI: 10.1021/acsenergylett.7b00865 ACS Energy Lett. 2018, 3, 7−13

Letter

ACS Energy Letters

Figure 3. SEM images of (a) 1T-MoS2 and (b) Co/1T-MoS2 and TEM images of exfoliated (c) 1T-MoS2 and (d) stacked Co/1T-MoS2 (refer to the Supporting Information for images of other samples).

indicating that redox chemistry between the cations and MoS2 did not occur during the intercalation process. Figure 3 shows scanning electron microscopy (SEM) images of 1T-MoS2 and Co/1T-MoS2 and transmission electron microscopy (TEM) images of the exfoliated 1T-MoS2 and flocculated samples of Co/1T-MoS2. The exfoliated 1T-MoS2 sheets are very thin, and upon flocculation, several sheets stack over one another as evident from the SEM images and TEM image contrast. SEM and TEM images of the other intercalated samples are shown in the Supporting Information (Figures S8 and S9). During the intercalation process we do not experimentally observe a difference in the sizes of the 1TMoS2 NS that stack to form the intercalated structures. Electrochemical experiments were carried out to understand the effect of the intercalated metals on the HER activity of 1TMoS2. Measurements were performed with a three-electrode system (see the Supporting Information for further information). Figure 4A exhibits polarization curves for the different samples. These data show the effect of the intercalated atoms on the overpotential (η) for HER associated with the different electrocatalysts (at a current density of 10 mA/cm2). The lowest overpotential in the MoS2 series studied was associated with Na/1T-MoS2 (η = 183 mV), and this material had an onset potential as low as 120 mV. This η value was ∼50 mV lower than the η value for pristine 1T-MoS2. The η value of 230 mV vs RHE for 1T-MoS2 obtained in our study is consistent with prior HER studies for this same material.22,23,26 An important conclusion that can be drawn from the data is that all the 1T-MoS2 materials with intercalated metal ions show a lower η than does intercalant-free 1T-MoS2. While Na/1TMoS2 exhibits the lowest η value, all of the other intercalated metal ion-1T-MoS2 samples exhibit η values that are within 12

retained upon stacking in three dimensions with intercalated metal cation. Additional evidence for the presence of the 1T phase and the super lattice structure was obtained from X-ray photoelectron spectroscopy (XPS). The Mo 3d region associated with MoS2 consists of two primary peaks at ∼229 and 232 eV corresponding to the 3d3/2 and 3d5/2 components of Mo4+ (Figure S4). Similar to the analysis of XPS for 2H- and 1TMoS2 in prior work, deconvolution of the peaks obtained in this study for MoS2 shows additional peaks with a shift to lower binding energy of ∼1.0 eV, which corresponds to the 1T phase.25,34 Analysis of the S 2p1/2 and 2p3/2 XPS data also shows that each feature is composed of two contributions, due to the presence of the 2H and 1T phases of MoS2. The presence of two phases gives rise to the super lattice structure. Deconvolution of both the Mo 3d and S 2p regions shows that ∼75% of the exfoliated MoS2 material is in the 1T phase. The intercalated samples also show a similar contribution of the 1T phase suggesting that metal ion intercalation does not alter the relative concentration of the 2H and 1T phases (see Figures S4−S6). Analysis of the samples with energy dispersive spectroscopy (EDS) shows that the amount of intercalant varies in the range of 1−2% when 10 mM divalent cation solutions are used to prepare samples and 4−5 atomic % is observed using a 20 mM solution of monovalent Na+ cation (Figures S7 and S10). We did not observe an increase in the intercalated cation concentration when the MoS2 sheets were exposed to higher cation solution concentrations. We believe this result is because the exfoliated sheets flocculate only when the charge on the sheet is neutralized by the intercalated cations. XPS of the Ni and Co intercalated samples (Figure S6b,c) showed the presence of Ni2+ and Co2+, respectively, 9

DOI: 10.1021/acsenergylett.7b00865 ACS Energy Lett. 2018, 3, 7−13

Letter

ACS Energy Letters

Figure 4. (A) Polarization plots of current density (j) vs V after iR correction showing enhanced catalytic activity of metal cation intercalated 1T-MoS2. The inset exhibits a magnified view emphasizing the differences in overpotential between the samples (at −10 mA cm−2). (B) Tafel plots for 1T-MoS2 and metal cation intercalated 1T-MoS2 samples. (C) Nyquist plot showing the decrease in charge-transfer resistance after intercalation. (D) Plots of scan rate vs current density used for ECSA calculations.

in our research by zeta potential measurements (see Table S1). It was found, for example, that upon intercalation with Na+, the zeta potential decreased from a value of −43 to −29 mV. In view of this prior assertion, it is possible that the removal of some of the charge from the surface of 1T-MoS2 due to intercalation plays a role in facilitating the electron transfer between the NS and protons during the HER. To investigate the effect that intercalation had on the charge-transfer resistance (Rct) of the samples, we carried out electrochemical impedance studies. Figure 4C exhibits data obtained at a η value of 200 mV. 1T-MoS2 exhibits an Rct value of ∼60 Ω, compared to a Rct = ∼10 kΩ for the nonmetallic 2H-MoS2 phase (Figure S12). The lowest Rct values are observed for the intercalated 1T-MoS2 samples (Rct varies from 25 to 55 Ω), consistent with the superior HER activity of these samples. The electrochemically active surface area (ECSA) both gives an estimation of the active sites and is proportional to the double-layer capacitance (Cdl). The slope of capacitive current for cyclic voltamograms (Figure S13) plotted versus scan rate (Figure 4D) gives half the value of Cdl. An increase in Cdl is associated with an increase in ECSA.41 The Cdl values are tabulated in Table 1. On the basis of these values, all the intercalated samples display a higher Cdl. Higher Cdl values imply that the capacitance on the catalyst has improved, suggesting that intercalated cations reside in the interlayer region. The improvement in the ECSA upon intercalation suggests that more sites are accessible for HER upon intercalation. We do not believe that the improved ECSA is due to changes in interlayer spacing, because all the intercalated

mV of the Na-intercalated samples and more than 40 mV lower than the η associated with 1T-MoS2 (see Table 1). Table 1. Summary of Electrochemical Data sample 1T-MoS2 Na/1T-MoS2 Ca/1T-MoS2 Ni/1T-MoS2 Co/1T-MoS2

η (mV) Tafel slope (mV/dec) 230 183 190 191 195

45 45 46 47 46

Rct (Ω)

Cdl (mF cm−2)

60 25 30 30 40

3.75 9.20 7.35 6.30 6.15

Figure 4B exhibits Tafel plots, derived from extrapolating the linear region of a plot of η vs log j. The Tafel slope for Na/1TMoS2 is 45 mV per decade, similar to the Tafel slope for 1TMoS2 (also 45 mV decade−1). Tafel plots for all the samples are provided in the Supporting Information (Figure S11). On the basis of prior studies, the relatively low Tafel slope values obtained for our MoS2 samples imply that a two-electrontransfer process takes place in the HER, consistent with a Volmer−Heyrowsky mechanism where the rate-limiting step in the process is the desorption of H2.40 The similar Tafel slopes for all the samples suggest that the intercalation process does not alter the HER mechanism that is associated with pristine 1T-MoS2. It is also important to mention that prior studies26 have suggested that removing some negative charge from exfoliated 1T-MoS2 NS (by exposing to iodine) also improves the catalytic HER activity of 1T-MoS2 (based on a lowering of η). This reduction in charge upon intercalation is emphasized 10

DOI: 10.1021/acsenergylett.7b00865 ACS Energy Lett. 2018, 3, 7−13

Letter

ACS Energy Letters

effective HER catalysts than 1T-MoS2 without intercalant. These calculations suggest that the introduction of DOS at the EF, whether from a transition metal or sp metal, favorably affects HER by a lowering of ΔGH from the value associated with intercalant-free 1T-MoS2. To summarize, we have used a solution-based technique to intercalate metal cations into the interlayer region of 1T-MoS2 to improve the efficiency of the material for the HER. Intercalation likely lowers the η for the HER at least in part by increasing the ECSA of the intercalated samples. Also, DFT calculations suggest that the intercalation of either a sp or dband metal cation into the 1T phase of MoS2 results in a change of ΔGH, relative to the ΔGH associated with the intercalant-free material, that can be related to enhanced activity for the HER.

samples have similar spacings. This improved ECSA is an interesting experimental observation considering that prior research of the 2H phase of MoS2 suggested that an increase in accessible reaction sites and lowering of the η for this phase occurred with an expansion in interlayer distance.42 It may be the case that a reduced charge on the sheets upon cation intercalation may allow more reaction site accessibility for the proton reduction cycle, resulting in the decreased overpotential and increased ECSA. Chronopotentiometric measurements carried out at 10 mA/ cm2 show that the intercalated samples are stable over a period of 24 h (see Figure S14). Raman spectroscopy of Na/1T-MoS2 after the stability test shows that the 1T structure is preserved (Figure S15). All the experimental data taken together show that the intercalated samples are stable and more active than intercalant-free 1T-MoS2. Finally, it is useful to revisit the electrochemical polarization results in view of DFT calculations. Specifically, we calculated the hydrogen atom adsorption free-energy change (ΔGH) for 1T-MoS2 in the absence and presence of the different intercalated cations. Prior studies have generally shown that ΔGH is an excellent descriptor to assess the electrochemical HER activity of catalytic surfaces.14,31−33 In particular, optimal catalytic activity is associated with a ΔGH ∼ 0.14,32 DFT calculations were performed in the Vienna Ab Initio Simulation Package (VASP)43 within the projector augmented wave method (PAW).44 The newly developed meta-GGA SCAN45 and its nonlocal van der Waals corrected form, SCAN +rVV10,46 were used as the exchange−correlation approximations (refer to the Supporting Information for detailed information including geometries). Figure 5 summarizes the



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsenergylett.7b00865. Experimental conditions, DFT calculations and structures, Figures S1−S17 as cited in the text, and zeta potentials (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Nuwan H. Attanayake: 0000-0001-7622-2337 Akila C. Thenuwara: 0000-0002-6146-9238 Yaroslav V. Aulin: 0000-0002-7370-1241 Daniel R. Strongin: 0000-0002-1776-5574 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by the Center for the Computational Design of Functional Layered Materials, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Basic Energy Sciences under Award # DESC0012575. A.P. (who performed the density functional calculations) was supported by the National Science Foundation under Grant No, DMR-1607868 (J.P.P.).



Figure 5. Free energies of hydrogen atom adsorption on the basal plane of pristine 1T-MoS2 as well as on the intercalated 1T-MoS2.

ABBREVIATIONS TMDs, transition-metal dichalcogenides; HER, hydrogen evolution reaction; DFT, density functional theory; DOS, density of states

computational results. The largest ΔGH value is obtained for 1T-MoS2 in the absence of intercalant, while the Na + intercalated sample is associated with the ΔGH value closest to zero. These calculations are consistent with the experimental finding that the highest and lowest overpotential are associated with 1T-MoS2 and Na/1T-MoS2. Figure S14 exhibits density of states (DOS) plots for the different 1T-MoS2 systems. The presence of the localized d-states of Co and Ni at the Fermi level (EF) are prominent features for Ni/1T-MoS2 and Co/1TMoS2 while delocalized sp states characterize the Na/1T-MoS2 and Ca/1T-MoS 2 DOS plots. Interestingly, while the contributions of the different intercalants to the DOS of 1TMoS2 are varied, all these intercalated systems are more



REFERENCES

(1) Walter, M. G.; Warren, E. L.; McKone, J. R.; Boettcher, S. W.; Mi, Q.; Santori, E. A.; Lewis, N. S. Solar Water Splitting Cells. Chem. Rev. 2010, 110, 6446−6473. (2) Gray, H. B. Powering the planet with solar fuel. Nat. Chem. 2009, 1, 7−7. (3) Lewis, N. S.; Nocera, D. G. Powering the planet: Chemical challenges in solar energy utilization. Proc. Natl. Acad. Sci. U. S. A. 2006, 103, 15729−15735. (4) Bard, A. J.; Fox, M. A. Artificial Photosynthesis: Solar Splitting of Water to Hydrogen and Oxygen. Acc. Chem. Res. 1995, 28, 141−145.

11

DOI: 10.1021/acsenergylett.7b00865 ACS Energy Lett. 2018, 3, 7−13

Letter

ACS Energy Letters (5) Bockris, J. O’M.; Reddy, A. K. N. Modern Electrochemistry; Springer US: New York, 1970. (6) Haldane, J. B. S. Daedalus or Science and the Future, A paper read to the Heretics, Cambridge, on February 4th, 1923; transcript, 1993. (7) Chhowalla, M.; Shin, H. S.; Eda, G.; Li, L.-J.; Loh, K. P.; Zhang, H. The chemistry of two-dimensional layered transition metal dichalcogenide nanosheets. Nat. Chem. 2013, 5, 263−275. (8) He, Z.; Que, W. Molybdenum disulfide nanomaterials: Structures, properties, synthesis and recent progress on hydrogen evolution reaction. Appl. Mater. Today 2016, 3, 23−56. (9) Reale, F.; Sharda, K.; Mattevi, C. From bulk crystals to atomically thin layers of group VI-transition metal dichalcogenides vapour phase synthesis. Appl. Mater. Today 2016, 3, 11−22. (10) Benck, J. D.; Hellstern, T. R.; Kibsgaard, J.; Chakthranont, P.; Jaramillo, T. F. Catalyzing the Hydrogen Evolution Reaction (HER) with Molybdenum Sulfide Nanomaterials. ACS Catal. 2014, 4, 3957− 3971. (11) Merki, D.; Hu, X. Recent developments of molybdenum and tungsten sulfides as hydrogen evolution catalysts. Energy Environ. Sci. 2011, 4, 3878−3888. (12) Chua, X. J.; Luxa, J.; Eng, A. Y. S.; Tan, S. M.; Sofer, Z.; Pumera, M. Negative Electrocatalytic Effects of p-Doping Niobium and Tantalum on MoS2 and WS2 for the Hydrogen Evolution Reaction and Oxygen Reduction Reaction. ACS Catal. 2016, 6, 5724−5734. (13) Chianelli, R. R.; Siadati, M. H.; De la Rosa, M. P.; Berhault, G.; Wilcoxon, J. P.; Bearden, R.; Abrams, B. L. Catalytic Properties of Single Layers of Transition Metal Sulfide Catalytic Materials. Catal. Rev.: Sci. Eng. 2006, 48, 1−41. (14) Hinnemann, B.; Moses, P. G.; Bonde, J.; Jørgensen, K. P.; Nielsen, J. H.; Horch, S.; Chorkendorff, I.; Nørskov, J. K. Biomimetic Hydrogen Evolution: MoS2 Nanoparticles as Catalyst for Hydrogen Evolution. J. Am. Chem. Soc. 2005, 127, 5308−5309. (15) Chia, X.; Eng, A. Y. S.; Ambrosi, A.; Tan, S. M.; Pumera, M. Electrochemistry of Nanostructured Layered Transition-Metal Dichalcogenides. Chem. Rev. 2015, 115, 11941−11966. (16) Jaramillo, T. F.; Jørgensen, K. P.; Bonde, J.; Nielsen, J. H.; Horch, S.; Chorkendorff, I. Identification of Active Edge Sites for Electrochemical H2 Evolution from MoS2 Nanocatalysts. Science 2007, 317, 100−102. (17) Li, G.; Zhang, D.; Qiao, Q.; Yu, Y.; Peterson, D.; Zafar, A.; Kumar, R.; Curtarolo, S.; Hunte, F.; Shannon, S.; et al. All The Catalytic Active Sites of MoS2 for Hydrogen Evolution. J. Am. Chem. Soc. 2016, 138, 16632−16638. (18) Li, H.; Tsai, C.; Koh, A. L.; Cai, L.; Contryman, A. W.; Fragapane, A. H.; Zhao, J.; Han, H. S.; Manoharan, H. C.; AbildPedersen, F.; et al. Activating and optimizing MoS2 basal planes for hydrogen evolution through the formation of strained sulphur vacancies. Nat. Mater. 2016, 15, 48−53. (19) Murphy, D. W.; Di Salvo, F. J.; Hull, G. W.; Waszczak, J. V. Convenient preparation and physical properties of lithium intercalation compounds of Group 4B and 5B layered transition metal dichalcogenides. Inorg. Chem. 1976, 15, 17−21. (20) Ambrosi, A.; Sofer, Z.; Pumera, M. Lithium Intercalation Compound Dramatically Influences the Electrochemical Properties of Exfoliated MoS2. Small 2015, 11, 605−612. (21) Chua, C. K.; Loo, A. H.; Pumera, M. Top-Down and BottomUp Approaches in Engineering 1 T Phase Molybdenum Disulfide (MoS2): Towards Highly Catalytically Active Materials. Chem. - Eur. J. 2016, 22, 14336−14341. (22) Ambrosi, A.; Sofer, Z.; Pumera, M. 2H [rightward arrow] 1T phase transition and hydrogen evolution activity of MoS2, MoSe2, WS2 and WSe2 strongly depends on the MX2 composition. Chem. Commun. 2015, 51, 8450−8453. (23) 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.

(24) Acerce, M.; Voiry, D.; Chhowalla, M. Metallic 1T phase MoS2 nanosheets as supercapacitor electrode materials. Nat. Nanotechnol. 2015, 10, 313−318. (25) Eda, G.; Yamaguchi, H.; Voiry, D.; Fujita, T.; Chen, M.; Chhowalla, M. Photoluminescence from Chemically Exfoliated MoS2. Nano Lett. 2011, 11, 5111−5116. (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 Lett. 2013, 13, 6222−6227. (27) Yin, Y.; Han, J.; Zhang, Y.; Zhang, X.; Xu, P.; Yuan, Q.; Samad, L.; Wang, X.; Wang, Y.; Zhang, Z.; et al. Contributions of Phase, Sulfur Vacancies, and Edges to the Hydrogen Evolution Reaction Catalytic Activity of Porous Molybdenum Disulfide Nanosheets. J. Am. Chem. Soc. 2016, 138, 7965−7972. (28) Chen, Z.; Leng, K.; Zhao, X.; Malkhandi, S.; Tang, W.; Tian, B.; Dong, L.; Zheng, L.; Lin, M.; Yeo, B. S.; et al. Interface confined hydrogen evolution reaction in zero valent metal nanoparticlesintercalated molybdenum disulfide. Nat. Commun. 2017, 8, 14548. (29) Thenuwara, A. C.; Cerkez, E. B.; Shumlas, S. L.; Attanayake, N. H.; McKendry, I. G.; Frazer, L.; Borguet, E.; Kang, Q.; Remsing, R. C.; Klein, M. L.; et al. Nickel Confined in the Interlayer Region of Birnessite: an Active Electrocatalyst for Water Oxidation. Angew. Chem., Int. Ed. 2016, 55, 10381−10385. (30) Thenuwara, A. C.; Shumlas, S. L.; Attanayake, N. H.; Aulin, Y. V.; McKendry, I. G.; Qiao, Q.; Zhu, Y.; Borguet, E.; Zdilla, M. J.; Strongin, D. R. Intercalation of Cobalt into the Interlayer of Birnessite Improves Oxygen Evolution Catalysis. ACS Catal. 2016, 6, 7739− 7743. (31) Tang, Q.; Jiang, D.-e. Mechanism of Hydrogen Evolution Reaction on 1T-MoS2 from First Principles. ACS Catal. 2016, 6, 4953−4961. (32) Nørskov, J. K.; Bligaard, T.; Logadottir, A.; Kitchin, J. R.; Chen, J. G.; Pandelov, S.; Stimming, U. Trends in the Exchange Current for Hydrogen Evolution. J. Electrochem. Soc. 2005, 152, J23−J26. (33) Tsai, C.; Chan, K.; Nørskov, J. K.; Abild-Pedersen, F. Theoretical insights into the hydrogen evolution activity of layered transition metal dichalcogenides. Surf. Sci. 2015, 640, 133−140. (34) Papageorgopoulos, C. A.; Jaegermann, W. Li intercalation across and along the van der Waals surfaces of MoS2(0001). Surf. Sci. 1995, 338, 83−93. (35) Heising, J.; Kanatzidis, M. G. Exfoliated and Restacked MoS2 and WS2: Ionic or Neutral Species? Encapsulation and Ordering of Hard Electropositive Cations. J. Am. Chem. Soc. 1999, 121, 11720− 11732. (36) Divigalpitiya, W. M. R.; Frindt, R. F.; Morrison, S. R. Inclusion Systems of Organic Molecules in Restacked Single-Layer Molybdenum Disulfide. Science 1989, 246, 369−371. (37) Zubavichus, Y. V.; Slovokhotov, Y. L.; Schilling, P. J.; Tittsworth, R. C.; Golub, A. S.; Protzenko, G. A.; Novikov, Y. N. X-ray absorption fine structure study of the atomic and electronic structure of molybdenum disulfide intercalation compounds with transition metals. Inorg. Chim. Acta 1998, 280, 211−218. (38) Dungey, K. E.; Curtis, M. D.; Penner-Hahn, J. E. Structural Characterization and Thermal Stability of MoS2 Intercalation Compounds. Chem. Mater. 1998, 10, 2152−2161. (39) Heising, J.; Kanatzidis, M. G. Structure of Restacked MoS2 and WS2 Elucidated by Electron Crystallography. J. Am. Chem. Soc. 1999, 121, 638−643. (40) Li, Y.; Wang, H.; Xie, L.; Liang, Y.; Hong, G.; Dai, H. MoS2 Nanoparticles Grown on Graphene: An Advanced Catalyst for the Hydrogen Evolution Reaction. J. Am. Chem. Soc. 2011, 133, 7296− 7299. (41) McCrory, C. C. L.; Jung, S.; Ferrer, I. M.; Chatman, S. M.; Peters, J. C.; Jaramillo, T. F. Benchmarking Hydrogen Evolving Reaction and Oxygen Evolving Reaction Electrocatalysts for Solar Water Splitting Devices. J. Am. Chem. Soc. 2015, 137, 4347−4357. 12

DOI: 10.1021/acsenergylett.7b00865 ACS Energy Lett. 2018, 3, 7−13

Letter

ACS Energy Letters (42) Gao, M.-R.; Chan, M. K. Y.; Sun, Y. Edge-terminated molybdenum disulfide with a 9.4-Å interlayer spacing for electrochemical hydrogen production. Nat. Commun. 2015, 6, 7493. (43) Hafner, J. Ab-initio simulations of materials using VASP: Density-functional theory and beyond. J. Comput. Chem. 2008, 29, 2044−2078. (44) Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B: Condens. Matter Mater. Phys. 1994, 50, 17953−17979. (45) Sun, J.; Remsing, R. C.; Zhang, Y.; Sun, Z.; Ruzsinszky, A.; Peng, H.; Yang, Z.; Paul, A.; Waghmare, U.; Wu, X.; et al. Accurate firstprinciples structures and energies of diversely bonded systems from an efficient density functional. Nat. Chem. 2016, 8, 831−836. (46) Peng, H.; Yang, Z.-H.; Perdew, J. P.; Sun, J. Versatile van der Waals Density Functional Based on a Meta-Generalized Gradient Approximation. Phys. Rev. X 2016, 6, 041005.

13

DOI: 10.1021/acsenergylett.7b00865 ACS Energy Lett. 2018, 3, 7−13