lammps模拟案例:晶体材料熔点的模拟
大家好,我是小马老师。
本文分享一个lammps模拟熔点的一个案例代码。
lammps模拟材料的熔点方法有很多,本文通过温度-势能曲线判断熔点。
使用npt系综控制升温,在融化过程中输出温度和势能,绘制温度-势能曲线,势能发生突变的温度即为熔点。
in文件代码:
#模型设置
units real
atom_style atomic
neighbor 1.0 bin
neigh_modify every 1 check yes
timestep 1.0
#建立模型
lattice fcc 5.58
region box block 0 5 0 5 0 5
create_box 1 box
create_atoms 1 box
mass 1 83.798
# 力场设置
pair_style lj/cut 7.0
pair_coeff 1 1 0.343 3.591
#输出设置
thermo 100
thermo_style custom step temp pe press vol density lx
dump 2 all custom 100 dump.xyz id type x y z
#能量最小化
fix 1 all box/relax iso 100
min_style cg
minimize 0.0 1.0e-8 1000000 10000000
unfix 1
#弛豫
fix 1 all momentum 1 linear 1 1 1 angular
fix 2 all npt temp 1 1 100 iso 100 100 1000
run 20000
unfix 2
unfix 1
#融化
fix 2 all npt temp 1 150 100 iso 100 100 1000
run 150000
unfix 2
绘制温度-势能曲线如下图所示,熔点大约120K左右:
公众号:lammps加油站