Charge Mediated Semiconducting-to-Metallic Phase Transition in

May 21, 2015 - School of Chemical Engineering, University of Adelaide, Adelaide, SA ...... Citation data is made available by participants in Crossref...
0 downloads 3 Views 1MB Size
Article pubs.acs.org/JPCC

Charge Mediated Semiconducting-to-Metallic Phase Transition in Molybdenum Disulfide Monolayer and Hydrogen Evolution Reaction in New 1T′ Phase Guoping Gao,† Yan Jiao,‡ Fengxian Ma,† Yalong Jiao,† Eric Waclawik,† and Aijun Du*,† †

School of Chemistry, Physics and Mechanical Engineering Faculty, Queensland University of Technology, Garden Point Campus, Brisbane, QLD 4001, Australia ‡ School of Chemical Engineering, University of Adelaide, Adelaide, SA 5005, Australia S Supporting Information *

ABSTRACT: The phase transition of single layer molybdenum disulfide (MoS2) from semiconducting 2H to metallic 1T and then to 1T′ phases, and the effect of the phase transition on hydrogen evolution reaction (HER) are investigated within this work by density functional theory. Experimentally, 2H-MoS2 has been widely used as an excellent electrode for HER and can get charged easily. Here we find that the negative charge has a significant impact on the structural phase transition in a MoS2 monolayer. The thermodynamic stability of 1T-MoS2 increases with the negative charge state, comparing with the 2H-MoS2 structure before phase transition and the kinetic energy barrier for a phase transition from 2H to 1T decreases from 1.59 to 0.27 eV when 4e− are injected per MoS2 unit. Additionally, 1T phase is found to transform into the distorted structure (1T′ phase) spontaneously. On their activity toward hydrogen evolution reaction, 1T′-MoS2 structure shows comparable hydrogen evolution reaction activity to the 2H-MoS2 structure. If the charge transfer kinetics is taken into account, the catalytic activity of 1T′-MoS2 is superior to that of 2H-MoS2. Our finding provides a possible novel method for phase transition of MoS2 and enriches understanding of the catalytic properties of MoS2 for HER.



identified as the active catalytic center for HER,30,34,35 single layer 2H-MoS2 with largely exposed sulfur edge sites is considered to be a promising candidate for HER.13,30,36,37 However, one of the limitations of 2H-MoS2 as fuel cell cathode material is its poor charge transfer properties.29,37,38 Two methods have been developed to increase the cathodic current: growing 2H-MoS2 on conductive supports, and transforming 2H-MoS2 to conductive 1T-MoS2. Experiments have demonstrated MoS2 grown on conductive supports,38−40 such as graphene,40 have showed markedly faster HER kinetics. However, conducting 1T-MoS2 nanosheets have also exhibited superior electrocatalytic performance for HER.18,41,42 This has been attributed to its metallic properties, where the elimination of Schottky barriers improved the charge transfer.21,41 Phase transition from 2H-MoS2 to 1T-MoS2, or stabilization of 1T-MoS2, has been realized by both chemical method6,43 and physical method.44 Li+ intercalation is the most studied chemical method to promote the phase transition.18,45−47 However, a previous study on the Li-intercalated 1T-MoS2 showed that only a quarter of the edge sites are active for HER.16 Hence, to achieve a higher HER efficiency, a new

