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

ecustxu · 发表于 Post on 2020-4-30 16:34:40

根据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的阴离子，推荐用什么基组计算能量和激发态？
谢谢大神！

冰释之川 · 发表于 Post on 2020-4-30 16:37:52

将不同基组算出来的单点能进行比较是没意义的。基组越大，算出来的能量越负，这很正常。

snljty · 发表于 Post on 2020-4-30 17:11:43

本帖最后由 snljty 于 2020-4-30 17:16 编辑

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

sun0615 · 发表于 Post on 2020-4-30 18:10:47

本帖最后由 sun0615 于 2020-4-30 18:19 编辑

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

hebrewsnabla · 发表于 Post on 2020-4-30 18:15:16

本帖最后由 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得到的

snljty · 发表于 Post on 2020-4-30 18:16:35

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

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

hebrewsnabla · 发表于 Post on 2020-4-30 18:57:43

本帖最后由 hebrewsnabla 于 2020-4-30 19:41 编辑

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

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

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

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

liyuanhe211 · 发表于 Post on 2020-4-30 19:44:53

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

liyuanhe211 · 发表于 Post on 2020-4-30 19:45:05

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

不会差这么多的

hebrewsnabla · 发表于 Post on 2020-4-30 19:52:45

我按照博文中的方法试了一下，并无问题，加ma和不加的区别为几个mH。

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

ecustxu · 发表于 Post on 2020-5-1 10:42:54

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

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

ecustxu · 发表于 Post on 2020-5-1 10:44:14

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

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

ecustxu · 发表于 Post on 2020-5-1 10:52:30

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

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

ecustxu · 发表于 Post on 2020-5-1 10:56:37

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

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

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

sobereva · 发表于 Post on 2020-5-5 08:57:41

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 用ma-qzvp基组计算碘离子（I-）单点能时比不加弥散函数大10倍，哪里出问题了？

浏览过的版块