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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】利用侯博的方法优化结构出现的问题（求验证）
61/1返回列表
查看: 2686  |  回复: 11
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

qiqi2926

铜虫 (小有名气)

[交流] 【求助】利用侯博的方法优化结构出现的问题（求验证） 已有5人参与

利用侯博的文章中优化的部分进行优化，照例对ENCUT，KPT，和SIGMA进行收敛测试，得到ENCUT=650eV，K=4*4*4，SIGMA=0.1。然后采用结构参数优化的复杂情况，因为体系为正交。
其中INCAR.relax：
SYSTEM = NH3BH3-pmn21
ENCUT = 650
ISTART = 0; ICHARG = 2
ISMEAR = 1; SIGMA = 0.1
NSW = 60; IBRION = 2
ISIF = 5
POTIM = 0.2
EDIFF = 1E-5; EDIFFG = -1E-3
PREC = Accurate

INCAR.static：
SYSTEM = NH3BH3-pmn21
ENCUT = 650
ISTART = 0; ICHARG = 2
ISMEAR = -5
PREC = Accurate


KPT：
Auto generation
0
Gamma
8 8 8
0.0 0.0 0.0

POTCAR为：
paw_GGA

但是计算得到的结果居然为：

5.041 98.6200 -62.534945
5.141 104.6000 -65.167658
5.241 110.8300 -67.385080
5.341 117.2900 -69.230747
5.441 124.0100 -70.744046
5.541 130.9700 -71.959857
5.641 138.1900 -72.909812
5.741 145.6700 -73.622417
5.841 153.4200 -74.123520
。。。
与理论期望不符。。。
望高手指点
回复此楼

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

VASP

» 猜你喜欢

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

valenhou001

至尊木虫 (职业作家)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢侯老师指点 2010-07-26 18:23:09
问题是你采用了ISIF=5的原因，原子的位置没有进行优化。 建议采用ISIF=4，你的体系是绝缘体吧？那在结构优化的过程中也可以采用ISMEAR =0 (Gaussian broadening), 或ISMEAR=-5。

[ Last edited by valenhou001 on 2010-7-26 at 18:13 ]
赞一下(2人) 回复此楼
2楼2010-07-26 18:07:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢侯老师指点 2010-07-27 09:25:57
没有必要测试不同ENCUT下的晶格常数。采用实验的晶格常数，计算不同ENCUT下，总能的收敛情况就足够了。不同ENCUT下，得到的晶格常数是略有变化的，是比较正常的。
ENCUT对Cu,Ｃr, S取到800~900似乎有点太高了。

研究此类材料的重点应该是它的磁性：磁序的能量高低等。

[ Last edited by valenhou001 on 2010-7-27 at 09:15 ]
赞一下(2人) 回复此楼
4楼2010-07-27 09:11:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢 2010-07-28 10:31:57
引用回帖:
Originally posted by identation at 2010-07-27 23:46:47:


看了您分享的那篇怎样做有意义的第一原理计算文章了,里面提到,不同ECUT下力的收敛难于总能的收敛,那如果进行几何优化,需不需要测试力的收敛呢?谢谢老师

既然你已经知道了一点原因，一般不需要对力进行ENCUT的收敛性测试。
赞一下(2人) 回复此楼
8楼2010-07-28 09:54:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢 2010-07-28 10:32:05
引用回帖:
Originally posted by nkleof at 2010-07-27 15:10:33:

谢谢侯老师指点，但是存在一个问题，就是如何确定晶格常数值作为后续计算的参数？由于不同的ENCUT导致的晶格常数相差还是不小的，可能所计算出来的性质不一定都一样吧？
我看见上次那篇 做有意义的DFT计算 文献 ...

主要的还是电子结构等的计算结果。晶格常数在实验值左右5%以内是可以完全接受的。
赞一下(2人) 回复此楼
9楼2010-07-28 09:57:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢侯老师知道 2010-07-29 20:40:33
结构优化中原子受力收敛和测试总能的变化对ENCUT的收敛是两码事。 要先保证ENCUT足够高，以能保证有足够的计算精度，然后再来做结构优化的。 力是通过总能来算出来的，原则上，测试力的变化对ENCUT的收敛难度要比测试总能对ENCUT的更严格。测试总能的变化对ENCUT的收敛，就能给力的收敛提供一个保证了。

[ Last edited by valenhou001 on 2010-7-29 at 16:02 ]
赞一下(2人) 回复此楼
12楼2010-07-29 15:59:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qiqi2926 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