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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】关于高斯的基组，审稿意见让用triple-zeta基组【完结】
231/3返回列表123下一页
查看: 3829  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhangguangping

木虫 (著名写手)

[交流] 【求助】关于高斯的基组，审稿意见让用triple-zeta基组【完结】 已有10人参与

在文章的计算中，我们体系包含H,S,C以及Au。使用的是LAN2DZ基组。审稿意见让用trip-zeta基组。LAN2DZ这个基组对H,S,C来说是一个价键DZ基组。对Au说是Los Alamos ECP加上DZ。
所以想请问一下：
（1）Los Alamos ECP是什么？是赝势吗？
（2）要想达到审稿人的要求，我们使用混合基组，对H,S,C使用cc-pVTZ基组,对Au使用CEP-121G基组是不是可以的？

谢谢大家的回复！

谢谢大家的讨论。这个帖子已经完结了。这个帖子还有更多的回帖和讨论的。但是现在只有30多个回帖了。问题的解决都在哪些回帖中。可能是小木虫的服务器更新的时候，将数据丢失了一些吧。

下面是我们计算时采用的输入文件：
#P Becke3LYP/GenECP opt=(z-matrix,Maxcyc=202) GFInput Nosymm scf=(Maxcyc=2000,conver=6)


体系的结构说明

C H S 0   
6-311G
****      
Au 0
S   1   1.00
      2.8090000              1.0000000        
S   1   1.00
      1.5950000              1.0000000        
S   1   1.00
      0.5327000              1.0000000        
S   1   1.00
      0.2826000              1.0000000        
S   1   1.00
      0.0598000              1.0000000        
P   1   1.00
      3.6840000              1.0000000        
P   1   1.00
      1.6660000              1.0000000        
P   1   1.00
      0.5989000              1.0000000        
P   1   1.00
      0.0977000              1.0000000        
P   1   1.00
      0.0279000              1.0000000        
D   1   1.00
      1.2870000              1.0000000        
D   1   1.00
      0.4335000              1.0000000        
D   1   1.00
      0.1396000              1.0000000
****

AU 0
AU-ECP     4     60
g potential
  5
1    622.6287956            -60.0000000        
2    136.2843607           -555.5292312        
2     33.1549781           -168.0019785        
2      9.9894895            -63.0399875        
2      3.0481312             -4.2516681        
s-g potential
  6
0    194.7374304              3.0000000        
1    351.5327447             38.6020880        
2    122.3270402            864.8370727        
2     32.0914617            374.9935520        
2      5.2451812            289.7910100        
2      4.4916223           -152.4532773        
p-g potential
  4
0    420.6158801              2.0000000        
1    109.4417815             73.8885625        
2     34.1714280            326.6729872        
2      5.9879750            126.5814591        
d-g potential
  5
0    219.2666158              3.0000000        
1    122.7297786             55.6793149        
2     63.1063369            449.1987335        
2     18.3684520            215.0269091        
2      4.4972844             64.0840995        
f-g potential
  5
0    108.5506037              4.0000000        
1     56.4795527             51.8065335        
2     29.2069159            231.2183113        
2      9.5440543            119.0047386        
2      2.8965118             15.3424188        

采用的都是TZ基组。

