24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2134  |  回复: 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 的主题更新
信息提示
请填处理意见