24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1240  |  回复: 8

小窗夜梦

金虫 (小有名气)

hero

[求助] 文献中频率不能重复

文献中用Gamess/mp2计算,我用实验室Gaussian/MP2,所用基组与文献中一样,但是其中频率重复不了,求高人指教!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

Tobethehero!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★
gmy1990: 金币+2 2013-09-29 06:41:22
如果有d,f,g基函数,可能是球函数和笛卡尔函数的差别;还有可能用了不同的冻芯设置
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2013-09-28 22:45:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luqing6879

木虫 (著名写手)

引用回帖:
2楼: Originally posted by beefly at 2013-09-28 22:45:58
如果有d,f,g基函数,可能是球函数和笛卡尔函数的差别;还有可能用了不同的冻芯设置

球函数和笛卡尔函数的差别能差多少呢?
cishan.goodweb.cn/慈善点击,点滴帮助。
3楼2013-09-29 02:39:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
gkf高: 金币+2 2013-09-29 23:22:01
相差的多么?gaussian的mp2有解析二阶导数,而gamess的只有解析1阶导数,需要数值微分,所以数值精度很有限;而且gamess的解析导数不能做高于一定角动量的基函数。如果只有能量数值导数求频率的话误差达到20cm-1也不奇怪
4楼2013-09-29 10:15:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小窗夜梦

金虫 (小有名气)

hero

引用回帖:
2楼: Originally posted by beefly at 2013-09-28 22:45:58
如果有d,f,g基函数,可能是球函数和笛卡尔函数的差别;还有可能用了不同的冻芯设置

有什么解决的办法没有?
Tobethehero!
5楼2013-09-29 10:43:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小窗夜梦

金虫 (小有名气)

hero

引用回帖:
4楼: Originally posted by virtualzx at 2013-09-29 10:15:02
相差的多么?gaussian的mp2有解析二阶导数,而gamess的只有解析1阶导数,需要数值微分,所以数值精度很有限;而且gamess的解析导数不能做高于一定角动量的基函数。如果只有能量数值导数求频率的话误差达到20cm-1也不 ...

The gold atom was treated with 19 explicit electrons using the (SBKJC)
85 effective core potential (ECP) with scalar relativistic
corrections augmented with a 3f/2g set of polarization functions
with f function exponents of 2.00, 0.84 and 0.31, and g function
exponents of 1.90 and 0.69, and one set of s and p diffuse
functions with exponents of 0.01. Spherical Gaussian functions
were used.

这个里面Au的基组是怎样定义的?
Tobethehero!
6楼2013-09-29 15:49:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小窗夜梦

金虫 (小有名气)

hero

引用回帖:
2楼: Originally posted by beefly at 2013-09-28 22:45:58
如果有d,f,g基函数,可能是球函数和笛卡尔函数的差别;还有可能用了不同的冻芯设置

首先非常感谢您的回复。因为我现在做的是要解释频移顺序,所以必须要算到正确的顺序,才行。文中用到的基组是这样自定义的:

The gold atom was treated with 19 explicit electrons using the (SBKJC)
85 effective core potential (ECP) with scalar relativistic
corrections augmented with a 3f/2g set of polarization functions
with f function exponents of 2.00, 0.84 and 0.31, and g function
exponents of 1.90 and 0.69, and one set of s and p diffuse
functions with exponents of 0.01. Spherical Gaussian functions
were used.

Au     0
SP   4   1.00
      7.4190000              0.2225460              0.0199240        
      4.0230000             -1.0860450             -0.2999970        
      1.6980000              1.1560390              0.7489190        
      0.6271000              0.5180610              0.5040230        
SP   1   1.00
      0.1515000              1.0000000              1.0000000        
SP   1   1.00
      0.0492500              1.0000000              1.0000000        
SP   1   1.00
      1.5020000              1.0000000              1.0000000
SP   1   1.00
         0.01                1.0000000              1.0000000        
D   3   1.00
      3.6300000             -0.0874020        
      1.9120000              0.4686340        
      0.8423000              0.6548050        
D   1   1.00
      0.3756000              1.0000000        
