Theoretical Investigation of the Intercalation ... - ACS Publications

Jun 2, 2017 - Tsinghua-Berkeley Shenzhen Institute (TBSI), Tsinghua University, ... School at Shenzhen, Tsinghua University, Shenzhen 518055, P. R. Ch...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF ARIZONA

Article

Theoretical Investigation of the Intercalation Chemistry of Li/Na Ions in Transition Metal Dichalcogenides Shaoxun Fan, Xiaolong Zou, Hongda Du, Lin Gan, Chengjun Xu, Wei Lv, Yan-Bing He, Quanhong Yang, Feiyu Kang, and Jia Li J. Phys. Chem. C, Just Accepted Manuscript • Publication Date (Web): 02 Jun 2017 Downloaded from http://pubs.acs.org on June 2, 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 21

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

Theoretical Investigation of the Intercalation Chemistry of Li/Na Ions in Transition Metal Dichalcogenides Shaoxun Fan,† Xiaolong Zou,∗,‡ Hongda Du,† Lin Gan,† Chengjun Xu,† Wei Lv,¶ Yan-Bing He,¶ Quan-Hong Yang,¶,‡ Feiyu Kang,†,¶,‡ and Jia Li∗,† Guangdong Provincial Key Laboratory of Thermal Management Engineering and Materials, Graduate School at Shenzhen, Tsinghua University, Shenzhen 518055, P. R. China, Tsinghua-Berkeley Shenzhen Institute (TBSI), Tsinghua University, Shenzhen 518055, P. R. China, and Engineering Laboratory for Functionalized Carbon Materials and Shenzhen Key Laboratory for Graphene-based Materials, Graduate School at Shenzhen, Tsinghua University, Shenzhen 518055, P. R. China E-mail: [email protected]; [email protected]



To whom correspondence should be addressed Guangdong Provincial Key Laboratory of Thermal Management Engineering and Materials, Graduate School at Shenzhen, Tsinghua University, Shenzhen 518055, P. R. China ‡ Tsinghua-Berkeley Shenzhen Institute (TBSI), Tsinghua University, Shenzhen 518055, P. R. China ¶ Engineering Laboratory for Functionalized Carbon Materials and Shenzhen Key Laboratory for Graphene-based Materials, Graduate School at Shenzhen, Tsinghua University, Shenzhen 518055, P. R. China †

1

ACS Paragon Plus Environment

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

Abstract Among various two-dimensional compounds, transition metal dichalcogenides (TMDs or MX2 ) are a group of materials attracting growing research interest for potential applications as battery electrodes. Here we systematically investigate the electrochemical performance of a series of MX2 (M = Mo, W, Nb, Ta; X = S, Se) upon Li/Na intercalation through first-principles calculations. MoX2 and WX2 were found to have lower voltages while those of NbX2 and TaX2 were higher than 1.5 V. By applying the rigidband model, we found that the energy gained for electrons to transfer from Li/Na to MX2 could serve as a descriptor for characterizing voltages of MX2 .The linear relation between the descriptor and voltages is useful for screening candidates for electrodes with desired voltage. Migration barriers for Li/Na ions were approximately 0.3 eV in MoX2 /WX2 and 0.5 eV in NbX2 /TaX2 . The low barriers suggest a reasonable rate performance when these TMDs are used as electrodes. By stacking different MX2 with distinct properties, TMDs heterostructures could be adopted to provide tunable electrochemical properties, including voltage, capacity and electronic conductivity while keeping barriers for Li/Na ions little changed. Thus this strategy offers another degree of freedom for rational design of layered electrode materials.

2

ACS Paragon Plus Environment

Page 2 of 21

Page 3 of 21

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

Introduction Since the discovery of graphene in 2004, 1 much effort has been devoted to searching for other two-dimensional (2D) materials that exhibit excellent physical and chemical properties. 2 In particular, 2D transition metal dichalcogenides (TMDs) have received considerable attention because of their interesting electronic, optical, and chemical properties. 3 2D TMDs are a group of layered compounds with a general chemical formula of MX2 , where M and X denote transition metal (e.g., Mo, W, etc.) and chalcogen (S, Se, Te) elements, respectively. Each layer of MX2 is three atoms thick, with one transition metal layer sandwiched between two chalcogen layers. These three layers are stacked together by weak van der Waals (vdW) interactions. 3,4 Due to the unique physical properties of 2D TMDs, many potential applications were explored, including logic integrated circuits, 5,6 photodetectors, 7,8 and light-emitting devices. 9,10 Besides their possible applications in electronics and optoelectronics, the layered structure of 2D TMDs also renders them promising applications for energy storage, such as in lithium-ion batteries (LIBs) and sodium-ion batteries (SIBs) 11–15 . Over the past decade, LIBs have become the most commonly used energy solution for portable electronic devices because of their high energy density. 16 However, the high cost of Li limits the use of LIBs in small- and mid-size energy storage applications. Compared to LIBs, SIBs are particularly suitable for grid-scale energy storage because of the natural abundance and easy accessibility of Na 17,18 . Nevertheless, LIBs and SIBs still suffer from low capacity and/or poor cycling stability of electrode materials. Further improvement of performance of LIBs and SIBs calls for search for other advanced electrode materials, including both cathode and anode materials. In this regard, 2D TMDs with tunable composition for both metals and chalcogenides have attracted considerable attention as promising electrode materials. 19 For the most extensively studied TMD, MoS2 , the electrochemical performance of Li/Na intercalation has been investigated both experimentally and theoretically. 20–23 Many MoS2 nanomaterials have been synthesized that show high capacity along with good cycling stability as LIB or SIB anode materials. 24,25 Regarding other TMD compounds, MoSe2 was 3

