lammps 案例:金球融化模拟代码练习

分享一个金球融化过程的lammps模拟案例。

代码已经注释以供参考。本例主要包含的知识点

(1)建模

设置晶格后建立一个15*15*15晶格的微信公众号怎么赚钱boxbox内全系统运维工资一般多少部填充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加油站。


                                            lammps 案例:金球融化模拟代码练习