Metal-Ions Intercalation Mechanism in Layered Anode From First-Principles Calculation.

Layered structure (MoS2) has the potential use as an anode in metal-ions (M-ions) batteries. Here, first-principles calculations are used to systematically investigate the diffusion mechanisms and structural changes of MoS2 as anode in lithium (Li)-, sodium (Na)-, magnesium (Mg)- and Zinc (Zn)-ions batteries. Li and Na ions are shown to be stored in the MoS2 anode material due to the strong adsorption energies (~-2.25 eV), in contrast to a relatively weak adsorption of Mg and Zn ions for the pristine MoS2. To rationalize the results, we evaluate the charge transfer from the M-ions to the MoS2 anode, and find a significant hybridization between the adsorbed atoms and S atoms in the MoS2 anode. Furthermore, the migration energy barriers of M ions are explored using first-principles with the climbing image nudged elastic band (CINEB) method, and the migration energy barrier is in the order of Zn > Mg > Li > Na ions. Our results combined with the electrochemical performance experiments show that Li- and Na-ions batteries have good cycle and rate performance due to low ions migration energy barrier and high storage capability. However, the MoS2 anode shows poor electrochemical performance in Zn- and Mg-ions batteries, especially Zn-ion batteries. Further analysis reveals that the MoS2 structure undergoes the phase transformation from 2H to 1T during the intercalation of Li and Na ions, leading to strong interaction between M ions and the anode, and thus higher electrochemical performance, which, however, is difficult to occur in Mg- and Zn-ions batteries. This work focuses on the theoretical aspects of M-ions intercalation, and our findings may stimulate the experimental work for the intercalation of multi-ions to maximize the capacity of anode in M-ions batteries.

Introduction

The rapid development of eco-friendly batteries will bring huge benefits to electrical vehicles and capacity devices (Manoj et al., 2018), and the rechargeable batteries with high energy density and long cycle life have attracted considerable attention in terms of improving the energy storage efficiency. Lithium (Li)-ions batteries represent a highly attractive and challenging alternative to rechargeable batteries. While many important achievements have been achieved for these batteries (Ju et al., 2019; Wang et al., 2020), the content of Li is limited in the Earth's crust, which will affect the extensive applications of Li-ions batteries (Schmuch et al., 2018; Sun et al., 2019), suggesting the necessity of the research of multivalent metal-ions (M-ions) batteries.

M (Li, Na, Mg, and Zn)-ions batteries have been demonstrated as very promising rechargeable batteries (Forsyth et al., 2019), and in order to enhance their capacity, many research efforts have been devoted to designing electrodes and solid electrolytes, including the nanostructured Si, transition metal oxides, and layered structure (Hu et al., 2019, 2020a). However, the ability to maintain high capacity and long cycle remains a bottleneck for M-ions batteries, especially anode materials, with a stable anode as the key component in M-ions batteries (Hu et al., 2020a). Therefore, it should be noted that finding an anode with fast ion and electron conduction is the biggest challenge in improving the performance of M-ions batteries.

Interestingly, due to their large surface-to-volume ratios, two-dimensional materials have currently become the research focus in nanostructured anode materials and have been successfully applied in M-ions batteries based on theoretical and experimental studies (Wang et al., 2018; Mohanapriya and Jha, 2019; Hu et al., 2020b). For example, layer anode materials, such as MoS2, WS2, borophene and graphene, were investigated for Li-ions batteries, and they exhibited good electrochemical performance as the anodes of Li (Na)-ions batteries (Xie et al., 2015). In particular, MoS2 nanoplate anodes in Li-ions batteries have been demonstrated to possess a capacity of 1,062 mAh g−1 (Cui et al., 2018). Previous investigations have attributed the high performance of the MoS2 anode behavior to its unique structural characteristics, such as the layered structure, which can provide more ions channels and storage compartments (Hu et al., 2014; Wang et al., 2017). However, the electrochemical performance of MoS2 is poor in divalent-ions (Mg, Zn) batteries. Thus far, few studies have been performed on the diffusion mechanisms and structural changes of the layer anode in M (Li, Na, Mg and Zn)-ions batteries to describe its different electrochemical behavior.

