计算化学公社

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

[综合交流] 同体系同泛函不同基组计算出的HOMO-LUMO gap差别很大

[复制链接 Copy URL]

171

帖子

0

威望

1853

eV
积分
2024

Level 5 (御坂)

15个丙烯单体构成的PP寡聚分子,使用gaussian16进行优化和振动分析,采用B3LYP泛函和D3(BJ)色散校正。
做了两次计算,仅基组不同。第一次用6-31G*,第二次用6-311G**.
计算获得的dos谱中,HOMO变化不大,但LUMO能量差距非常明显,导致gap差别也非常明显。
请教一下可能的原因?我做了多次测试,猜测可能不是构象的原因。哪种结果比较可信呢?
下图为6-31G*


下图为6-311G**

171

帖子

0

威望

1853

eV
积分
2024

Level 5 (御坂)

11#
 楼主 Author| 发表于 Post on 2022-11-8 18:20:31 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-8 17:28
可以,不过加大一点基组有益无害

好的 谢谢老师的详细指教!

1万

帖子

0

威望

8983

eV
积分
20731

Level 6 (一方通行)

10#
发表于 Post on 2022-11-8 17:28:08 | 只看该作者 Only view this author
rugals 发表于 2022-11-8 10:17
缺陷是靠对聚丙烯分子链的一些修饰进行引入。比如链节中间共聚一些特殊结构,或者搞上去一些侧链。
预想 ...

可以,不过加大一点基组有益无害
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

171

帖子

0

威望

1853

eV
积分
2024

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2022-11-8 17:17:07 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-8 16:47
那样的话主要取决于缺陷能不能算准。聚丙烯的带隙这么宽,带隙差个1~2eV也不影响本征聚丙烯完全测不出导 ...

缺陷是靠对聚丙烯分子链的一些修饰进行引入。比如链节中间共聚一些特殊结构,或者搞上去一些侧链。
预想的是它能够在带隙中引入一些能级。这样的话是否可以用刚才我提到的模拟和计算策略进行分析?

1万

帖子

0

威望

8983

eV
积分
20731

Level 6 (一方通行)

8#
发表于 Post on 2022-11-8 16:47:36 | 只看该作者 Only view this author
rugals 发表于 2022-11-8 08:44
谢谢老师!
我是想评价一下聚合物材料的导电性。将它类比于半导体的能带结构分析其绝缘和电荷陷阱相关的 ...

那样的话主要取决于缺陷能不能算准。聚丙烯的带隙这么宽,带隙差个1~2eV也不影响本征聚丙烯完全测不出导电性的事实,但是缺陷是可以导致导电性明显变化的。所以应该以你要研究的缺陷为准来选择计算方法
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

171

帖子

0

威望

1853

eV
积分
2024

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2022-11-8 15:44:21 | 只看该作者 Only view this author
本帖最后由 rugals 于 2022-11-8 15:50 编辑
wzkchem5 发表于 2022-11-8 15:14
如果忽略B3LYP的误差的话,你这个计算结果表明聚丙烯被光激发以后,电子很可能会直接被激发到真空里。此 ...

谢谢老师!
我是想评价一下聚合物材料的导电性。将它类比于半导体的能带结构分析其绝缘和电荷陷阱相关的一些性质。

我准备用MD模拟出物理性质稳定的团簇(或者在模拟大体系MD后抠一个团簇)后再去做dos。这样的话是否可以还用Gaussian+B3LYP/6-31G*的组合?或者必须用第一性原理?

171

帖子

0

威望

1853

eV
积分
2024

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2022-11-8 15:35:33 | 只看该作者 Only view this author
sobereva 发表于 2022-11-8 12:31
本来MO能量就是虚构的概念,尤其是空轨道的物理意义更差,绝非基组越大结果越好。甚至基组过大时由于泛函的 ...

谢谢老师!

1万

帖子

0

威望

8983

eV
积分
20731

Level 6 (一方通行)

5#
发表于 Post on 2022-11-8 15:14:08 | 只看该作者 Only view this author
rugals 发表于 2022-11-7 12:18
谢谢老师,这篇博文我看过了。我认为是基组的原因。但我没法确定哪个结果更合理一些。

因为按照博文的 ...

如果忽略B3LYP的误差的话,你这个计算结果表明聚丙烯被光激发以后,电子很可能会直接被激发到真空里。此时取决于你计算HOMO-LUMO gap的目的:(1)如果你就是要算气相里的孤立分子,那应该取大基组的结果,甚至应该加一些弥散函数再算,但是这个结果是不能和体相的带隙结果相比的;(2)如果你要算体相聚合物带隙的话,建议用周期性的程序,或者用embedded cluster model可能也行,主要是不让电子激发到真空里去。
但是B3LYP的误差不能简单忽略,因为B3LYP虽然算带隙比较合适,但是算电离能,尤其是大分子的电离能,有系统误差(因为自相互作用误差的缘故),一旦电离能不准,真空能级和LUMO的相对位置就会不准。
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

5万

eV
积分
120101

管理员

公社社长

4#
发表于 Post on 2022-11-8 12:31:23 | 只看该作者 Only view this author
本来MO能量就是虚构的概念,尤其是空轨道的物理意义更差,绝非基组越大结果越好。甚至基组过大时由于泛函的离域性误差问题反倒算出更没化学意义的LUMO
怎么计算得当看实际具体目的
北京科音自然科学研究中心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

171

帖子

0

威望

1853

eV
积分
2024

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2022-11-7 19:18:33 | 只看该作者 Only view this author
本帖最后由 rugals 于 2022-11-7 19:21 编辑

谢谢老师,这篇博文我看过了。我认为是基组的原因。但我没法确定哪个结果更合理一些。

因为按照博文的意思,B3LYP/6-31G*很适合算聚合物。但是6-311G**比它级别更高一些。

1万

帖子

0

威望

8983

eV
积分
20731

Level 6 (一方通行)

2#
发表于 Post on 2022-11-7 18:35:16 | 只看该作者 Only view this author
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

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

GMT+8, 2025-8-13 22:16 , Processed in 0.163720 second(s), 24 queries , Gzip On.

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