在5月22-26日期间,本论坛仅限等级≥level 4的成员发新的主题贴,带来的不便请谅解!(此期间请勿给管理员Sobereva发私信问学术问题)

计算化学公社

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

[新手求助] 单点能计算时conver=6和conver=8的结果能否比较

[复制链接 Copy URL]

25

帖子

1

威望

556

eV
积分
601

Level 4 (黑子)

各位大佬好!小弟最近为了DIAS分析交了一部分单点能任务,有一部分出现了l502报错。按照sob老师博文(http://sobereva.com/61)里的方法:“G09/16对单点计算默认是SCF=tight,相当于SCF=conver=8,然而此时对密度矩阵的收敛要求有点太严了,往往不容易达到,而且达到的时候能量通常都早已经收敛到极高精度了。对于计算单点能的目的可放心降到conver=6,相当于把收敛标准放宽100倍,通常收敛时能量变化已经非常小了。”加了scf(novaracc,conver=6)后确实收敛了。于是小弟有了另一个顾虑:其他部分的单点能计算用默认设置自动收敛了,也就是 conver=8,其结果能量值能否与用conver=6的结果放在一起进行比较?


如需要统一标准,以下两个方法是否均可行呢?
1)对其他部分同一加上scf(novaracc,conver=6)指令进行计算
2)对默认设定下难以收敛的任务,以此运算结果的波函数为初猜,进一步进行conver=8的运算


第一次发帖,如果有不恰当的地方还请各位海涵,感谢指教!

6万

帖子

99

威望

6万

eV
积分
127029

管理员

公社社长

2#
发表于 Post on 2025-4-8 06:40:58 | 只看该作者 Only view this author
对于普通泛函算单点能,可以比较

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
czhang8 + 3 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心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

4334

帖子

4

威望

9868

eV
积分
14282

Level 6 (一方通行)

MOKIT开发者

3#
发表于 Post on 2025-4-8 07:07:35 | 只看该作者 Only view this author
正常收敛限的任务 去读取scf(conver=6)任务的轨道,收敛后再比较。如果是开壳层体系,注意检查波函数稳定性。

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
czhang8 + 3 谢谢

查看全部评分 View all ratings

自动做多参考态计算的程序MOKIT

25

帖子

1

威望

556

eV
积分
601

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2025-4-8 07:11:32 | 只看该作者 Only view this author
感谢各位!确实是开壳层的计算,用的b3lyp-d3 和def2-tzvp,我去学习一下怎么检查波函数稳定性。另外还想请问一下,加了XQC/QC才达到收敛的任务,和默认算法的任务,能量之间可以比较吗?

4334

帖子

4

威望

9868

eV
积分
14282

Level 6 (一方通行)

MOKIT开发者

5#
发表于 Post on 2025-4-8 12:56:06 | 只看该作者 Only view this author
czhang8 发表于 2025-4-8 07:11
感谢各位!确实是开壳层的计算,用的b3lyp-d3 和def2-tzvp,我去学习一下怎么检查波函数稳定性。另外还想请 ...

可以。
自动做多参考态计算的程序MOKIT

14

帖子

1

威望

311

eV
积分
345

Level 3 能力者

6#
发表于 Post on 2026-4-30 09:32:44 | 只看该作者 Only view this author
sobereva 发表于 2025-4-8 06:40
对于普通泛函算单点能,可以比较

老师我在根据mcpb.py对一个从酶促反应里挖出来的含过渡金属的开壳团簇体系(200+原子),含有较高负电(-5)进行一系列gaussian计算,小模型可以收敛(opt,fre计算),但用来计算mk-resp电荷的大模型试了很多方法scf都无法收敛,如果我的目的只是单纯为了得到这个体系里金属的参数跑md,可不可以在计算大模型的mk-resp电荷时使用conver=5?观察了scf log其实只要十来步就能达到这个能级但后面就开始震荡了

6万

帖子

99

威望

6万

eV
积分
127029

管理员

公社社长

