Effect of Intercalated Metals on the Electrocatalytic Activity of 1T-MoS2

Nov 9, 2017 - Department of Chemistry, Temple University, 1901 North 13th Street, ..... (42) It may be the case that a reduced charge on the sheets up...
0 downloads 0 Views 1016KB Size
Subscriber access provided by READING UNIV

Letter

Effect of Intercalated Metals on the Electrocatalytic Activity of 1T-MoS for the Hydrogen Evolution Reaction 2

Nuwan H. Attanayake, Akila C. Thenuwara, Abhirup Patra, Yaroslav V. Aulin, Thi Tran, Himanshu Chakraborty, Eric Borguet, Michael L. Klein, John P Perdew, and Daniel R. Strongin ACS Energy Lett., Just Accepted Manuscript • DOI: 10.1021/acsenergylett.7b00865 • Publication Date (Web): 09 Nov 2017 Downloaded from http://pubs.acs.org on November 9, 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.

ACS Energy Letters 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 17

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 Energy Letters

Effect of Intercalated Metals on the Electrocatalytic Activity of 1TMoS2 for the Hydrogen Evolution Reaction. Nuwan H. Attanayake1,3, Akila C. Thenuwara1,3, Abhirup Patra2,3, Yaroslav V. Aulin1,3, Thi M. Tran1 Himanshu Chakraborty1,3,4, Eric Borguet1,3, Michael L. Klein 1,3,4, John P. Perdew2,3 and Daniel R. Strongin1,3,* [1]

Department of Chemistry, Temple University, 1901 N. 13th Street, Philadelphia, PA 19122 (USA)

[2]

Department of Physics, Temple University 1925 N. 12th Street, Philadelphia, PA 19122 (USA)

[3]

Center for Computational Design of Functional Layered Materials (CCDM)

[4]

Institute for Computational Molecular Science, Temple University, 1925 N. 12th Street, Philadelphia,

Pennsylvania 19122 (USA)

1

ABSTRACT: We show that intercalation of cations (Na+, Ca2+, Ni2+, and Co2+) into the interlayer region

2

of 1T-MoS2 is an effective strategy to lower the overpotential for the hydrogen evolution reaction (HER).

3

In acidic media the onset potential for 1T-MoS2 with intercalated ions is lowered by 60 mV relative to

4

pristine 1T-MoS2 (180 mV). Density functional theory (DFT) calculations show a lowering in the Gibbs

5

free energy for H-adsorption (GH) on these intercalated structures relative to intercalant-free 1T-MoS2.

6

The DFT calculations suggest that Na+ intercalation results in a GH close to zero. Consistent with

7

calculation, experiments show that the intercalation of Na+ ions into the interlayer region of 1T-MoS2

8

results in the lowest overpotential for the HER.

9 ACS Paragon Plus Environment

ACS Energy Letters

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

10

Page 2 of 17

TOC GRAPHICS

11 12

2

ACS Paragon Plus Environment

Page 3 of 17

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 Energy Letters

13

The release of CO2 into the atmosphere through anthropogenic activities, such as the combustion of

14

fossil fuels has severe consequences for both flora and fauna. A proposed solution to this growing problem

15

would be to use green energy from renewable energy sources such as wind and solar.1-3 Capture and

16

storage of solar energy as a renewable energy source would be one of the most important steps toward

17

turning away from the use of fossil fuels. Storing solar energy in batteries, however, is not currently a

18

plausible solution due to their low energy densities. Hence, there is significant interest in developing

19

methods to store solar energy in the form of chemical fuels. In this context, splitting water into hydrogen

20

and oxygen using solar energy would be one such method to harness solar energy and store it in chemical

21

bonds (i.e., hydrogen).4-6 Solar-induced water splitting has generally been investigated in two ways: (1)

22

the direct photochemical splitting of water over suitable semiconductor catalyst and (2) the use of photo-

23

voltaics to harness solar energy and to then split water with an electrolyzer.1 One crucial step in this latter

