🔬 mesoscale physics

Position-dependent tight-binding model for Li impurities in monolayer and bilayer graphene

This paper presents a semi-empirical, position-dependent tight-binding model for lithium impurities in monolayer and bilayer graphene, calibrated against density-functional theory to accurately reproduce electronic band structures and elucidate how localized Li-induced perturbations modify low-energy properties through symmetry breaking and intervalley mixing.

Original authors: Hernan Aguirre, Hernan L. Calvo, Eduardo M. Perassi

Published 2026-07-21

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Hernan Aguirre, Hernan L. Calvo, Eduardo M. Perassi

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a world where the smallest building blocks of matter are like a giant, invisible trampoline made of carbon atoms. This is graphene: a sheet of carbon so thin it's just one atom deep, yet incredibly strong and conductive. It's the superstar of the material science world, promising to revolutionize everything from super-fast computers to the batteries that power our electric cars. But here's the catch: to make these batteries work, we need to stuff them with lithium ions (tiny, charged lithium atoms) to store energy. The problem is, we don't fully understand how these lithium ions dance around on or inside the graphene sheets. Do they just sit there? Do they change the way electricity flows? To answer this, scientists use a mix of super-computer simulations and clever math models. They want to know exactly how a single lithium ion messes with the "traffic" of electrons moving through the graphene, because that traffic is what makes our devices run.

In this study, a team of researchers from Argentina decided to build a new, smarter map to track this dance. Instead of treating lithium ions as static, boring dots that only sit in specific, perfect spots, they created a "position-dependent" model. Think of it like a video game where the terrain changes dynamically based on where the player stands. They tested this map on two different stages: a single layer of graphene (monolayer) and a double-layer sandwich (bilayer). By comparing their new map against ultra-detailed computer simulations (called Density Functional Theory), they found that their model is incredibly accurate. It successfully predicts how the energy levels of the electrons shift and how the "roads" for electricity open or close depending on exactly where the lithium ion is hovering.

Here is what they discovered: The lithium ion acts like a very shy, localized guest at a party. Its influence is strongest right next to it and fades away quickly as you move further out. On a single layer of graphene, the lithium ion mostly just changes the "height" of the energy floor for the nearby carbon atoms, and you can actually ignore the complex "hand-holding" (hopping) between the lithium and the carbon to get a good picture of what's happening. However, on the double-layer graphene, the story changes. The lithium ion is squeezed between two sheets, and now those complex hand-holds become essential; you can't ignore them without losing the plot.

The researchers also noticed something fascinating about the shape of the lithium's influence. On a single layer, the lithium's effect stretches out vertically, like a tall, thin egg (a prolate shape). But when it's sandwiched between two layers, the second layer squishes it down, making it look more like a flat pancake (an oblate shape). This happens because the top layer of graphene acts like a lid, confining the lithium's electric field.

Crucially, the paper rules out the idea that lithium ions interact strongly with each other over long distances. Their simulations show that as you spread the lithium ions further apart (making the "supercell" bigger), the interaction between them drops off rapidly. This means that in a battery with very few lithium ions (a dilute system), each ion acts mostly on its own, without worrying about its neighbors. The team is very confident in these results because their model matches the heavy-duty computer simulations almost perfectly, especially near the "Dirac points"—the special spots in graphene where electrons behave like massless particles.

While the model is a powerful tool for understanding the basics, the authors are careful to note that it is a simulation framework, not a physical experiment they performed in a lab. They suggest that this new, efficient way of modeling could be a great starting point for future studies on how lithium moves through graphene, which is a key step toward building better, faster-charging batteries. They haven't built a new battery yet, but they've handed the engineers a much better blueprint for how the ingredients behave.

Technical Summary: Position-Dependent Tight-Binding Model for Li Impurities in Monolayer and Bilayer Graphene

Problem Statement
Lithium-ion batteries (LIBs) rely on the efficient transport and storage of Li ions within graphitic electrodes. While graphene and bilayer graphene (BLG) are promising candidates due to their high conductivity and storage capacity, understanding the microscopic interaction between Li ions and the graphene electronic structure is critical. Specifically, there is a need to characterize how Li adsorption (on monolayer graphene, MLG) and intercalation (within AB-stacked BLG) modify the low-energy electronic properties, such as band structures near the Fermi level. Existing models often restrict impurities to static, highly symmetric sites, failing to capture the continuous nature of Li diffusion and the specific symmetry breakings induced by arbitrary impurity positions.

