The Role of Intrinsic Defects in Electrocatalytic ... - ACS Publications

Jan 4, 2017 - ... of Technical Physics, Chinese Academy of Science, 200083 Shanghai, China ... *E-mail: [email protected]. Phone: 86-0571-86875622., ...
1 downloads 0 Views 994KB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Article

The Role of Intrinsic Defects in Electrocatalytic Activity of Monolayer VS Basal Planes for the Hydrogen Evolution Reaction 2

Yang Zhang, Xiaoshuang Chen, Yan Huang, Chong Zhang, Feng Li, and Haibo Shu J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b11987 • Publication Date (Web): 04 Jan 2017 Downloaded from http://pubs.acs.org on January 4, 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.

The Journal of Physical Chemistry C 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

The Journal of Physical Chemistry

The Role of Intrinsic Defects in Electrocatalytic Activity of Monolayer VS2 Basal Planes for the Hydrogen Evolution Reaction Yang Zhang, †, ‡ Xiaoshuang Chen ‡,*, Yan Huang,‡ Chong Zhang,§ Feng Li,§ Haibo Shu,§,* †

National Synchrotron Radiation Laboratory, University of Science and Technology of China, 230029 Hefei, China ‡ National Laboratory for Infrared Physics, Shanghai Institute of Technical Physics, Chinese Academy of Science, 200083 Shanghai, China § College of Optical and Electronic Technology, China Jiliang University, 310018 Hangzhou, China ABSTRACT: Two-dimensional VS2 nanomaterials have emerged as highly efficient and inexpensive electrocatalysts for the hydrogen evolution reaction (HER), and the further improvement of their HER performance depends on the understanding of the catalytic mechanism and activity in various pristine and defective structures. Here structural stability, electronic properties, and HER activity of monolayer VS2 nanosheets with various intrinsic point defects are studied by using first-principles calculations. Compared to the most studied 2H-phase MoS2 basal plane, both 2H- and 1T-phase VS2 basal planes exhibit superior catalytic activity due to their metallic properties. With the introduction of intrinsic point defects onto VS2 basal planes, we find that there are four types of stable defects in 2H phase (i.e., Sad, Svac, Vad, and VS) and three types of stable defects in 1T phase (i.e., Sad, Svac, and Vad), respectively. Moreover, the formation of Svac, Vad, and VS structures in 2H phase and Vad in 1T phase can enhance the HER activity of basal planes, which implies that the synthesis of VS2 nanosheets at the V-rich condition facilitates to achieve high HER performance. The HER activity of pristine and defective VS2 structures can be well understood by a Fermi-abundance model that is also suitable to describe a broad class of electrocatalytic HER systems. This work provides a deep insight into the HER activity of single-layer VS2 and the guidance for synthesizing highly active electrocatalysts in transition-metal dichalcogenides.

1 Environment ACS Paragon Plus

The Journal of Physical Chemistry

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

1. INTRODUCTION Electrochemical water splitting through the hydrogen evolution reaction (HER) is a sustainable method to generate molecular hydrogen (H2),1-3 which has been proposed as one of the most promising clean energy solutions. Employing appropriate electrocatalysts is critical for the HER because they can accelerate the conversion from proton-electron pairs to H2 at low overpotentials (η).4 Platinum (Pt)based materials are still regarded as highly efficient catalysts due to their outstanding HER performance with high exchange current density and small Tafel slope.5 However, the high cost and scarcity of Ptbased electrocatalysts have limited their widespread use.6,7 Thus, the development of cost-effective and efficient catalysts with earth-abundant elements is strongly desirable. Recently, two-dimensional (2D) transition metal dichalcogenides (TMDs), such as MoS2 and WS2,8-12 have emerged as an alternative to Pt-based catalysts for the HER. The HER active sites have been identified from the exposed metallic edges of MoS2 and WS2 nanosheets but their less conducting basal planes are catalytically inert,13-15 which becomes a main limitation for the further improvement of HER performance. In order to explore efficient methods to achieve highly active TMDs-based catalysts, considerable efforts have been made during the past few years. Two feasible schemes have been proposed: (i) increasing the concentration of exposed edges and enhancing the conductivity of basal plane in typical semiconducting TMD materials (e.g., MoS2),16-20 (ii) finding metallic TMD materials that facilitates the HER on their basal planes.21 As a typical layered TMD material, vanadium disulfide (VS2) is an analogous of MoS2 in structure,22 in which sandwich-like S-V-S layers are interacted by van der Waals forces. Unlike semiconducting MoS2, the single-layer VS2 in 2H phase exhibits metallic characteristic2325

which makes it potentially suitable for the HER electrocatalysts.26 Motivated by unique electronic

properties of VS2, a variety of methods have been developed to fabricate high-quality VS2 monolayer and nanosheets.26-28 The synthesized VS2 nanomaterials were found to have the 2H phase in monolayer structures and the 1T phase in their bulk-like structures,29 respectively. From these experimental and theoretical studies,19,20 the phase transition from 2H to 1T can leads to the significant improvement for the HER activity of MoS2 nanostructures. Despite the large potential as the HER electrocatalysts, there are few studies about the electrocatalytic activity of 2D-VS2 with both two phases. On the other hand, engineering intrinsic defects is an important way to tune electronic properties of 2D-TMD materials, consequently leading to the change in their HER activity. For example, sulphur (S) vacancies are a universal intrinsic defect in 2D-MS2 (M = Mo, W, and Sn) nanomaterials, and the inert basal plane of these materials can be activated by introducing the S vacancies.30-32 In addition to the vacancy defects, other types of intrinsic point defects, such as antisites, holes, substitution atoms, and adadtoms, have also been identified in experiment.33,34 However, how the intrinsic defects affect the catalytic activity of 2D-VS2 remains unclear. In this work, we present a systematic investigation on the 2 Environment ACS Paragon Plus

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

