计算化学公社

标题: SCF（自洽场）收敛主要影响的是能量吗？ [打印本页]

作者
Author: diaolanxinyu 时间: 2015-4-21 22:26
标题: SCF（自洽场）收敛主要影响的是能量吗？
请问：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)关键词是不是只影响能量？

谢谢！

作者
Author: sobereva 时间: 2015-4-21 22:51
默认收敛标准是
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矩阵（进而影响频率）等等。

作者
Author: diaolanxinyu 时间: 2015-4-21 23:18
恩，知道了，谢谢sobereva老师！

作者
Author: nkallwar 时间: 2015-4-28 11:21

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收敛标准影响更大些？

作者
Author: sobereva 时间: 2015-4-28 16:19

nkallwar 发表于 2015-4-28 11:21
所以，在做单点能计算的时候，设置SCF=conver=6，大致还是很可信的？看输出文件里一个scf迭代完成 ...

是的。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)