24小时热门版块排行榜

返回列表

yaoyao-4921

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 328.5
帖子: 7
在线: 32.2小时
虫号: 2387285
注册: 2013-03-29
专业: 凝聚态物性 II ：电子结构

[求助] 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
不知道我哪里错了
不知道我说清楚没

回复此楼