请教计算光学性质时独立粒子近似和考虑局域场近似的区别

bushiyueliang

首先非常感谢您能打开这个帖子，我刚开始做光学性质的计算，有看到基于独立粒子近似和考虑局域场近似，但是没搞明白他们的区别以及如何选择，只知道后者是在前者的基础上INCAR中加入LRPA=.TRUE.
所以想请教一下大佬们什么情况下用哪种，一般计算下哪种更好，再次感谢

万里云

光学性质和逆介电矩阵的第(0, 0)个元素有关。计算的时候，如果把整个矩阵都算出来再取逆，就考虑了局域场效应。直接取第(0, 0)个元素的逆就没有考虑局域场效应。

RPA是算极化率的一种近似，算完极化率才能算介电矩阵。

光学性质计算量都不小，求矩阵逆的耗时可以忽略，所以建议一直考虑晶体场效应。

至于算得准不准，PBE+RPA给实验算足够了。一般理论计算需要HSE+RPA。以光学为卖点的计算建议GW+BSE。

如果做实验的要求定量符合测量结果，建议跑路。

bushiyueliang

万里云 发表于 2023-9-14 17:07
光学性质和逆介电矩阵的第(0, 0)个元素有关。计算的时候，如果把整个矩阵都算出来再取逆，就考虑了局域场效 ...

谢谢，有一点不明白，把整个矩阵都算出来再取逆这是依靠哪一步完成的？是考虑局域场效应吗？我在网上看到的教程是在独立粒子近似的条件下计算完，再用考虑局域场近似处理上一步的结果，这和您说的是一样的吗？我刚开始学习，对这方面不了解

万里云

bushiyueliang 发表于 2023-9-14 18:07
谢谢，有一点不明白，把整个矩阵都算出来再取逆这是依靠哪一步完成的？是考虑局域场效应吗？我在网上看到 ...

计算流程是先做scf得到基态电荷密度，再算bands得到价带和导带，然后基于RPA算微观极化率，接着从微观极化率算微观介电函数，最后取微观介电矩阵的逆矩阵的第(0, 0)个元素。

bushiyueliang

万里云 发表于 2023-9-14 18:51
计算流程是先做scf得到基态电荷密度，再算bands得到价带和导带，然后基于RPA算微观极化率，接着从微观极 ...

非常感谢！查了一下午也没您说的明白！

bushiyueliang

万里云 发表于 2023-9-14 18:51
计算流程是先做scf得到基态电荷密度，再算bands得到价带和导带，然后基于RPA算微观极化率，接着从微观极 ...

大佬您好，请问您有推荐的计算氧化物光学性质的实例吗？我试着用IPA和RPA的方法做计算，但是不会处理计算生成的结果，想问问您有没有详细的计算实例可以一步一步比对的，谢谢

万里云

bushiyueliang 发表于 2023-9-15 10:27
大佬您好，请问您有推荐的计算氧化物光学性质的实例吗？我试着用IPA和RPA的方法做计算，但是不会处理计算 ...

看vasp的官方教程。计算结果在vasprun.xml里。

柠檬鬼lemon

您好，我最近也在研究光学性质，我在官网上看到考虑局域场近似的方法是读取上一步独立例子近似的WAVCAR然后ALGO=CHI，请问在考虑局域场效应时是否是这样计算就可以了。此时LPRA默认=TRUE，若将LRPA=FALSE给出的解释是“为了包括RPA之外的局部场效应，即来自DFT交换和相关的贡献，必须指定LRPA=F"那么如果我现在需要用杂化泛函计算光学性质，是否应该将LRPA=F呀

柠檬鬼lemon

万里云 发表于 2023-9-14 17:07
光学性质和逆介电矩阵的第(0, 0)个元素有关。计算的时候，如果把整个矩阵都算出来再取逆，就考虑了局域场效 ...

您好，请问我现在在用杂化泛函+RPA做光学性质，计算结束后用vaspkit生成数据，因为材料是各向同性所以三个方向的数据是相同的，我就用xx方向做了图，但是得到的却是三条线，这是怎么回事呀，请问这三条数据线我应该怎么选择呀

tjuphysfang

万里云 发表于 2023-9-14 17:07
光学性质和逆介电矩阵的第(0, 0)个元素有关。计算的时候，如果把整个矩阵都算出来再取逆，就考虑了局域场效 ...

博主您好，请问这部分理论可以在那些地方学习啊？我刚开始做光学性质的计算，对于上面帖子提到的内容理论背景不是很懂。