计算化学公社

标题: Multiwfn 输出定域化轨道的cube文件经计算后轨道电子数非0,1,2 [打印本页]

作者
Author: 赵小壮 时间: 2024-3-26 09:48
标题: Multiwfn 输出定域化轨道的cube文件经计算后轨道电子数非0,1,2
本帖最后由赵小壮于 2024-3-26 09:50 编辑

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

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

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

复制代码

作者
Author: sobereva 时间: 2024-3-26 14:53
单个原子根本没有轨道定域化的事，毫无意义

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)