lammps高熵合金熔炼过程分子动力学模拟

原创 公众号 LAMMPS交流站javascript:void(0);

   大家好，前面小编讲解了一期层合结构的建模，小编本期基于之前的层合结构，设置5种金属建立高熵合金模型，进行一个高熵合金熔炼过程的分子动力学模拟。

   不足之处请各位指出，希望能给读者提供一定的帮助，同时希望专业前辈提出不足，小编会和大家共同学习。

  

01

基本思路

1.建立初始构型

2.能量最小化

 

3.保温固熔及保温固凝

 

02

完整代码

  

完整in文件如下：

# Ni包覆Bi:1能量最小化、2弛豫、3升温融化、4保温、5降温、6弛豫
#------------------start-------------------
#读取模型数据
units      metal
boundary        p p p                        #固定边界条件
atom_style      atomic
read_data       cenghe.data                #读取建模数据
#------------------------------------------
#热力学信息输出
timestep        0.001                        #时间步长
thermo      500                     #100步输出一次
thermo_style    custom step temp vol pe     #输出类型：pe-势能，ke-动能
thermo_modify   lost ignore flush yes 
#------------------------------------------
#近邻信息
neighbor        1.0 bin                      #近邻原子距离
neigh_modify    every 1 delay 0 check yes
#------------------------------------------
#设置势函数
pair_style      lj/cut 10.0             #lj势/截断半径
pair_coeff    * * 0.03246 2.8932          #epsilon:ev  sigma:A
#------------------------------------------  
#温度初始化
velocity      all create 300 898955 dist gaussian    #原子速度调控温度
#------------------------------------------
#1能量最小化    
dump            1 all atom 5 1min.xyz        #能量最小化输出文件
min_style       cg
minimize        1e-15 1e-15 1000 1000
undump          1
reset_timestep  0
#------------------------------------------
#2nvt系综固熔
dump      2 all atom 100 2nvt_hot.xyz  #升温过程输出文件
fix             2 all nvt temp 2000 2000 0.01 #nvt保温：2000
run          10000                 #升温步长
undump          2
unfix           2
reset_timestep  0                            #重置时间步长
#------------------------------------------ 
#3nvt系综固凝
dump      3 all atom 100 3nvt_cool.xyz  #升温过程输出文件
fix             3 all nvt temp 300 300 0.01 #nvt升温：300-2000
run          10000                 #升温步长
#------------------end---------------------

文中所用到的data文件已在之前教程层合结构建模给出。

03

运行结果

  在cmd窗口运行in文件，将生成的轨迹文件拖到OVITO软件中，并在OVITO软件中进行渲染，就产生了如下结果：

 

注意

1.文中所需的data文件及势文件已上传至QQ群-lammps交流站。

更多案例请关注公众号lammps交流站