| 查看: 841 | 回复: 2 | |||
1455878431新虫 (小有名气)
|
[求助]
怎么用rpa算吸附能? 已有1人参与
|
|
请问各位大佬,如何用vasp的rpa算吸附能?能给出incar以及相关参数吗 发自小木虫IOS客户端 |
» 猜你喜欢
材料专业求调剂
已经有5人回复
一志愿天津大学化学工艺专业(081702)315分求调剂
已经有10人回复
344求调剂
已经有7人回复
一志愿985,本科211,0817化学工程与技术319求调剂
已经有7人回复
271材料工程求调剂
已经有4人回复
材料专硕英一数二306
已经有4人回复
294求调剂材料与化工专硕
已经有4人回复
一志愿中海洋材料工程专硕330分求调剂
已经有4人回复
070303一志愿西北大学学硕310找调剂
已经有9人回复
0703化学调剂,求各位老师收留
已经有9人回复
|
rpa 发自小木虫IOS客户端 |
2楼2019-11-05 15:44:27
【答案】应助回帖
|
手册里有讨论: 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













回复此楼