QE中CeO2计算提示:not orthogonal operation错误求助

风飞 · 发表于 Post on 2021-1-11 22:17:02

老师，您好，
使用QE 优化CeO2的111面，提示说：not orthogonal operation，出错求助
我使用的赝式是，QE中推荐的赝式：

Ce.GGA-PBE-paw-v1.0.UPF (3.18 MB, 下载次数 Times of downloads: 5)
我首先进行变胞优化，输入文件如下：

CeO2-TJ.in (1.4 KB, 下载次数 Times of downloads: 8)
然后切111面，输入文件如下：

CeO2-TJ-111-1.in (1.68 KB, 下载次数 Times of downloads: 7) ，部分输入如图

结果出现错误，出错的out文件如下：

CeO2-TJ-111-1.out (494.23 KB, 下载次数 Times of downloads: 2) ，出错提示如图

提示：not orthogonal operation
我尝试加入U_projection_type='ortho-atomic'/‘norm-atomic'/‘atomic’/‘pseudo’ /cell_dofree='all' 也没有解决
于是我尝试用ps赝式

Ce.pbe-spdfn-kjpaw_psl.1.0.0.UPF (2.99 MB, 下载次数 Times of downloads: 1) ：Ce.pbe-spdfn-kjpaw_psl.1.0.0.UPF ，其他参数不变
但是在变胞计算的时候，发现得到的晶胞参数(4.3746)和实验(5.4976)相差的比较的大。

请问该怎样解决呢？

此外，请问QE 中的starting_ns_eigenvalue参数怎么设置呢？QE中starting_ns_eigenvalue的介绍如下：

(在DFT+U运行的第一次迭代中，它覆盖了原子种类ityp的ispin分量的ns占用矩阵的第m个特征值。对于noncolin情况，ispin索引运行到npol。lmax值由应用Hubbard U的最大角动量数给出。保留未设置的特征值不变。)

卡开发发 · 发表于 Post on 2021-1-12 01:22:19

1、和U的projector是否正交无关，可能和对称性有关。比如密度网格使得对称性破缺。可能能够通过调整密度截断或者试试不要用对称性。
2、一般情况下可以不用设置，不过如果收敛到错误的电子态或因d、f等占据问题SCF无法收敛可以试着调整下初始占据，这里指定磁量子数、alpha或beta轨道，元素种类和占据数来设置这个占据矩阵。

风飞 · 发表于 Post on 2021-1-12 09:06:31

嗯，好的，谢谢老师

zyj19831206 · 发表于 Post on 2021-1-12 11:01:27

卡开发发发表于 2021-1-12 01:22
1、和U的projector是否正交无关，可能和对称性有关。比如密度网格使得对称性破缺。可能能够通过调整密度截 ...

卡老师，QE有身边么比较好论坛，或者比较全的手册是哪个网址？

卡开发发 · 发表于 Post on 2021-1-12 12:53:17

zyj19831206 发表于 2021-1-12 11:01
卡老师，QE有身边么比较好论坛，或者比较全的手册是哪个网址？

我不太清楚诶，就是靠着官方的一些lecture的资料和“科学搜索”。
不过没有关系，可以先试着用一下，我觉得平面波方法的程序逻辑上大体都是差不多的，有什么问题或者想法论坛上一块讨论分享一下也不错。

风飞 · 发表于 Post on 2021-1-19 16:40:42

卡开发发发表于 2021-1-12 01:22
1、和U的projector是否正交无关，可能和对称性有关。比如密度网格使得对称性破缺。可能能够通过调整密度截 ...

老师，您好，最终，我是通过不使用对称性解决的

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

[Quantum ESPRESSO] QE中CeO2计算提示:not orthogonal operation错误求助

评分 Rate

评分 Rate

浏览过的版块