在以上LAMMPS脚本中，需要特别注意的是，命令“set type 1 type/ratio 2 0.2 7603”执行的结果与随机数7603有关，相同的随机数通过该命令生成的结果是一样的因此我们在生成下一个模型时需要提前对随机数进行更改。后面四行通过sed命令将拿到的随机用来修改FeNiCrCoCu.in文件，这样就可以保证在循环运行FeNiCrCoCu.in脚本时，生成的结构都是不同的。本期文章，我们会详细解读上一期贴出的脚本，以便让读者更深的了解LAMMPS程序的运行过程。
上一期我们介绍了对于格点原子无序的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
最低0.47元/天 解锁文章
1570
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
