如何使用cp2k计算aimd体系输出静电势和费米能级

fly1992

程俊老师在PRL 119, 016801 (2017)这篇文献中计算了aimd过程中的电极电势，根据程俊老师的方法：

后面三项都可以按照之前的文献，直接拿到具体的值，需要计算的主要是前两项。
想请教一下各位老师，文献里说第一项指的是界面模型的费米能，第二项指块体水的静电势，并且两个都可以从轨迹文件中直接读取。

这两个能值否在社长的Multiwfn设置inp文件进行输出？费米能能否选择Multiwfn里的-3然后选择-1，最后在output文件里找到？静电势该如何设置？

此外，这两个值能不能对平衡后的某几个时间点的结构进行静态计算得到？
谢谢各位老师？