ACS Paragon Plus Environment

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

found to be a good anode candidate for both LIBs and SIBs, with high reversible specific capacity and excellent cycling stability. 26,27 Both 2H- and 3R-NbS2 were synthesized as electrode materials for LIBs with specific capacities of about 169.5 and 169.0 mAh/g at 0.05 C rate, respectively. 28 Meanwhile, 2H-NbS2 also showed a high specific capacity of 143.6 mAh/g in SIBs. 29 These examples highlight the potential of finding ideal electrode materials in the vast choices of layered TMDs. However, there is no systematic understanding of the electrochemical performance of TMDs to guide experimental efforts to find novel electrode materials. In this work, employing first-principles calculations, we systematically investigated the electrochemical performance of a series of bulk MX2 (M = Mo, W, Nb, Ta; X = S, Se) TMDs upon Li/Na intercalation. MoX2 and WX2 are found to be good anode materials for LIBs or SIBs due to the low voltages and good electronic conductivities after intercalation. Moreover, Li and Na ions have the lowest energy barriers for the migration in WX2 , which is beneficial for the rate performance of batteries. As for NbX2 and TaX2 , the voltages are larger than 1.5 V when the Li/Na are intercalated. And the barriers for the migration of Li/Na ions are higher in NbX2 and TaX2 . Importantly, a linear dependence between the energy gained for electrons to transfer into empty states of MX2 (Wf illing ) and the voltage of the Li/Na intercalation has been established, making Wf illing a simple descriptor for the voltages of MX2 used for LIBs or SIBs. TMDs heterostructure composed of MX2 with distinct properties, are found to be able to provide tunable electrochemical performance in electronic conductivity, voltage, and capacity, offering another strategy in rational design of electrode materials.

Computational Methods All calculations were performed using density functional theory (DFT) implemented in the Vienna ab initio Simulation Package. 30 The projected augmented wave potential 31 and gen-

4

ACS Paragon Plus Environment

Page 4 of 21

Page 5 of 21

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

eralized gradient approximation of the Perdew–Burke–Ernzerhof functional 32 were used to describe the electron–ion interaction and exchange-correlation energy, respectively. The cutoff energy was set to 650 eV. We used 2 × 2 × 1 and 2 × 2 × 2 supercells for 2H-, and 1T-phases bulk TMDs, respectively, so that both systems consisted of two layers of TMDs including eight transition metal and sixteen chalcogen atoms. The Γ-centered 7 × 7 × 5 grid and 17 × 17 × 11 k-point samplings were used for geometry optimization and electronic structure calculations, respectively. All structures were relaxed without any symmetry restriction until the total forces acted on every atom were less than 0.01 eV/˚ A. To accurately describe vdW interactions, we used the DFT-D3 empirical correction method, 33 which has been proven reliable for describing long-range vdW interactions. The climbing-image nudged elastic band method 34 was used to determine the energy barrier and minimum energy path for the migration of Li/Na ions in the interlayer of TMDs.

Results and Discussion In TMDs, each of the transition metal atoms is coordinated to six chalcogen atoms with a triangular prismatic or octahedral structure. Based on the coordination structure of the transition metal and the stacking sequence of MX2 layers, TMDs can be classified into several phases including 1T and 2H, 3 which are the most common phases. Because other phases such as 3R possess similar intercalation environments for Li/Na ions to that of the 1T or 2H phase, only these two phases are considered in our work. Figure 1 presents the structures of the 1T and 2H phases of TMDs. In the 1T phase, transition metal atoms are surrounded by chalcogen atoms forming triangular prisms and MX2 layers are stacked in AA sequence. In the 2H phase, transition metal atoms are octahedrally coordinated, with the MX2 layers stacked in AB sequence. Both 2H- and 1T-MX2 contain two kinds of intercalation sites for Li/Na ions: octahedral sites and tetrahedral sites, as shown in Figure 1. The octahedral site is more favorable

5

ACS Paragon Plus Environment

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

