24小时热门版块排行榜    

查看: 1097  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 小窗夜梦 的主题更新
信息提示
请填处理意见