避开坑！用LAMMPS构建双层石墨烯（AA/AB堆垛）的3个常见错误与修正

避开坑！用LAMMPS构建双层石墨烯（AA/AB堆垛）的3个常见错误与修正

在分子动力学模拟中，石墨烯的双层结构建模看似简单，却暗藏诸多细节陷阱。许多研究者在初次尝试用LAMMPS构建AA或AB堆垛的双层石墨烯时，往往会遇到原子重叠、层间距异常或周期性边界条件导致的图像失真等问题。本文将聚焦三个最易被忽视的关键错误，通过对比错误与正确的代码实现，帮助您快速定位问题根源。

1. 层间距设置：被低估的"隐形杀手"

1.1 错误现象：看似合理的参数导致原子穿透

初学者常直接复制单层石墨烯的晶格参数，仅简单修改z方向基矢长度。例如以下典型错误代码：

lattice custom 1.421 a1 3 0 0 a2 0 1.732 0 a3 0 0 3.0 &
    basis 0 0 0 &
    basis 0.333 0 0 &
    basis 0.5 0.5 0 &
    basis 0.833 0.5 0

表面看a3向量的z分量已扩大，但实际上会导致：

• 上层原子穿透下层晶格
• 实际层间距偏离理论值(≈3.35Å)
• 后续模拟中非物理的原子相互作用

1.2 修正方案：精确控制垂直距离

正确的AA堆垛应确保：

1. 基矢a3的z分量需包含层间距缓冲
2. basis原子的z坐标需精确偏移

lattice custom 1.421 a1 3 0 0 a2 0 1.732 0 a3 0 0 6.714 &
    basis 0 0 0 &
    basis 0.333 0 0 &
    basis 0.5 0.5 0 &
    basis 0.833 0.5 0 &
    basis 0 0 0