INTRODUCTION MoS 2 is one of the most studied transition metal dichalcogenides (TMD)1,2 and is considered to be a promising electrocatalyst for the hydrogen evolution reaction (HER). Several prototypes of MoS2 have been reported, including 2H, 1T, and 1T′.3−8 Among them, 2H9,10 and 1T3,11,12 have attracted the most extensive attention from both theory calculations and experiments in recent years.13−15 Due to the localization behavior of the d-band of transition metal,16,17 the electronic properties of the two phases are dramatically different:8 the 2H-MoS2 phase is semiconductive, while the 1T-MoS2 phase is metallic.18 Generally, the 1T phase is unstable under ambient condition and transforms back to the thermodynamically stable 2H phase17 or 1T′ (distorted 1T) phase.8 Previous studies showed that 1T phase could be stabilized by electron donation via ion intercalation,14,19,20 and the 1T and 2H phases can coexist in the same TMD nanosheets.11,21−23 However, the mechanism for the phase transition is not elucidated yet. Due to the weak van der Waals interaction between layers, MoS2 can be exfoliated into monolayer structure, which exhibits unique and fascinating electronic properties2,24−29 that could be used in a variety of potential applications, including HER,30 energy storage,31 high-performance electronics,32 and sensitive photodetectors.33 Since sulfur edge sites of MoS2 were © 2015 American Chemical Society

Received: May 15, 2015 Revised: May 21, 2015 Published: May 21, 2015 13124

DOI: 10.1021/acs.jpcc.5b04658 J. Phys. Chem. C 2015, 119, 13124−13128

Article

The Journal of Physical Chemistry C

the 1T-MoS2 by 0.80 eV per three-atom-MoS2 unit. If the Mo atoms of 1T-MoS2 distort along the y direction slightly, it is 1T′-MoS2 form (see Figure 1c,f), and the band gap of 1T′MoS2 is 0.006 eV (see Figure 1i). The energy of 1T′ phase is 0.15 eV lower than 1T phase per MoS2 unit. In the following, we investigate the phase transition between 2H and 1T MoS2 and then to 1T′ phase. The phase transition involves one of the S atoms moving from one pyramidal position to the other pyramidal position of the unit cell (see Figure S1). Figure 2a shows the pathway of the transition of 2H to 1T under charge states from 0e− to 4e− per MoS2. The relative energy of 1T-MoS2 gradually decreases from positive (0.80 and 0.31 eV in neutral state and 1e−, respectively) to negative (−0.05, −0.10, and −0.30 eV in charge state of 2e−, 3e−, and 4e−, respectively). This indicates that charge injection can stabilize the 1T-MoS2; in the case where more than 1e− is injected, 1T-MoS2 is more stable than 2H-MoS2. The phase transition energy barrier also decreases from 1.59 eV in neutral to 1.11, 0.71, 0.54, and 0.27 eV in charge states of 1e−, 2e−, 3e−, and 4e−, respectively. The general trend of MoS2 phase transition is that negative charge inducing 2H-MoS2 transform into 1T-MoS2 easily. In charge state of 4e−, not only is the 1TMoS2 form more stable than 2H-MoS2 by 0.30 eV, but also the barrier of 0.27 eV is easily overcome. Intercalation lithium (Li) or sodium (Na) on MoS2, which are studied extensively in recent experiments14,31,42,45 and theoretical calculations,56−58 introduce excessive electrons into MoS2 and promotes the phase transition. We intercalate a layer of lithium (Li) and sodium (Na) atoms on the top of twodimensional MoS2 surface to compare the ion’s effect on the transition of 2H-MoS2 to 1T-MoS2 with charge effect discussed above. The relative energies of Li and Na intercalated 1T-MoS2 are 0.16 and 0.53 eV, respectively. Meanwhile, the barrier for the Li and Na intercalated MoS2 phase transition from 2H to 1T is only 1.00 and 1.25 eV, respectively. By comparing these results with that of phase transition in the different states, ion intercalation could enhance the stability of 1T-MoS2 and reduce the barrier of the transition for 2H-MoS2 to 1T-MoS2 to some extent. However, charge injection is a superior method to promote the phase transition. To explore the mechanism for the phase transition, we calculate the different charge densities of negatively charged, Li intercalated, and Na intercalated MoS2 (see Figure S2). In all three cases, the charge transfer direction is similar, which indicates that the mechanism of ion intercalation and negative charge injection are the same. Charge mainly accumulates on the S atom induced in the transition, while it is depleted in the area among S and Mo. These all indicate that extra electrons weaken the Mo−S bonds. Meanwhile, there is a small amount of electrons to be transferred into the nonbonding d-orbitals of Mo. According to ligand field theory and refs 16 and 59, the nonbonding d-orbitals of 2H-MoS2 (group 6 TMD) split into three d-orbitals: dz2 (filled), dx2−y2,xy (empty), and dxz,yz (empty). The extra electron entrancing into nonbonding d-orbitals of MoS2 makes it an isoelectronic of group 7 TMD, whose nonbonding d-orbitals prefer to split into two d-orbitals: dyz, xz, xy (partially filled) and dz2,x2−y2 (empty) and exhibit a metallic character. The increase of d-electron count destabilizes the 2H-MoS2. Therefore, the extra electron, from charge injection or ion intercalation, promotes the phase transition by both weakening the Mo−S bond and decreasing the relative energy of 1T-MoS2.

