计算化学公社

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

[新手求助] 求助:用高斯优化分子体系得出的能隙值大小顺序与实验不符

[复制链接 Copy URL]

27

帖子

0

威望

311

eV
积分
338

Level 3 能力者

老师您好,最近我模拟了如下图的两个分子,首先使用高斯对两分子进行优化(b3lyp/6-31g(d)),并得出HOMO,LUMO和能隙值,但两分子能隙值大小顺序和文献中相反并且差距较大(如下图),我检查了分子结构并无错误,优化完的二面角也与实验数据非常相似,并且泛函和实验中用的一样。所以我疑惑到底是哪里出了错误?分子体系及优化完的结构如下。(我给通讯作者发邮件想要一下输入文件但一直没有回)

紧接着我用优化完的文件计算了吸收光谱(CAM-b3lyp/6-311g(d,p)),计算结果可以和实验吻合。

202308111210148752..png (25.86 KB, 下载次数 Times of downloads: 13)

202308111210148752..png

202308111200587808..png (45.61 KB, 下载次数 Times of downloads: 14)

202308111200587808..png

21

帖子

0

威望

319

eV
积分
340

Level 3 能力者

2#
发表于 Post on 2023-8-11 13:14:50 | 只看该作者 Only view this author
https://doi.org/10.1016/j.dyepig.2023.111322
“In order to better analyze the photoelectric properties of the three dyes, the electrostatic potential (ESP) of CS-70–72 was calculated using the DFT (B3LYP/6-31G*) level of theory.”
"To further research into the optimized ground-state molecular geometries of CS-70–72, the density functional theory (DFT) calculations were performed with the Gaussian 09 program package at the B3LYP/6-311G(d,p) level."
会不会是看错了,原文优化基态用的基组是6-311G(d,p)

27

帖子

0

威望

311

eV
积分
338

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-8-11 14:13:09 | 只看该作者 Only view this author
w5326178 发表于 2023-8-11 13:14
https://doi.org/10.1016/j.dyepig.2023.111322
“In order to better analyze the photoelectric propert ...

你好,基组我确实没有用实验上的6-311G(d,p),是由于一开始我用这个基组优化了cs-71这个分子,能隙值为
2.11,和实验差的更多,另一方面是有点费时间。但我认为即使基组不一样也只会导致数值上有所不同,不会导致计算结论与实验结论不同吧?

1万

帖子

0

威望

7402

eV
积分
18170

Level 6 (一方通行)

4#
发表于 Post on 2023-8-11 20:20:21 | 只看该作者 Only view this author
学物理不秃头 发表于 2023-8-11 07:13
你好,基组我确实没有用实验上的6-311G(d,p),是由于一开始我用这个基组优化了cs-71这个分子,能隙值为
...

仔细阅读http://sobereva.com/543,了解为什么实验测得的所谓HOMO、LUMO能级不是真正的HOMO、LUMO能级,以及为什么和计算对不上很正常
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

5万

帖子

99

威望

5万

eV
积分
112492

管理员

公社社长

5#
发表于 Post on 2023-8-12 03:34:04 | 只看该作者 Only view this author
描述得很混乱
“泛函和实验中用的一样”,难道你说的是实验是理论计算?理论计算一般不能叫做实验
如果文献里就是做的理论计算,重复不出来这种问题下文该说的都说了
谈谈重复不出来计算化学文献里的数据的可能原因
http://sobereva.com/678http://bbs.keinsci.com/thread-38911-1-1.html

如果你说的实验就是一般意义的实验,HOMO、LUMO能级是根本没法通过实验测定的,4L提到的我的博文里明确说了
北京科音自然科学研究中心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!

27

帖子

0

威望

311

eV
积分
338

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2023-8-12 09:28:45 | 只看该作者 Only view this author
wzkchem5 发表于 2023-8-11 20:20
仔细阅读http://sobereva.com/543,了解为什么实验测得的所谓HOMO、LUMO能级不是真正的HOMO、LUMO能级, ...

好的,谢谢老师

27

帖子

0

威望

311

eV
积分
338

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2023-8-12 10:26:05 | 只看该作者 Only view this author
sobereva 发表于 2023-8-12 03:34
描述得很混乱
“泛函和实验中用的一样”,难道你说的是实验是理论计算?理论计算一般不能叫做实验
如果文 ...

不好意思老师,我表达的有些错误。文献中确实是一般的实验,但他除了对设计的分子进行实验测试之外,还进行了部分理论计算,我上面所说的泛函就是他理论计算优化时的泛函。但是文献中说HOMO能级是通过所谓的循环伏安法测得,我之前看了您的那篇博文明白了他的错误性。他的能隙是通过公式E0-0 = 1241/λonset计算的(我看了某些文章这个公式计算的应该是指的光学能隙吧?)。但是我疑惑的还是对于这两个分子为什么我这个理论上的HOMO-LUMO gap能隙值和实验上这个E0-0大小顺序相反?不应该是有正相关性?课题组的师姐说对应的不好就不行了,但她们也不知道哪错了。

1万

帖子

0

威望

7402

eV
积分
18170

Level 6 (一方通行)

8#
发表于 Post on 2023-8-12 15:48:54 | 只看该作者 Only view this author
学物理不秃头 发表于 2023-8-12 03:26
不好意思老师,我表达的有些错误。文献中确实是一般的实验,但他除了对设计的分子进行实验测试之外,还进 ...

确实有正相关性,但是这个正相关性是统计上的关系,不是必然联系。比如假设你研究了100个化合物,把E0-0对HOMO-LUMO gap作图,会发现有一定的线性,但是那是一大堆散点落在一根直线周围的那种线性,这些点并没有严格落在那条直线上。所以你任意取两个散点,就会有较小概率发现这两个点里面HOMO-LUMO gap的化合物E0-0反倒小
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

27

帖子

0

威望

311

eV
积分
338

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2023-8-12 18:47:53 | 只看该作者 Only view this author
wzkchem5 发表于 2023-8-12 15:48
确实有正相关性,但是这个正相关性是统计上的关系,不是必然联系。比如假设你研究了100个化合物,把E0-0 ...

感谢感谢老师,我明白了,祝您学术长青

本版积分规则 Credits rule

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

GMT+8, 2024-11-26 22:53 , Processed in 0.165983 second(s), 24 queries , Gzip On.

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