计算化学公社

标题: 使用Multiwfn处理Gaussian 和 ORCA 计算得到同一个分子的 Chole Cele差异很大 [打印本页]

作者
Author: qinsc 时间: 2025-6-22 08:40
标题: 使用Multiwfn处理Gaussian 和 ORCA 计算得到同一个分子的 Chole Cele差异很大
本帖最后由 qinsc 于 2025-6-22 14:25 编辑

各位老师好，我使用Multiwfn的电子空穴分析功能处理ORCA （6.0.1）计算计算结果，得到一个分子的 Chole Cele，发现形状有些奇怪，和之前的类似分子结果不一样，遂又使用Gaussian（16）算了一下，Multiwfn处理可以得到正常的图像，想问一下我在使用ORCA的时候是哪里没注意到吗？
ORCA结果如下
(, 下载次数 Times of downloads: 75) (, 下载次数 Times of downloads: 70)

Gaussian结果如下
(, 下载次数 Times of downloads: 75) (, 下载次数 Times of downloads: 70)

可以看到gaussian的结果更合理。

分子结构如下
(, 下载次数 Times of downloads: 2)
ORCA inp 文件如下
(, 下载次数 Times of downloads: 4)
Gaussian 输入文件如下
(, 下载次数 Times of downloads: 6)

作者
Author: sobereva 时间: 2025-6-22 09:08
ORCA和Gaussian根本没这个计算功能，明明是用Multiwfn按照下文算的。描述问题时不要有误导性
使用Multiwfn做空穴-电子分析全面考察电子激发特征
http://sobereva.com/434（http://bbs.keinsci.com/thread-10775-1-1.html）

作者
Author: wal 时间: 2025-6-22 11:15
先把打印阈值统一

作者
Author: qinsc 时间: 2025-6-22 14:23

sobereva 发表于 2025-6-22 09:08
ORCA和Gaussian根本没这个计算功能，明明是用Multiwfn按照下文算的。描述问题时不要有误导性
使用Multiwfn ...

卢老师已经修改描述，之前说的不够清晰。

作者
Author: qinsc 时间: 2025-6-22 14:45

wal 发表于 2025-6-22 11:15
先把打印阈值统一

谢谢老师，高斯在设定IOp(9/40=4)后结果还是一样，这里主要还是Multiwfn处理ORCA的输出文件得到的结果不太符合预期，我也试过更改tprint的值，结果都一样。
奇怪的一点是除了CholeCele的图像外，电子空穴分布的图像是非常接近的，所以是不是CholeCele所依赖的一些输出ORCA 有些问题呢？

作者
Author: wal 时间: 2025-6-22 17:51

qinsc 发表于 2025-6-22 14:45
谢谢老师，高斯在设定IOp(9/40=4)后结果还是一样，这里主要还是Multiwfn处理ORCA的输出文件得到的结果不 ...

不太清楚，不过我看手册感觉cholecele只是平滑化处理了一下
如果更改打印阈值不能产生可见的变化，也可能是tda近似引起的

作者
Author: sobereva 时间: 2025-6-23 05:01

qinsc 发表于 2025-6-22 14:45
谢谢老师，高斯在设定IOp(9/40=4)后结果还是一样，这里主要还是Multiwfn处理ORCA的输出文件得到的结果不 ...

ORCA默认用TDA，要用Gaussian与之对比先把TDA还是TD统一了再说。并且仔细阅读Multiwfn手册3.21.A节ORCA users:下面IMPORTANT NOTE:部分的文字
并且注意wB97X-D3和wB97XD相差并不仅限于色散校正部分

作者
Author: qinsc 时间: 2025-6-23 12:57

wal 发表于 2025-6-22 17:51
不太清楚，不过我看手册感觉cholecele只是平滑化处理了一下
如果更改打印阈值不能产生可见的变化，也可 ...

好的，谢谢老师，我再研究一下TDA和TD

作者
Author: qinsc 时间: 2025-6-23 14:13

sobereva 发表于 2025-6-23 05:01
ORCA默认用TDA，要用Gaussian与之对比先把TDA还是TD统一了再说。并且仔细阅读Multiwfn手册3.21.A节ORCA u ...

谢谢卢老师回答，我看到了 IMPORTANT NOTE，所以ORCA在计算时使用了默认的TDA（但TDA FALSE我也试过，结果是一样的）。
刚刚验证了Gaussian使用TDA，结果仍是符合预期的。
我也使用过其他泛函比如CAM-B3LYP进行对比，结果还是Gaussian更合理，ORCA不符合预期，两者差别很大。
另外Gaussian和ORCA的结果经过Multiwfn处理后空穴电子分布图是几乎相同的，所以不确定CholeCele是波函数的问题还是后续处理的问题。

作者
Author: qinsc 时间: 2025-6-23 14:14

wal 发表于 2025-6-22 17:51
不太清楚，不过我看手册感觉cholecele只是平滑化处理了一下
如果更改打印阈值不能产生可见的变化，也可 ...

老师，我用Gaussian的TDA也试了一下，结果和TD的一样，我感觉应该不是TDA造成的。

作者
Author: wal 时间: 2025-6-23 14:26

qinsc 发表于 2025-6-23 14:14
老师，我用Gaussian的TDA也试了一下，结果和TD的一样，我感觉应该不是TDA造成的。

除了泛函，TDA，打印阈值之外，我能看到的差异就剩ORCA默认用的RI近似了。
你可以关掉RI试试，如果还是不能跟Gaussian对应上，那就等sob老师解答吧

作者
Author: qinsc 时间: 2025-6-23 18:28

wal 发表于 2025-6-23 14:26
除了泛函，TDA，打印阈值之外，我能看到的差异就剩ORCA默认用的RI近似了。
你可以关掉RI试试，如果还是 ...

还真是，这个给忘了，谢谢老师提醒，我再试试。

作者
Author: qinsc 时间: 2025-6-23 22:41

wal 发表于 2025-6-23 14:26
除了泛函，TDA，打印阈值之外，我能看到的差异就剩ORCA默认用的RI近似了。
你可以关掉RI试试，如果还是 ...

关掉RI还是不行，我再研究一下

作者
Author: sobereva 时间: 2025-6-24 05:06

qinsc 发表于 2025-6-23 22:41
关掉RI还是不行，我再研究一下

对比Gaussian和ORCA计算得到的波函数文件里的轨道（看等值面图）、组态系数（可以用《使用Multiwfn便利地查看所有激发态中的主要轨道跃迁贡献》http://sobereva.com/529的做法对比）
记住ORCA默认用的是球谐型基函数，注意看下文。能统一的都尽可能统一再比
谈谈不同量子化学程序计算结果的差异问题
http://sobereva.com/573（http://bbs.keinsci.com/thread-20122-1-1.html）

如果死活弄不清楚，把两个程序计算设置尽可能全都一致情况下的Gaussian的输入输出文件、fch文件，ORCA的输入输出文件、molden文件，全都压缩传到网盘上贴链接

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