24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2674  |  回复: 8

kewei09

银虫 (小有名气)

[交流] 【求助】td计算三重态问题已有1人参与

用td算三重态,过程没有出错但是每个态的振子强度f都是0,在算单重态的时候就没有这个现象,总觉得输入有什么问题,以下是输入文件。
#P b3lyp/6-31+g* td(direct,triplets,nstates=10) test pop=full scrf=(pcm,read)

输出文件里
Excited State   1:   Triplet-A      1.1271 eV 1100.00 nm  f=0.0000
Excited State   2:   Triplet-A      1.9968 eV  620.92 nm  f=0.0000
……
十个态都是0啊

另外问一下如果我想计算发射光谱,除了上面这个方法还有什么好方法可以试一下

[ Last edited by lei0736 on 2009-11-13 at 22:00 ]
回复此楼

» 猜你喜欢

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

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

Jasminer

铁杆木虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
正常,f=0说明没有吸收
这不是算发射光谱
2楼2009-09-17 19:12:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
loovfnd(金币+5,VIP+0):作图的能力太强了 9-17 21:06
你的输入没错,但是你的思想有错
楼主很显然对td计算根本就没有好好研究,就匆匆忙忙开始了计算。
从稳定的基态,垂直跃迁到单重态的第一激发态,这是可以的。因此“在算单重态的时候就没有这个现象”。你这是算的吸收谱。
从稳定的基态,一般来说很难直接跃迁到第一激发态的叁重态。因此“每个态的振子强度f都是0”。你这算的什么谱都不是。
如果要得到发射谱,应该从第一激发态的单重态的最低能态,垂直跃迁到基态。此时基态的电子结构发生了变化,而分子结构尚未发生变化。请自己研究一下“垂直跃迁规则”
要用td算发射谱,你要计算第一激发态的单重态的最低能态。然后用第一激发态的单重态的最低能态的构型,计算td

[ Last edited by yjcmwgk on 2009-9-17 at 19:41 ]
3楼2009-09-17 19:39:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
erylingjet(金币+2,VIP+0):感谢交流! 9-18 09:31
这个问题和吸收、发射没有关系。

在非相对论近似下,单态-三态的跃迁是禁戒的(这是个常识,可以读任何一本分子光谱或量子力学的教科书),故f=0。

标量相对论也是如此。

只有考虑自旋-轨道耦合,f才可能不为0。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2009-09-17 22:46:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

psfan

铁虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
erylingjet(金币+1,VIP+0):感谢交流 9-18 09:31
引用回帖:
Originally posted by beefly at 2009-9-17 22:46:
这个问题和吸收、发射没有关系。

在非相对论近似下,单态-三态的跃迁是禁戒的(这是个常识,可以读任何一本分子光谱或量子力学的教科书),故f=0。

标量相对论也是如此。

只有考虑自旋-轨道耦合,f才可能 ...

对于lz的问题,如果考虑自旋-轨道耦合该用什么方法算呢?
5楼2009-09-18 07:27:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
erylingjet(金币+2,VIP+0):感谢交流! 9-18 09:31
比较容易的方法是用adf的so-tddft。

也可以用gamess的casscf或mcqdpt2做旋轨耦合计算,不过输入更复杂,而且计算量很大。molcas,molpro和columbus等程序也有类似的功能。我推荐用molpro,因为输入最简单,计算最快。不论选哪个,用之前一定要懂点群论(除非分子是C1群)。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
6楼2009-09-18 07:51:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

psfan

铁虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
目前理论水平下算过渡金属等重原子的旋轨耦合效应结果如何?和实验值符合的好么?
7楼2009-09-18 12:05:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abbott

金虫 (著名写手)

不要用QQ问我东西

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+1,VIP+0): 9-23 21:19
triplets,nstates=10

lz,  计算得到的f=0,是由于算法的问题,
g03不能求解一个矩阵方程导致的.....

比如ADF 就能计算三重态的F值
Chemistry[]==[]Chem[]is[]try!!!
8楼2009-09-23 14:14:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

blueybz

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by yjcmwgk at 2009-09-17 19:39:30:
你的输入没错,但是你的思想有错
楼主很显然对td计算根本就没有好好研究,就匆匆忙忙开始了计算。
从稳定的基态,垂直跃迁到单重态的第一激发态,这是可以的。因此“在算单重态的时候就没有这个现象”。你这是算 ...

可能楼主用的是 优化好的三重态结构计算的,可以是磷光发射。f=0,因为跃迁是禁阻的。强度为0
9楼2011-04-03 00:55:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kewei09 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见