| 查看: 8272 | 回复: 36 | |||||||||
[交流]
【求助】关于激发态优化及发射光谱的计算
|
|||||||||
|
有3个问题想问下各位虫友。 1、用gaussian计算激发态的时候nstates和root都是默认的，后来我优化激发态构型，然后再次做TD计算发射光谱也都是默认的态，实际上我想要观察的是Excited State 3的变化，请问有没有指定Excited State3得出的结果会不一样吗？比如优化激发态构型得到的稳定构型会不一样吗？ 2、我主要是想计算发射光谱。见论坛有人提到，在对激发态构型进行优化时，第一个循环给出的激发能就是吸收谱，最后一步给出的激发能就是荧光谱。我想问下这里说的吸收谱是指紫外吸收光谱吗？荧光谱是指荧光吸收谱吗？ 3、我开始计算物质A，能正常结束。 计算A与B的复合物的时候（物质A是笼型的，然后物质B是个小分子，进到笼内被吸附），结构优化正常，无虚频，但是计算激发态的时候报错了，重新调整了结构，优化了再计算激发态，还是收敛失败。 输入文件是： # opt td b3lyp/6-31g(d,p) geom=connectivity 报错信息是： >>>>>>>>>> Convergence criterion not met. SCF Done: E(RB3LYP) = -4648.00956095 A.U. after 129 cycles Convg = 0.1628D-05 -V/T = 2.0075 Convergence failure -- run terminated. Error termination via Lnk1e in /home/soft/g09/l502.exe at Sun Mar 27 13:08:44 2011. 是不是不能计算2个分子也就是复合物的激发态啊？ 希望各位高手能帮个忙~谢谢~ [ Last edited by ben_ladeng on 2011-3-31 at 06:00 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 高斯及相关软件的使用 | 量子化学-第一性原理 | 气态分子化学动力学 | 研究生 |
| 计算化学 | 不懂计算 |
» 猜你喜欢
华中农业大学植物科学技术学院李国田老师课题组劳动聘用制科研助理招聘公告
已经有20人回复
-
微纳hub，微纳加工的宝库
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有172人回复
-
推荐一款思路新颖的论文写作工具Flowing
已经有0人回复
-
醋酸盐体系溶胶凝胶法合成富锂锰基正极
已经有1人回复
-
牛磺酸钠
已经有1人回复
-
Renlab@monash全奖博士机会，最后三天！
已经有0人回复
-
微纳加工里的 lift-off 工艺
已经有0人回复
-
今年面上没有消息，想来是没了，散个金攒攒人品，为明年祈福
已经有52人回复
-
合肥工业大学854仪器技术综合专业课考研
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- td做激发态优化出错
已经有19人回复
- gaussian 计算得到的荧光光谱选择其最大发射光谱是看f最大对应的吗？
已经有14人回复
- gaussian 激发态优化后的chk文件和荧光计算后的chk文件有何不同
已经有5人回复
- [求助]：高斯中的CIS是否能够优化开壳层的激发态结构，例如双重态，四重态的激发态
已经有8人回复
- 【求助】用TDDFT优化出来激发态之后的布居分析
已经有20人回复
- 【求助】单重激发态的高斯算法
已经有29人回复
- 【求助】Gaussian09 优化激发态出错
已经有8人回复
- 【求助】如何计算有机小分子的三重态能级（T1-S0）？
已经有12人回复
» 抢金币啦！回帖就可以得到:
-
求助2013年-2011年，武汉哪所学校的实验室合成过单价3万/克（或毫升）的剧毒成品？
+5/2385
-
助理研究员——吉林大学-苏州国家实验室光电半导体智能设计与应用团队
+1/488
-
2027年博士生（硕博连读/直博）吉林大学-苏州国家实验室光电半导体智能设计与应用团队
+1/485
-
诚征女友
+1/261
-
华东理工大学 生物工程学院 生物反应器工程国家重点室 博士后招聘
+1/79
-
华东理工大学 生物反应器工程全国重点实验室 博士后招聘
+1/79
-
山东第一医科大学医工交叉团队2025年招收纳米医学、药物化学等领域副教授、博士后
+1/74
-
Postdoctoral fellow in Medicinal Chemistry, McGill University
+1/73
-
西湖大学物理学、光学、电子信息方向博士生招生（长期有效）
+5/35
-
坐标南京
+1/30
-
北京师范大学珠海校区 李朝强团队招多名硕士博士
+1/12
-
西安交通大学补亚忠课题组招收申请考核制博士生
+2/8
-
四川大学周加境课题组招博士研究生/研究助理（自组装材料/生物医学方向）
+1/6
-
多本Scopus收录国际英文期刊征稿，全免 APC！含医学期刊
+1/6
-
国家自然科学基金申报者发榜前心理状态调查，麻烦帮忙填一下哇！
+1/6
-
大家有research gate的账号吗？有偿借用
+1/4
-
宁波诺丁汉大学新能源与环境材料技术研究中心招收2027年春季入学博士生
+1/4
-
分不清 EI 会议 CA 和 JA？90% 投稿人踩过的检索大坑
+1/4
-
浙江大学王宇课题组诚聘博士后/博士研究生/科研助理（招聘长期有效）
+1/3
-
中科院大连化物所诚聘博士后/科研助理/项目聘用人员
+1/2
2楼2011-03-31 08:45:37
5楼2011-03-31 10:36:11
ar02hi3v(金币+1): 谢谢啊，是我理解能力差，呵呵~ 2011-03-31 20:29:23
|
好吧，你说的这个不是我的思路，我现在对自己的表达能力表示强烈怀疑。我的思路是下面这样的： 1、优化基态构型，#P opt freq b3lyp/6-31g(d,p) geom=connectivity 2、TD单点，得到吸收谱(#P td(NStates=6) b3lyp/6-31g(d,p) geom=connectivity) 之所以要算这一步是为了要确定计算出来的哪个激发态是你需要优化的 3、以基态结构为初始结构用TD优化激发态结构，这个过程中将会计算很多步垂直激发能，其中第一步所给出的垂直激发态能将和上一步的TD单点给出的完全相同，即吸收谱，而最后一步给出的则是荧光谱。（#P Opt td(NStates=6,root=3) b3lyp/6-31g(d,p) geom=connectivity) [ Last edited by pwzhou on 2011-3-31 at 15:32 ] |
14楼2011-03-31 15:28:37
|
关于第2个问题，这是原帖地址http://muchong.com/bbs/viewthread.php?tid=1872007 帖子的第3楼有说到，我没看明白：“第二步，激发态结构优化（以优化好的基态结构为初始构型），在激发态结构优化的第一个循环给出的激发能就是吸收谱，最后一步给出的激发能就是荧光谱。”这里指的是在激发态构型优化的输出文件里找的这两个吗？那么再次做TD后出来的激发态波长是什么谱？谢谢~ |
3楼2011-03-31 09:28:36
★ ★ ★
ar02hi3v(金币+5): 谢谢指教~ 2011-03-31 10:37:37
375642546(金币+3): 2011-04-08 16:55:56
ar02hi3v(金币+5): 谢谢指教~ 2011-03-31 10:37:37
375642546(金币+3): 2011-04-08 16:55:56
|
这么讲可能会更清楚一点，在优化好的基态构型上做TD得到的是吸收谱，在优化好的激发态构型上做TD得到的是荧光谱。其实好好看分子光谱的书，这个还是很容易理解的。原帖中说的激发态结构优化（以优化好的基态结构为初始构型），这样你第一步实际上就是在基态构型上做的TD，所以是吸收谱；当结构优化完成后，最后的那个结构就是激发态的结构，在此结构上做的TD就是荧光谱。 如果你要优化的是第三激发态的结构，那就必须指定root=3，否则你得到的是第一激发态的结构。 你那个不收敛是SCF不收敛，增大循环次数就是了scf(maxcycle=500) |
4楼2011-03-31 10:00:35
6楼2011-03-31 11:23:44
7楼2011-03-31 11:30:17
ar02hi3v(金币+1): 感谢积极参与讨论~ 2011-03-31 14:00:15
|
这个纯粹就是多此一举，画蛇添足的做法。 用TD优化完的结构当然就是激发态结构了，在这个结构上做的TD当然是荧光谱。在此基础上为什么还要对激发态进行CIS优化？要知道CIS和TD是优化激发态的两种不同方法，利用CIS和TD进行优化都可以得到激发态的结构，而且TD本身的计算结果要比CIS要好，为啥还要多此一举再优化一个CIS方法下的激发态结构？还非得在CIS优化好的激发态结构上再做TD？难道这个比在TD优化好的激发态结构上做TD得到的结果还要好？ 以前在G03的时候，TD没有解析梯度，无法有效的完成激发态的优化，所以通常的做法就是DFT优化基态结构，在基态结构上做TD得到吸收谱；然后用CIS优化激发态结构，最后在CIS优化的激发态结构上做TD得到荧光谱。但是在TD可以优化激发态结构的基础上，根本没必要再用CIS。 |
8楼2011-03-31 11:46:36
9楼2011-03-31 12:03:56
10楼2011-03-31 12:07:21
11楼2011-03-31 14:08:04
12楼2011-03-31 14:12:00
ar02hi3v(金币+5): 今天愚人节，全成版主啥的了，不过我的金币是真的哦，哈哈~ 2011-04-01 00:18:34
|
TD激发态优化过程中每一步都要对所得到的结构进行垂直激发能的计算的，优化的最后一步不就是在激发态结构上做垂直激发能计算吗？为何还要再单独算一步（当然你可以再算一步，但是这一步的结果和优化最后一步的结果是完全一致的，你可以测试一下） 每个激发态对应的轨道跃迁、跃迁偶极距以及谐振强度都是不一样的，假如你计算得到的第一个或者前几个激发态的谐振强度都是0，也就说这些态是暗态，而对应的吸收谱的跃迁则是第二或者更高的激发态，那么这时候你就要根据计算结果和你所关注的问题来选择你要优化的是那个激发态了。 总之就是具体问题具体分析，有些基础概念一定要搞清楚，最重要的不是明白怎么算，而是要明白为什么要这么算？ [ Last edited by pwzhou on 2011-3-31 at 19:02 ] |
16楼2011-03-31 19:01:32
|
我整理了一下你的思路，你看对不对。 1、优化基态构型，我用的# opt freq b3lyp/6-31g(d,p) geom=connectivity 2、做TD，得到吸收谱(我写的# opt td b3lyp/6-31g(d,p) geom=connectivity) 这里电脑应该是先进行opt，再进行TD的吧？现在想想，好像我这步的opt有点多余了。 3、优化激发态构型，# opt cis/6-31g(d,p) geom=connectivity cis=direct 这是以前算的，我用的是g09的，不知道你说的用TD-DTF的要怎么输入，麻烦给个具体的，还有就是细节方面有什么注意的吗？比如这一步的基组和下一步的基组要一致之类的。（要写上指定哪个态，这个我记住了）。 4、再次做TD得到发射谱，# b3lyp/6-31g(d,p) geom=connectivity td 问的有点啰嗦，希望大虾耐心赐教~ |
13楼2011-03-31 14:38:06
15楼2011-03-31 16:57:55
17楼2011-04-03 00:19:21
18楼2011-04-03 10:08:29
19楼2011-04-03 21:03:45
20楼2011-04-08 16:00:25
21楼2011-04-09 14:07:09
22楼2011-04-23 19:29:17
23楼2011-04-24 00:17:01
24楼2011-05-27 15:27:23
25楼2011-05-27 15:31:04
26楼2012-05-03 11:06:37
27楼2012-05-12 09:30:29
28楼2012-05-12 09:37:30
29楼2012-07-05 16:53:55
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
我想问一个问题，我算激发态优化时用的如下方法#p hf/6-31+g(d) cis(singlets,nstates=5) scrf(pcm,solvent=CH3CH2OCH2CH3 ) opt optcyc=300 maxdisk=20GB 出现了下面的结果 Root 1 : ******************** Change is ******************** Root 2 : ******************** Change is -177.750956439561200 Root 3 : ******************** Change is -109.476836921007000 Root 4 : ******************** Change is -6.854283804947429 Root 5 : ******************** Change is -58.866471450604870 这是怎么回事呢？请求你的帮助，我都算了一个多月了，还没出来，还不报错。 |
30楼2012-08-16 08:50:52
31楼2012-12-14 10:03:27
32楼2013-03-04 09:17:21
33楼2013-03-04 09:23:57
34楼2014-04-21 11:57:07
35楼2014-10-05 09:27:59
36楼2014-10-05 09:29:45
37楼2015-11-23 20:02:32