Figure 1: Schematics of the structures of (a) 1T- and (b) 2H-MX2 ; (c) octahedral and (d) tetrahedral intercalation sites for 2H-MX2 . Transition metal atoms, chalcogen atoms and intercalated Li/Na atoms are represented by purple, yellow and green spheres, respectively. for lithiation than the tetrahedral site in MoS2 . 22 This is the case for all TMDs studied here upon Li or Na ion intercalation, most likely because of the higher coordination number and approximately three times larger space for inserted ions in octahedral sites than in tetrahedral ones. Therefore, we consider that all Li/Na ions intercalated only at octahedral sites to evaluate their stability and voltage profiles. The volume changes of TMDs during Li/Na intercalation to form the fully lithiated/sodiated state are listed in Table 1. These values indicate that Li or Na intercalation will tend to markedly expand the lattice in c direction, and the expansion caused by sodiation can be one to two times larger than that by lithiation because the ion radius of Na (1.02 ˚ A) is considerably larger than that of Li (0.76 ˚ A). Moreover, transition metal sulfides or selenides with transition metals in the same column of the periodic table have similar lattice constants before and after ion intercalation, especially in a and b directions. For pristine 1T-MX2 , the lattice constants are also similar to those of 2H-MX2 . The layered structure of the TMDs effectively mitigates the volume expansion issue commonly faced by electrode materials. 6

ACS Paragon Plus Environment

Page 6 of 21

Page 7 of 21

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

Table 1: Lattice constants and volumes of deintercalated, and Li- and Na-ion intercalated MX2 and their maximum capacities.

2H phase MoS2 MoSe2 WS2 WSe2 NbS2 NbSe2 TaS2 TaSe2 1T phase MoS2 MoSe2 WS2 WSe2 NbS2 NbSe2 TaS2 TaSe2

a,b/˚ A

MX2 c/˚ A

V/˚ A3

a,b/˚ A

LiMX2 c/˚ A

V/˚ A3

a,b/˚ A

NaMX2 c/˚ A

V/˚ A3

Capacity/ mAh g−1

6.3 6.6 6.3 6.6 6.6 6.9 6.6 6.9

12.3 13.0 12.4 13.0 11.8 12.4 12.2 12.7

427.2 487.3 432.2 488.7 451.1 509.9 464.4 522.8

6.5 6.8 6.5 6.8 6.7 7.0 6.7 7.0

13.1 14.0 13.3 14.1 12.6 13.4 12.6 13.4

478.7 552.7 480.8 552.8 486.1 564.1 486.5 562.6

6.6 6.9 6.6 6.8 6.8 7.1 6.7 7.1

14.5 15.3 14.7 15.9 14.1 14.8 14.1 14.8

543.0 626.2 545.4 629.6 558.5 639.6 554.9 638.7

166.8 105.2 107.6 78.1 170.0 106.4 108.9 78.8

6.4 6.6 6.4 6.6 6.7 6.9 6.7 6.9

12.1 13.0 12.1 13.1 11.7 12.4 12.0 12.5

422.8 485.8 431.4 489.4 449.2 509.1 467.8 520.4

6.5 7.0 6.7 7.0 6.7 7.0 6.7 7.0

12.7 12.9 12.3 12.9 12.6 13.3 12.7 13.4

470.0 548.3 479.2 547.7 491.0 566.6 490.4 565.3

6.6 7.1 6.8 7.1 6.8 7.1 6.9 7.0

14.2 14.4 13.7 14.4 14.1 14.9 14.1 15.1

541.4 625.2 545.1 621.9 566.4 648.6 573.4 646.8

166.8 105.2 107.6 78.1 170.0 106.4 108.9 78.8

To provide insight into the electrochemical properties of the Li/Na intercalation process into MX2 , eight different intercalation states, from (Li/Na)0.125 MX2 to (Li/Na)MX2 , were chosen to study the phase transition, voltages, and ionic/electronic transport in the Li/Na intercalated states. By comparing the total energies of intercalated states in different phases (Figure 2 and Figure S1), phase transformations could be identified in the Li/Na intercalated MoX2 and WX2 from 2H phase to 1T phase, which is consistent with reported results. 20,21,35–37 For NbX2 and TaX2 , no such transformation occurred, which can be understood by analyzing the energy acquired for charge transfer as detailed below. Convex hulls 38 (Figure S2 and S3, Supporting Information) were constructed by calculating formation enthalpies of intermediate states to determine energetically stable states and thus voltage profiles of various TMDs. The equation used to derive voltages is expressed

7

ACS Paragon Plus Environment

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 21

Figure 2: Total energies of 1T- and 2H-MoS2 in different (a) lithiated and (b) sodiated states. Red and black points represent 1T and 2H phases, respectively. as: 39 V =−

