24小时热门版块排行榜    

查看: 709  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 11:51 by xiangy672
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见