Nanowires Catalyze the Hydrogen Evolution ... - ACS Publications

Subscriber access provided by CORNELL UNIVERSITY LIBRARY

Article

Template-Grown MoS2 Nanowires Catalyze Hydrogen Evolution Reaction: Ultra-low Kinetic Barriers with High Active Site Density Chongyi Ling, Yixin Ouyang, Li Shi, Shijun Yuan, Qian Chen, and Jinlan Wang ACS Catal., Just Accepted Manuscript • Publication Date (Web): 14 Jun 2017 Downloaded from http://pubs.acs.org on June 14, 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 Catalysis 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 16

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 Catalysis

Template-Grown MoS2 Nanowires Catalyze Hydrogen Evolution Reaction: Ultra-low Kinetic Barriers with High Active Site Density Chongyi Ling1, Yixin Ouyang1, Li Shi1, Shijun Yuan1, Qian Chen1 and Jinlan Wang*1,2 1

School of Physics, Southeast University, Nanjing 211189, China

2

Synergetic Innovation Center for Quantum Effects and Applications (SICQEA), Hunan Normal

University, Changsha 410081, China

KEYWORDS: MoS2 nanowires; Volmer-Tafel mechanism; low barriers; hydrogen evolution reaction; electrochemical water splitting

ABSTRACT: Molybdenum disulfide (MoS2) is considered as one of the most promising lowcost catalyst for hydrogen evolution reaction (HER). So far, the limited active sites and high kinetic barriers for H2 evolution still impede its practical application in electrochemical water splitting. In this work, on basis of comprehensive first-principles calculations, we predict that the recently produced template-grown MoS2 nanowires (NWs) on Au(755) surfaces hold both ultralow kinetic barriers for H2 evolution and ultra-high active site density simultaneously. The calculated kinetic barrier of H2 evolution through the Tafel mechanism is only 0.49 eV on the Mo edges, making the Volmer-Tafel mechanism operative and Tafel slope can be as low as 30

ACS Paragon Plus Environment

1

ACS Catalysis

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 16

mV/decade. Through substituting Au(755) substrate with non-noble metals, such as Ni(755) and Cu(755), the activity will be maintained. This work provides a possible way to achieve the ultrahigh HER activity of MoS2 based catalysts.

Introduction Molybdenum disulfide (MoS2), as a promising low-cost catalyst for hydrogen evolution reaction (HER), has attracted enormous attention since 2005.1 However, because of the limited active site density and high kinetic barriers for H2 evolution, the catalytic performance of MoS2 is still far from the requirements to commercial application. Generally, there are two strategies to enhance the activity of a catalyst: i) increasing the active site density and ii) improving the intrinsic activity of each active site. As the active sites of MoS2 are limited at their open edges,2 increasing the active site density has become a widely used strategy, including exposing more edge sites through preparing defect-rich, mesoporous or vertically aligned MoS23-7, and activating the inert basal planes by introducing defects or doping8-10. On the contrary, enhancing the intrinsic activity has gained little attention. In fact, the available Tafel slope (b) limit for 2.303RT

MoS2 based catalysts is as high as 40 mV/decade ( b = (1+α)F ≈40 mV/dec), suggesting the Volmer-Heyrovsky mechanism is operative.11-19 Theoretical studies have also shown that the evolution of H2 through Tafel reaction on MoS2 based catalysts needs to conquer a relatively high kinetic barrier (~1.50 eV).20-22 In comparison, the so-far best HER catalyst, platinum (Pt), can catalyze HER through Volmer-Tafel mechanism with an energy barrier as low as ~0.80 eV, leading to the Tafel slope only 30 mV/decade ( b =

2.303RT 2F

≈30 mV/dec).23,

24

The high H2

evolution barriers of MoS2 based catalysts are believed to be the reason why they present poorer catalytic activity than Pt.25 Therefore, designing MoS2 based HER catalysts with low kinetic


2

Page 3 of 16

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 Catalysis

