计算化学公社

标题: 求助，做BSSE修正 [打印本页]

作者
Author: 尚艳磊 时间: 2018-8-24 09:11
标题: 求助，做BSSE修正
我的体系是CH3OCH2CH3+NO2的体系，用的是G4MP2方法。审稿人的意见是让我做BSSE修正。但是我试过用关键词： g4mp2 counterpoise=2 出现了错误，如下：Electronic spatial extent (au):  <R**2>=          638.4379
Charge=             0.0000 electrons
Dipole moment (field-independent basis, Debye):
X=             1.2190 Y=             0.5225 Z=          -0.2919  Tot=             1.3580
Quadrupole moment (field-independent basis, Debye-Ang):
XX=          -23.7319 YY=          -24.1468 ZZ=          -26.6264
XY=             0.1498 XZ=          -0.1988 YZ=          -0.1286
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
XX=             1.1032 YY=             0.6882 ZZ=          -1.7914
XY=             0.1498 XZ=          -0.1988 YZ=          -0.1286
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=          -116.1076  YYY=          -5.4843  ZZZ=          -0.1770  XYY=          -38.9257
  XXY=          -3.6315  XXZ=             0.1891  XZZ=          -39.8883  YZZ=          -3.5138
  YYZ=          -0.0480  XYZ=          -0.0125
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=          -449.2186 YYYY=          -364.3927 ZZZZ=          -49.2588 XXXY=          18.8610
XXXZ=             1.1661 YYYX=          25.6754 YYYZ=          11.3076 ZZZX=          -0.8933
ZZZY=             8.6255 XXYY=          -133.7674 XXZZ=          -83.7661 YYZZ=          -69.0277
XXYZ=             4.9816 YYXZ=             0.1597 ZZXY=             9.7496
N-N= 1.334505959366D+02 E-N=-7.193087510918D+02  KE= 1.928719664286D+02
CPIOFr:  IOpCl= 1 IRwI=-1 IRwCP=  731 ICalc=  1 LCPTot=    190887 Len1MO=    629720 IndFrg=    1450327
LenMC>LenMMx in CPIOFr.
Error termination via Lnk1e in /home-lustre/shangyl/Gaussian/g16/l122.exe at Thu Aug 23 11:57:27 2018.
Job cpu time:    0 days  1 hours 32 minutes 42.0 seconds.
Elapsed time:    0 days  0 hours  5 minutes 22.4 seconds.
File lengths (MBytes):  RWF= 9864 Int=    0 D2E=    0 Chk=    19 Scr=    1

已经算出来的结果表明，在G4MP2的复合方法中B3LYP和CCSD是完成了的。但是MP2就失败了。
求助前辈们，我该怎么做？

作者
Author: yflchx 时间: 2018-8-24 10:18
算这种反应，不需要考虑BSSE吧？审稿人的出发点是什么?

作者
Author: 尚艳磊 时间: 2018-8-24 10:51

yflchx 发表于 2018-8-24 10:18
算这种反应，不需要考虑BSSE吧？审稿人的出发点是什么?

审稿人的出发点是觉得 reaction complex 需要把它的能量算的清楚

作者
Author: sobereva 时间: 2018-8-24 16:47
热力学组合方法根本不能和counterpoise联用，如果非要用，得自己手动做G4MP2的各个能量计算步骤，但原理上并不合适
另外，用这个热力学组合方法算复合物的热力学量也绝对不是恰当的选择，因为此方法用B3LYP做优化，而众所周知B3LYP不带色散校正时优化复合物是很不合适的

作者
Author: 尚艳磊 时间: 2018-8-24 19:30

sobereva 发表于 2018-8-24 16:47
热力学组合方法根本不能和counterpoise联用，如果非要用，得自己手动做G4MP2的各个能量计算步骤，但原理上 ...

是呀，Sob老师。我贴一段审稿人的话吧。
G4MP2 uses a small basis set for geometry optimisations which makes BSSE corrections necessary ... "G4MP2 theory is based on G4 theory$^6$. It uses B3LYP/6-31G(2df,p) optimized geometries for a series of single point energy calculations at higher levels of theory. In a previous study$^5$, this basis set was found to perform well for most geometries with some exceptions hydrogen bonded complexes and Jahn-Teller distorted systems."
他认为我需要做BSSE校正。但是G4MP2又不能做。我试了B3LYP/6-31G(2df,p)做BSSE。其实结果并不好啊。
Sob老师，给些建议吧

作者
Author: sobereva 时间: 2018-8-24 19:38

尚艳磊发表于 2018-8-24 19:30
是呀，Sob老师。我贴一段审稿人的话吧。
G4MP2 uses a small basis set for geometry optimisations whi ...

当前就根本不该用G4(MP2)，一方面没好处，另一方面不灵活
优化和单点都手动计算，干什么都方便
前面强调了，对于弱相互作用体系，用B3LYP超级不适合，当前问题关键不在于优化用的基组，在于理论方法烂

弱相互作用计算常识看
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413
DFT-D色散校正的使用
http://sobereva.com/210
乱谈DFT-D
http://sobereva.com/83
大体系弱相互作用计算的解决之道
http://sobereva.com/214
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

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