24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1527  |  回复: 9

钢金术师

金虫 (小有名气)

[求助] Mg的K点测试 已有2人参与

有哪位虫友做过单胞mg的K点测试,结果如何?我做单胞mg的K点测试发现K点一直不收敛,输入的INCAR 和测试的结果如下:
SYSTEM= Mg
ISTART= 0
ICHARG= 2
ISMEAR= 1
SIGMA= 0.2
EDIFF = 1.00E-05
EDIFFG = -0.01
PREC = Accurate
LWAVE = .FALSE.
LCHARG = .FALSE.
ENCUT = 550
IBRION = 1
NSW = 100
ISPIN = 1
NELM = 100
POTIM= 0.5
ISIF= 3
GGA= PE
结果:
7 *4 *4                -6.196808
8* 5*5                 -6.171068
9 *6*6                 -6.159724
10 *7*7                -6.180129
11* 8*8               -6.179336
12* 9*9               -6.170015
13 *10*10           -6.165933
14 *11 *11          -6.165145
15 *12 *12          -6.170765
16 *13*13           -6.172224
17 *14*14           -6.170126
18 *15*15           -6.169038
19 *16*16           -6.167489
赝势选用的是2个电子的mg的赝势。
测试的结果一支也没有收敛到0.01ev,请各路高手帮忙分析一下。
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangjiang12

铁虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2014-06-08 18:02:28
Mg是HCP结构吧,建议用ISIF = 6,你做收敛性测试之前,最好先优化出一个结构,然后用那个优化出来的结构做收敛性测试。小建议,PREC = NORMAL就可以了,ENCUT=300就可以了,没必要那么大
2楼2014-03-11 10:48:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dyc_2008

捐助贵宾 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2014-06-08 18:02:49
收敛分绝对收敛和相对收敛,你说的这个是绝对收敛。

相对收敛的话,要求要低得多。比如计算层错能,涉及到两个胞的能量相减。两个胞各自的能量如果按你上面的标准,是没有收敛的。但是他们的差值,却收敛了。换句话说,即使用更高密度的K点,也不会改变层错能。原因在于:相减的时候,两个胞的能量误差抵消了。

另:要是你看看文献里面的K点密度,其实也不那么高,相对收敛就行。还有就是晶格常数a,误差在5%都是可以接受的。
一布福田,二修菩提
3楼2014-03-11 14:28:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dyc_2008

捐助贵宾 (著名写手)

单胞的话  15 *12 *12     足够,我用的就差不多是这个密度
一布福田,二修菩提
4楼2014-03-11 14:30:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钢金术师

金虫 (小有名气)

引用回帖:
3楼: Originally posted by dyc_2008 at 2014-03-11 14:28:55
收敛分绝对收敛和相对收敛,你说的这个是绝对收敛。

相对收敛的话,要求要低得多。比如计算层错能,涉及到两个胞的能量相减。两个胞各自的能量如果按你上面的标准,是没有收敛的。但是他们的差值,却收敛了。换句 ...

但是一般在做k点测试的时候 不是应该最终会归一化的吗,误差应该达到0.01ev才对,但是这个mg却一直在波动
5楼2014-03-11 16:08:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钢金术师

金虫 (小有名气)

引用回帖:
4楼: Originally posted by dyc_2008 at 2014-03-11 14:30:35
单胞的话  15 *12 *12     足够,我用的就差不多是这个密度

但是这个不收敛的原因在哪里啊,我怕我的参数不对,我又做了单胞钛的,收敛的很好,只有这个mg的很不一样啊
6楼2014-03-11 16:09:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzjust

新虫 (小有名气)

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2014-06-08 18:02:58
我最近算Fe和Mg系,也是出现这样的情况。
总是害怕自己的参数设置不正确,而且确实是和别的文献中的计算结果差距较大(都用CASTEP)。
楼主现在的问题解决了吗?
7楼2014-06-03 08:46:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钢金术师