barrier for Tafel reaction to improve the intrinsic activity is highly desirable, which has not yet been arrived. Actually, the high kinetic barrier for H2 evolution through Tafel reaction (Ea) on MoS2 based catalysts originates from the large H-H distance between two neighboring adsorbed H atoms (dH-H).21 Decreasing the dH-H is thus expected to reduce the Ea and improve the intrinsic activity of MoS2. Very recently, MoS2 nanowire (NW) arrays have been successfully grown on an Au(755) surface template,26 where the produced NWs were aligned well along the substrate steps with high density. Moreover, through a careful analysis on the bonding characteristics, we found this structure contains a very short dH-H at its Mo edges. Meanwhile, the rather narrow width of ~0.67 nm signifies the ultra-high density of exposed edge sites of MoS2. Therefore, such a MoS2 nanostructure may present extremely low Ea for HER and have ultra-high active site density simultaneously. In this work, we systematically evaluate the HER performance of MoS2 NWs@Au(755), within the framework of density functional theory (DFT). It is found the kinetic barrier for H2 evolution on their Mo edges via the Tafel reaction is indeed greatly reduced as compared with other MoS2 based catalysts, which is as low as 0.49 eV. The ultra-low Ea endow MoS2 NWs@Au(755) with H2 evolution rates about 16 orders faster than other MoS2 based catalysts. The HER performance of MoS2 NWs@Cu(755) and Ni(755) are further investigated considering the high cost of Au, which also present ultra-high HER activity. Therefore, these ultra-narrow MoS2 NWs are expected to present the smallest Tafel slope and achieve greatest HER performance among MoS2 based catalysts ever reported. Computational Details


3

ACS Catalysis

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 16

The DFT calculations were performed through the projector augmented wave (PAW) method27 as implemented in the Vienna ab initio simulation package (VASP).28,

29

The

generalized gradient approximation in the Perdew-Burke-Ernzerhof (PBE) form30, 31 and a cutoff energy of 350 eV for plane-wave basis set were adopted. The convergence threshold was 10-5 eV and 0.03 eV/Å for energy and force, respectively. To avoid the interaction between two periodic units, a vacuum space at least 10 Å was used. The model of MoS2 NW@Au(755) was composed of a three layer Au(755) substrate with (1 × 10) surface unit cells and a MoS2 NW with 9 times of unit cell along the growth direction. For MoS2 NW@Ni(755) or MoS2 NW@Cu(755), the model was constructed by a three layer Ni(755) or Cu(755) substrate with (1 × 5) surface unit cells and a MoS2 NW with 4 times of unit cell along the growth direction (more details are included in Figure S1 in Supporting Information). To reduce the computational burden, the bottom two layers were fixed for all the three systems. DFT-D232 calculations were used to describe the weak interaction. The minimum energy pathway for H2 evolution is determined by using the climbing nudged elastic band method,33, 34 and each transition state is further confirmed by vibrational analysis. The HER performance was evaluated by calculating the reaction free energy for hydrogen adsorption (∆GH), which can be obtained by35, 36: ∆GH = ∆EH + ∆EZPE – T∆S. The ∆EH is the hydrogen adsorption energy, while the ∆EZPE and ∆S are the zero-point energy difference and the entropy difference between the adsorbed and gas phase, respectively. For the case of MoS2 NWs, the ∆EZPE – T∆S is calculated to be 0.29 eV. Results and Discussions To uncover how to reduce the dH-H of MoS2 based catalysts, analysis on bonding characters of various catalysts are carried out. Generally, the S atoms bond to other atoms by forming sp3 hybrid orbitals in MoS2. These four sp3 hybrid orbitals form a tetrahedral structure. On the


4

Page 5 of 16

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 Catalysis

basal planes of both 2H and 1T-MoS2, S atoms bond to three neighboring Mo atoms. Due to the directional properties of covalent bonds, the orientation of S-H bonds (OS-H) have only one direction, which is perpendicular to the S-S orientation (OS-S) as shown in Figure 1a. As a result, the shortest dH-H is equal to the distance between two adjacent S atoms, which exceeds 3.0 Å. For the case of Mo edge in MoS2 sheet, each S atom interacts with two Mo atoms. Thereby, the adsorbed H atoms have two orientations, both of which are perpendicular to the OS-S, resulting in two distinct structures as displayed in Figure 1b and 1c. Therefore, the shortest dH-H is still equal to the distance between two neighboring S atoms, larger than 3.0 Å. Under these circumstances, significant distortions of the surface geometry are needed to allow the dH-H to be small enough for a bond formation.21 The corresponding barriers are thus high (~1.50 eV) and the Tafel reaction is inhibited.


