计算化学公社

标题: 使用Gaussian09w刚性扫描Ac2分子键长获得解离曲线报错 [打印本页]

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

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

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

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

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