approach, which can not only transform 2H-MoS2 to 1T-MoS2 easily but also maintain the number of active sites per unit MoS2, should be explored. Converting single layer MoS2 from 2H to 1T phase by physical methods, including strain effect and charge injection, are possible ways to obtain metallic 1T-MoS2 without sacrificing the number of active sites. For example, Voiry et al. demonstrated that 2H-WS2 could transform to 1TWS2 with improved catalytic activity by strain effect.48 Another physical factor, charge injection, which could be found in the phase transition of MoS2 in Li ion battery,49 remains unexplored. Therefore, in-depth theoretical understandings of the HER activity on 1T-MoS2 or 1T′-MoS2 are highly desired, in addition to the phase transition process. Experimentally, 2H-MoS2 has been widely used as an excellent electrode for HER and can get charged easily. Within this study, we investigate the effect of charge injection on the conversion of 2H-MoS2 to 1T-MoS2 by using the nudged elastic band (NEB) method50,51 as implemented in the Vienna ab initio simulation package.52,53 Then HER properties on distorted 1T′-MoS2 are studied by density functional theory (DFT) calculations based on Dmol3 software.54,55 More details on computational methods can be seen in the Supporting Information. Our investigations predict that an easy physical method, i.e., charge injection, can transform 2H-MoS2 to 1T′MoS2 with high density of active edges toward HER. These active sites in 1T′-MoS2 engender superior electrocatalytic HER activity compared to the well-developed 2H-MoS2.



RESULTS AND DISCUSSION MoS2 possesses a variety of polytypic structures, including 2H, 1T, and 1T′. The geometry structures and their electronic band structures are shown in Figure 1. The space-group of 2H-MoS2 is P6̅m2, where the Mo site is in a trigonal prism coordination (Figure 1a,d). The space-group of 1T-MoS2 is P3̅m1, with the Mo in octahedral coordination (Figure 1b,e). Figure 1g,h shows that the electronic properties of 2H-MoS2 and 1T-MoS2 are quite different: the band gap of 2H-MoS2 is 1.74 eV, while 1TMoS2 is metallic. The optimized 2H-MoS2 is more stable than

Figure 1. Top view and side view of 2H, 1T, and 1T′ 4 × 4 × 1 MoS2 and their band structures. Cyan, Mo; yellow, S. 13125

DOI: 10.1021/acs.jpcc.5b04658 J. Phys. Chem. C 2015, 119, 13124−13128

Article

The Journal of Physical Chemistry C

Figure 2. (a) Minimum energy pathways of the phase transition from 2H to 1T in different charge states or with ion intercalated. All the energies are relative to the total energy of 2H-MoS2 in corresponding charge states or with corresponding ion intercalated. (b) Spontaneous atom relaxation from 1T to new 1T′ phase.