5

ACS Catalysis

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 16

Figure 1. Models of S atoms bond to other atoms in MoS2 (a) basal plane, (b) & (c) Mo edge as well as (d), (e) & (f) NWs grown on Au(755) surfaces. The cyan, yellow, golden and white balls represent the Mo, S, Au and H atoms, respectively. The cyan spindle structures describe the sp3 hybrid orbitals of S atoms. The red and blue arrows represent the orientations of S-S and S-H, respectively. For the newly synthesized MoS2 NW@Au(755), each S atom of Mo-edge interacts with a Mo atom and an Au atom when a H atom adsorbs. Thus, the adsorbed H atoms also have two orientations, resulting in three distinct structures (Figure 1d to 1f). However, unlike the former two cases where the OS-H is perpendicular to the OS-S, the OS-H in the Mo edge of NWs is nearly parallel to the OS-S. Accordingly, when two H atoms adsorb along different OS-H (Figure 1d), the H-H distance will be significantly reduced, which is about half of that in 2H-MoS2 Mo edge or 1T-MoS2 basal plane. Therefore, such a MoS2 nanostructure may possess a very low energy barrier for H2 evolution via Tafel reaction. Moreover, this template-grown NW presents ultranarrow width of 0.67nm, which provides abundant exposed edge sites as well. Therefore, we systematically study the HER performance of MoS2 NW@Au(755).

Figure 2. (a) Structure of MoS2 nanowire grown on Au(755) surface and labels of different sites for hydrogen adsorption. (b) Calculated free energy diagram for hydrogen evolution on different kinds of S sites. (c) Minimum-energy pathways of the Tafel reaction on SMo-edge sites and the


6

Page 7 of 16

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 Catalysis

insets are the initial state (IS), transition state (TS) and final state (FS). The golded, yellow and cyan balls represent the Au, S and Mo atoms, respectively. We first perform a qualitative description on the HER performance of MoS2 NW@Au(755). The equilibrium structure is presented in Figure 2a, and there are four kinds of sites for hydrogen adsorption, labeled as SMo-edge, Ssurf, SS-edge and SS-edge(Au), respectively. The calculated ∆GH of different sites are illustrated in Figure 2b, where the SMo-edge and SS-edge(Au) sites are highly active with ∆GH of 0.15 and 0.20 eV, respectively. Besides, the ∆GH of Ssurf or SS-edge site is either positive (0.66 eV) or relatively negative (-0.30 eV), indicating that both of them are not suitable for HER. These results are in accordance with that of 2H-MoS2 sheet, whose Mo edge, S edge and basal plane own modulate, relatively negative and too positive ∆GH, respectively.21 Besides, S atoms on both SMo-edge and SS-edge(Au) sites present metallic nature (Figure S2a and S2b), However, a question remains unresolved: can the intrinsic activity of MoS2 NW@Au(755) be largely improved as we expected? To answer this question, kinetic energy barriers (Ea) for H2 evolution on SMo-edge site of MoS2 NW@Au(755) is further explored. As shown Figure 2c, the Ea of Tafel reaction on SMo-edge of MoS2 NW@Au(755) is as low as 0.49 eV, which is not only greatly reduced as compared with that of other MoS2 based catalysts (~1.5 eV), but also much lower than that of Pt(111) surface (~0.8 eV). According to the Arrhenius formula (k = A·exp(Ea/RT)), the rate of H2 production on MoS2 NW@Au(755) SMo-edge site is estimated to be amazingly increased by about 9 orders faster than that of 1T MoS2 (catalyze HER through Volmer-Heyrovsky mechanism with a Heyrovsky reaction barrier of ~1.0 eV)20 at room temperature. It is even about 5 orders faster than that of Pt(111) with hydrogen coverage lower than one monolayer, whose barrier is about 0.8 eV.


7

ACS Catalysis

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 16

