24小时热门版块排行榜    

查看: 1229  |  回复: 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的回帖

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的回帖
查看全部 3 个回答

毛毛虫_

银虫 (小有名气)

楼主,这个问题解决了吗?我现在也出现了这个问题
努力爬行的小虫子·····
2楼2019-04-19 09:46:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 +5 晴天加油 2026-08-26 6/300 2026-08-27 04:49 by littlezl
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +3 a089 2026-08-26 3/150 2026-08-26 23:21 by young90
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +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
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见