The Journal of Physical Chemistry

geometry and stability of various potential intrinsic defects in monolayer VS2 basal plane with both 2H and 1T phases, and explore their HER activity by using density-functional theory (DFT) calculations. Our calculated results show that there are four types of stable defects in 2H phase (i.e., Sad, Svac, Vad, and VS) and three types of stable defects in 1T phase (i.e., Sad, Svac, and Vad) in allowed chemical potential ranges, respectively. We find that not all stable intrinsic defects contribute to the improvement of HER activity in VS2 basal planes, and highly active VS2 catalysts are easily achieved at the V-rich condition. Based on the analysis of electronic structures, we propose a Fermi-abundance model to understand the defect dependence of HER activity in monolayer VS2. 2. COMPUTATIONAL DETAILS The single-layer VS2 in 2H and 1T phases are shown in Figure 1a. The optimized lattice constants are a=b= 3.166 Å for 2H-VS2 and a=b= 3.174 Å for 1T-VS2, which is in good agreement with previous theoretical studies.25,29 Based on the optimized single-layer VS2 sheets, the defect models are created by introducing intrinsic point defects onto the 2D-VS2 supercell which has the in-plane lattice constants of 16.452 Å×15.830 Å in 2H phase and 16.491 Å×15.870 Å in 1T phase, respectively. Here we consider ten possible point defects in 2H-VS2, including of monosulfur vacancy (Svac), disulfur vacancy (S2vac), S adatom (Sad), S substitution at the Mo site (SV), disulfur substitution at the Mo site (S2V), V vacancy (Vvac), V adatom (Vad), one S atom replaced by one V atom (VS), two S atoms replaced by one V atom (VS2), and two S atoms replaced by two V atoms (V2S2), and six point defects in 1T-VS2, including of Svac, Sad, SV, Vvac, Vad, and VS, respectively. For the Sad and Vad structures, the optimal adsorption site is obtained by the energy comparison of various possible configurations. The optimized atomic models of various defective structures have been presented in Figures S1 and S2 of Supporting Information (SI). To evaluate the stability of various defective structures, their formation energies have been calculated by the following equation, E f = ED − ET + ∆ni ∑ µi i

(1)

where ED and ET are total energies of monolayer VS2 with and without intrinsic defects, respectively. µi is the chemical potential of atomic species i (i = V and S), and ∆ni is the difference of the number of V and S atoms between perfect and defective structures. To maintain thermodynamic equilibrium (µVS2 = µV +2µS) of VS2, the allowable value of µS is µS(bulk) –∆Hf/2 < µS < µS(bulk), where the upper (lower) limit corresponds to S-rich (V-rich) condition and ∆Hf is the heat of formation. Here ∆Hf is defined as ∆Hf = EVS2 – EV – ES, where EVS2, EV, and ES are the energies of VS2 monolayer, V atom in bulk, S atom in bulk, respectively. The computed ∆Hf is 1.30 eV and 1.29 eV for 2H- and 1T-phase VS2, respectively. The catalyst-supported HER is a multistep electrochemical process. The first step is the hydrogen adsorption onto a VS2 catalyst (i.e., Volmer reaction), which is described by H+ + e- + *  H*, where 3 Environment ACS Paragon Plus

The Journal of Physical Chemistry

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 prefix of * denotes the adsorption site of the catalyst. The subsequent step is the release of H2 molecules from the catalyst following either the Heyrovsky (H* + H+ + e-  H2 + *) or Tafer (2H*  H2 + 2*) reaction. To insight into the HER pathway and activity, the change of Gibbs free energies for the H adsorption (∆GH) on various VS2 catalysts have been calculated in the acid solution using the standard hydrogen electrode (CHE) model as follows,35 ∆GH = ∆EH + ∆EZPE – T∆S

(2)

where ∆EH is the adsorption energy for adding one H atom onto VS2 catalysts and is defined as ∆EH = EVS2+nH – EVS2+(n-1)H – 1/2H2, ∆EZPE and ∆S are zero-point energies (ZPE) and entropy differences between the adsorbed system and gas phase, respectively and T is temperature which refers to the room temperature (T = 298.15 K). EZPE is defined as EZPE = 1/ 2∑ hωi , where ω is the vibration frequency i=N

and N is the number of H independent freedom degree in adsorbed system or H2 gas molecule. The ∆S term is defined as ∆S = 1/2SH2 – SH*, where SH* is the vibration entropy of adsorbed H and defined as

S H * = k B ∑ ln(1 − e− hωi / kBT ) , where kB and T are Boltzmann constant and temperature respectively,36,37 i=N