Furthermore, if the angle of 1T-MoS2’s lattice vectors a and b are changed from 60° to 90° and supercell is enlarged from 1 × 1 × 1 cell to 1 × 2 × 1, the Mo atoms of 1T-MoS2 relaxes along the y direction, and the lattice transforms into the distorted 1T′-MoS2 form after cell relaxation with no barrier (see Figure 2b). The relaxation energy from 1T to 1T′ phase is −0.15 eV per MoS2 unit. Therefore, 1T phase would transform into the distorted structure (1T′ phase) spontaneously. Finally, we calculate the catalytic properties of 1T′-MoS2 on HER by Dmol3. The overall HER pathway includes two steps: Volmer step (bonding of hydrogen on the catalyst (H*) from initial state (H+ + e− + *)) and Heyrovsky60 or Tafel61 step (releasing of molecular hydrogen (1/2H2 + *), where * donates the catalyst). Experiments62 have proven that Volmer step is the rate-determining step of HER on MoS2. Therefore, the Gibbs free-energy of the adsorption atomic hydrogen (ΔG0H) is a key quantity to describe the HER activity of the catalyst. The ΔG0H of an ideal catalyst for HER should be close to 0 (ΔG0H → 0).63 In addition, previous DFT studies show that the activities sites of MoS2 mainly are Mo-edge with 50% S coverage30,37,64,65 rather than the basal plane.66,67 As the Mo of 1T′-MoS2 forms a 1D zigzag chain, there are two different Mo-edges determined by the directions. To simulate the two different Mo-edges, we build a model of 7Mo × 4Mo slab, and both the (101̅0) and (1̅010) Mo-edges are covered by 50% S (see Figure S3). The ΔG0H of the two 1T′-MoS2 Mo-edges was calculated under two coverages, 0.25 and 0.50 ML, as shown in Figure 3. The binding energy of hydrogen adsorbed on (10̅ 10) Mo-edge at 0.25 ML is too strong with a free energy of −0.421 eV, while the binding energy of hydrogen adsorbed on (1̅010) Mo-edge at 0.50 ML has a free energy of 0.187 eV. On the (101̅0) Moedge, the free energy of hydrogen is close to zero at both 0.25 ML (0.062 eV) and 0.50 ML (0.059 eV). The free energy of hydrogen on (1010̅ ) Mo-edge of 2H-MoS2 at 0.50 ML is 0.015 eV. These results indicate that the catalytic activity on (101̅0) Mo-edge of 1T′-MoS2 is comparable to the catalytic activity of 2H-MoS2. However, the poor charge transfer in semiconducting 2H-MoS2 significantly limited the kinetics of HER as suggested by experiments.39 If taking the charge transfer into

Figure 3. Free energy diagram of HER processing on (101̅0) (left panel) and (1̅010) (right panel) of 1T′-MoS2 under different coverage at equilibrium potential. Cyan, Mo; yellow, S; white, H.

account, the catalytic performance of 1T′-MoS2 exhibits superior catalytic activity, largely attributed to its conductivity.



CONCLUSIONS

In conclusion, we report that injected negative charge can significantly enhance the stability of 1T-MoS2 and promote the phase transition from semiconducting 2H-MoS2 to metallic 1TMoS2, which then further distorts into the 1T′-MoS2 phase. Remarkably, 1T′-MoS2 demonstrates comparable catalytic activity to the 2H-MoS2 phase, based on free energy calculations. Taking the charge transfer into account, the metallic 1T′-MoS2 would show superior catalytic activity compared to semiconducting 2H-MoS2. Our results suggest a practical method for the transition of 2H-MoS2 to 1T-MoS2 and enrich the information on catalytic properties of 1T′-MoS2 for HER application. 13126

DOI: 10.1021/acs.jpcc.5b04658 J. Phys. Chem. C 2015, 119, 13124−13128

Article

The Journal of Physical Chemistry C



