24小时热门版块排行榜    

查看: 605  |  回复: 2

huilaoshu999

金虫 (著名写手)

[交流] 【求助】请教关于溶质与溶剂分子间的静电作用能的问题? 已有2人参与

------------------------------------------------------------------------------
Polarizable Continuum Model (PCM)
=================================
Model                : PCM.
Atomic radii         : UA0 (Simple United Atom Topological Model).
Polarization charges : Total charges.
Charge compensation  : None.
Solution method      : Matrix inversion.
Cavity               : GePol (RMin=0.200 OFac=0.890).
                        Default sphere list used, NSphG=    3.
                        Tesserae with average area of 0.200 Ang**2.
Solvent              : Acetonitrile, Eps=  36.640000.
------------------------------------------------------------------------------
Using symmetry in molecular cavity generation.
One-electron integrals computed using PRISM.
NBasis=    63 RedAO= T  NBF=    46    17
NBsUse=    63 1.00D-06 NBFU=    46    17
Harris functional with IExCor=  205 diagonalized for initial guess.
ExpMin= 4.38D-02 ExpMax= 4.17D+03 ExpMxC= 6.27D+02 IAcc=2 IRadAn=         4 AccDes= 0.00D+00
HarFok:  IExCor= 205 AccDes= 0.00D+00 IRadAn=         4 IDoV=1
ScaDFX=  1.000000  1.000000  1.000000  1.000000
Initial guess orbital symmetries:
       Occupied  (A1) (A1) (A1) (A1) (A1) (A1) (E) (E) (A1) (E)
                 (E)
       Virtual   (E) (E) (A1) (A1) (E) (E) (A1) (E) (E) (A1) (E)
                 (E) (A1) (E) (E) (A1) (A1) (A1) (A1) (E) (E) (E)
                 (E) (A1) (A1) (E) (E) (A1) (E) (E) (A1) (A1) (E)
                 (E) (A1) (E) (E) (A1) (E) (E) (E) (E) (A1) (E)
                 (E) (E) (E) (A1) (A1) (A1) (A1) (A1)
The electronic state of the initial guess is 1-A1.
Requested convergence on RMS density matrix=1.00D-04 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-02.
Requested convergence on             energy=1.00D-02.
No special actions if energy rises.
Keep R1 integrals in memory in canonical form, NReq=     2767344.
Error on total polarization charges =  0.00935
SCF Done:  E(RHF) =  -131.940813624     A.U. after    5 cycles
             Convg  =    0.9612D-04             -V/T =  2.0018
             S**2   =   0.0000
--------------------------------------------------------------------
Variational PCM results
=======================
                     (a.u.) =    -131.929015
                     (a.u.) =    -131.940814
Total free energy in solution:
  with all non electrostatic terms            (a.u.) =    -131.932196
--------------------------------------------------------------------
(Polarized solute)-Solvent               (kcal/mol) =      -7.40
--------------------------------------------------------------------
Cavitation energy                        (kcal/mol) =       9.36
Dispersion energy                        (kcal/mol) =      -4.20
Repulsion energy                         (kcal/mol) =       0.25
Total non electrostatic                  (kcal/mol) =       5.41
--------------------------------------------------------------------
Partition over spheres:
Sphere  on Atom  Surface  Charge   GEl     GCav    GDR
    1      C1       2.32   0.004   -0.02    0.67   -0.04
    2      C2      64.67  -0.263   -2.53    6.01   -2.56
    3      N6      23.88   0.242   -4.71    2.68   -1.35
Added spheres:     7.51   0.008   -0.14    0.00    0.00
--------------------------------------------------------------------

**********************************************************************

            Population analysis using the SCF density.

**********************************************************************

Orbital symmetries:
       Occupied  (A1) (A1) (A1) (A1) (A1) (A1) (E) (E) (A1) (E)
                 (E)
       Virtual   (A1) (A1) (E) (E) (E) (E) (A1) (E) (E) (A1) (A1)
                 (E) (E) (A1) (E) (E) (A1) (A1) (A1) (E) (E) (E)
                 (E) (A1) (A1) (A1) (E) (E) (E) (E) (A1) (A1) (A1)
                 (E) (E) (E) (E) (A1) (E) (E) (E) (E) (A1) (E)
                 (E) (E) (E) (A1) (A1) (A1) (A1) (A1)
