Multiwfn 输出定域化轨道的cube文件经计算后轨道电子数非0,1,2

赵小壮

大家好，我在计算一个第四周期单原子的波函数文件后，输出该原子1s轨道的电子密度的cube文件，然后计算该轨道上的电子数目时，发现1s轨道上电子数目有15.3个

附件是我的轨道文件，我在定域化轨道后通过主功能6 26确实把其它轨道占据数清零了，然后我的check cube文件电子数目的code是没有问题的，用过很多次了。
并且以上过程和code用在第二周期时候都是正常的，算出来的电子数目也都是1.9999这个级别。但是算第三周期和第四周期时候就不对了。请问这是怎么一回事

1. 
   
2. Multiwfn one.molden.input << EOF > out.txt
   
3. 6
   
4. 26
   
5. 2,3,4,5,6,7,8,9,10,11
   
6. 0.0
   
7. q
   
8. -1
   
9. 5
   
10. 1
    
11. 3
    
12. 2
    
13. EOF
    
14. 
    
15. mv density.cub density_1s.cub
    

复制代码

sobereva

单个原子根本没有轨道定域化的事，毫无意义

周期表越靠后的原子，原子核附近内核轨道密度变化越剧烈，越不可能靠均匀格点准确积分，哪怕格点间距设得非常小。要积分这种情况应当用Multiwfn主功能100的子功能4通过下文介绍的原子中心格点来积分
密度泛函计算中的格点积分方法
http://sobereva.com/69