计算化学公社

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

[Gaussian/gview] 有关stable关键词的疑问

[复制链接 Copy URL]

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 2238663 于 2020-9-24 11:02 编辑

各位老师好!我在对一个配合物进行电荷分布计算时,采用pop关键词,但由于SCF难以收敛,参考了Sob老师的《解决SCF不收敛问题的方法》一文。
我在单点分析任务中采用了SCF=conver=6,文中提到这样放宽收敛对单点能的计算基本不产生影响,可以大胆使用,但是我用了stable测出我的波函数不稳定……
对于电荷分布的计算,使用该关键词结果是否还可信呢?

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

2#
发表于 Post on 2020-9-24 10:30:24 | 只看该作者 Only view this author
所以和stable关键词有什么关系

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

3#
发表于 Post on 2020-9-24 10:34:44 | 只看该作者 Only view this author
当前问题和stable关键词没有丝毫联系
conver=6的波函数一般来说可用
北京科音自然科学研究中心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

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2020-9-24 11:02:22 | 只看该作者 Only view this author
sobereva 发表于 2020-9-24 10:34
当前问题和stable关键词没有丝毫联系
conver=6的波函数一般来说可用

我加了stable关键词 说我的波函数不稳定 我不知道电荷分布结果是否可信

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

5#
发表于 Post on 2020-9-24 11:03:47 | 只看该作者 Only view this author
2238663 发表于 2020-9-24 11:02
我加了stable关键词 说我的波函数不稳定 我不知道电荷分布结果是否可信

不稳定就stable=opt变稳定咯

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2020-9-24 11:47:13 | 只看该作者 Only view this author
喵星大佬 发表于 2020-9-24 11:03
不稳定就stable=opt变稳定咯

恩…体系比较大 是直接用坐标计算单点任务 不涉及结构优化

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

7#
发表于 Post on 2020-9-24 12:05:43 | 只看该作者 Only view this author
2238663 发表于 2020-9-24 11:47
恩…体系比较大 是直接用坐标计算单点任务 不涉及结构优化

本来stable=opt和结构优化也没有任何联系
北京科音自然科学研究中心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

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

8#
发表于 Post on 2020-9-24 12:28:39 | 只看该作者 Only view this author
2238663 发表于 2020-9-24 11:47
恩…体系比较大 是直接用坐标计算单点任务 不涉及结构优化

stable=opt的opt指的是优化波函数使其变得稳定,和单独写的opt关键字(结构优化)没有任何关系
而且并不是说不设scf=conver=6就一定能稳定,不要把不稳定怪罪到把conver放松到6上,不稳定性主要是因为你的初猜和最后的稳定波函数有定性差别,比如空间对称性不同,或该破缺的对称性(空间或自旋对称性)没有破缺
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
本团队长期招收研究生,有意者可私信联系

4289

帖子

4

威望

9543

eV
积分
13912

Level 6 (一方通行)

MOKIT开发者

9#
发表于 Post on 2020-9-24 12:47:35 | 只看该作者 Only view this author
2238663 发表于 2020-9-24 11:02
我加了stable关键词 说我的波函数不稳定 我不知道电荷分布结果是否可信

波函数不稳定,要看是哪种类型的不稳定,如果是内部不稳定(internal instability),布局分析不可信。

如果你用了RDFT,且显示的是RHF -> UHF不稳定,那就不好讲了。此时若你稳定的UDFT波函数能量只比RDFT能量低一点点,且<S**2>非常接近于理想值,那么两个波函数的布局分析结果是很接近的,是可信的。若你稳定的UDFT波函数能量比RDFT能量低很多,那RDFT结果不怎么可信。
自动做多参考态计算的程序MOKIT

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2020-9-24 13:17:36 | 只看该作者 Only view this author
先谢谢各位老师的指导了。对波函数稳定性的测试是internal instability,请教一下:我下一步使用stable=opt后,波函数稳定后,布局分析产生的电荷值可信吗?

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

11#
发表于 Post on 2020-9-24 14:25:26 | 只看该作者 Only view this author
2238663 发表于 2020-9-24 13:17
先谢谢各位老师的指导了。对波函数稳定性的测试是internal instability,请教一下:我下一步使用stable=opt ...

只能说是不会因为波函数稳定性的问题导致不可信,因为决定电荷是否可信的还有其他因素,比如你的体系是否有强相关,是否容易受到自相互作用误差的影响,等等。
强相关体系有instability的概率比较大,所以反过来讲,出现instabililty的时候建议检查一下体系是否有强相关特征,如果有明显的强相关,就需要用多参考态方法,DFT算出来的无论如何是不可信的。
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
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 01:03 , Processed in 0.182146 second(s), 20 queries , Gzip On.

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