学习交流：VASP计算激发态构型时INCAR的设置等问题

石不转

参考文章（https://zhuanlan.zhihu.com/p/454376223）通过改变电子占据情况来描述激发态的产生，我仿照文章设置了适用于自己体系的INCAR，成功计算出了激发态构型。之前的计算失败有两方面原因：一是FERWE标签名拼写错误；二是忽视了参考文章中一些非关键但也有影响的参数。谨以此帖记录学习过程中遇到的问题和解决方案，希望前辈们多多指点其中的错误和纰漏，先在此一并致谢。

我使用的INCAR除了MAGMOM、FERWE、FERDO三项之外，其余均与参考文章中相同，相同部分就不再赘述：
MAGMOM设置为 54*0（无磁性体系， NIONS=54）

1. MAGMOM = 54*0

复制代码

FERWE、FERDO设置为适用于自己体系多个K点的形式（NKPTS=120, NBANDS=320）。
由于K点数量多，FERWE、FERDO较冗长（目前没有找到简写的方式），以下仅放一小段作为示例。不同K点占据情况之间相隔了3个空格。

1. FERWE = 217*1.0 1*0.0 1*1.0 101*0.0   217*1.0 1*0.0 1*1.0 101*0.0
   
2. FERDO = 218*1.0 102*0.0   218*1.0 102*0.0   218*1.0 102*0.0

复制代码

盈朕

你好，这篇参考文章（https://zhuanlan.zhihu.com/p/454376223）文章原文找不到了，您还有吗
FERWE与FERDO具体如何设置有教程吗

石不转

盈朕 发表于 2022-10-20 19:35
你好，这篇参考文章（https://zhuanlan.zhihu.com/p/454376223）文章原文找不到了，您还有吗
FERWE与FERDO ...

您好。FERWE与FERDO的设置可参考VASP官网。官网原话中“The occupancies must be specified for all bands and k points”且“The occupancies must be between 0 and 1”。具体设置方法可参照以下该文章作者的原话（如有冒犯或侵权请老师告知我删除）：

grep NKPTS OUTCAR 查看具体有多少个K点，假定有五个，且最初你一个K点设置如下，
FERWE = 11*1.0 1*0.0 1*1.0 15*0.0
FERDO = 11*1.0 1*1.0 1*0.0 15*0.0
现在你重复五次，即可设置成功。
FERWE = 11*1.0 1*0.0 1*1.0 15*0.0  11*1.0 1*0.0 1*1.0 15*0.0 11*1.0 1*0.0 1*1.0 15*0.0 11*1.0 1*0.0 1*1.0 15*0.0 11*1.0 1*0.0 1*1.0 15*0.0
FERDO = 11*1.0 1*1.0 1*0.0 15*0.0  11*1.0 1*1.0 1*0.0 15*0.0 11*1.0 1*1.0 1*0.0 15*0.0 11*1.0 1*1.0 1*0.0 15*0.0 11*1.0 1*1.0 1*0.0 15*0.0

盈朕

石不转 发表于 2022-11-3 11:28
您好。FERWE与FERDO的设置可参考VASP官网。官网原话中“The occupancies must be specified for all band ...

谢谢老师，这样计算完激发态构型，结果应该在哪看，用激发态构型进行其他计算，应该用哪个文件的内容呢。一个小请求老师，可以加个联系方式私下交流吗

石不转

盈朕 发表于 2022-11-3 15:57
谢谢老师，这样计算完激发态构型，结果应该在哪看，用激发态构型进行其他计算，应该用哪个文件的内容呢。 ...

您好。以下陋见供您参考，欢迎批评指正：
1）计算完成后可在输出的EIGENVAL文件中查看电子占据情况，如电子占据与FERWE、FERDO设置相同，则改变电子占据成功。
2）本人尚未使用激发态构型进行其他计算。因刚接触不久，尚未明确一步计算后应提取哪些信息进行下一步计算，在此一同求问，期待高手的回复，先行谢过。
3）本人才疏学浅，私下交流恐误人前途，不如公开讨论或可获得论坛内高手指点。望见谅。

盈朕

石不转 发表于 2022-11-3 21:52
您好。以下陋见供您参考，欢迎批评指正：
1）计算完成后可在输出的EIGENVAL文件中查看电子占据情况，如 ...

谢谢您的指点，望此后在此论谈多多向您以及坛内高手请教