24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1233  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 小窗夜梦 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿南航,数一英一学硕317求调剂!! +5 Acaciad 2026-04-04 5/250 2026-04-05 12:31 by 搏击518
[考研] 一志愿北京化工大学,初试成绩350求调剂 +9 沿岸?贝壳 2026-04-04 14/700 2026-04-05 01:09 by 沿岸?贝壳
[考研] 材料专硕322分 +11 哈哈哈吼吼吼哈 2026-04-02 11/550 2026-04-04 23:37 by 永字号
[考研] 材料调剂 +11 一样YWY 2026-04-02 13/650 2026-04-04 23:10 by 无际的草原
[考研] 278求调剂 +14 范婷娜 2026-04-04 15/750 2026-04-04 22:15 by lqwchd
[考研] 材料调剂 +18 一样YWY 2026-04-02 19/950 2026-04-04 22:14 by hemengdong
[考研] [调剂信息]085408光电信息 求调剂 总分291分数一英一 +3 iz11az 2026-04-02 3/150 2026-04-04 19:09 by 蓝云思雨
[考研] 一志愿北交大材料工程总分358 +6 cs0106 2026-04-03 6/300 2026-04-04 11:20 by w_xuqing
[考研] 322求调剂 +6 FZAC123 2026-04-03 6/300 2026-04-03 22:23 by 科研小专家
[考研] 学硕288调剂!!! +3 小王xw123 2026-04-03 3/150 2026-04-03 21:20 by 啵啵啵0119
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
[考研] 一志愿厦门大学材料工程专硕354找调剂!!! +8 贝呗钡钡 2026-03-30 8/400 2026-04-03 09:41 by hypershenger
[考研] 279求调剂 +5 傅文秋 2026-04-02 5/250 2026-04-02 18:10 by 笔落锦州
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 一志愿北交大材料工程总分358 +8 cs0106 2026-04-01 9/450 2026-04-02 10:36 by 不吃魚的貓
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 00:21 by 百秒光年
[考研] 0710生物学求调剂! +6 叙述文 2026-03-31 6/300 2026-04-01 09:39 by JourneyLucky
[考研] 339求调剂 +5 zjjkt 2026-03-31 5/250 2026-04-01 09:18 by JourneyLucky
[考研] 323分 食品与营养调剂 +3 嘿ooo 2026-03-31 3/150 2026-03-31 09:38 by longlotian
信息提示
请填处理意见