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

技术摘要：单层与双层石墨烯中锂杂质的位置相关紧束缚模型

问题陈述
锂离子电池（LIBs）依赖于锂离子在石墨化电极中的高效传输与存储。虽然石墨烯和双层石墨烯（BLG）因其高导电性和高储能容量而成为极具前景的候选材料，但理解锂离子与石墨烯电子结构之间的微观相互作用至关重要。具体而言，需要表征锂吸附（在单层石墨烯 MLG 上）和插层（在 AB 堆叠的双层石墨烯 BLG 中）如何改变低能电子特性，例如费米能级附近的能带结构。现有的模型通常将杂质限制在静态、高度对称的位点，无法捕捉锂扩散的连续性以及由任意杂质位置引起的特定对称性破缺。

方法论
作者开发了一种半经验的、位置相关的紧束缚（TB）模型，用于描述 MLG 和 BLG 中的锂杂质。该方法整合了以下要素：

• 第一性原理基准测试： 使用 Quantum ESPRESSO 软件包（PBE 泛函，PAW 赝势）进行密度泛函理论（DFT）计算，以生成参考电子能带结构。采用了各种超胞尺寸（$n=1$ 到 $4$）来模拟不同的锂稀释构型，并将锂离子放置在特定的高对称及中间位置（$P_1$至$P_4$）以及不同的垂直高度。
• 模型框架： 该 TB 模型基于 Slonczewski-Weiss-McClure (SWM) $\pi$ 能带框架构建。它引入了一个单一的“伪轨道”来代表锂离子与碳 $p_z$ 轨道的相互作用。
• 位置相关性： 与静态模型不同，该框架引入了连续的位置相关性。锂-碳相互作用通过随距离按各向异性高斯分布衰减的在位能量修正（$\epsilon$）和跳跃项（$\gamma$）来进行建模。这使得该模型能够处理超胞内任意的锂位置，而不会增加拟合参数的数量。
• 拟合程序： 使用 Levenberg–Marquardt 优化算法将模型参数（$\alpha = \{\epsilon_{Li}, \epsilon_0, \gamma_0, \sigma, \sigma_z, \tau, \tau_z\}$）拟合到 DFT 能带结构，特别优先考虑狄拉克点（超胞的 $\Gamma$ 点）附近的准确性。
• 解析分解： 为了解释结果，作者采用了划分方案，将有效哈密顿量分解为直接修正（局部势能和对称性破缺）和间接修正（通过虚跃迁介导的耦合）。

主要结果

• 模型准确性： 拟合后的 TB 模型能够准确重现 MLG 和 BLG 在各种超胞尺寸和锂位置下的 DFT 低能电子能带。它成功捕捉了由杂质引起的对称性破缺、谷间混合以及带隙开启。
• 局域化与超胞尺寸： 研究发现锂离子引入的扰动是强局域化的。随着超胞尺寸的增加（锂浓度降低），能带结构的修改能标随之减小，因为越来越多的碳原子保持未受扰动状态。
• 单层石墨烯 (MLG)：
  • 对于 MLG，主要效应源于对碳在位能量的局部修改。
  • 用于在位修正的拟合参数对超胞尺寸的依赖较弱。
  • 至关重要的是，研究发现对于 MLG 而言，锂-碳跳跃项可以忽略不计。一个忽略这些跳跃项的简化模型在重现低能电子结构方面，其准确度与全模型几乎相当。
  • MLG 中杂质势的空间分布呈长球形（沿 z 轴拉长），这与开放的表面环境一致。
• 双层石墨烯 (BLG)：
  • 与 MLG 不同，在 BLG 中，锂-碳跳跃项是必不可少的，若忽略这些项会导致精度显著下降。
  • BLG 中杂质势的空间分布呈扁球形（沿 z 轴压缩），这归因于垂直静电约束以及石墨烯层间的对称屏蔽作用。
  • 该模型成功重现了复杂的能带分裂和带隙开启现象，这些现象高度依赖于相对于锂位置的特定子晶格对称性。
• 相互作用能： 相互作用能（$E_{int}$）的计算表明，MLG 中的锂-锂相互作用比 BLG 中更强。BLG 中第二层的存在更有效地屏蔽了锂-锂相互作用，导致相互作用能在超胞尺寸增加时衰减得更快。

意义与主张
本文声称提供了一个高效、通用的框架，用于研究石墨烯基材料中的稀释锂杂质。通过保持与杂质位置或超胞尺寸无关的恒定拟合变量数量，该模型提供了一个“轻量级且通用”的工具，用于输运模拟。作者指出，该框架为未来研究锂扩散及杂质诱导的输运现象（例如杂质附近的电流诱导力/电迁移）提供了有用的起点。这项工作并非提出新的实验装置，而是建立了一个理论工具，用以解释电子-离子相互作用，并弥合第一性原理计算与大规模输运模拟之间的鸿沟。