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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】利用侯博的方法优化结构出现的问题（求验证）
121/2返回列表12下一页
查看: 2874  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 化学工程与技术324调剂 +23 孙常华 2026-04-09 25/1250 2026-04-11 00:07 by 骑牛渡寒江
[考研] 材料类284调剂 +40 想换手机不想解� 2026-04-08 48/2400 2026-04-10 23:28 by 314126402
[考研] 吉大计算机技术331分，英语六级，求调剂 +3 峰峰021116 2026-04-09 3/150 2026-04-10 20:01 by chemisry
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 环境专硕调剂 +16 会说话的肘子 2026-04-06 16/800 2026-04-10 10:30 by asy1wn
[考研] 一志愿电子科技大学085600材料与化工 329分求调剂 +14 Naiko 2026-04-04 14/700 2026-04-09 16:56 by luoyongfeng
[考研] 一志愿武理车辆 281 求调剂 +5 上岸研究生. 2026-04-07 5/250 2026-04-09 15:56 by only周
[硕博家园] 有没有学校材料专业收跨调(一志愿085410) +5 momo(上岸版) 2026-04-06 8/400 2026-04-09 15:07 by only周
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 286求调剂 +20 Faune 2026-04-06 20/1000 2026-04-07 11:33 by 诗与自由
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +6 白易辰 2026-04-06 7/350 2026-04-06 22:46 by 875465
[考研] 085405软件工程301分求调剂，专硕可跨专业，四六级已过 +3 静静想想 2026-04-05 3/150 2026-04-06 15:23 by nepu_uu
[考研] 求调剂 +10 Hll胡 2026-04-04 10/500 2026-04-05 20:09 by nepu_uu
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 考研调剂 +3 mcbbc 2026-04-04 3/150 2026-04-05 10:03 by barlinike
信息提示
关闭
请填处理意见
关闭