EAx1 M X2 − (EAx2 M X2 + (x1 − x2 )EA ) (x1 − x2 )e

(1)

where A represents Li or Na, EA and EAx M X2 represent Li/Na bulk energy per atom and total energy of Ax MX2 per unit formula, respectively. Figure 3 shows calculated average voltages as a function of Li/Na concentration in various MX2 . As seen from Figure 3, for all cases, the voltages of lithiated states are higher than those of sodiated states. Moreover, the voltages of lithiated/sodiated MoX2 and WX2 are below 1 V at first, and become higher after the phase transition, while voltages of NbX2 and TaX2 exceed 1.5 V, consistent with experimental results. 28 The relatively low voltages of MoX2 and WX2 mean that they are suitable for anode materials. Specifically, the voltage of WX2 is about 0.3 V lower than that of MoX2 , so WX2 may be capable of storing more energy as anode materials than MoX2 . If matched with a cathode that has a voltage of 4 V, WX2 can provide a volumetric energy density that is 250 Wh/L (approximately 10%) higher than that of MoX2 . The migration barrier of ions is the dominant factor that determines the rate performance of batteries. Figure S4 shows the typical migration path of an intercalated ion from one octahedral site to an adjacent one through a tetrahedrally coordinated intermediate site. The migration barriers of a Li/Na ion in different 2H-A0.125 MX2 along this path are shown in Table S1. For Li ion, it has the lowest energy barriers of around 0.25 eV in WS2 and WSe2 . Although the barriers become higher for NbX2 and TaX2 , they are below 0.6 eV in all 8

ACS Paragon Plus Environment

Page 9 of 21

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 3: Average voltages as a function of (a) Li and (b) Na concentration in Nb/TaX2 and (c) Li and (d) Na concentration in Mo/WX2 . considered TMDs, which could ensure reasonable rate performance at room temperature. 40 Regarding Na ions, its migration barriers are slightly higher in MX2 and WX2 than the barriers of Li ions while those in NbX2 and TaX2 are slightly lower compared to lithium ions. The lowest barriers of 0.27 and 0.28 eV occur in WS2 and WSe2 too. In comparison, in-plane migration barriers of Li and Na ions in graphite are about 0.48 and 0.29 eV, respectively. 41 These results suggest that WS2 and WSe2 have the highest Li-ion conductivities among considered 2D materials. Moreover, the barriers may also be influenced by the concentration of inserted ions. The migration barrier of the lithium ion in Li0.031 MoS2 is calculated to be 0.44 eV, which is close to the reported result. 23 The migration barriers of Li ions in Li0.25 MoS2 and Li0.5 MoS2 are 0.23 eV and 0.19 eV, respectively. It can be concluded that as the ion concentration increases, the migration barrier would become even lower, which is probably because the intercalation of more ions can expand the interlayer spacing more, thus allowing ions to migrate more easily. Rate performance may also be influenced by the electronic properties of materials, espe-

9

ACS Paragon Plus Environment

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

Figure 4: Density of states of (a) WS2 , (b) NbS2 , (c) Li0.125 WS2 , (d) Li0.125 NbS2 , (e) LiWS2 , and (f) LiNbS2 . The Fermi level was set to zero. LiWS2 is in 1T phase while others are in 2H phases. cially for semiconducting TMDs. Therefore, density of states (DOS) were also calculated to understand the change of electronic structures and conducting behaviors of TMDs during ion intercalation. Because TMDs with transition metals in the same column of the periodic table have similar behaviors upon intercalation, WS2 and NbS2 are chosen as examples in the following discussion. Figure 4 shows the DOS of WS2 and NbS2 in different lithiated states. Pristine WS2 has a bandgap of approximately 0.9 eV, suggesting its semiconducting behavior. Once Li ions are intercalated into the vdW gap, the system becomes n-type doped and turns into a conductor. It will maintain this conducting behavior after the phase transition until fully lithiated when it turns into semiconductor again as shown in Figure 4e. In contrast, NbS2 shows a different trend, transforming from a metal in the pristine form to a semiconductor after being fully lithiated. This phenomenon can be understood by considering the rigid-band model. The electronic structure of MX2 changes little when lithiated, especially in the case of low ion concentration. The intercalation of Li ions introduces additional electrons from Li-2s orbitals into the empty bands of the MX2 substrate. Thus, the Fermi level of the lithiated MX2 will move up. For WS2 , the Fermi level moves from

10

ACS Paragon Plus Environment

Page 10 of 21

