用Stuttgart基组选择不同核数赝势计算单点能，差异很大，这是为啥？

Tanya · 发表于 Post on 2018-9-2 20:39:08

本帖最后由 Tanya 于 2018-9-3 08:30 编辑

本人初学Gaussian中。

先用6-31g*/SDD混合基组，B3LYP方法，对我们实验室常用的化合物“三[双(三甲基硅基)胺基]镧”进行了结构优化，然后再用6-311+g**和5个不同的赝势基组计算结构优化后的单点能，发现差别挺大的，结果在下表中。

赝势基组	Energy （Hartree）	Dipole oment（Ddbye）
ECP28MWB-ANO	RB3LYP Energy -3056.36540455	2.0897
ECP28MWB-SEG	RB3LYP Energy -3056.36613588	2.0942
MWB28	RB3LYP Energy -3056.36613588	2.0942
MWB46	RB3LYP Energy -2652.27969067	2.1697
MWB47	UB3LYP Energy -2646.64291569	1.1488

注：

1、ECP28MWB-ANO和ECP28MWB-SEG是从网站上下载的赝势基组；

2、MWB28、MWB46和MWB47是Gaussian内置的Stuttgart基组关键词。

3、除La外，其他元素用6-311+g**

4、方法B3LYP，加上了GD3。

从计算结果可以看出，ECP28MWB-SEG和MWB28完全相同，说明内置的La小核y赝势基组MWB28与网站上的ECP28MWB-SEG完全一样，并且与同样是小核的ECP28MWB-ANO结果差别不大。但是MWB28，MWB46，MWB47之间的差别就很大了，我该相信哪个结果呢？

三[双(三甲基硅基)胺基]镧的结构：

file:///C:/Temp/msohtmlclip1/01/clip_image008.png

以下是计算结果截图：

ECP28MWB-ANO：

ECP28MWB-SEG：

MWB28：

MWB46：

MWB47：

19970306 · 发表于 Post on 2020-5-23 09:21:50

19970306 发表于 2020-5-22 21:46
您好，很抱歉打搅您，请问一下ECP28MWB-ANO和基组库里的ANO-RCC有什么区别吗

找到了

19970306 · 发表于 Post on 2020-5-22 21:46:52

您好，很抱歉打搅您，请问一下ECP28MWB-ANO和基组库里的ANO-RCC有什么区别吗

Tanya · 发表于 Post on 2018-9-4 09:59:26

sobereva 发表于 2018-9-4 03:00
在量化里，绝对能量是没有意义的，只有相对能量才是化学上感兴趣的，比较不同级别计算结果，也都应当比相 ...

嗯，这个是理解的。
这段时间一直在看你的博客和帖子，各种比较

。所以想看看自己新买的电脑计算不同的基组花多少时间，也算是个学习过程吧。

sobereva · 发表于 Post on 2018-9-4 03:00:50

Tanya 发表于 2018-9-3 22:36
做过很长时间分析测试工作，即使用不同方法测试，正确结果应该很接近，否则测试方法就有问题。不同方法结果 ...

在量化里，绝对能量是没有意义的，只有相对能量才是化学上感兴趣的，比较不同级别计算结果，也都应当比相对能量对应的化学量（诸如势垒、反应能量变化、激发能等）

Tanya · 发表于 Post on 2018-9-3 22:36:27

做过很长时间分析测试工作，即使用不同方法测试，正确结果应该很接近，否则测试方法就有问题。不同方法结果不一样，就有点心慌慌。看来这种思维方式不能带到量化中来。

sobereva · 发表于 Post on 2018-9-3 20:04:23

Tanya 发表于 2018-9-3 08:21
谢谢Sob老师的回答。
还有疑问：
同一个赝势基组，虽然内核电子数不同。但都有研究给出了专门的赝势， ...

不同基组间的能量不能横向对比这是量化计算最基本常识之一
赝势下计算的能量只是计算过程中直接考虑的电子的能量，显然大核和小核赝势会因为电子数不同而能量有巨大差异

Tanya · 发表于 Post on 2018-9-3 08:23:47

Accelerator 发表于 2018-9-3 04:14
想知道这些截图是什么可视化软件出来的

用了WebMO的免费版，可以从这里下载：https://www.webmo.net/#opennewwindow

单机使用挺方便的，可以直接用Web方式访问服务器，提交任务，可以在线看任务进行情况。
我的服务器放在学校办公室，通过学校的VPN方式，还可以直接在家里提交任务。

Tanya · 发表于 Post on 2018-9-3 08:21:20

sobereva 发表于 2018-9-2 23:50
计算中考虑的内核电子数不同，单点能自然有巨大差异，属于常识

谢谢Sob老师的回答。
还有疑问：
同一个赝势基组，虽然内核电子数不同。但都有研究给出了专门的赝势，不就是为了更准确么，差别这么大也算正常？
分子中一共有98个原子，只有La用了不同的基组，能量就差别就达到15%（MWB46相对MWB28），一般相对误差大约多少算是基本合理呢？
是不是不同的基组之间不能做横向比较？

Accelerator · 发表于 Post on 2018-9-3 04:14:15

想知道这些截图是什么可视化软件出来的

sobereva · 发表于 Post on 2018-9-2 23:50:33

计算中考虑的内核电子数不同，单点能自然有巨大差异，属于常识

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

[Gaussian/gview] 用Stuttgart基组选择不同核数赝势计算单点能，差异很大，这是为啥？

评分 Rate

浏览过的版块