计算化学公社

标题: 加了gen关键词，但坐标后面无机组函数对结果有什么影响？ [打印本页]
作者
Author: shengwenyang    时间: 2020-9-23 16:59
标题: 加了gen关键词，但坐标后面无机组函数对结果有什么影响？
各位老师好，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)，请问老师，这代表了什么呢？这样会对结果产生什么影响？
作者
Author: shengwenyang    时间: 2020-9-23 19:01
我在这以下帖子找到了答案，谢谢大家。
1. http://bbs.keinsci.com/thread-19095-1-1.html
2. http://sobereva.com/51
作者
Author: ldatea    时间: 2020-9-23 20:47
本帖最后由 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

  3. 1

  5. 0 1
  6. H 0. 0. 0.
  7. 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

  3. 1

  5. 0 1
  6. H 0. 0. 0.
  7. H 0. 0. 0.742
  8.      
复制代码
会得到
  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. ****
复制代码

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


作者
Author: sobereva    时间: 2020-9-24 10:45
那叫“基组”

gen 6-311+G(2d,2p)这种写法没有任何意义
实际用了什么基函数用gfinput关键词一看便知
作者
Author: shengwenyang    时间: 2020-9-24 16:23
ldatea 发表于 2020-9-23 20:47
很显然你没有抓住问题的关键。这两帖子和你的情况有明显区别。
G09的手册的基组部分写了：
6-31G&#8224 ...

谢谢您，其实我是知道导致的基组是不一样的，很明显我也列出了log日志的标准基组。确实上面的帖子跟我的问题不完全一样，但是我得到了一些启发。谢谢。
作者
Author: shengwenyang    时间: 2020-9-24 16:24
sobereva 发表于 2020-9-24 10:45
那叫“基组”

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

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



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3