计算化学公社

标题: 使用混合基组出现价电子数多于基函数的警示应该怎么解决 [打印本页]

作者
Author: Priscilla 时间: 2023-9-8 12:43
标题: 使用混合基组出现价电子数多于基函数的警示应该怎么解决
如题，菜鸟第一次使用混合基组，对I使用了def2svp的基组，对其他C和H使用了6-31G*基组，监控log文件出现如下的输出：

Warning!  I  atom 1 has 53 valence electrons but only  26 basis functions. This is less than a minimal basis set!
Warning!  I  atom 2 has 53 valence electrons but only  26 basis functions.
This is less than a minimal basis set!
Warning!  I  atom 3 has 53 valence electrons but only  26 basis functions.
This is less than a minimal basis set!
Warning!  I  atom 4 has 53 valence electrons but only  26 basis functions.
This is less than a minimal basis set!

的警示，但是高斯依然在迭代计算，没有遇到过这种情况，请问大佬们这样的警示是什么原因，需要怎么做呢？

作者
Author: hebrewsnabla 时间: 2023-9-8 14:09
说明没有给 I 指定赝势。

作者
Author: sobereva 时间: 2023-9-8 15:15
写输入文件没经验的话还不如直接都用def2-SVP，省得genecp的时候写错

作者
Author: Priscilla 时间: 2023-9-8 15:58

sobereva 发表于 2023-9-8 15:15
写输入文件没经验的话还不如直接都用def2-SVP，省得genecp的时候写错

请问老师我这样写错在哪，不是说def2系列并不需要额外指定赝势吗，我就没写......
C H 0
6-31G*
****
N Cu I 0
Def2SVP
****

作者
Author: sobereva 时间: 2023-9-8 22:14

Priscilla 发表于 2023-9-8 15:58
请问老师我这样写错在哪，不是说def2系列并不需要额外指定赝势吗，我就没写......
C H 0
6-31G*

不需要定义赝势显然是你把基组名直接写在route section的情况

作者
Author: wanwan5339 时间: 2023-11-22 18:05
本帖最后由 wanwan5339 于 2023-11-22 18:15 编辑

hebrewsnabla 发表于 2023-9-8 14:09
说明没有给 I 指定赝势。

你好，我也是出现他这种情况，我是Gaussian里没有def2svp基组，我是在网站上复制粘贴在分子坐标最下面的，还需要指定I赝势吗，而且我的log文件最后出现Error on total polarization charges = 0.05443，这个应该是说我的结构不收敛吧？
请问这样的话我怎么解决呢，下面是我的关键词
%chk=diatrizoatea.chk
%nprocshared=12
%mem=2GB
# opt fre M062X/gen scrf
我的虚频消除打错了freq打成了fre，估计这个也是有影响的
我是不是既要加赝势，也要根据sob老师之前提到解决scf不收敛的方法加关键词呢

作者
Author: hebrewsnabla 时间: 2023-11-22 19:39

wanwan5339 发表于 2023-11-22 18:05
你好，我也是出现他这种情况，我是Gaussian里没有def2svp基组，我是在网站上复制粘贴在分子坐标最下面的 ...

阅读 http://sobereva.com/60

，而且我的log文件最后出现Error on total polarization charges = 0.05443，这个应该是说我的结构不收敛吧？

这不是报错。

作者
Author: wanwan5339 时间: 2023-11-22 21:38

hebrewsnabla 发表于 2023-11-22 19:39
阅读 http://sobereva.com/60

我看了，我感觉自己也没写错，都是按格式来的。我是真的不知道如何解决

而且我用的就是第二个自定义基组的格式输入的，因为我def2svp这个基组不是自带的我是复制进去的，我如果加赝势的话，关键词就不能是gen了吧，是不是用gencp，然后分子序号之后复制基组，然后在复制赝势，请大佬解答，我是自学的，太多不理解的了，请多多包涵

作者
Author: hebrewsnabla 时间: 2023-11-22 21:43

wanwan5339 发表于 2023-11-22 21:38
我看了，我感觉自己也没写错，都是按格式来的。我是真的不知道如何解决
而且我用的就是第二个自定义 ...

你说的没有问题。

作者
Author: wanwan5339 时间: 2023-11-22 21:49
本帖最后由 wanwan5339 于 2023-11-22 21:57 编辑

hebrewsnabla 发表于 2023-11-22 21:43
你说的没有问题。

大佬，感谢。我加赝势的时候是直接写46MDF还是把我下面这张图的内容复制进去呢（图上的内容是在Stuttgart系列赝势的开发者的网站上找到的）？
我还想请问您，Error on total polarization charges = 0.05443，我也阅读了一些帖子，帖子上大部分回复都是说阅读sob老师的一篇博文，我是不是也要根据博文上的解决办法加一些关键词呢？

作者
Author: hebrewsnabla 时间: 2023-11-22 21:57

wanwan5339 发表于 2023-11-22 21:49
大佬，感谢。我加赝势的时候是直接写46MDF还是把我下面这张图的内容复制进去呢（图上的内容是在Stuttgart ...

你没有贴图。按 http://sobereva.com/60 第6节的方法写。

不需要。

作者
Author: wanwan5339 时间: 2023-11-22 21:58
本帖最后由 wanwan5339 于 2023-11-22 22:00 编辑

hebrewsnabla 发表于 2023-11-22 21:57
你没有贴图。按 http://sobereva.com/60 第6节的方法写。

不需要。

好的，非常感谢

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