24

effort is the hydrogen evolution reaction (HER) step:

25 26

2H+ + 2e- → H2

(1)

27 28

Precious metals such as Pt are still the most efficient HER catalysts, but they are not ideal for mass scale

29

H2 production via water electrolysis due to their scarcity and high cost. Toward the goal of economic

30

feasibility there is a significant research effort designed to find efficient and earth abundant electrocata-

31

lytic materials that can catalyze the HER at the lowest possible overpotential. 2D transition metal dichal-

32

cogenides (TMD) have shown significant promise for the HER in acidic media.7-9 Within this material

33

class, WS2 and MoS2 are the most studied.10-12 This research effort has been inspired in part by the excel-

34

lent catalytic activity of nano and micro sized MoS2 structures in petroleum refining industries.13 TMD

35

are layered structures with van der Waals forces binding the sheets together. The activity of the TMD

3

ACS Paragon Plus Environment

ACS Energy Letters

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 17

36

varies depending on structure and morphology. For example, bulk MoS2 was thought to be inactive to-

37

wards HER until nano structures were realized and found to be excellent HER catalysts.8, 14-15 Three pol-

38

ymorphs of MoS2 have been identified. The polymorphic semiconducting 2H type with a trigonal pris-

39

matic coordination to six sulfur atoms was shown to have active sites on the edges on under coordinated

40

Mo atoms. The basal planes of these materials were found in general to be inactive towards HER,14, 16

41

while edge sites and sulfur anion vacancies provided sites for increased HER activity.17 Furthermore, it

42

was shown that introduction of lattice strain in the layered sheets increased the HER activity. 18 Also, of

43

note is that the interplanar conductivity of the TMD is 2200 times less than the conductivity along the

44

basal plane.10 The paucity of active sites and low out of plane conductivity are key issues that limit the

45

electrocatalytic activity of bulk 2H phase of MoS2 (i.e., 2H-MoS2) for the HER.

46

Toward the goal of making MoS2 a more efficient HER catalyst, it has been shown that the 2H-MoS2

47

phase can be transformed into the polymorphic metallic 1T type MoS 2 (via the intercalation of Li ions19-

48

23

49

stable nano-sheets (NS) of 1T-MoS2, with lattice distortions.24-26 These distorted metallic 1T-MoS2 NS

50

show good electrochemical activity relative to the 2H phase on both the edge sites and the basal plane.26-

51

27

52

hough this charge has been proposed to reduce their electrocatalytic efficiency.7, 24

) with octahedral coordination of S around the Mo. Lithium intercalation followed by exfoliation yields

The exfoliated NS bear a surface negative charge that enables a stable colloidal solution in water, alt-

53

A hypothesis tested in this research is that the electrocatalytic activity of the layered 1T-MoS2 NS for

54

the HER can be enhanced by controlling its electronic and structural properties via the intercalation of

55

metal cations. Prior studies have already shown examples where the electrocatalytic behavior for water

56

splitting of particular layered materials are sensitive to the nature of metal/cations in the interlayer re-

57

gion.28-30 To test this hypothesis for 1T-MoS2 we determined the electrocatalytic HER activity of this

58

layered material individually intercalated with Co2+, Ni2+, Ca2+, and Na+. Density Functional Theory

59

(DFT) calculations were employed to gain an insight into the Gibbs free energy of adsorption for hydrogen 4

ACS Paragon Plus Environment

Page 5 of 17

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 Energy Letters

60

(ΔGH) on each of the intercalated samples. The value of ΔGH has been widely used as a descriptor for

61

evaluating catalysts for HER.14, 31-33 In particular, an improved catalytic activity is observed for ΔGH val-

62

ues close to zero.

63

To prepare 1T-MoS2 samples we used a well-established method that yields lithium intercalated bulk

64

2H-MoS2 by using strong reducing agents such as LiBH4 and n-butyl lithium.19 The intercalated material,

65

