应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】下面基函数的大小排序谁知道
191/2返回列表12下一页
查看: 3384  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

loovfnd

至尊木虫 (著名写手)

[交流] 【求助】下面基函数的大小排序谁知道

1. aug-cc-pVDZ，
2. aug-cc-pVTZ，
3. aug-cc-pVQZ，
4. aug-cc-pV5Z，
5. 6-311++g（3df）
回复此楼

» 猜你喜欢
莫强求！
1楼 2009-08-24 09:05:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

我来做个总结

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
loovfnd(金币+2,VIP+0):多谢交流！ 8-24 14:00
对于第二周期元素C，几个常用基组的基组规模测试结果：
                        basis functions          primitive gaussians
sto-3g                       5                          15
3-21g                        9                          25
6-31g*                      15                          28
6-31+g**                    19                          32
6-311++g**                  22                          36
6-311++g(3df)               39                          58
6-311++g(3d2f,2pd)          46                          68
sddall                       8                          16
tzvp                        19                          35
aug-cc-pvdz                 23                          43
aug-cc-pvtz                 46                          73
aug-cc-pvqz                 80                         127
aug-cc-pv5z                127                         209
aug-cc-pv6z                189                         326
===gjf测试文件开始===
%kjob L302
# hf/sto-3g

a

0 5
C

--Link1--
%kjob L302
# hf/3-21g

a

0 5
C

--Link1--
%kjob L302
# hf/6-31g*

a

0 5
C

--Link1--
%kjob L302
# hf/6-31+g(d,p)

a

0 5
C

--Link1--
%kjob L302
# hf/6-311++g**

a

0 5
C

--Link1--
%kjob L302
# hf/6-311++g(3df)

a

0 5
C

--Link1--
%kjob L302
# hf/6-311++g(3d2f,2pd)

a

0 5
C

--Link1--
%kjob L302
# hf/sddall

a

0 5
C

--Link1--
%kjob L302
# hf/tzvp

a

0 5
C

--Link1--
%kjob L302
# hf/aug-cc-pvdz

a

0 5
C

--Link1--
%kjob L302
# hf/aug-cc-pvtz

a

0 5
C

--Link1--
%kjob L302
# hf/aug-cc-pvqz

a

0 5
C

--Link1--
%kjob L302
# hf/aug-cc-pv5z

a

0 5
C

--Link1--
%kjob L302
# hf/aug-cc-pv6z

a

0 5
C

===gjf文件结束===
赞一下(12人) 回复此楼
10楼2009-08-24 13:50:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
loovfnd(金币+3):感谢帮助
你可以加入一个命令%kjob L??? 1
在合适的Link下，添加一个kjob命令。
来测试你的机组规模

我曾经测试过一个计算任务的基组规模，原始输入文件如下：
===gjf文件开始===
%kjob L302 1
%mem=1gb
# pbe1pbe/genecp gfprint

PcCePc-D4d, sddall for Ce, 6-311+g(3d2f) for the four inner N, and 6-31g(d) for other atoms.