In this work, the adsorption and diffusion mechanisms of Li, Na, Mg, and Zn ions on the ordinary MoS2 (2H) structure were explored based on the density functional theory (DFT) and experimental analysis. In addition, the storage capability of Li, Na, Mg, and Zn ions in the MoS2 material were predicted and verified by experiments. Moreover, the structural deformation and electronic properties of MoS2 during ions intercalation were also investigated. Our theoretical and experimental results show that when used as the anode material of M-ions (Li and Na) batteries, MoS2 had significantly less variation in volume during the cycling process, while in Mg and Zn-ions batteries, MoS2 exhibited poor electrochemical performance due to the high migration energy barrier and low adsorption energy caused by changes in the structural properties based on first-principles. Furthermore, the MoS2 anode was found to undergo a phase transition from 2H to 1T during the intercalation of Li and Na ions, which was difficult to occur during the intercalation of Mg and Zn ions. This phase transition contributes to improving the performance of the MoS2 anode in M-ions batteries. Our results facilitate the understanding of the mechanisms of ions diffusion and structural changes of layer materials and provide useful information for designing high-performance anode materials, especially multivalent M-ions batteries.

Computational and Experimental Methods

Computational Methods

First-principles were used to describe the ions behavior in the anode based on density functional theory (DFT) with the Vienna ab initio simulation package (VASP) code (Kresse and Furthmüller, 1996). In addition, Perdewe-Burkee-Ernzerhof (PBE) generalized gradient approximation and the projected augmented wave (PAW) method were used to describe the ion-electron interactions in our systems (Perdew and Yue, 1986; Filippi et al., 1996). In this study, the plane-wave cutoff energy was set to 450 eV, and van der Waals corrections (optPBE-vdW) were adopted during structural optimization for the layer materials, and the vdWs interactions were described exactly by using DFT-D3 correction method of Grimme's scheme (Grimme et al., 2011). During the optimization, the Brillouin zone was represented by Monkhorst-Pack (MP), and the k-point mesh of 8 × 8 × 2 was used. Finally, the ion migration energies were acquired using the climbing-image nudged elastic band (CINEB) method (Henkelman et al., 2000; Yao et al., 2017). In our calculation, the structural optimization was considered complete when the force convergence criterion was <0.03 eV Å−1 and the total energy per unit cell was within 10−5 eV. CI-NEB calculations were performed with linear interpolating 5 images between the initial and final states of the diffusion paths, and spring constants is set as −5. The geometry and energy of the images were then relaxed until the largest norm of the force orthogonal to the path was <0.03 eV Å−1. To estimate the adsorption energy of M-ions on the MoS2 anode material, the adsorption energy was calculated by the equation: Eb = (Etotal − EMoS − nEM)/n, where Etotal is the ground energy of M-ions adsorbed on MoS2; EMoS, the ground energy of MoS2; EM, the chemical potential of M (Li, Na, Mg, or Zn) atoms; n, the number of M atoms (Tian et al., 2020).

Experimental Methods

In a typical synthesis, 1.0 mmol ammonium molybdate [(NH4)6Mo7O24] and 28 mmol thiourea (CH4N2S) were mixed with 60 mL deionized water containing 1.0 g PVP, followed by stirring for 2 h and transferring the solution into a 100 mL Teflon-lined stainless-steel autoclave at 200°C for 24 h. After cooling down to room temperature, the black precipitate was obtained by centrifugation, washed several times with water and ethanol, dried at 80°C overnight, and collected as the MoS2 material (Li and Peng, 2018).

For the electrochemical measurement, CR2032 cells were assembled in an argon-filled glove box by mixing the active samples, super-P and polyvinylidene fluoride at a weight ratio of 7: 2: 1. Next, the slurry was coated on Cu foil and dried at 60°C under vacuum for 10 h. After cutting into 12 mm discs, the working electrodes were obtained. Meanwhile, lithium (Li) metal or metal sodium (Na) was used as the anode, and the separators were commercial polypropylene (Celgard 2500 membrane) for Li-ions batteries and glass fiber (Whatman) for Na-ions batteries. For Li-ions batteries, the electrolyte was obtained by dissolving 1 M LiPF6 in ethylene carbonate (EC) and dimethyl carbonate (DMC) with 1:1 vol %. For Na-ions batteries, the electrolyte was prepared by dissolving 1 M NaClO4 and 5% fluoroethylene carbonate (FEC) in EC and DMC with 1:1 vol %. The electrochemical performances were recorded on Land cell test station (CT2001A) within the potential range of 0–3 V (vs Li/Li+).

Results and Discussions

