| 查看: 7226 | 回复: 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 ] |
» 收录本帖的淘帖专辑推荐
高斯及相关软件的使用 | 量子化学-第一性原理 | 气态分子化学动力学 | 研究生 |
计算化学 | 不懂计算 |
» 猜你喜欢
求助立方晶体In(OH)3的CIF卡片
已经有2人回复
大湾区大学刘天辉课题组招收2026级博士生(与中山大学物理学院联合培养)
已经有22人回复
物理化学论文润色/翻译怎么收费?
已经有173人回复
PbS量子点紫外吸收
已经有0人回复
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
[电子教材]天津大学物理化学(第七版)上、下册
已经有2人回复
阴离子交换膜电解二氧化碳还原
已经有0人回复
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
中国科学院成都山地灾害与环境研究所2026年特别研究助理(博士后)招聘启事
已经有0人回复
九江学院2026年最新高层次人才招聘公告
已经有0人回复
深圳信息职业技术学院2026年教师招聘公告(最新)
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐,对您同样有帮助:
td做激发态优化出错
已经有19人回复
gaussian 计算得到的荧光光谱选择其最大发射光谱是看f最大对应的吗?
已经有14人回复
gaussian 激发态优化后的chk文件和荧光计算后的chk文件有何不同
已经有5人回复
[求助]:高斯中的CIS是否能够优化开壳层的激发态结构,例如双重态,四重态的激发态
已经有8人回复
【求助】用TDDFT优化出来激发态之后的布居分析
已经有20人回复
【求助】单重激发态的高斯算法
已经有29人回复
【求助】Gaussian09 优化激发态出错
已经有8人回复
【求助】如何计算有机小分子的三重态能级(T1-S0)?
已经有12人回复
» 抢金币啦!回帖就可以得到:
高级研发工程师/研发工程师 – 锂离子电池正极材料回收与制备-香港纳米及先进材料研发
+1/283
山东理工大学材料学院泰山学者青年专家招收2026年调剂硕士研究生
+1/76
哈工大 化工学院招26级博士生(钙钛矿太阳能电池/传感器方向)
+1/33
中山大学公卫学院招收2026年4月申请-考核制的博士生补录取(化学生物学、蛋白质组学)
+3/25
26届计算机、电子信息类、电科、控制、通信考研T j信息pp骏:74+08+76+6+20/
+1/14
【青岛大学】2026年生物与医药申请考核制博士生招生(含少数民族骨干人才)
+1/12
26届计算机、电子信息类、电科、控制、通信考研T j信息pp骏:74+08+76+6+20/
+1/9
中国科学院宁波材料所 极端环境高分子工程材料团队 诚聘博士后
+1/8
温州大学李新华教授课题组2026年招收硕士研究生信息
+1/7
look 26届计算机、电子信息类、电科、控制、通信考研T j信息醇:10+33+77+47+03
+1/7
湖北工业大学2026年博士研究生(生物医药专博+轻工技术与工程学博)
+1/5
26计算机、电子信息类、电科、控制、通信考研T j醇:10+33+77+47+03
+1/5
诚邀有志之士加入江苏大学环境学院合成生物学课题组!(长期接收学硕专硕、硕博连读)
+1/4
11408 372分 求调剂 有顶会论文国创竞赛六级
+1/4
140300设计,一志愿985,354求调剂
+1/4
【全奖, 免GRE】美国科罗拉多大学科罗拉多斯普林斯分校----固体力学
+1/3
环境纳米技术@水质风险控制 课题组接收研究生报考/调剂公告
+1/3
燕山大学-意大利米兰理工大学国际化科研团队招收硕士研究生
+1/3
新加坡南洋理工大学- 光电/ 智能传感/ 脑机接口方向 博士后
+1/1
青岛科技大学高分子学院--让我们一起做有趣的研究
+1/1
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













回复此楼



