计算化学公社

标题: Gaussian或Gaussview如何实现SOMO [打印本页]

作者
Author: diaolanxinyu 时间: 2014-11-25 22:23
标题: Gaussian或Gaussview如何实现SOMO
大家好：
请教大家几个问题；

对于电子贴敷或自由基类，文献中提到“the singly occupied molecular orbitals (SOMOs) of the electron attached states”，
我用gaussview画出的HOMO，LUMO以及spin density图都与文献中的SOMO不一致，我的问题是：
1. SOMO与HOMO,LUMO的区别是什么？
2. SOMO可以用Gamess-US实现，请问Gaussian和view怎么实现呢？

谢谢大家！:)

作者
Author: sobereva 时间: 2014-11-25 22:34
HOMO、LUMO是对于闭壳层体系而言的。SOMO是对于开壳层体系而言的，指单占据的轨道，如图所示
(, 下载次数 Times of downloads: 98)
只要是支持ROHF或RODFT的程序就都能得到SOMO。

比如#P ROB3LYP/6-31G*计算，然后把chk文件转换为fch文件，载入进Multiwfn，进主功能0，文本窗口就会提示哪个是SOMO，然后在图形窗口的轨道列表里点击相应编号即可看到SOMO图形。

作者
Author: diaolanxinyu 时间: 2014-11-25 22:48
谢谢sobereva老师！我这就去试试！:)

作者
Author: diaolanxinyu 时间: 2014-11-25 23:43
sobereva老师您好：

附件fig1是文献中用Gamess-US作出的丙酮的SOMO图，fjg2是我用Multiwfn做出的SOMO图，fig3是用view做出的HOMO图，跟文献中不太一样，是不是因为8.2ev的问题？

作者
Author: diaolanxinyu 时间: 2014-11-25 23:47
这是一篇相关文献。谢谢sobereva老师！

作者
Author: diaolanxinyu 时间: 2014-11-26 00:28
附件是用#p rob3lyp/6-311+g*计算的fchk文件

是不是文献中的SOMO图只能用Gamess-US画了？

谢谢！

作者
Author: sobereva 时间: 2014-11-26 08:21
对于画SOMO，不存在“只能用GAMESS-US”的情况
什么程序用RO算出来的都是SOMO，画出来都是一样的。
Multiwfn和gview画出来的也是完全一样的，只要你把isovalue设为相同（Multiwfn的isovalue比gview默认的大）

和文中的差异可能在于：理论方法不同，基组不同（文中还额外以even-temper方式又加了s和p弥散），等值面数值不同，绘制的SOMO不是一般意义的SOMO。最后一点最有可能，基态下的SOMO能量不可能那么高，我没仔细看他的文章，他所谓的SOMO大抵只是某个空轨道上后来占了个电子后的情况，因此实际上画的是某个虚轨道而已。

作者
Author: diaolanxinyu 时间: 2014-11-26 12:04
非常感谢sobereva老师，您的耐心细致讲解，使我收获很大，我会再仔细阅读文献，并继续学习Multiwfn这个强大的软件！

谢谢您！

作者
Author: ter20 时间: 2014-11-28 11:02

sobereva 发表于 2014-11-25 22:34
HOMO、LUMO是对于闭壳层体系而言的。SOMO是对于开壳层体系而言的，指单占据的轨道，如图所示

只要是支持 ...

现在文献中关于SOMO的用法比较混乱，其实个人认同Sob的观点，只有在RO框架下才有SOMO的概念，HOMO其实是特指双占据。但是现在不少用unrestricted方法做的有关自由基的计算都会用SOMO的概念，其实对应的都是HOMO(α)

作者
Author: zhangyina 时间: 2016-6-28 22:21

sobereva 发表于 2014-11-25 22:34
HOMO、LUMO是对于闭壳层体系而言的。SOMO是对于开壳层体系而言的，指单占据的轨道，如图所示

只要是支持 ...

是使用rob3lyp做单点能计算吗？那优化的过程也必须用rob3lyp来做吗？

作者
Author: sobereva 时间: 2016-6-29 02:09

zhangyina 发表于 2016-6-28 22:21
是使用rob3lyp做单点能计算吗？那优化的过程也必须用rob3lyp来做吗？

是。优化过程用RO还是U都可以，U速度更快结果更接近实际。

作者
Author: guoy14iccas 时间: 2017-6-21 11:35
这篇文章中也是计算含自由基分子的SOMO轨道的。没有SOMO计算经验，想请教各位老师。本文中用的是UCAM-B3LYP/6-31G(d,p)计算的SOMO。然后给出了自旋密度分布。在网上找到说是自旋密度的定义就是alpha电子密度减去beta电子的密度？这个密度怎么算？附件中还给出自旋污染S^^2？这个怎么算？

作者
Author: sobereva 时间: 2017-6-21 13:59

guoy14iccas 发表于 2017-6-21 11:35
这篇文章中也是计算含自由基分子的SOMO轨道的。没有SOMO计算经验，想请教各位老师。本文中用的是UCAM-B3LYP ...

谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353

U是非限制性开壳层计算，此时根本就没有SOMO。SOMO是对RO计算而言的。

作者
Author: guoy14iccas 时间: 2017-6-21 15:27

sobereva 发表于 2017-6-21 13:59
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353

谢谢sob老师。

作者
Author: 虎王 时间: 2021-10-3 15:42

sobereva 发表于 2014-11-25 22:34
HOMO、LUMO是对于闭壳层体系而言的。SOMO是对于开壳层体系而言的，指单占据的轨道，如图所示

只要是支持 ...

除了RO以外，UDFT应该也能得到SOMOs。

作者
Author: sobereva 时间: 2021-10-3 22:07

虎王发表于 2021-10-3 15:42
除了RO以外，UDFT应该也能得到SOMOs。

UKS没有一般意义的SOMO，因为本身就没法说哪个算是singly occupied，尤其是自旋极化显著的时候

顶多是做双正交化之后再说SOMO
用于非限制性开壳层波函数的双正交化方法的原理与应用
http://sobereva.com/448（http://bbs.keinsci.com/thread-11377-1-1.html）

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