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

风飞 · 发表于 Post on 2021-1-14 15:19:45

本帖最后由风飞于 2021-1-14 15:21 编辑

老师，
我变胞优化了Fe的bulk结构，输入文件如下

Fe.in (930 Bytes, 下载次数 Times of downloads: 6)
然后搭建了Fe（111）面，但始终不收敛，请问该怎样修改呢？
输入文件

Fe-111.in (1.18 KB, 下载次数 Times of downloads: 21)
输出文件太大了，传不了, 优化的过程中，一步也没有跑出来，就提示说不收敛
偃式用的是QE推荐的
部分输入文件

部分输出文件

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

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

&control
nstep=2000
/

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

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

风飞 · 发表于 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 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-15 00:01:16

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

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

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

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

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

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

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

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

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

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

应当测试确定。

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

[Quantum ESPRESSO] QE 优化Fe(111)面，不收敛，请指教