24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4385  |  回复: 19

1244464635

金虫 (正式写手)

引用回帖:
8楼: Originally posted by c_chi at 2014-06-03 14:56:30
谢谢你的指点哈,我这里还有几点疑问
1. 我觉得在做结构优化的时候可以不加磁性,然后在算能量和其他性质的时候再加磁性参数,这样计算量会小一些
2. 我是用psl生成的赝势,有两种kjpaw和rrkjus,应该是分别对应 ...

1.这个我还真的不清楚,你有试过么?要是磁性对于结构优化影响不大的话,确实是可以再relax时不加磁性。
2.对于赝势我了解得很少,我用的都是赝势库里的赝势,我用的Mn的赝势和你用的一样,这个(van=vanderbilt ultrosoft)就是超软的啊。我看过一些关于赝势的帖子,大家貌似都认为用自己生成的赝势需经过严格测试才行。
3.超软相比于PAW势确实老了些。与其他用vasp的同行交流之后我发现PWSCF相对于vasp计算速度确实慢一些。
4.mixing-mode的选择应该只是与收敛速度关系较大。
5.我是直接将electron-maxstep设置得大一些。我个人感觉前几个scf(尤其是第一个)能在400步内收敛就可以接受,主要是看后面的scf收敛速度。
6.加U对于结构和能量还是有较大影响的,尤其是对于有4f电子的体系,不过大家对于3d电子的研究好像没那么多。
7.AIMD=ab inito molecule dynamic?这个我没做过..
8.可以试试调小mixing-beta,我最小用过0.01。degauss不能再加大了,不过在relax快要结束时应当适当地减小degauss。
9.upscale:在结构优化快要完成时,conv-thr会变得很小,调节这个参数可使conv-thr不至于变得过小;
starting-ns-eigenvalue:设置原子轨道的初始占据态。其他的好像没用过。
对于pwscf我也有很多不懂得东西,好多东西都是自己瞎试,也不知理解的是否到位。大家以后继续讨论哈。
11楼2014-06-03 16:12:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1244464635

金虫 (正式写手)

引用回帖:
9楼: Originally posted by c_chi at 2014-06-03 15:05:31
问一下,你有没有做过电催化?...

没做过啊
12楼2014-06-03 16:12:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)

引用回帖:
11楼: Originally posted by 1244464635 at 2014-06-03 16:12:12
1.这个我还真的不清楚,你有试过么?要是磁性对于结构优化影响不大的话,确实是可以再relax时不加磁性。
2.对于赝势我了解得很少,我用的都是赝势库里的赝势,我用的Mn的赝势和你用的一样,这个(van=vanderbilt u ...

恩,我也是刚刚接触pwscf不久,大家相互学习,也期待大牛们的指点哈。
我有试过结构优化不加磁性,一般没什么问题,而且在结构优化的时候可以降低K点和cutoff。
赝势我是从官网上下的赝势库pslibrary.1.0.0.tar.gz生成的。
对于加U我是担心加了之后就更慢了。
AIMD是ab inito molecule dynamic,如果做结构优化都这么慢,那AIMD就是在难以忍受了,要绝望了。。pwscf的PAW势又不太成熟。
upscale不是有默认值吗?是要将这个值调大?
starting_ns_eigenvalue是在加U的时候使用的吧?应该如何调整比较合适?
13楼2014-06-03 16:41:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1244464635

金虫 (正式写手)

引用回帖:
13楼: Originally posted by c_chi at 2014-06-03 16:41:32
恩,我也是刚刚接触pwscf不久,大家相互学习,也期待大牛们的指点哈。
我有试过结构优化不加磁性,一般没什么问题,而且在结构优化的时候可以降低K点和cutoff。
赝势我是从官网上下的赝势库pslibrary.1.0.0.tar. ...

upscale默认值为10,一般将其调小。
starting_ns_eigenvalue确实是在加U时用的,根据原子的价电子数,可以给其某些轨道预设占据态为1。譬如说Ce3+只有一个4f电子,那么我就可以将其f轨道第七个eigenvalue设为1。
14楼2014-06-03 17:26:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)

引用回帖:
14楼: Originally posted by 1244464635 at 2014-06-03 17:26:17
upscale默认值为10,一般将其调小。
starting_ns_eigenvalue确实是在加U时用的,根据原子的价电子数,可以给其某些轨道预设占据态为1。譬如说Ce3+只有一个4f电子,那么我就可以将其f轨道第七个eigenvalue设为1。...

好的,谢谢!有问题再交流
15楼2014-06-04 19:32:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pangrui1985

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
c_chi: 金币+5, ★★★很有帮助 2014-07-04 13:50:47
打开氧的PP,有这一句
Suggested minimum cutoff for wavefunctions:  47. Ry
因此,你的截断远远不够,这是不收敛的原因
16楼2014-07-04 09:46:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)

引用回帖:
16楼: Originally posted by pangrui1985 at 2014-07-04 09:46:50
打开氧的PP,有这一句
Suggested minimum cutoff for wavefunctions:  47. Ry
因此,你的截断远远不够,这是不收敛的原因

好的,谢谢!
17楼2014-07-04 13:50:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
18楼2015-03-21 09:47:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

周虫虫

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★
c_chi: 金币+3, 好的,谢谢 2016-05-19 03:11:34
我用PWscf也遇到了同样的问题,调了很多参数都不行,至今也没有彻底解决这个问题,但还是有一点经验:
1、截断能的设置可以稍微高一些,可以把更多的平面波包含进来
2、赝势选用冻结核心半径小的,超软的不行可以换模守恒的试一下
3、将混合因子 mixing_beta值降低,这个值是更新电荷密度的时候新的电荷密度所占的权重,降低新电荷密度的权重可以在之前的电荷密度基础上缓慢靠近最小总能对应的电荷密度
4、降低精度,即将总能收敛阈值调高
5、用Langevin动力学方法,这种方法总能计算不需要达到很高精度,但其他性质计算比较好
6、增加原子分子数目,原子数过少,单个原子的运动对体系影响很大
以上是我的一些经验,欢迎批评讨论~
扎根学术
19楼2016-05-18 22:38:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

applezju

木虫 (正式写手)

楼主,我的计算也是电子步很难收敛,600步都不收敛。能否麻烦您看看是什么参数设置的有问题,谢谢
http://muchong.com/bbs/viewthread.php?tid=10401175&fpage=1
混,混混
20楼2016-06-13 21:37:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 c_chi 的主题更新
信息提示
请填处理意见