24小时热门版块排行榜    

查看: 1449  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 小窗夜梦 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +5 kulium 2026-08-21 6/300 2026-08-21 18:15 by alongwaytogo
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 14/700 2026-08-21 16:51 by applepetty
[基金申请] 人气不行了 +7 fansofjerry 2026-08-21 7/350 2026-08-21 14:13 by Vivilian
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 18/900 2026-08-21 13:46 by 医学老男孩
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 今天基金会出结果吗?20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[基金申请] 科研孤儿太难了 +15 我4大白菜 2026-08-20 16/800 2026-08-21 11:07 by yming1319
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳今天,20号变了 +4 archvillain 2026-08-20 4/200 2026-08-21 08:10 by gatelove
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
信息提示
请填处理意见