计算化学公社

标题: 求助，黄铁矿（110）面几何优化一直失败，找不到解决办法。 [打印本页]

作者
Author: yassuo 时间: 2022-3-31 20:52
标题: 求助，黄铁矿（110）面几何优化一直失败，找不到解决办法。
最近一直在优化黄铁矿110表面，切出来的面已经放在附件，选择几何优化一直显示SCF不收敛。修改期间参照网上意见修改过cutoff energy（精度从粗到fine都试过）、K点（试过3x3x1,2x2x1）、density mixing中的charge（数值从0.1 改到0.5）、增加smearing值（目前加到0.3），增加empty band（目前是30%），增加SCF循环次数（目前已经加到500），但是都没有收敛，显示的问题也都是SCF不收敛的问题。同样的参数切黄铁矿111面可以收敛，所以非常困惑，目前已经不知道该修改什么参数才能收敛了

，希望有大佬能指点迷津，万分感谢！！！！

作者
Author: 卡开发发 时间: 2022-3-31 23:05
你看看你的晶体结构对称性是否为Pnnm？如果是这个对称性，那么这种FeS2没有磁性，按照非自旋极化算就行，也许这样更容易收敛。

作者
Author: yassuo 时间: 2022-4-1 15:14
本帖最后由 yassuo 于 2022-4-1 15:15 编辑

卡开发发发表于 2022-3-31 23:05
你看看你的晶体结构对称性是否为Pnnm？如果是这个对称性，那么这种FeS2没有磁性，按照非自旋极化算就行，也 ...

非常感谢回复。
我的晶体结构显示的对称性是P1，但是我也按照您提供的方法进行了尝试，在coarse精度下可以收敛。目前我想在此基础上进行更高精度的优化，但是在相同参数的情况下选择medium精度还是会出现SCF不收敛的情况，请问这种情况还有比较有效的解决办法吗？

作者
Author: 卡开发发 时间: 2022-4-1 16:39

yassuo 发表于 2022-4-1 15:14
非常感谢回复。
我的晶体结构显示的对称性是P1，但是我也按照您提供的方法进行了尝试，在coarse精度下可 ...

晶体结构你找一下对称性看看，判断一下是哪种结构，再去对照这个体系磁序情况。降低精度不会有任何帮助，应该给予足够的精度，如果结构没问题，先查磁序，如果不行再考虑调整迭代参数。

作者
Author: yassuo 时间: 2022-4-1 19:56
本帖最后由 yassuo 于 2022-4-1 20:03 编辑

卡开发发发表于 2022-4-1 16:39
晶体结构你找一下对称性看看，判断一下是哪种结构，再去对照这个体系磁序情况。降低精度不会有任何帮助， ...

您好，感谢回复。
请问关于晶体结构的分类和对应磁序有推荐的参考资料吗？对这方面知识我还比较欠缺，需要补一下。另外还有一个疑问：我对pyrite（111）和pyrite（100）两个面进行几何优化时都加了自旋极化并且都收敛了，是否需要去掉自旋极化重新提交计算呢，这对最终计算出的能量影响大吗？

作者
Author: 卡开发发 时间: 2022-4-1 20:25

yassuo 发表于 2022-4-1 19:56
您好，感谢回复。
请问关于晶体结构的分类和对应磁序有推荐的参考资料吗？对这方面知识我还比较欠缺，需 ...

这方面也只好看固体物理，至于磁序我一般还是参考Materials Project的网站，但也会找文献做确认。我不确定你的块体结构是否有磁性，如果确实没有，那么算出来不是基态的稳定结构，也就没有实际意义了。

作者
Author: yassuo 时间: 2022-4-1 21:10

卡开发发发表于 2022-4-1 20:25
这方面也只好看固体物理，至于磁序我一般还是参考Materials Project的网站，但也会找文献做确认。我不确 ...

不知道我理解的对不对：如果没有磁性，需要去掉自旋极化重新计算吗?

作者
Author: 卡开发发 时间: 2022-4-1 21:24

yassuo 发表于 2022-4-1 21:10
不知道我理解的对不对：如果没有磁性，需要去掉自旋极化重新计算吗?

对，所以课题要做好初步的调研再下手。

作者
Author: yassuo 时间: 2022-4-1 22:01

卡开发发发表于 2022-4-1 21:24
对，所以课题要做好初步的调研再下手。

好勒，非常感谢指导

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