计算化学公社

标题: MRCC的DFT计算与Gaussian不符 [打印本页]

作者
Author: liyaohx 时间: 2020-1-2 15:14
标题: MRCC的DFT计算与Gaussian不符
本帖最后由 liyaohx 于 2020-1-2 15:19 编辑

这是我MRCC的输入文件：

basis=aug-cc-pVDZ
calc=pw91
scftype=uhf
mem=50000MB
dfbasiscor=none
dfbasis_scf=none
mult=8
charge=-1
verbosity=3
scfiguess=ao
unit=angs
geom=xyz
5

V                -0.37577200 -0.61692500 0.00001200
V                -2.06758100 0.02363600 -0.00001000
C                -0.58122200 1.38919100 0.00000200
C                0.72022300 1.21999200 0.00002000
Cr                2.30679700 -0.08372700 -0.00000800

计算结果：
MRCC：-3010.478169
Gaussian：-3008.835230
感觉结果很诡异，不知道为什么结果偏差这么大？是我输入文件写的不对吗？跪求大神能指点一下。

作者
Author: Warm_Cloud 时间: 2020-1-2 17:50
开壳层体系波函数很容易跑到不同的极小值

作者
Author: sobereva 时间: 2020-1-2 19:05
鉴于MRCC给出的能量更低，应当对Gaussian的情况用stable做一下波函数稳定性测试。也可以考虑把MRCC的molden文件载入Multiwfn，然后导出成fch，再用unfchk转化成chk，从而在Gaussian计算时用guess=read从中读取初猜波函数，看看能否收敛到比之前能量更低的态

但在此之前，建议先让UHF的能量对上，然后再说DFT的事。

作者
Author: yflchx 时间: 2020-1-2 20:38
你用DFT，积分格点一致吗？MRCC的DFT默认用adaptive integration grid，而Gaussian不是。

1. 确保积分格点一致。

2. 确保最后收敛的电子态一致。MRCC里用UHF有时需要指定occ。（很多程序的SCF不像Gaussian那样优越，需要研究者更多的参与其中。）

3.你看看MRCC的论坛，integ.f模块作者刚重新写过，adaptive integration grid对某些体系存在问题。

作者
Author: liyaohx 时间: 2020-1-2 20:58
感谢各位老师的回复，我这边再按老师们的建议尝试一下

作者
Author: liyaohx 时间: 2020-1-3 09:36
本帖最后由 liyaohx 于 2020-1-3 10:03 编辑

我对Gaussian使用stable=opt，检验结果显示波函数是稳定的，用Mutiwfn转换mrcc产生的molden到fchk的时候提示说缺少基函数信息，无法转换。这个体系的UHF我试了，得到的能量跟Gaussian在一个量级，不会出现偏差这么大的情况。我又试了下N2和Mn的原子簇，发现DFT都能和Gaussian对上，就只有Cr的这个体系不行。我这个体系没有对称性，我看了程序输出文件显示的占据数信息是对的，α/β=45/38

作者
Author: liyaohx 时间: 2020-1-3 11:13
本帖最后由 liyaohx 于 2020-1-3 11:20 编辑

找到问题了，就是adaptive integration grid的锅，设置agrid=LD0110-LD0974就正常了

，请问各位老师知道怎么让mrcc显示T1 诊断值和调节CC的shift值吗，之前一直用的molpro，molpro里面这两个都可以

作者
Author: yflchx 时间: 2020-1-3 12:03

liyaohx 发表于 2020-1-3 11:13
找到问题了，就是adaptive integration grid的锅，设置agrid=LD0110-LD0974就正常了，请问各位老师知 ...

我一直用LNO-CCSD(T)，T1和D1自动输出，你看看是不是都这样？
CCSD iteration has converged in 17 steps.
5229.01 MB memory was utilized by CCSD

T1 diagnostic: 1.497E-02 D1 diagnostic: 7.118E-02

至于CC的shift调节，这个我几乎用不到，也没关心过。好好看看手册，或上MRCC论坛问问。

作者
Author: chrinide 时间: 2020-1-3 12:42
MRCC的DFT默认的求积格点有点问题，对过渡金属比较多的体系问题比较大，这个体系对DFT格点比较敏感，你做以下设置，就可以得到正确的结果

agrid=LD0110-LD0590
grtol=16

***FINAL KOHN-SHAM ENERGY: -3008.8348612999907346 [AU]

就和G对上了

作者
Author: zjxitcc 时间: 2020-1-3 12:46

liyaohx 发表于 2020-1-3 11:13
找到问题了，就是adaptive integration grid的锅，设置agrid=LD0110-LD0974就正常了，请问各位老师知 ...

T1诊断值高斯也能给出，做个CCSD计算就好了。用MOLPRO和MRCC有点杀鸡用牛刀···

作者
Author: sobereva 时间: 2020-1-3 18:15

liyaohx 发表于 2020-1-3 09:36
我对Gaussian使用stable=opt，检验结果显示波函数是稳定的，用Mutiwfn转换mrcc产生的molden到fchk的时候提 ...

“提示说缺少基函数信息”是因为后缀不是.molden，导致没有当做molden输入文件来读入

作者
Author: liyaohx 时间: 2020-1-7 16:22
感谢各位老师的回复

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