计算化学公社

标题: k点密度对收敛性的影响 [打印本页]

作者
Author: jzjy1994 时间: 2017-9-16 16:55
标题: k点密度对收敛性的影响
请问一下各位老师k点密度是否会对LDA+U的单点能自洽的收敛性产生影响，因为最近所做的LDA+U体系在较为稀疏的k（0.055左右的密度）点时候收敛性很不错，大概跑个30几次就收敛了，然后k点一取密就变得极难收敛了（大约0.03）然后在计算能带时k点是可以取的比较密的，这种情况是因为什么造成的呢

作者
Author: 卡开发发 时间: 2017-9-16 19:15
数值的稳定性具体原因我不大好说，有很多因素会有影响，比如截断能、密度混合的方式等。比较难收敛可以考虑降低density mixing的参数，如果还不行可以考虑读取DFT的密度作为初猜来做，或者U小的先算收敛再慢慢增加。

作者
Author: jzjy1994 时间: 2017-9-17 09:45

卡开发发发表于 2017-9-16 19:15
数值的稳定性具体原因我不大好说，有很多因素会有影响，比如截断能、密度混合的方式等。比较难收敛可以考虑 ...

谢谢老师，我不太明白为什么k点增加会使其难收敛，另外低k点密度下计算的结构和能量是不是不可信，另外就是使一些没有自旋的体系带上自旋竟然都能帮助收敛。。。

作者
Author: jzjy1994 时间: 2017-9-17 09:46

卡开发发发表于 2017-9-16 19:15
数值的稳定性具体原因我不大好说，有很多因素会有影响，比如截断能、密度混合的方式等。比较难收敛可以考虑 ...

=。=sorry说错了，是给一个体系中没自旋的原子带上自旋都能帮助收敛

作者
Author: 卡开发发 时间: 2017-9-17 11:34

jzjy1994 发表于 2017-9-17 09:46
=。=sorry说错了，是给一个体系中没自旋的原子带上自旋都能帮助收敛

你要是用CASTEP的话这还不是最诡异的，最诡异的是初猜波函数好像是随机构造的，如果采用不同的random seed的话有时候相同的参数也会出现收敛和不收敛的情况。你说这种设置初始磁矩的方式有，但大部分体系应该用不着这样设置，+U确实数值稳定性会稍微差一些。

作者
Author: jzjy1994 时间: 2017-9-18 10:06

卡开发发发表于 2017-9-17 11:34
你要是用CASTEP的话这还不是最诡异的，最诡异的是初猜波函数好像是随机构造的，如果采用不同的random see ...

=。=是啊，老师，有时候有一个随机数效果就特别好，请教一下老师如何使得castep能使用和某一次模拟相同的随机数

作者
Author: 卡开发发 时间: 2017-9-18 11:15

jzjy1994 发表于 2017-9-18 10:06
=。=是啊，老师，有时候有一个随机数效果就特别好，请教一下老师如何使得castep能使用和某一次模拟相同的 ...

castep输出文件会打印出某次计算使用的seed ID，输入文件也支持给定seed ID，具体你可以研究一下Keywords，太简单我就不详细说了。

作者
Author: jzjy1994 时间: 2017-9-18 11:20

卡开发发发表于 2017-9-18 11:15
castep输出文件会打印出某次计算使用的seed ID，输入文件也支持给定seed ID，具体你可以研究一下Keywords ...

=0=老师，我能查到某次计算的seed ID，但是在其输入文件中输入相应的id貌似会报错，我有查关键词，上面说随机数只能选0或者大于0或者小于0，=。=seedID很长的一串。输进去就报错。

作者
Author: jzjy1994 时间: 2017-9-18 11:22

jzjy1994 发表于 2017-9-18 11:20
=0=老师，我能查到某次计算的seed ID，但是在其输入文件中输入相应的id貌似会报错，我有查关键词，上面说 ...

00或者说老师是我的关键词输错了。。。我之前输入的关键词是rand_seed

作者
Author: 卡开发发 时间: 2017-9-18 11:32

jzjy1994 发表于 2017-9-18 11:20
=0=老师，我能查到某次计算的seed ID，但是在其输入文件中输入相应的id貌似会报错，我有查关键词，上面说 ...

RAND_SEED应该是一串数字才是，手头现在没有castep所以没办法帮你具体测试。

作者
Author: jzjy1994 时间: 2017-9-18 11:36
RAND_SEED (.param)
Keyword type
Integer (expert)

Description
This keyword controls the initialization of random number seeds. Available options are:

0 - the seed for the random number generation is selected pseudorandomly.
> or < 0 - this value is used as a seed for the random number generator.
Default
0

Example
RAND_SEED : -25
See Also:
CASTEP keyword glossary
CASTEP parameters keywords

您好，老师，这是castep里面help给的提示。

作者
Author: 卡开发发 时间: 2017-9-18 12:03

jzjy1994 发表于 2017-9-18 11:36
RAND_SEED (.param)
Keyword type
Integer (expert)

是啊，按说你应该在.param文件当中其中一行加上
RAND_SEED : 你上一次计算CASTEP文件中打印的SEED_ID
这个ID应该是一串数。

作者
Author: jzjy1994 时间: 2017-9-18 14:27

卡开发发发表于 2017-9-18 12:03
是啊，按说你应该在.param文件当中其中一行加上
RAND_SEED : 你上一次计算CASTEP文件中打印的SEED_ID
...

=。=是啊，我再找一下资料看看吧，谢谢老师了，浪费了您很多时间

欢迎光临计算化学公社 (http://ccc.keinsci.com/)