24小时热门版块排行榜    

查看: 2512  |  回复: 20
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

清风凌云007

金虫 (小有名气)

[求助] 用G09优化激发态及荧光光谱计算 已有1人参与

我要算一个体系的荧光光谱,看了好多小木虫的帖子越看越糊涂了,还看了G09的英文说明,在SCRF这有说明(http://www.gaussian.com/g_tech/g_ur/k_scrf.htm),我只算气态的,总结了下是不是需要step1、2、4、5、6?step7是否需要计算?希望计算高手们、前辈们给指点迷津,我这样算到底对不对,如果不对应该怎么算呢,现在真是非常着急,都不想放假了,先在此谢过各位大神!!我的计算步骤如下(现在在进行第三步),
1、# opt freq b3lyp/6-31g(d)             基态分子结构优化及频率计算
2、# td(nstates=20) b3lyp/6-31g(d)           吸收光谱计算
3、# opt td(nstates=20,root=1,direct) b3lyp/6-31g(d)            单重态第一激发态优化(优化基础应该用优化基态结果还是吸收光谱计算结果?????)
4、# freq td(nstates=20,root=1,direct) b3lyp/6-31g(d)           单重态第一激发态优化结构频率计算
5、# td(nstates=20) b3lyp/6-31g(d)            单重态第一激发态发射
回复此楼

» 收录本帖的淘帖专辑推荐

经验 自己感觉好的东西

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清风凌云007

金虫 (小有名气)

送红花一朵
引用回帖:
4楼: Originally posted by Jasminer at 2015-01-24 05:32:20
你的步骤没错

多谢多谢
8楼2015-01-24 13:50:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答

清风凌云007

金虫 (小有名气)

怎么没人回复呢 呜呜  各位高手请进贴回复吧 非常感谢
2楼2015-01-23 19:33:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
5楼2015-01-24 10:36:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
清风凌云007: 金币+50, ★★★很有帮助, 希望大神继续指导我的回帖,多谢呀 2015-01-24 14:00:50
手册上的例子还有我们在论坛里讨论的都是在溶剂中的,需要考虑溶剂的平衡化和非平衡化的作用,
还有在TDDFT和DFT过程中,溶剂的作用。所以需要1-7,确切的说是1.3.4.5.6.7这几部。
如果说,你只计算在空气中的发射谱,那么这个很简单,下面的步骤
1:# opt freq b3lyp/6-31g(d)             基态分子结构优化及频率计算
2:# opt td(nstates=20,root=1,direct) b3lyp/6-31g(d)            单重态第一激发态优化(优化基础应该用优化基态结果还是吸收光谱计算结果?????
ps:针对你的问题,是这样的,这是激发态的优化,如果说计算发射,我们要看从s1-s0的发射,基于s0
优化的结构优化s1

在进行完上述步骤之后,在优化s1的输出文件中,会有如下的数据:
Excitation energies and oscillator strengths:

Excited State   1:      Singlet-A      1.7417 eV  711.84 nm  f=0.0001  <S**2>=0.000
     167 ->168        -0.70045
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -2006.02553586   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:      Singlet-A      2.5732 eV  481.82 nm  f=0.1422  <S**2>=0.000
     167 ->169        -0.69630

Excited State   3:      Singlet-A      2.5993 eV  476.99 nm  f=0.0028  <S**2>=0.000
     165 ->168         0.44525
     166 ->168         0.52944
默认的是3个,因为默认的nstates=3.看第一个td的波长就是发射的荧光的波长

ps:很多人,还在优化了s1的结构后,又做了一个td,其实没有必要,当然你以可以进行,计算出来的结果,和我们优化s1出来的结构的第一td是一样的。

综上所述,计算空气中的荧光,只需要2步,优化s0,基于s0优化s1,然后从s1输出文件找第一个td波长。
over

» 本帖已获得的红花(最新10朵)

It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
6楼2015-01-24 10:53:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 19:34 by lXgU0RvNTDC9
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 19:34 by lXgU0RvNTDC9
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 19:22 by lXgU0RvNTDC9
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 18:58 by lXgU0RvNTDC9
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 18:46 by lXgU0RvNTDC9
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +3 process2012 2026-08-23 4/200 2026-08-23 12:54 by huixian257
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见