计算化学公社

标题: QE中CeO2计算提示:not orthogonal operation错误求助 [打印本页]

作者
Author: 风飞 时间: 2021-1-11 22:17
标题: QE中CeO2计算提示:not orthogonal operation错误求助
老师，您好，
      使用QE 优化CeO2的111面，提示说：not orthogonal operation，出错求助
我使用的赝式是，QE中推荐的赝式： (, 下载次数 Times of downloads: 5)
我首先进行变胞优化，输入文件如下： (, 下载次数 Times of downloads: 8)
然后切111面，输入文件如下： (, 下载次数 Times of downloads: 7) ，部分输入如图
(, 下载次数 Times of downloads: 37)
结果出现错误，出错的out文件如下： (, 下载次数 Times of downloads: 2) ，出错提示如图
(, 下载次数 Times of downloads: 35)
提示：not orthogonal operation
我尝试加入U_projection_type='ortho-atomic'/‘norm-atomic'/‘atomic’/‘pseudo’  /cell_dofree='all'  也没有解决
于是我尝试用ps赝式 (, 下载次数 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的介绍如下：
(, 下载次数 Times of downloads: 30)
(在DFT+U运行的第一次迭代中，它覆盖了原子种类ityp的ispin分量的ns占用矩阵的第m个特征值。对于noncolin情况，ispin索引运行到npol。lmax值由应用Hubbard U的最大角动量数给出。保留未设置的特征值不变。)
(, 下载次数 Times of downloads: 44)

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

作者
Author: 风飞 时间: 2021-1-12 09:06
嗯，好的，谢谢老师

作者
Author: zyj19831206 时间: 2021-1-12 11:01

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

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

作者
Author: 卡开发发 时间: 2021-1-12 12:53

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

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

作者
Author: 风飞 时间: 2021-1-19 16:40

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

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

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