Perodic NBO与vasp计算单原子催化剂体系，NBO无法输出可视化Cube文件

csedulikai

本人已经安装好NBO，是在Gethub上面下载的,vasp也已经编译好了，是6.3.0版本。目前测试过网上的范例Si体系没有问题，能算出可视化轨道，但是在计算Co-N4单原子体系的时候，在运行NBO程序时想要得到cube文件进行作图，但是一直报错：There are a different number of AOs and NBOs in visualization subroutine. Something is wrong   
尝试过修改投影的basis.inp没有效果，也尝试了Generalization of Natural Bond Orbital Analysis to Periodic Systems:Applications to Solids and Surfaces via Plane-Wave DensityFunctional Theory这篇文章删除扰乱的s,p的方法，也没有奏效，以下是我的basis.inp文件
****
Co 0
S   2  1.00
6.50   0.154
1.80   0.535
P   2  1.00
5.20   0.162
1.60   0.570
D   2  1.00
1.20   0.600
0.35   0.400
****
C 0
S   3  1.00
10.00   0.15
2.50   0.55
0.70   0.35
P   2  1.00
1.80   0.60
0.50   0.40
****
N 0
S   3  1.00
12.00   0.15
3.00   0.55
0.80   0.35
P   2  1.00
2.20   0.60
0.60   0.40
****
O 0
S   3  1.00
14.00   0.15
3.50   0.55
1.00   0.35
P   2  1.00
2.50   0.60
0.70   0.40
****

以下是我的nbo.config

#####################################################################
#####Configuration file for the periodic NBO code of JRS and BDD#####
#####################################################################

#####NBO search parameters####
0.80   #Occupancy cutoff for one-center NBOs
0.95   #Occupancy cutoff for two-center NBOs
#####Visualization output control parameters####
     T   #Control over printing of .cube files for visualization.
     F   #density - Whether density (T) or wavefunctions (F) are visualized.
  1 2   #vis_start vis_end - Start and end of NBOs to print .cube files for
40 40 40   #mesh - Number of points along each lattice vectors to use in .cube files
2 2 2   #box_int - Number of unit cell to use for .cube file. See READ_ME.txt for guidance
  -0.50   #origin_fact - Shift of the origin for .cube file. See READ_ME.txt for guidance
以下是我的体系