and SH2 is the vibration entropy of H2 gas molecule obtained from the standard thermodynamic database.38 The detailed data of various energies and entropies have been listed in Table S1 of SI. All DFT calculations are performed using the projector augmented wave (PAW) method39 as implemented in the Vienna ab initio Simulation Package (VASP).38 The electronic exchange-correlation energy is treated by generalized-gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE).41,42 The kinetic energy cutoff for the plane-wave expansion is set to 400 eV. The k-point sampling in the Brillouin zone is implemented by the Monkhorst-Pack scheme with the grids of 16×16×1 and 4×4×1for VS2 primitive cells and supercells, respectively. For the geometry optimization, the convergence criteria of electronic energy and forces acting on each atom are 10-3 eV and 10-2 eV/Å, respectively. In the calculation of electronic structures, the DFT-GGA method can lead to the underestimation of band gap, thus the GGA+U method have been considered as a comparison. We have investigated the Hubbard parameter Ueff = U – J from 0 eV to 4 eV to compute electronic structure of pristine and defective VS2 structures. The climbing-image nudged elastic band (cNEB) method43 is employed to search transition states for the transition of VS2 monolayer from 2H to 1T phase. In the calculation, the number of image is set to 6 and the spring constant is -5. Simulated scan tunneling microscopy (STM) images are obtained using the Tersoff-Hamann approximation based on calculated charge densities, and the STM tip height is set to ~ 1.8 Å in the simulation.

4 Environment ACS Paragon Plus

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

The Journal of Physical Chemistry

Figure 1 (a) Atomic structures and barrier (in eV) of phase transition from 2H to 1T phase in monolayer VS2. (b), (c) spin-polarized band structures of 2H- and 1T-phase VS2 monolayers. The solid and dash lines denote spin-up and spindown bands, respectively.

3. RESULTS AND DISCUSSION We first investigate structural and electronic properties of prefer single-layer VS2, which can provide a reference for the study of defective structures. Figure 1a shows atomic structures and energetics of 2Hand 1T-phase VS2 monolayers. The 2H-VS2 monolayer indicates relatively higher stability than the 1T structure, which agrees well with previous studies.29 However, there is a very small energy difference (~0.04 eV) between 2H and 1T structures. Moreover, the calculated 2H-to-1T phase-transition barrier is 0.65 eV that is far lower than that of monolayer MoS2 (~1.80 eV).44 The result suggests the possibility of phase coexistence during the synthesis of single-layer VS2 using regular experimental methods, such as the chemical intercalation and chemical vapor deposition (CVD). The electronic band structures show that both 2H and 1T structures are metallic with ferromagnetic (FM) coupling (see Figures 1b and 1c), which is in good agreement with previous reports.23-25, 45 Therefore, the phase transition does not change the metallic characteristic of VS2. The calculated density of states (DOS) indicates that the magnetism of monolayer VS2 systems mainly originates from spin polarization of V-3d states (see Figures S3 and S4), resulting in that the magnetic moment of V atom is 0.5 µB in 2H phase and 4.5µB in 1T phase, respectively. Owing to the difference of FM coupling between 2H and 1T phases, there is a different distribution for electronic states around the Fermi level. The distribution of occupied states near Fermi level in 1T phase is larger than that of 2H phase. Such a difference may induce different HER activities in 2H- and 1T-phase VS2 basal planes, which will be discussed blow. It needs to be noted that the DFTGGA method underestimates the band gap. In order to confirm the metallic characteristic of 2H- and 5 Environment ACS Paragon Plus

The Journal of Physical Chemistry

1T-phase VS2 monolayers, their DOS are calculated using the GGA+U method. It can be found that both 2H- and 1T-phase VS2 monolayers are still metallic with the change of Ueff value from 0 eV to 4 eV (see Figure S5). Furthermore, we also examine electronic structure of defective VS2 structures (Svac), and calculated DOS shows that the Svac structure in both 2H and 1T phases can maintain its metallic characteristic with different Ueff values (see Figure S6). The above results suggest that metallic pristine and defective VS2 structures calculated by the DFT-GGA method are valid. To explore the defect effect on the HER activity of monolayer VS2 basal planes, it needs to firstly understand the stability of various intrinsic point defects on VS2 basal planes. The formation energies (Ef) of sixteen intrinsic point defects in 2H-VS2 as a function of S chemical potential difference (∆µS, ∆µS = µS–µS(bulk)) are plotted in Figure 2b. It can be found that Svac is the most stable point defect in the most of allowed ∆µS range, and Sad becomes the most stable defect under the S-rich regime, and Vad and VS indicate higher stability under the V-rich growth condition. The similar result has also been found in other transition-metal sulfides, such as MoS2 and SnS.33,46 For other intrinsic point defects, such as SV, S2V, S2vac,VS2, and V2S2, they present relatively larger formation energies. According to the defect concentration C ∝ exp(−E f / kBT ) , a large Ef value means that the defect is difficult to be formed on VS2 basal planes during the experimental synthesis. As shown in Figure 2b, the stability of intrinsic defects on 1T basal plane is very similar to the case of 2H phase: Svac is the most stable defect for most of allowed range of ∆µS, and Sad and Vad become the most stable intrinsic defect under the S-rich and Vrich conditions, respectively. 8

6

4

S2V SV

V2S2

(b)

VS2

6

SV

4

S2vac

Vvac

Ef (eV)

(a)

Ef (eV)

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

Vvac

VS

2

2

VS 0

0

Svac

Vad -1.2 Mo-rich

Sad -2

-0.8

-0.4

∆µS (eV)

0.0 S-rich

Svac

Vad -1.2 Mo-rich

-0.8

-0.4

∆µS (eV)

Sad 0.0 S-rich

Figure 2 Formation energies of various intrinsic point defects on (a) 2H and (b) 1T VS2 basal planes as a function of S chemical potential difference ∆µS.

