利用Siesta求解金属功函数遇到若干问题

Puying

最近试图用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中绘制得到的。

Puying

做这种工作还是用基于平面波的第一性原理程序好。。。。不要用siesta