应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 求助，帮忙计算一下
51/1返回列表
查看: 1130  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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, 小女先谢过
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-03-15 15:47:51
已阅   关注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这个不需要
赞一下(1人) 回复此楼
前途光明，出路难觅！
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 的主题更新
51/1返回列表
信息提示
关闭
请填处理意见
关闭