手机版

LAMMPS模拟计算：从输入文件到结果分析的全流程经验

发布时间：2026-07-07 来源：科研学术网

小中大

字号：

LAMMPS是材料科学领域使用最广泛的MD软件——免费、开源、并行效率高、力场支持全面。但LAMMPS的学习曲线陡峭，输入脚本语法复杂，很多新手在第一步就被劝退。本文从实际项目出发，分享LAMMPS模拟的全流程经验，重点是”怎么设参数才靠谱”。

一、LAMMPS输入文件结构

1.1 四大模块

一个完整的LAMMPS输入文件包含：

1.2 单位制选择

单位制	长度	时间	能量	温度	适用
metal	Å	ps	eV	K	金属/材料(最常用)
real	Å	fs	kcal/mol	K	生物/有机
si	m	s	J	K	工程尺度(极少用)
lj	σ	τ	ε	1	无量纲

经验：90%的材料模拟用metal单位。注意metal单位的时间步长是ps，如果用real单位则是fs——别搞混了。

二、结构构建经验

2.1 LAMMPS内置建模

lammps

复制

# 创建FCC Cu超胞
lattice        fcc 3.615
region         box block 0 20 0 20 0 20
create_box     1 box
create_atoms   1 box

2.2 外部工具建模

工具	用途	输出格式
Atomsk	晶体/缺陷/界面建模	LAMMPS data
OVITO	可视化+结构编辑	多格式
VMD	生物分子建模	LAMMPS data
MS	图形界面建模	导出data文件
Packmol	液体/混合体系	坐标文件

经验：Atomsk是LAMMPS建模的神器——一行命令创建位错、晶界、层错等缺陷结构：

bash

复制

atomsk --create fcc 3.615 Cu Cu.lmp
atomsk initial.lmp -dislocation 0.5*box 0.5*box 1 0 0 0.5 0.5 final.lmp

三、力场设置经验

3.1 常用力场速查

lammps

复制

3.2 力场文件路径

踩坑经验：LAMMPS找不到力场文件是最常见的报错。解决方法：

用绝对路径：pair_coeff * * /home/user/potentials/Cu_u3.eam.alloy Cu
或在运行时指定：lmp -in in.run -var potfile Cu_u3.eam.alloy

四、MD模拟参数详解

4.1 时间步长

lammps

复制

timestep       0.001    # metal单位: 0.001 ps = 1 fs

体系	步长(fs)	说明
金属(无H)	1-2	标准选择
含C/O	1.0	通用
含H	0.5	C-H振动
高温	0.5	剧烈振动
ReaxFF	0.25	化学反应

4.2 恒温器/恒压器

lammps

复制

# NVT (Nose-Hoover)
fix 1 all nvt temp 300 300 0.1

# NPT (Nose-Hoover)
fix 1 all npt temp 300 300 0.1 iso 0.0 0.0 1.0

# Berendsen (平衡阶段)
fix 1 all temp/berendsen 300 300 0.1
fix 2 all press/berendsen iso 0.0 0.0 1.0

经验：

Tdamp/Pdamp参数：0.1 ps(metal单位100步)是标准选择
太小(0.01)→温度震荡；太大(1.0)→响应慢
NVT用Nose-Hoover；NPT用Nose-Hoover或Parrinello-Rahman

4.3 输出设置

lammps

复制

# 热力学信息输出
thermo         100              # 每100步输出
thermo_style   custom step temp pe ke etotal press vol density

# 轨迹输出
dump           1 all custom 1000 dump.lammpstrj id type x y z vx vy vz

# 重启文件
restart        10000 restart.*.bin

五、典型模拟流程

5.1 金属熔化-凝固

lammps

复制

参考价：3000-5000元

5.2 拉伸模拟

lammps

复制

参考价：2000-4000元

5.3 热导率计算(NEMD)

lammps

复制

# 分为热源/热库/中间区
region         hot block INF INF INF INF 0 2
region         cold block INF INF INF INF 18 20
group          hot region hot
group          cold region cold

# NVE + 温度重标
fix            1 all nve
fix            2 hot temp/rescale 1 310 310 1 1
fix            3 cold temp/rescale 1 290 290 1 1

# 计算温度梯度
compute        chunks all chunk/atom bin/1d z 0.5
fix            4 all ave/chunk 100 1000 100000 chunks temp file temp_profile.txt

run            1000000      # 1 ns

参考价：4000-8000元

六、后处理分析

6.1 RDF计算

lammps

复制

compute        rdf all rdf 100
fix            1 all ave/time 100 100 10000 c_rdf file rdf.txt mode vector

6.2 MSD计算

lammps

复制

compute        msd all msd
fix            1 all ave/time 100 1 100 c_msd[4] file msd.txt

经验：扩散系数D = MSD斜率/6。在LAMMPS中也可以直接用：

lammps

复制

compute        1 all msd
fix            2 all ave/time 1 1 1000 c_1[4] file msd.txt
# 后处理: D = slope(msd vs t) / 6

6.3 用OVITO分析

OVITO是LAMMPS轨迹分析的最佳工具：

可视化原子轨迹
计算配位数、CNA(公共近邻分析)
位错识别(DXA)
应变张量分析
制作高质量动画

七、并行与加速

7.1 并行设置

bash

复制

# MPI并行
mpirun -np 32 lmp -in in.run

# 按原子分区
# LAMMPS自动按空间分区，不需要手动设置

经验：

金属体系：1核/500-1000原子
ReaxFF：1核/200-500原子（计算量大）
不要过度并行——核数超过原子数/100后效率下降

7.2 GPU加速

lammps

复制

# LAMMPS GPU包
package        gpu force/neigh 1
pair_style     gpu/lj/cut 10.0

经验：GPU加速对短程力(LJ/EAM)效果3-5倍，对长程力(PME)效果不大。

八、项目报价参考

项目类型	体系大小	模拟时间	参考价
结构优化+RDF	<5000原子	100-500ps	1500-3000元
熔化/凝固模拟	<10000原子	500ps-1ns	3000-5000元
拉伸/压缩模拟	<5000原子	200-500ps	2000-4000元
扩散系数	<10000原子	1-5ns	3000-6000元
热导率(NEMD)	<10000原子	1-5ns	4000-8000元
ReaxFF反应模拟	<3000原子	100-500ps	5000-10000元
完整材料研究	不限	1-10ns	5000-15000元