计算化学公社

 找回密码 Forget password
 注册 Register
Views: 3586|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[Quantum ESPRESSO] 关于QE声子计算不收敛的问题的求助

[复制链接 Copy URL]

8

帖子

0

威望

71

eV
积分
79

Level 2 能力者

最近在QE计算声子谱的ph.x这一步骤的时候,一直会出现不收敛的情况






而且这个收敛情况为什么会由大到小,然后再由小到大,最后固定在一个数值一直不变做重复计算,不知道该怎么做了。

432

帖子

11

威望

3424

eV
积分
4076

Level 6 (一方通行)

2#
发表于 Post on 2022-8-21 11:02:43 | 只看该作者 Only view this author
声子计算你的收敛限设置了多少?tr2_th 默认是1.d^-12,你的收敛限离这个还差挺远
试试减小alpha_mix(1),默认0.7,并且增大nmix_ph到16
自由发挥,野蛮生长

78

帖子

0

威望

3308

eV
积分
3386

Level 5 (御坂)

3#
发表于 Post on 2022-8-21 11:25:51 | 只看该作者 Only view this author
可以检查看看结构优化过程和声子计算前一步的自洽计算是否收敛。

不知道大家比较过没有,QE计算声子谱比VASP慢了很多很多,如果开自旋极化和DFT+U就速度更慢了

8

帖子

0

威望

71

eV
积分
79

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2022-8-21 15:24:57 | 只看该作者 Only view this author
丁越 发表于 2022-8-21 11:02
声子计算你的收敛限设置了多少?tr2_th 默认是1.d^-12,你的收敛限离这个还差挺远
试试减小alpha_mix(1), ...

感谢 找到问题了 之前没有设置alpha_mix(1)的值 默认0.7太大了  减小之后发现正常收敛计算了

8

帖子

0

威望

71

eV
积分
79

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2022-8-21 15:27:47 | 只看该作者 Only view this author
waitingseven 发表于 2022-8-21 11:25
可以检查看看结构优化过程和声子计算前一步的自洽计算是否收敛。

不知道大家比较过没有,QE计算声子谱比 ...

vasp对于超胞的计算会出现一些不收敛的问题,导致虚频,想看看QE的计算情况如何,而且QE还可以续算,vasp倒是断了就没了,服务器提交时间有上限,不能用太长的时间去算vasp的声子谱

78

帖子

0

威望

3308

eV
积分
3386

Level 5 (御坂)

6#
发表于 Post on 2022-8-22 09:24:05 | 只看该作者 Only view this author
Dawn_4D 发表于 2022-8-21 15:27
vasp对于超胞的计算会出现一些不收敛的问题,导致虚频,想看看QE的计算情况如何,而且QE还可以续算,vasp ...

VASP的虚频是GGA泛函的原因,换到LDA结果会改善,采用DFPT的QE和VASP都不支持DFT-D3,我觉得QE有时候结果好一些是因为赝势种类多,总可以选出来一套结果没有虚频的赝势,对于声子谱这种计算量大的任务,VASP速度比QE快的太多了

1

帖子

0

威望

19

eV
积分
20

Level 1 能力者

7#
发表于 Post on 2024-5-11 17:26:55 | 只看该作者 Only view this author
丁越 发表于 2022-8-21 11:02
声子计算你的收敛限设置了多少?tr2_th 默认是1.d^-12,你的收敛限离这个还差挺远
试试减小alpha_mix(1), ...

您好  我按照您的方法,增大nmix_ph,但是nmix_ph默认最大是5,但是官网说可以有更大的取值(8-20),请问我应该如何修改使其能够设置更大的值呢

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 14:48 , Processed in 0.268249 second(s), 30 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list