计算化学公社

标题: 结构优化后单点能计算使用混合基组结果异常 [打印本页]

作者
Author: 天东耳日 时间: 2020-11-30 20:32
标题: 结构优化后单点能计算使用混合基组结果异常
各位老师，我在计算离子液体与汞的相互作用，结构优化时使用pbe0/6-311g**+SDD赝势，结构优化后没有虚频，能量为-1127.6966295 au。结构优化后再进行单点能计算，用混合基组，对阳离子与汞用def2TZVP，阴离子用ma-def2，弥散定义拷贝自340文，计算没有报错，能量为-4813.9520782，且出现大量虚频。请各位老师帮忙看看，是不是混合基组输入哪里写错了，应该怎么改正。结构优化输出文件，能量输入输出文件在附件中。PS：我计算单独的离子液体（没有加汞）单点能时用的混合基组写法与这个是一致的，计算结果就没有问题。

作者
Author: biogon 时间: 2020-11-30 20:37
可能是贋势未读入

作者
Author: zjxitcc 时间: 2020-11-30 21:49
def2-TZVP对Hg有赝势的，你既没写读赝势的关键词，也没写赝势数据，当然结果是错的了。

作者
Author: RAL 时间: 2020-11-30 21:53
感觉单点能计算的基组输入没有问题，可能结构优化的时候写法有问题所以Hg的部分基函数没有读入，从输出文件中看不出来。
PS1: 可以用#p获得更详细的输出。
PS2: 更换基组后势能面形状会改变，低水平下的极小点结构在高水平下一般就不是极小点了，此时计算频率徒增耗时，所得热力学结果没有意义。

作者
Author: 天东耳日 时间: 2020-12-1 08:55
谢谢各位老师，是赝势没有读入，重写了输入文件就可以了

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