7#
发表于 Post on 2026-4-30 11:18:05 | 只看该作者 Only view this author
ZiyuSong 发表于 2026-4-30 09:32
老师我在根据mcpb.py对一个从酶促反应里挖出来的含过渡金属的开壳团簇体系(200+原子),含有较高负电(- ...

至少=6
该怎么解决怎么解决
解决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

14

帖子

1

威望

311

eV
积分
345

Level 3 能力者

8#
发表于 Post on 2026-4-30 23:14:48 | 只看该作者 Only view this author
sobereva 发表于 2026-4-30 11:18
至少=6
该怎么解决怎么解决
解决SCF不收敛问题的方法

好的,谢谢sob老师

14

帖子

1

威望

311

eV
积分
345

Level 3 能力者

9#
发表于 Post on 2026-4-30 23:15:03 | 只看该作者 Only view this author
sobereva 发表于 2026-4-30 11:18
至少=6
该怎么解决怎么解决
解决SCF不收敛问题的方法

好的,谢谢sob老师

14

帖子

1

威望

311

eV
积分
345

Level 3 能力者

10#
发表于 Post on 2026-5-13 04:58:35 | 只看该作者 Only view this author
sobereva 发表于 2026-4-30 11:18
至少=6
该怎么解决怎么解决
解决SCF不收敛问题的方法

老师我发现这个体系优化时scf在第一步优化前的收敛过程中rmsdp可以收敛得比较顺利但maxdp在缓慢下降的趋势中出现了多次震荡(在e-3-e-4)而且在默认圈数耗尽前也没法收敛,我现在想的解决办法一是增大默认圈数,二是减少macconv提早进入xqc,不知道这样合理吗?或者先设置成conver=6进行优化再用默认值优化一遍?大概哪种解决方式更合理?

6万

帖子

99

威望

6万

eV
积分
127029

管理员

公社社长

11#
发表于 Post on 2026-5-13 16:22:40 | 只看该作者 Only view this author
ZiyuSong 发表于 2026-5-13 04:58
老师我发现这个体系优化时scf在第一步优化前的收敛过程中rmsdp可以收敛得比较顺利但maxdp在缓慢下降的趋 ...

增大默认圈数有用的前提是能看到有收敛趋势

没事甭用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

14

帖子

1

威望

311

eV
积分
345

Level 3 能力者

12#
发表于 Post on 2026-5-13 22:47:01 | 只看该作者 Only view this author
sobereva 发表于 2026-5-13 16:22
增大默认圈数有用的前提是能看到有收敛趋势

没事甭用QC,博文里明确说了缺点和风险,博文里那么多更值 ...

是这样的老师,我确实仔细看了那个博文里的所有方法,适合我这个体系的我也加上去了。现在的的泛函本身没有弥散(6311g**),所以在scf里就只加了fermi没有使用novaracc,noincfock,int=acc2e=12因为软件是g16).然后因为体系里有过渡金属ni且观察到gap较小,所以提升了vshift=500,隐式溶剂使用了ether是对应一篇文献里的参数。这个体系是用mcpb.py截出来的酶促反应中心含有个卟啉大环,我冻结了除alpha carbon之外的几个外周的o,我附上了我用的routine,希望老师能抽空看一眼对于我这个体系还有什么更好的解决办法,因为目前尝试的所有没有qc的方法确实没法收敛

3rd_opt.gjf

11.31 KB, 下载次数 Times of downloads: 3

6万

帖子

99

威望

6万

eV
积分
127029

管理员

公社社长

13#
发表于 Post on 2026-5-14 01:09:02 | 只看该作者 Only view this author
ZiyuSong 发表于 2026-5-13 22:47
是这样的老师,我确实仔细看了那个博文里的所有方法,适合我这个体系的我也加上去了。现在的的泛函本身没 ...

弥散函数是对于基组而言的不是对于泛函而言的

QC容易收敛到不稳定波函数,如果非要用这个,最后记得做波函数稳定性测试

别乱用Guess=Mix,根本不是这种场合用的

也可以尝试PBE0,由于gap更大因此SCF可能比TPSSh更容易收敛
北京科音自然科学研究中心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

14

帖子

1

威望

311

eV
积分
345

Level 3 能力者

14#
发表于 Post on 2026-5-15 04:00:03 | 只看该作者 Only view this author
sobereva 发表于 2026-5-14 01:09
弥散函数是对于基组而言的不是对于泛函而言的

QC容易收敛到不稳定波函数,如果非要用这个,最后记得做 ...

好的,谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2026-5-24 07:28 , Processed in 0.226906 second(s), 25 queries , Gzip On.

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