Page 11 of 21

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 5: Voltages as a function of (a) interlayer binding energy and (b) Wf illing for different MX2 . Orange and green labels represent 2H and 1T phases. Red and Black plots represent voltages of Na and Li intercalation, respectively. The dashed lines in (b) are linear fitting of the data. the valence band maximum to the conduction band when lithiated, markedly enhancing its electronic conductivity. For NbS2 , the Fermi level moves into a gap when fully lithiated; this extreme situation is potentially harmful to electron transport. More importantly, these analyses show that the intercalation process is mainly determined by charge transfer from Li/Na to TMDs, regardless of whether they are metallic or semiconducting. To gain a deeper understanding of the intercalation, we divide it into two steps: the expansion of the interlayer spaces between TMD layers when intercalated and the interaction between intercalated Li/Na atoms and the TMDs. The space expansion between TMD layers is controlled by the interlayer interaction. As a direct reflection of the energy change of intercalation, voltage might be influenced by this interaction because a stronger interaction means more energy will be required to expand the layers, thus leading to a lower voltage. We calculated the interlayer binding energies of each MX2 and compared them with their first stage voltages. As seen in Figure 5a, the interlayer binding energies of MSe2 are larger than those of MS2 , and this may be partly responsible for the voltage trends. However, the interlayer binding energies of metallic (Ta/Nb)X2 are similar to those of (Mo/W)X2 , while their voltages are much higher. In addition, the relationships between voltages and interlayer binding energies for TMDs are not monotonic. All these suggest that the difference in interlayer binding energies of the TMDs cannot account for their change of voltage. 11

ACS Paragon Plus Environment

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 21

As shown above, the main effect of intercalation can be described by electron transfer from Li/Na to the TMDs. The voltage of intercalation should be strongly correlated with the energy change of this process. Because the rigid-band model is reliable in dilute cases, we used the state-filling model 42 to calculate the energy acquired when an extra electron transferred into an empty state above the Fermi level of pristine MX2 described by:

Wf illing =

Z E′

Eρ(E)dE

Ef

(2)

where E is the energy that referenced to the vacuum level by comparing local potentials of bulk MX2 and their corresponding slabs, ρ(E) is the DOS of pristine MX2 and E ′ is derived from

R E′ Ef

ρ(E)dE = 1, which assumes that there is one extra electron transferred into MX2

per supercell. The amount of charge transfer may not be precisely one, as previously shown for 2D blue phosphorus. 43 The exact values are calculated by Bader analysis and the results are shown in Table S2. The amount of charge transfer is above 0.8 in all cases. Actually, their influence on Wf illing is small, as shown in SI, and can be approximated as one here. The lower the Wf illing is, the larger the energy gain during intercalation. We take the lithiated state of A0.125 MX2 as the dilute limit, which keeps the electronic structure of MX2 in the lithiated state similar to that of the pristine one. Figure 5b plots the corresponding voltages of each MX2 upon lithiation as a function of calculated Wf illing . Intriguingly, the voltage displays a good linear dependence on Wf illing , which confirms that electron transfer plays an important role in the Li/Na intercalation process of 2H phase TMDs. This simple relation should be applicable to all TMDs or even other layered materials, as shown in Figure 5 for NbX2 and TaX2 in their 1T phases. Such a linear correlation provides guidance to explore the electrochemical performance of other materials. Besides, Wf illing could help us to understand the phase transition in TMDs. As seen in Figure 2, the energy of 1T-MoS2 decreases faster and becomes lower than that of 2H-MoS2 as lithium/sodium concentration increases. Based on the discussions above, this trend can be partly attributed to a lower Wf illing value of

12

ACS Paragon Plus Environment

Page 13 of 21

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

about -5.1 eV for 1T-MoS2 compared to -4.2 eV for its 2H phase, which leads to a larger energy decrease when lithiated. It should be noted that the deviation of 1T-MoX2 and WX2 from the linear dependence originates from the distortion of the intercalated compound, as shown in Figure S5. For other four 1T-MX2 , their Wf illing are higher than those of their 2H phases which rationalizes the stability of their 1T phases. TMD heterojunctions have attracted much research interest recently because various combinations of different 2D TMDs could endow them with unique properties useful in electronics or other applications. 44,45 Such heterostructures can also be used as battery electrodes, and could potentially provide another degree of freedom to tune material properties. Complementary electrochemical performance incorporating merits of different TMDs could be potentially achieved in this way. For example, MoS2 has a relatively high capacity, while WS2 has a lower voltage than MoS2 and migration barrier of Li/Na ion in it is the lowest among all TMDs studied here. Therefore, we designed a MoS2 -WS2 heterostructure to investigate whether this combination is promising for LIBs. As shown in Table 2, the lattice constants, voltage ranges and capacities of the heterostructure all lie between those of its two components. Meanwhile, the migration barrier of Li ions in 2H-(MoS2 )0.50 (WS2 )0.50 is 0.25 eV, the same as in WS2 . Besides this semiconductor–semiconductor heterostructure, we also examined metal– semiconductor (NbS2 -MoS2 ) and metal–metal (NbS2 -TaS2 ) heterostructures. There are no phase transition found in these two materials during intercalation. Similar to the case of MoS2 -WS2 , their properties like lattice constants (Table 2) and voltage profiles (Figure S6) are also between those of their components. Interestingly, the metal–semiconductor heterostructure was always a metal during the whole intercalation process (Figure S7), which indicates an overall higher electronic conductivity than those of the other TMDs we studied. Given the vast number of layered materials with distinct properties, these results suggest that hetero-stacking could serves as an effective strategy to improve the electrochemical performance of electrodes based on 2D materials.

