lammps中文教程13
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可有效降低出错概率,是⼀个谁⽤谁知道,⽤了都说好的命令。
评论
春风得意马蹄疾

干货推荐

【代码分享】费米面嵌套(Fermi surface nesting)函数计算
张阳小号1
651 2

CP2K计算能带结构与DOS图
齐振宏
310 2

有髓神经纤维超微结构
Mr弘🔬
56 5

顶级期刊子刊的含金量怎么样?
九妹
41 1

Codex + Zotero + Obsidian,研究生真的可以构建自己的个人论文知识库
海山了
117 2

提取文献数据点,你只知道GetData?
木若林溪
557 1