0 1
H                 -1.68395303   -7.28940820    2.78817738
C                 -1.82074393   -6.23510460    2.56314038
C                 -0.70999313   -5.44022282    2.28518040
C                 -3.12142436   -5.69634512    2.56314038
C                 -0.93505629   -4.09455913    1.98912040
H                  0.29735008   -5.84415908    2.29309142
C                 -3.34477749   -4.34885941    2.28518040
H                 -3.96365537   -6.34512458    2.78817738
C                 -2.23410588   -3.55647518    1.98912040
C                 -0.03161546   -2.98966237    1.68726442
H                 -4.34270278   -3.92218626    2.29309142
C                 -2.09165504   -2.13636604    1.68726442
N                  1.28923817   -3.11249628    1.74366744
N                 -0.76494518   -1.84674103    1.45160441
N                 -3.11249628   -1.28923817    1.74366744
C                  2.13636604   -2.09165504    1.68726442
C                 -2.98966237    0.03161546    1.68726442
N                  1.84674103   -0.76494518    1.45160441
C                  3.55647518   -2.23410588    1.98912040
N                 -1.84674103    0.76494518    1.45160441
C                 -4.09455913    0.93505629    1.98912040
C                  2.98966237   -0.03161546    1.68726442
C                  4.34885941   -3.34477749    2.28518040
C                  4.09455913   -0.93505629    1.98912040
C                 -2.13636604    2.09165504    1.68726442
C                 -5.44022282    0.70999313    2.28518040
C                 -3.55647518    2.23410588    1.98912040
N                  3.11249628    1.28923817    1.74366744
C                  5.69634512   -3.12142436    2.56314038
H                  3.92218626   -4.34270278    2.29309142
C                  5.44022282   -0.70999313    2.28518040
N                 -1.28923817    3.11249628    1.74366744
C                 -6.23510460    1.82074393    2.56314038
H                 -5.84415908   -0.29735008    2.29309142
C                 -4.34885941    3.34477749    2.28518040
C                  2.09165504    2.13636604    1.68726442
C                  6.23510460   -1.82074393    2.56314038
H                  6.34512458   -3.96365537    2.78817738
H                  5.84415908    0.29735008    2.29309142
C                  0.03161546    2.98966237    1.68726442
C                 -5.69634512    3.12142436    2.56314038
H                 -7.28940820    1.68395303    2.78817738
H                 -3.92218626    4.34270278    2.29309142
C                  2.23410588    3.55647518    1.98912040
N                  0.76494518    1.84674103    1.45160441
H                  7.28940820   -1.68395303    2.78817738
C                  0.93505629    4.09455913    1.98912040
H                 -6.34512458    3.96365537    2.78817738
C                  3.34477749    4.34885941    2.28518040
C                  0.70999313    5.44022282    2.28518040
C                  3.12142436    5.69634512    2.56314038
H                  4.34270278    3.92218626    2.29309142
C                  1.82074393    6.23510460    2.56314038
H                 -0.29735008    5.84415908    2.29309142
H                  3.96365537    6.34512458    2.78817738
H                  1.68395303    7.28940820    2.78817738
H                  3.96365537   -6.34512458   -2.78817738
C                  3.12142436   -5.69634512   -2.56314038
C                  3.34477749   -4.34885941   -2.28518040
C                  1.82074393   -6.23510460   -2.56314038
C                  2.23410588   -3.55647518   -1.98912040
H                  4.34270278   -3.92218626   -2.29309142
C                  0.70999313   -5.44022282   -2.28518040
H                  1.68395303   -7.28940820   -2.78817738
C                  0.93505629   -4.09455913   -1.98912040
C                  2.09165504   -2.13636604   -1.68726442
H                 -0.29735008   -5.84415908   -2.29309142
C                  0.03161546   -2.98966237   -1.68726442
N                  3.11249628   -1.28923817   -1.74366744
N                  0.76494518   -1.84674103   -1.45160441
N                 -1.28923817   -3.11249628   -1.74366744
C                  2.98966237    0.03161546   -1.68726442
C                 -2.13636604   -2.09165504   -1.68726442
N                  1.84674103    0.76494518   -1.45160441
C                  4.09455913    0.93505629   -1.98912040
N                 -1.84674103   -0.76494518   -1.45160441
C                 -3.55647518   -2.23410588   -1.98912040
C                  2.13636604    2.09165504   -1.68726442
C                  5.44022282    0.70999313   -2.28518040
C                  3.55647518    2.23410588   -1.98912040
C                 -2.98966237   -0.03161546   -1.68726442
C                 -4.34885941   -3.34477749   -2.28518040
C                 -4.09455913   -0.93505629   -1.98912040
N                  1.28923817    3.11249628   -1.74366744
C                  6.23510460    1.82074393   -2.56314038
H                  5.84415908   -0.29735008   -2.29309142
C                  4.34885941    3.34477749   -2.28518040
N                 -3.11249628    1.28923817   -1.74366744
C                 -5.69634512   -3.12142436   -2.56314038
H                 -3.92218626   -4.34270278   -2.29309142
C                 -5.44022282   -0.70999313   -2.28518040
C                 -0.03161546    2.98966237   -1.68726442
C                  5.69634512    3.12142436   -2.56314038
H                  7.28940820    1.68395303   -2.78817738
H                  3.92218626    4.34270278   -2.29309142
C                 -2.09165504    2.13636604   -1.68726442
C                 -6.23510460   -1.82074393   -2.56314038
H                 -6.34512458   -3.96365537   -2.78817738
H                 -5.84415908    0.29735008   -2.29309142
C                 -0.93505629    4.09455913   -1.98912040
N                 -0.76494518    1.84674103   -1.45160441
H                  6.34512458    3.96365537   -2.78817738
C                 -2.23410588    3.55647518   -1.98912040
H                 -7.28940820   -1.68395303   -2.78817738
C                 -0.70999313    5.44022282   -2.28518040
C                 -3.34477749    4.34885941   -2.28518040
C                 -1.82074393    6.23510460   -2.56314038
H                  0.29735008    5.84415908   -2.29309142
C                 -3.12142436    5.69634512   -2.56314038
H                 -4.34270278    3.92218626   -2.29309142
H                 -1.68395303    7.28940820   -2.78817738
H                 -3.96365537    6.34512458   -2.78817738
Ce                 0.00000000    0.00000000    0.00000000

C H 0
6-31g(d)
****
101 0
6-311+g(3d2f)
****
14 0
6-311+g(3d2f)
****
18 0
6-311+g(3d2f)
****
20 0
6-311+g(3d2f)
****
45 0
6-311+g(3d2f)
****
70 0
6-311+g(3d2f)
****
74 0
6-311+g(3d2f)
****
76 0
6-311+g(3d2f)
****
13 0
6-31g(d)
****
15 0
6-31g(d)
****
32 0
6-31g(d)
****
28 0
6-31g(d)
****
69 0
6-31g(d)
****
71 0
6-31g(d)
****
88 0
6-31g(d)
****
84 0
6-31g(d)
****
Ce 0
sddall
****

