干货分子动力学 干货详情

lammps中文教程13

春风得意马蹄疾33
分子动力学LAMMPS能量最小化结构弛豫共轭梯度法最速下降法原子模拟计算材料学初始模型优化原子重叠

1.13 minimize命令
进⾏lammps模拟,如果初始模型不合理、原⼦间距离过近,容易导致体系能量过⾼,原⼦速度过快致使
原⼦⻜出box。
前⽂介绍过“delete_atoms”命令可删除重叠或距离过近原⼦,但在⼤多数模型中,并不需要删除原⼦,有
些模型也不允许删除原⼦,这时候就需要⽤到minimize命令。
1. minimize命令
minimize命令通过⼒的计算排开距离过近的原⼦,使体系的能量降到⼀个较低的⽔平,便于后续的弛豫
操作。
minimize语法为:
minimize etol ftol maxiter maxeval
后⾯四个参数为minimize运⾏停⽌条件,其中etol、ftol分别为能量和⼒的容差,maxiter、maxeval为迭
代最⼤步数。
在最⼩化过程中,达到四个条件中的任意⼀个时,最⼩化即停⽌运⾏。
⼒或者能量已经降低到⼀个较低⽔平,或者⼒或能量还没有降低到⼀个低值,但设定的迭代步数已经⽤
完,都会停⽌最⼩化。
2. min_style命令
minimize最⼩化采⽤的算法由min_style命令确定,min_style语法为:
min_style style
style = cg or hftn or sd or quickmin or fire or fire/old or spin or spin/cg or
spin/lbfgs
最常⽤的为共轭梯度法(cg)、最速下降法(sd)。

3. 使⽤举例
min_style cg
minimize 1.0e-4 1.0e-6 100 1000
如果已经进⾏了最⼩化,但是后续弛豫仍然会报错,可查看log⽂件,排查minimize运⾏停⽌的条件,适
当调⼩能量和⼒的容差或者加⼤迭代步数。
如以上语句可改为:
min_style cg
minimize 1.0e-15 1.0e-15 10000 10000
mimimize可有效降低出错概率,是⼀个谁⽤谁知道,⽤了都说好的命令。

评论

0/1000发布评论
全部评论

春风得意马蹄疾

关注
TA的主页

干货推荐