24小时热门版块排行榜    

查看: 976  |  回复: 5

461018249

木虫 (小有名气)

[求助] 紫外计算

我计算紫外的步骤是
1 对分子结构进行优化 # B3LYP/6-31g* opt  freq
2 用上述优化后的分子结构 #p td rb3lyp/6-31+g(d) geom=connectivity
得到的应该就是紫外谱图了
得到的结果是
Excitation energies and oscillator strengths:

Excited State   1:   Singlet-A      2.1808 eV  568.53 nm  f=1.1775
     149 ->150         0.62205
This state for optimization and/or second-order correction.
Total Energy, E(RPA) =  -1658.19079878   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:   Singlet-A      2.7537 eV  450.24 nm  f=0.2086
     147 ->150        -0.17227
     148 ->150         0.62405
     149 ->151         0.15662

Excited State   3:   Singlet-A      2.9378 eV  422.04 nm  f=0.0061
     147 ->150         0.68151
     148 ->150         0.14129
问题1. 568.53 nm是不是就是最大吸收波长?如果是,为什么比实验值小了一百多纳米呀?反倒是568.53 nm * f=1.1775 的值和实验值比较接近 但是我没见过这样的计算方法呀  他们之间应该是正比的关系的呀  如果不是的话  那应该是什么
2. 我加大nstates=10,得出的结果和上面的差不多 就是多了几个states而已 到底什么样的结果是可靠的  还是我计算的问题
请各位高手赐教
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

【答案】应助回帖

感谢参与,应助指数 +1
461018249(金币+10): ★★★很有帮助 3 2012-02-16 15:29:06
b3lyp倾向于高估激发能(也就是低估激发波长)
td计算对泛函很敏感
可以采用bp86去计算
参考J Phys Chem A 2011, 115, 13811

还有啊,楼主居然打算把波长和振子强度乘到一块儿?搞啥妖蛾子呢?这么做有啥物理意义?理论计算不是瞎搞
振子强度是吸光率的标示。光子捕获效率(Light Harvesting Efficiency,LHE)可以通过这个公式计算LHE=1-[10^(-f)]
参考J Mater Chem 2012, 22, 568

话说你们单位有没有专门搞量子化学计算的人啊?你跟他合作一下吧。量子化学计算不是黑匣子计算。不懂这行儿的人,他们算的东西,十之八九不靠谱
2楼2012-02-16 13:43:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
本帖仅楼主可见
3楼2012-02-16 13:53:24
已阅   申请QC强帖   回复此楼   编辑   查看我的主页

gongxd325

至尊木虫 (文坛精英)

【答案】应助回帖

感谢参与,应助指数 +1
要是你的计算结果和实验结果完全一致,那肯定有问题
你算的是气态单个分子,实验测的可不是
一个人做点好事并不难难的是一辈子做好事!!!
4楼2012-02-16 14:26:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

461018249

木虫 (小有名气)

引用回帖:
: Originally posted by gongxd325 at 2012-02-16 14:26:41:
要是你的计算结果和实验结果完全一致,那肯定有问题
你算的是气态单个分子,实验测的可不是

谢谢你的回复  可是差了100多是不是太多了点
5楼2012-02-16 15:22:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

461018249

木虫 (小有名气)

引用回帖:
: Originally posted by yjcmwgk at 2012-02-16 13:43:05:
b3lyp倾向于高估激发能(也就是低估激发波长)
td计算对泛函很敏感
可以采用bp86去计算
参考J Phys Chem A 2011, 115, 13811

还有啊,楼主居然打算把波长和振子强度乘到一块儿?搞啥妖蛾子呢?这么做有啥物 ...

谢谢你 我看了你推荐了文章  怎么都说要计算这三个呢geometries, electronic structures and absorption spectra   第一个应该是基态构型吧  电子结构是什么构型  请赐教  还有里面有提到经过一个CAM-B3lyp的校正  这是怎么校正的 谢谢回复
6楼2012-02-16 15:28:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 461018249 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-29 21:35 by 4qO8zhfjr2AX
[基金申请] 面上意见出来了 +8 黄鸟于飞Chao 2026-08-29 16/800 2026-08-29 20:45 by 黄鸟于飞Chao
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-29 17:22 by 4FFAWE8HcgUD
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 14:03 by jCd0dEvKHShX
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:02 by jCd0dEvKHShX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 11:49 by jCd0dEvKHShX
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
信息提示
请填处理意见