lammps多孔结构

原创   公众号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交流站