计算化学公社

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

[Gaussian/gview] TDDFT计算六甲基苯碳正离子荧光光谱该用什么方法和基组

[复制链接 Copy URL]

24

帖子

0

威望

81

eV
积分
105

Level 2 能力者

各位老师好,
    我想计算六甲基苯碳正离子的荧光光谱。我阅读了老师的采用TDDFT计算光谱性质的文章:http://sobereva.com/314。但在实际操作中遇到了问题。
    我的计算方法及流程如下:1. 将六甲基碳正离子进行优化,优化方法为 rb3lyp/dgdzvp,检查后得到没有虚频的优化结构;2.在所获得的优化结构基础上加上关键词TD进行TDDFT计算。
    遇到的问题:1.报错L502.exe;2.我还尝试了 wb97xd方法以及6-31g基组的相互组合,但报了同样的错误。
    求助目的:我想请教各位老师,我是不是方法使用错了,或者TDDFT计算上设置遗漏或者有错误?碳正离子带正电荷charge给了+1或者给为0我都尝试了但都无法正常结束计算。如下是我的结构文档:

%nprocshared=12
%mem=12GB
# opt rb3lyp/dgdzvp geom=connectivity td

B6p

1 1
C                  1.27065100   -0.71939500   -0.24291800
C                  0.00002100   -1.48994600   -0.18090600
C                 -1.27062900   -0.71942400   -0.24290100
C                 -1.26049600    0.65275700   -0.07787900
C                 -0.00000500    1.31863600    0.07728100
C                  1.26049000    0.65278700   -0.07787800
C                  2.54337800    1.45968400   -0.05974700
C                 -0.00003800    2.79337600    0.35508000
C                 -2.54340800    1.45960200   -0.05975800
C                 -2.53589600   -1.50686600   -0.44135100
C                  0.00003200   -2.32406200    1.16746400
C                  2.53591800   -1.50683800   -0.44136700
H                  0.00002400   -2.23917800   -0.98599400
H                  2.47085100    2.34509600   -0.69694800
H                  2.78268000    1.80382800    0.95303600
H                  3.39382700    0.87664000   -0.41048700
H                  0.88686000    3.11124500    0.90329400
H                 -0.00029300    3.34338600   -0.59808400
H                 -0.88676400    3.11110500    0.90366400
H                 -2.47081400    2.34521800   -0.69666200
H                 -2.78292200    1.80342400    0.95308800
H                 -3.39376300    0.87662700   -0.41084700
H                 -2.34033800   -2.57348500   -0.56499000
H                 -3.05947600   -1.16181800   -1.34071900
H                 -3.23045200   -1.37784200    0.39622600
H                 -0.88418500   -2.96137500    1.21329200
H                 -0.00011200   -1.65052800    2.02743000
H                  0.88440200   -2.96115500    1.21342800
H                  3.05962300   -1.16169000   -1.34061800
H                  3.23038200   -1.37793800    0.39631400
H                  2.34035100   -2.57343900   -0.56514300

1 2 1.0 6 2.0 12 1.0
2 3 1.0 11 1.0 13 1.0
3 4 2.0 10 1.0
4 5 1.5 9 1.0
5 6 1.5 8 1.0
6 7 1.0
7 14 1.0 15 1.0 16 1.0
8 17 1.0 18 1.0 19 1.0
9 20 1.0 21 1.0 22 1.0
10 23 1.0 24 1.0 25 1.0
11 26 1.0 27 1.0 28 1.0
12 29 1.0 30 1.0 31 1.0
13
14
15
16
17


    请老师帮助,谢谢!

8

帖子

0

威望

1182

eV
积分
1190

Level 4 (黑子)

2#
发表于 Post on 2020-2-22 11:10:00 | 只看该作者 Only view this author
计算流程没有什么问题。

以后不要说“报错L502.exe”,参见http://bbs.keinsci.com/thread-4829-1-1.html,此文给出常见的 Gaussian FAQ。
“提问时请不要:仅说/仅截图 lXXX.exe 报错;仅按照 lXXX.exe 报错解决问题。问“L1234.exe报错该怎么解决”这种话没人能回答。Gaussian报错时会给出出错link名称,但哪怕是最容易确定的L502是SCF不收敛都有例外情况,其他的link报错几乎都需要此行之前的具体报错信息才能提供解决方案。同时也切勿仅按照Link数寻找解决方案,例如本文后面说“L801,Excessive mixing of core and valence orbitals 错误”应当使用Full选项解决,但从未说“L801错误”都应该用Full解决,一定要看具体报错信息。”

