计算化学公社

标题: TDDFT计算六甲基苯碳正离子荧光光谱该用什么方法和基组 [打印本页]

作者
Author: bingomini 时间: 2020-2-22 08:37
标题: TDDFT计算六甲基苯碳正离子荧光光谱该用什么方法和基组
各位老师好，
我想计算六甲基苯碳正离子的荧光光谱。我阅读了老师的采用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

请老师帮助，谢谢！

作者
Author: clwater 时间: 2020-2-22 11:10
计算流程没有什么问题。

以后不要说“报错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够了。

作者
Author: bingomini 时间: 2020-2-22 11:35

clwater 发表于 2020-2-22 11:10
计算流程没有什么问题。

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

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

作者
Author: sobereva 时间: 2020-2-22 17:44
dgdzvp比起大小相仿佛的6-31G*、def2-SV(P)根本没什么额外好处，还非常非主流，用这种东西很容易招致审稿人问你为什么用这个、怎么论证合理性，到时候很麻烦

作者
Author: bingomini 时间: 2020-2-27 17:48

sobereva 发表于 2020-2-22 17:44
dgdzvp比起大小相仿佛的6-31G*、def2-SV(P)根本没什么额外好处，还非常非主流，用这种东西很容易招致审稿人 ...

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

作者
Author: bingomini 时间: 2020-2-27 19:01

sobereva 发表于 2020-2-22 17:44
dgdzvp比起大小相仿佛的6-31G*、def2-SV(P)根本没什么额外好处，还非常非主流，用这种东西很容易招致审稿人 ...

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

作者
Author: sobereva 时间: 2020-2-28 02:58

bingomini 发表于 2020-2-27 17:48
谢谢老师解答，我采用了6-31G(d,p)计算出来了。还想请教老师，DFT或者Multiwfn可以计算荧光量子产率吗？

DFT和Multiwfn不是并列的关系

目前没有靠谱的方法能定量算荧光量子产率，因为牵扯的因素太多，去公社论坛首页google框搜索，之前有很多讨论

作者
Author: sobereva 时间: 2020-2-28 02:58

bingomini 发表于 2020-2-27 19:01
老师，我想请教：用Multiwfn做出来的荧光光谱谱图中，发射波长对应的纵坐标“Molecular adsorption coeff ...

不能

作者
Author: bingomini 时间: 2020-2-28 08:34

sobereva 发表于 2020-2-28 02:58
DFT和Multiwfn不是并列的关系

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

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

作者
Author: sobereva 时间: 2020-3-1 23:20

bingomini 发表于 2020-2-28 08:34
谢谢老师指正。我采用了TDDFT方法计算出多甲基苯碳正离子的荧光性质，然后采用Multiwfn处理获得了该物质 ...

发光强度

作者
Author: bingomini 时间: 2020-3-2 15:09

sobereva 发表于 2020-3-1 23:20
发光强度

谢谢老师！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)