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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】利用侯博的方法优化结构出现的问题（求验证）
51/1返回列表
查看: 2886  |  回复: 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(金币+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的回帖
查看全部 12 个回答

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的回帖

nkleof

银虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by valenhou001 at 2010-07-26 18:07:30:
问题是你采用了ISIF=5的原因，原子的位置没有进行优化。 建议采用ISIF=4，你的体系是绝缘体吧？那在结构优化的过程中也可以采用ISMEAR =0 (Gaussian broadening), 或ISMEAR=-5。

[ Last edited by valenhou00 ...

跟帖请问一下，在优化时遇到所得晶格常数随ENCUT值增大而增大的情况（从ENCUT=300到800eV，晶格常数从9.836涨到9.852，而实验值在9.814），问题可能出在哪里？采用的优化方法如下：
INCAR --》
SYSTEM = CuCr2S4 cubic cell spinel

PREC = High
ENCUT = 900
LREAL = A
IALGO = 48
NPAR = 8
NSIM = 8
ISTART = 0
ICHARG = 2

NELMIN = 4
LORBIT = 12
EDIFF = 1E-5
ISMEAR = -5
LWAVE = .False.

NSW = 60
IBRION = 2
ISIF = 3

ISPIN = 2
MAGMOM = 2*0 4*4 8*0

KPOINTS --》
Auto
0
M
13  13  13
0   0   0

k点网格做过收敛测试，从此往后更密的网格造成的总能变化在0.00015eV/cell，每原胞14个原子。

体系是准金属性的，但是在费米面处的态密度值很小，且是磁性体系。是否有可能是ISMEAR=-5造成的问题？
赞一下(1人) 回复此楼
3楼2010-07-26 19:20:38
已阅   回复此楼   关注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的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 药学求调剂 +12 喽哈加油 2026-04-14 14/700 2026-04-15 23:49 by 陈皮皮
[考研] 急需调剂 +5 绝不放弃22 2026-04-15 5/250 2026-04-15 23:11 by Tide man
[考研] 271求调剂 +32 2261744733 2026-04-11 33/1650 2026-04-15 22:03 by noqvsozv
[考研] 297，工科调剂? +10 河南农业大学-能 2026-04-14 10/500 2026-04-15 21:50 by noqvsozv
[考研] 化学070300 求调剂 +23 哈哈哈^_^ 2026-04-12 23/1150 2026-04-14 16:30 by zhouxiaoyu
[考研] 农学0904 312求调剂 +4 Say Never 2026-04-11 4/200 2026-04-14 09:10 by zs92450
[考研] 求调剂，985材料与化工348分 +9 涵竹刘 2026-04-11 14/700 2026-04-13 22:26 by 涵竹刘
[考研] 生物学调剂 +11 小冉要努力 2026-04-10 13/650 2026-04-13 11:46 by 电化学及催化
[考研] 一志愿085802 323分求调剂 +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[考研] 2本，初试303，0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 一志愿厦大0856，306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 085400 328分 求调剂 +10 喂你一个大橙子 2026-04-09 14/700 2026-04-11 19:53 by lqspecial
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 一志愿东北大学控制工程085406数二英二385，求调剂 +8 Ezra_Zhang 2026-04-09 8/400 2026-04-11 09:15 by 猪会飞
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 344求调剂 +7 丶风雪夜归人丶 2026-04-09 7/350 2026-04-10 12:05 by pengliang8036
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 本科西工大 0856 324求调剂 +10 wysyjs25 2026-04-09 11/550 2026-04-10 08:37 by 5268321
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 332，085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
信息提示
关闭
请填处理意见
关闭