24小时热门版块排行榜    

查看: 5966  |  回复: 29

lltyy

金虫 (小有名气)

[求助] 请教各位计算牛人有关HOMO/LUMO能极差的问题,万急!

小弟希望模拟紫外吸收光谱,即在基态优化的结构基础上算了一次TD-DFT,基态优化的log文件在GV中可以看到各个能级的能量,E(LUMO)- E(HOMO)能得到一个能极差。但是在TD计算后的log文件中,也能找到HOMO->LUMO的一个跃迁,也显示一个能量!这两个能量差不一样!有没有哪位大侠能帮我解释一下这是为什么啊!不理解TD-DFT的内涵!毕业论文万急!小弟先拜谢了!

[ 来自科研家族 木虫量化 ]
回复此楼

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

Gaussian 碳材料及二维材料 计算化学

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
lltyy: 金币+5 2012-05-07 12:52:37
引用回帖:
4楼: Originally posted by lihb734 at 2012-05-07 08:55:15:
首先确保你的计算方案一致的,否则那只能是见鬼了。

完全赞同lihb的言论
先检查你的计算基组一样吗?再检查你的计算方法一致吗?还有你的分子构型没变吧?
如果三者回答都是yes,那么你可以写信给高斯公司骂人了
5楼2012-05-07 09:37:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

purplez

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
lltyy: 金币+5 2012-05-08 21:34:36
引用回帖:
10楼: Originally posted by zhangmt at 2012-05-07 15:14:15:
简单的说说,关于Gaussian计算中输出的LUMO轨道的能量(LUMO波函数的本征值),有好多误解在里面。因为一般而言,大家会从Koopman定理的角度认知,认为轨道能量即是以电子距核无穷远处为能量零点,则电子轨道的能 ...

9楼已经很专业了
简单的说就是:DFT的基态HOMO,LUMO和常见的轨道能量定义其实有很大的差别(其实HF的轨道定义更接近我们一般理解的轨道能量)

另外,基态用DFT计算的HOMO-LUMO只能定性的解释反应能力,和激发态的吸收有关,但是能量上并不同(否则,要TD-DFT作甚?)

你就直接去TD-DFT的结果去和实验结果对比就好了

» 本帖已获得的红花(最新10朵)

23楼2012-05-08 20:42:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

procedure

铜虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
lltyy: 金币+5 2012-05-06 13:24:46
TDDFT计算的HOMO到LUMO跃迁能量为垂直跃迁能量,而DFT得到homo能量和Lumo能量差与垂直跃迁得到的能量不一致
2楼2012-05-06 13:12:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

blueybz

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
lltyy: 金币+5 2012-05-08 13:08:34
基态优化的log文件在GV中可以看到各个能级的能量,E(LUMO)- E(HOMO),这个是在HOMO有两个电子,而LUMO没有电子的情况下的能量差值。TD计算后的log文件中,HOMO->LUMO的一个跃迁,是在HOMO有1个电子,LUMO一个电子情况下的差值。
15楼2012-05-08 09:14:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

lltyy

金虫 (小有名气)

引用回帖:
2楼: Originally posted by procedure at 2012-05-06 13:12:39:
TDDFT计算的HOMO到LUMO跃迁能量为垂直跃迁能量,而DFT得到homo能量和Lumo能量差与垂直跃迁得到的能量不一致

照这么说我怎么觉得垂直跃迁能量应该更大啊!但是实际上TDDFT算出来的比较小!
3楼2012-05-06 13:24:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
lltyy: 金币+5 2012-05-07 12:52:42
首先确保你的计算方案一致的,否则那只能是见鬼了。
前途光明,出路难觅!
4楼2012-05-07 08:55:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lltyy

金虫 (小有名气)

引用回帖:
4楼: Originally posted by lihb734 at 2012-05-07 08:55:15:
首先确保你的计算方案一致的,否则那只能是见鬼了。

# opt=tight ub3lyp/genecp scrf=(cpcm,solvent=thf) geom=connectivity



C N O H Cl 0
6-31G(d,p)
****
Re 0
lanl2dz
****

Re 0
lanl2dz
------------------------------------------------------------------------------------
# td=(singlets,nstates=80) ub3lyp/genecp scrf=(cpcm,solvent=thf)
geom=connectivity



C N O H Cl 0
6-31G(d,p)
****
Re 0
lanl2dz
****

Re 0
lanl2dz

这计算方案应该是一致的吧?求指教!
6楼2012-05-07 12:59:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lltyy

金虫 (小有名气)

引用回帖:
5楼: Originally posted by yjcmwgk at 2012-05-07 09:37:51:
完全赞同lihb的言论
先检查你的计算基组一样吗?再检查你的计算方法一致吗?还有你的分子构型没变吧?
如果三者回答都是yes,那么你可以写信给高斯公司骂人了

还真是都一样!不懂了!
7楼2012-05-07 13:02:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

【答案】应助回帖


御剑江湖: 金币+1, 谢谢,系统太呆笨了,呵呵 2012-05-08 19:53:52
电荷多重度给一下
“emuch真变态,只给上述回复还不让”
前途光明,出路难觅!
8楼2012-05-07 14:24:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

引用回帖:
7楼: Originally posted by lltyy at 2012-05-07 13:02:33:
还真是都一样!不懂了!

太好了
写信骂高斯公司吧
9楼2012-05-07 14:28:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

★ ★ ★ ★ ★
zhou2009: 金币+5 2012-05-08 17:41:32
简单的说说,关于Gaussian计算中输出的LUMO轨道的能量(LUMO波函数的本征值),有好多误解在里面。因为一般而言,大家会从Koopman定理的角度认知,认为轨道能量即是以电子距核无穷远处为能量零点,则电子轨道的能量即等于将电子从处于该轨道移到无穷远处的能量差值。
然而,事实上,Gaussian在计算LUMO轨道能量时,由于该处事实上并不存在电子,因此也就不存在该处那个“虚拟”的“电子”对于内层(或同层)电子的排斥作用,因此,如果事实上真的将一个电子从无穷远处移到此轨道,则新加入的电子对内层电子的排斥同样导致内层电子的轨道能量发生变化,内层电子会因为多了外加的排斥作用而更加贴近内核,势能降低而动能增加,而内层电子对新加的电子的排斥则小于原先预计的值。整个分子体系的能量变化绝不仅仅是LUMO的轨道能。换到楼主的激发态的体系中,内层一个电子激发到高层的轨道,比如LUMO或LUMO+n等等,由于内层离去个电子,原本的电子之间的排斥减小了,内层电子的结构会变化得更加松散,内层的各个轨道的能级也会变化,相应地,对于要处于LUMO层的那个激发电子的排斥也会增加,对于原本闭壳层的分子而言,激发一个电子同时还损失一个配对能。所以,整个分子体系能量的变化也绝不仅仅是LUMO-HOMO那么简单——不然的话,还要TD计算干什么?
所以,LUMO轨道的能量值,仅仅就是个示意,定性地看看,在有限的化学经验而言,才有一点点意义。凭它是无法处理定量的光谱计算的。。
类似地推论一下,即可知道,HUMO轨道能量值不等于第一电离能,LUMO轨道能量值也不等于电子亲和势。
以上是个人理解,说错的地方,请专家指点。

» 本帖已获得的红花(最新10朵)

一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
10楼2012-05-07 15:14:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lltyy 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601求调剂 +3 Du.11 2026-03-16 3/150 2026-03-16 20:42 by 无际的草原
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
信息提示
请填处理意见