计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First-principles) › QE 优化Fe(111)面，不收敛，请指教

QE 优化Fe(111)面，不收敛，请指教

查看数: 10253 | 评论数: 8 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

风飞

发布时间: 2021-1-14 15:19

正文摘要:

本帖最后由风飞于 2021-1-14 15:21 编辑老师，我变胞优化了Fe的bulk结构，输入文件如下

回复 Reply

卡开发发 发表于 Post on 2021-1-15 10:32:48

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

应当测试确定。

风飞 发表于 Post on 2021-1-15 00:50:39

本帖最后由风飞于 2021-1-15 01:04 编辑

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

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

lcq845402558 发表于 Post on 2021-1-15 00:14:59

这么大的截断能呀。。我自己算都不舍得用这么高的精度。羡慕。

风飞 发表于 Post on 2021-1-15 00:01:16

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

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

卡开发发 发表于 Post on 2021-1-14 23:52:33

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

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

风飞 发表于 Post on 2021-1-14 23:11:36

卡开发发发表于 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呢？

卡开发发 发表于 Post on 2021-1-14 22:39:10

scf收敛标准没必要1e-13，无特定情况Fe做DFT+U也不是必要的。

rabbitkiller 发表于 Post on 2021-1-14 19:50:04

QE默认优化走200步，可能不够，改成2000步试试……

&control
nstep=2000
/

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register