计算化学公社

标题: 求助：sobMECP运行runMECP.sh报错，Problem with the ab initio Job. [打印本页]

作者
Author: 小学徒 时间: 前天 22:00
标题: 求助：sobMECP运行runMECP.sh报错，Problem with the ab initio Job.
各位老师好，最近在使用sobMECP程序寻找五重态到三重态的势能面交叉点（MECP）时遇到问题，想请教一下可能的原因和解决方案。
问题描述：
runfirst.sh 正常结束，log结尾显示Normal termination。
运行 ./runMECP.sh 后，程序报错，显示
“ProgFile ExistsFirst Input OK
Problem with the ab initio Job.
Problem with Fortran program
Step Number 0 -- MECP not yet converged
An error has occurred,possibly in the Gaussian Job”

ReportFile中只有 “ERROR” 一词，没有其他具体错误信息。之前运行过CH2的例子，可以MECP收敛。想问一下各位老师这可能是什么原因？如何进一步定位错误根源？如何解决？是否初猜不合适？Gaussian版本：G16 A.03，直接命令行运行。有关初猜结构选择：自旋态：五重态 → 三重态，选择了两个自旋态之间更稳定的过渡态TS的几何结构作为初猜。DFT计算表明，交叉点可能位于中间体I和II之间，是否将初猜结构选为两个自旋态之间更稳定的中间体II?
感谢各位老师指教！

作者
Author: sobereva 时间: 前天 22:16
没你的输入文件没法进行测试
sobMECP的脚本是完全透明的，直接从源代码里搜大写的ERROR结合实际情况就能判断是因为什么触发的。跟初猜结构没关系

作者
Author: 小学徒 时间: 前天 22:36

sobereva 发表于 2026-8-13 22:16
没你的输入文件没法进行测试
sobMECP的脚本是完全透明的，直接从源代码里搜大写的ERROR结合实际情况就能判 ...

老师您好，这是我的输入文件，Input_Tail没有内容

作者
Author: sobereva 时间: yesterday 01:02
我可以正常跑下去

作者
Author: sobereva 时间: yesterday 01:20
我又看了下，问题在于你用了scf=xqc，常规SCF迭代没收敛后自动用了QC，导致SCF Done:总共输出了两次，使得提取出的ab_initio内容格式错乱。甭用scf=xqc

作者
Author: 小学徒 时间: yesterday 10:04
本帖最后由小学徒于 2026-8-14 12:06 编辑

好的好的，非常感谢老师的帮助，想问一下老师除了将scf=xqc的关键词删除，老师在我上传的输入文件上还有什么其他的更改吗？

作者
Author: 小学徒 时间: yesterday 11:56
老师您好，我尝试了Input_Header_A和Input_Header_B的“#n uwb97xd/def2svp force guess(read)”命令行，运行runfirst.sh后显示正常结束，但是运行 ./runMECP.sh 后，程序报错，显示如下错误，并未出现4楼的界面。ReportFile中只有 “ERROR” 一词，没有其他具体错误信息。想问一下各位老师这可能是什么原因？如何解决？感谢各位老师指教！

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