计算化学公社

标题: guess=read使用时报错chk文件无基组信息，是哪里出错了？ [打印本页]

作者
Author: ganlanchina 时间: 2024-11-16 20:53
标题: guess=read使用时报错chk文件无基组信息，是哪里出错了？
本帖最后由 ganlanchina 于 2024-11-25 21:26 编辑

在络合物三线态优化过程中使用关键词 #p opt freq def2svp em=gd3bj upbe1pbe无法收敛。
根据sob老师意见，使用小基组优化#p opt freq 6-31g(d) em=gd3bj upbe1pbe顺利收敛，得到chk文件另存为1116.chk；
后希望读取此优化下的波函数作为初猜，继续在def2svp下优化：
%chk=1116.chk
%nprocshared=32
%mem=100GB
#p opt freq def2svp em=gd3bj upbe1pbe guess=read

（输入文件名为1116.gjf）
然而运行很快报错：
Initial guess from the checkpoint file: "1116/1116.chk"
Basis set data is not on the checkpoint file.
Error termination via Lnk1e in /public/home/ganlan/apprepo/gaussian/16-hy/scripts/../app/g16/l401.exe at Sat Nov 16 18:51:11 2024.

输出文件1116.log附于附件。
本人有机化学出生，纯自己摸索了一个月的计算化学入门的，基础很不扎实。感觉应该是犯了很低级的错误。各位前辈，能否指点一下是哪里出错了，怎么改进？

作者
Author: ganlanchina 时间: 2024-11-16 21:47
def2svp下输入文件的坐标不是6-31g(d)优化完最后一帧的坐标
而是6-31g(d)那个优化的输入文件的坐标
这可能是报错的原因吗

作者
Author: sobereva 时间: 2024-11-17 02:22
从输出文件可见，明明你写的是%chk=1116/1116.chk，而不是你帖子里说的%chk=1116.chk
确保chk的路径写对了，并且确保确实chk文件是6-31G*任务直接产生的。
既然用def2-SVP进一步优化，显然结构也应当用6-31G*优化完的。应写上geom=allcheck从chk里面读取坐标

另外，6-31G*和def2-SVP大小几乎没多大差异，除非你还要给氢加极化（当前体系并不需要），否则直接用6-31G*的结果就完了

作者
Author: ganlanchina 时间: 2024-11-18 08:48
感谢老师回复，后面确认了一下路径确实有问题。感觉是超算slurm系统的客服帮忙编写的批量提交脚本有问题。

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