计算化学公社

标题: 计算电子亲和能（VEA）时遇到一些问题，求解答 [打印本页]

作者
Author: Marc0 时间: 2024-5-22 15:05
标题: 计算电子亲和能（VEA）时遇到一些问题，求解答
尝试用高斯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是否有必要加同样的关键词重算？

作者
Author: sobereva 时间: 2024-5-22 15:15
该说的都说了，没有可补充的
解决SCF不收敛问题的方法
http://sobereva.com/61

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

看具体什么关键词

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