D   1   1.00
      0.1544000              1.0000000   
F   1   1.00
        2.00                 1.0000000
F   1   1.00
        0.84                 1.0000000
F   1   1.00
        0.31                 1.0000000
G   1   1.00
        1.90                 1.0000000
G   1   1.00
        0.69                 1.0000000     
****

AU     0
AU-ECP     4     60
g-ul potential
  1
1      4.3876300            -10.7235800        
s-ul potential
  3
0      1.5563600              6.3561200        
2      3.7159300           -364.4403900        
2      4.0679200            428.1975300        

这个基组是我在https://bse.pnl.gov/bse/portal——Au(Gaussian94)——SBKJC VDZ ECP的基础上修改的。不知对不对,能否请您帮忙看看?在此,先谢谢您了!
Tobethehero!
7楼2013-09-29 16:46:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

ecp不全
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
8楼2013-09-29 22:31:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小窗夜梦

金虫 (小有名气)

hero

引用回帖:
8楼: Originally posted by beefly at 2013-09-29 22:31:39
ecp不全

不好意思,我贴错了。应该是:

****

AU     0
AU-ECP     4     60
g-ul potential
  1
1      4.3876300            -10.7235800        
s-ul potential
  3
0      1.5563600              6.3561200        
2      3.7159300           -364.4403900        
2      4.0679200            428.1975300        
p-ul potential
  3
0      1.1879800              4.4151800        
2      3.0155100           -136.5755000        
2      3.5958800            194.2053500        
d-ul potential
  2
0     35.2500000              8.8819800        
2      5.0230700             86.7661200        
f-ul potential
  1
0      1.6888100              6.2160300        

您能帮我看看我用的基组自定义,对了吗?

另外,您对简谐振动频率和非简谐振动频率了解吗,可否介绍一下资料?我对这块不是太懂
Tobethehero!
9楼2013-09-30 11:18:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小窗夜梦 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿厦大生物学332求调剂 +5 池池池池池池 2026-04-08 5/250 2026-04-08 21:23 by zhyzzh
[考研] 307分材料专业求调剂 +12 Hll胡 2026-04-05 12/600 2026-04-08 16:33 by luoyongfeng
[考研] 336求调剂,一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 277求调剂 +4 考研调剂lxh 2026-04-06 6/300 2026-04-08 10:40 by 逆水乘风
[考研] 22408 一志愿双一流人工智能300分 四六级,数据分析国奖 +4 zzfeng123 2026-04-06 6/300 2026-04-07 21:02 by zzfeng123
[考博] 博士申请 +3 IQwQl 2026-04-05 3/150 2026-04-07 20:31 by greychen00
[考研] 277求调剂 数一104分 +9 瓶子PZ 2026-04-05 14/700 2026-04-07 17:52 by 蓝云思雨
[考研] 生物学363调剂求助 +7 fanzhang6666 2026-04-06 9/450 2026-04-07 17:37 by lijunpoly
[考研] 302分求调剂 一志愿安徽大学085601 +12 zyx上岸! 2026-04-04 12/600 2026-04-07 02:09 by BruceLiu320
[考研] 华南理工0703化学,总分336求调剂 +10 手机用户 2026-04-03 11/550 2026-04-06 22:20 by michael2011
[考研] 调剂 一志愿吉林大学357分 +5 .Starry. 2026-04-04 5/250 2026-04-06 09:28 by cql1109
[考研] 材料专硕283求调剂 +17 试试看呗 2026-04-04 18/900 2026-04-06 09:24 by 286640313
[考研] 考研生物学考A区211,初试322,科目生化和生物综合,求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[考研] 生物学308分求调剂(一志愿华东师大) +8 相信必会光芒万 2026-04-05 10/500 2026-04-05 12:19 by Hdyxbekcb
[考研] 考研调剂 +3 mcbbc 2026-04-04 3/150 2026-04-05 10:03 by barlinike
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 22:08 by 啵啵啵0119
[考研] 工科341分调剂 +3 洛多罗 2026-04-03 3/150 2026-04-03 14:20 by 1753564080
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
信息提示
请填处理意见