求助，做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就失败了。
求助前辈们，我该怎么做？

yflchx

算这种反应，不需要考虑BSSE吧？审稿人的出发点是什么?

尚艳磊

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

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

sobereva

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

尚艳磊

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老师，给些建议吧

sobereva

尚艳磊 发表于 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