计算化学公社

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

[Gaussian/gview] 求助gaussian 16中计算能量转移EET得到负值的问题和计算公式

[复制链接 Copy URL]

21

帖子

0

威望

781

eV
积分
802

Level 4 (黑子)

使用gaussian 16中的EET部分计算了激发态给受体分子之间的电子耦合,请问为什么会有负值出现呢?查阅了qchem中对VET耦合值的描述公式如下图,请问高斯中的计算公式是什么?
两种计算得到的耦合值区别是什么?




202212010924266619..png (45.21 KB, 下载次数 Times of downloads: 11)

202212010924266619..png

202211302034265437..png (19.73 KB, 下载次数 Times of downloads: 10)

202211302034265437..png

5万

帖子

99

威望

5万

eV
积分
112393

管理员

公社社长

2#
发表于 Post on 2022-12-1 09:33:53 | 只看该作者 Only view this author
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,并清楚、准确反映出帖子具体内容,避免有任何歧义和含糊性,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “gaussian 16中计算能量转移EET” 改了,以后务必注意
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

5万

帖子

99

威望

5万

eV
积分
112393

管理员

公社社长

3#
发表于 Post on 2022-12-1 09:35:54 | 只看该作者 Only view this author
相关公式看DOI: 10.1002/wcms.1392
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

89

帖子

0

威望

1713

eV
积分
1802

Level 5 (御坂)

4#
发表于 Post on 2022-12-1 10:06:24 | 只看该作者 Only view this author
gaussian给出的J的符号不一定对,这个需要自己再判断一下
当然,一般用这个算速率的话就不用考虑这个问题

21

帖子

0

威望

781

eV
积分
802

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-12-1 10:46:51 | 只看该作者 Only view this author
sobereva 发表于 2022-12-1 09:35
相关公式看DOI: 10.1002/wcms.1392

谢谢老师的回复,我以后一定注意!

21

帖子

0

威望

781

eV
积分
802

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2022-12-1 11:04:10 | 只看该作者 Only view this author
本帖最后由 liushu 于 2022-12-1 11:05 编辑
风起~ 发表于 2022-12-1 10:06
gaussian给出的J的符号不一定对,这个需要自己再判断一下
当然,一般用这个算速率的话就不用考虑这个问题

老师您好,我在10.1002/ange.202113425的支持文献中看到如下关于能量转移速率的公式,想模仿计算kET来证明主客体间能够发生TTET过程。请问您说的gaussian计算J光谱重叠积分的不准确性,会使得该公式下的VET也会有不准确性吗?

energy transfer.png (200.48 KB, 下载次数 Times of downloads: 10)

energy transfer.png

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

7#
发表于 Post on 2022-12-1 15:17:47 | 只看该作者 Only view this author
电子耦合矩阵元本来就不一定是正的,因为和两个态的相对相位有关。之所以有的程序或者文章只报正的值,是因为他们要么取了绝对值,要么通过适当选取两个态的相对相位来使得电子耦合矩阵元恒非负。Q-Chem属于后一种做法,看一下你截图里这几个公式的推导就懂了。如果只是为了计算能量转移速率的话,只有电子耦合矩阵元的模平方是重要的,所以符号无所谓。
不过我不确定”电子耦合“和”电子耦合矩阵元“是不是一个概念,还是说前者是后者的绝对值。所以不好说”电子耦合“是否应当出现负值,这个还请专业人士补充。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

89

帖子

0

威望

1713

eV
积分
1802

Level 5 (御坂)

8#
发表于 Post on 2022-12-1 15:33:28 | 只看该作者 Only view this author
liushu 发表于 2022-12-1 11:04
老师您好,我在10.1002/ange.202113425的支持文献中看到如下关于能量转移速率的公式,想模仿计算kET来证 ...

我说的J就是你文章中的V_{ET},用marcus公式算的话,不用管这个符号
关于符号的定义,可以看看Spano的综述,

101

帖子

0

威望

2095

eV
积分
2196

Level 5 (御坂)

9#
发表于 Post on 2023-12-12 20:50:47 | 只看该作者 Only view this author
您好,请教一下,如何用gaussian16计算两分子间不同态之间的电子耦合?非常感谢

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 18:23 , Processed in 0.187390 second(s), 23 queries , Gzip On.

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