计算化学公社

标题: 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/) Powered by Discuz! X3.3