24小时热门版块排行榜    

查看: 1395  |  回复: 4

phoeby_ice

金虫 (小有名气)

[求助] 求助,帮忙计算一下

小女正用高斯计算一分子的PA,结果总是很大(>500kcal/mol),一般300多为合理,求助各位大神帮忙算下,我哪里弄错了?
PA=H(ArO-)+H(H+)-H(ArOH)       H(H+)=6.197kj/mol(文献值)
即分子失去一个氢离子的能量

代表上式中ArOH的结构优化输入文件:
%chk=ra.chk
%mem=800MW
%nprocshared=8
# opt freq b3lyp/gen geom=connectivity pseudo=read

ra

0 1
C                 -2.30570368   -1.77035906    0.00026300
C                 -0.93290268   -1.77035906    0.00026300
C                 -0.21096468   -0.54523806    0.00026300
C                 -0.92883468    0.67862794    0.00063300
C                 -2.35025068    0.64925494    0.00078700
C                 -3.02175568   -0.54872006    0.00049700
H                 -2.86643468   -2.71670506    0.00018300
C                  1.21076232   -0.51569506    0.00000000
H                 -2.89749168    1.60371094    0.00086900
C                  1.16677632    1.90342194    0.00047500
C                  1.88278832    0.68152994    0.00001600
O                  3.31237433    0.71593612   -0.00022061
H                  3.62516914    1.03220751   -0.85094490
O                 -0.19889857   -2.99760715   -0.00003331
H                 -0.81101858   -3.73710776    0.00689173
O                 -4.45134477   -0.58299862    0.00062961
O                 -0.20608468    1.90335394    0.00074800
O                  1.92211427   -1.75621044    0.00002599
C                  1.95177175    3.22833028    0.00022579
C                  1.26780682    4.44433325    0.00066890
C                  3.34651314    3.21270119   -0.00033170
C                  1.97850423    5.64438754   -0.00012737
H                  0.16821686    4.45631575    0.00033114
C                  4.05757713    4.41303358   -0.00013233
H                  3.88574400    2.25439065   -0.00050109
C                  3.37384399    5.62878274   -0.00016840
H                  1.43946494    6.60289318   -0.00041217
H                  5.15726557    4.40040214   -0.00015657
O                  4.10224325    6.85936551   -0.00090787
H                  4.18849101    7.18123984    0.89940221
C                 -4.92699326   -0.65421452    1.34732438
H                 -5.99590185   -0.70251120    1.34596448
H                 -4.53104280   -1.52887632    1.81964840
H                 -4.60993976    0.21545659    1.88402657

1 2 1.5 6 1.5 7 1.0
2 3 1.5 14 1.0
3 4 1.5 8 1.0
4 5 1.5 17 1.0
5 6 1.5 9 1.0
6 16 1.0
7
8 11 1.0 18 2.0
9
10 11 2.0 17 1.0 19 1.0
11 12 1.0
12 13 1.0
13
14 15 1.0
15
16 31 1.0
17
18
19 20 1.5 21 1.5
20 22 1.5 23 1.0
21 24 1.5 25 1.0
22 26 1.5 27 1.0
23
24 26 1.5 28 1.0
25
26 29 1.0
27
28
29 30 1.0
30
31 32 1.0 33 1.0 34 1.0
32
33
34

C H 0
6-31g*
****
O 0
6-31+g*
****



优化结构后用MP2计算能量。

如果有算出上述理想的数值,请把输入输出文件发给我吧~
whatever, 小女先谢过
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

viger87

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
要不要考虑下加溶剂化模型?
有花堪折直须折,末待花落空折枝!
2楼2012-03-15 22:49:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

phoeby_ice

金虫 (小有名气)

引用回帖:
2楼: Originally posted by viger87 at 2012-03-15 22:49:51:
要不要考虑下加溶剂化模型?

谢谢你的意见,别人直接能够算出来,应该不是溶剂化的问题吧,主要是我很困惑这个问题,不想就这么绕开
3楼2012-03-16 09:57:32
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

【答案】应助回帖


感谢参与,应助指数 +1
gmy1990: 金币+1 2012-03-16 12:24:58
geom=connectivity pseudo=read这个不需要
前途光明,出路难觅!
4楼2012-03-16 10:34:50
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

phoeby_ice

金虫 (小有名气)

引用回帖:
4楼: Originally posted by lihb734 at 2012-03-16 10:34:50:
geom=connectivity pseudo=read这个不需要

试过了,删掉也算不对
还有别的方法么
5楼2012-03-22 17:18:56
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 phoeby_ice 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +13 羊腰板 2026-08-21 13/650 2026-08-23 09:25 by lfy8008
[教师之家] 跳槽后在研项目怎么办? +4 简单化xn 2026-08-22 8/400 2026-08-23 07:20 by ilovexiaobin
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见