24小时热门版块排行榜

返回列表

vigaryang

木虫 (正式写手)

应助: 16 (小学生)
金币: 3466.3
散金: 677
红花: 6
帖子: 471
在线: 495.7小时
虫号: 227097
注册: 2006-03-23
性别: GG
专业: 应用化学

[求助] 溶液中的相互作用能量

我想计算溶质A和溶质B在溶剂环境下的相互作用能量，所用方法是用溶剂化（IEFPCM方法）条件下得到的复合物AB的能量减去溶剂化条件下单体A和单体B各自的能量。我在计算时所用的AB、A、B各自的“能量”数据来源于用GaussView打开优化输出文件时summary选项中提供的E(RB+HF-LYP)数据。
令人郁闷的是，通过上述方法计算得到的相互作用能居然是个大于0的正数。而以前在气相中计算得到的是负数，几十个KJ/mol的样子。

我想请教一下，上述结果的错误缘由是什么。是不是不该采用GaussView提供的数据？采用IEFPCM模型得到的log文件中给出了好几个能量值，不知道该选哪一个：

SCF Done:  E(RB+HF-LYP) =  -2995.19728131（也即GaussView显示的数据）    A.U. after 8 cycles
         Convg  = 0.9331D-08          -V/T =  2.0068
         S**2 = 0.0000
--------------------------------------------------------------------
Variational PCM results
=======================
                  (a.u.) = -2995.136943
                  (a.u.) = -2995.197281
Total free energy in solution:
  with all non electrostatic terms          (a.u.) = -2995.176521
--------------------------------------------------------------------
(Polarized solute)-Solvent             (kcal/mol) =    -37.86
--------------------------------------------------------------------
Cavitation energy                      (kcal/mol) =    25.59
Dispersion energy                      (kcal/mol) =    -13.37
Repulsion energy                      (kcal/mol) =    0.81
Total non electrostatic                (kcal/mol) =    13.03

回复此楼