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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】利用侯博的方法优化结构出现的问题（求验证）
121/2返回列表12下一页
查看: 2872  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

ApolloYang

木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢提示 2010-07-27 16:26:45
引用侯博士的话
“对六角晶格, E是 a,c的函数，根据体积V与a,c的关系，可以得到E也是a, V的函数，或者E是V, c的函数或 E是V, c/a的函数。明白了这个几个关系之后就能大致明白那个例子里面所展示的内容。考虑的是计算不同V下面的总能，对给定一个V，初始化一个a, c/a，保持体积不变去优化优化a, c/a。计算了一组V下的总能，就有了E~V的数据，然后采用状态方程拟合，得到平衡体积V0，然后再初始化该体积下的a, c/a再进行固定体积的优化得到平衡态下的a0, c0/a0。”
赞一下(3人) 回复此楼
6楼2010-07-27 15:13:19
已阅   回复此楼   关注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的回帖

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

nkleof

银虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by valenhou001 at 2010-07-27 09:11:41:
没有必要测试不同ENCUT下的晶格常数。采用实验的晶格常数，计算不同ENCUT下，总能的收敛情况就足够了。不同ENCUT下，得到的晶格常数是略有变化的，是比较正常的。
ENCUT对Cu,Ｃr, S取到800~900似乎有点太高了。
...

谢谢侯老师指点，但是存在一个问题，就是如何确定晶格常数值作为后续计算的参数？由于不同的ENCUT导致的晶格常数相差还是不小的，可能所计算出来的性质不一定都一样吧？
我看见上次那篇 做有意义的DFT计算 文献里面，晶格常数也是会收敛的吧？为什么在这里会出现ENCUT 增大后，晶格常数持续增大的情况？
赞一下(1人) 回复此楼
5楼2010-07-27 15:10:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by valenhou001 at 2010-07-27 09:11:41:
没有必要测试不同ENCUT下的晶格常数。采用实验的晶格常数，计算不同ENCUT下，总能的收敛情况就足够了。不同ENCUT下，得到的晶格常数是略有变化的，是比较正常的。
ENCUT对Cu,Ｃr, S取到800~900似乎有点太高了。
...

看了您分享的那篇怎样做有意义的第一原理计算文章了,里面提到,不同ECUT下力的收敛难于总能的收敛,那如果进行几何优化,需不需要测试力的收敛呢?谢谢老师
赞一下(1人) 回复此楼
7楼2010-07-27 23:46:47
已阅   回复此楼   关注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的回帖

qiqi2926

铜虫 (小有名气)
引用回帖:
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 ...

谢谢侯老师。。。
回复此楼
10楼2010-07-28 11:23:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qiqi2926 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北理工298英一数二已上岸，感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 求调剂 +9 翩翩一书生 2026-04-09 9/450 2026-04-10 11:43 by wp06
[考研] 材料工程302分求调剂 +18 zyx上岸！ 2026-04-04 18/900 2026-04-10 10:07 by 314126402
[考研] 青岛科技大学材料学院，环境学院调剂补录4月10日以前都可以 +3 1青科大。 2026-04-09 5/250 2026-04-10 09:58 by 翩翩一书生
[考研] 求调剂 +15 张zic 2026-04-05 16/800 2026-04-10 08:12 by kangsm
[考研] 生物学求调剂 一志愿沪9，326分 +7 刘墨墨 2026-04-06 7/350 2026-04-10 08:11 by kangsm
[考研] 289求调剂 +4 L1ttleTiger 2026-04-04 4/200 2026-04-09 19:53 by xiayizhi
[考研] 一志愿中国科学院上海有机所，有机化学356分找调剂 +11 Nadiums 2026-04-09 11/550 2026-04-09 18:04 by lijunpoly
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 复试调剂，一志愿郑州大学材料与化工289分 +31 硕星赴 2026-04-08 31/1550 2026-04-09 16:54 by Delta2012
[考研] 1U盾记得记得就 +9 sanjin020722 2026-04-08 10/500 2026-04-09 14:11 by 诗与自由
[考研] 求调剂，262机械专硕 +6 嗯yyl 2026-04-08 6/300 2026-04-09 12:01 by zhouyuwinner
[考研] 求调剂 +3 猪肉墩粉条cc 2026-04-08 4/200 2026-04-09 10:05 by 猪肉墩粉条cc
[考研] 求助071001调剂！！！ +7 黄守松 2026-04-05 8/400 2026-04-09 09:07 by 徐良白眉大侠
[考研] 化学308分求调剂 +21 你好明天你好 2026-04-07 23/1150 2026-04-08 22:32 by 凯凯要变帅
[考研] 求调剂 +9 月@163.com 2026-04-07 11/550 2026-04-08 14:48 by qlm5820
[考研] 生物学学硕求调剂：351分一志愿南京师范大学生物学专业 +6 …～、王…～ 2026-04-06 7/350 2026-04-06 18:54 by macy2011
[考研] 考研调剂 +5 美丽的youth_ 2026-04-04 6/300 2026-04-06 06:57 by houyaoxu
[考研] 309求调剂 +4 快乐的小白鸽 2026-04-04 5/250 2026-04-04 15:55 by cql1109
信息提示
关闭
请填处理意见
关闭