最近试图用siesta求解金属的功函数。阅读了一些文献，很多文献都是用一种名为“bulk plus band lineup method”的方法计算的，包括siesta的教程ppt也指明了这种方法，主要是用来计算VBO和CBO(band offset)，当然也可以用来计算金属的功函数。此种方法包含两项，一个是bulk项：对金属bulk进行计算后获得其planar average electrostatic potential 和Fermi level，并求其差；另一项是针对slab相的计算，针对slab的计算时为了获得其真空能级和slab中的静电势。我的理解是这两项算完以后，将bulk的静电势和slab里面的静电势放在同一个水平线上，便获得了真空能级和费米能级的相对位置，求差便可得到金属的功函数。主要有如下个问题：
1.上述做法是否正确？是否体现了bulk plus band lineup method 的思想？
2.siesta 计算完成后其out文件中有一栏能量项（如下），最后的Fermi= xxx 是不是指的费米能级；如果不是那么feimi能级具体该如何计算呢？
siesta: Final energy (eV):
siesta: Band Struct. = -1620.094511
siesta: Kinetic = 43092.791084
siesta: Hartree = 9082.228050
siesta: Eldau = 0.000000
siesta: Eso = 0.000000
siesta: Ext. field = 0.000000
siesta: Exch.-corr. = -6273.012046
siesta: Ion-electron = -54736.661070
siesta: Ion-ion = -29487.373850
siesta: Ekinion = 0.000000
siesta: Total = -38322.027832
siesta: Fermi = -1.699580
3.我对Pb和Cu进行了计算，并且画出了bulk和slab的静电势能曲线，但是按照上述方法得不到正确的值，和实验值相差十万八千里，不知道是哪里出了问题。还请各位老师不吝赐教
图像是在输入fdf文件中添加 SaveElectrostaticPotential T 后产生xxx.VH文件，再利用macroave工具对其进行处理得到.PAV文件和.MAV文件后在Excel中绘制得到的。