Methodology
The authors developed a semi-empirical, position-dependent tight-binding (TB) model to describe Li impurities in MLG and BLG. The approach integrates the following elements:

First-Principles Benchmarking: Density Functional Theory (DFT) calculations were performed using the Quantum ESPRESSO package (PBE functional, PAW pseudopotentials) to generate reference electronic band structures. Various supercell sizes ( $n=1$ to $4$) were used to model different Li dilution configurations, with Li ions placed at specific high-symmetry and intermediate positions ( $P_1$ to $P_4$ ) and varying vertical heights.
Model Framework: The TB model is built upon the Slonczewski-Weiss-McClure (SWM) $\pi$ -band framework. It introduces a single "pseudoorbital" for the Li ion to represent its interaction with the carbon $p_z$ orbitals.
Position Dependence: Unlike static models, this framework incorporates continuous position dependence. The Li-C interaction is modeled via onsite energy corrections ( $\epsilon$ ) and hopping terms ( $\gamma$ ) that decay with distance according to anisotropic Gaussian profiles. This allows the model to handle arbitrary Li positions within the supercell without increasing the number of fitting parameters.
Fitting Procedure: A Levenberg–Marquardt optimization algorithm was used to fit the model parameters ( $\alpha = \{\epsilon_{Li}, \epsilon_0, \gamma_0, \sigma, \sigma_z, \tau, \tau_z\}$ ) to the DFT band structures, specifically prioritizing accuracy near the Dirac points ( $\Gamma$ point of the supercell).
Analytical Decomposition: To interpret the results, the authors employed a partitioning scheme to separate the effective Hamiltonian into direct corrections (local potential and symmetry breaking) and indirect corrections (mediated coupling via virtual transitions).

Key Results

Model Accuracy: The fitted TB model accurately reproduces the DFT electronic bands near the Fermi level for both MLG and BLG across various supercell sizes and Li positions. It successfully captures symmetry breakings, intervalley mixing, and band-gap openings induced by the impurity.
Localization and Supercell Size: The perturbation introduced by the Li ion is found to be strongly localized. As the supercell size increases (lowering Li concentration), the energy scale of the band structure modifications decreases, as an increasing number of carbon atoms remain unperturbed.
Monolayer Graphene (MLG):
- For MLG, the dominant effect arises from local modifications to the carbon onsite energies.
- The fitting parameters for onsite corrections show weak dependence on supercell size.
- Crucially, the Li–C hopping terms were found to be negligible for MLG. A simplified model neglecting these hopping terms reproduces the low-energy electronic structure with nearly the same accuracy as the full model.
- The spatial profile of the impurity potential in MLG is prolate (elongated along the z-axis), consistent with the open surface environment.
Bilayer Graphene (BLG):
- In contrast to MLG, Li–C hopping terms are essential in BLG and cannot be neglected without significant loss of precision.
- The spatial profile of the impurity potential in BLG is oblate (compressed along the z-axis), attributed to vertical electrostatic confinement and symmetric screening between the graphene layers.
- The model successfully reproduces the complex band splitting and gap openings in BLG, which depend heavily on the specific sublattice symmetry relative to the Li position.
Interaction Energies: Calculations of interaction energies ( $E_{int}$ ) indicate that Li–Li interactions are stronger in MLG than in BLG. The presence of the second layer in BLG screens Li–Li interactions more effectively, causing the interaction energy to decay more rapidly with increasing supercell size.

Significance and Claims
The paper claims to provide an efficient, general framework for studying dilute Li impurities in graphene-based materials. By maintaining a constant number of fitting variables independent of the impurity position or supercell size, the model offers a "lightweight yet general" tool for transport simulations. The authors state that this framework constitutes a useful starting point for future investigations into Li diffusion and impurity-induced transport phenomena, such as determining current-induced forces (electromigration) near impurities. The work does not propose new experimental setups but rather establishes a theoretical tool to interpret electronic-ionic interactions and bridge the gap between first-principles calculations and large-scale transport simulations.

More like this