计算电子亲和能（VEA）时遇到一些问题，求解答

Marc0

尝试用高斯09计算三苯基胺的电子亲和能，方法如下：
1. 优化分子构象，关键词：#p opt b3lyp/6-31g(d)
2. 优化完成后计算基态分子的单点能
3. 读取第一步优化后的.chk文件，直接另存为一个.com文件，计算该状态下阴离子的单点能。
步骤2和3均采用 m062x/ma-TZVP 的计算档次，仅仅是电荷数与自旋多重度的区别（前者是0 1后者是-1 2）。以步骤3为例，输入文件如下（参考了博文：给Ahlrichs的def2系列基组加弥散的方法 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)，用的是笨办法直接在输入文件中自定义基组）：


%nprocshared=10
%mem=20GB
#p m062x/gen


Anion opt


-1 2

C                 -0.33359036    1.38161140   -0.00120493
C                  0.37678851    2.28835252   -0.80487572
C                  0.05187099    3.64340114   -0.79289022
C                 -0.99450701    4.11690160    0.00161203
C                 -1.70877083    3.21631863    0.79468784
C                 -1.37913027    1.86232987    0.80385798
H                  1.18187626    1.92409366   -1.43511850
H                  0.61295750    4.33001607   -1.42160720
H                 -1.24966295    5.17277356    0.00270368
H                 -2.52109061    3.56982271    1.42453539
H                 -1.92872879    1.16970872    1.43341554
N                  0.00033833   -0.00021140   -0.00245012
C                 -1.02981670   -0.98000324   -0.00120511
C                 -2.17001076   -0.81723993   -0.80492650
C                 -0.92400406   -2.12610644    0.80360265
C                 -3.18175919   -1.77548907   -0.79312111
H                 -2.25621081    0.06213080   -1.43524015
C                 -1.93232198   -3.08801295    0.79410107
H                 -0.04932476   -2.25626114    1.43289534
C                 -3.06930871   -2.91862269    0.00110351
H                 -4.05684585   -1.63231657   -1.42181903
H                 -1.83288069   -3.96857599    1.42358219
H                 -3.85665624   -3.66700476    0.00197623
C                  1.36385934   -0.40211194   -0.00112356
C                  2.30386818    0.26506821    0.80125387
C                  1.79305906   -1.47298537   -0.80215929
C                  3.64123150   -0.12659699    0.79189078
H                  1.97942090    1.08925380    1.42857155
C                  3.12907320   -1.86909410   -0.79042607
H                  1.07470831   -1.98980088   -1.43054302
C                  4.06317836   -1.19763912    0.00133074
H                  4.35411598    0.40191076    1.41961199
H                  3.44254896   -2.70007412   -1.41711573
H                  5.10514888   -1.50468613    0.00225489



-H     0
S   3   1.00
     34.0613410              0.60251978E-02
      5.1235746              0.45021094E-01
      1.1646626              0.20189726
S   1   1.00
      0.32723041             1.0000000
S   1   1.00
      0.10307241             1.0000000
P   1   1.00
      0.8000000              1.0000000
****
-C     0
S   6   1.00
  13575.3496820              0.22245814352E-03
   2035.2333680              0.17232738252E-02
    463.22562359             0.89255715314E-02
    131.20019598             0.35727984502E-01
     42.853015891            0.11076259931
     15.584185766            0.24295627626
S   2   1.00
      6.2067138508           0.41440263448
      2.5764896527           0.23744968655
S   1   1.00
      0.57696339419          1.0000000
S   1   1.00
      0.22972831358          1.0000000
S   1   1.00
      0.95164440028E-01            1.0000000
P   4   1.00
     34.697232244            0.53333657805E-02
      7.9582622826           0.35864109092E-01
      2.3780826883           0.14215873329
      0.81433208183          0.34270471845
P   1   1.00
      0.28887547253          1.0000000
P   1   1.00
      0.10056823671          1.0000000
D   1   1.00
      1.09700000             1.0000000
D   1   1.00
      0.31800000             1.0000000
F   1   1.00
      0.76100000             1.0000000
S   1   1.00
  0.317214800093E-01   1.000
P   1   1.00
  0.335227455700E-01   1.000
****
-N     0
S   6   1.00
  19730.8006470              0.21887984991E-03
   2957.8958745              0.16960708803E-02
    673.22133595             0.87954603538E-02
    190.68249494             0.35359382605E-01
     62.295441898            0.11095789217
     22.654161182            0.24982972552
S   2   1.00
      8.9791477428           0.40623896148
      3.6863002370           0.24338217176
S   1   1.00
      0.84660076805          1.0000000
S   1   1.00
      0.33647133771          1.0000000
S   1   1.00
      0.13647653675          1.0000000
P   4   1.00
     49.200380510            0.55552416751E-02
     11.346790537            0.38052379723E-01
      3.4273972411           0.14953671029
      1.1785525134           0.34949305230
P   1   1.00
      0.41642204972          1.0000000
P   1   1.00
      0.14260826011          1.0000000
D   1   1.00
      1.65400000             1.0000000
D   1   1.00
      0.46900000             1.0000000
F   1   1.00
      1.09300000             1.0000000
S   1   1.00
  0.454921789167E-01   1.000
P   1   1.00
  0.475360867033E-01   1.000
****


我的问题是：
1. 计算步骤是否正确？所选择的计算档次是否合理？
2. 步骤3计算阴离子能量时报错，没有收敛。
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UM062X) =  -749.631641773     A.U. after  129 cycles
            NFock=128  Conv=0.15D-05     -V/T= 2.0056
参考博文：谈谈不同量子化学程序计算结果的差异问题 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)，加了关键词  SCF(novaracc,noincfock) int=acc2e=12三管齐下，结果仍然是不收敛。想知道还可以尝试什么办法？
3. 假设计算阴离子能量，通过加入一些关键词最终收敛，而步骤2计算基态能量时没有加关键词，两个能量之间可以直接相减作为VEA的值吗？步骤2是否有必要加同样的关键词重算？

sobereva

该说的都说了，没有可补充的
解决SCF不收敛问题的方法
http://sobereva.com/61

gjf内容别直接贴，直接作为附件上传，否则空行数目都乱了，本来Gaussian对空行就很敏感

看具体什么关键词