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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 结构优化与静态计算的能量相差1.4eV，怎么可能呢？
341/1返回列表
查看: 3552  |  回复: 33
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lijuanshu

铜虫 (正式写手)

[交流] 结构优化与静态计算的能量相差1.4eV，怎么可能呢？

我用VASP优化结构得到的能量是-306.54153eV，把优化后的CONTCAR拷成POSCAR做静态计算得到的能量却是-307.95863eV，它们的能量相差了1.4eV，这也太不可思议了吧，还请大家分析下原因啊，不胜感激。
回复此楼

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

第一性原理

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2013-03-14 21:38:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Lathander

新虫 (小有名气)
★
lijuanshu(金币+1): 谢谢参与
体积变化了，basis set不一样了，wavefunction不一样了，等等
赞一下 回复此楼
2楼2013-03-15 07:07:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
2楼: Originally posted by Lathander at 2013-03-15 07:07:18
体积变化了，basis set不一样了，wavefunction不一样了，等等

静态计算时，原子的坐标只是发生很小的移动，我感觉都可以忽略不计。怎样才能使结构优化与静态计算的能量一致呢？
赞一下 回复此楼
3楼2013-03-15 08:43:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
我用杂质元素A参杂B，A有两个间隙位置，A在间隙位置1的结构优化的能量是-306.54153eV，静态计算的能量是-307.95863eV；A在间隙位置2的结构优化的能量是-307.92116eV，静态计算的能量是-305.71585eV，我现在都搞不懂到底哪个位置更稳定了，还请大家帮我分析下，如何解决这种问题。
赞一下 回复此楼
4楼2013-03-15 09:01:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)
★ ★
lijuanshu(金币+1): 谢谢参与
ben_ladeng: 金币+1, 谢谢 2013-03-16 14:24:06
当然是
优化后的CONTCAR拷成POSCAR做静态计算得到的能量.
应该是位置1稳定。静态计算的能量是-307.95863eV。两个计算的参数设置要一样···
赞一下 回复此楼
5楼2013-03-15 09:13:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
5楼: Originally posted by Illusionist at 2013-03-15 09:13:09
当然是
优化后的CONTCAR拷成POSCAR做静态计算得到的能量.
应该是位置1稳定。静态计算的能量是-307.95863eV。两个计算的参数设置要一样···

参数的设置都是一样的， 但是我算的有些结构的结构优化和静态计算的能量就是几乎相同的，不知道为什么？以下是我结构优化的INCAR：
SYSTEM =
ENCUT =
ISTART = 0
ICHARG = 2
ISMEAR = 1
SIGMA = 0.2
NSW = 2000
IBRION = 2
ISPIN = 2
NELM = 2000
POTIM = 0.2
EDIFF = 1E-4
IALGO = 48
LREAL = Auto
GGA = 91
VOSKOWN = 1
下面是我静态计算的INCAR：
SYSTEM =
ENCUT =
ISTART = 0
ICHARG = 2
ISMEAR = 1
SIMGA = 0.2
ISPIN = 2
NELM = 2000
EDIFF = 1E-4
IALGO = 48
LREAL = Auto
GGA = 91
VOSKOWN = 1
不知道在设置上是否有些问题，还请指出，不胜感激。
赞一下 回复此楼
6楼2013-03-15 09:54:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)
ENCUT最好自己设定。ISMEAR 取值。
你设定这么多，最好了解每一项的意义。

最好先看看侯老师的一些资料·照着做一做··


个人看法，不一定正确··还请包涵·
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : VASP的个人经验手册-侯柱峰2004版.pdf
    • 2013-03-15 14:00:57, 245.92 K
7楼2013-03-15 14:01:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
7楼: Originally posted by Illusionist at 2013-03-15 14:01:26
ENCUT最好自己设定。ISMEAR 取值。
你设定这么多，最好了解每一项的意义。

最好先看看侯老师的一些资料·照着做一做··


个人看法，不一定正确··还请包涵·

这个手册我是看过的，我觉得这些设置应该是没问题啊，但是能量怎么差这么多啊
赞一下 回复此楼
8楼2013-03-15 16:20:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongjijun

木虫 (正式写手)
★
lijuanshu(金币+1): 谢谢参与
增加K点，提高ecut试试。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
9楼2013-03-16 13:58:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongjijun

木虫 (正式写手)
增加K点，提高ecut试试。

