LAMMPS是材料体系分子动力学模拟的主流工具，特别适合金属材料、聚合物、纳米材料等大规模体系的模拟。结合多个材料模拟项目经验，聊聊LAMMPS计算中的关键设置与调优经验。
1. 为什么材料体系用LAMMPS更多
GROMACS在生物分子领域是绝对主流，但用到材料体系(金属、陶瓷、聚合物熔体等)时，LAMMPS的优势就体现出来了。最核心的原因是：LAMMPS对力场的支持更灵活，几乎任何原子间相互作用势都可以通过pair_style来配置，而GROMACS的力场体系主要针对生物分子设计，对金属材料支持有限。
我做铝合金沉淀析出模拟的时候，最初尝试用GROMACS，发现很难找到合适的力场来描述Al-Mg-Si体系中原子的相互作用。后来换用LAMMPS，用EAM(Embedded Atom Method)势，问题迎刃而解。这个经历让我意识到：选软件不能只看名气，要看它是否适合你的体系类型。
2. 力场选择：几种常用势的取舍
LAMMPS里的力场通过pair_style命令指定。金属材料最常用的是EAM和MEAM势：EAM势计算量适中，对FCC和BCC金属的熔点、弹性常数等性质描述较好;MEAM势是EAM的改进版，考虑了方向性键合，对共价键较强的材料(如硅、碳)更合适。
我做铁素体钢中碳扩散的模拟时，用的是EAM势(Fe-C的嵌入原子势有现成的库可用，比如”Fe_C.eam”)。但后来做石墨烯力学性能的模拟时，EAM完全不适用(EAM假设金属键，石墨烯是共价键)，换成AIREBO势才得到合理结果。力场选错，后面所有的计算都是浪费时间。
一个实用建议：LAMMPS的potentials目录下有大量预定义的势文件，先看看有没有适合你体系的，不要自己从头写参数。我从头写过一次势函数参数，花了两周时间调参，最后发现potentials目录下早就有人做好了更好的版本——这个教训让我后来先搜库再动手。
3. 输入脚本：in文件的编写逻辑
LAMMPS用in文件(输入脚本)来定义整个模拟流程。一个标准的in文件包含以下部分：初始化(units、dimension、boundary)→ 定义原子(read_data或lattice + create_atoms)→ 设置力场(pair_style、pair_coeff、bond_style等)→ 定义输出(thermo、dump)→ 运行(run)。
我写in文件有个习惯：把参数集中放在文件开头，用variable命令定义，后面引用变量名。这样需要调整参数的时候，只需要改文件开头的variable行，不用在几百行脚本里逐个找。比如：
variable T equal 300
variable dt equal 0.001
variable nsteps equal 1000000
units metal
timestep ${dt}
fix 1 all nvt temp ${T} ${T} 0.1
run ${nsteps}
这个写法还有一个好处：用同一个in文件跑不同温度或者不同步长的时候，改一个variable就够了，不容易改漏。
4. 积分算法与系综选择
LAMMPS通过fix命令来设定积分算法和统计系综。NVE(微正则系综)用fix nve;NVT(正则系综)用fix nvt;NPT(等压等温系综)用fix npt。选择哪个系综取决于你要模拟的物理场景。
我做纳米压痕模拟的时候，用NPT系综(压头区域的z方向用fix npt允许体积变化，样品底部固定)。但如果做理想晶体中的点缺陷扩散，NVE系综其实更合适——因为它严格守恒能量和动量，扩散系数的计算结果更干净。系综选错，统计结果可能有系统性偏差。
5. 大规模并行：性能优化的实际经验
LAMMPS的并行效率在同类软件里算很好的，但要用好也需要一些经验。核心参数是处理器网格划分(processors命令)：LAMMPS把模拟盒子划分成若干子域，每个CPU核心负责一个子域。划分方式直接影响通信开销。
我的经验是：让processors的划分尽可能接近立方体的子域形状。如果模拟盒子在x方向很长，但processors设成了Px=1. Py=8. Pz=8.子域就会在x方向很长、y和z方向很薄，子域之间的通信效率会下降。用”processors * * *”让LAMMPS自动优化划分，通常能得到不错的效果。如果还嫌慢，可以试试用”bond/improved”或者”neigh_modify”调整邻居列表更新频率——这两个参数对计算效率的影响有时候比增加CPU核心数还明显。