Figure 3 shows atomic structures and simulated STM images of these stable intrinsic point defects on VS2 basal planes. For the Svac structures, introducing an S vacancy (Svac) onto both 2H- and 1T-phase VS2 basal planes brings three nearest-neighboring V atoms with the dangling bonds (see Figure 3a,e). The coulomb repulsion effect induces the increase of atomic distance among these V atoms from ~3.15 6 Environment ACS Paragon Plus

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

The Journal of Physical Chemistry

Å to ~3.86 Å. The simulated STM images display that there are bright spots at the sites of three nearestneighboring V atoms approaching the S vacancy in 2H phase (see Figure 3a), but there is a dark region around the S vacancy in 1T phase (see Figure 3e). The difference of STM images means that the S vacancy should have different effects on the HER activity of VS2 basal plane in two phases. For the Sad structures, the S adatom prefers to occupy at the top of V atom and bonds with one V atom and two S atoms in 2H phase (see Figure 3b), while the S adatom prefers to situate at the top of S atom, forming SS dimer in 1T phase (see Figure 3f). The simulated STM images indicate a highly bright spot at the site of S adatom in both 2H and 1T phases due to the localization of unpaired electrons at the defect, implying that the S adatom has a high activity to trap H atoms during the HER. In the Vad structures, the V adatom prefers to adsorb at the hollow site regardless of the crystal phase. It is seen from Figures 3c,g that simulated STM images present a triangle-like bright spot in both 2H and 1T phases, which arises from unpaired electrons of V adatom delocalized to its adjacent V atoms. For the VS structure (see Figure 3d), the V substitution atom causes a local structural distortion in 2H basal plane, leading to the elongation of V-V bonds and the constriction of V-S bonds. Hence, the simulated STM image displays a high bright spot at the site of V substitution atom and the dark region at its nearest-neighboring V atoms. The above results reflect that these intrinsic defects lead to a significant change in electronic states of VS2 basal planes at the vicinity of Fermi level, which may have a significant effect on their HER activity.

Figure 3 Side and top views of atomic structures and simulated STM images of stable intrinsic point defects on (a)-(d) 2H and (e)-(g) 1T VS2 basal planes. The simulated STM images are obtained at a voltage of -1.0 V using the TersoffHamann approximation.

7 Environment ACS Paragon Plus

The Journal of Physical Chemistry

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

Now we turn to investigate the HER activity of perfect and defective VS2 basal planes. The change of Gibbs free energy for hydrogen adsorption (∆GH) is an important parameter to measure the intrinsic catalytic activity of a catalyst. An ideal catalyst should be able to facilitate HER with zero for ∆GH (∆GH = 0). Generally, the ideal case cannot be achieved due to the binding interaction of H with the catalyst. Therefore, a catalyst with lower ∆GH would have better HER performance. Based on this principle, the free-energy diagram through either Volmer-Heyrovsky or Volmer-Tafel pathway has been firstly calculated to reveal the HER mechanism on VS2 basal planes. As is shown in Figure 4a (VolmerHeyrovsky pathway), the free energy change (∆GH) on both 2H- and 1T-phase pristine basal planes is the first Volmer step, which serves as the rate-limiting step for the overall HER at equilibrium potential. In contrast, as shown in Figure 4b (Volmer-Tafel pathway), the ∆GH of the Tafel step on two pristine basal planes is larger than that of the Volmer step, and the overall ∆GH of the Volmer-Tafel pathway is also greater than that of the Volmer-Heyrovsky pathway. Therefore, the Volmer-Heyrovsky pathway should be energetically favorable for the HER of VS2-based catalysts and the ∆GH of the Volmer step is thus used to evaluate the HER activity of VS2 catalysts. The result is similar to the case of heteroatomdoped graphene HER electrocatalysts.47 It is seen from Figure 4a that the calculated ∆GH of pristine VS2 basal plane in 2H and 1T phases is 0.26 eV and -0.16 eV respectively, which is far lower than that of monolayer MoS2 basal plane (∆GH =1.89 ~1.92 eV)25,48 and closes to that of Pt (∆GH =-0.09 eV).49 The lower ∆GH of VS2 basal planes originate from their metallic electronic structures which facilitate to the charge transfer from the catalytic surfaces to the H adatoms. The result suggests a large potential of monolayer VS2 as the highly efficient HER electrocatalysts.

Figure 4 Schematic illustration and the corresponding free-energy diagram of the HER following (a) the VolmerHeyrovsky pathway and (b) the Volmer-Tafel pathway on 2H- and 1T-phase pristine VS2 basal planes.

8 Environment ACS Paragon Plus

Page 9 of 16

Based on the Volmer-Heyrovsky pathway, the free-energy diagrams of various defective VS2 structures in 2H and 1T phases have been compared. For the comparison, the data of pristine VS2 and Pt catalysts have been provided here. We find that the introduction of most of intrinsic point defects more or less increases the surface activity of VS2 basal planes except for 1T-Vad and 1T-Svac structures (see Figure 5a), resulting in calculated ∆GH more negative. For example, the calculated ∆GH of Vad, VS, Svac, and Sad structures in 2H phase is 0.17, 0.46, 0.55, and 0.73 eV lower than that of pristine 2H basal plane, respectively. Similarly, the calculated ∆GH of 1T-Sad structure is 0.41 eV lower than that of the pristine structure (see Figure 5b). Overall, the defective structures formed under V-rich condition, such as Vad and VS in 2H phase and Vad in 1T phase, have the relatively smaller absolute values of ∆GH (< 0.26 eV). Therefore, the V-rich growth condition should facilitate the synthesized VS2 catalysts with high HER activity. We find that the formation of Svac defect in 2H and 1T basal phases leads to different change trends for the ∆GH relative to the pristine ones. As shown in Figures 5a and 5b, the ∆GH shows the reduction of ~0.55 eV in 2H basal plane but has the increase of ~0.42 eV in 1T basal plane with the formation of Svac defect. In addition, introducing Sad defect onto VS2 basal planes can induce the significant reduction of ∆GH (> 0.4 eV) in both 2H and 1T phase, which is related to the high activity of S adatom to trap H atom. These results agree well with the prediction from the simulated STM images. 0.4

