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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » Mg的K点测试
101/1返回列表
查看: 1853  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

钢金术师

金虫 (小有名气)

[求助] 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，请各路高手帮忙分析一下。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-03-10 20:06:51
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 钢金术师 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 05:44 by qL1cU46YIdK0
[基金申请] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 05:43 by qL1cU46YIdK0
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 05:43 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 8/400 2026-11-06 05:22 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 05:21 by qL1cU46YIdK0
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 13/650 2026-11-06 04:49 by qL1cU46YIdK0
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 04:47 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 PYZmcHr5kSaq 2026-11-05 7/350 2026-11-06 02:59 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:40 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 01:29 by qL1cU46YIdK0
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-05 4/200 2026-11-05 23:39 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