24小时热门版块排行榜    

查看: 928  |  回复: 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的回帖

461018249

木虫 (小有名气)

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

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

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode +8 documentary 2026-08-10 8/400 2026-08-10 19:20 by loufangrui
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 为什么网上很多人说本周 12号出结果 +3 瞬息宇宙 2026-08-10 3/150 2026-08-10 18:27 by Tide man
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +14 Lanmanbaby 2026-08-09 22/1100 2026-08-10 14:57 by wwncly
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见