原创 公众号LAMMPS交流站
点击上方蓝字关注我们
大家好,我是YJ学长,本文的模型是小编在lammps模拟Ni的熔化过程中偶然生成了一帧多孔结构,小编觉得图挺好看的,因此做一期推送共同来讨论这一现象。
因能力有限,这种偶然生成的结构欢迎专业前辈能提出建议。
首先是本文建立的Ni金属团簇模型。
01
完整代码
完整in文件如下:
# Ni:1能量最小化、2弛豫、
#------------------start-------------------
#读取模型数据
units metal
atom_style atomic
read_data Ni.data #读取建模数据
#------------------------------------------
#热力学信息输出
thermo 100 #100步输出一次
thermo_style custom step temp vol pe #输出类型:pe-势能,ke-动能
neighbor 2.0 bin #近邻原子距离
#------------------------------------------
#设置势函数
pair_style eam #势类型
pair_coeff * * Ni.lammps.eam #选择用该势原子
#------------------------------------------
#1能量最小化
thermo 100
dump 1 all atom 5 1min.xyz #能量最小化输出文件
min_style cg
minimize 1e-15 1e-15 1000 1000
undump 1
reset_timestep 0
#------------------------------------------
#温度初始化
velocity all create 300 89895 #原子速度调控温度
#------------------------------------------
#2nvt系综弛豫
fix 1 all nvt temp 298 298 0.01 #弛豫参数
dump 1 all atom 100 2nvt_equ.xyz #弛豫输出文件
run 10000 #弛豫步数
02
运行结果
接着,在cmd窗口运行in文件,将生成的xyz文件拖到OVITO软件中,并在OVITO软件中进行渲染,就产生了如下结果:
弛豫过程出现的一帧多孔结构
更多案例代码,尽在公众号LAMMPS交流站