计算化学公社

标题: 为什么在Multiwfn中计算自旋布居时，Alpha pop.的值出现了负值，怎么理解？ [打印本页]

作者
Author: yinbaoqi 时间: 2019-11-25 19:19
标题: 为什么在Multiwfn中计算自旋布居时，Alpha pop.的值出现了负值，怎么理解？
如图，请问有人遇到这种情况吗?该如何理解

作者
Author: sobereva 时间: 2019-11-25 19:45
如果你用了弥散函数，Mulliken或其变体算的轨道成份以及布居数就可能出现负值，这是方法自身的缺陷。用这类方法时绝对不能带弥散函数
如果是其它情况，上传量化程序的输入文件以及给Multiwfn用的输入文件

作者
Author: yinbaoqi 时间: 2019-11-25 20:13
谢谢社长！确实用了弥散函数，但我看对自旋密度的分布图好像没有影响，因为含有非金属元素就用了6-311G*。

作者
Author: sobereva 时间: 2019-11-26 23:40

yinbaoqi 发表于 2019-11-25 20:13
谢谢社长！确实用了弥散函数，但我看对自旋密度的分布图好像没有影响，因为含有非金属元素就用了6-311G*。

自旋密度是实空间函数，因此不怕弥散函数。弥散函数对于多数实空间函数的分布特征影响也比较小

作者
Author: yinbaoqi 时间: 2019-11-27 15:40
本帖最后由 yinbaoqi 于 2019-11-27 16:05 编辑

sobereva 发表于 2019-11-26 23:40
自旋密度是实空间函数，因此不怕弥散函数。弥散函数对于多数实空间函数的分布特征影响也比较小

太感谢您的回复！我尝试了SCPA的方法好像不会涉及弥散函数的这个问题，我试图用SCPA的方法分析自旋布局，金属用的是HSE06/SDD，但我发现用Mulliken和SCPA做布局分析的时候得到的结果差别很大，如图Cu4的自旋布局0.43uB的是用SCPA方法，而Mulliken方法得到的Cu4只有0.14uB。体系只有一个单电子的我应该用哪种方法更可信？

作者
Author: sobereva 时间: 2019-11-28 02:20

yinbaoqi 发表于 2019-11-27 15:40
太感谢您的回复！我尝试了SCPA的方法好像不会涉及弥散函数的这个问题，我试图用SCPA的方法分析自旋布局， ...

有弥散函数下SCPA结果同样无意义
此文对原因有明确体现
分子轨道成分的计算（http://sioc-journal.cn/Jwk_hxxb/CN/abstract/abstract340458.shtml）

必须有弥散函数的情况，应当用Hirshfeld或Becke方法计算

作者
Author: yinbaoqi 时间: 2019-11-28 09:49

sobereva 发表于 2019-11-28 02:20
有弥散函数下SCPA结果同样无意义
此文对原因有明确体现
分子轨道成分的计算（http://sioc-journal.cn/J ...

谢谢社长！我去尝试一下

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