Figure 3. Structures of (a) MoS2 NW@Ni(755) and (b) MoS2 NW@Cu(755) as well as the labels of different sites for hydrogen adsorption. (c) Calculated free energy diagram for hydrogen evolution on different S sites of MoS2 NW@Ni(755) and MoS2 NW@Cu(755), respectively represented by solid and dash lines. (d) Calculated ∆GH as a function of number of electron S atom gains (Ne) for MoS2 NW@Ni(755) (blue dots and line) and MoS2 NW@Cu(755) (red quadrangles and line), respectively. The R2 is respectively 0.99 and 0.98, indicating the high linear correlation between the Ne and ∆GH. The blue, bronze, yellow, cyan and white balls represent the Ni, Cu, S, Mo and H atoms, respectively. Despite the prominent HER performance, the high expense of Au makes such a catalyst commercially impracticable, alternative non-noble metal substrates for the growth of MoS2 NWs are very necessary. For the template growth of MoS2 NWs, the atomic staircases of Au(755) surface is crucial.26 Furthermore, the special Au-S (edge) bonds lead to the greatly reduced H-H distance between two adjacent adsorbed H on S edge, which finally results in the ultra-low energy barrier of Tafel reaction. As non-noble metals, Ni and Cu have similar crystal structure to Au (face-centered-cubic structure). Naturally, the (755) surfaces of Ni and Cu also contain periodical atomic staircases which are analogous to that of Au(755) surface. Moreover, Ni(755) and Cu(755) surfaces have been used as templates to prepare nanomaterials experimentally.37, 38 Moreover, we computed the binding energies of MoS2 NWs on different substrates (Details can be found in SI). As shown in Table S1, Eb for MoS2 NW@Au(755), Cu(755) and Ni(755) are -


8

Page 9 of 16

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 Catalysis

0.92, -2.55 and -3.97 eV, respectively, indicating that MoS2 NWs on Ni(755) and Cu(755) is more stable than that on Au(755). Therefore, the experimental synthesis of MoS2 NWs on Ni(755) and Cu(755) surfaces and their excellent performance for HER can be expected.

Figure 4. Minimum-energy pathways of the Tafel reaction on different active sites: (a) or (b) SSedge

and (c) or (d) SMo-edge sites of MoS2 NW@Cu(755) or MoS2 NW@Ni(755), respectively. The

insets are the corresponding structures of initial state (IS), transition state (TS) and final state (FS). The blue, bronze, yellow and cyan balls represent the Ni, Cu, S and Mo atoms, respectively. Similar to MoS2 NW@Au(755), both MoS2 NW@Ni(755) and MoS2 NW@Cu(755) contain four different kinds of S sites for hydrogen adsorption: SMo-edge, Ssurf, SS-edge and SSedge(Ni/Cu)

as shown in Figure 3a and 3b. For MoS2 NW@Ni(755), both the edge sites, SMo-edge and

SS-edge present ultra-high HER activity with the calculated ∆GH of 0.19 and -0.12 eV, respectively. Other two kinds of sites (Ssurf and SS-edge(Ni)) have relatively positive ∆GH (0.83 and


9

ACS Catalysis

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 16

0.54 eV, respectively), therefore they are not active sites for HER (Figure 3c). Similar tendency has also been observed for MoS2 NW@Cu(755): SMo-edge and SS-edge are highly active, while Ssurf and SS-edge(Cu) are inert for HER. The calculated ∆GH for these four kinds of sites are 0.01, -0.16, 0.69 and 0.42 eV, respectively (see Figure 3c). The different HER activities of different sites can be ascribed to the different number of electrons that S atom gains (Ne). As shown in Figure 3d, the calculated ∆GH present nearly perfect linear correlation with the Ne of corresponding S atoms and a larger Ne will lead to a more positive ∆GH, which is in accordance with our previous studies.39, 40 Similar to MoS2 NW@Au(755), MoS2 NW@Ni(755) and Cu(755) can also ensure efficient charge transfer due to the metallic nature of S atoms on the active sites (Figure S2c to S2f). The kinetic energy barriers for H2 evolution on both kinds of edge sites of MoS2 NW@Ni(755) and MoS2 NW@Cu(755) are further explored. We first calculate the free energy diagrams of the second H* adjacent to the first one. As displayed in Figure S3, ∆GH tends to increase upon the adsorption of the second H, consistent with previous studies.35, 40 Except for SMo-edge site of NW@Ni(755) whose ∆GH of the second H adsorption is a little high, the rest three kinds of sites still possess nearly zero ∆GH. We then compute the Tafel reaction barriers on these active sites. As shown in Figure 4a and 4b, the Ea of Tafel reaction for SS-edge of both NWs keep costly, 1.80 and 1.74 eV for MoS2 NW@Cu(755) and NW@Ni(755), respectively, which is similar to other MoS2 based catalysts and suggesting that the Tafel-Heyrovsky is operative. Nevertheless, similar to that of MoS2 NW@Au(755), the Ea of Tafel reaction on SMo-edge of NW@Cu(755) and Ni(755) are very low, just 0.64 and 0.48 eV, respectively, as shown in Figure 4c and 4d. To determine the underlying mechanism of HER on these two kinds of sites, the energy barriers of Volmer reaction and Heyrovsky reaction are also calculated referring to


