24小时热门版块排行榜    

查看: 1218  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思:在爱情中有过遗憾吗? +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错,点赞! +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题,请指教! 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
请填处理意见