[ Last edited by zhangguangping on 2011-9-1 at 10:12 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
弘德明志博学笃行
1楼 2011-03-01 10:08:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
建议使用 aug-cc-pVTZ for H, S, C
aug-cc-pVTZ-pp for Au + Stuttgart small core relativistic ECP for Au (60mwb)


Lan2DZ 是多年前的东西了
赞一下(1人) 回复此楼
更痛苦的是梦没醒路已经不能走了
2楼2011-03-01 10:50:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
Originally posted by mchen10 at 2011-03-01 03:50:24:
建议使用 aug-cc-pVTZ for H, S, C
aug-cc-pVTZ-pp for Au + Stuttgart small core relativistic ECP for Au (60mwb)


Lan2DZ 是多年前的东西了

aug-cc-pVTZ for H, S, C
这一个我知道是在我的那个基础上加上极化和弥散了。
但是
aug-cc-pVTZ-pp for Au + Stuttgart small core relativistic ECP for Au (60mwb)
这一个是什么基组啊（加了相对论?）？没有查到啊。他的关键字是什么呢？请指教！
赞一下 回复此楼
弘德明志博学笃行
3楼2011-03-01 11:16:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by zhangguangping at 2011-03-01 11:16:34:
aug-cc-pVTZ for H, S, C
这一个我知道是在我的那个基础上加上极化和弥散了。
但是
aug-cc-pVTZ-pp for Au + Stuttgart small core relativistic ECP for Au (60mwb)
这一个是什么基组啊（加了相对论?）？ ...

genecp手动输入

basis set 和 ecp 可以在 emsl basis set exchange找到，选gaussian94格式的下载就行了

里面有海量的basis set 挑合适的用。。
赞一下(1人) 回复此楼
更痛苦的是梦没醒路已经不能走了
4楼2011-03-01 11:41:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
Originally posted by mchen10 at 2011-03-01 04:41:53:
genecp手动输入

basis set 和 ecp 可以在 emsl basis set exchange找到，选gaussian94格式的下载就行了

里面有海量的basis set 挑合适的用。。

谢谢你的指点，这个对我来说难度有点高。
#P Becke3LYP/Gen opt=z-matrix  Nosymmetry
scf=(Maxcyc=2000 conver=6)
.
.
.
C H S 0   
cc-pVTZ
****      
Au 0      
CEP-121G
****

Au 0      
CEP-121G

我先使用这个试试
估计你的那个基组对我来说内存都是问题。谢谢你了。如果需要再请教你！
赞一下 回复此楼
弘德明志博学笃行
5楼2011-03-01 11:53:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
Originally posted by mchen10 at 2011-03-01 04:41:53:
genecp手动输入

basis set 和 ecp 可以在 emsl basis set exchange找到，选gaussian94格式的下载就行了

里面有海量的basis set 挑合适的用。。

看来你是一个高斯高手。呵呵呵！
赞一下 回复此楼
弘德明志博学笃行
6楼2011-03-01 11:56:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
你也可以尝试使用2008年开发的LANL2TZ系列
赞一下(1人) 回复此楼
前途光明，出路难觅！
7楼2011-03-01 15:01:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
Originally posted by lihb734 at 2011-03-01 08:01:48:
你也可以尝试使用2008年开发的LANL2TZ系列

为什么我用五楼的那个设置，计算过程中报告警告如下呢:
  Basis read from rwf:  (5D, 7F)
No pseudopotential information found on rwf file.
Integral buffers will be    131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned off.
   916 basis functions,  1755 primitive gaussians,  1039 cartesian basis functions
   324 alpha electrons      324 beta electrons
       nuclear repulsion energy     15996.3487521324 Hartrees.



Warning!  Au atom   30 has 79 valence electrons but only  31 basis functions.
This is less than a minimal basis set!
赞一下 回复此楼
弘德明志博学笃行
8楼2011-03-01 16:07:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
审稿人是让你非金属元素用TZ还是金属元素？感觉Au用DZ已经可以了，不过你用TZ的话计算量要大很多。
赞一下(1人) 回复此楼
Anewday,anewhour,anewminute,anewpeople.
9楼2011-03-01 16:21:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
Originally posted by qzhaosdu at 2011-03-01 09:21:57:
审稿人是让你非金属元素用TZ还是金属元素？感觉Au用DZ已经可以了，不过你用TZ的话计算量要大很多。

就是啊。没有说。仅仅说是DZ有点小，建议用TZ算一下。呵呵呵
赞一下 回复此楼
弘德明志博学笃行
10楼2011-03-01 16:33:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangguangping 的主题更新
231/3返回列表123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 68gx6bh3ei 2026-07-09 3/150 2026-07-10 06:33 by pbuc7it0i8
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ph7bxbho8f 2026-07-09 3/150 2026-07-10 05:53 by j4frgf6zcw
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[基金申请] b口会评 +6 当空照 2026-07-08 6/300 2026-07-09 20:49 by tonyliu0401
[基金申请] 信息学部会评时间 +3 xuel2011 2026-07-06 3/150 2026-07-09 19:05 by mollyzhang_2003
[基金申请] 祈祷青基必中 50+3 hadengry 2026-07-09 9/450 2026-07-09 17:37 by 文杰猪猪
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
[论文投稿] MSER送审了还被拒稿 +3 S778950906 2026-07-08 5/250 2026-07-09 15:46 by S778950906
[基金申请] 关于如何从代码看上不上会 +11 且听虎啸 2026-07-09 14/700 2026-07-09 15:27 by libeny
[有机交流] 硝基还原 20+3 暮年飞雪 2026-07-07 5/250 2026-07-09 13:30 by tfang
[基金申请] 化学口会评？应用基础研究和基础研究，希望能分开 +3 Tide man 2026-07-08 3/150 2026-07-09 11:36 by dede5556
[基金申请] 今年E04面上 +15 布布和一二 2026-07-05 21/1050 2026-07-09 10:25 by 学员fOzRO9
[硕博家园] 大龄残疾硕士的一点执念 +14 3776737374 2026-07-06 22/1100 2026-07-09 06:07 by lovechizi
[基金申请] 生命口会评 +4 wuke100666 2026-07-07 4/200 2026-07-08 16:38 by 逍遥不够自在
[教师之家] 网络举报学术不端成为了男女之间情杀的手段了，真荒唐。。。 +3 zju2000 2026-07-04 3/150 2026-07-07 21:00 by hyx552200
[有机交流] Suzuki偶联中硼酸酯原料掉下的硼酸酯是否可以自身偶联，生成B2Pin2杂质 500+3 天下大同 2026-07-06 3/150 2026-07-07 10:31 by martyr_oyx
[论文投稿] 刑法学论文投稿求助 +3 高源678 2026-07-04 5/250 2026-07-06 12:20 by 平天尺
[论文投稿] Journal of Environmental Chemical Engineering +3 雪花天涯 2026-07-05 3/150 2026-07-06 09:00 by bobvan
[基金申请] 最后一年，祈求好运 +3 pyw1992 2026-07-06 3/150 2026-07-06 08:53 by 苏_frp
信息提示
关闭
请填处理意见
关闭