10

Page 11 of 16

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 Catalysis

previous works.20, 24 20 H2O molecules are used to simulate the solution effect and the results are presented in Figure 5. Similar to the cases of Pt(111) and 1T MoS2,20,

24

the Ea of Volmer

reaction are very small for both NW@Ni(755) and Cu(755), just 0.17 and 0.06 eV, respectively (Figure 5a and 5b). However, for Heyrovsky reaction, the Ea are 1.29 and 1.31 eV for NWs@Ni(755) and Cu(755) respectively, demonstrating that Volmer-Tafel mechanism is operative and Tafel reaction is the rate-determining step. Correspondingly, the Tafel slope can be as low as 30 mV/decade according to the equation (b =

2.303RT 2F

≈30 mV/dec), which will be the

lowest record for MoS2 based HER catalysts. Therefore, not only the cost will be greatly reduced, but also the ultra-high HER performance can be maintained by substituting Au(755) with Ni(755) or Cu(755) as the templates for the growth of MoS2 NWs. The experimental realization of MoS2 NW@Cu(755) and NW@Ni(755) is thus desirable. Closer examination of the structures reveals that the ultra-low Ea of Tafel reaction for these NWs indeed originates from the greatly reduced H-H distance between two neighboring adsorbed H atoms on SMo-edge sites. As clearly displayed in Figure 6, the barrier increases linearly with the enlargement of H-H distance. For 2H-MoS2 Mo edge or 1T-MoS2 basal plane, the shortest dH-H is normally larger than 3.0 Å. To form the H-H bond, a significant distortion of the surface geometry is thus needed, leading to a high kinetic barrier for H2 evolution, e.g., ~1.5eV. For Mo edges of MoS2 NWs, the dH-H between two neighboring H atoms is only ~1.50 Å, respectively, nearly half of that on 2H-MoS2 Mo edge or 1T-MoS2 basal plane. Correspondingly, the Ea is dramatically reduced by nearly 1.0 eV as compared with that of other MoS2 based catalysts, making the H2 evolution through the Volmer-Tafel route possible, which eventually improves the HER activity. It should be noted that Ea and dH-H of different systems present close


11

ACS Catalysis

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 16

linear correlation, suggesting that dH-H may be served as a simple descriptor for intrinsic HER activity of catalysts.

Figure 5. Minimum-energy pathways of Volmer and Heyrovsky reaction on SMo-edge sites of MoS2 NW@Ni(755) (a and c, respectively) and NW@Cu(755) (b and d, respectively). The inserts are the corresponding initial state (IS), transition state (TS) and final state (FS). The blue, bronze, yellow, cyan, red and white balls and lines represent the Ni, Cu, S, Mo, O and H atoms, respectively.


12

Page 13 of 16

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 Catalysis

Figure 6. Calculated activation barriers (Ea) for H2 desorption through the Tafel mechanism as a function of the H-H distance (dH-H), where the R2 is 0.95, indicating the high linear correlation between the Ea and dH-H. The purple pentagon, red stars, blue triangles, orange dot and green quadrangle represent the cases of MoS2 NWs@Ni(755), NWs@Cu(755), Mo edge and Pt(111) surface, respectively. The detailed information of two middle points is presented in Figure S4. Conclusions In summary, we have predicted the newly template-grown MoS2 nanowire on Au(755) and its analogous MoS2 NW@Ni(755) and MoS2 NW@Cu(755), owns both abundant active sites and greatly improved intrinsic HER activity as compared to other MoS2 based catalysts. The free energy for hydrogen adsorption is close to ~ 0 eV on both edges and more importantly, the kinetic barrier of H2 evolution through the Tafel mechanism can be as low as 0.48 eV on the Mo edges. The ultra-low kinetic barrier of MoS2 NW@Au(755) and MoS2 NW@Ni(755) suggests that the H2 evolution rate is at least 16 orders faster than other MoS2 based catalysts and even 5 orders faster than that on the Pt(111) surface. The Volmer-Tafel route is thus possible and a Tafel slope as low as 30 mV/decade is expected for MoS2 based catalysts. The ultra-low kinetic


