计算化学公社

标题: 求助：VASP计算出的磁态分布与文献不符 [打印本页]

作者
Author: yuzj41d 时间: 2022-5-17 17:18
标题: 求助：VASP计算出的磁态分布与文献不符
本帖最后由 yuzj41d 于 2022-5-17 17:30 编辑

在对CeO2含氧缺陷体系计算时，模拟了上表面O缺陷的情况如图所示
在文献DOI: 10.1103/PhysRevB.79.193401 中，O缺陷的第二近邻Ce原子也就是 4、6、12三个Ce两个处于+3价时是最稳定的，文中是6、12。
于是我分别进行了以下模拟分别的用默认铁磁、6 12、 6 8、分别为Ce3+进行模拟
但是结果磁态分别分布在8 9、6 9 、8 9上也就是说基本至少有一个是第一近邻，我百思不得其解。请各位大佬解惑。
(, 下载次数 Times of downloads: 17)
(, 下载次数 Times of downloads: 21) (, 下载次数 Times of downloads: 15) (, 下载次数 Times of downloads: 9)

三次磁态分布（部分OUTCAR截图）

(, 下载次数 Times of downloads: 14) (, 下载次数 Times of downloads: 15) (, 下载次数 Times of downloads: 18)

作者
Author: yuzj41d 时间: 2022-5-19 11:04
希望大佬解答

作者
Author: 乘风万里 时间: 2022-5-19 16:03
本帖最后由乘风万里于 2022-5-19 16:17 编辑

DFT+U的自洽结果与初始的 occupation matrix有关，如果不设，很容易收敛到亚稳态，文献[1]中对算UO2有详细的讨论。文献[2][3]讨论了怎么设置occupation matrix能算出低能态。
Zunger 18年计算NiO, CoO的文章[4]引了以上的文献，并按照以上文献设置occupation matrix，以防算出亚稳态。
我自己算LixCoO2的时候，LMAXMIX设不设，自洽后的局域磁矩都会变化，虽然VASP wiki上说LMAXMIX参数只是加快收敛。甚至按说明书d电子体系设LMAXMIX=4，反而自洽出一个亚稳态，不设LMAXMIX时，自洽结果的能量却更低，应该也是跟初始占据有关。

1. PHYSICAL REVIEW B 79, 235125 (2009)
2. PHYSICAL REVIEW B 82, 195128 (2010)
3. Phys. Chem. Chem. Phys., 2014, 16, 21016—21031
4. PHYSICAL REVIEW B 97, 035107 (2018)

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