24小时热门版块排行榜    

查看: 1180  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 小窗夜梦 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 化学工程321分求调剂 +7 大米饭! 2026-03-15 7/350 2026-03-16 10:25 by 了了了了。。
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考研] 311求调剂 +3 26研0 2026-03-15 3/150 2026-03-15 09:12 by JourneyLucky
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 云南财经大学信息学院计算机学硕专硕学位点 +3 zjptai 2026-03-10 5/250 2026-03-14 01:23 by 飞行琦
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 纺织、生物、化学、材料相关专业招生了 +4 耶耶业 2026-03-09 7/350 2026-03-12 19:05 by Equinoxhua
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
请填处理意见