24小时热门版块排行榜    

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

xtdut

银虫 (小有名气)

[交流] 【求助】如何计算有机小分子的三重态能级(T1-S0)?

看了论坛里很多大牛关于类似问题的回复,但是还是不知道在Gaussian中如何具体操作。拿最简单的吡啶环为例来说吧:
1.先优化基态的最低能量构型:
# opt b3lyp/6-31+g(d)
Charge =  0 Multiplicity = 1
2.接着计算三重态的所谓垂直激发能:
# td=(triplets,nstates=3,root=1) b3lyp/6-31+g(d)
Charge =  0 Multiplicity = 1
3.然后优化此三重态的最低能量构型:
# opt b3lyp/6-31+g(d)
Charge =  0 Multiplicity = 3
问题是第二步计算了三个三重态,此步优化的是不是这三个三重态呢?还是应该在第二步中只计算一个三重态,然后在此步中优化此三重态也就是我们需要的T1态呢?还有这一步中自旋多重度为3没有错误吧?引申一下,若是在第二步中计算了十个单重态(也即十个吸收峰,假设它们均存在),在此步中优化单重态的最低能量构型如何计算,自旋多重度也是3对吗?
4.最关键一步,就是计算从最低能量构型的三重态回落到基态的过程,不知道如何计算?
我知道对各位大牛来说,这是个最简单计算发射光谱(不论是荧光光谱还是磷光光谱)的事情,原理也就是F-C原理搞清楚了,计算应该没有问题,可是对我们菜鸟来说,这些具体的操作才是问题的关键,希望各位大牛能够耐心指点一下,而不是一句F-C原理了事,谢谢大家了!
另外相关问题是:在计算类似问题时,我们这些非专业的人也就是用普通带XP的计算机算算而已,分子中要是有3到7个环左右(而这在有机合成中是非常常见的)时,可能光优化构型就得好几天,这中途要是碰上什么事情(停电、误操作等)停机了,就只能从头开始算,有没有什么办法能够接着上次终止的地方继续计算的,如何操作?我知道对各位专业人士来说这种事情无法想象,可国情如此,只好这么不专业地蛮干了。
回复此楼

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

分子模拟要术 不懂计算 量子化学-第一性原理 Gaussian
量化学者

» 猜你喜欢

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

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

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

jiewei

专家顾问 (著名写手)

小才

引用回帖:
Originally posted by cjhztt at 2009-11-8 20:13:
收益匪浅,xtdut算的(T1-S0)发射光谱应该是磷光吧,我想请教一下,计算一个分子,如何确定应该算荧光还是磷光啊,也就是说从何得知一个化合物是有荧光性质还是磷光性质?
望不吝赐教

你看它最后结果出来的谐振强度,比如说你算磷光的时候谐振强度为0那么这个分子具有荧光性质
我想做个科学家,可我又喜欢写诗;我想做一个诗人,可我却读的理工科..
8楼2009-11-08 23:16:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

ZDBWHZ

金虫 (正式写手)


xtdut(金币+1,VIP+0):谢谢您的回复,麻烦您能不能稍微耐心点教教怎么直接读取文件的最后一个结构计算?读取哪一个文件,是.out文件,还是.chk文件,亦或是.rwf文件呢?接着怎么知道那个是最后一个结构呢?最后是用GView打开算还是用GUASSIAN打开算?知道对您这种高手这不是问题,可咱是个菜鸟,麻烦您能不能说详细点,非常感谢! 11-6 13:20
关于停电的事情,直接读取文件的最后一个结构计算啊。
2楼2009-11-06 11:52:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuy069

铁杆木虫 (正式写手)

★ ★ ★ ★ ★
xtdut(金币+5,VIP+0):非常感谢您的回复! 11-6 13:11
首先你的3个步骤都是正确的。
第二步计算td,gaussian默认三个态,如果是对称性高的分子建议用 Largest Abelian subgroup 后的Nop 个数。为什么计算一个态需要用3个态或更多的态呢,那是因为求解这个态,需要用其它的态来对它就解。第2步,自旋多重度还是写1,因为垂直激发是基于基态的构型,所以电荷和自旋多重度都是基态的。
第3步,优化激发态的构型,电荷和自旋多重度当然需要激发态的,所以电荷和自旋多重度为0 3.
要想计算有机小分子的三重态能级(T1-S0),即是算激发态的发射,但是gaussian中没这么直接的,我们可以转换一下思路:要算发射,基于F-C原理,这时用优化好的激发态的构型,但是电荷和自旋多重度写0 1,那这时的构型相当于垂直时的基态构型,这时再做td,从输出文件中找到Excited State   1:      Triplet-A      2.6936 eV  460.28 nm  f=0.0000  =2.000 那么2.6936 eV  就是你所要的激发能。
3楼2009-11-06 12:06:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xtdut

银虫 (小有名气)

谢谢回复,金币已送出!

引用回帖:
Originally posted by wuy069 at 2009-11-6 12:06:
首先你的3个步骤都是正确的。
第二步计算td,gaussian默认三个态,如果是对称性高的分子建议用 Largest Abelian subgroup 后的Nop 个数。为什么计算一个态需要用3个态或更多的态呢,那是因为求解这个态,需要用其 ...

太感谢您了,还是有点不太明白,麻烦您耐心给指导一下:
也就是说第四步应该这样写:
# td=(singlets,nstates=3,root=1) b3lyp/6-31+g(d)
0 1
还是应该这样写:
# td=(triplets,nstates=3,root=1) b3lyp/6-31+g(d)
0 1
不太明白此时该计算单重态的发射还是三重态的发射?另外还想请教您一下,我在第三步优化时优化了三个三重态,还是只优化了第一个三重态(也即T1态),如果我想知道第二或第三个三重态的能级(也即T2或者T3时),是不是只需要在第四步中将root改为等于2或者3即可呢(当然nstates的值也应该相应地改大点)?
4楼2009-11-06 13:10:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 284求调剂 +21 让我上岸吧阿西 2026-04-09 21/1050 2026-04-15 21:53 by noqvsozv
[考研] 通信工程求调剂!!! +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 310求调剂 +16 666真好 2026-04-11 18/900 2026-04-15 13:28 by 黑科技矿业
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 26药学专硕105500求调剂 +6 喽哈加油 2026-04-13 6/300 2026-04-14 16:40 by zhouxiaoyu
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
[考研] 305求调剂 +8 玛卡巴卡boom 2026-04-11 8/400 2026-04-14 09:04 by pengliang8036
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 一志愿东北大学控制工程085406数二英二385,求调剂 +8 Ezra_Zhang 2026-04-09 8/400 2026-04-11 09:15 by 猪会飞
[考研] 309求调剂 +14 wdhw 2026-04-10 15/750 2026-04-10 21:06 by zhouxiaoyu
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
信息提示
请填处理意见