查看: 2951  |  回复: 9
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

faqianliu

银虫 (小有名气)

[交流] 【求助】关于TD-DFT计算紫外光谱及采用分子轨道进行解释的几个问题

关于TD-DFT计算紫外光谱及采用分子轨道进行解释的几个问题,请高手指正:

部分问题选自Adv. Mater. 2008, 20, 1870–1875
论文解析了三个晶体的结构,以TD-DFT结论来解释紫外光谱信息,从而解释材料的三阶非线性。

以配合物1为例,其结构参见附件图1


附件图2为其紫外光谱曲线以及分子轨道解释
两条紫外光谱曲线分别为配合物1的实验值和使用gaussian03 TD-DFT计算的理论值

其中论文中peak1的量子化学解释为HOMO-2到LUMO的跃迁,并指出为(3ps至5dw的跃迁)

提供了部分TD-DFT结果,并给出了相应轨道的轨道系数,参见附图3


现有几个疑问请教各位:
问题1:根据excited state 8 的结果,peak1应为四种跃迁的组合,为何论文中只认定是101-104的跃迁,是不是因为101-104的跃迁系数最大?

问题2:判定为101-104跃迁后,根据后面的轨道系数,如何判定为3ps至5dw的跃迁?从101-104的轨道系数好像看不出S的3P轨道和W的5d轨道的特殊之处,是不是通过图形看的?

问题3:论文中计算的紫外光谱均处于500-200nm之间,为什么我采用TD-DFT计算的光谱范围为4000以下,其中有2XXXnm的,也有1XXXnm的,利用nstate=20才计算到500nm左右,利用nstate=50计算到300nm左右,汗死,我看该论文提供的第一激发态就处于500nm以下,能不能加入某些参数使gaussian只计算500-200nm之间的紫外光谱?

问题4:论文中称该计算是在DMF溶剂中进行,并给出了DMF的介电常数为36.7,但是如何在输入文件中体现,应该是scrf=(solvent=dmf 36.7)?请高手指正?最好给出命令行,非常感谢!

问题5:我计算的体系为开壳层多金属配合物体系,除避免在解释时使用HOMO及LUMO轨道外,在计算时有没有特别需要注意的问题?我计算的体系自旋多重度为2,计算结果如下:
Excited State   1:   ?Spin  -A      0.2640 eV 4695.86 nm  f=0.0000
    100B ->105B       -0.10199
    100B ->106B       -0.21311
    100B ->107B        0.80119
    100B ->108B        1.05220
    101B ->107B       -0.17708
    101B ->108B       -0.23320
其中所有的激发态为什么会出现问号?


本人为做晶体的,对量子化学不懂,刚接触gaussian半个月左右,看了论坛很多资料,这几个问题也无法解释,还请各位高手不要见笑,再谢!

[ Last edited by faqianliu on 2009-11-2 at 13:28 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuxiaocunde

铁杆木虫 (著名写手)

古德里安将军

图片看不到
显示是:
此图片来自QQ空间
未经允许不准使用
新征途
3楼2009-11-02 12:39:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

faqianliu

银虫 (小有名气)

lei0736(金币+0,VIP+0):还是看不到 居然放在QQ空间 11-2 11:38
终于把图片搞定了,哪位高手给指导一下

[ Last edited by faqianliu on 2009-11-2 at 11:23 ]
2楼2009-11-02 11:12:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

faqianliu

银虫 (小有名气)

引用回帖:
Originally posted by liuxiaocunde at 2009-11-2 12:39:
图片看不到
显示是:
此图片来自QQ空间
未经允许不准使用

重新整理了一下

[ Last edited by faqianliu on 2009-11-2 at 13:28 ]
4楼2009-11-02 13:07:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+5,VIP+0): 11-2 14:34
问题1:根据excited state 8 的结果,peak1应为四种跃迁的组合,为何论文中只认定是101-104的跃迁,是不是因为101-104的跃迁系数最大?
是,且平方以后相对更大。

问题2:判定为101-104跃迁后,根据后面的轨道系数,如何判定为3ps至5dw的跃迁?从101-104的轨道系数好像看不出S的3P轨道和W的5d轨道的特殊之处,是不是通过图形看的?
有经验者从轨道系数就可以看出来,当然从图形看更简单直接。

问题3:论文中计算的紫外光谱均处于500-200nm之间,为什么我采用TD-DFT计算的光谱范围为4000以下,其中有2XXXnm的,也有1XXXnm的,利用nstate=20才计算到500nm左右,利用nstate=50计算到300nm左右,汗死,我看该论文提供的第一激发态就处于500nm以下,能不能加入某些参数使gaussian只计算500-200nm之间的紫外光谱?
你的体系算出来是多少就是多少。至于有没有参数“使gaussian只计算500-200nm之间的紫外光谱”,同求……不过我没找到过。

问题4:论文中称该计算是在DMF溶剂中进行,并给出了DMF的介电常数为36.7,但是如何在输入文件中体现,应该是scrf=(solvent=dmf 36.7)?请高手指正?最好给出命令行,非常感谢!
命令行看说明书scrf一节。但该文献用的哪种溶剂模型应有具体说明,或许在supporting里面。

问题5:我计算的体系为开壳层多金属配合物体系,除避免在解释时使用HOMO及LUMO轨道外,在计算时有没有特别需要注意的问题?我计算的体系自旋多重度为2,计算结果如下:
Excited State   1:   ?Spin  -A      0.2640 eV 4695.86 nm  f=0.0000
    100B ->105B       -0.10199
    100B ->106B       -0.21311
    100B ->107B        0.80119
    100B ->108B        1.05220
    101B ->107B       -0.17708
    101B ->108B       -0.23320
其中所有的激发态为什么会出现问号?
开壳层体系我见过多种说法,有说gaussian的td不能正确处理的,有说应该选用RO基态波函数的,有说无影响的。我倾向于RO方法,但无论如何gaussian的td方法是不适于处理开壳层体系的,推荐软件orca、nwchem等。也有人说可以用一个闭壳层的激发态代替开壳层作为参考态进行td计算,具体怎么做、结果如何我不清楚。
问号是因为基态为开壳层时软件辨别不了激发后的自旋态。

[ Last edited by Jasminer on 2009-11-2 at 14:19 ]
5楼2009-11-02 13:44:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +6 Isak0320 2026-10-12 8/400 2026-10-12 20:09 by caixiaoye1
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见