(b)

Pristine 0.2 +

0.0

H +e

-

Vad

-0.2

Pt VS

-0.4

Svac

-0.6

Sad

0.4

Svac

0.2 1/2 H2

+

∆GH (eV)

(a)

∆GH (eV)

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

The Journal of Physical Chemistry

0.0

H +e

-

Pt

-0.2

Vad

1/2 H2

Pristine

-0.4 Sad -0.6

Reaction coordinate

Reaction coordinate

Figure 5 Calculated free-energy diagrams for the hydrogen evolution on pristine and defective VS2 basal planes in (2) 2H phase and (b) 1T phase. For the comparison, ∆GH of Pt has also been included here.

In order to understand the inherent mechanism of HER activity, electronic structures of pristine and defective VS2 structures in 2H and 1T phases have been calculated. In the HER process, the adsorption and desorption of H are strongly related to its binding strength on a catalyst that substantially depends on electronic structures of the catalyst. The previous studies have demonstrated that the d-band center (εd)50,51 is an efficient descriptor for the activity of transition-metal-based materials and it is defined as, +∞

εd

∫ = ∫

ED ( E ) dE

−∞ +∞ −∞

D ( E ) dE

9 Environment ACS Paragon Plus

(3)

The Journal of Physical Chemistry

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

where D(E) is d states of transition metal in catalysts. Figure 6a indicates ∆GH of various VS2 catalysts as a function of εd. However, we find that ∆GH does not display a regular trend with εd. In other words, the d-band center is not suitable to describe the HER activity of VS2 catalysts. Considering that closer the electronic state is to the Fermi level, the greater its contribution to bonding interaction. Therefore, the reactivity of a catalyst should be determined by both its density of states (D(E)) and a weight factor (w(E)) that quantifies the contribution of every electronic state to the surface reaction.52 For the HER, the activity of an electrocatalyst relies on the abundance of occupied states near the Fermi level because the reaction refers to the charge transfer from the catalyst to the H adatoms. The magnitude of occupied states near the Fermi level (FA) is obtained by the sum of the weight contribution of occupied states as follows (see Figure 6b), FA = ∫

EF

−∞

(4)

D( E ) w( E )dE

where EF is the energy of Fermi level. From previous studies, the derivative of the Fermi-Dirac function, 52 − fT′(E − EF ) , has been validated as a good weight function. The Fermi-Dirac function is defined as,

fT ( E − E F ) =

1 exp[( E − EF ) / kT ] + 1

(5)

Figure 6 (a) ∆GH of pristine and defective VS2 structures as a function of d-band center (εd). (b) The abundance of occupied states (FA) is defined as a weight sum of occupied states of density, E D( E ) w( E )dE , where w(E) is the weight



F

−∞

function that is applied by using the derivative of the Fermi-Dirac function, − fT′(E − EF ) . (c) ∆GH of pristine and defective VS2 structures as a function of DF using kT = 0.4 eV. The overlapped shadow region indicates the optimal range for highly active HER catalysts.

10 Environment ACS Paragon Plus

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

The Journal of Physical Chemistry

Using the − fT′(E − EF ) as the weight function can ensure that reaches a maximum at EF and descends to zero with the increase of |E–EF|, as shown in Figure 6b. It can be found that the spreading of the w(E) is sensitive to the magnitude of kT that reaches the optimal value at 0.4 eV. Based on above analysis, we propose a Fermi-abundance (DF) model to describe the HER activity and it is defined as follows, EF

DF

∫ = ∫

ED( E ) fT′ ( E − EF )dE

−∞ EF −∞

(6)

D( E ) fT′ ( E − EF )dE

Considering that electronic states of perfect and defective VS2 near the Fermi level are mainly from V atoms, here D(E) refers to the projected DOS of V atoms in various VS2 structures (see Figures S7 and S8 in SI). Figure 6c shows that ∆GH of pristine and defective VS2 structures as a function of DF. The ∆GH values of various VS2 catalysts indicate an approximately linear trend in relation to their DF values, suggesting that DF can become a good descriptor for the HER activity. A smaller DF value represents the distribution of few occupied states at the vicinity of Fermi level, resulting in that the catalyst has a weak surface activity to capture H atoms (e.g., 1T-Svac). In contrast, the catalyst with a larger DF value will have a strong surface activity to trap H (e.g., 1T-Sad). The calculated DOS indicates that the 1T-Sad structure has greater occupied states than that of the 1T-Svac structure near the Fermi level (see Figure S8). However, a good catalyst should balance the H adsorption/desorption strength based on the Sabatier principle. Therefore, the DF values of the VS2 catalysts located in the range from -0.6 to -0.5 eV should have higher HER activity, such as 2H-Vad and 2H-VS structures. The result reflects that the HER activity of these VS2 catalysts is actually decided by the occupied states near the Fermi level. To examine the rationality of Fermi-abundance model, we also calculated the DF value of Pt(111) surface (see Figure S9). The DF of Pt(111) surface is -0.54, thus Pt(111) surface has a good HER activity. Such a result agrees well with theoretical and experimental reports on the HER activity of Pt catalysts.35,49 Therefore, the Fermi-abundance model should be of generality and be also used to describe the HER activity of other type of electrocatalysts.

