计算化学公社

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

[Gaussian/gview] SCF(自洽场)收敛主要影响的是能量吗?

[复制链接 Copy URL]

117

帖子

0

威望

1126

eV
积分
1243

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
请问:SCF(自洽场)收敛主要影响的是能量吗?
当在Gaussian计算中们加入SCF(CONVER=6)时,会有如下自洽场的收敛标准信息(密度矩阵)

Requested convergence on RMS density matrix=1.00D-06 within 500 cycles.
Requested convergence on MAX density matrix=1.00D-04.
Requested convergence on             energy=1.00D-04.
No special actions if energy rises.
SCF Done:  E(UHF) =  -114.960376997     A.U. after    6 cycles
            NFock=  6  Conv=0.54D-06     -V/T= 2.0025




问题一:SCF(CONVER=6)关键词是不是影响能量E的值,是不是可以理解成CONVER=6,能量精确到小数点后六位,等于8能量精确到小数点后八位?


问题二:SCF(CONVER=6)关键词是不是只影响能量?

谢谢!

5万

帖子

99

威望

5万

eV
积分
112353

管理员

公社社长

2#
发表于 Post on 2015-4-21 22:51:38 | 只看该作者 Only view this author
默认收敛标准是
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.

用scf(conver=N)  可指定RMS density matrix收敛限是1E-N,MAX density matrix和energy收敛限都是1E(-N+2)。例如下面是N=11
Requested convergence on RMS density matrix=1.00D-11 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-09.
Requested convergence on             energy=1.00D-09.

所以显然=6不是精确到后N位。默认是N=8,这个标准对波函数收敛精度要求较高,此时达到收敛限时,通常能量已经收敛得很精确了,因此如果只关心能量,可以设为=6让收敛标准适当降低。

不只影响能量,还影响波函数、受力、Hessian矩阵(进而影响频率)等等。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

117

帖子

0

威望

1126

eV
积分
1243

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2015-4-21 23:18:29 | 只看该作者 Only view this author
恩,知道了,谢谢sobereva老师!

67

帖子

0

威望

1039

eV
积分
1106

Level 4 (黑子)

4#
发表于 Post on 2015-4-28 11:21:44 | 只看该作者 Only view this author
sobereva 发表于 2015-4-21 22:51
默认收敛标准是
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested c ...

      所以,在做单点能计算的时候,设置SCF=conver=6,大致还是很可信的?看输出文件里一个scf迭代完成的过程里,后面的能量变化已经很小了...
     不过在结构优化里,由于scf收敛标准对受力、频率影响比较大,所以scf=conver=6得到的结果可信度就差些?因此如果用scf=conver=6做完结构优化最好再去掉这个关键词再优化一下比较好么?
    这么看起来,似乎就是结构优化对基组、方法的依赖性小一些,通常在较低的水平下就能得到基本可靠的结构,但反而受scf收敛标准影响更大些?

5万

帖子

99

威望

5万

eV
积分
112353

管理员

公社社长

5#
发表于 Post on 2015-4-28 16:19:19 | 只看该作者 Only view this author
nkallwar 发表于 2015-4-28 11:21
所以,在做单点能计算的时候,设置SCF=conver=6,大致还是很可信的?看输出文件里一个scf迭代完成 ...

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

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 18:44 , Processed in 0.191942 second(s), 21 queries , Gzip On.

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