计算化学公社

标题: 用ma-qzvp基组计算碘离子（I-）单点能时比不加弥散函数大10倍，哪里出问题了？ [打印本页]

作者
Author: ecustxu 时间: 2020-4-30 16:34
标题: 用ma-qzvp基组计算碘离子（I-）单点能时比不加弥散函数大10倍，哪里出问题了？
根据sob老师博文 http://sobereva.com/509中的方法，计算碘负离子的单点能，gaussian 软件，输入文件如下：# m062x genecp

Title Card Required

-1 1
I 0.00000000 0.00000000 0.00000000

@C:\ma-qZVP.txt

@C:\def2-ECP.txt

得到的碘负离子的能量是 -2443.21525680 a.u. 不加弥散函数，直接用qzvp基组得到的是-297.72395997 a.u ，显然后者更合理。但sob老师说阴离子体系要加弥散函数，不知道哪里出问题了，请老师指教。
还有我不加弥散函数算的I3-, I5-的UV曲线比较接近实验值，加了弥散函数明显就不对了。对于I的阴离子，推荐用什么基组计算能量和激发态？
谢谢大神！

作者
Author: 冰释之川 时间: 2020-4-30 16:37
将不同基组算出来的单点能进行比较是没意义的。基组越大，算出来的能量越负，这很正常。

作者
Author: snljty 时间: 2020-4-30 17:11
本帖最后由 snljty 于 2020-4-30 17:16 编辑

qzvp是什么，不记得def系列有叫QZVP的基组。你说的是def2-QZVP么。
你这两个数值差的太多，怎么看着像一个忘了写赝势另一个一个用了赝势了的感觉？
只看一个能量没有什么值更合理之说，你实际要算的也是和别的体系对比算相对能量。(比如HI分子和H+I-对比，只是举例。)

作者
Author: sun0615 时间: 2020-4-30 18:10
本帖最后由 sun0615 于 2020-4-30 18:19 编辑

是否是 spin orbital coupling 的问题? (如果是跟实验值比较有误差)

作者
Author: hebrewsnabla 时间: 2020-4-30 18:15
本帖最后由 hebrewsnabla 于 2020-4-30 18:19 编辑

snljty 发表于 2020-4-30 17:11
qzvp是什么，不记得def系列有叫QZVP的基组。你说的是def2-QZVP么。
你这两个数值差的太多，怎么看着像一个 ...

def系列还真有QZVP

但是sob的ma-QZVP是给def2-QZVP加ma得到的

作者
Author: snljty 时间: 2020-4-30 18:16

hebrewsnabla 发表于 2020-4-30 18:15
def系列还真有QZVP

谢谢，Gaussian里面有么，记得以前看说明书没有呀

作者
Author: hebrewsnabla 时间: 2020-4-30 18:57
本帖最后由 hebrewsnabla 于 2020-4-30 19:41 编辑

snljty 发表于 2020-4-30 18:16
谢谢，Gaussian里面有么，记得以前看说明书没有呀

G16手册里有是有，但是算了个例子发现给我算的是def2QZVP。（迷惑

orca里面是有的，不过只支持H-Kr。

看了文献和turbomole基组库发现QZVP和def2QZVP似乎是一回事。

作者
Author: liyuanhe211 时间: 2020-4-30 19:44
两个文件里都把使用的赝势和基组打印出来，然后上传，一定有至少其中一个不符合你“只加了弥散”的描述

作者
Author: liyuanhe211 时间: 2020-4-30 19:45

冰释之川发表于 2020-4-30 16:37
将不同基组算出来的单点能进行比较是没意义的。基组越大，算出来的能量越负，这很正常。

不会差这么多的

作者
Author: hebrewsnabla 时间: 2020-4-30 19:52
我按照博文中的方法试了一下，并无问题，加ma和不加的区别为几个mH。

你检查检查log，看看赝势有没有读进去。

作者
Author: ecustxu 时间: 2020-5-1 10:42

snljty 发表于 2020-4-30 17:11
qzvp是什么，不记得def系列有叫QZVP的基组。你说的是def2-QZVP么。
你这两个数值差的太多，怎么看着像一个 ...

谢谢您。确实发现计算中赝势没有读进去。

作者
Author: ecustxu 时间: 2020-5-1 10:44

sun0615 发表于 2020-4-30 18:10
是否是 spin orbital coupling 的问题? (如果是跟实验值比较有误差)

不太懂，我再好好查查，感谢

作者
Author: ecustxu 时间: 2020-5-1 10:52

liyuanhe211 发表于 2020-4-30 19:44
两个文件里都把使用的赝势和基组打印出来，然后上传，一定有至少其中一个不符合你“只加了弥散”的描述

谢谢您，发现是我的赝势没有读进去。我把两个@文件中间的空行删掉了，好像是读进去了。您看对吗？

作者
Author: ecustxu 时间: 2020-5-1 10:56

hebrewsnabla 发表于 2020-4-30 19:52
我按照博文中的方法试了一下，并无问题，加ma和不加的区别为几个mH。

你检查检查log，看看赝势有没有读 ...

感谢感谢，检查了下确实是赝势没有读进去。我把两个@文件中间的空行删了后算了下，发现可以了。不知道是我电脑版本问题还是sob老师在这个地方多加了空行。我这小菜鸟只知道照葫芦画瓢

作者
Author: sobereva 时间: 2020-5-5 08:57

ecustxu 发表于 2020-5-1 10:56
感谢感谢，检查了下确实是赝势没有读进去。我把两个@文件中间的空行删了后算了下，发现可以了。不知道是 ...

博文没有问题
不同系统下和换行符有关的差异而已

作者
Author: ecustxu 时间: 2020-5-6 00:05

sobereva 发表于 2020-5-5 08:57
博文没有问题
不同系统下和换行符有关的差异而已

明白了，谢谢sob老师

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