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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】关于高斯的基组，审稿意见让用triple-zeta基组【完结】
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
131/2返回列表12下一页
查看: 3329  |  回复: 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的回帖

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的回帖

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的回帖

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的回帖

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

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

谢谢你的指点！
回复此楼
弘德明志博学笃行
13楼2011-03-01 16:52:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
Originally posted by qzhaosdu at 2011-03-01 09:45:42:
你是不是所有的元素都用的lanl2dz啊？你可以非金属用6-311之类的，然后 Au用lanl2dz

我们Au的个数不是很多。所以TZ能够接受。如果用赝势的话，只计算19个电子就行。
赞一下 回复此楼
弘德明志博学笃行
14楼2011-03-01 17:05:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
Originally posted by qzhaosdu at 2011-03-01 09:45:42:
你是不是所有的元素都用的lanl2dz啊？你可以非金属用6-311之类的，然后 Au用lanl2dz

cc-pTZ本身就包含着极化函数，并且极化的还不少。可能会加重计算量。
赞一下 回复此楼
弘德明志博学笃行
16楼2011-03-01 18:01:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
Originally posted by 求学之子 at 2011-03-01 10:16:38:
DZ主要是对第一周期的过渡金属，审稿人叫你用那个是因为那个基组对第二三周期的过渡金属算出来的结果更好，这都是大家默认的，这也是基组适用范围的问题，不想退稿就按审稿人说的做比较好！

谢谢指教！
回复此楼
弘德明志博学笃行
17楼2011-03-01 18:05:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangguangping 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