Lithium Intercalation in Graphene/MoS2 Composites: First-Principles Insights

As a storage material for Li-ion batteries, graphene/molybdenum disulfide (Gr/MoS2) composites have been intensively studied in experiments. But the relevant theoretical works from first-principles are lacking. In the current work, van-der-Waals-corrected density functional theory calculations are performed to investigate the interaction of Li in Gr/MoS2 composites. Three interesting features are revealed for the intercalated Gr/Li(n)/MoS2 composites (n = 1 to 9). One is the reason for large Li storage capacity of Gr/MoS2: due to the binding energies per Li atom increase with the increasing number of intercalated Li atoms. Secondly, the band gap opening of Gr is found, and the band gap is enlarged with the increasing number of intercalated Li atoms, up to 160 meV with nine Li; hence these results suggest an efficient way to tune the band gap of graphene. Thirdly, the Dirac cone of Gr always preserve for different number of ionic bonded Li atoms.

💡 Research Summary

The authors present a comprehensive first‑principles investigation of lithium intercalation in graphene/molybdenum disulfide (Gr/MoS₂) heterostructures, a material class that has attracted considerable experimental attention for lithium‑ion battery (LIB) anodes but has lacked systematic theoretical insight. Using density functional theory (DFT) with the Perdew‑Burke‑Ernzerhof (PBE) generalized‑gradient approximation and a Grimme D3 van‑der‑Waals correction, they model a 4 × 4 graphene supercell stacked on a 3 × 3 MoS₂ slab, separated by a 13 Å vacuum region to avoid spurious interactions. Lithium atoms are placed in the interlayer space between graphene and MoS₂, and the number of Li atoms (n) is varied from 1 to 9, representing realistic loadings observed experimentally.

Structural optimization reveals that Li atoms preferentially occupy sites that allow simultaneous coordination to carbon atoms of graphene and sulfur atoms of MoS₂, forming mixed Li–C and Li–S ionic bonds. The binding energy per Li atom, defined as E_b = (E_Gr/MoS₂ + n E_Li − E_Gr/Li_n/MoS₂)/n, increases monotonically from ~1.85 eV for a single Li to ~2.31 eV for nine Li. This trend indicates a cooperative stabilization mechanism: as more Li atoms are inserted, the electrostatic environment becomes more favorable, and charge transfer from Li to the host layers intensifies. Bader charge analysis quantifies an average transfer of ~0.85 e⁻ per Li atom, confirming that Li behaves essentially as a cation and that the donated electrons are delocalized over both graphene and MoS₂.

Electronic‑structure calculations show that pristine graphene’s gapless Dirac cone is preserved in the composite, but the presence of Li induces a modest band gap opening in the graphene-derived bands. The gap widens with Li concentration, reaching a maximum of ~160 meV for the nine‑Li configuration. Importantly, the Dirac point itself remains intact; the linear dispersion is shifted but not destroyed, indicating that the ionic nature of the Li–C interaction does not break the sublattice symmetry of graphene. Density‑of‑states (DOS) analysis corroborates this picture: the Fermi level moves into the conduction band as Li concentration increases, reflecting a transition toward metallic behavior. Charge‑density difference plots (Δρ) reveal electron accumulation around Li and depletion on neighboring C and S atoms, illustrating the formation of a built‑in electric field that both stabilizes the intercalated Li and modulates the electronic bands.

The authors connect these theoretical findings to experimental observations. The increasing per‑Li binding energy rationalizes the exceptionally high specific capacities (>600 mAh g⁻¹) reported for Gr/MoS₂ anodes, as each additional Li becomes energetically more favorable to host. The tunable graphene band gap suggests a dual‑functionality: beyond energy storage, Gr/MoS₂ could serve as a platform for graphene‑based electronic devices where a controllable gap is desirable. Moreover, the preservation of the Dirac cone implies that high carrier mobility is retained, which is advantageous for rapid charge‑discharge kinetics in LIBs.

Limitations of the study are acknowledged. Calculations are performed at 0 K on static structures, neglecting finite‑temperature effects, Li diffusion barriers, and explicit electrolyte interactions that are critical in real batteries. The authors propose future work involving ab‑initio molecular dynamics, nudged elastic band calculations for Li migration pathways, and explicit modeling of solid‑electrolyte interphase (SEI) formation to capture the full complexity of operating cells.

In summary, this paper provides three key insights: (1) the per‑Li binding energy in Gr/MoS₂ rises with Li loading, explaining the material’s high storage capacity; (2) intercalated Li opens a modest, tunable band gap in graphene while preserving its Dirac cone, offering a route to engineer graphene’s electronic properties; and (3) the Dirac cone’s robustness ensures that high electronic conductivity is maintained, supporting fast ion transport. These results not only deepen the fundamental understanding of Li‑graphene/MoS₂ interactions but also guide the rational design of next‑generation high‑capacity, high‑rate LIB anodes and multifunctional 2D electronic heterostructures.