24小时热门版块排行榜    

查看: 2307  |  回复: 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 德尚中行 2026-08-27 5/250 2026-08-29 16:16 by manplx
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 14:23 by jCd0dEvKHShX
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 13:51 by jCd0dEvKHShX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:40 by jCd0dEvKHShX
[基金申请] 中青基了要发朋友圈吗? +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见