计算化学公社

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

[Gaussian/gview] 加了gen关键词,但坐标后面无机组函数对结果有什么影响?

[复制链接 Copy URL]

12

帖子

0

威望

71

eV
积分
83

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师好,Sobereva老师帖子说,在输入文件中加入gen关键词代表坐标部分后面读入机组函数。这次我错误操作在坐标前面部分写入机组函数也加入了gen关键词,如:M062X/gen 6-311+G(2d,2p),坐标部分后面并无其他机组函数 。我查看log文件看到不一样的标准机组: 不加gen关键词是:Standard basis: 6-311+G(d,p) (5D, 7F),这应该是读取了坐标前面的机组函数;但是加了gen的关键词是:Standard basis: 6-31+(d',p') (6D, 7F),请问老师,这代表了什么呢?这样会对结果产生什么影响?

12

帖子

0

威望

71

eV
积分
83

Level 2 能力者

2#
 楼主 Author| 发表于 Post on 2020-9-23 19:01:49 | 只看该作者 Only view this author
我在这以下帖子找到了答案,谢谢大家。
1. http://bbs.keinsci.com/thread-19095-1-1.html
2. http://sobereva.com/51

115

帖子

1

威望

1122

eV
积分
1257

Level 4 (黑子)

3#
发表于 Post on 2020-9-23 20:47:09 | 只看该作者 Only view this author
本帖最后由 ldatea 于 2020-9-25 01:52 编辑
shengwenyang 发表于 2020-9-23 19:01
我在这以下帖子找到了答案,谢谢大家。
1. http://bbs.keinsci.com/thread-19095-1-1.html
2. http://sob ...

很显然你没有抓住问题的关键。这两帖子和你的情况有明显区别。
G09的手册的基组部分写了:
6-31G†:Gaussian 09 还包含 George Petersson 及合作者的 6-31G†和 6-31G‡基组,定
义为完备基组方法[Petersson88, Petersson91]的一部分。它们可以用关键词 6-31G(d')和
6-31G(d',p')调用,并可以添加单个或两个弥散函数;还可以添加 f 函数:例如,
6-31G(d'f),等

你那种奇怪的关键词导致的结果就是用了6-31+(2d',2p') (6D, 7F)
这和6-31+g(2d,2p)是有区别的。
  1. #p M062X/gen 6-311+G(2d,2p) gfinput

  2. 1

  3. 0 1
  4. H 0. 0. 0.
  5. H 0. 0. 0.742
复制代码
会得到
  1. Leave Link  202 at Wed Sep 23 20:34:29 2020, MaxMem=   104857600 cpu:               0.0 elap:               0.2
  2. (Enter /mnt/c/users/wyn/AppData/g16/l301.exe)
  3. Standard basis: 6-31+(2d',2p') (6D, 7F)
  4. Ernie: Thresh=  0.10000D-02 Tol=  0.10000D-05 Strict=F.
  5. AO basis set in the form of general basis input (Overlap normalization):
  6.       1 0
  7. S   3 1.00       0.000000000000
  8.       0.1873113696D+02  0.3349460434D-01
  9.       0.2825394365D+01  0.2347269535D+00
  10.       0.6401216923D+00  0.8137573261D+00
  11. S   1 1.00       0.000000000000
  12.       0.1612777588D+00  0.1000000000D+01
  13. P   1 1.00       0.000000000000
  14.       0.1500000000D+01  0.1000000000D+01
  15. P   1 1.00       0.000000000000
  16.       0.3750000000D+00  0.1000000000D+01
  17. ****
  18.       2 0
  19. S   3 1.00       0.000000000000
  20.       0.1873113696D+02  0.3349460434D-01
  21.       0.2825394365D+01  0.2347269535D+00
  22.       0.6401216923D+00  0.8137573261D+00
  23. S   1 1.00       0.000000000000
  24.       0.1612777588D+00  0.1000000000D+01
  25. P   1 1.00       0.000000000000
  26.       0.1500000000D+01  0.1000000000D+01
  27. P   1 1.00       0.000000000000
  28.       0.3750000000D+00  0.1000000000D+01
  29. ****
复制代码

  1. #p M062X 6-31+G(2d,2p) gfinput

  2. 1

  3. 0 1
  4. H 0. 0. 0.
  5. H 0. 0. 0.742
  6.      
复制代码
会得到
  1. Leave Link  202 at Wed Sep 23 20:41:42 2020, MaxMem=   104857600 cpu:               0.0 elap:               0.2
  2. (Enter /mnt/c/users/wyn/AppData/g16/l301.exe)
  3. Standard basis: 6-31+G(2d,2p) (6D, 7F)
  4. Ernie: Thresh=  0.10000D-02 Tol=  0.10000D-05 Strict=F.
  5. AO basis set in the form of general basis input (Overlap normalization):
  6.       1 0
  7. S   3 1.00       0.000000000000
  8.       0.1873113696D+02  0.3349460434D-01
  9.       0.2825394365D+01  0.2347269535D+00
  10.       0.6401216923D+00  0.8137573261D+00
  11. S   1 1.00       0.000000000000
  12.       0.1612777588D+00  0.1000000000D+01
  13. P   1 1.00       0.000000000000
  14.       0.2200000000D+01  0.1000000000D+01
  15. P   1 1.00       0.000000000000
  16.       0.5500000000D+00  0.1000000000D+01
  17. ****
  18.       2 0
  19. S   3 1.00       0.000000000000
  20.       0.1873113696D+02  0.3349460434D-01
  21.       0.2825394365D+01  0.2347269535D+00
  22.       0.6401216923D+00  0.8137573261D+00
  23. S   1 1.00       0.000000000000
  24.       0.1612777588D+00  0.1000000000D+01
  25. P   1 1.00       0.000000000000
  26.       0.2200000000D+01  0.1000000000D+01
  27. P   1 1.00       0.000000000000
  28.       0.5500000000D+00  0.1000000000D+01
  29. ****
复制代码

这个两个基组的定义就是不一样的,至少极化函数是不一样的。


评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
snljty + 3 正解

查看全部评分 View all ratings

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2020-9-24 10:45:33 | 只看该作者 Only view this author
那叫“基组”

gen 6-311+G(2d,2p)这种写法没有任何意义
实际用了什么基函数用gfinput关键词一看便知
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

12

帖子

0

威望

71

eV
积分
83

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2020-9-24 16:23:12 | 只看该作者 Only view this author
ldatea 发表于 2020-9-23 20:47
很显然你没有抓住问题的关键。这两帖子和你的情况有明显区别。
G09的手册的基组部分写了:
6-31G&#8224 ...

谢谢您,其实我是知道导致的基组是不一样的,很明显我也列出了log日志的标准基组。确实上面的帖子跟我的问题不完全一样,但是我得到了一些启发。谢谢。

12

帖子

0

威望

71

eV
积分
83

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2020-9-24 16:24:23 | 只看该作者 Only view this author
sobereva 发表于 2020-9-24 10:45
那叫“基组”

gen 6-311+G(2d,2p)这种写法没有任何意义

谢谢sobereva老师提醒,是我的疏忽,不严谨,确实是基组。

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 17:23 , Processed in 0.168495 second(s), 21 queries , Gzip On.

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