计算化学公社

标题: 求推荐高精度小分子计算方法 [打印本页]

作者
Author: gauss98 时间: 2017-4-12 14:30
标题: 求推荐高精度小分子计算方法
本人打算研究甲烷燃烧，乙烯高温裂解等问题
需要算大量自由基反应过程

我用了 B3LYP-D3， M062x -D3，配上 6-311G**, 6-311++g**, 6-311++(2d,p), def2tzvp, def2tzvpp等，用 CO， N2， O2， CH4分子算键长和原子化能来检验
总得不到一致的好结果，
求老师们推荐较好精度的计算方法和检验方法，计算需要求反应速度，谢谢！

作者
Author: 冰释之川 时间: 2017-4-12 16:27
这么小的分子是不是可以CCSD优化起来了

，算能量用CCSD(T)

作者
Author: nanski 时间: 2017-4-12 17:54
最终是为了算反应速率的话，还是要保证能量和自由能比较准就好，可以测试下Gn方法比如G3，G4，或者其他CBS方法，甚至是W1之类。
检验的话，就用一个你研究的物种的转化反应，用不同方法算一下比较下反应自由能垒，这应该是研究反应速率时最关心的问题吧。
这些小分子很多有数据可以参照，拿标准自由能做个对比也是个不错的办法
另外DFT如果结果和上述方法差别不大的话也是可以接受的，毕竟都是主族小分子，但是算自由能的时候一定要谨慎，时刻记得看一看每个结果是否合理。如果自由能不太好的话，可能要考虑非简谐矫正等办法

作者
Author: sobereva 时间: 2017-4-12 18:05
高精度计算就根本别考虑那些DFT，基组小于def2TZVPP也没有考虑的必要。

算热力学数据，算得动就用W1，是Gaussian支持的精度最高的热力学组合方法。算不动就降到CBS-APNO，还算不动就降到G4，还算不动就降到CBS-QB3。

作者
Author: gauss98 时间: 2017-4-12 18:15

sobereva 发表于 2017-4-12 18:05
高精度计算就根本别考虑那些DFT，基组小于def2TZVPP也没有考虑的必要。

算热力学数据，算得动就用W1，是 ...

谢谢！
有相关方法的输入文件模板么？ test之类的，以前没做过。

作者
Author: sobereva 时间: 2017-4-12 18:18

gauss98 发表于 2017-4-12 18:15
谢谢！
有相关方法的输入文件模板么？ test之类的，以前没做过。

# W1
就完了

作者
Author: gauss98 时间: 2017-4-12 21:15

sobereva 发表于 2017-4-12 18:18
# W1
就完了

CBS等方法都试过通过了，W1方法始终出现错误
Gaussian 09:  ES64L-G09RevD.01 24-Apr-2013
            12-Apr-2017
******************************************
%nprocshared=8
Will use up to 8 processors via shared memory.
----
# W1
----
QPErr --- An ambiguous keyword was detected.
# W1
'
Last state= "GCL"
TCursr=    1232 LCursr=       2
Error termination via Lnk1e in /home/gauss/g09/l1.exe at Wed Apr 12 20:23:00 2017.
Job cpu time:    0 days  0 hours  0 minutes  0.3 seconds.
File lengths (MBytes):  RWF=    12 Int=    0 D2E=    0 Chk=    2 Scr=    2

输入文件是

%nprocshared=8
# W1

N2

0 1
N                   0.       0.       0.
N                   0.       0.       1.12

计算时弹出Error: segmentation violation
rax 0000000000000000, rbx ffffffffffffffff, rcx ffffffffffffffff
rdx 0000000000007f62, rsp 00007ffffff5e538, rbp 00007ffffff5eab0
rsi 000000000000000b, rdi 0000000000007f62, r8  00002aaaaaabc180
r9  0000000000000000, r10 00007ffffff5e2c0, r11 0000000000000202
r12 0000000000000000, r13 0000000000000000, r14 00007ffffff5eaf8
r15 00000000000003e6
  --- traceback not available

作者
Author: sobereva 时间: 2017-4-12 22:18

gauss98 发表于 2017-4-12 21:15
CBS等方法都试过通过了，W1方法始终出现错误
Gaussian 09: ES64L-G09RevD.01 24-Apr-2013
...

我写错了，应该是W1U、W1BD、W1RO三种之一
用W1RO就可以，和W1的原始版本对应。W1BD和W1U都是后来的变体，但精度不增反降。

作者
Author: gauss98 时间: 2017-4-13 21:03

sobereva 发表于 2017-4-12 22:18
我写错了，应该是W1U、W1BD、W1RO三种之一
用W1RO就可以，和W1的原始版本对应。W1BD和W1U都是后来的变 ...

谢谢！
测试了下，CBS-APNO似乎不错

但是这些方法都是一系列组合方法，开始用优化，如何在该级别下求过渡态？

我先用 DFT优化了个过渡态，作为输入用CBS，但还是自动优化成反应物了。

谢谢

作者
Author: sobereva 时间: 2017-4-14 00:13

gauss98 发表于 2017-4-13 21:03
谢谢！
测试了下，CBS-APNO似乎不错

用热力学组合方法的时候同时用opt(TS,calcfc,noeigen)，这样自动优化时就会往过渡态优化

欢迎光临计算化学公社 (http://bbs.keinsci.com/)