应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 怎么用rpa算吸附能？
31/1返回列表
查看: 1014  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

1455878431

新虫 (小有名气)

[求助] 怎么用rpa算吸附能？ 已有1人参与

请问各位大佬，如何用vasp的rpa算吸附能？能给出incar以及相关参数吗

发自小木虫IOS客户端
回复此楼

» 猜你喜欢
1楼 2019-11-05 15:40:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

123dingyc
rpa

发自小木虫IOS客户端
2楼2019-11-05 15:44:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zmyxmc

新虫 (初入文坛)

【答案】应助回帖

手册里有讨论：
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1455878431 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