(14) Wang, X.; Shen, X.; Wang, Z.; Yu, R.; Chen, L. Atomic-Scale Clarification of Structural Transition of MoS2 upon Sodium Intercalation. ACS Nano 2014, 11394−11400. (15) Zheng, Y.; Jiao, Y.; Jaroniec, M.; Qiao, S. Z. Advancing the Electrochemistry of the Hydrogen-Evolution Reaction through Combining Experiment and Theory. Angew. Chem., Int. Ed. 2015, 54, 52−65. (16) Chhowalla, M.; Shin, H. S.; Eda, G.; Li, L. J.; Loh, K. P.; Zhang, H. The Chemistry of Two-Dimensional Layered Transition Metal Dichalcogenide Nanosheets. Nat. Chem. 2013, 5, 263−275. (17) Bissessur, R.; Kanatzidis, M. G.; Schindler, J.; Kannewurf, C. Encapsulation of Polymers into MoS2 and Metal to Insulator Transition in Metastable MoS2. J. Chem. Soc., Chem. Commun. 1993, 1582−1585. (18) 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. (19) Kan, M.; Wang, J. Y.; Li, X. W.; Zhang, S. H.; Li, Y. W.; Kawazoe, Y.; Sun, Q.; Jena, P. Structures and Phase Transition of a MoS2 Monolayer. J. Phys. Chem. C 2014, 118, 1515−1522. (20) Yan, S.; Qiao, W.; He, X.; Guo, X.; Xi, L.; Zhong, W.; Du, Y. Enhancement of Magnetism by Structural Phase Transition in MoS2. Appl. Phys. Lett. 2015, 106, 012408. (21) Kappera, R.; Voiry, D.; Yalcin, S. E.; Branch, B.; Gupta, G.; Mohite, A. D.; Chhowalla, M. Phase-Engineered Low-Resistance Contacts for Ultrathin MoS2 Transistors. Nat. Mater. 2014, 13, 1128− 1134. (22) Eda, G.; Fujita, T.; Yamaguchi, H.; Voiry, D.; Chen, M.; Chhowalla, M. Coherent Atomic and Electronic Heterostructures of Single-Layer MoS2. ACS Nano 2012, 6, 7311−7317. (23) Yamaguchi, H.; et al. Spatially Resolved Photoexcited ChargeCarrier Dynamics in Phase-Engineered Monolayer MoS2. ACS Nano 2014, 9, 840−849. (24) Radisavljevic, B.; Radenovic, A.; Brivio, J.; Giacometti, V.; Kis, A. Single-Layer MoS2 Transistors. Nat. Nanotechnol. 2011, 6, 147−150. (25) Lee, C.; Yan, H.; Brus, L. E.; Heinz, T. F.; Hone, J.; Ryu, S. Anomalous Lattice Vibrations of Single-and Few-Layer MoS2. ACS Nano 2010, 4, 2695−2700. (26) Zong, X.; Yan, H.; Wu, G.; Ma, G.; Wen, F.; Wang, L.; Li, C. Enhancement of Photocatalytic H2 Evolution on CdS by Loading MoS2 as Cocatalyst under Visible Light Irradiation. J. Am. Chem. Soc. 2008, 130, 7176−7177. (27) Lopez-Sanchez, O.; Lembke, D.; Kayci, M.; Radenovic, A.; Kis, A. Ultrasensitive Photodetectors Based on Monolayer MoS2. Nat. Nanotechnol. 2013, 8, 497−501. (28) Yu, J. H.; Lee, H. R.; Hong, S. S.; Kong, D.; Lee, H. W.; Wang, H.; Xiong, F.; Wang, S.; Cui, Y. Vertical Heterostructure of TwoDimensional MoS2 and WSe2 with Vertically Aligned Layers. Nano Lett. 2015, 15, 1031−1035. (29) Wang, H.; Yuan, H.; Sae Hong, S.; Li, Y.; Cui, Y. Physical and Chemical Tuning of Two-Dimensional Transition Metal Dichalcogenides. Chem. Soc. Rev. 2015, 44, 2664−2680. (30) 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. (31) Hwang, H.; Kim, H.; Cho, J. MoS2 Nanoplates Consisting of Disordered Graphene-like Layers for High Rate Lithium Battery Anode Materials. Nano Lett. 2011, 11, 4826−4830. (32) Radisavljevic, B.; Whitwick, M. B.; Kis, A. Integrated Circuits and Logic Operations Based on Single-Layer MoS2. ACS Nano 2011, 5, 9934−9938. (33) Yin, Z.; Li, H.; Li, H.; Jiang, L.; Shi, Y.; Sun, Y.; Lu, G.; Zhang, Q.; Chen, X.; Zhang, H. Single-Layer MoS2 Phototransistors. ACS Nano 2011, 6, 74−80. (34) Byskov, L.; Nørskov, J.; Clausen, B.; Topsøe, H. Edge Termination of MoS2 and CoMoS Catalyst Particles. Catal. Lett. 2000, 64, 95−99.

