计算化学公社

标题: 准分子密度计算问题 [打印本页]

作者
Author: masgils 时间: 2023-7-7 15:40
标题: 准分子密度计算问题
各位老师好：

我使用Multiwfn读取了 examples/atomwfn/C .wfn 作为输入，计算了坐标 0,0,0 处的电子密度（原子核），得到Density of all electrons结果是0.1169734707E+03
操作步骤为：
Multiwfn "C .wfn"
1 # 主功能1 Output all properties at a point
0,0,0 # 原子核坐标
2 # 单位选择了Angstrom

然后我以只含有一个碳原子的pdb文件作为输入，计算原子核处的电子密度，由于此时载入的是pdb文件，计算得到的是准分子密度，得到Density of electrons结果是0.1213720000E+03
操作步骤为：
Multiwfn test.pdb
1 # 主功能1 Output all properties at a point
-10.992,-8.796,60.617 # 原子核坐标
2 # 单位选择了Angstrom

我理解的Multiwfn准分子密度计算是先读取了examples/atomwfn目录中原子wfn文件（这里的wfn应该是做过球对称化处理之后的？），将GTF的中心放置到待计算结构的原子核处进行电子密度的计算。所以我预期直接用单原子wfn文件算电子密度和用单原子pdb算准分子密度在原子核处应该得到一样的结果，想请问下我得到不同结果的原因？

我用的Multiwfn版本为 Version 3.8(dev), release date: 2022-Dec-18
附件是我使用的2个输入文件，以及Multiwfn在屏幕上打印出的结果

作者
Author: sobereva 时间: 2023-7-8 10:21
第二种情况Multiwfn屏幕上没提示examples/atomwfn里的wfn文件，自然跟它们毫无关系
此时准分子密度是Multiwfn内置的球对称化的自由原子密度叠加产生的，这些密度在Multiwfn手册附录3 Detail of built-in atomic densities节里写明了

作者
Author: masgils 时间: 2023-7-8 11:03
谢谢社长解答！

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