24小时热门版块排行榜    

查看: 2900  |  回复: 43

刘仕晨

木虫 (正式写手)

引用回帖:
28楼: Originally posted by gemucai at 2013-03-28 16:56:00
啊,你用的是新版的pwscf?这么智能了
如果换成以前,他会告诉你说nr3设置得不正确,你就要一个个试,现在给你自动调整,挺好的。
有些数字因为约数比较少,开发者可能觉得不太保险,不能保证FFT格点上都有原子,所 ...

没有最新版,是4.3.2的,我刚刚上面说的那个问题,第一个是没有设置nr1 nr2 nr3 ,我就是想看看两个结果,然后设置了,增加了nr3,却发现能量高了一点点
31楼2013-03-28 17:00:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
30楼: Originally posted by gemucai at 2013-03-28 17:00:39
不知道你有没有做过ecut的测试,到最后会有幅度较小的震荡,这是正常的,这时候我们认为它已经收敛。...

我现在就是在做ecutwfc测试,在从90,慢慢增加,步数为10,在增加到140的时候,出现了这个对称性丢失警告,所以想增加nr1 nr2 nr3继续算……
32楼2013-03-28 17:02:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
32楼: Originally posted by 刘仕晨 at 2013-03-28 17:02:47
我现在就是在做ecutwfc测试,在从90,慢慢增加,步数为10,在增加到140的时候,出现了这个对称性丢失警告,所以想增加nr1 nr2 nr3继续算……...

每隔10Ry算一次,差别这么小,可以说能量收敛了。到最后会有一定的小起伏。
33楼2013-03-28 18:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
33楼: Originally posted by gemucai at 2013-03-28 18:09:25
每隔10Ry算一次,差别这么小,可以说能量收敛了。到最后会有一定的小起伏。...

不是的,当ecutwfc=150时,总能为-167.74382003 Ry,当ecutwfc=160时,总能量为-167.78218062 Ry,相差0.52ev,我觉得是不是还需要继续加大ecutwfc?正常应该收敛到相差10mev么,就是奇怪,居然要这么大的ecut
34楼2013-03-29 13:26:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
34楼: Originally posted by 刘仕晨 at 2013-03-29 13:26:13
不是的,当ecutwfc=150时,总能为-167.74382003 Ry,当ecutwfc=160时,总能量为-167.78218062 Ry,相差0.52ev,我觉得是不是还需要继续加大ecutwfc?正常应该收敛到相差10mev么,就是奇怪,居然要这么大的ecut...

是不是结构有问题
35楼2013-03-29 14:41:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
34楼: Originally posted by 刘仕晨 at 2013-03-29 13:26:13
不是的,当ecutwfc=150时,总能为-167.74382003 Ry,当ecutwfc=160时,总能量为-167.78218062 Ry,相差0.52ev,我觉得是不是还需要继续加大ecutwfc?正常应该收敛到相差10mev么,就是奇怪,居然要这么大的ecut...

可以尝试继续增加k点,但我估计意义不大。

我之前也碰到过类似的情况,pwscf的赝势质量实在是堪忧。我当时算的结果,也是有发散的趋势:
http://muchong.com/bbs/viewthread.php?tid=3934098

所以我觉得,不必拘泥于这些细节,你的ecut已经够高了。
36楼2013-03-29 14:44:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
35楼: Originally posted by cj4566 at 2013-03-29 14:41:16
是不是结构有问题...

我上面提供了我的IN文件的,就是ecut不一样,结构应该没有问题,我用xcrysden看过了,应该问题不大……你可以看看
37楼2013-03-29 14:57:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
36楼: Originally posted by gemucai at 2013-03-29 14:44:40
可以尝试继续增加k点,但我估计意义不大。

我之前也碰到过类似的情况,pwscf的赝势质量实在是堪忧。我当时算的结果,也是有发散的趋势:
http://muchong.com/bbs/viewthread.php?tid=3934098

所以我觉得,不必 ...

我再继续加大看看,因为我以前用vasp老板都是要求到mev级别的,K点是问题不大,因为我作了K点的收敛性测试的,两个相邻的只有0.5mev,所以现在只在增加ecut。我这还只是自洽,不知道这么大的ecut,跑费米面能不能跑的动呢……我PC机先学习的,I7-3770,8G
38楼2013-03-29 15:00:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
38楼: Originally posted by 刘仕晨 at 2013-03-29 15:00:02
我再继续加大看看,因为我以前用vasp老板都是要求到mev级别的,K点是问题不大,因为我作了K点的收敛性测试的,两个相邻的只有0.5mev,所以现在只在增加ecut。我这还只是自洽,不知道这么大的ecut,跑费米面能不能跑 ...

总之我觉得pwscf的赝势质量不行,不过有时候不管收敛判据,死马当活马医也能出一些结果。
39楼2013-03-29 15:02:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)

引用回帖:
39楼: Originally posted by gemucai at 2013-03-29 15:02:21
总之我觉得pwscf的赝势质量不行,不过有时候不管收敛判据,死马当活马医也能出一些结果。...

是的,之前我感觉ecut=30差不多了,就这样算的,费米面什么的都画出来了,不过后来发现都错了,你能给我一点你选官网的势的经验么,比如你算一个体系,你一般选那些势,选谁的,或者有一些你的经验?还有我们画band图,一定要用example里的gnuplot软件么?
40楼2013-03-29 16:49:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cj4566 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 G2rhSvqvhNzl 2026-11-10 4/200 2026-11-11 08:34 by 1Nb0IVlmV6v0
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 4/200 2026-11-11 08:12 by 1Nb0IVlmV6v0
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 07:03 by Z7qfCsdIiD1T
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 9/450 2026-11-11 04:54 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:23 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-11 04:21 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:10 by Z7qfCsdIiD1T
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:07 by Z7qfCsdIiD1T
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-11 01:24 by Z7qfCsdIiD1T
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 22:39 by Z7qfCsdIiD1T
[考博] 申请环境,微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +4 yexuqing 2026-11-04 6/300 2026-11-10 11:13 by susizheng
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售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 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
信息提示
请填处理意见