ASSOCIATED CONTENT

S Supporting Information *

Computational details, schematic diagram of 2H-MoS2 to 1TMoS2 phase transition, different charge densities of negatively charged, Li intercalated, and Na intercalated MoS2, relaxation energy, geometric structures of 1T′-MoS2 stripes for HER, and differential hydrogen adsorption energy on 1T′-MoS2. The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.5b04658.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We acknowledge generous grants of high-performance computer time from computing facility at Queensland University of Technology and Australian National Facility. A.D. greatly appreciates the Australian Research Council QEII Fellowship (DP110101239) and financial support of the Australian Research Council under Discovery Project (DP130102420).



REFERENCES

(1) Feldman, Y.; Wasserman, E.; Srolovitz, D. J.; Tenne, R. HighRate, Gas-Phase Growth of MoS2 Nested Inorganic Fullerenes and Nanotubes. Science 1995, 267, 222−225. (2) Splendiani, A.; Sun, L.; Zhang, Y.; Li, T.; Kim, J.; Chim, C.; Galli, G.; Wang, F. Emerging Photoluminescence in Monolayer MoS2. Nano Lett. 2010, 10, 1271−1275. (3) Wypych, F.; Schö llhorn, R. 1T-MoS2, a New Metallic Modification of Molybdenum Disulfide. J. Chem. Soc., Chem. Commun. 1992, 1386−1388. (4) Wang, H.; Skeldon, P.; Thompson, G.; Wood, G. Synthesis and Characterization of Molybdenum Disulphide Formed from Ammonium Tetrathiomolybdate. J. Mater. Sci. 1997, 32, 497−502. (5) Wilson, J.; Yoffe, A. The Transition Metal Dichalcogenides Discussion and Interpretation of The Observed Optical, Electrical and Structural Properties. Adv. Phys. 1969, 18, 193−335. (6) Schonfeld, B.; Huang, J.; Moss, S. Anisotropic Mean-Square Displacements (MSD) in Single-Crystals of 2H-and 3R-MoS2. Acta Crystallogr., Sect. B: Struct. Sci. 1983, 39, 404−407. (7) Song, I.; Park, C.; Choi, H. C. Synthesis and Properties of Molybdenum Disulphide: from Bulk to Atomic Layers. RSC Adv. 2015, 5, 7495−7514. (8) Qian, X.; Liu, J.; Fu, L.; Li, J. Quantum Spin Hall Effect in TwoDimensional Transition Metal Dichalcogenides. Science 2014, 346, 1344−1347. (9) Beal, A.; Hughes, H. Kramers-Kronig Analysis of the Reflectivity Spectra of 2H-MoS2, 2H-MoSe2 and 2H-MoTe2. J. Phys. C: Solid State Phys. 1979, 12, 881−890. (10) El-Mahalawy, S.; Evans, B. The Thermal Expansion of 2HMoS2, 2H-MoSe2 and 2H-WSe2 between 20 and 800 DegreesC. J. Appl. Crystallogr. 1976, 9, 403−406. (11) Eda, G.; Yamaguchi, H.; Voiry, D.; Fujita, T.; Chen, M.; Chhowalla, M. Photoluminescence from Chemically Exfoliated MoS2. Nano Lett. 2011, 11, 5111−5116. (12) Wypych, F.; Weber, T.; Prins, R. Scanning Tunneling Microscopic Investigation of 1T-MoS2. Chem. Mater. 1998, 10, 723− 727. (13) 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. 13127