LixMoS2, was then exfoliated in water to produce metallic monolayer NS of MoS2.34 Prior studies have

66

shown that these NS can be restacked by adding cations to a NS containing solution.35

67

To prepare intercalated samples for this study, solutions containing either Na+, Ca2+, Ni2+ or Co2+ ions

68

were added separately to a solution containing 2 mg/mL 1T-MoS2 NS. This addition led to a precipitated

69

phase. The solution was centrifuged and the precipitate was washed three times with deionized (DI) water

70

(see SI). Samples individually prepared with Na2+, Ca2+, Ni2+, and Co2+ will be referred to hereafter as

71

Na/1T-MoS2, Ca/1T-MoS2, Ni/1T-MoS2 and Co/1T-MoS2, respectively.

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

X-ray diffraction (XRD) was used to determine the interlayer spacing of the various samples. Figure 1

73

exhibits diffractograms for bulk 2H-MoS2 and 1T-MoS2 NS that were collected for analysis upon drying, 5

ACS Paragon Plus Environment

ACS Energy Letters

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 17

74

and 1T-MoS2 NS that immediately precipitated from solution after individual exposure to Co2+, Ni2+, Ca2+,

75

and Na+. Inspection of the data shows that bulk 2H-MoS2 exhibits a (002) reflection associated with an

76

interlayer spacing of 6.36 Ȧ.36-37 Exfoliated MoS2 sheets that stacked upon drying exhibit an interlayer

77

spacing of 12.87 Ȧ, due to a water double layer that persists in the interlayer region.37 The comparatively

78

weak Bragg reflections associated with 1T-MoS2 NS (Figure 1 and SI Figure S2 for complete diffracto-

79

grams) after exfoliation and drying suggest that this material likely has poor long range order (relative to

80

the other samples). An analysis of the XRD does show that the interlayer distance is a function of the

81

intercalated metal cation. Based on an analysis of the (001) reflection, the presence of Na+ results in the

82

smallest interlayer distance (12.33 Å), while Ni2+ results in the greatest interlayer distance (12.73 Å).

83

To confirm that the intercalated samples maintained the 1T phase, Raman spectroscopy was carried out

84

(Figure 2). The distinct phonon modes at 154 (J1), 219 (J2) and 327 cm-1 (J3), which are characteristic of

85

metallic 1T-MoS2, are attributed to the formation of a super lattice structure.35 The presence of the 1T

86

super lattice strucuture was also observed in selective area elecron diffraction (SAED) patterns of both

87

1T-MoS2 and intercalated 1T-MoS2 samples (SI Figure S3).27, 38-39 The broadening of 404 (A1g) and 380

88

cm-1 (E12g) modes, as well as suppressed intensity of E12g are also typical for 1T phase. These results

89

suggest that the metallic nature of the 1T NS is retained upon stacking in 3D with intercalated metal

90

cation.

6

ACS Paragon Plus Environment

Page 7 of 17

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 91 27 28 92 29 30 93 31 32 33 94 34 35 95 36 37 96 38 39 40 97 41 42 98 43 44 45 99 46 47100 48 49101 50 51 52102 53 54103 55 56 104 57 58 59 60

ACS Energy Letters

Figure 2. Raman spectra of 2H-MoS2 (black),1T-MoS2 (green), and metal ion intercalated 1TMoS2 samples. 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+ (SI Figure S4). Similar to the analysis of XPS for 2H- and 1T-MoS2 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, that corresponds to the 1T phase.25, 34 Analysis of the S 2p1/2 and 2p3/2 XPS 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 phase (see SI Figure S4-S6). Analysis of the samples with energy dispersive spectroscopy (EDS) suggests that the amount of intercalant varies in the range of 1-2% when 10 mM divalent cation solutions and 4-5 atomic % is observed with the 20 mM solution of monovalent Na+ cation (Figure S7 and S10). We did not observe an increase in the intercalated cation concentration when the 7

ACS Paragon Plus Environment

