计算化学公社

 找回密码 Forget password
 注册 Register
Views: 408|回复 Reply: 13
打印 Print 上一主题 Last thread 下一主题 Next thread

[波函数分析求助] 使用Multiwfn处理Gaussian 和 ORCA 计算得到同一个分子的 Chole Cele差异很大

[复制链接 Copy URL]

65

帖子

1

威望

535

eV
积分
620

Level 4 (黑子)

本帖最后由 qinsc 于 2025-6-22 14:25 编辑

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

Gaussian结果如下


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

分子结构如下
structure.xyz (9.63 KB, 下载次数 Times of downloads: 2)
ORCA inp 文件如下
wb97xdtddft.inp (125 Bytes, 下载次数 Times of downloads: 4)
Gaussian 输入文件如下
structure.gjf (8.97 KB, 下载次数 Times of downloads: 6)

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

2#
发表于 Post on 2025-6-22 09:08:47 | 只看该作者 Only view this author
ORCA和Gaussian根本没这个计算功能,明明是用Multiwfn按照下文算的。描述问题时不要有误导性
使用Multiwfn做空穴-电子分析全面考察电子激发特征
http://sobereva.com/434http://bbs.keinsci.com/thread-10775-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

410

帖子

5

威望

1630

eV
积分
2140

Level 5 (御坂)

鸩羽

3#
发表于 Post on 2025-6-22 11:15:10 | 只看该作者 Only view this author
先把打印阈值统一
某不知名实验组从苞米地里长出来的计算选手

65

帖子

1

威望

535

eV
积分
620

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2025-6-22 14:23:12 | 只看该作者 Only view this author
sobereva 发表于 2025-6-22 09:08
ORCA和Gaussian根本没这个计算功能,明明是用Multiwfn按照下文算的。描述问题时不要有误导性
使用Multiwfn ...

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

65

帖子

1

威望

535

eV
积分
620

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2025-6-22 14:45:30 | 只看该作者 Only view this author
wal 发表于 2025-6-22 11:15
先把打印阈值统一

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

410

帖子

5

威望

1630

eV
积分
2140

Level 5 (御坂)

鸩羽

6#
发表于 Post on 2025-6-22 17:51:40 | 只看该作者 Only view this author
qinsc 发表于 2025-6-22 14:45
谢谢老师,高斯在设定IOp(9/40=4)后结果还是一样,这里主要还是Multiwfn处理ORCA的输出文件得到的结果不 ...

不太清楚,不过我看手册感觉cholecele只是平滑化处理了一下
如果更改打印阈值不能产生可见的变化,也可能是tda近似引起的
某不知名实验组从苞米地里长出来的计算选手

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

7#
发表于 Post on 2025-6-23 05:01:40 | 只看该作者 Only view this author
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相差并不仅限于色散校正部分
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

65

帖子

1

威望

535

eV
积分
620

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2025-6-23 12:57:19 | 只看该作者 Only view this author
wal 发表于 2025-6-22 17:51
不太清楚,不过我看手册感觉cholecele只是平滑化处理了一下
如果更改打印阈值不能产生可见的变化,也可 ...

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

65

帖子

1

威望

535

eV
积分
620

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2025-6-23 14:13:16 | 只看该作者 Only view this author
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是波函数的问题还是后续处理的问题。

65

帖子

1

威望

535

eV
积分
620

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2025-6-23 14:14:04 | 只看该作者 Only view this author
wal 发表于 2025-6-22 17:51
不太清楚,不过我看手册感觉cholecele只是平滑化处理了一下
如果更改打印阈值不能产生可见的变化,也可 ...

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

410

帖子

5

威望

1630

eV
积分
2140

Level 5 (御坂)

鸩羽

11#
发表于 Post on 2025-6-23 14:26:18 | 只看该作者 Only view this author
qinsc 发表于 2025-6-23 14:14
老师,我用Gaussian的TDA也试了一下,结果和TD的一样,我感觉应该不是TDA造成的。

除了泛函,TDA,打印阈值之外,我能看到的差异就剩ORCA默认用的RI近似了。
你可以关掉RI试试,如果还是不能跟Gaussian对应上,那就等sob老师解答吧
某不知名实验组从苞米地里长出来的计算选手

65

帖子

1

威望

535

eV
积分
620

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2025-6-23 18:28:12 | 只看该作者 Only view this author
wal 发表于 2025-6-23 14:26
除了泛函,TDA,打印阈值之外,我能看到的差异就剩ORCA默认用的RI近似了。
你可以关掉RI试试,如果还是 ...

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

65

帖子

1

威望

535

eV
积分
620

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2025-6-23 22:41:49 | 只看该作者 Only view this author
wal 发表于 2025-6-23 14:26
除了泛函,TDA,打印阈值之外,我能看到的差异就剩ORCA默认用的RI近似了。
你可以关掉RI试试,如果还是 ...

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

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

14#
发表于 Post on 2025-6-24 05:06:58 | 只看该作者 Only view this author
qinsc 发表于 2025-6-23 22:41
关掉RI还是不行,我再研究一下

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

如果死活弄不清楚,把两个程序计算设置尽可能全都一致情况下的Gaussian的输入输出文件、fch文件,ORCA的输入输出文件、molden文件,全都压缩传到网盘上贴链接
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-12 21:00 , Processed in 0.218635 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list