使用Gaussian09w刚性扫描Ac2分子键长获得解离曲线报错

小算符

各位老师：
       使用Gaussian09w对Ac2分子键长做刚性扫描，使用PBE0泛函，SARC-DKH2基组，基组定义从Basic set网站下载，使用DKH2处理相对论效应，使用了nosymm指令，键长从2.0A以步长0.1A扫描到5.0A，程序在l1模块出现2070错误。请问为什么会出现这种情况，如何解决？（输入文件和输出文件已上传至附件）

sobereva

r1 2.000 30 0.1前头必须空一行，仔细看
详谈使用Gaussian做势能面扫描
http://sobereva.com/474（http://bbs.keinsci.com/thread-12660-1-1.html）

zjxitcc

补充一点，这么做计算，结果基本上是完全错误的。简单来说有两点错误：（1）使用了RDFT；（2）扫描方向从短到长。（除非你是故意想获得一条错误的RPBE0解离曲线，你手上有另一条非常准确的解离曲线，想用RPBE0的错误曲线来衬托你高精度计算结果合理）。

你需要在键长最远的地方算一个UDFT单点，检验波函数稳定性，看自旋布居是否正确，必要时使用片段组合波函数做SCF初始猜测。然后读取这个结构的波函数 往键长缩短的方向做刚性扫描。如果做更细心一点的计算，双原子分子体系通常的做法是每个键长下的结构单独计算，都尝试片段组合波函数、检验波函数稳定性。