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

原创 公众号 LAMMPS交流站icon-default.png?t=LA92javascript: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交流站

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值