金虫 (小有名气)

引用回帖:
7楼: Originally posted by zzjust at 2014-06-03 08:46:05
我最近算Fe和Mg系,也是出现这样的情况。
总是害怕自己的参数设置不正确,而且确实是和别的文献中的计算结果差距较大(都用CASTEP)。
楼主现在的问题解决了吗?

我用的VASP ,对于Mg,K点采用G方法(如果是M方法,奇数和偶数分开测试),测试到最后会收敛的,只不过K点数量比较多
8楼2014-06-03 10:45:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我叫小殊子

铜虫 (小有名气)

【答案】应助回帖

引用回帖:
7楼: Originally posted by zzjust at 2014-06-03 08:46:05
我最近算Fe和Mg系,也是出现这样的情况。
总是害怕自己的参数设置不正确,而且确实是和别的文献中的计算结果差距较大(都用CASTEP)。
楼主现在的问题解决了吗?

在算Mg和Fe结构优化时老是出现错误,请问一下是什么原因呢?

[ 发自小木虫客户端 ]
甘肃正天合律师事务所专职律师
9楼2014-06-04 01:47:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钢金术师

金虫 (小有名气)

引用回帖:
9楼: Originally posted by 我叫小殊子 at 2014-06-04 01:47:04
在算Mg和Fe结构优化时老是出现错误,请问一下是什么原因呢?
...

具体什么错误啊?
10楼2014-06-04 08:58:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 钢金术师 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料求调剂 一志愿哈工大总分298分,前三科223分 +5 dongfang59 2026-03-27 5/250 2026-03-28 04:53 by wxiongid
[考研] 085701环境工程求调剂 +9 多久上课 2026-03-27 9/450 2026-03-28 03:58 by fmesaito
[考研] 0856材料化工调剂 总分330 +10 zhubinhao 2026-03-27 10/500 2026-03-28 03:34 by fmesaito
[考研] 【求调剂】085601材料工程专硕 | 总分272 | +5 脚滑的守法公民 2026-03-27 5/250 2026-03-27 22:50 by lcl13500..
[考研] 266分,求材料冶金能源化工等调剂 +5 哇呼哼呼哼 2026-03-27 7/350 2026-03-27 21:41 by lbsjt
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +7 18419759900 2026-03-25 8/400 2026-03-27 16:38 by 18419759900
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +9 vv迷 2026-03-22 9/450 2026-03-27 15:59 by 不吃魚的貓
[考研] 一志愿北化085600材料专硕275|有文章专利|求调剂 +3 Micky11223 2026-03-25 3/150 2026-03-27 14:52 by 醉在风里
[考研] 化学308分求调剂 +8 你好明天你好 2026-03-23 9/450 2026-03-27 14:01 by 杨光于青云
[考研] 求调剂 +3 刘柯@ 2026-03-24 4/200 2026-03-27 11:28 by shangxh
[论文投稿] Journal of Mechanical Science and Technology +3 Russ_ss 2026-03-25 5/250 2026-03-27 10:49 by 陆小果画大饼
[考研] 325求调剂 +5 李嘉图·S·路 2026-03-23 5/250 2026-03-27 00:42 by wxiongid
[考研] 336材料求调剂 +7 陈滢莹 2026-03-26 9/450 2026-03-27 00:20 by wxiongid
[考研] 325求调剂 +3 Aoyijiang 2026-03-23 3/150 2026-03-26 20:46 by 不吃魚的貓
[考研] 一志愿厦门大学化学学硕307求调剂 +8 y7czhao 2026-03-26 8/400 2026-03-26 19:51 by 不吃魚的貓
[考研] 281求调剂 +3 亚克西good 2026-03-26 5/250 2026-03-26 19:48 by 不吃魚的貓
[考研] B区考研调剂 +4 yqdszhdap- 2026-03-22 5/250 2026-03-25 08:51 by baoball
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
信息提示
请填处理意见