4. CONCLUSIONS In summary, we have performed systematic DFT calculations to insight into the role of intrinsic point defects in HER activity of monolayer VS2 catalysts. The calculation of formation energies indicates that there are four types of stable intrinsic defects in 2H phase (i.e., Sad, Svac, Vad, and VS) and three types of stable intrinsic defects in 1T phase (i.e., Sad, Svac, and Vad), respectively. The calculated results of ∆GH suggest that the formation of Svac, Vad, and VS structures in 2H phase and the Vad structure in 1T phase can improve the HER activity of VS2 basal planes, which implies that the synthesis of VS2 should be maintained at the V-rich regime facilitating the catalysts with high HER performance. Based on the 11 Environment ACS Paragon Plus

The Journal of Physical Chemistry

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

analysis of electronic structures, we find that the HER activity of pristine and defective VS2 structures can be well understood by a Fermi-abundance model, which is also applicable to describe a broad class of electrocatalysts for the HER.  ASSOCIATED CONTENT  Supporting Information Atomic structures of intrinsic point defects in 2H and 1T VS2 basal planes, DFT energies and entropies for calculating ∆GH of pristine and defective VS2 structures, electronic structure of 2H and 1T VS2 monolayers, DOS of pristine VS2 monolayers and Svac defective structures calculated using GGA+U method, projected DOS of V atoms in pristine and stable defective VS2 structures, and projected DOS of Pt(111) Surface. This material is available free of charge via the internet at http://pubs.acs.org.  AUTHOR INFORMATION Corresponding Authors *E-mail: [email protected] (Haibo Shu), Phone: 86-0571-86875622; [email protected] (Xiaoshuang Chen), Phone: 86-021-35052062. Notes The authors declare no competing financial interest.  ACKNOWLEDGMENT This work was supported in part by the National Natural Science Foundation of China (Grant no. 11404309, 51402275, 11334008, and 61290301) and the Fund of Shanghai Science and Technology Foundation (Grant no. 13JC1408800). Computational resources from the Shanghai Supercomputer Center are acknowledged.

REFERENCES AND NOTES (1) Turner, J. A. Sustainable Hydrogen Production. Science 2004, 305, 972-974. (2) Merki, D.; Hu, X. Recent Developments of Molybdenum and Tungsten Sulfides as Hydrogen Evolution Catalysts. Energy Environ. Sci. 2011, 4, 3878-3888. (3) Laursen, A. B.; Kegnaes, S.; Dahl, S.; Chorkendorff, I. Molybdenum Sulfides – Efficient and Viable Materials for Electro- and Photoelectrocatalytic Hydrogen Evolution. Energy Environ. Sci. 2012, 5, 5577-5591. (4) Jiao, Y.; Zheng, Y.; Jaroniec, M.; Qiao, S. Z. Design of Electrocatalysts for Oxygen- and Hydrogen-Involving Energy Conversion Reactions. Chem. Soc. Rev. 2015, 44, 2060-2086.

12 Environment ACS Paragon Plus

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

The Journal of Physical Chemistry

(5) Esposito, D. V.; Hunt, S. T.; Stottlemyer, A. L.; Dobson, K. D.; McCandless, B. E.; Birkmire, R. W.; Chen, J. G. Low-Cost Hydrogen-Evolution Catalysts Based on Monolayer Platinum on Tungsten Monocarbide Substrates. Angew. Chem. Int. Ed. 2010, 49, 9859-9862. (6) Durst, J.; Siebel, A.; Simon, C.; Hasché, F.; Herranz, J.; Gasteiger, H. A. New Insights into the Electrochemical Hydrogen Oxidation and Evolution Reaction Mechanism. Energy Environ. Sci. 2014, 7, 2255-2260. (7) Vesborg, P. C.; Jaramillo, T. F. Addressing the Terawatt Challenge: Scalability in the Supply of Chemical Elements for Renewable Energy. RSC Adv. 2012, 2, 7933-7947. (8) Wang, H.; Lu, Z.; Xu, S.; Kong, D.; Cha, J. J.; Zheng, G.; Hsu, P.-C.; Yan, K.; Bradshaw, D.; Prinz, F. B.; et al. Electrochemical Tuning of Vertically Aligned MoS2 Nanofilms and Its Application in Improving Hydrogen Evolution Reaction. Proc. Nat. Acad. Sci. USA 2013, 110, 19701-19706. (9) Yu, Y.; Huang, S.-Y.; Li, Y.; Steinmann, S. N.; Yang, W.; Cao, L. Layer-Dependent Electrocatalysis of MoS2 for Hyrogen Evolution. Nano Lett. 2014, 14, 553-558. (10)

Voiry, D.; Yang, J.; Chhowalla, M. Recent Strategies for Improving the Catalytic Acitivity of 2D TMD

Nanosheets Toward the Hydrogen Evolution Reaction. Adv. Mater. 2016, 28, 6197-6206. (11)

