计算化学公社

 找回密码 Forget password
 注册 Register
Views: 879|回复 Reply: 11
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 使用混合基组出现价电子数多于基函数的警示应该怎么解决

[复制链接 Copy URL]

70

帖子

0

威望

310

eV
积分
380

Level 3 能力者

如题,菜鸟第一次使用混合基组,对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!

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



831

帖子

1

威望

7190

eV
积分
8041

Level 6 (一方通行)

2#
发表于 Post on 2023-9-8 14:09:16 | 只看该作者 Only view this author
说明没有给 I 指定赝势。

5万

帖子

99

威望

5万

eV
积分
112393

管理员

公社社长

3#
发表于 Post on 2023-9-8 15:15:51 | 只看该作者 Only view this author
写输入文件没经验的话还不如直接都用def2-SVP,省得genecp的时候写错
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

70

帖子

0

威望

310

eV
积分
380

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2023-9-8 15:58:16 | 只看该作者 Only view this author
sobereva 发表于 2023-9-8 15:15
写输入文件没经验的话还不如直接都用def2-SVP,省得genecp的时候写错

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

5万

帖子

99

威望

5万

eV
积分
112393

管理员

公社社长

5#
发表于 Post on 2023-9-8 22:14:06 | 只看该作者 Only view this author
Priscilla 发表于 2023-9-8 15:58
请问老师我这样写错在哪,不是说def2系列并不需要额外指定赝势吗,我就没写......
C H 0
6-31G*

不需要定义赝势显然是你把基组名直接写在route section的情况
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

45

帖子

0

威望

135

eV
积分
180

Level 3 能力者

6#
发表于 Post on 2023-11-22 18:05:55 | 只看该作者 Only view this author
本帖最后由 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不收敛的方法加关键词呢

831

帖子

1

威望

7190

eV
积分
8041

Level 6 (一方通行)

7#
发表于 Post on 2023-11-22 19:39:05 | 只看该作者 Only view this author
wanwan5339 发表于 2023-11-22 18:05
你好,我也是出现他这种情况,我是Gaussian里没有def2svp基组,我是在网站上复制粘贴在分子坐标最下面的 ...

阅读 http://sobereva.com/60

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


这不是报错。

45

帖子

0

威望

135

eV
积分
180

Level 3 能力者

8#
发表于 Post on 2023-11-22 21:38:49 | 只看该作者 Only view this author
hebrewsnabla 发表于 2023-11-22 19:39
阅读 http://sobereva.com/60

我看了,我感觉自己也没写错,都是按格式来的。我是真的不知道如何解决
而且我用的就是第二个自定义基组的格式输入的,因为我def2svp这个基组不是自带的我是复制进去的,我如果加赝势的话,关键词就不能是gen了吧,是不是用gencp,然后分子序号之后复制基组,然后在复制赝势,请大佬解答,我是自学的,太多不理解的了,请多多包涵

831

帖子

1

威望

7190

eV
积分
8041

Level 6 (一方通行)

9#
发表于 Post on 2023-11-22 21:43:21 | 只看该作者 Only view this author
wanwan5339 发表于 2023-11-22 21:38
我看了,我感觉自己也没写错,都是按格式来的。我是真的不知道如何解决
而且我用的就是第二个自定义 ...

你说的没有问题。

45

帖子

0

威望

135

eV
积分
180

Level 3 能力者

10#
发表于 Post on 2023-11-22 21:49:57 | 只看该作者 Only view this author
本帖最后由 wanwan5339 于 2023-11-22 21:57 编辑

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

微信图片_20231122215654.png (12.14 KB, 下载次数 Times of downloads: 16)

微信图片_20231122215654.png

831

帖子

1

威望

7190

eV
积分
8041

Level 6 (一方通行)

11#
发表于 Post on 2023-11-22 21:57:06 | 只看该作者 Only view this author
wanwan5339 发表于 2023-11-22 21:49
大佬,感谢。我加赝势的时候是直接写46MDF还是把我下面这张图的内容复制进去呢(图上的内容是在Stuttgart ...

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

不需要。

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
wanwan5339 + 4 谢谢

查看全部评分 View all ratings

45

帖子

0

威望

135

eV
积分
180

Level 3 能力者

12#
发表于 Post on 2023-11-22 21:58:28 | 只看该作者 Only view this author
本帖最后由 wanwan5339 于 2023-11-22 22:00 编辑
hebrewsnabla 发表于 2023-11-22 21:57
你没有贴图。按 http://sobereva.com/60 第6节的方法写。

不需要。

好的,非常感谢

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 18:22 , Processed in 0.307039 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list