Our structures calculated have been fully relaxed using first-principles calculation. Figure 1 shows the optimized crystal structure of MoS2 (supercell), and the optimized lattice constants are a = b = 3.18 Å and c = 15.12 Å for unit-cells. In addition, the theoretical XRD pattern was calculated based on DFT, which was consistent with the experimental results shown in Figure 1C. In Figure 1D, it can be seen that there are four possible sites in the MoS2 crystal structure. The large adsorption energy of ions in anode plays a fundamental role in providing a high ion storage energy in M-ions batteries, and the determination of suitable adsorption sites is a premise for the first-principles prediction of adsorption energies based on DFT. Therefore, the adsorption of M ions (Li, Na, Mg, and Zn) on MoS2 was investigated in our work.

Figure 1

The MoS2 structure drawings viewed along (A) c-axis and (B) the stereo perspective. (C) The XRD pattern of the MoS2 sample. (D) Adsorption site of M-ions.

The greater the negative adsorption energy, the better the thermodynamics and the more favorable for the adsorption thermodynamically. In our calculation, the 3 × 3 × 1 supercell for MoS2 was used to characterize the adsorption energy, corresponding to MxMoS2. Herein, as shown in Figure 1D, the different adsorption sites of metal on the MoS2 had been considered, including top site where metal atom sits directly above Mo (Mo top), hollow site above the center of Mo-S hexagon (center Mo-S hexagon), bridge site at the middle of Mo-S bond (Mo-S bridge) and the site directly above S atom (S top). And the adsorption energy had been shown in Table 1. Furthermore, the adsorption energy of Li and Na ions was found to be lower than that of Mg and Zn ions, implying the storage of Li and Na ions in the MoS2 anode material, thus resulting in higher specific capacity for Li- and Na- ions batteries in Figure 2A. In Figure 2B, the adsorption energies of M ions (Li, Na, Mg and Zn) were seen to increase gradually with the increase of M-ions concentration, and the large coulomb repulsion became increasingly apparent in the adjacent positively charged M-ions due to the high metal adsorption concentration. It noted that the Na curve decreases the fastest among others with higher adsorption concentration. And it may caused by the strong coulomb repulsion between Na neighboring positively charged, the large Na ion radius and adsorption energy. In addition, the adsorption energy of M-ions was lower than the cohesive energy of metal, such as Li metals (−2.01 eV), in a real battery system to ensure a positive discharge potential. When the number of atom adsorption increases, the adsorption energy of M-ions was larger than the cohesive energy of metal, which may form clusters. In Figure 2B, the adsorption energy of Zn-ions was shown to be larger than that of the other ions, suggesting that a small amount of Zn-ions can be stored in the MoS2 structure, leading to the low voltage and capacity of Zn-ions batteries. As a typical Li-ions batteries, we predict the capacity of MoS2 for Li batteries theoretically. From the Figure 2B, it is found that one Li atom adsorbed on a unit cell of MoS2, corresponding to storage capacity of ~687 mAh g−1, the adsorption energy becomes smaller than the chemical potential.

Table 1

The adsorption energies (eV) with different sites in MoS2 structure.

Sites	Li	Na	Mg	Zn
Mo top	−1.24	−0.61	−0.31	−0.28
S top	−1.08	−0.29	−0.17	−0.12
Mo-S bridge	−1.16	−0.55	−0.28	−0.26
Center Mo-S hexagon	−2.54	−2.02	−0.97	−0.71

Figure 2

(A) The adsorption energies center Mo-S hexagon sites and (B) the adsorption energies with M-ions concentration.

To study the physical origin of ions adsorption and anode performance, the charge density differences and density of states (DOS) of the adsorption structures were calculated and shown in Figures 3, 4. Figure 3 shows the charge density differences obtained by subtracting the total electron densities of MoS2 and isolated metal (Li, Na, Mg, and Zn) atom in center of Mo-S hexagon sites from that of the MxMoS2 structure. The isovalue was set as 0.03 eÅ−3, charge depletion was in green and accumulation in red. the distribution was similar in the charge density difference between Li ions and Na ions. Additionally, the charge rehybridization during the intercalation of Mg and Zn ions was obviously greater than that of Li and Na ions. Furthermore, a considerable alteration can be observed in the charge accumulation region between Li/Mg ions and S atoms, because these atoms are closer to one layer of sulfurs in the MoS2 structure. However, in Figure 3B, the transferred electrons were shown to be largely localized for Na ions, thus reducing the energy consumption during Na-ions diffusion.

Figure 3

The charge density difference and bader charge transfer of M-ions in MoS2 structure. (A) Li-ion, (B) Na-ion, (C) Mg-ion, (D) Zn-ion.

Figure 4

The projected density of states in MoS2 structure during Li (A), Na (B), Mg (C), and Zn (D) intercalation.

