VASP计算frequency dependent 介电函数矩阵问题求解

灰飞的旋律

有一些关于vasp计算光学性质的问题想各位帮我解答下，刚接触光学性质计算这块，也看了手册上各个参数的解释还有官方的教程各种示例，网络资料，然而还是有些问题没想明白
1. 在计算frequency dependent介电函数矩阵时，也就是LOPTICS=.TRUE.时，哈密顿量的精确对角化（ALGO=EXACT）是必要的吗？手册里是建议当计算频率依赖介电函数矩阵时，因为NABNDS=默认值2~3倍，这样空带数就超过了30%，建议加上exact，那不加会对介电函数矩阵有什么影响呢，结果还可信？
2. 我的体系时3个原子的原胞，在进行这种exact计算时，内存竟然占到了60G，感觉有点不可思议啊，请问大家都是这样的吗？计算个小原胞也这么耗内存？补充下我的参数设置吧，我是基于PBE泛函的基态计算结果然后增加空带数进行exact计算也就是独立粒子近似IPA计算，用的一个节点，28核，64G内存，体系是六方晶系二维单层结构，上下各10A真空层。
INCAR如图1：

，KPOINTS图2，

3.  采用密度泛函微扰理论DFPT（LEPSILON = .TRUE.）算出来的静态介电张量static dielectic tensor和频率依赖的介电函数矩阵之间是什么关系，如果想算光吸收系数只能用频率依赖的是吧？静态介电张量的用处是什么呢？
4.  如果我想计算个超胞体系的光学性质，总共27个原子。其中9个第四周期金属原子，其余均为非金属，还是基于PBE基态计算，采用RPA近似计算频率依赖的介电函数矩阵，大概估测得多大内存才能满足呢？
5. 还有个就是关于下面这个官方教程的理解，step2和3都没有设置LRPA的值，那为啥2是IPA，3就是RPA呢，该怎么理解呢

因为刚开始接触，有些问题可能比较浅薄，还请不吝赐教或者指点一下深入学习了解的资料，方向啥的

小毛毛

您好！我的体系计算出来的outcar一直默认的是RPA（独立粒子近似），怎么去选择独立粒子近似呢？两者的区别是什么啊！谢谢谢谢

灰飞的旋律

小毛毛 发表于 2019-5-20 11:40
您好！我的体系计算出来的outcar一直默认的是RPA（独立粒子近似），怎么去选择独立粒子近似呢？两者的区别 ...

有个叫LRPA的参数，你看看手册是咋说的

小毛毛

我尝试过将其设置为T或者F，得出来的结果都是RPA（随机相近似），是否可以给我参考一下您的输入文件。我们可以加QQ吗？我的QQ是1579041114

卡开发发

小毛毛 发表于 2019-5-20 20:32
我尝试过将其设置为T或者F，得出来的结果都是RPA（随机相近似），是否可以给我参考一下您的输入文件。我们 ...

我想可能是ALGO=CHI

小毛毛

ALGO应该是可以在IPA和DPA条件下计算频率依赖介电函数吧

黑挖掘机

请问计算介电函数第二步时，ALGO=Exact，任务被自动停止，日志出现killed with signal 9 ，这是什么问题呢？当使用ALGO = Normal时能正常计算，基于此计算基础，进行第三步时，ALGO=CHI，OUTCAR里没有INVERSE MACROSCOPIC和HEAD OF MICROSCOPIC的信息，这是什么问题呢？请问是否需要严格执行以上计算参数？不胜感激

灰飞的旋律

黑挖掘机 发表于 2021-11-5 01:31
请问计算介电函数第二步时，ALGO=Exact，任务被自动停止，日志出现killed with signal 9 ，这是什么问题呢 ...

ALGO=Exact，任务被自动停止，日志出现killed with signal 9 ，，内存不够用了吧

DCL

请问楼主，明白以上问题了吗？我才接触十分迷惑，还希望楼主能指点一二

Dominic

有人试过官方给出的Si例子吗，我用两个服务器测试，还是官方给的输入文件，算出来Si在所有能量范围内的吸收都是零，有人遇到过这种问题吗？

学不会计算化学

可不可以加个qq，讨论一下计算光学性质遇到的问题

tjuphysfang

卡开发发 发表于 2019-5-20 22:50
我想可能是ALGO=CHI

wiki中SiC的例子分别用了当ALGO=CHI时 LRPA=T和LRPA=F的结果但是无论在INCAR中设置T还是F，OUTCAR结果都是(including local field effects in RPA (Hartree))，这是为什么呀 卡卡老师。

卡开发发

tjuphysfang 发表于 2024-10-29 15:25
wiki中SiC的例子分别用了当ALGO=CHI时 LRPA=T和LRPA=F的结果但是无论在INCAR中设置T还是F，OUTCAR结果都 ...

要么INCAR和OUTCAR传来看看？

tjuphysfang

卡开发发 发表于 2024-10-29 20:13
要么INCAR和OUTCAR传来看看？

问题解决了，卡卡老师，但是为什么当LRPA = F配合 ALGO = CHI的时候，output会提醒下面的话呢？
这边我重新跑了一下程序，还没算完 但是output中会有这个警告。
-----------------------------------------------------------------------------
|                                                                             |
|           W    W    AA    RRRRR   N    N  II  N    N   GGGG   !!!           |
|           W    W   A  A   R    R  NN   N  II  NN   N  G    G  !!!           |
|           W    W  A    A  R    R  N N  N  II  N N  N  G       !!!           |
|           W WW W  AAAAAA  RRRRR   N  N N  II  N  N N  G  GGG   !            |
|           WW  WW  A    A  R   R   N   NN  II  N   NN  G    G                |
|           W    W  A    A  R    R  N    N  II  N    N   GGGG   !!!           |
|                                                                             |
|      You calculate the electron-hole interaction using AEXX=0.              |
|      This is hardly what you want to do.                                    |
|      Maybe you have forgotten to set AEXX in the INCAR file.                |
|    ( For LHFCALC=.TRUE. the default is  AEXX=0.25, but if LHFCALC is not    |
|      set the default is AEXX=0.0 )                                          |
|                                                                             |
-----------------------------------------------------------------------------

卡开发发

tjuphysfang 发表于 2024-10-30 17:29
问题解决了，卡卡老师，但是为什么当LRPA = F配合 ALGO = CHI的时候，output会提醒下面的话呢？
这边我 ...

看上去似乎用LRPA=F的时候会自动设置成LHFCALC=T。