Yu, S.; Kim, J.; Yoon, K. R.; Jung, J.-W.; Oh, J.; Kim. I.-D. Rational Design of Efficient Electrocatalysts

for Hydrogen Evolution Reaction: Single Layers of WS2 Nanoplates Anchored to Hollow Nitrogen-Doped Carbon Nanofibers. ACS Appl. Mater. Interfaces, 2015, 7, 28116-28121. (12)

Xu, K.; Wang, F.; Wang, Z.; Zhan, X.; Wang, Q.; Cheng, Z.; Safdar, M.; He, J. Component-Controllable

WS2(1-x)Se2x Nanotubes for Efficient Hydrogen Evolution Reaction. ACS Nano 2014, 8, 8468-8476. (13)

Jaramillo, T. F.; Jøergensen, 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. (14)

Karunadasa, H. I.; Montalvo, E.; Sun, Y.; Majda, M.; Long, J. R.; Chang, C. J. A Molecular MoS2 Edge

Site Mimic for Catalytic Hydrogen Generation. Science 2012, 335, 698-702. (15)

Voriy, D.; Salehi, M.; Silva, R.; Fujita, T.; Chen, M.; Asefa, T.; Shenoy, V. B.; Eda, G.; Chhowalla, M.

Conducting MoS2 Nanosheets as Catalyst for Hyrogen Evolution Reaction. Nano Lett. 2013, 13, 6222-6227. (16)

Kibsgaard, J.; Chen, Z.; Reinecke, B. N.; Jaramillo, T. F. Engineering the Surface Structure of MoS2 to

Preferentially Expose Active Edge Sites for Electrocatalysis. Nat. Mater. 2012, 11, 963-969. (17)

Zhao, X.; Ma, X.; Sun, J.; Li, D.; Yang, X. Enhanced Catalytic Activities of Surfactant-Assisted

Exfoliated WS2 Nanodots for Hydrogen Evolution. ACS Nano 2016, 10, 2159-2166. (18)

Deng, J.; Li, H.; Xiao, J.; Tu, Y.; Deng, D.; Yang, H.; Li, J.; Ren, P.; Bao, X. Triggering the

Electrocatalytic Hydrogen Evolution Activity of the Inert Two-Dimensional MoS2 Surface via Single-Atom Metal Doping. Energy Environ. Sci. 2015, 8, 1594-1601. (19)

Tang, Q.; Jiang, D. E. Mechanism of Hydrogen Evolution Reaction on 1T-MoS2 from First Principles.

ACS Catal. 2016, 6, 4953-4961. (20)

Lukowski, M. A.; Daniel, A. S.; Meng, F.; Forticaux, A.; Li, L.; Jin, S. Ehanced Hydrogen Evolution

Catalysis from Chemial Exfoliated Metallic MoS2 Nanosheets. J. Am. Chem. Soc. 2013, 135, 10274-10277. (21)

Pan, H. Metal Dichalcogenides Monolayers: Novel Catalysts for Electrochemical Hydrogen Production.

Sci. Rep. 2014, 4, 5348.

13 Environment ACS Paragon Plus

The Journal of Physical Chemistry

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

(22)

Page 14 of 16

Tang, C.; Zhang, H. Two-Dimensional Transition Metal Dichalcogenide Nanosheet-Based Composites.

Chem. Soc. Rev. 2015, 44, 2713-2731. (23)

Ma, Y.; Dai, Y.; Guo, M.; Niu, C.; Zhu, Y.; Huang, B. Evidence of the Exitence of Magnetism in Pristine

VX2 Monolayers (X = S, Se) and Their Strain-Induced Tunable Magnetic Properties. ACS Nano 2012, 6, 1695-1701. (24)

Jing, Y.; Zhou, Z.; Cabrera, C. R.; Chen, Z. F. Metallic VS2 Monoalyer: A Promising 2D Anode Material

for Lithium Ion Batteries. J. Phys. Chem. C 2013, 117, 25409-25413. (25)

Fan, X.; Wang, S.; An, Y.; Lau, W. Catalytic Activity of MS2 Monolayer for Electrochemical Hydrogen

Evolution. J. Phys. Chem. C 2016, 120, 1623-1632. (26)

Liang, H.; Shi, H.; Zhang, D.; Ming, F.; Wang, R.; Zhuo, J.; Wang, Z. Solution Growth of Vertical VS2

Nanoplate Arrays for Electrocatalytic Hyrogen Evolution. Chem. Mater. 2016, 28, 5587-5591. (27)

Feng, J.; Sun, X.; Wu, C.; Peng, L.; Lin, C.; Hu, S.; Yang, J.; Xie, Y. Metallic Few-Layered VS2 Ultathin

Nanosheets: High Two-Dimensional Conductivity for In-Plane Supercapacitors. J. Am. Chem. Soc. 2011, 133, 17832-17838. (28)

Fang, W.; Zhao, H.; Xie, Y.; Fang, J.; Xu, J.; Chen, Z. Facile Hydrothermal Synthesis of VS2/graphene

Nanocomposites with Superior High-Rate Capability as Lithium-Ion Battery Cathodes. ACS Appl. Mater. Interfaces 2015, 7, 13044-13052. (29)

Zhang, H.; Liu, L.-M.; Lau, W.-M. Dimensional-Dependent Phase Transition and Magnetic Properties of

VS2. J. Mater. Chem. A 2013, 1, 10821-10828. (30)

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. (31)

