LAMMPS实操系列（二）: 大量FCC-CoCrCuFeNi高熵合金建模与最稳定结构筛选

上一期我们介绍了对于格点原子无序的FCC-CoCrCuFeNi高熵合金理论模拟中的建模思想，以及LAMMPS的程序实现。本期文章，我们会详细解读上一期贴出的脚本，以便让读者更深的了解LAMMPS程序的运行过程。

下面就是程序实现过程的详细解读与分析：

units      metal      #定义单位系统为metal

boundary      p p p        #定义三个方向都为周期性边界条件

atom_style    atomic       #定义原子类型为atomic

timestep      0.001        #设定时间步长为1 fs

#构建模型

lattice       fcc 3.56       #定义fcc晶格常数

region        box block 0 10 0 10 0 10   #定义模拟盒子大小

create_box    5 box        #定义模拟盒子中最多可以包含5种原子

label       loop       #定义脚本循环起始标签为loop

variable      a loop 1000      #定义跟踪循环次数变量a为1至1000     

log        log.$a       #定义记录标准输出文件的文件名以次数结尾

create_atoms   1 box   #现在模拟盒子中将格点原子用第1种原子Fe代替

set    type 1 type/ratio 2 0.2 7603  #将1/5的Fe原子换成第2种原子Ni

set    type 1 type/ratio 3 0.5 7603  #将1/2的Fe原子换成第3种原子Cr

set    type 1 type/ratio 4 0.5 7603  #将1/2的Fe原子换成第4种原子Co

set    type 3 type/ratio 5