计算化学公社

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

[Gaussian/gview] 求助:Search did not lower the energy significantly. 怎么解决?

[复制链接 Copy URL]

235

帖子

0

威望

3374

eV
积分
3609

Level 5 (御坂)

输入文件如下(因为不收敛,所以使用了qc)。谢谢!!!

%chk=C18d9h.chk
%mem=50GB
%nproc=24
#opt=(maxcyc=200) wB97XD/def2TZVP scf=(maxcyc=500,qc)

C18d9h

0 1   
C             1.15341512737721    6.88190633624316    0.00000000000000
C            -1.15341512737721    6.88190633624316    0.00000000000000
C            -3.54003528830724    6.01324647973673    0.00000000000000
C            -5.30717110012025    4.53044452986863    0.00000000000000
C            -6.57706516638981    2.33092237312197    0.00000000000000
C            -6.97764175146146    0.05913708935555    0.00000000000000
C            -6.53661368895820   -2.44206472543734    0.00000000000000
C            -5.38319667185486   -4.43983972107969    0.00000000000000
C            -3.43760646315421   -6.07238545881840    0.00000000000000
C            -1.26989406626957   -6.86136690299060    0.00000000000000
C             1.26989406626957   -6.86136690299060    0.00000000000000
C             3.43760646315421   -6.07238545881840    0.00000000000000
C             5.38319667185486   -4.43983972107969    0.00000000000000
C             6.53661368895820   -2.44206472543734    0.00000000000000
C             6.97764175146146    0.05913708935555    0.00000000000000
C             6.57706516638981    2.33092237312197    0.00000000000000
C             5.30717110012025    4.53044452986863    0.00000000000000
C             3.54003528830724    6.01324647973673    0.00000000000000
Be            5.72500790882509    6.66250025708779   -1.48556119378526




输出文件中出现下面的这样的错误
Search did not lower the energy significantly.
No lower point found -- switch to steepest descent.
Search did not lower the energy significantly.
Restarting incremental Fock formation.
Search did not lower the energy significantly.
Restarting incremental Fock formation.
Search did not lower the energy significantly.
No lower point found -- run aborted.
Error termination via Lnk1e in /Share/home/05141/g09/l508.exe at Tue Sep  8 07:39:55 2020.
Job cpu time:       2 days 20 hours 29 minutes  1.7 seconds.


6万

帖子

99

威望

6万

eV
积分
125128

管理员

公社社长

2#
发表于 Post on 2020-9-8 08:52:28 | 只看该作者 Only view this author
用qc完全是多余的。18碳环体系用wB97XD收敛实在太容易了,除非你的初始结构太烂
不要以为用qc就能解决问题

较大体系记得用#P,省得难收敛还蒙在鼓里
北京科音自然科学研究中心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

235

帖子

0

威望

3374

eV
积分
3609

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2020-9-8 09:11:33 | 只看该作者 Only view this author
sobereva 发表于 2020-9-8 08:52
用qc完全是多余的。18碳环体系用wB97XD收敛实在太容易了,除非你的初始结构太烂
不要以为用qc就能解决问题 ...

谢谢社长。最初就只有  #opt wB97XD/def2TZVP  这些关键词,但几种结构都不收敛,所以才使用了qc

334

帖子

1

威望

6153

eV
积分
6507

Level 6 (一方通行)

4#
发表于 Post on 2020-9-8 09:35:27 | 只看该作者 Only view this author
顺便看了眼结构,记得之前看社长的帖子,C18环键长是1.3和1.2 Angstrom 左右的交替的么, 而你这个输入文件里的结构是2.5和2.3 A 键长,差这么多是有意而为之的么。键长差了这么远,不好收敛很正常。还是说这个坐标其实原本的单位是bohr呀。

235

帖子

0

威望

3374

eV
积分
3609

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2020-9-8 09:43:06 | 只看该作者 Only view this author
paramecium86 发表于 2020-9-8 09:35
顺便看了眼结构,记得之前看社长的帖子,C18环键长是1.3和1.2 Angstrom 左右的交替的么, 而你这个输入文件 ...

这个还真没注意,先是用社长文章里的结构优化的,然后加了其他原子。多谢提醒!!

136

帖子

0

威望

2207

eV
积分
2343

Level 5 (御坂)

6#
发表于 Post on 2023-3-14 20:51:17 | 只看该作者 Only view this author
sobereva 发表于 2020-9-8 08:52
用qc完全是多余的。18碳环体系用wB97XD收敛实在太容易了,除非你的初始结构太烂
不要以为用qc就能解决问题 ...

老师,请问Search did not lower the energy significantly. l508 这是什么原因导致的报错啊?我没有加QC,做的柔性扫描键长,同样报这个错误。

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

7#
发表于 Post on 2023-3-14 21:04:15 | 只看该作者 Only view this author
嘤嘤嘤 发表于 2023-3-14 13:51
老师,请问Search did not lower the energy significantly. l508 这是什么原因导致的报错啊?我没有加QC ...

有没有加其他类似关键词,例如xqc?
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

136

帖子

0

威望

2207

eV
积分
2343

Level 5 (御坂)

8#
发表于 Post on 2023-3-14 21:09:18 | 只看该作者 Only view this author
wzkchem5 发表于 2023-3-14 21:04
有没有加其他类似关键词,例如xqc?

老师,加了。以下是我的输入文件:# um062x/6-31g* guess(mix,always) opt=modredundant scf=(xqc,maxconventionalcycles=64)。柔性扫描40步,0.05埃

6万

帖子

99

威望

6万

eV
积分
125128

管理员

公社社长

9#
发表于 Post on 2023-3-14 21:31:34 | 只看该作者 Only view this author
嘤嘤嘤 发表于 2023-3-14 21:09
老师,加了。以下是我的输入文件:# um062x/6-31g* guess(mix,always) opt=modredundant scf=(xqc,maxcon ...

常规SCF迭代不收敛时自动切换到qc,等于还是用了qc
北京科音自然科学研究中心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

136

帖子

0

威望

2207

eV
积分
2343

Level 5 (御坂)

10#
发表于 Post on 2023-3-14 22:08:17 | 只看该作者 Only view this author
sobereva 发表于 2023-3-14 21:31
常规SCF迭代不收敛时自动切换到qc,等于还是用了qc

老师,可是不加之前他是502不收敛的。之前加了这一对关键词还是不收敛。# um062x/6-31g* guess(mix,always) opt=modredundant pop=always nosymm SCF=NoIncFock int=superfine。不知道为啥,就很简单的柔性扫描键长一直不收敛。

6万

帖子

99

威望

6万

eV
积分
125128

管理员

公社社长

11#
发表于 Post on 2023-3-14 22:54:03 | 只看该作者 Only view this author
嘤嘤嘤 发表于 2023-3-14 22:08
老师,可是不加之前他是502不收敛的。之前加了这一对关键词还是不收敛。# um062x/6-31g* guess(mix,alway ...

解决SCF不收敛问题的方法
http://sobereva.com/61
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 04:29 , Processed in 0.228963 second(s), 20 queries , Gzip On.

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