13

ACS Paragon Plus Environment

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

Table 2: Lattice constants, voltage ranges and capacities of the heterostructure. MoS2 -WS2 NbS2 -TaS2 NbS2 -MoS2

a,b / ˚ A 6.3 6.6 6.5

c/˚ A 12.3 11.9 12.0

Voltage range / V 0.95 - 2.81 1.93 - 2.85 1.29 - 2.47

Capacity / mAh g−1 137.2 139.5 168.4

Conclusions To summarize, we have performed a series of DFT calculations to investigate the electrochemical properties of a group of TMDs intended for use as novel electrode materials. It was found that MoX2 and WX2 can be used as anode materials for LIBs and SIBs because they have relatively low voltages and high electronic conductivities during Li/Na intercalation. As for NbX2 and TaX2 , their voltages versus lithium/sodium are higher than those of MoX2 and WX2 . Our calculations indicated that the voltages during Li/Na intercalation were dominated by charge transfer from Li/Na to MX2 and the energy gained for the extra electron to transfer into an empty state above the Fermi level of pristine MX2 (Wf illing ) is a good descriptor for predicting the voltages of layered materials. Moreover, Li/Na ions have lowest migration barriers in WS2 and WSe2 , suggesting these materials should provide batteries with favorable rate performance. Based on the above assessment, heterostructures were investigated and shown to be good candidates for electrodes as it could combine merits of different TMDs. This may be a promising approach to design electrodes based on 2D materials with high electrochemical performance.

Supporting Information Available Electrochemical performance of TMDs, including migration barriers for Li/Na ions in various 2H-MX2 ; amount of charge transferred; formation enthalpy and convex hulls of MX2 upon Li/Na intercalation; total energies of 1T and 2H phase TMDs; schematic of migratin path of Li/Na ions in MX2 ; schematic structure of 1T-Li0.125 MoS2 ; voltage profiles of three kinds of heterosturctures and their components; density of states of NbS2 -MoS2 heterostructure in 14

ACS Paragon Plus Environment

Page 14 of 21

Page 15 of 21

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

different lithation states. This material is available free of charge via the Internet at http://pubs.acs.org/.

Acknowledgement This work was supported by the Ministry of Science and Technology of China (Grant No. 2014CB932400), the National Natural Science Foundation of China (Grant No. 51232005), Shenzhen Projects for Basic Research (Grant No. KQCX20140521161756227) and the National Program for Thousand Young Talents of China. Guangzhou Supercomputing Center is also acknowledged for allowing the use of computational resources including TIANHE-2.

References (1) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666–669. (2) Butler, S. Z.; Hollen, S. M.; Cao, L.; Cui, Y.; Gupta, J. A.; Guti´errez, H. R.; Heinz, T. F.; Hong, S. S.; Huang, J.; Ismach, A. F. et al. Progress, Challenges, and Opportunities in Two-Dimensional Materials Beyond Graphene. ACS Nano 2013, 7, 2898–2926. (3) Wang, Q. H.; Kalantar-Zadeh, K.; Kis, A.; Coleman, J. N.; Strano, M. S. Electronics and Optoelectronics of Two-Dimensional Transition Metal Dichalcogenides. Nat. Nanotechnol. 2012, 7, 699–712. (4) Lv, R.; Robinson, J. A.; Schaak, R. E.; Sun, D.; Sun, Y.; Mallouk, T. E.; Terrones, M. Transition Metal Dichalcogenides and Beyond: Synthesis, Properties, and Applications of Single- and Few-Layer Nanosheets. Acc. Chem. Res. 2015, 48, 56–64.

15

ACS Paragon Plus Environment

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

