计算化学公社

标题: QE 优化Fe(111)面，不收敛，请指教 [打印本页]

作者
Author: 风飞 时间: 2021-1-14 15:19
标题: QE 优化Fe(111)面，不收敛，请指教
本帖最后由风飞于 2021-1-14 15:21 编辑

老师，
   我变胞优化了Fe的bulk结构，输入文件如下 (, 下载次数 Times of downloads: 6)
   然后搭建了Fe（111）面，但始终不收敛，请问该怎样修改呢？
输入文件 (, 下载次数 Times of downloads: 21)
输出文件太大了，传不了,  优化的过程中，一步也没有跑出来，就提示说不收敛
偃式用的是QE推荐的
部分输入文件
(, 下载次数 Times of downloads: 25)
部分输出文件
(, 下载次数 Times of downloads: 51)

作者
Author: rabbitkiller 时间: 2021-1-14 19:50
QE默认优化走200步，可能不够，改成2000步试试……

&control
nstep=2000
/

作者
Author: 卡开发发 时间: 2021-1-14 22:39
scf收敛标准没必要1e-13，无特定情况Fe做DFT+U也不是必要的。

作者
Author: 风飞 时间: 2021-1-14 23:11

卡开发发发表于 2021-1-14 22:39
scf收敛标准没必要1e-13，无特定情况Fe做DFT+U也不是必要的。

嗯，我最开始的收敛标准是10d-8，mixing_beta=0.3，结果不收敛，然后看了手册中关于Fe计算的例子，然后就慢慢提高到了10d-13 ，mixing_beta=0.05, 结果还是不收敛，
我想计算Fe Co Ni的111面，请问是不是无特殊情况，这些都不需要加DFT+U呢？
此外，对于这些元素的计算，在什么情况下需要做DFT+U呢？

作者
Author: 卡开发发 时间: 2021-1-14 23:52

风飞发表于 2021-1-14 23:11
嗯，我最开始的收敛标准是10d-8，mixing_beta=0.3，结果不收敛，然后看了手册中关于Fe计算的例子， ...

1e-8都难收敛，1e-13岂不是更难收敛？
DFT+U确实能消除一定的自作用，但U值选取也比较讲究（虽然这样也不见得就能保证得到好的结果），相应也会引起电子局域化导致SCF确实会难收敛一些，对于纯金属体系如果不是要得到好的能带结构，必要性不是特别大。不要啥都想往上加，先权衡这些参数的意义和计算成本。

作者
Author: 风飞 时间: 2021-1-15 00:01

卡开发发发表于 2021-1-14 23:52
1e-8都难收敛，1e-13岂不是更难收敛？
DFT+U确实能消除一定的自作用，但U值选取也比较讲究（虽然这样也 ...

嗯，好的谢谢老师，我记在我的小本本上了，我再去试试

作者
Author: lcq845402558 时间: 2021-1-15 00:14
这么大的截断能呀。。我自己算都不舍得用这么高的精度。羡慕。

作者
Author: 风飞 时间: 2021-1-15 00:50
本帖最后由风飞于 2021-1-15 01:04 编辑

lcq845402558 发表于 2021-1-15 00:14
这么大的截断能呀。。我自己算都不舍得用这么高的精度。羡慕。

这个是贗势里推荐的, 请问你计算Fe元素的时候，用的截断能是多少呢？

作者
Author: 卡开发发 时间: 2021-1-15 10:32

风飞发表于 2021-1-15 00:50
这个是贗势里推荐的, 请问你计算Fe元素的时候，用的截断能是多少呢？

应当测试确定。

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