24小时热门版块排行榜    

查看: 2059  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

小范范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的回帖

小范范1989

木虫 (著名写手)

引用回帖:
7楼: Originally posted by 清风凌云007 at 2015-01-24 13:49:13
多谢大神指教,我看了两个帖子,还是有一点不明白,就是怎么做发射光谱图呢,要是作图是要在优化激发态基础上做TD 吗?也就是我的最后一步
5、# td(nstates=20) b3lyp/6-31g(d)            单重态第一激发态发射
...

按照我们这个计算方式,计算的只是一个波长,
如果说非要做出谱线,那么就采用下面的步骤:
1:那你计算td的文件的输出文件.log或者是.out,放到高斯view中,
点击results,然后选择UV-VIS,会看到一个图像,如下图一,但是这个图像时我们计算td时候全部的,我们要看的是第一个td对应的发射。所以这个图像不是。
2:在图像的中,在x轴上,找到我们输出文件第一个td对应的波长,假如说是350nm。就在350nm处,右键。--save data-(可以设置开始,结束,还有间隔范围)---保存数据后,使用origin作图。如图二
3:图3才是我们的荧光光谱。
用G09优化激发态及荧光光谱计算
3.png


用G09优化激发态及荧光光谱计算-1
1.png


用G09优化激发态及荧光光谱计算-2
2.jpg

It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
12楼2015-01-24 16:00:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

清风凌云007

金虫 (小有名气)

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

清风凌云007

金虫 (小有名气)

哎呀呀呀
3楼2015-01-23 19:56:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

你的步骤没错

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

4楼2015-01-24 05:32:20
已阅   回复此楼   关注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的回帖

清风凌云007

金虫 (小有名气)

送红花一朵
引用回帖:
6楼: Originally posted by 小范范1989 at 2015-01-24 10:53:34
手册上的例子还有我们在论坛里讨论的都是在溶剂中的,需要考虑溶剂的平衡化和非平衡化的作用,
还有在TDDFT和DFT过程中,溶剂的作用。所以需要1-7,确切的说是1.3.4.5.6.7这几部。
如果说,你只计算在空气中的发射 ...

多谢大神指教,我看了两个帖子,还是有一点不明白,就是怎么做发射光谱图呢,要是作图是要在优化激发态基础上做TD 吗?也就是我的最后一步
5、# td(nstates=20) b3lyp/6-31g(d)            单重态第一激发态发射

还有就是画图要用chk文件吗?只有out或log文件应该怎么画图呢?在我的输入文件(如下)中输入%chk=2.gif.chk就是生成了chk文件吧,但是我打不开,看了好多论坛上的帖子有的说要转化成fchk,但是我按照说的转化不了
%chk=2.chk
%mem=100MW
%nproc=1
# opt b3lyp/6-31+g
7楼2015-01-24 13:49:13
已阅   回复此楼   关注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的回帖

Jasminer

铁杆木虫 (著名写手)

引用回帖:
7楼: Originally posted by 清风凌云007 at 2015-01-24 13:49:13
多谢大神指教,我看了两个帖子,还是有一点不明白,就是怎么做发射光谱图呢,要是作图是要在优化激发态基础上做TD 吗?也就是我的最后一步
5、# td(nstates=20) b3lyp/6-31g(d)            单重态第一激发态发射
...

作图不需chk,任一能作展宽的软件都可以,比如说GV。

[ 发自小木虫客户端 ]
9楼2015-01-24 14:21:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清风凌云007

金虫 (小有名气)

引用回帖:
9楼: Originally posted by Jasminer at 2015-01-24 14:21:47
作图不需chk,任一能作展宽的软件都可以,比如说GV。
...

就是说在GaussView中打开out或log文件,在results中选哪个选项呢?
图片是我算了TD吸收后点UV-VIS后出现的图,别的就不会了
用G09优化激发态及荧光光谱计算-3
360截图20150124144958116.jpg

10楼2015-01-24 14:52:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 清风凌云007 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 265求调剂 +10 梁梁校校 2026-03-17 10/500 2026-03-23 21:17 by 一切OK
[考研] 一志愿华东理工大学081700,初试分数271 +3 kotoko_ik 2026-03-23 4/200 2026-03-23 20:47 by 我的船我的海
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 306求调剂 +9 chuanzhu川烛 2026-03-18 9/450 2026-03-23 13:17 by luoyongfeng
[考研] 北科281学硕材料求调剂 +8 tcxiaoxx 2026-03-20 9/450 2026-03-23 12:16 by tcxiaoxx
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 311求调剂 +6 冬十三 2026-03-18 6/300 2026-03-22 20:18 by edmund7
[考研] 315分,诚求调剂,材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 384求调剂 +3 子系博 2026-03-22 4/200 2026-03-22 11:04 by 搏击518
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
信息提示
请填处理意见