24小时热门版块排行榜    

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

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的回帖

yongleli

木虫 (正式写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by abraham89 at 2013-04-06 21:37:09
多谢回复。也就是说我要不停地尝试泛函和基组,然后选择比较适合的。那分子结构呢,我已经优化过了,还怎么调节呢...

1.换方法优化

2.根据实验结果手动设置
5楼2013-04-07 03:42:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qwertybobo

银虫 (小有名气)

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

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

你这法二略有一丝牛啊~
6楼2013-04-13 23:32:32
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 abraham89 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 17:31 by 4FFAWE8HcgUD
[基金申请] 面上意见出来了 +6 黄鸟于飞Chao 2026-08-29 12/600 2026-08-29 17:20 by 黄鸟于飞Chao
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 11:49 by jCd0dEvKHShX
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:40 by jCd0dEvKHShX
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +14 晴天加油 2026-08-26 15/750 2026-08-29 11:06 by biily
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
信息提示
请填处理意见