24小时热门版块排行榜    

查看: 1048  |  回复: 8

wynli

铁虫 (小有名气)

[求助] 结合能比文献上的大很多

# B3LYP counterpoise=2

Title Card Required

1 1 1 1 0 1
Cu                    0.00011   0.00009  -0.00002  1
C                     0.8911    1.62704   0.76534  2
C                    -1.44342   0.5651   -1.27456  1
C                     1.35789  -1.10946  -0.97599  1
C                    -0.80559  -1.08269   1.48515  1
O                     2.12517  -1.73651  -1.52736  1
O                    -1.2611   -1.69454   2.32423  1
O                     1.39469   2.54638   1.19791  2
O                    -2.25915   0.88436  -1.99467  1

C     0
S   6   1.00
   4232.6100000              0.0020290        
    634.8820000              0.0155350        
    146.0970000              0.0754110        
     42.4974000              0.2571210        
     14.1892000              0.5965550        
      1.9666000              0.2425170        
S   1   1.00
      5.1477000              1.0000000        
S   1   1.00
      0.4962000              1.0000000        
S   1   1.00
      0.1533000              1.0000000        
P   4   1.00
     18.1557000              0.0185340        
      3.9864000              0.1154420        
      1.1429000              0.3862060        
      0.3594000              0.6400890        
P   1   1.00
      0.1146000              1.0000000        
D   1   1.00
      0.7500000              1.0000000        
****
O     0
S   6   1.00
   7816.5400000              0.0020310        
   1175.8200000              0.0154360        
    273.1880000              0.0737710        
     81.1696000              0.2476060        
     27.1836000              0.6118320        
      3.4136000              0.2412050        
S   1   1.00
      9.5322000              1.0000000        
S   1   1.00
      0.9398000              1.0000000        
S   1   1.00
      0.2846000              1.0000000        
P   4   1.00
     35.1832000              0.0195800        
      7.9040000              0.1241890        
      2.3051000              0.3947270        
      0.7171000              0.6273750        
P   1   1.00
      0.2137000              1.0000000        
D   1   1.00
      0.8500000              1.0000000        
****
Cu     0
S   6   1.00
337200.0000000              0.0002800        
  50072.9000000              0.0021700        
  11373.4000000              0.0112200        
   3239.8200000              0.0446210        
   1071.9700000              0.1383600        
    395.0990000              0.3125500        
S   2   1.00
    158.3990000              0.4106000        
     67.3591000              0.2145300        
S   1   1.00
     22.2983000              1.0000000        
S   1   1.00
      9.3935700              1.0000000        
S   1   1.00
      2.5784800              1.0000000        
S   1   1.00
      0.9640800              1.0000000        
S   1   1.00
      0.1133030              1.0000000        
S   1   1.00
      0.0407910              1.0000000        
P   3   1.00
   2245.2900000              0.0022600        
    532.1060000              0.0184100        
    172.1950000              0.0864600        
P   3   1.00
     65.3239000              0.2545000        
     27.0551000              0.4341900        
     11.7435000              0.3339000        
P   1   1.00
      4.6938200              1.0000000        
P   2   1.00
      1.9066700              0.5693600        
      0.7114450              0.2470200        
P   1   1.00
      0.2325980              1.0000000        
P   1   1.00
      0.0692990              1.0000000        
D   3   1.00
     48.5439000              0.0309700        
     13.5490000              0.1619600        
      4.5162800              0.3781200        
D   1   1.00
      1.4732900              1.0000000        
D   1   1.00
      0.4148750              1.0000000        
D   1   1.00
      0.1491000              1.0000000        
F   3   1.00
      8.0267373              0.1737857        
      2.8785728              0.5973381        
      1.2326390              0.3929396        
****

输入文件是这样的,最后结果是58.3kcal/mol,文献里算出的是17.01kcal/mol
差别很大啊,我跟文献唯一不同的是我算的第二个配体CO的结合能。文献里算的是最后一个配体的结合能。(这里四个配体的结合能我都算了,都比17.01大很多。) 怎么个情况
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gujing136

禁虫 (文坛精英)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+3 2012-09-24 19:01:27
话说你用DFT的貌似精度比较低。
可以用稍微高点的基组么?
叫老猫的兔子
2楼2012-09-24 17:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

不可理喻_

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
帮顶下,不是很清楚!
加油,GOGO!!!
3楼2012-09-24 19:17:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★
gmy1990: 金币+2 2012-09-24 20:31:00
如果输入文件不是重新编辑过,那么用的是STO-3G基组
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2012-09-24 19:54:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wynli

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by beefly at 2012-09-24 19:54:51
如果输入文件不是重新编辑过,那么用的是STO-3G基组

后边加gen
是吧
5楼2012-09-24 20:02:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CKX

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
弱弱地问下你文件后面的C,O,Cu的机组信息哪来的?正确吗?为什么不直接用6-311++Gdp?

[ 发自手机版 http://muchong.com/3g ]
6楼2012-09-25 07:51:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gengcy

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
问题没理解,是按照文献计算方法计算得到的结果不一样,还是和文献中的实验值比较后结果不一样?
分析
1.我相信文献不会用楼主给的那种方法来算,基组问题,太不可靠。
2.没仔细看构型,如果b3lyp算的是分子间作用,那么肯定和实验值差距非常大。
3.出现金属原子,是不是该赝势?不过,通常b3lyp算金属,比较麻烦,Cu应该还好,要是Fe......你懂的
7楼2012-09-25 09:44:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

哈哈大姑爷

至尊木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wynli: 金币+15 2012-11-10 10:44:07
如果输入文件不是重新编辑过,那么用的是STO-3G基组
GO ON
8楼2012-09-25 10:08:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wynli

铁虫 (小有名气)

引用回帖:
8楼: Originally posted by 哈哈大姑爷 at 2012-09-25 10:08:30
如果输入文件不是重新编辑过,那么用的是STO-3G基组

什么意思
9楼2012-09-25 10:13:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wynli 的主题更新
信息提示
请填处理意见