24小时热门版块排行榜    

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

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的回帖
查看全部 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-31 8/400 2026-09-02 09:37 by 昵称告罄
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +3 瞬息宇宙 2026-09-01 5/250 2026-09-02 08:01 by anjeeshine
[文学芳草园] 梦想 +5 myrtle 2026-08-26 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
信息提示
请填处理意见