计算化学公社

标题: Gauss计算结构分子开壳层体系判别？ [打印本页]

作者
Author: lao7 时间: 2017-4-20 12:28
标题: Gauss计算结构分子开壳层体系判别？
本帖最后由 lao7 于 2017-4-20 12:35 编辑

读了《使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图》一文，文中提到对于CDA的计算，开壳层和闭壳层体系处理方法不同，问题是：1）如何根据高斯计算结果，判定分子体系是开壳层体系还是闭壳层体系呢？2）文中提到需要对复合物，分子（A和B）分别进行再计算，才能进行CDA处理。那么如果计算方法没有用到HF计算，是否闭壳层需要pop=NO。仅仅写成，例如：“B3LYP/6-31G** Opt Freq”即可。

作者
Author: 冰释之川 时间: 2017-4-20 14:08
(1) 看你的输入坐标里的自旋多重度设置，自旋多重度为1是，一般是闭壳层单重态（当然如果是自旋极化单重态的话，则是开壳层的）
(2)没用到post-HF就直接用复合物和片段的单点fchk文件就可以了。你想要折腾，就用它们的输出文件*.out,不过此时对片段计算必须使用pop=full关键词，对复合物必须用pop=full IOp(3/33=1)关键词。因为Multiwfn要从复合物的输出文件中读取以原子中心高斯函数为基的重叠矩阵，这是要写IOp(3/33=1)的原因。由于Multiwfn还要读取复合物及片段的轨道系数，所以必须写pop=full。对于你优化过的复合物，做单点就是“B3LYP/6-31G**”完事了。

作者
Author: sobereva 时间: 2017-4-20 19:51
不是根据高斯计算结果判断体系是开壳层还是闭壳层，而是在高斯计算之前你就理应知道是开壳层还是闭壳层，懂得基本化学和理论化学常识就很容易判断。

不做后HF计算显然不需要pop=NO。把博文里的例子完整做一遍就知道了。

作者
Author: lao7 时间: 2017-4-20 20:53

sobereva 发表于 2017-4-20 19:51
不是根据高斯计算结果判断体系是开壳层还是闭壳层，而是在高斯计算之前你就理应知道是开壳层还是闭壳层，懂 ...

博文里面的几个例子的附件，怎么无法下载呢？谢谢！

作者
Author: sobereva 时间: 2017-4-21 04:55

lao7 发表于 2017-4-20 20:53
博文里面的几个例子的附件，怎么无法下载呢？谢谢！

博文里的例子文件都是Multiwfn自带的。其它链接也都正常。

作者
Author: lao7 时间: 2017-4-22 21:00

sobereva 发表于 2017-4-21 04:55
博文里的例子文件都是Multiwfn自带的。其它链接也都正常。

谢谢！我演练一遍。

作者
Author: Somnus荣小荣 时间: 2020-4-16 15:40

sobereva 发表于 2017-4-20 19:51
不是根据高斯计算结果判断体系是开壳层还是闭壳层，而是在高斯计算之前你就理应知道是开壳层还是闭壳层，懂 ...

老师，请问一下，如果是计算之前是在单态开壳层下进行，过程中会跑到单态闭壳层下吗，

作者
Author: sobereva 时间: 2020-4-19 00:47

Somnus荣小荣发表于 2020-4-16 15:40
老师，请问一下，如果是计算之前是在单态开壳层下进行，过程中会跑到单态闭壳层下吗，

如果某个结构下基态电子结构本来就是闭壳层的，即便是U，也会自然而然得到和闭壳层计算相同的结果

作者
Author: maxuexiang 时间: 2020-11-2 16:31

sobereva 发表于 2020-4-19 00:47
如果某个结构下基态电子结构本来就是闭壳层的，即便是U，也会自然而然得到和闭壳层计算相同的结果

sob老师，有没有可能某个结构根本不存在开壳层单重态呢

作者
Author: zjxitcc 时间: 2020-11-2 16:54

maxuexiang 发表于 2020-11-2 16:31
sob老师，有没有可能某个结构根本不存在开壳层单重态呢

无需用词“某个”。大部分的有机分子基态都是闭壳层单重态的。

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