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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】结构优化和频率计算后的能量问题
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
131/2返回列表12下一页
查看: 966  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

aaacjj

铜虫 (初入文坛)

[交流] 【求助】结构优化和频率计算后的能量问题

请教各位虫友，在对一个含一个氯四苯环和一个汞原子进行优化和频率计算后，出来的能量比别人的小，比如别人的是-615.36325a.u，而我的计算出来的确是-1117.45668644a.u，相差的比较大，是什么问题啊？
我采用的是B3LYP方法，C,H,Cl用的是6-31G*，汞用的是Lanl2dz基组
回复此楼

» 猜你喜欢
1楼2009-11-11 10:04:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfan2192

金虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+1,VIP+0): 11-11 13:09
这个问题 好像有好多人讨论
你比较的是绝对能量吧
应该统计相对能量（比如说原子键合能）
赞一下(11人) 回复此楼
相逢是缘！
2楼2009-11-11 11:22:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaacjj

铜虫 (初入文坛)
引用回帖:
Originally posted by zhangfan2192 at 2009-11-11 11:22:
这个问题 好像有好多人讨论
你比较的是绝对能量吧
应该统计相对能量（比如说原子键合能）

绝对能量指的就是高斯最后直接算出来那个数值咯？相对能量这个不太懂啊，能不能详细的说下或者有什么资料提供下，小虫在此谢过
赞一下 回复此楼
3楼2009-11-11 12:49:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+1,VIP+0):感谢交流 11-11 19:29
很明显，你的赝势写法写错了，你比别人少了大概10个电子，所以能量差了
那么远。。。
赞一下(11人) 回复此楼
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭，其实这是修养不足学识浅薄的一种体现，也是可恶的偏执和愚蠢的自以
4楼2009-11-11 13:10:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaacjj

铜虫 (初入文坛)
引用回帖:
Originally posted by zhangmt at 2009-11-11 13:10:
很明显，你的赝势写法写错了，你比别人少了大概10个电子，所以能量差了
那么远。。。

少了10个电子？晕哦，不是吧~~~
我的赝势写法是按照自定义基组来写的，请这位虫友具体说明下好吗？该怎么改才行

[ Last edited by aaacjj on 2009-11-11 at 13:19 ]
赞一下 回复此楼
5楼2009-11-11 13:17:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+1,VIP+0):感谢交流 11-11 19:29
我猜的哦。不然不会差那么多啊——以后尽量贴出来你怎么写的，明眼人一看就知道了。应该是：
Hg 0
lanl2dz
****

Hg 0
lanl2dz

我猜你没有写后面第二遍的定义吧。所以就会少了10个内核层电子的。
赞一下(11人) 回复此楼
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭，其实这是修养不足学识浅薄的一种体现，也是可恶的偏执和愚蠢的自以
6楼2009-11-11 14:03:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaacjj

铜虫 (初入文坛)
引用回帖:
Originally posted by aaacjj at 2009-11-11 13:17:

少了10个电子？晕哦，不是吧~~~
我的赝势写法是按照自定义基组来写的，请这位虫友具体说明下好吗？该怎么改才行

[ Last edited by aaacjj on 2009-11-11 at 13:19 ]

贴上说明部分哦，