DOI: 10.1021/acs.jpcc.5b04658 J. Phys. Chem. C 2015, 119, 13124−13128

Article

The Journal of Physical Chemistry C

(55) Delley, B. An All-Electron Numerical Method for Solving the Local Density Functional for Polyatomic Molecules. J. Chem. Phys. 1990, 92, 508−517. (56) Chen, X.; He, J.; Srivastava, D.; Li, J. Electrochemical Cycling Reversibility of LiMoS2 Using First-Principles Calculations. Appl. Phys. Lett. 2012, 100, 263901. (57) Chen, X.; Chen, Z.; Li, J. Critical Electronic Structures Controlling Phase Transitions Induced by Lithium Ion Intercalation in Molybdenum Disulphide. Chin. Sci. Bull. 2013, 58, 1632−1641. (58) Enyashin, A. N.; Seifert, G. Density-Functional Ftudy of LixMoS2 Intercalates (0 × 1). Comput. Theor. Chem. 2012, 999, 13−20. (59) Kertesz, M.; Hoffmann, R. Octahedral vs. Trigonal-Prismatic Coordination and Clustering in Transition-Metal Dichalcogenides. J. Am. Chem. Soc. 1984, 106, 3453−3460. (60) Skulason, E.; Karlberg, G. S.; Rossmeisl, J.; Bligaard, T.; Greeley, J.; Jonsson, H.; Norskov, J. K. Density Functional Theory Calculations for the Hydrogen Evolution Reaction in an Electrochemical Double Layer on the Pt(111) Electrode. Phys. Chem. Chem. Phys. 2007, 9, 3241−3250. (61) Skúlason, E.; Tripkovic, V.; Björketun, M. E.; Gudmundsdottir, S.; Karlberg, G.; Rossmeisl, J.; Bligaard, T.; Jónsson, H.; Nørskov, J. K. Modeling the Electrochemical Hydrogen Oxidation and Evolution Reactions on the Basis of Density Functional Theory Calculations. J. Phys. Chem. C 2010, 114, 18182−18197. (62) Gao, M.-R.; Liang, J.-X.; Zheng, Y.-R.; Xu, Y.-F.; Jiang, J.; Gao, Q.; Li, J.; Yu, S.-H., An Efficient Molybdenum Disulfide/Cobalt Diselenide Hybrid Catalyst for Electrochemical Hydrogen Generation. Nat. Commun. 2015, 6. (63) 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. (64) Tsai, C.; Chan, K.; Abild-Pedersen, F.; Norskov, 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, 13156−13164. (65) Huang, Y.; Nielsen, R. J.; Goddard Iii, W. A.; Soriaga, M. P. The Reaction Mechanism with Free Energy Barriers for Electrochemical Dihydrogen Evolution on MoS2. J. Am. Chem. Soc. 2015, 137, 6692− 6698. (66) Wang, H.; Zhang, Q.; Yao, H.; Liang, Z.; Lee, H.-W.; Hsu, P.-C.; Zheng, G.; Cui, Y. High Electrochemical Selectivity of Edge versus Terrace Sites in Two-Dimensional Layered MoS2 Materials. Nano Lett. 2014, 14, 7138−7144. (67) Pandey, M.; Vojvodic, A.; Thygesen, K. S.; Jacobsen, K. W. Two-Dimensional Metal Dichalcogenides and Oxides for Hydrogen Evolution: A Computational Screening Approach. J. Phys. Chem. Lett. 2015, 6, 1577−1585.