13

ACS Catalysis

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 16

barriers are originated from the greatly reduced H-H distance between two neighboring adsorbed hydrogen atoms, suggesting that H-H distance may be served as a simple descriptor for intrinsic HER activity of the electrochemical catalysts. This work suggests that template-grown of MoS2 nanowire catalysts provide a promising way to achieve the greatest HER performance of MoS2 based catalysts, which is even better than the noble metal catalyst Pt. ASSOCIATED CONTENT Supporting Information. The detail information of the models we used, the partial density of states of S atom p orbital at different sites of MoS2@Au(755), Cu(755) and Ni(755), free energy diagram for the adsorption of the second H adjacent to the first one and the minimum-energy pathways of the Tafel reaction on S edges of MoS2 NW@Cu(755) Ni(755) as well as the binding energies of MoS2 NWs with different substrates are included in the Supporting Information. AUTHOR INFORMATION Corresponding Author *Email: [email protected]

ACKNOWLEDGMENT This work is supported by the National Science Funding (Grants 21525311, 21373045), the Ministry of Science and Technology (Grant 2017YFA0204800), Jiangsu 333 project (BRA2016353), and the Scientific Research Foundation of Graduate School of Southeast University (YBJJ1670) in China. The authors thank the computational resources at the SEU and National Supercomputing Center in Tianjin.


14

Page 15 of 16

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 Catalysis

REFERENCES (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)

B. Hinnemann, P. G. Moses, J. Bonde, K. P. Jørgensen, J. H. Nielsen, S. Horch, I. Chorkendorff and J. K. Nørskov, J. Am. Chem. Soc., 2005, 127, 5308-5309. T. F. Jaramillo, K. P. Jørgensen, J. Bonde, J. H. Nielsen, S. Horch and I. Chorkendorff, Science, 2007, 317, 100-102. J. Kibsgaard, Z. Chen, B. N. Reinecke and T. F. Jaramillo, Nat. Mater., 2012, 11, 963969. J. Xie, H. Zhang, S. Li, R. Wang, X. Sun, M. Zhou, J. Zhou, X. W. Lou and Y. Xie, Adv. Mater., 2013, 25, 5807-5813. Y. Tan, P. Liu, L. Chen, W. Cong, Y. Ito, J. Han, X. Guo, Z. Tang, T. Fujita, A. Hirata and M. W. Chen, Adv. Mater., 2014, 26, 8023-8028. D. Kong, H. Wang, J. J. Cha, M. Pasta, K. J. Koski, J. Yao and Y. Cui, Nano Lett., 2013, 13, 1341-1347. J. Shi, D. Ma, G.-F. Han, Y. Zhang, Q. Ji, T. Gao, J. Sun, X. Song, C. Li and Y. Zhang, ACS Nano, 2014, 8, 10196-10204. H. Li, C. Tsai, A. L. Koh, L. Cai, A. W. Contryman, A. H. Fragapane, J. Zhao, H. S. Han, H. C. Manoharan and F. Abild-Pedersen, Nat. Mater., 2015, 15, 48-53. J. Deng, H. Li, J. Xiao, Y. Tu, D. Deng, H. Yang, H. Tian, J. Li, P. Ren and X. Bao, Energy Environ. Sci., 2015, 8, 1594-1601. G. Gao, Q. Sun and A. Du, J. Phys. Chem. C, 2016, 120, 16761-16766. Y. Li, H. Wang, L. Xie, Y. Liang, G. Hong and H. Dai, J. Am. Chem. Soc., 2011, 133, 7296-7299. D. Gopalakrishnan, D. Damien and M. M. Shaijumon, ACS Nano, 2014, 8, 5297-5303. Y. Yu, S. Y. Huang, Y. Li, S. N. Steinmann, W. Yang and L. Cao, Nano Lett., 2014, 14, 553-558. J. Kibsgaard, T. F. Jaramillo and F. Besenbacher, Nat. Chem., 2014, 6, 248-253. W. Xiao, P. Liu, J. Zhang, W. Song, Y. P. Feng, D. Gao and J. Ding, Adv. Energy. Mater., 2016, 10.1002/aenm.201602086. D. Deng, K. S. Novoselov, Q. Fu, N. Zheng, Z. Tian and X. Bao, Nat. Nanotechnol., 2016, 11, 218-230. Q. Lu, Y. Yu, Q. Ma, B. Chen and H. Zhang, Adv. Mater., 2015, 28, 1917-1933. Y. Yan, B. Xia, Z. Xu and X. Wang, ACS Catal, 2014, 4, 1693-1705. Q. Gong, L. Cheng, C. Liu, M. Zhang, Q. Feng, H. Ye, M. Zeng, L. Xie, Z. Liu and Y. Li, ACS Catal, 2015, 5, 2213-2219. Q. Tang and D.-e. Jiang, ACS Catal, 2016, 6, 4953-4961. C. Tsai, K. Chan, J. K. Nørskov and F. Abild-Pedersen, Surf. Sci., 2015, 640, 133-140. Y. Ouyang, C. Ling, Q. Chen, Z. Wang, L. Shi and J. Wang, Chem. Mater., 2016, 28, 4390-4396. B. Conway and L. Bai, J. Electroanal. Chem., 1986, 198, 149-175. E. Skulason, G. S. Karlberg, J. Rossmeisl, T. Bligaard, J. Greeley, H. Jonsson and J. K. Norskov, Phys. Chem. Chem. Phys., 2007, 9, 3241-3250. Z. W. Seh, J. Kibsgaard, C. F. Dickens, I. Chorkendorff, J. K. Norskov and T. F. Jaramillo, Science, 2017, 355. H. Xu, S. Liu, Z. Ding, S. J. Tan, K. M. Yam, Y. Bao, C. T. Nai, M. F. Ng, J. Lu, C. Zhang and K. P. Loh, Nat. Commun., 2016, 7, 12904.


