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

2026/6/16 16:29:34

避开坑！用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堆垛应确保：

基矢a3的z分量需包含层间距缓冲
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.5 & # 第二层原子 basis 0.333 0 0.5 & basis 0.5 0.5 0.5 & basis 0.833 0.5 0.5

关键参数对照表：

参数	错误值	修正值	物理意义
a3(z)	3.0	6.714	包含两倍层间距缓冲
basis z偏移	无	0.5	精确半周期定位

提示：使用Ovito可视化时，开启"Coordination Analysis"可快速检测原子重叠

2. AB堆垛的位移陷阱：1/3还是0.5？

2.1 典型错误：直接平移导致的错位

AB堆垛需要上层原子位于下层六元环中心，但以下代码存在严重问题：

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.5 & # 错误位移 basis 0.833 0.5 0.5

这种简单xy平面0.5偏移会导致：

破坏原始六方对称性
实际堆垛方式既非AA也非AB
力学性能计算结果失真

2.2 正确实现：严格遵循晶体学位移

AB堆垛的黄金法则是上层原子应在下层六元环中心，对应位移量为(1/3, 1/3)。修正代码如下：

basis 0 0 0 & basis 0.333 0 0 & basis 0.5 0.5 0 & basis 0.833 0.5 0 & basis 0.167 0.333 0.5 & # 正确AB位移 basis 0.5 0.333 0.5 & basis 0.667 0.833 0.5 & basis 0.833 0.833 0.5

验证技巧：

在Ovito中启用"Polyhedral Template Matching"
检查相邻层原子投影是否形成AB特征图案
测量最近邻碳原子间距应为≈2.46Å

3. 周期性边界条件的"幽灵原子"问题

3.1 错误现象：边缘原子的异常行为

当模拟盒尺寸与晶格常数不成整数倍时，会出现：

region box block 0 20 0 20 -5 5 # 任意尺寸 create_box 1 box region graphene block 0 20 0 20 0 1 create_atoms 1 region graphene

这会导致：

截断的六元环结构
边缘原子缺失键合
应力分布异常

3.2 专业解决方案：晶格对齐技巧

正确做法应确保：

模拟盒尺寸是晶格常数的整数倍
使用round函数自动适配

variable Lx equal 20 variable Ly equal 20 variable x equal ${Lx}/3 # 3是x方向晶格参数 variable y equal ${Ly}/1.732 variable xbox equal round(v_x) variable ybox equal round(v_y) region box block 0 ${xbox} 0 ${ybox} -5 5

进阶技巧：

使用lattice命令的origin选项微调位置
结合displace_atoms进行后期调整
采用delete_atoms去除边缘悬键

4. 验证与调试实战指南

4.1 Ovito诊断三板斧

配位数分析：理想石墨烯配位数应为3

ovito.modifiers.CoordinationAnalysisModifier(cutoff=1.8)

键序参数：检测sp²杂化完整性

ovito.modifiers.BondAngleAnalysisModifier()

层间距离分布：确认堆垛准确性

4.2 LAMMPS内置检查命令

在in文件中添加这些诊断命令：

compute coord all coord/atom 1.8 dump coord all custom 100 coord.dump id type c_coord run 0

典型问题与解决速查表：

现象	可能原因	解决方案
原子重叠报警	层间距设置错误	检查a3向量和basis z偏移
能量异常高	AB堆垛位移不准确	确认上层原子位于六元环中心
应力分布不对称	周期性边界条件不匹配	调整模拟盒尺寸为晶格整数倍
温度漂移	边缘原子未正确处理	使用fix shake或删除边缘原子