[ 发自手机版 http://muchong.com/3g ]

发重了，抱歉!

[ Last edited by gongjijun on 2013-3-16 at 16:28 ]
赞一下 回复此楼
10楼2013-03-16 14:48:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
10楼: Originally posted by gongjijun at 2013-03-16 14:48:23
增加K点，提高ecut试试。



发重了，抱歉!

但是我想所有的计算都尽量保持参数的一致，像K点和ENCUT我都是保持一致的，前面一些的计算结果符合的很好，只是这两个出现了问题
赞一下 回复此楼
11楼2013-03-17 09:27:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongjijun

木虫 (正式写手)
引用回帖:
11楼: Originally posted by lijuanshu at 2013-03-17 09:27:36
但是我想所有的计算都尽量保持参数的一致，像K点和ENCUT我都是保持一致的，前面一些的计算结果符合的很好，只是这两个出现了问题...

把这个优化好的结构作为初始结构再优化一次看看。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
12楼2013-03-17 17:29:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
12楼: Originally posted by gongjijun at 2013-03-17 17:29:34
把这个优化好的结构作为初始结构再优化一次看看。
...

好的，我试试这种方法，谢谢啊
赞一下 回复此楼
13楼2013-03-18 07:59:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
★
lijuanshu(金币+1): 谢谢参与
在优化和计算之前要做收敛测试
赞一下 回复此楼
14楼2013-03-18 11:39:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

随便注册

银虫 (初入文坛)
★
lijuanshu(金币+1): 谢谢参与
期待后续，楼主找出问题来了没？
赞一下 回复此楼
15楼2013-03-18 14:47:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
14楼: Originally posted by liliangfang at 2013-03-18 11:39:09
在优化和计算之前要做收敛测试

如何测试，还希望详细说明啊
赞一下 回复此楼
16楼2013-03-18 15:07:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

npuouwen

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
15楼: Originally posted by 随便注册 at 2013-03-18 14:47:26
期待后续，楼主找出问题来了没？

求问简单二元合金系的能量可能达到5万ev不呀。。。
赞一下 回复此楼
17楼2013-03-18 15:37:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
16楼: Originally posted by lijuanshu at 2013-03-18 15:07:02
如何测试，还希望详细说明啊...

候柱峰有个中文的教程，你可以搜索一下，里面介绍如何测试ENCUT，SIGMA，k等
赞一下 回复此楼
18楼2013-03-18 20:49:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
17楼: Originally posted by npuouwen at 2013-03-18 15:37:01
求问简单二元合金系的能量可能达到5万ev不呀。。。...

这么大啊，我的都只是几百个eV，我也不是很清楚，可能不同体系的能量不同吧
赞一下 回复此楼
19楼2013-03-19 08:04:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
12楼: Originally posted by gongjijun at 2013-03-17 17:29:34
把这个优化好的结构作为初始结构再优化一次看看。
...

我再次优化出现的问题就更严重了，两个体系的情况分别如下：
A：第一次：-307.92116 eV； 第二次：-304.02535 eV
B：第一次：-306.54153 eV； 第二次：-307.96147 eV
我看到这些数据都要崩溃了，这都什么情况啊
赞一下 回复此楼
21楼2013-03-19 08:22:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
15楼: Originally posted by 随便注册 at 2013-03-18 14:47:26
期待后续，楼主找出问题来了没？

仍在探索
回复此楼
22楼2013-03-19 08:22:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你做的是什么？体材料？金属半导体绝缘体？~
赞一下 回复此楼
23楼2013-03-19 13:36:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
23楼: Originally posted by Illusionist at 2013-03-19 13:36:58
你做的是什么？体材料？金属半导体绝缘体？~

研究的元素在金属间化合物晶界处的偏聚情况
赞一下 回复此楼
24楼2013-03-19 14:02:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
24楼: Originally posted by lijuanshu at 2013-03-19 14:02:54
研究的元素在金属间化合物晶界处的偏聚情况...

• For semiconductors or insulators use the tetrahedron method (ISMEAR=-5), if the cell is too large (or if you use only a
single or two k-points) use ISMEAR=0 in combination with a small SIGMA=0.05.
• For relaxations in metals always use ISMEAR=1 or ISMEAR=2 and an appropriate SIGMA value (the entropy term should
be less than 1 meV per atom). Mind: Avoid to use ISMEAR>0 for semiconductors and insulators, since it might cause
problems.
For metals a sensible value is usually SIGMA= 0.2 (which is the default).
• For the calculations of the DOS and very accurate total energy calculations (no relaxation in metals) use the tetrahedron
method (ISMEAR=-5).


对于半导体而已，我遇到过有时候ismear=1和-5计算能量相差较大。·····
赞一下 回复此楼
25楼2013-03-19 14:48:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

第七敌法师

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
应该是你设置的问题，你试试把那些不必要的参数去掉，系统提醒你出问题时在加入里面。可能我们算得方向不一样。我感觉，你结构优化时的INCAR文件里，有的参数根本就是多余的。
赞一下 回复此楼
26楼2013-03-19 15:24:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
25楼: Originally posted by Illusionist at 2013-03-19 14:48:36
• For semiconductors or insulators use the tetrahedron method (ISMEAR=-5), if the cell is too large (or if you use only a
single or two k-points) use ISMEAR=0 in combination with a small SIGM ...

我计算的时候ISMEAR都是取的是1啊
赞一下 回复此楼
27楼2013-03-19 15:45:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
26楼: Originally posted by 第七敌法师 at 2013-03-19 15:24:40
应该是你设置的问题，你试试把那些不必要的参数去掉，系统提醒你出问题时在加入里面。可能我们算得方向不一样。我感觉，你结构优化时的INCAR文件里，有的参数根本就是多余的。

我个人认为这些参数设置应该没有问题吧，这些参数我也都查了VASP的手册
赞一下 回复此楼
28楼2013-03-19 15:46:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

沐浴晨光_7月

银虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主你ISPIN=2说明你考虑了自旋，如果这样的话，再次做静态可以尝试把CHGCAR和WAVECAR都加上，这样一来，可以快点，也可以计算的更准确。
赞一下 回复此楼
29楼2013-03-22 10:36:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
29楼: Originally posted by 沐浴晨光_7月 at 2013-03-22 10:36:35
楼主你ISPIN=2说明你考虑了自旋，如果这样的话，再次做静态可以尝试把CHGCAR和WAVECAR都加上，这样一来，可以快点，也可以计算的更准确。

“再次做静态可以尝试把CHGCAR和WAVECAR都加上“，你的意思是把结构优化的CHGCAR和WAVECAR都拷到静态计算的文件夹下面吗？但是我在静态计算的时候设置ISTART = 0，ICHARG =2，就是拷过来也没用吧
赞一下 回复此楼
30楼2013-03-22 12:47:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

沐浴晨光_7月

银虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
30楼: Originally posted by lijuanshu at 2013-03-22 12:47:14
“再次做静态可以尝试把CHGCAR和WAVECAR都加上“，你的意思是把结构优化的CHGCAR和WAVECAR都拷到静态计算的文件夹下面吗？但是我在静态计算的时候设置ISTART = 0，ICHARG =2，就是拷过来也没用吧...

嗯，但可以把ISTART注销，采用默认值，拷贝过来就有用了。你可以看看手册里的ISTART说明。
赞一下 回复此楼
31楼2013-03-22 12:56:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
31楼: Originally posted by 沐浴晨光_7月 at 2013-03-22 12:56:01
嗯，但可以把ISTART注销，采用默认值，拷贝过来就有用了。你可以看看手册里的ISTART说明。...

当WAVECAR存在时，ISTART = 1，那ICHARG也是取默认值吗
赞一下 回复此楼
32楼2013-03-22 13:06:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

沐浴晨光_7月

银虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
32楼: Originally posted by lijuanshu at 2013-03-22 13:06:43
当WAVECAR存在时，ISTART = 1，那ICHARG也是取默认值吗...

这些vasp都有啊，你要自己仔细看看咯
赞一下 回复此楼
33楼2013-03-22 15:38:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宇ll

银虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
请问你的问题最后解决了吗？？
赞一下 回复此楼
34楼2018-08-07 12:16:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
lijuanshu20楼
2013-03-19 08:08   回复  
引用回帖:
18楼: Originally posted by liliangfang at 2013-03-18 20:49:18 候柱峰有个中文的教程，你可以搜索一下，里面介绍如何测试ENCUT，SIGMA，k等...

相关版块跳转 我要订阅楼主 lijuanshu 的主题更新
341/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