(35) Bollinger, M. V.; Lauritsen, J. V.; Jacobsen, K. W.; Nørskov, J. K.; Helveg, S.; Besenbacher, F. One-Dimensional Metallic Edge States in MoS2. Phys. Rev. Lett. 2001, 87, 196803. (36) 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. (37) Tsai, C.; Abild-Pedersen, F.; Nørskov, J. K. Tuning the MoS2 Edge-Site Activity for Hydrogen Evolution via Support Interactions. Nano Lett. 2014, 14, 1381−1387. (38) Chen, W.; Santos, E. J. G.; Zhu, W.; Kaxiras, E.; Zhang, Z. Tuning the Electronic and Chemical Properties of Monolayer MoS2 Adsorbed on Transition Metal Substrates. Nano Lett. 2013, 13, 509− 514. (39) 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. (40) Chang, K.; Chen, W. l-Cysteine-Assisted Synthesis of Layered MoS2/Graphene Composites with Excellent Electrochemical Performances for Lithium Ion Batteries. ACS Nano 2011, 5, 4720−4728. (41) 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. (42) Wang, H.; Lu, Z.; Kong, D.; Sun, J.; Hymel, T. M.; Cui, Y. Electrochemical Tuning of MoS2 Nanoparticles on Three-Dimensional Substrate for Efficient Hydrogen Evolution. ACS Nano 2014, 8, 4940− 4947. (43) Lin, Y.-C.; Dumcenco, D. O.; Huang, Y.-S.; Suenaga, K. Atomic Mechanism of the Semiconducting-to-Metallic Phase Transition in Single-Layered MoS2. Nat. Nanotechnol. 2014, 9, 391−396. (44) Enyashin, A. N.; Yadgarov, L.; Houben, L.; Popov, I.; Weidenbach, M.; Tenne, R.; Bar-Sadan, M.; Seifert, G. New Route for Stabilization of 1T-WS2 and MoS2 Phases. J. Phys. Chem. C 2011, 115, 24586−24591. (45) Wang, H.; et al. Electrochemical Tuning of Vertically Aligned MoS2 Nanofilms and its Application in Improving Hydrogen Evolution Reaction. Proc. Natl. Acad. Sci. U.S.A. 2013, 110, 19701− 19706. (46) Chia, X.; Ambrosi, A.; Sedmidubský, D.; Sofer, Z.; Pumera, M. Precise Tuning of the Charge Transfer Kinetics and Catalytic Properties of MoS2 Materials via Electrochemical Methods. Chem. Eur. J. 2014, 20, 17426−17432. (47) Chia, X.; Ambrosi, A.; Sofer, Z.; Luxa, J.; Pumera, M. Catalytic and Charge Transfer Properties of Transition Metal Dichalcogenides Arising from Electrochemical Pretreatment. ACS Nano 2015, 9, 5164− 5179. (48) Voiry, D.; Yamaguchi, H.; Li, J.; Silva, R.; Alves, D. C.; Fujita, T.; Chen, M.; Asefa, T.; Shenoy, V. B.; Eda, G. Enhanced Catalytic Activity in Strained Chemically Exfoliated WS2 Nanosheets for Hydrogen Evolution. Nat. Mater. 2013, 12, 850−855. (49) Acerce, M.; Voiry, D.; Chhowalla, M. Metallic 1T phase MoS2 Nanosheets as Supercapacitor Electrode Materials. Nat. Nanotechnol. 2015, 10, 313−318. (50) 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−9904. (51) Henkelman, G.; Jónsson, H. Improved Tangent Estimate in the Nudged Elastic Band Method for Finding Minimum Energy Paths and Saddle Points. J. Chem. Phys. 2000, 113, 9978−9985. (52) Kresse, G.; Furthmüller, J. Efficiency of ab-initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15−50. (53) 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−11186. (54) Delley, B. From Molecules to Solids with the DMol3 Approach. J. Chem. Phys. 2000, 113, 7756−7764. 13128

DOI: 10.1021/acs.jpcc.5b04658 J. Phys. Chem. C 2015, 119, 13124−13128