Ce 0
sddall
===gjf文件结束===

[ Last edited by yjcmwgk on 2009-8-24 at 09:35 ]
赞一下(2人) 回复此楼
2楼2009-08-24 09:19:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
%KJob LN [M]
指示程序在执行Link N的第M个事件之后停止。例如%KJob L502 2表示执行模块502两次后停止。如果M值省略，默认值为1。

引自高斯说明书zork中译版
赞一下(1人) 回复此楼
3楼2009-08-24 09:37:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
loovfnd(金币+3):感谢帮助
yjcmwgk(金币+1,VIP+0):看看7楼和8楼吧，你和我都说话太草率了，不过，你说的更接近正确答案。谢谢你 8-24 13:23
引用回帖:
Originally posted by loovfnd at 2009-8-24 09:05:
1. aug-cc-pVDZ，
2. aug-cc-pVTZ，
3. aug-cc-pVQZ，
4. aug-cc-pV5Z，
5. 6-311++g（3df）

从大到小排列
43251
赞一下(8人) 回复此楼
据权威报纸调查，84%的青少年犯罪者有玩电子游戏的经历，故应该禁止经营电子游戏。据我们所知，100%的青少年犯罪者有穿鞋的经历，故制鞋厂应该停业。
4楼2009-08-24 10:48:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loovfnd

至尊木虫 (著名写手)
★
yjcmwgk(金币+1,VIP+0):看看7楼和8楼吧，我和hairan都说话草率了 8-24 13:21
哪个正确呢？

yjcmwgk   4>3>2>1>5
hairan    4>3>2>5>1

[ Last edited by yjcmwgk on 2009-8-24 at 12:53 ]
赞一下(6人) 回复此楼
莫强求！
5楼2009-08-24 11:10:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+2,VIP+0):感谢交流 8-24 13:21
从形式上看，6-311++g（3df）是包含了极化和弥散函数的triple-zeta基组，因此和aug-cc-pVTZ相当。但是一般情况下比后者的函数数量要少一些，特别是对于过渡金属，非常明显。此外，6-311++g（3df）是分块收缩的，计算相关能的时候没有收敛极限。
赞一下(7人) 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
6楼2009-08-24 12:50:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
loovfnd(金币+2):感谢
对于第一周期元素的测试：
6-311++g(3df) for H atom：
4 basis functions, 6 primitive gaussians
aug-cc-pvdz for H atom:
9 basis functions, 11 primitive gaussians
对于第二周期元素的测试：
6-311++g(3df) for C atom：
39 basis functions, 58 primitive gaussians
aug-cc-pvdz for C atom:
23 basis functions, 43 primitive gaussians
对于第三周期元素的测试：
6-311++g(3df) for Si atom:
47 basis functions, 72 primitive gaussians
aug-cc-pvdz for Si atom:
27 basis functions, 83 primitive gaussians
对于第四周期元素的测试：
6-311++g(3df) for Fe atom:
81 basis functions, 129 primitive gaussians
aug-cc-pvdz for Si atom:
超出Gaussian03D.01版的aug-CC-pVDZ基组的范围了
=========
结论：你很难说这两个基组的规模谁比谁大
我当时说aug-cc-pvdz比6-311++g(3df)规模大，这句话实在是太草率了。I am sorry!
赞一下(2人) 回复此楼
7楼2009-08-24 13:18:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
测试文件如下，大家可以自己运行一下看看：
===gjf文件开始===
%kjob L302
# hf/6-311++g(3df)

a

0 2
H

--Link1--
%kjob L302
# hf/aug-cc-pvdz

a

0 2
H

--Link1--
%kjob L302
# hf/6-311++g(3df)

a

0 5
C

--Link1--
%kjob L302
# hf/aug-cc-pvdz

a

0 5
C

--Link1--
%kjob L302
# hf/6-311++g(3df)

a

0 5
Si

--Link1--
%kjob L302
# hf/aug-cc-pvdz

a

0 5
Si

--Link1--
%kjob L302
# hf/6-311++g(3df)

a

0 3
Fe

--Link1--
%kjob L302
# hf/aug-cc-pvdz

a

0 3
Fe

===gjf文件结束===

[ Last edited by yjcmwgk on 2009-8-24 at 13:25 ]
赞一下 回复此楼
8楼2009-08-24 13:20:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqx

木虫 (正式写手)

我来钻个牛角尖

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
loovfnd(金币+3,VIP+0):感谢你的牛角尖 8-24 13:49
1，2，3，4的顺序是绝对的；5和1不一定，和具体体系可能有关系：如果体系只有氢原子，则5>1，如氢分子，二者大小分别为18，8; 如果体系中氢原子和碳原子个数比例为>=4:1，则5>=1; 如在CH4分子中，5=1 (二者的基函数个数均为59).
感兴趣的虫子可参考我在chem8的一个贴子：
http://hk.chem8.org/bbs/viewthread.php?tid=29962&highlight=
赞一下(11人) 回复此楼
9楼2009-08-24 13:38:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 loovfnd 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