The electronic state is 1-A1.
Alpha  occ. eigenvalues --  -15.58977 -11.30022 -11.27522  -1.24896  -1.03740
Alpha  occ. eigenvalues --   -0.69516  -0.62497  -0.62497  -0.57522  -0.48287
Alpha  occ. eigenvalues --   -0.48287
Alpha virt. eigenvalues --    0.06679   0.09087   0.11881   0.11881   0.14621
Alpha virt. eigenvalues --    0.14621   0.15710   0.18266   0.18266   0.20642
Alpha virt. eigenvalues --    0.29575   0.29850   0.29850   0.30087   0.36654
Alpha virt. eigenvalues --    0.36654   0.40985   0.47406   0.69593   0.85120
Alpha virt. eigenvalues --    0.85120   0.91827   0.91827   0.99891   1.06998
Alpha virt. eigenvalues --    1.13497   1.16367   1.16367   1.24274   1.24274
Alpha virt. eigenvalues --    1.36343   1.43002   1.65732   1.68865   1.68865
Alpha virt. eigenvalues --    1.76572   1.76572   1.97844   2.20130   2.20130
Alpha virt. eigenvalues --    2.26970   2.26970   2.51998   2.58696   2.58696
Alpha virt. eigenvalues --    3.11613   3.11613   3.27249   3.67329   4.52141
Alpha virt. eigenvalues --    4.59286   5.05432
          Condensed to atoms (all electrons):
              1          2          3          4          5          6
     1  C    4.787052   0.138084  -0.050464  -0.050464  -0.050464   0.754286
     2  C    0.138084   5.396979   0.403581   0.403581   0.403581  -0.052259
     3  H   -0.050464   0.403581   0.422474  -0.020543  -0.020543  -0.002631
     4  H   -0.050464   0.403581  -0.020543   0.422474  -0.020543  -0.002631
     5  H   -0.050464   0.403581  -0.020543  -0.020543   0.422474  -0.002631
     6  N    0.754286  -0.052259  -0.002631  -0.002631  -0.002631   6.888662
Mulliken atomic charges:
              1
     1  C    0.471969
     2  C   -0.693548
     3  H    0.268125
     4  H    0.268125
     5  H    0.268125
     6  N   -0.582796
Sum of Mulliken charges=   0.00000
Atomic charges with hydrogens summed into heavy atoms:
              1
     1  C    0.471969
     2  C    0.110828
     3  H    0.000000
     4  H    0.000000
     5  H    0.000000
     6  N   -0.582796
Sum of Mulliken charges=   0.00000
Electronic spatial extent (au):  =   161.6983
Charge=     0.0000 electrons
Dipole moment (field-independent basis, Debye):
    X=     0.0000    Y=     0.0000    Z=    -5.1515  Tot=     5.1515
Quadrupole moment (field-independent basis, Debye-Ang):
   XX=   -17.7371   YY=   -17.7371   ZZ=   -21.9228
   XY=     0.0000   XZ=     0.0000   YZ=     0.0000
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
   XX=     1.3952   YY=     1.3952   ZZ=    -2.7905
   XY=     0.0000   XZ=     0.0000   YZ=     0.0000
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=     0.0000  YYY=     0.8857  ZZZ=   -17.6082  XYY=     0.0000
  XXY=    -0.8857  XXZ=    -2.8911  XZZ=     0.0000  YZZ=     0.0000
  YYZ=    -2.8911  XYZ=     0.0000
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=   -23.9495 YYYY=   -23.9495 ZZZZ=  -184.9841 XXXY=     0.0000
XXXZ=     0.0000 YYYX=     0.0000 YYYZ=    -1.4541 ZZZX=     0.0000
ZZZY=     0.0000 XXYY=    -7.9832 XXZZ=   -31.1307 YYZZ=   -31.1307
XXYZ=     1.4541 YYXZ=     0.0000 ZZXY=     0.0000
N-N= 5.886190824948D+01 E-N=-4.255290960902D+02  KE= 1.317065684090D+02
Symmetry A'   KE= 1.269810328033D+02
Symmetry A"   KE= 4.725535605678D+00
这个gaussian03计算结果文件中,静电作用力从哪里看?我只看到了空穴能、散射-排斥作用能和非静电作用能总和啊,溶剂化能无法计算啊 !
回复此楼

» 猜你喜欢

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

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

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

hzfish

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+1):谢谢交流! 2010-04-23 22:11
分子坐标下加关键字:scfvac
2楼2010-04-22 22:33:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+2):谢谢交流! 2010-04-23 22:11
空穴能、散射、排斥作用能三项加起来不就是那个非静电作用项嘛。
(Polarized solute)-Solvent那个就是静电项啊。
两个加起来就是△G(sol)啊。
也就是
Total free energy in solution with all non electrostatic terms (溶液中的能量) 这一项  减去
     (气相中的能量)   之差

[ Last edited by zhangmt on 2010-4-22 at 23:28 ]
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
3楼2010-04-22 23:26:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huilaoshu999 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见