Le, D.; Rawal, T. B.; Rahman, T. S. Single-Layer MoS2 with Sulfur Vacancies: Structure and Catalytic

Application. J. Phys. Chem. C 2014, 118, 5346-5351. (32)

Lin, L.; Miao, N.; Wen, Y.; Zhang, S.; Ghosez, P.; Sun, Z.; Allwood, D. A. Sulfur-Depleted

Monolayered Molybdenum Disulfide Nanocrystals for Superelectochemical Hydrogen Evolution Reaction. ACS Nano 2016, 10, 8929-8937. (33)

Hong, J.; Hu, Z.; Probert, M.; Li, K.; Lv, D.; Yang, X.; Gu, L.; Mao, N.; Feng, Q.; Xie, L.; et al.

Exploring Atomic Defects in Molybdenum Disulphide Monolayers. Nat. Comm. 2015, 6, 6293. (34)

Najmaei, S.; Yuan, J.; Zhang, J.; Ajayan, P.; Lou, J. Synthesis and Defect Investigation of Two-

Dimensional Molybdenum Disulfide Atomic Layers. Acc. Chem. Res. 2015, 48, 31-40. (35)

Greeley, J.; Jaramillo, T. F.; Bonde, J.; Chorkendorff, I.; Nørskov, J. K. Computational High-Throughput

Screening of Electrocatalytic Materials for Hydrogen Evolution. Nat. Mater. 2006, 5, 909-913. (36)

Reuter, K.; Scheffler, M. Composition, Structure, and Stability of RuO2(110) as a Function of Oxygen

Pressure. Phys. Rev. B 2001, 65, 035406. (37)

Shu, H. B.; Tao, X.-M.; Ding, F. What are the Active Carbon Species During Graphene Chemical Vapor

Deposition Growth? Nanoscale 2015, 7, 1627-1634. (38)

Cramer, C. J. Essentials of Computational Chemistry Theories and Models, 2nd ed.; John Wiley & Sons,

Ltd.: West Sussex, England, 2004.

14 Environment ACS Paragon Plus

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

The Journal of Physical Chemistry

(39)

Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method.

Phys. Rev. B 1999, 59, 1758. (40)

Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953-17979.

(41)

Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for ab Initio Total-Energy Calculations Using a

Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169-11185. (42)

Payne, M. C.; Teter, M. P.; Allan, D. C.; Arias, T. A.; Joannopoulos, J. D. Iterative Minimization

Technologies for Ab Initio Total-Energy Calcualtions: Molecular Dynamics and Conjugate Gradients. Rev. Mod. Phys. 1992, 64, 1045-1097. (43)

Henkelman, G.; Uberuaga, B. P.; Jónsson, H. A Climbing Image Nudged Elastic Band Method for

Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901. (44)

Shu, H.; Li, F.; Hu, C.; Liang, P.; Cao, D.; Chen, X. S. The Capacity Fading Mechanism and

Improvement of Cycling Stability in MoS2-based Anode Materials for Lithium-Ion Batteries. Nanoscale 2016, 8, 2918-2926. (45)

Wasey, A. H. M. A.; Chakrabarty, S.; Das, G. P. Quantum Size Effects in Layered VX2 (X=S, Se)

Materials: Manifestation of Metal to Semimetal or Semiconductor Transition. J. Appl. Phys. 2015, 117, 064313. (46)

Vidal, J.; Lany, S.; d´Avezac, M.; Zunger, A.; Zakutayev, A.; Francis, J.; Tate, J. Band-Structure, Optical

Properties, and Defect Physics of the Photovoltaic Semiconductor SnS. Appl. Phys. Lett. 2016, 1, 16130. (47)

Jiao, Y.; Zheng, Y.; Davey, K.; Qiao, S.-Z. Activity Origin and Catalyst Design Principles for

Electrocatalytic Hydrogen Evolution on Heteroatom-Doped Graphene. Nat. Comm. 2012, 100, 032104. (48)

Tsai, C.; Chan, K.; Pedersenb, F. A.; Nørskov, J. K. Active Edge Sites in MoSe2 and WSe2 Catalysts for

the Hydrogen Evolution Reaction: a Density Functional Study. Phys. Chem. Chem. Phys. 2014, 16, 1315613164. (49)

Ouyang, Y.; Ling, C.; Chen, Q.; Wang, Z.; Shi, L.; Wang, J. L. Activating Inert Basal Planes of MoS2 for

Hydrogen Evolution Reaction through the Formation of Different Intrinsic Defects. Chem. Mater. 2016, 28, 4390-4396. (50)

Hammer, B.; Nørskov, J. K. Electronic Factors Determining the Reactivity of Metal Surfaces. Surf. Sci.

1995, 343, 211-220. (51)

Li, F.; Shu, H.; Hu, C.; Shi, Z.; Liu, X.; Liang, P.; Chen, X. S. Atomic Mechanism of Electrocatalytically

Active Co-N Complexes in Graphene Basal Plane for Oxygen Reduction Reaction. ACS Appl. Mater. Interfaces 2015, 7, 27405-27413. (52)

Huang, B.; Xiao, L.; Lu, J.; Zhang, L. Spatially Resolved Quantification of the Surface Reactivity of

Solid Catalysts. Angew. Chem. Int. Ed. 2016, 55, 6239-6243.

15 Environment ACS Paragon Plus

The Journal of Physical Chemistry

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

TOC Graphic

16 Environment ACS Paragon Plus

Page 16 of 16