计算化学公社

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

[ORCA] ORCA和Gaussian中TDDFT三重态计算结果相差大

[复制链接 Copy URL]

39

帖子

1

威望

3984

eV
积分
4043

Level 6 (一方通行)

跳转到指定楼层 Go to specific reply
楼主
各位老师们好,我最近在ORCA和Gaussian中分别用M062X泛函计算了蒽分子的三重激发态能级,但是两个软件的结果相差很大。但是单重态的结果却是相近的,所以比较疑惑,请各位看看是哪里出了问题。谢谢!

An-orca.log

79.52 KB, 下载次数 Times of downloads: 110

An-gaussian.log

45.46 KB, 下载次数 Times of downloads: 40

538

帖子

1

威望

5777

eV
积分
6335

Level 6 (一方通行)

2#
发表于 Post on 2021-1-25 09:15:38 | 只看该作者 Only view this author
我觉得这差别也不至于用“很大”来描述。真要比的话把ORCA的RIJCOSX去掉再比一比。这么小的体系没必要用RIJCOSX。

39

帖子

1

威望

3984

eV
积分
4043

Level 6 (一方通行)

3#
 楼主 Author| 发表于 Post on 2021-1-25 09:51:09 | 只看该作者 Only view this author
niobium 发表于 2021-1-25 09:15
我觉得这差别也不至于用“很大”来描述。真要比的话把ORCA的RIJCOSX去掉再比一比。这么小的体系没必要用RIJ ...

谢谢回复!我去掉RIJCOSX之后重新算了一遍,发现相差还是有点大。主要是Gaussian计算的三重态能量普遍比ORCA的高出接近0.2eV,这我感觉已经不在误差范围内了。

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

4#
发表于 Post on 2021-1-25 22:33:22 | 只看该作者 Only view this author
计算更多的根试试,并且只让Gaussian计算三重态试试
北京科音自然科学研究中心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

39

帖子

1

威望

3984

eV
积分
4043

Level 6 (一方通行)

5#
 楼主 Author| 发表于 Post on 2021-1-26 10:20:51 | 只看该作者 Only view this author
sobereva 发表于 2021-1-25 22:33
计算更多的根试试,并且只让Gaussian计算三重态试试

谢谢sob老师回复!我又计算了20个态,同时让Gaussian只计算三重态,结果还是没有改善。

An-gaussian.log

50.33 KB, 下载次数 Times of downloads: 14

An-orca.log

111.94 KB, 下载次数 Times of downloads: 25

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

6#
发表于 Post on 2021-1-26 19:22:53 | 只看该作者 Only view this author
zhouoh 发表于 2021-1-26 10:20
谢谢sob老师回复!我又计算了20个态,同时让Gaussian只计算三重态,结果还是没有改善。

原因不好说。换个泛函比如PBE0,或者都用TDA再比比试试。始终不开RI
北京科音自然科学研究中心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

26

帖子

0

威望

339

eV
积分
365

Level 3 能力者

7#
发表于 Post on 2022-2-9 22:20:16 | 只看该作者 Only view this author
想问下最后题主问题解决了吗,遇到了相同的问题

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

8#
发表于 Post on 2022-2-10 00:09:45 | 只看该作者 Only view this author
huahua0621 发表于 2022-2-9 15:20
想问下最后题主问题解决了吗,遇到了相同的问题

有可能是ORCA 4.2.1的bug,我用ORCA 5.0.2测楼主的输入文件,结果只和高斯差0.003eV。
可以下载ORCA 5.0.2试试,或者等几天ORCA 5.0.3出来以后下载5.0.3
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

26

帖子

0

威望

339

eV
积分
365

Level 3 能力者

9#
发表于 Post on 2022-2-10 22:16:24 | 只看该作者 Only view this author
wzkchem5 发表于 2022-2-10 00:09
有可能是ORCA 4.2.1的bug,我用ORCA 5.0.2测楼主的输入文件,结果只和高斯差0.003eV。
可以下载ORCA 5.0 ...

我用的是5.0.2,和singlet一样的输入只是改了MULT和坐标前的多重度,T1计算的结果甚至比S1的绝热能都高 算的QD-NEVPT2

4291

帖子

4

威望

9565

eV
积分
13936

Level 6 (一方通行)

MOKIT开发者

10#
发表于 Post on 2022-2-10 22:42:37 | 只看该作者 Only view this author
huahua0621 发表于 2022-2-10 22:16
我用的是5.0.2,和singlet一样的输入只是改了MULT和坐标前的多重度,T1计算的结果甚至比S1的绝热能都高  ...

你这问题与贴主所问完全不是同一个问题,你的问题复杂多了,而且高斯也没有NEVPT2给你对比,不能用“遇到了相同的问题”来描述。

QD-NEVPT2要讲出你的计算过程(计算流程,例如先算啥,再算啥,怎么选活性轨道,SA-CASSCF做完之后是否轨道符合自己的期望?SA-CASSCF用了几个根做平均?QD用了几个态?所用基组?),你看这细节是TDDFT的四五倍,你一个都没讲
自动做多参考态计算的程序MOKIT

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

11#
发表于 Post on 2022-2-11 02:21:10 | 只看该作者 Only view this author
huahua0621 发表于 2022-2-10 15:16
我用的是5.0.2,和singlet一样的输入只是改了MULT和坐标前的多重度,T1计算的结果甚至比S1的绝热能都高  ...

T1比S1能量高的分子是存在的(虽然极其罕见),算出这样的结果并不一定代表算错了。我记得听过Aspuru-Guzik讲过这种体系的计算,你可以搜一下他的文章
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

26

帖子

0

威望

339

eV
积分
365

Level 3 能力者

12#
发表于 Post on 2022-2-16 20:47:07 | 只看该作者 Only view this author
zjxitcc 发表于 2022-2-10 22:42
你这问题与贴主所问完全不是同一个问题,你的问题复杂多了,而且高斯也没有NEVPT2给你对比,不能用“遇到 ...

DFT优化,QD-NEVPT2做单点,用DFT的轨道做的初猜,先调的CAS轨道,15个根做平均,然后算的QD-NEVPT2,吸收光谱基本和实验吻合,读了单重态的轨道去算的三重态,QD也是15个根,基组是def2-svp 对Ir用的SDD 算三重态也用的15个根,轨道也基本和之前一致,就是能量差距很大

4291

帖子

4

威望

9565

eV
积分
13936

Level 6 (一方通行)

MOKIT开发者

13#
发表于 Post on 2022-2-16 21:06:31 | 只看该作者 Only view this author
huahua0621 发表于 2022-2-16 20:47
DFT优化,QD-NEVPT2做单点,用DFT的轨道做的初猜,先调的CAS轨道,15个根做平均,然后算的QD-NEVPT2,吸收 ...

你可以另开一贴,晒出你的分子,你的计算结果,你期望的结果(或实验值),以便详细分析、讨论
自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

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

GMT+8, 2026-2-27 20:15 , Processed in 0.172272 second(s), 23 queries , Gzip On.

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