应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】结构优化中，如何选择POTIM？
北京石油化工学院2026年研究生招生接收调剂公告
81/1返回列表
查看: 5325  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

Prettyswan

木虫 (小有名气)

[交流] 【求助】结构优化中，如何选择POTIM？ 已有6人参与

说明手册上说：
引用回帖:
Within all relaxation algorithms (IBRION=1,2 and 3) the parameter POTIM should be supplied in the INCAR file. For IBRION0, the forces are scaled internally before calling the minimization routine. Therefore for relaxations, POTIM has no physical meaning and serves only as a scaling factor. For many systems, the optimal POTIM is around 0.5. Because the Quasi-Newton algorithm and the damped algorithms are sensitive to the choice of this parameter, use IBRION=2, if you are not sure how large the optimal POTIM is.

这里推荐在结构优化中使用IBRION=2和POTIM=0.5.
在侯柱峰编写的《vasp软件包使用入门指南》中《§6.2结构参数的优化》中“一次成型”一节写道：
引用回帖:
SYSTEM = Mg-hex
ENCUT = 250
ISTART = 0; ICHARG = 2
ISMEAR = 1; SIGMA = 0.2
NSW = 60; IBRION = 2
ISIF = 3
POTIM = 0.2
EDIFF = 1E-6; EDIFFG = -1E-3
PREC = Accurate

那么应该如何选择POTIM的值呢？
回复此楼

» 收录本帖的淘帖专辑推荐

VASP

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-12-16 23:22:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★
sunyang1988(金币+2):谢谢帮助 2010-12-18 12:15:15
我一般使用默认值POTIM=0.5。而且IBRION=2的算法比较可信。
驰豫的时候POTIM的设置只会影响运行的效率，不会对结果又明显的影响。
不用管它，默认值很不错。

[ Last edited by stractor on 2010-12-18 at 10:53 ]
赞一下(1人) 回复此楼
2楼2010-12-18 10:52:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)
引用回帖:
Originally posted by stractor at 2010-12-18 10:52:01:
我一般使用默认值POTIM=0.5。而且IBRION=2的算法比较可信。
驰豫的时候POTIM的设置只会影响运行的效率，不会对结果又明显的影响。
不用管它，默认值很不错。

[ Last edited by stractor on 2010-12-18 at 1 ...

结构驰豫时，POTIM既然只是一个没有物理含义的缩放系数，那么肯定对于结果没有影响。
不过，POTIM与计算效率的影响是怎么样的呢？
赞一下 回复此楼
3楼2010-12-19 08:42:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangguiyang

新虫 (小有名气)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指教，期待你更好的表现 2010-12-21 08:19:50
对于IBRION=2，POTIM不同值对效率影响不大的，就用0.5。如果对于初始结构严重偏离基态的，应先用较小的POTIM,例如0.1,0.05，甚至0.01。力减少到一定程度后，再用0.5.
赞一下(2人) 回复此楼
4楼2010-12-21 06:03:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swx0789

铜虫 (小有名气)

LZ,请教你个问题

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
youzhizhe:鼓励交流。 2010-12-25 15:31:03
youzhizhe(金币+1):忘记金币了 ，不好意思。鼓励交流。 2010-12-25 15:31:21
我现在也在优化晶体的结构参数。看了后博士的中文入门手册，虽然讲的很详细，当有个地方还不是很明白。就他例子的关于Mg的晶格参数的优化（P25）。
我取10个点，循环完后，拟合得到最佳的a值。他后面说“在该体积下，重复上面第2和3步，得到平衡状态下的其他晶胞参数。这一步就是：在得到了E(V)曲线后，通过拟合状态方程得到平衡状态下的体积，计算出上面脚本中的变量$i的值，并改变$i的循环值，再运行run_cell计算一次，得到其他的结构参数c和位置u.。”
对这句话我不太理解，安我自己的意思，通过拟合上面E（V）曲线得到平衡状态下的体积，计算出变量$i,然后在该值下迟豫，就可以得到b，c。而不是候博士上面说的“改变$i的循环值，再运行run_cell计算一次”。
不知道大家怎么理解他这句话的。
赞一下(3人) 回复此楼
5楼2010-12-23 08:24:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zhang668(金币+3):多谢指教！ 2010-12-23 18:51:30
引用回帖:
Originally posted by swx0789 at 2010-12-23 08:24:33:
我现在也在优化晶体的结构参数。看了后博士的中文入门手册，虽然讲的很详细，当有个地方还不是很明白。就他例子的关于Mg的晶格参数的优化（P25）。
我取10个点，循环完后，拟合得到最佳的a值。他后面说“在该体积 ...

“改变$i的循环值，再运行run_cell计算一次”
此时$i只附了一个值，就是前面由EOS拟合得到的平衡体积V_0计算(V_0 = 0.5*sqrt(3.0) a^3 * c/a,  c/a的值取自run_cell所设置的，即c/a的值与前一步中run_cell中的相同。这样得到了一个平衡体积下初始的a, c/a)。这一次run_cell就是固定体积(将体积固定在所拟合出来的V_0)，根据run_cell所设置的初始a, c/a进行固定体积的优化。

E(a,c) -----> E(V, a)------->E(V_0,a)------> V_0下最优的a, c的值。
赞一下(3人) 回复此楼
6楼2010-12-23 18:47:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swx0789

铜虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by valenhou001 at 2010-12-23 18:47:48:


“改变$i的循环值，再运行run_cell计算一次”
此时$i只附了一个值，就是前面由EOS拟合得到的平衡体积V_0计算(V_0 = 0.5*sqrt(3.0) a^3 * c/a,  c/a的值取自run_cell所设置的，即c/a的值与前一步中run_cell中 ...

谢谢斑竹的回答，我理解了这个意思。当我还有点疑问，vasp 好像是在开始十步优化的时候，都是固定角度和c/a的值进行优化的。万一我优化的是一个三斜晶系，我该怎么办。我还是安你的这个步骤：E（a，b，c)->E(v,a)->E（V_0,a)->V_0下的最佳a，b，c，alpha，beta，gama吗？不知道这个有通用步骤没有？
因为我还想优化结构：a！＝b！＝c，三个角度都为90度。如果还是按Mg那种方法的话，b/a的值将会始终不变。

