分享一个金球融化过程的lammps模拟案例。
(1)建模
设置晶格后建立一个15*15*15晶格的微信公众号怎么赚钱box,box内全系统运维工资一般多少部填充Au原子。
在bo系统运维包括哪些内容x内设置一个球形的区域(region),使用delete_atoms命初始化电脑时出现问题令删除球形区域外的Au原子,实现金球的建模。
(2)势函数设置
pair_styleeam
pair_coeff**Au_u3.eam
本例只含有1种原子,使用eam势,设置方式比较简单,Au_u3.eam后不需要加原子列表。
(3)nvt系综下的升温
fix 2 all nvt temp 298 2400 0.1
温度由298K升温至2400K,升温时间=模拟步数*timestep=150ps,因此可得升温速率=(2400-298)K/150ps。
设置不同的模拟步数可实现不同升温速率下的模拟。
in源代码如下:
# 金球融化模拟
#模拟基本参数
dimension3
unitsmetal
atom_styleatomic
boundaryppp
#设置晶格参数
latticefcc4.0782
regionboxblock015015015
create_box1box
#在box内填充Au原子
create_atoms1box
mass1196.97
#设置一个球形区域,取球形的外侧
regionsphsphere7.07.07.04.0sideout
#删除球形外侧区域原子,只保留球内原子
delete_atomsregionsph
#设置ema势
pair_styleeam
pair_coeff**Au_u3.eam
#近邻列表参数定义
neighbor3.0bin
#每20步更新一次近邻列表
neigh_modifyevery20delay0checkno
#能量最小化
minimize1.0e-81.0e-810000100000
min_stylecg
#模拟步长
timestep0.001
reset_timestep0
#热力学信息输出
thermo5000
#自定义输出:步数、势能、动能、 温度、体积、压力
thermo_stylecustomsteppeketempvolpress
#初始化温度298K
velocityallcreate29839849momyesrotyesdistgaussian
#1000步保存一次模拟结果
dump1dallcustom1000dump.meltidtypexyz
#nvt系综下由298K->2400K,实现融化
fix2allnvttemp29824000.1
#运行150000步
run150000
END
陆续推出lammps基础教程,敬请扫码关注微信公众号:lammps加油站。
发表评论