To further study the electronic structure of MoS2 during the intercalation of M-ions, the DOS was investigated (Figure 4), with the Fermi level set to zero. In Figures 4A,B, the band gap of MoS2 was seen to vanish during the adsorption of Li and Na-ions in contrast to the existence of the band gap in the adsorption of Mg and Zn-ions on MoS2. For comparison, the DOS of pure MoS2 had be repeated (Figure 5), which was in accordance with the reported previously (Hao et al., 2018; Chen et al., 2020). These results indicated that the semiconductor MoS2 may be transformed into metal during the insertion of Li and Na ions, while the properties of the MoS2 semiconductor remained unchanged during the insertion of Mg and Zn ions (Figure 4D), resulting in low electronic conductivity in divalent-ions batteries. Additionally, the conduction bands (CBM) were dominated by the Mo orbitals based on the DOS results. Meanwhile, the adsorbed ions (Li and Na) showed a discernible common peak between −3 and −1.5 eV with the S orbitals of MoS2, suggesting the hybridization between the adsorbed atoms and the MoS2 anode. However, the atoms adsorbed by Zn ions did not exhibit such characteristics. This is consistent with charge transfer analysis, where Li and Na loses electrons while MoS2 gains electrons with an itinerant feature. The electronic structures of MoS2 adsorbed by Li and Na ions were metallic, which can ensure a good electrical conduction, while the MoS2 adsorbed by Zn ions was shown to have semi-conductivity according to the electronic structure, leading to poor electrical conduction in Zn-ions batteries.

Figure 5

The projected density of states in MoS2 structure.

For the MoS2 structure, the migration pathways of M-ions and the corresponding migration energy barriers were calculated using the CINEB method and shown in Figure 6. In Figure 6A the migration pathways of M-ions (Li, Na, Mg and Zn) were seen to move from the center Mo-S hexagon site to the adjacent Mo top site. Based on our DFT results, the migration energy barriers of M-ions were estimated and shown in Figure 6B, which were 0.22, 0.13, 0.47, and 0.61 eV for Li, Na, Mg and Zn ions, respectively. The energy barrier of Na ions was lower than that of the other ions and in the order of Zn>Mg>Li>Na, which was consistent with previous reports (Shu et al., 2016; Sun et al., 2018). It is worth noting that the migration barrier was lower than that of Li-ions in Na-ions, possibly due to the distribution of charge density differences. Besides, the divalent nature significantly induced charge rehybridization during the intercalation of Mg or Zn-ions, leading to sluggish mobility. These results agreed with the electronic structure analysis. Generally, the migration of ions inside the anode can directly describe the charging and discharging rates (C rate) for M-ions batteries, and thus the migration barrier of ions is always considered as a desirable design parameter.

Figure 6

The M-ions migration pathways (A) and corresponding migration energy (B) in MoS2 structure.

Furthermore, the electrochemical experiments were performed for M-ions batteries, and the electrochemical performances of Li- and Na-ions batteries with the MoS2 anode at different current densities were shown in Figures 7A–C. Even at the high current density of 2 A g−1 for Li-ions batteries and 10 A g−1 for Na-ions batteries, the discharge platform was still in good condition and could ensure effective ion/electron transmission. In Figure 7C, the capacity of the MoS2 electrode was seen to remain at a high value when the current density returned to 0.1 A g−1 for both Li- and Na-ions batteries, indicating the good rate capability. The long cycle performance of Li- and Na-ions batteries is shown in Figure 7D, with a good coincidence observed for the charge capacity and discharge capacity. In Figure 8, it can be seen that, compared with Li- and Na-ions batteries with MoS2, Mg- and Zn-ions batteries with the MoS2 anode, especially for the Zn-ions batteries, had very poor electrochemical performance. All these results were consistent with the theoretical calculation results.

Figure 7

Galvanostatic charge-discharge curves of MoS2 at varied current density for (A) Li-ions battery anode and (B) Na-ions battery anode. (C) Rate properties of the MoS2 anode for Li- and Na-ions battery. (D) Cycling properties of MoS2 anode for Li- and Na-ions battery.

Figure 8

Galvanostatic charge-discharge curves of MoS2 for (A) Mg and (B) Zn-ions battery.

Structural Transformation