(5) Radisavljevic, B.; Radenovic, A.; Brivio, J.; Giacometti, V.; Kis, A. Single-Layer MoS2 Transistors. Nat. Nanotechnol. 2011, 6, 147–150. (6) Kim, S.; Konar, A.; Hwang, W.-S.; Lee, J. H.; Lee, J.; Yang, J.; Jung, C.; Kim, H.; Yoo, J.-B.; Choi, J.-Y. et al. High-Mobility and Low-Power Thin-Film Transistors Based on Multilayer MoS2 Crystals. Nat. Commun. 2012, 3, 1011. (7) Lopez-Sanchez, O.; Lembke, D.; Kayci, M.; Radenovic, A.; Kis, A. Ultrasensitive Photodetectors Based on Monolayer MoS2 . Nat. Nanotechnol. 2013, 8, 497–501. (8) Roy, K.; Padmanabhan, M.; Goswami, S.; Sai, T. P.; Ramalingam, G.; Raghavan, S.; Ghosh, A. Graphene-MoS2 Hybrid Structures for Multifunctional Photoresponsive Memory Devices. Nat. Nanotechnol. 2013, 8, 826–830. (9) Baugher, B. W. H.; Churchill, H. O. H.; Yang, Y.; Jarillo-Herrero, P. Optoelectronic Devices Based on Electrically Tunable p-n Diodes in a Monolayer Dichalcogenide. Nat. Nanotechnol. 2014, 9, 262–7. (10) Zhang, Y. J.; Oka, T.; Suzuki, R.; Ye, J. T.; Iwasa, Y. Electrically Switchable Chiral Light-Emitting Transistor. Science 2014, 344, 725–728. (11) Pumera, M.; Sofer, Z.; Ambrosi, A. Layered Transition Metal Dichalcogenides for Electrochemical Energy Generation and Storage. J. Mater. Chem. A 2014, 2, 8981. (12) Ko, Y. N.; Choi, S. H.; Park, S. B.; Kang, Y. C. Hierarchical MoSe2 Yolkshell Microspheres with Superior Na-ion Storage Properties. Nanoscale 2014, 6, 10511. (13) Jing, Y.; Zhou, Z.; Cabrera, C. R.; Chen, Z. Graphene, Inorganic Graphene Analogs and Their Composites for Lithium Ion Batteries. J. Mater. Chem. A 2014, 2, 12104. (14) Jing, Y.; Ortiz-Quiles, E. O.; Cabrera, C. R.; Chen, Z.; Zhou, Z. Layer-by-Layer Hybrids of MoS2 and Reduced Graphene Oxide for Lithium Ion Batteries. Electrochim. Acta 2014, 147, 392–400. 16

ACS Paragon Plus Environment

Page 16 of 21

Page 17 of 21

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

(15) Zhao, X.; Zhang, X.; Wu, D.; Zhang, H.; Ding, F.; Zhou, Z. Ab Initio Investigations on Bulk and Monolayer V2 O5 as Cathode Materials for Li-, Na-, K- and Mg-Ion Batteries. J. Mater. Chem. A 2016, 4, 16606–16611. (16) Croguennec, L.; Palacin, M. R. Recent Achievements on Inorganic Electrode Materials for Lithium-Ion Batteries. J. Am. Chem. Soc. 2015, 137, 3140–3156. (17) Kundu, D.; Talaie, E.; Duffort, V.; Nazar, L. F. The Emerging Chemistry of Sodium Ion Batteries for Electrochemical Energy Storage. Angew. Chemie Int. Ed. 2015, 54, 3431–3448. (18) Yang, J.; Zhou, X.; Wu, D.; Zhao, X.; Zhou, Z. S-Doped N-Rich Carbon Nanosheets with Expanded Interlayer Distance as Anode Materials for Sodium-Ion Batteries. Adv. Mater. 2017, 29, 1604108. (19) Morales, J.; Santos, J.; Tirado, J. L. Electrochemical Studies of Lithium and Sodium Intercalation in MoSe2 . Solid State Ionics 1996, 83, 57–64. (20) Xiong, F.; Wang, H.; Liu, X.; Sun, J.; Brongersma, M.; Pop, E.; Cui, Y. Li Intercalation in MoS2 : In Situ Observation of Its Dynamics and Tuning Optical and Electrical Properties. Nano Lett. 2015, 15, 6777–6784. (21) Gao, P.; Wang, L.; Zhang, Y.; Huang, Y.; Liu, K. Atomic-Scale Probing of the Dynamics of Sodium Transport and Intercalation-Induced Phase Transformations in MoS2 . ACS Nano 2015, 9, 11296–11301. (22) Mortazavi, M.; Wang, C.; Deng, J.; Shenoy, V. B.; Medhekar, N. V. Ab Initio Characterization of Layered MoS2 as Anode for Sodium-Ion Batteries. J. Power Sources 2014, 268, 279–286. (23) Li, Y.; Wu, D.; Zhou, Z.; Cabrera, C. R.; Chen, Z. Enhanced Li Adsorption and

17

ACS Paragon Plus Environment

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