ACS Energy Letters

105

1 2 3 106 4 5 107 6 7 108 8 9 10109 11 12 110 13 14 15111 16 17112 18 19 113 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 17

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 and S6c) showed the presence of Ni2+ and Co2+, respectively, 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/1TMoS2. 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

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 SI for images of other samples)

8

ACS Paragon Plus Environment

Page 9 of 17

114

1 2 3 115 4 5 116 6 7 8 117 9 10118 11 12 119 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 Energy Letters

intercalated samples are shown in the SI (Figure S8-S9). During the intercalation process we do not experimentally observe a difference in the sizes of the 1T-MoS2 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 1T-MoS2. Measurements were performed with a three electrode system (see SI for further information). Figure 4A exhibits polarization curves for the different samples. These data show the effect

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.

9

ACS Paragon Plus Environment

ACS Energy Letters

120

1 2 3 121 4 5 122 6 7 8 123 9 10124 11 12 13125 14 15126 16 17 18127 19 20128 21 22 23129 24 25130 26 27 28131 29 30132 31 32 133 33 34 35134 36 37135 38 39 136 40 41 42137 43 44138 45 46 47139 48 49140 50 51 141 52 53 54142 55 56143 57 58 59 60

Page 10 of 17

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/1T-MoS2 exhibits the lowest  value, all of the other intercalated metal ion-1T- MoS2 samples exhibit  values that are within 12 mV of the Naintercalated samples and more than 40 mV lower than the  associated with 1T-MoS2 (see Table 1). Figure 4B exhibits Tafel plots, derived from extrapolating the linear region of a plot of  vs log j. The Tafel slope for Na/1T-MoS2 is 45 mV per decade, similar to the Tafel slope for 1T-MoS2 (also 45 mV decade-1). Tafel plots for all the samples are provided in the SI (Figure S11). Based on prior studies, the relatively low Tafel slope values obtained for our MoS2 samples imply that a two-electron-transfer 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 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 mV 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 10

ACS Paragon Plus Environment

Page 11 of 17

144

1 2 3 145 4 5 146 6 7 8 147 9 10148 11 12 13149 14 15150 16 17151 18 19 20152 21 22153 23 24 154 25 26 27155 28 29156 30 31 157 32 33 34158 35 36 159 37 38 39160 40 41161 42 43 44162 45 46163 47 48 164 49 50 51165 52 53 54 55 56 57 58 59 60

ACS Energy Letters

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 non-metallic 2H-MoS2 phase (figure SI 12). 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 SI13) plotted vs scan rate (Figure 4D) gives half the value of Cdl. The ECSA is proportional to the Cdl values (an increase in Cdl is associated with an increase in ECSA).41 The Cdl values are tabulated in Table 1. Based on 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, since all the intercalated 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 SI figure S14). Raman spectroscopy of Na/1T-MoS2 after the stability test shows that the 1T structure is preserved (SI figure S15). All the experimental data taken together show that the intercalated samples are stable and more active than intercalant-free 1T-MoS2.

11

ACS Paragon Plus Environment

ACS Energy Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24166 25 26 167 27 28 29168 30 31169 32 33 34170 35 36171 37 38 172 39 40 41173 42 43174 44 45 175 46 47 48176 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 17

Table 1: Summary of electrochemical data Sample

η (mV)

Tafel Rct (Ω) Slope (mV/dec)

Cdl (mF cm-2)

1TMoS2

230

45

60

3.75

Na/1TMoS2

183

45

25

9.20

Ca/1TMoS2

190

46

30

7.35

Ni/1TMoS2

191

47

30

6.30

Co/1TMoS2

195

46

40

6.15

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 non-local van der Waals corrected form, SCAN+rVV1046 were used as the exchange-correlation approximations (Refer to the SI for detailed information including geometries). Figure 5 summarizes the 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.

12

ACS Paragon Plus Environment

