| 查看: 1916 | 回复: 6 | ||
[求助]
求助有关高斯计算激发态能量准确度问题
|
|
最近在用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 |
回复此楼» 猜你喜欢
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有117人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
分子动力学模拟合作需求
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教：高斯如何计算单重态（激发态）过渡态
已经有7人回复
- 急求用高斯09怎么计算电子激发态
已经有7人回复
- 从高斯激发态优化输出文件中提取原子坐标
已经有9人回复
- 高斯计算激发态并行问题
已经有6人回复
- [求助]：高斯中的CIS是否能够优化开壳层的激发态结构，例如双重态，四重态的激发态
已经有8人回复
- 用高斯09计算的tddft方法优化激发态出现问题如何解决？
已经有15人回复
- 【求助】高斯09 激发态过渡态问题
已经有9人回复
- 【求助】高斯09优化激发态几何后算频率出错。求教改正方案。
已经有4人回复
- 【求助】请教高手，如何用高斯09的tddft方法计算激发态
已经有7人回复
- 【求助】高斯菜鸟求助激发态问题
已经有24人回复
- 【求助】单重激发态的高斯算法
已经有29人回复
stalart
木虫 (正式写手)
- 应助: 284 (大学生)
- 金币: 1541.4
- 红花: 9
- 帖子: 488
- 在线: 129小时
- 虫号: 1986267
- 注册: 2012-09-09
- 专业: 理论和计算化学
2楼2013-04-06 08:15:24
yongleli
木虫 (正式写手)
- 应助: 93 (初中生)
- 金币: 3769.1
- 散金: 369
- 红花: 24
- 帖子: 764
- 在线: 101.3小时
- 虫号: 303595
- 注册: 2006-12-03
- 专业: 原子和分子物理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
abraham89: 金币+10, ★★★★★最佳答案 2013-04-07 21:46:42
gmy1990: 金币+2 2013-04-07 23:05:20
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
4楼2013-04-06 21:37:09
yongleli
木虫 (正式写手)
- 应助: 93 (初中生)
- 金币: 3769.1
- 散金: 369
- 红花: 24
- 帖子: 764
- 在线: 101.3小时
- 虫号: 303595
- 注册: 2006-12-03
- 专业: 原子和分子物理
5楼2013-04-07 03:42:37
qwertybobo
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 280.7
- 散金: 18
- 帖子: 201
- 在线: 117.8小时
- 虫号: 915153
- 注册: 2009-11-28
- 专业: 高分子物理与高分子物理化
6楼2013-04-13 23:32:32
7楼2013-04-14 20:24:29
