求助：CP2K使用HLE17泛函计算苯分子吸附体系的能带的若干问题与审稿意见

HydroKagami

老师们好，我使用CP2K在HLE17泛函下计算了某个碳基Fe单原子体系的slab以及其吸附若干有机气体分子的能带。
用Multiwfn解析BAND.bs，slab材料本身的gap为0.009eV，对HCHO、丙酮、苯分子吸附的gap为0.192、0.362、0.392eV。
但是根据吸附能数据，苯分子的吸附能是最差的，并且不同于含氧的分子有明确的结合位点产生化学吸附。
我同时计算了原子电荷（Hirshfeld电荷）显示苯分子的电荷转移量只有0.059e，远小于另外两种有机分子，又通过IGMH分析发现苯分子的确主要以弱相互作用的形式吸附在slab上。
有3个审稿人对这个反常gap打开的机制产生了疑问，并且基于另一个审稿人要求的对PBE泛函计算的测试，我发现PBE计算苯分子吸附体系的能带，其gap只有0.007eV，和PBE计算的纯slab的gap几乎一致，而另外两种有机气体在PBE计算下的带隙仍保持同一趋势。
我本来试图说明是苯分子吸附诱导的slab形变导致的带隙打开，我又补充计算了HLE17下，对苯分子吸附后体系删去苯分子，只对因苯分子吸附导致了形变的slab计算能带，但是这个结果显示gap为0.009eV，也不能说明我的猜想。
我现在还做了对Fe和苯的PDOS和OPDOS分析（HLE17泛函下结合对角化给出的molden），显示Fe的PDOS和苯的PDOS几乎没有重叠，但HOCO的等值面在苯上，LUCO的等值面在slab上。请问我应该如何说明这个带隙打开的机制？还是说HLE17的计算在这个体系中出现了问题，把苯的gap也纳入了？