个人没有使用过dgdzvp基组,无法评论。对这种简单有机正离子体系,常建议使用6-31G*或者def2-SVP基组,如果有进一步计算资源,可使用6-311G**、 def-TZVP等。以及“6-31G、6-311G不加极化函数根本没法用!”,可参见http://sobereva.com/336

方法我觉得B3LYP够了。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

24

帖子

0

威望

81

eV
积分
105

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2020-2-22 11:35:51 | 只看该作者 Only view this author
clwater 发表于 2020-2-22 11:10
计算流程没有什么问题。

以后不要说“报错L502.exe”,参见http://bbs.keinsci.com/thread-4829-1-1.htm ...

感谢您的指正与回答。一定注意。
我试试看更换基组。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2020-2-22 17:44:05 | 只看该作者 Only view this author
dgdzvp比起大小相仿佛的6-31G*、def2-SV(P)根本没什么额外好处,还非常非主流,用这种东西很容易招致审稿人问你为什么用这个、怎么论证合理性,到时候很麻烦
北京科音自然科学研究中心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

24

帖子

0

威望

81

eV
积分
105

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2020-2-27 17:48:07 | 只看该作者 Only view this author
sobereva 发表于 2020-2-22 17:44
dgdzvp比起大小相仿佛的6-31G*、def2-SV(P)根本没什么额外好处,还非常非主流,用这种东西很容易招致审稿人 ...

谢谢老师解答,我采用了6-31G(d,p)计算出来了。还想请教老师,DFT或者Multiwfn可以计算荧光量子产率吗?

24

帖子

0

威望

81

eV
积分
105

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2020-2-27 19:01:20 | 只看该作者 Only view this author
sobereva 发表于 2020-2-22 17:44
dgdzvp比起大小相仿佛的6-31G*、def2-SV(P)根本没什么额外好处,还非常非主流,用这种东西很容易招致审稿人 ...

老师,我想请教:用Multiwfn做出来的荧光光谱谱图中,发射波长对应的纵坐标“Molecular adsorption coefficient”峰值能够用来衡量荧光产率吗?不同物质的荧光产率在理论计算的哪个物理量可以比较或者计算呢?谢谢老师!

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

7#
发表于 Post on 2020-2-28 02:58:17 | 只看该作者 Only view this author
bingomini 发表于 2020-2-27 17:48
谢谢老师解答,我采用了6-31G(d,p)计算出来了。还想请教老师,DFT或者Multiwfn可以计算荧光量子产率吗?

DFT和Multiwfn不是并列的关系

目前没有靠谱的方法能定量算荧光量子产率,因为牵扯的因素太多,去公社论坛首页google框搜索,之前有很多讨论
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

8#
发表于 Post on 2020-2-28 02:58:32 | 只看该作者 Only view this author
bingomini 发表于 2020-2-27 19:01
老师,我想请教:用Multiwfn做出来的荧光光谱谱图中,发射波长对应的纵坐标“Molecular adsorption coeff ...

不能
北京科音自然科学研究中心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

24

帖子

0

威望

81

eV
积分
105

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2020-2-28 08:34:01 | 只看该作者 Only view this author
sobereva 发表于 2020-2-28 02:58
DFT和Multiwfn不是并列的关系

目前没有靠谱的方法能定量算荧光量子产率,因为牵扯的因素太多,去公社 ...

谢谢老师指正。我采用了TDDFT方法计算出多甲基苯碳正离子的荧光性质,然后采用Multiwfn处理获得了该物质的荧光光谱。我想请问所获得的荧光光谱的纵坐标物理意义是什么呢?比如紫外可见光谱的纵坐标是吸光系数,能反应物质对这个波长光的吸收性质。我想知道荧光的纵坐标的物理意义是什么,还请老师指教,谢谢!C:\Users\bingo\Desktop

DISLIN_1.PNG (33.45 KB, 下载次数 Times of downloads: 28)

DISLIN_1.PNG

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

10#
发表于 Post on 2020-3-1 23:20:37 | 只看该作者 Only view this author
bingomini 发表于 2020-2-28 08:34
谢谢老师指正。我采用了TDDFT方法计算出多甲基苯碳正离子的荧光性质,然后采用Multiwfn处理获得了该物质 ...

发光强度
北京科音自然科学研究中心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

24

帖子

0

威望

81

eV
积分
105

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2020-3-2 15:09:30 | 只看该作者 Only view this author

谢谢老师!

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 14:05 , Processed in 0.289745 second(s), 24 queries , Gzip On.

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