计算化学公社

 找回密码 Forget password
 注册 Register
Views: 181|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[辅助/分析程序] 求助:sobMECP运行runMECP.sh报错,Problem with the ab initio Job.

[复制链接 Copy URL]

13

帖子

0

威望

301

eV
积分
314

Level 3 能力者

各位老师好,最近在使用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?
感谢各位老师指教!



202608132153024084..png (11.97 KB, 下载次数 Times of downloads: 1)

202608132153024084..png

202608132201418426..png (84.63 KB, 下载次数 Times of downloads: 1)

势能面交叉示例

势能面交叉示例

202608132149318981..png (6.22 KB, 下载次数 Times of downloads: 1)

202608132149318981..png

Job0_B.log

90.23 KB, 下载次数 Times of downloads: 5

Job0_A.log

91.51 KB, 下载次数 Times of downloads: 2

ab_initio

5.23 KB, 下载次数 Times of downloads: 2

ReportFile

6 Bytes, 下载次数 Times of downloads: 0

6万

帖子

99

威望

6万

eV
积分
128952

管理员

公社社长

2#
发表于 Post on 前天 22:16 | 只看该作者 Only view this author
没你的输入文件没法进行测试
sobMECP的脚本是完全透明的,直接从源代码里搜大写的ERROR结合实际情况就能判断是因为什么触发的。跟初猜结构没关系
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

13

帖子

0

威望

301

eV
积分
314

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 前天 22:36 | 只看该作者 Only view this author
sobereva 发表于 2026-8-13 22:16
没你的输入文件没法进行测试
sobMECP的脚本是完全透明的,直接从源代码里搜大写的ERROR结合实际情况就能判 ...

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

geom

3.29 KB, 下载次数 Times of downloads: 2

Input_Header_A

100 Bytes, 下载次数 Times of downloads: 3

Input_Header_B

101 Bytes, 下载次数 Times of downloads: 1

6万

帖子

99

威望

6万

eV
积分
128952

管理员

公社社长

4#
发表于 Post on yesterday 01:02 | 只看该作者 Only view this author
我可以正常跑下去

202608140101593750..png (54.23 KB, 下载次数 Times of downloads: 51)

202608140101593750..png
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
128952

管理员

公社社长

5#
发表于 Post on yesterday 01:20 | 只看该作者 Only view this author
我又看了下,问题在于你用了scf=xqc,常规SCF迭代没收敛后自动用了QC,导致SCF Done:总共输出了两次,使得提取出的ab_initio内容格式错乱。甭用scf=xqc
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

13

帖子

0

威望

301

eV
积分
314

Level 3 能力者

6#
 楼主 Author| 发表于 Post on yesterday 10:04 | 只看该作者 Only view this author
本帖最后由 小学徒 于 2026-8-14 12:06 编辑

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

13

帖子

0

威望

301

eV
积分
314

Level 3 能力者

7#
 楼主 Author| 发表于 Post on yesterday 11:56 | 只看该作者 Only view this author
老师您好,我尝试了Input_Header_A和Input_Header_B的“#n uwb97xd/def2svp force guess(read)”命令行,运行runfirst.sh后显示正常结束,但是运行 ./runMECP.sh 后,程序报错,显示如下错误,并未出现4楼的界面。ReportFile中只有 “ERROR” 一词,没有其他具体错误信息。想问一下各位老师这可能是什么原因?如何解决?感谢各位老师指教!

202608141153204916..png (12.97 KB, 下载次数 Times of downloads: 1)

202608141153204916..png

ab_initio

5.3 KB, 下载次数 Times of downloads: 0

geom

3.35 KB, 下载次数 Times of downloads: 2

Input_Header_A

92 Bytes, 下载次数 Times of downloads: 0

Input_Header_B

93 Bytes, 下载次数 Times of downloads: 0

Job0_A.log

89.56 KB, 下载次数 Times of downloads: 0

Job0_B.log

89.78 KB, 下载次数 Times of downloads: 0

MECP.f

16.71 KB, 下载次数 Times of downloads: 2

ReportFile

6 Bytes, 下载次数 Times of downloads: 0

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-8-15 13:39 , Processed in 0.235165 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list