隐式溶剂模型通过NELECT加电势的问题

zquan · 发表于 Post on 2022-9-14 20:47:21

我有一个问题想请教一下各位大佬，很多文献在计算固液界面体系时用VASPSOL计算隐式溶剂化模型时，加入背景电荷(NELECT)，可以加上背景电荷。我没太搞明白这个的操作过程，我理解的是，加入背景电荷后，通过VASPSOL计算得到体系的费米能级，把费米能级减去SHE的参考电势(-0.46eV)就可以认为是电极电势，然后拟合一下背景电荷的个数和电极电势的关系就能确定对于已知电势的体系需要加多少NELECT了，不知道我的理解对不对，还请各位大佬赐教

卡开发发 · 发表于 Post on 2022-9-15 06:07:30

1、理论上说应该是μ=-Ef，这是Koopmans近似。但实际上三维周期边界并没有真正意义的真空能级作为参考，因此我们认为远离表面的位置（也就是真空层1/2附近）的电势平稳段为真空能级Evac，此时μ=Evac-Ef。然后再是φ和μ之间的常数关系。
2、φ和q应该未必是完全线性的，可以不断调整q来得到你想要得到的φ。同时应该考虑加入电荷引入的修正，你可以参考FERMI_SHIFT的关键字。
3、上述方法有一定局限性，必须得把上下表面做成对称的。

zquan · 发表于 Post on 2022-9-15 15:40:03

卡开发发发表于 2022-9-15 06:07
1、理论上说应该是μ=-Ef，这是Koopmans近似。但实际上三维周期边界并没有真正意义的真空能级作为参考，因 ...

感谢老师

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] 隐式溶剂模型通过NELECT加电势的问题

浏览过的版块