Page 13 of 17

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25177 26 27178 28 29179 30 31 32180 33 34181 35 36 182 37 38 39183 40 41184 42 43 44185 45 46186 47 48 49187 50 51188 52 53 189 54 55 56 57 58 59 60

ACS Energy Letters

Figure 5. Free energies of hydrogen atom adsorption on the basal plane of pristine 1TMoS2 as well as on the intercalated 1T-MoS2. 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/1T- MoS2 while delocalized sp states characterize the Na/1TMoS2 and Ca/1T- MoS2 DOS plots. Interestingly, while the contributions of the different intercalants to the DOS of 1T- MoS2 are varied, all these intercalated systems are more effective HER catalysts than 1TMoS2 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 d-band metal cation into the 1T phase of MoS2 results

13

ACS Paragon Plus Environment

ACS Energy Letters

190

1 2 3 191 4 5 6 192 7 8 193 9 10 11194 12 13 14195 15 16 17196 18 19 197 20 21 22198 23 24 25199 26 27200 28 29 30201 31 32202 33 34 203 35 36 37204 38 39205 40 41 42206 43 44 45207 46 47208 48 49 50 209 51 52 53 54210 55 56 57211 58 59 60

Page 14 of 17

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. ASSOCIATED CONTENT Supporting Information. Experimental Conditions, DFT calculations and structure. Figures S1-S17 as cited in the text. Zeta potentials. AUTHOR INFORMATION Corresponding Author * Daniel R. Strongin. Email: [email protected] Author Contributions All authors have given approval to the final version of the manuscript.

ACKNOWLWEDGEMENTS 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 # DE-SC0012575. A.P. was supported by the National Science Foundation under Grant No. DMR-1607868 (J.P.P.). ABBREVIATIONS TMD, Transition metal dichalcogenide; HER, Hydrogen Evolution Reaction; DFT, Density Functional Theory; DOS, Density of States.

REFERENCES 14

ACS Paragon Plus Environment

Page 15 of 17

212

1 213 2 3 214 4 215 5 216 6 217 7 218 8 9 219 10220 11221 12222 13223 14 15224 16225 17226 18227 19228 20 229 21 22230 23231 24232 25233 26234 27 28235 29236 30237 31238 32239 33 240 34 35241 36242 37243 38244 39245 40 41246 42247 43248 44249 45250 46 251 47 48252 49253 50254 51255 52 256 53 54257 55 56 57 58 59 60

ACS Energy Letters

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. 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. 5. John O'M Bockris; Reddy, A. K. N. Modern Electrochemistry.; Springers US: 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 twodimensional 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. App. Mater. Today 2016, 3, 23-56. 9. Reale, F.; Sharda, K.; Mattevi, C. From bulk crystals to atomically thin layers of group VItransition metal dichalcogenides vapour phase synthesis. App. 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. 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.; Abild-Pedersen, 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.

15

ACS Paragon Plus Environment

ACS Energy Letters

258

1 259 2 3 260 4 261 5 262 6 263 7 264 8 9 265 10266 11267 12268 13269 14 15270 16271 17272 18273 19274 20 275 21 22276 23277 24278 25279 26280 27 28281 29282 30283 31284 32285 33 286 34 35287 36288 37289 38290 39291 40 41292 42293 43294 44295 45296 46 297 47 48298 49299 50300 51301 52 302 53 54303 55304 56 57 58 59 60

Page 16 of 17

21. Chua, C. K.; Loo, A. H.; Pumera, M. Top-Down and Bottom-Up 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. Nano 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, 1038110385. 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. 16

ACS Paragon Plus Environment

Page 17 of 17

305

1 306 2 3 307 4 308 5 309 6 310 7 311 8 9 312 10313 11314 12315 13316 14 15317 16318 17319 18320 19321 20 322 21 22323 23324 24 25325 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 Energy Letters

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. Journal of the American Chemical Society 2015, 137, 4347-4357. 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 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 first-principles 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.

17

ACS Paragon Plus Environment