24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1949  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上和省基金B类撒花 +17 花田半亩~白 2026-04-21 17/850 2026-04-23 07:30 by john198197
[考研] 有没有学校收留 +3 蒋昌鹏qtj 2026-04-20 3/150 2026-04-22 20:25 by 学员JpLReM
[考研] 312求调剂 +3 山河似你温柔 2026-04-22 3/150 2026-04-22 20:17 by 学员JpLReM
[考博] 华师大读博 +3 xq83 2026-04-22 5/250 2026-04-22 10:42 by xq83
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +5 yexuqing 2026-04-19 5/250 2026-04-22 10:01 by easeheart
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
请填处理意见