| 查看: 994 | 回复: 2 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
1455878431新虫 (小有名气)
|
[求助]
怎么用rpa算吸附能？ 已有1人参与
|
||
|
请问各位大佬，如何用vasp的rpa算吸附能？能给出incar以及相关参数吗 发自小木虫IOS客户端 |
回复此楼» 猜你喜欢
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有14人回复
-
重要来源：本周末出结果
已经有7人回复
-
结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。
已经有4人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有11人回复
-
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
filecode
已经有15人回复
【答案】应助回帖
|
手册里有讨论： https://cms.mpi.univie.ac.at/vasp/vasp/general_recipe_calculate_ACFDT_RPA_total_energies.html 简而言之，对半导体体系、保持k点和encut不变前提下是四步 1. 正常semilocal泛函下作SCF，得到WAVECAR 2. 拿1的WAVECAR，作单步非自洽的HF计算，得到总能量为EHF 3. 拿1的WAVECAR作为输入，作一步ALGO=exact精确对角化得到包含全部空轨道的WAVECAR，同时打开LOPTICS得到WAVEDER。后者不是必须，但本身体量也小，能加快RPA关联对k点的收敛。参考Harl等2010年的PRB。 4. 拿3算出来的WAVECAR和WAVEDER作为输入，作ALGO=acfdt，算出RPA关联能，EcRPA 最后RPA能量是EHF+EcRPA。算A在B上的吸附能无非是，E(ads) = ERPA(A+B) - ERPA(A) - ERPA(B)，跟DFT一样只不过换成了RPA。我没做过slab的，但如果你拿比较大的slab或k点算的话计算量会比较可怕。 |
3楼2019-11-15 10:52:11
|
rpa 发自小木虫IOS客户端 |
2楼2019-11-05 15:44:27
