计算化学公社

标题: spin density and SOMO plot之间的关系问题 [打印本页]

作者
Author: wull 时间: 2019-3-12 12:43
标题: spin density and SOMO plot之间的关系问题
求助
为什么算出的三重态的SOMO和spin density给出的结果完全不相关呢？
是不是这个优化的三重态不稳定不正确，需要stable=opt来获得稳定正确的波函数再优化呢？
谢谢各位啦

作者
Author: sobereva 时间: 2019-3-12 14:47
100%是你的计算有问题，或者对SOMO的理解有问题
我不晓得你怎么算的，如果是用RO方式算的，所有单占据轨道的密度总和直接对应自旋密度

建议参看
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353（http://bbs.keinsci.com/thread-4422-1-1.html）
用于非限制性开壳层波函数的双正交化方法的原理与应用
http://sobereva.com/448（http://bbs.keinsci.com/thread-11377-1-1.html）

对于UHF、UKS方式的计算，又要通过轨道来讨论自旋密度的本质构成，用Multiwfn产生的双正交化轨道讨论是最好选择

作者
Author: wull 时间: 2019-3-12 17:36

sobereva 发表于 2019-3-12 14:47
100%是你的计算有问题，或者对SOMO的理解有问题
我不晓得你怎么算的，如果是用RO方式算的，所有单占据轨道 ...

非常感谢sob老师，原因找到了，是我没有对U计算的结果做双正交化处理。

Multiwfn非常强大！

作者
Author: ljh123 时间: 2024-9-7 21:48

sobereva 发表于 2019-3-12 14:47
100%是你的计算有问题，或者对SOMO的理解有问题
我不晓得你怎么算的，如果是用RO方式算的，所有单占据轨道 ...

sob老师请问一下，对于这个贴的情况，确实如老师所说是因为没有正交化所导致的。但是如果我是想讨论单电子所在的轨道所在，那我可以认为两个单电子在阿尔法轨道的HOMO1和HOMO2轨道吗？而自旋密度这个概念我是否可以理解为仅仅是所有轨道自旋的贡献，而不是单电子所在的位置呢？

作者
Author: sobereva 时间: 2024-9-7 23:14

ljh123 发表于 2024-9-7 21:48
sob老师请问一下，对于这个贴的情况，确实如老师所说是因为没有正交化所导致的。但是如果我是想讨论单电 ...

对于自旋极化明显的情况，你只能说单电子几乎在两个双正交化轨道上（你可以按照双正交化轨道能量进行编号），做了双正交化之后就没有一般意义的HOMO了，HOMO是对于分子轨道而言的

自旋密度直接体现了三维空间中单电子的分布情况，下文解释得很清楚了
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353（http://bbs.keinsci.com/thread-4422-1-1.html）

作者
Author: ljh123 时间: 2024-9-8 00:27

sobereva 发表于 2024-9-7 23:14
对于自旋极化明显的情况，你只能说单电子几乎在两个双正交化轨道上（你可以按照双正交化轨道能量进行编号 ...

谢谢sob老师的回复，那么讨论开壳层体系的单电子分布轨道的时候，可以不用正交化讨论，而是根据阿尔法HOMO轨道形状来认为是单电子占据轨道吗？我有这样的疑问是因为，通过RO方法或者正交化方法得到的单电子占据轨道图和UDFT得到阿尔法HOMO轨道图形状差别非常大。就如该贴贴主的问题（贴主给出的SOMO1和SOMO2应该实际上是UDFT计算得到的2个阿尔法HOMO轨道图，只是文献把这个概念用错了），可以看到贴主给出的SOMO1，SOMO2（实际是阿尔法HOMO1，阿尔法HOMO2)，和自旋密度图差别很大。

作者
Author: sobereva 时间: 2024-9-8 10:58

ljh123 发表于 2024-9-8 00:27
谢谢sob老师的回复，那么讨论开壳层体系的单电子分布轨道的时候，可以不用正交化讨论，而是根据阿尔法HOM ...

“不用正交化讨论” 语义不明
我该说的都说了，仔细把下文认真阅读、理解了
用于非限制性开壳层波函数的双正交化方法的原理与应用
http://sobereva.com/448（http://bbs.keinsci.com/thread-11377-1-1.html）

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