Undoubtedly, the MoS2 anode underwent a phase transition from 2H to 1T during Li ions intercalation, which has been confirmed by theoretical calculations and experimental analysis (Wang et al., 2013; Du et al., 2016; Zhu et al., 2019). Figures 9A,B shows the 2H and 1T structures of MoS2, with a trigonal prism shape for 2H and an octahedron for 1T in the Mo coordination structure. In order to describe the phase transition in the MoS2 anode during M-ions intercalation, the energy difference between the 2H and 1T structures of MoS2 was calculated under different M-ions concentrations based on first-principles, and shown in Figure 9C. The energy of 2H was lower than that of the 1T structure at the early intercalation stage of M-ions, while the 1T structure obtained relatively more energy at a higher M-ions concentration, suggesting the occurrence of the transition from 2H to 1T. For Li-ions batteries, Li ions were adsorbed in the tetrahedral coordination center of the S-S, leading to the charge transfer from Li ions, which was confirmed by the above electronic structure analysis (Figures 3, 4). In addition, the Mo atom was found to be surrounded by six S atoms after lithiation, implying that the Mo coordination structure may be transformed from 2H to 1T due to the intercalation of Li ions and the transfer of electrons. The 1T-type MoS2 possessed a higher conductivity than the 2H-type MoS2 due to its disordered octahedral structure. Therefore, this conversion of MoS2 can improve the electronic conductivity. Moreover, the 1T structure of MoS2 was more stable in energy when intercalating Li ions. In Figure 9, the phase conversion of 2H to 1T was also shown to be faster in Na ions than in the other M-ions, corresponding to a relatively lower concentration of Na-ions, which was caused by the larger ionic radius and electron transfer in Na-ions batteries. Therefore, the total charge transfer from M-ions to MoS2 was calculated and shown in Figure 9D. However, in the case of a high concentration of Zn-ions and an abundant intercalation of Zn2+, the energy was still higher in 2H than in the 1T structure, indicating that it is difficult to convert 2H to 1T in the MoS2 structure. In Figure 9D, the Mo coordination structure showed no significant change during the intercalation of Zn-ions, with a small amount of charge transfer from Zn-ions to MoS2. The aforementioned results suggest that the phase conversion from 2H to 1T is beneficial to improve the performance of the MoS2 anode.

Figure 9

The 2H (A) and 1T (B) for MoS2 anodes structure. (C) Energy difference between the 2H and 1T structure of MoS2 with Metal ions concentration and (D) Charge transfer with M-ions concentration in 2H and 1T MoS2 structure.

Conclusion

In summary, the intercalation mechanism of metal ions (Li, Na, Mg, and Zn) and the intrinsic properties of MoS2 as an anode material in M-ions batteries were investigated by experiments and first-principles calculations, and the ground state properties, geometrical and electronic structures, as well as the intercalation mechanism of M-ions in MoS2 were explored by DFT. Based on the calculation of adsorption energies, the Li and Na ions were stored in the MoS2 anode material due to the low adsorption energies, which may result in higher specific capacity for Li- and Na-ions batteries, in contrast to a weak adsorption strength of MoS2 for Mg and Zn ions, which is not conducive to the storage of anode. These results were also confirmed by the electrochemical performance experiments. Moreover, the density of states and charge density differences were investigated to explore the intercalation mechanism of M-ions. Our results suggest the hybridization between the Li/Mg-ions and the MoS2 anode and a considerable alteration in the charge accumulation region due to the amount of charge transfer. However, such characteristics were not observed in the atoms adsorbed by Zn ions. Furthermore, the migration of M-ions was used to describe the charging and discharging rates (C rate) for M-ions batteries, and the energy barrier was shown to be in the order of Zn>Mg>Li>Na. Mg and Zn-ions had poor electrochemical performance due to weak adsorption and high energy barrier of ions, which also was confirmed by the electrochemical performance experiments. Interestingly, the phase conversion from 2H to 1T was found to occur during the intercalation of Li and Na ions, which may induce the high adsorption capacity and electron transfer, resulting in the high performance of the anode. However, the phase transformation failed to occur in Mg and Zn-ions batteries. In the future, we may improve the performance of multivalent-ions batteries by regulating the phase transformation of MoS2. For example, the intercalation of molecules between the layers can increase the interlayer spacing and facilitate the intercalation of ions and the transfer of electrons, thereby making it more prone to phase change. Overall, the intercalation mechanism of M-ions can help optimize the design of layered structure for the high performance of M-ion batteries, and shed light on more efficient battery technologies.

Data Availability Statement

The original contributions generated for the study are included in the article/supplementary material, further inquiries can be directed to the Corresponding author.

Author Contributions

BBC developed the idea for the study. JBZ did the analyses and wrote the paper. XDL, JJZ, HL, BWH, JQZ, and SMJ conceived of the study, designed the study, and collected the data. All authors analyzed the data and were involved in writing the manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.