计算化学公社

标题: 求助：电子亲和势和垂直电离势计算后得到的值为0，请问是什么原因？ [打印本页]

作者
Author: Trivers 时间: 2021-2-2 14:12
标题: 求助：电子亲和势和垂直电离势计算后得到的值为0，请问是什么原因？
各位大神：
   我参考了计算化学公社和思想家公社门口的博文后（下附链接），计算了一个中性有机分子的电子亲和势和垂直电离势，结果发现算出来的电子亲和势和垂直电离势都为0，请问是什么原因导致的，应该如何调整呢？
   我的计算方法如下：
   1、优化中性分子基态下的结构并计算单点能。由于分子不大且易于优化，所以直接用了B3LYP-D3BJ搭配def2TZVP基组进行优化，优化得到的能量直接作为单点能。同时为了得到热力学方面的能量，我也同时做了频率的计算。关键词如下：
         #p opt freq B3LYP/def2TZVP EM=GD3BJ scale=0.9590
   2、用优化后的结构进行电子亲和势和垂直电离势的计算。以电子亲和势为例，关键词如下：
         #p B3LYP/def2TZVP EM=GD3BJ geom=allcheck
         电荷量和自旋多重度分别为-1 和 2
   3、比较步骤1和步骤2得到的能量（log文件中搜索HF，取最末尾一个HF的值），结果发现二者完全相同。

      参考链接：
      http://bbs.keinsci.com/thread-1423-1-1.html
      http://bbs.keinsci.com/thread-6516-1-1.html
      http://bbs.keinsci.com/thread-1390-1-1.html
      http://sobereva.com/543

作者
Author: Trivers 时间: 2021-2-2 14:55
我好像知道是什么原因了……因为计算电子亲和势的时候用来geom=allcheck，所以读取了基态优化结果中包括电荷量、自旋多重度在内的所有信息并进行单点计算，所以得到的值和基态优化的能量必然完全相同……

作者
Author: snljty 时间: 2021-2-2 15:07

Trivers 发表于 2021-2-2 14:55
我好像知道是什么原因了……因为计算电子亲和势的时候用来geom=allcheck，所以读取了基态优化结果中包括电 ...

顺便，对于已经优化好的结构，D3色散校正加不加对电离能和电子亲合能计算没有影响。计算单点的水平通常建议比优化高一些。而且电子亲合能不加弥散函数用这种方法一定算不准。

作者
Author: Trivers 时间: 2021-2-2 16:28

snljty 发表于 2021-2-2 15:07
顺便，对于已经优化好的结构，D3色散校正加不加对电离能和电子亲合能计算没有影响。计算单点的水平通常建 ...

谢谢您的回复，所以是不是在计算涉及电子亲和势的内部重整化能的过程中，所有的单点能都应该用弥散函数进行计算？

作者
Author: sobereva 时间: 2021-2-3 05:09

Trivers 发表于 2021-2-2 16:28
谢谢您的回复，所以是不是在计算涉及电子亲和势的内部重整化能的过程中，所有的单点能都应该用弥散函数进 ...

是

如果是有机体系，用B3LYP算非常不理想。至少应当用M06-2X/ma-TZVP算单点

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