24小时热门版块排行榜    

查看: 1226  |  回复: 2

fangyongxinxi

新虫 (正式写手)

[求助] 为什么vasp能一直向高能量方向优化? 已有1人参与

为什么vasp能一直向高能量方向优化?
自己打算对一个体系以同样参数,进行结构多次优化,ISIF=3,其中晶格常数abc如下,后面的数字是经过优化得到的,第一列是初始数据(别人文献中查找的),第二列是第二次优化,第三列是第三次优化结果,依次类推。能量:第一次优化的能量 -55.64ev, 第二次,-53.53ev ,第三次,-52.22ev,第四次,-51.44ev。能量随着优化次数逐渐升高。

a 5.75 5.62 5.56 5.52 5.51
b 4.06 3.97 3.92 3.89 3.88
c 4.06 3.97 3.92 3.89 3.88


求高手解答:
1. 为什么会出现优化向高能量方向发生这样的情况?(一般情况是如果一个体系达到能量和受力标准,就会停止优化,比如MS,为什么vaso能一直优化?一种解释是这个体系有亚稳态存在,但能量趋势如何解释?自己将尝试用同样的方法优化一些简单体系,比如cu ,si 等,看是否有同样结果)
2. 就这个体系来说,自己应该选哪个作为性质计算的初始结构?第一次优化的,还是后面的,因为第一次优化的能量最小。

后面自己又进行了尝试:1. 首先设定ISIF=2,优化,然后再用isif=3优化,得到结果与上面直接用ISIF=3结果一样;2. 首先用ISIF=4,再用ISIF=3,得到结果和上面直接用ISIF=3一样。所以后面得出结论:对于本体系,直接用ISIF=3优化一次就行了。这样做对吗?

谢谢。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

毛毛虫_

银虫 (小有名气)

楼主,这个问题解决了吗?我现在也出现了这个问题
努力爬行的小虫子·····
2楼2019-04-19 09:46:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

1. 为什么会出现优化向高能量方向发生这样的情况?(一般情况是如果一个体系达到能量和受力标准,就会停止优化,比如MS,为什么vaso能一直优化?一种解释是这个体系有亚稳态存在,但能量趋势如何解释?自己将尝试用同样的方法优化一些简单体系,比如cu ,si 等,看是否有同样结果)

> VASP 就如你想的,就像一般的軟體 (比如 MS) 一樣,結構優化是會往能量低的方向前進
ISIF=3 通常在 1~2 次就會找到最低能量,之後再怎麼重覆做,結果都會是一樣,會發生像你遇到的情況,我能想到以下幾點

(1)  你的計算沒有 "收斂"
VASP 裡面有個 electron loop convergence criterion,叫 NELM。default value 是 60
如果你electron loop 到達這次數,VASP 會?娭仆V梗瑢е掠嬎悴皇諗浚浣Y果是不可信的
(2)  你的 VASP 沒有正確編輯
有可能你使用的 math kernel library 有 bug,或是編輯上有參數沒有給正確,導致計算結果是不可信
(3) 你 INCAR 內可能有加了些特殊的參數,導致計算不收斂,或是收斂到不正確的地方

在沒有看到 INCAR, POSCAR, POTCAR, KPOINTS 前,很難判斷發生什麼事
但能肯定的是,正常的 VASP 計算,是不會往能量高的地方去收斂

2. 就这个体系来说,自己应该选哪个作为性质计算的初始结构?第一次优化的,还是后面的,因为第一次优化的能量最小。
> 要先確認計算是"正確"
3楼2019-04-23 14:08:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fangyongxinxi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +14 爱看书的可乐 2026-08-26 16/800 2026-08-26 19:21 by skyjudy
[基金申请] 我不理解! +9 Edward_pc 2026-08-26 14/700 2026-08-26 18:03 by hsding2026
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见