计算化学公社

标题: 隐式溶剂模型通过NELECT加电势的问题 [打印本页]

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

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

作者
Author: zquan 时间: 2022-9-15 15:40

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

感谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)