24小时热门版块排行榜    

查看: 2300  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

abraham89

新虫 (初入文坛)

[求助] 求助有关高斯计算激发态能量准确度问题

最近在用Gaussian 09计算一些激发态,为了验证计算准确度,我算了丙酮的单重态和三重态,算出来E(S1)=426kJ/mol, E(T1)=366kJ/mol, 而一些书籍中这两个值分别为355和326,差距有点大。请问各位高手,我的计算有问题吗?多谢!
我的输入文件如下:

%chk=acetone MD.chk
%nproc=2
# b3lyp/6-311+g(d,p) opt freq

acetone

0 1
C                  1.29068400   -0.61278400   -0.00415300
C                  0.00002500    0.18521400   -0.00010300
O                  0.00000000    1.39672400    0.00001800
C                 -1.29069600   -0.61279700    0.00411500
H                  1.30715800   -1.32008000    0.83447400
H                  2.14646800    0.06086100    0.05926800
H                  1.35991700   -1.20668200   -0.92622400
H                 -2.14639800    0.06098500   -0.06148200
H                 -1.36077400   -1.20459700    0.92773200
H                 -1.30644400   -1.32207800   -0.83306700
1 2 1.0 5 1.0 6 1.0 7 1.0
2 3 2.0 4 1.0
3
4 8 1.0 9 1.0 10 1.0
5
6
7
8
9
10

--Link1--
%chk=acetone MD
%nproc=2
# b3lyp/6-311+g(d,p) TD=(singlets, Nstates=3, Root=1) Geom=check Guess=read

acetone singlets

0 1

--Link1--
%chk=acetone MD
%nproc=2
# b3lyp/6-311+g(d,p) TD=(Triplets, Nstates=3, Root=1) Geom=check Guess=read

acetone triplets

0 1
回复此楼

» 猜你喜欢

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

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

zyhappy

新虫 (小有名气)

引用回帖:
5楼: Originally posted by yongleli at 2013-04-07 03:42:37
1.换方法优化

2.根据实验结果手动设置...

您是传说中的李永乐?
7楼2013-04-14 20:24:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

stalart

木虫 (正式写手)

【答案】应助回帖


gmy1990: 金币+1 2013-04-07 23:05:35
计算值与文献实验值之间的差距实际上就是拟合的好坏,通过尝试不同的基组来拟合,验证,最终预报。
2楼2013-04-06 08:15:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongleli

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
abraham89: 金币+10, ★★★★★最佳答案 2013-04-07 21:46:42
gmy1990: 金币+2 2013-04-07 23:05:20
1.TDDFT的结果,依赖于泛函和结构。
结构的改变对结果影响比泛函大。

2.如果结构已经用CCSD(T)-f12优化好了,
并且利用Monte Carlo之类的办法确定了优化到的确实是全局能量最小点,
则Nstates越多越精确。

3.楼主这个case,槽点很多。
1. 结构优化了吗?精度够吗?
2.大家伙都知道,B3LYP算TDDFT误差很大。而且TDDFT的理论误差
达到0.3eV,也就是调节泛函、开大基组,仍然可能有
7 kcal/mol的误差。楼主的case里差别在十几个kcal/mol,
属于正常误差范围之内。如要精确计算光谱,需要先做一番
test,再选出适合自己的体系的泛函和基组的组合。
不是随随便便就能成功的。
当然,楼主也可以手动调一下分子结构,直到TDDFT计算出来的数值
和文献值符合为止。
3楼2013-04-06 09:21:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abraham89

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by yongleli at 2013-04-06 09:21:55
1.TDDFT的结果,依赖于泛函和结构。
结构的改变对结果影响比泛函大。

2.如果结构已经用CCSD(T)-f12优化好了,
并且利用Monte Carlo之类的办法确定了优化到的确实是全局能量最小点,
则Nstates越多越精确。

...

多谢回复。也就是说我要不停地尝试泛函和基组,然后选择比较适合的。那分子结构呢,我已经优化过了,还怎么调节呢
4楼2013-04-06 21:37:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我不理解! +5 Edward_pc 2026-08-26 8/400 2026-08-26 15:21 by climap
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 系统查不到 +9 董八千 2026-08-26 9/450 2026-08-26 12:30 by yanmofang
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 范进中举一文的中心思想 +8 炎黄贵胄 2026-08-22 9/450 2026-08-26 07:47 by WASM
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见