24小时热门版块排行榜    

查看: 431  |  回复: 1

yaoyao-4921

铜虫 (初入文坛)

[求助] Gaussian中内层分子的能级对上 已有1人参与

问题1:我计算了单独的LiClO4分子和单独的Cl原子的结构优化(原子也什么又不优化的)
应用全电子基组
LiClO4输入
# opt freq  pop=nbo B3LYP/6-311++G(3df,2pd)
nosymmetry

LiClO4

0 1
Li      2.548799992    0.948199987    0.899699986
Cl     -0.139500007   -0.176799998   -0.169400007
O       1.316900015   -0.750299990   -0.712700009
O      -0.084899999    1.297199965   -0.076399997
O      -0.426099986   -0.744099975    1.164299965
O      -1.206699967   -0.571799994   -1.112800002

Cl原子的输入
# opt freq  pop=nbo B3LYP/6-311++G(3df,2pd)
nosymmetry

Cl

0 2
Cl      0.0000000      0.0000000      0.0000000

LiClO4的分子能级
Alpha  occ. eigenvalues -- -101.90206 -19.15834 -19.15834 -19.15817 -19.15816
Alpha  occ. eigenvalues --   -9.76254  -7.52639  -7.52602  -7.52544  -2.04425
Alpha  occ. eigenvalues --   -1.30431  -1.07912  -1.07210  -1.04524  -0.66833
Alpha  occ. eigenvalues --   -0.57994  -0.56035  -0.55138  -0.42707  -0.40606
Alpha  occ. eigenvalues --   -0.40200  -0.39717  -0.38557  -0.32597  -0.31192
Alpha  occ. eigenvalues --   -0.31058
Alpha virt. eigenvalues --   -0.05979  -0.00461  -0.00103   0.00317   0.00719
Cl的原子轨道
Alpha  occ. eigenvalues -- -101.60520  -9.51563  -7.29071  -7.26812  -7.26812
Alpha  occ. eigenvalues --   -0.85118  -0.41373  -0.35896  -0.35896
Alpha virt. eigenvalues --    0.02441   0.06610   0.07110   0.07110   0.22762
  Beta  occ. eigenvalues -- -101.59839  -9.50902  -7.26852  -7.26419  -7.26419
  Beta  occ. eigenvalues --   -0.79767  -0.34390  -0.34390
  Beta virt. eigenvalues --   -0.24949   0.02779   0.07061   0.07062   0.07327
我相说对应的-101点多的都是Cl原子1s的能级,可是Cl原子和LiClO4中为什么有这么大的区别呢?能级的单位是Hartree的单位吧,就是说有0.3*27.2116=8.1eV的能量差
问题2:我算了一个孤立H原子的能级# opt freq  pop=nbo b3lyp/6-311+G(d)
Alpha  occ. eigenvalues --   -0.32170
Alpha virt. eigenvalues --    0.24504   2.14455
  Beta virt. eigenvalues --    0.02946   0.46427   2.56407
我想说,根据量子力学的计算,H原子的1s能级不是-13.6eV吗?这里明显不是啊,而且原子的能级和这个能级上有没有电子占据没有什么关系吧,就像解薛定谔方程那样,可是H的1s(Alpha)-0.3217     1s(beta)0.02946   差别这么大。。。
问题3:我计算水分子的溶剂化能
# opt freq B3lyp/6-31+g** scrf(read,pcm)
nosymmetry

water

0 1
O        0.3641    0.8512    0.0044
H        1.3027    0.5862    0.0046
H        0.1254    1.7967    0.0075

RADII=UAHF
EPS=78.3
DENSITY=1.0
scfvac
用Gaussian03算的,输出中/Delta
<psi(f)|H+V(f)/2|psi(f)>                     (a.u.) =     -76.446321
Total free energy in solution:
  with all non electrostatic terms            (a.u.) =     -76.616426
--------------------------------------------------------------------
(Unpolarized solute)-Solvent             (kcal/mol) =      -6.97
(Polarized solute)-Solvent               (kcal/mol) =      -8.65
Solute polarization                      (kcal/mol) =       0.88
Total electrostatic                      (kcal/mol) =      -7.76
--------------------------------------------------------------------
Cavitation energy                        (kcal/mol) =       4.45
Dispersion energy                        (kcal/mol) =    -154.93
Repulsion energy                         (kcal/mol) =      43.74
Total non electrostatic                  (kcal/mol) =    -106.74
DeltaG (solv)                            (kcal/mol) =    -114.51
我的是-114.51Kcal/mol         文献中是-6.3Kcal/mol
不知道我哪里错了
不知道我说清楚没
回复此楼

» 猜你喜欢

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

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

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

gauss98

禁虫 (正式写手)

★ ★
感谢参与,应助指数 +1
yaoyao-4921: 金币+1, DFT不是更准确吗,是不是DFT的能级不能比较,只能在同一个计算中比较不同的能级 2014-09-13 12:43:39
gmy1990: 金币+1 2014-09-13 18:16:40
本帖内容被屏蔽

2楼2014-09-13 08:30:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yaoyao-4921 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 Ocx3OKOzGU4q 2026-08-03 7/350 2026-08-04 13:22 by jzdxtiUmJ6PL
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 8/400 2026-08-04 13:15 by jzdxtiUmJ6PL
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 jRl3mE6ddZGq 2026-08-03 11/550 2026-08-04 12:58 by jzdxtiUmJ6PL
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 纯娱乐,不喜欢勿喷 +5 Tide man 2026-08-04 8/400 2026-08-04 11:43 by everxc
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 6/300 2026-08-04 11:34 by MWRYsWbJw4DE
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 7/350 2026-08-04 11:31 by MWRYsWbJw4DE
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 4/200 2026-08-04 11:26 by MWRYsWbJw4DE
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
信息提示
请填处理意见