[ Last edited by swx0789 on 2010-12-25 at 15:22 ]
赞一下(1人) 回复此楼
7楼2010-12-25 14:56:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2008130252

银虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
6楼: Originally posted by valenhou001 at 2010-12-23 18:47:48
“改变$i的循环值，再运行run_cell计算一次”
此时$i只附了一个值，就是前面由EOS拟合得到的平衡体积V_0计算(V_0 = 0.5*sqrt(3.0) a^3 * c/a,  c/a的值取自run_cell所设置的，即c/a的值与前一步中run_cell中的 ...

侯老师您好，看来了您的讲解我还是不大明白。这一次是改变哪个变量进行计算呢？也就是 i 代表了哪个值
赞一下 回复此楼
8楼2014-12-11 18:37:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Prettyswan 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600 286分 材料求调剂 +10 麻辣鱿鱼 2026-03-27 11/550 2026-03-29 21:52 by nanaliuyun
[考研] 277跪求调剂 +6 1915668 2026-03-27 10/500 2026-03-29 16:03 by 王亮_大连医科大
[考研] 一志愿北京理工大学本科211材料工程294求调剂 +8 mikasa的围巾 2026-03-28 8/400 2026-03-29 12:48 by 无际的草原
[考研] 295求调剂 +4 wei-5 2026-03-26 4/200 2026-03-28 23:20 by 小木虫tim
[考研] 学硕274求调剂 +9 Li李鱼 2026-03-26 9/450 2026-03-28 21:42 by bymhappy
[考研] 压国家一区线，求导师收留，有恩必谢！ +7 迷人的哈哈 2026-03-28 7/350 2026-03-28 16:47 by 催化大白
[考研] 复试调剂 +3 raojunqi0129 2026-03-28 3/150 2026-03-28 15:27 by 落睿可思
[考研] 0703本科郑州大学求调剂 +3 nhj_ 2026-03-25 3/150 2026-03-28 13:24 by Iveryant
[考研] 286求调剂 +12 PolarBear11 2026-03-26 12/600 2026-03-28 12:14 by zllcz
[考研] 0856，材料与化工321分求调剂 +12 大馋小子 2026-03-27 13/650 2026-03-28 10:56 by self2008
[考研] 0703化学求调剂，各位老师看看我！！！ +5 祁祺祺 2026-03-25 5/250 2026-03-27 21:44 by 东方猪猪
[考研] 考研调剂 +4 Sanmu-124 2026-03-26 4/200 2026-03-27 17:49 by kiokin
[考研] 085600材料与化工调剂 +10 A-哆啦Z梦 2026-03-23 16/800 2026-03-27 15:13 by caszguilin
[考研] 0856调剂 +5 求求让我有书读� 2026-03-26 6/300 2026-03-27 15:12 by caszguilin
[考研] 化学308分求调剂 +8 你好明天你好 2026-03-23 9/450 2026-03-27 14:01 by 杨光于青云
[考研] 325求调剂 +3 Aoyijiang 2026-03-23 3/150 2026-03-26 20:46 by 不吃魚的貓
[考研] 生物学 296 求调剂 +4 朵朵- 2026-03-26 6/300 2026-03-26 19:01 by 不吃魚的貓
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-25 3/150 2026-03-25 19:03 by Ainin_
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
信息提示
关闭
请填处理意见
关闭