Diffusion on MoS2 Zigzag Nanoribbons by Edge Effects: A Computational Study. J. Phys. Chem. Lett. 2012, 3, 2221–2227. (24) Sun, P.; Zhang, W.; Hu, X.; Yuan, L.; Huang, Y. Synthesis of Hierarchical MoS2 and Its Electrochemical Performance as an Anode Material for Lithium-ion Batteries. J. Mater. Chem. A 2014, 2, 3498–3504. (25) Hu, Z.; Wang, L.; Zhang, K.; Wang, J.; Cheng, F.; Tao, Z.; Chen, J. MoS2 Nanoflowers with Expanded Interlayers as High-Performance Anodes for Sodium-Ion Batteries. Angew. Chemie Int. Ed. 2014, 53, 12794–12798. (26) Wang, H.; Wang, X.; Wang, L.; Wang, J.; Jiang, D.; Li, G.; Zhang, Y.; Zhong, H.; Jiang, Y. Phase Transition Mechanism and Electrochemical Properties of Nanocrystalline MoSe2 as Anode Materials for the High Performance Lithium-Ion Battery. J. Phys. Chem. C 2015, 119, 10197–10205. (27) Yang, X.; Zhang, Z.; Fu, Y.; Li, Q. Porous Hollow Carbon Spheres Decorated with Molybdenum Diselenide Nanosheets as Anodes for Highly Reversible Lithium and Sodium Storage. Nanoscale 2015, 7, 10198–10203. (28) Liao, Y.; Park, K. S.; Singh, P.; Li, W.; Goodenough, J. B. Reinvestigation of the Electrochemical Lithium Intercalation in 2H- and 3R-NbS2 . J. Power Sources 2014, 245, 27–32. (29) Liao, Y.; Park, K. S.; Xiao, P.; Henkelman, G.; Li, W.; Goodenough, J. B. Sodium Intercalation Behavior of Layered Nax NbS2 (0 x 1). Chem. Mater. 2013, 25, 1699– 1705. (30) Kresse, G.; Furthm¨ uller, 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.

18

ACS Paragon Plus Environment

Page 18 of 21

Page 19 of 21

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

(31) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector AugmentedWave Method. Phys. Rev. B 1999, 59, 1758–1775. (32) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. (33) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (34) Henkelman, G.; Uberuaga, B. P.; Jonsson, H. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901. (35) Wang, X.; Shen, X.; Wang, Z.; Yu, R.; Chen, L. Atomic-Scale Clarification of Structural Transition of MoS2 upon Sodium Intercalation. ACS Nano 2014, 8, 11394–11400. (36) 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. (37) Cheng, Y.; Nie, A.; Zhang, Q.; Gan, L.-Y.; Shahbazian-Yassar, R.; Schwingenschlogl, U. Origin of the Phase Transition in Lithiated Molybdenum Disulfide. ACS Nano 2014, 8, 11447–11453. (38) Mayo, M.; Griffith, K. J.; Pickard, C. J.; Morris, A. J. Ab Initio Study of Phosphorus Anodes for Lithium- and Sodium-Ion Batteries. Chem. Mater. 2016, 28, 2011–2021. (39) Ong, S. P.; Chevrier, V. L.; Hautier, G.; Jain, A.; Moore, C.; Kim, S.; Ma, X.; Ceder, G. Voltage, Stability and Diffusion Barrier Differences Between Sodium-Ion and LithiumIon Intercalation Materials. Energy Environ. Sci. 2011, 4, 3680.

19

ACS Paragon Plus Environment

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

(40) Ceder, G.; Hautier, G.; Jain, A.; Ong, S. Recharging Lithium Battery Research with First-Principles Methods. MRS Bulletin 2011, 36, 185–191. (41) Wang, Z.; Ratvik, A. P.; Grande, T.; Selbach, S. M. Diffusion of Alkali Metals in the First Stage Graphite Intercalation Compounds by vdW-DFT Calculations. RSC Adv. 2015, 5, 15985–15992. (42) Liu, Y.; Wang, Y. M.; Yakobson, B. I.; Wood, B. C. Assessing Carbon-Based Anodes for Lithium-Ion Batteries: A Universal Description of Charge-Transfer Binding. Phys. Rev. Lett. 2014, 113, 028304. (43) Zhang, J.-J.; Dong, S. Prediction of Above 20 K Superconductivity of Blue Phosphorus Bilayer with Metal Intercalations. 2D Mater. 2016, 3, 035006. (44) Gong, Y.; Lin, J.; Wang, X.; Shi, G.; Lei, S.; Lin, Z.; Zou, X.; Ye, G.; Vajtai, R.; Yakobson, B. I. et al. Vertical and In-Plane Heterostructures from WS2 /MoS2 Monolayers. Nat. Mater. 2014, 13, 1135–1142. (45) Hong, X.; Kim, J.; Shi, S.-f.; Zhang, Y.; Jin, C.; Sun, Y.; Tongay, S.; Wu, J.; Zhang, Y.; Wang, F. Ultrafast Charge Transfer in Atomically Thin MoS2 /WS2 Heterostructures. Nat. Nanotechnol. 2014, 9, 682–686.

20

ACS Paragon Plus Environment

Page 20 of 21

Page 21 of 21

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

Graphical TOC Entry

21

ACS Paragon Plus Environment