15

ACS Catalysis

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

(27) (28) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40)

Page 16 of 16

P. E. Blöchl, Phys. Rev. B, 1994, 50, 17953-17979. G. Kresse and J. Furthmüller, Phys. Rev. B, 1996, 54, 11169-11186. G. Kresse and D. Joubert, Phys. Rev. B, 1999, 59, 1758-1775. J. P. Perdew, J. Chevary, S. Vosko, K. A. Jackson, M. R. Pederson, D. Singh and C. Fiolhais, Phys. Rev. B, 1992, 46, 6671-6687. J. P. Perdew and Y. Wang, Phys. Rev. B, 1992, 45, 13244-13249. S. Grimme, J. Comput. Chem., 2006, 27, 1787-1799. G. Henkelman and H. Jónsson, J. Chem. Phys., 2000, 113, 9978-9985. G. Henkelman, B. P. Uberuaga and H. Jónsson, J. Chem. Phys., 2000, 113, 9901-9904. J. K. Nørskov, T. Bligaard, A. Logadottir, J. R. Kitchin, J. G. Chen, S. Pandelov and U. Stimming, J. Electrochem. Soc., 2005, 152, J23. J. Greeley, T. F. Jaramillo, J. Bonde, I. Chorkendorff and J. K. Nørskov, Nat. Mater., 2006, 5, 909-913. K. Ogawa, K. Nakanishi and H. Namba, Surf. Sci., 2004, 566-568, 406-413. E. Rokuta, Y. Hasegawa, A. Itoh, K. Yamashita, T. Tanaka, S. Otani and C. Oshima, Surf. Sci., 1999, 427, 97-101. C. Ling, L. Shi, Y. Ouyang and J. Wang, Chem. Mater., 2016, 28, 9026-9032. C. Ling, L. Shi, Y. Ouyang, Q. Chen and J. Wang, Adv. Sci., 2016, 3, 1600180.


16

Nanowires Catalyze the Hydrogen Evolution ... - ACS Publications

Recommend Documents