%chk=c16h10-cl-hg.chk
%mem=20MW
%nproc=1
# opt freq=noraman b3lyp/geom=connectivity gen pseudo=read
Title Card Required
0 2
C
C                  1              B1
C                  2              B2    1              A1
C                  3              B3    2              A2    1              D1
C                  4              B4    3              A3    2              D2
C                  1              B5    2              A4    3              D3
C                  3              B6    2              A5    1              D4
C                  4              B7    3              A6    2              D5
C                  8              B8    4              A7    3              D6
C                  7              B9    3              A8    2              D7
C                  9             B10    8              A9    4              D8
C                 11             B11    9             A10    8              D9
C                 12             B12   11             A11    9             D10
C                 13             B13   12             A12   11             D11
C                 14             B14   13             A13   12             D12
C                 15             B15   14             A14   13             D13
H                 16             B16   15             A15   14             D14
H                 15             B17   14             A16   13             D15
H                  7             B18    3             A17    2             D16
H                  1             B19    2             A18    3             D17
H                  2             B20    1             A19    6             D18
H                  6             B21    1             A20    2             D19
H                 10             B22    7             A21    3             D20
H                 13             B23   12             A22   11             D21
Cl                12             B24   11             A23    9             D22
Hg                11             B25    9             A24    8             D23

   B1             1.39284800
   B2             1.40084111
   B3             1.41912588
   B4             1.41916471
   B5             1.39291258
   B6             1.43967017
   B7             1.43332053
   B8             1.41911862
   B9             1.35659936
   B10            1.40085481
   B11            1.39284152
   B12            1.39291500
   B13            1.40078014
   B14            1.43970354
   B15            1.35658151
   B16            1.10092562
   B17            1.10092811
   B18            1.10090311
   B19            1.10046946
   B20            1.10023090
   B21            1.10023190
   B22            1.10090377
   B23            1.10025826
   B24            1.76000000
   B25            2.59975995
   A1           120.38056419
   A2           119.26682426
   A3           119.97568755
   A4           120.72891530
   A5           122.05557592
   A6           120.01018893
   A7           120.01361062
   A8           121.31109807
   A9           119.26868222
   A10          120.37817715
   A11          120.72955231
   A12          120.38027165
   A13          122.05711711
   A14          121.31297141
   A15          121.09630657
   A16          117.59523536
   A17          117.58920632
   A18          119.63893281
   A19          120.05639298
   A20          120.06065827
   A21          121.09061840
   A22          120.05105654
   A23          119.63788467
   A24          138.95642064
   D1             0.00000000
   D2             0.00000000
   D3             0.00000000
   D4           179.99694778
   D5           179.98991667
   D6             0.00670983
   D7          -179.99555390
   D8           179.99727543
   D9             0.00000000
   D10            0.00000000
   D11           -0.00394707
   D12         -179.98919216
   D13         -180.00000000
   D14         -179.97796207
   D15           -0.00437896
   D16            0.02035327
   D17          179.98424563
   D18         -179.99266950
   D19          179.97711317
   D20         -179.97954966
   D21         -179.99387444
   D22         -179.99116621
   D23         -179.98730937

1 2 1.5 6 1.5 20 1.0
2 3 1.5 21 1.0
3 4 1.5 7 1.5
4 5 1.5 8 1.5
5 6 1.5 16 1.5
6 22 1.0
7 10 2.0 19 1.0
8 9 1.5 14 1.5
9 10 1.5 11 1.5
10 23 1.0
11 12 1.5
12 13 1.5 25 1.0
13 14 1.5 24 1.0
14 15 1.5
15 16 2.0 18 1.0
16 17 1.0
17
18
19
20
21
22
23
24
25
26

C H Cl 0
6-31G*
****
Hg 0
lanl2dz
****

Hg 0
lanl2dz
请各位虫友看看啊
赞一下 回复此楼
7楼2009-11-11 15:15:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfan2192

金虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+1,VIP+0): 11-11 19:11
引用回帖:
Originally posted by aaacjj at 2009-11-11 12:49:

绝对能量指的就是高斯最后直接算出来那个数值咯？相对能量这个不太懂啊，能不能详细的说下或者有什么资料提供下，小虫在此谢过

比如说AB的能量为EAB
A原子的能量为EA  B原子为Eb
平均结合能就为EAB-EA-EB
其他的依次类推
赞一下(11人) 回复此楼
相逢是缘！
8楼2009-11-11 16:24:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaacjj

铜虫 (初入文坛)
引用回帖:
Originally posted by zhangfan2192 at 2009-11-11 16:24:

比如说AB的能量为EAB
A原子的能量为EA  B原子为Eb
平均结合能就为EAB-EA-EB
其他的依次类推

谢谢这位虫友的耐心解答，其实我现在要算的就是你所说的这个公式，但由于计算出来的结果偏差比较大，导致结合能的数值偏差也很大，所以就卡住了
赞一下 回复此楼
9楼2009-11-11 17:01:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
你把单个原子能量相加，总能量都要低于-1000 au，说明你用的基组和别人不一样。例如，如果Cl改用lanl2dz，总能量就会高~550 au.
赞一下(1人) 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
10楼2009-11-12 00:00:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 aaacjj 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