24小时热门版块排行榜    

查看: 624  |  回复: 2
当前主题已经存档。

penghcp

木虫 (小有名气)

可爱小猫咪

[交流] 【求助】计算溶剂效应(pcm模型)能量问题,高手请进

各位大虾,小弟在用pcm计算时发现输出文件少了一项能量,具体情况如下:
关键字:#P B3LYP/6-31G(d) scf=tight int=finegrid SCRF=(PCM,Read,Solvent=Water)

输出文件关于溶剂化能量如下


Variational PCM results
=======================
                     (a.u.) =    -155.029051
                     (a.u.) =    -155.039159
Total free energy in solution:
  with all non electrostatic terms            (a.u.) =    -155.030980
--------------------------------------------------------------------
(Polarized solute)-Solvent               (kcal/mol) =      -6.34
--------------------------------------------------------------------
Cavitation energy                        (kcal/mol) =      11.83
Dispersion energy                        (kcal/mol) =      -7.15
Repulsion energy                         (kcal/mol) =       0.45
Total non electrostatic                  (kcal/mol) =       5.13
--------------------------------------------------------------------

当我查看高斯使用指南时发现他例子中的能量输出比我多了一项(Total Electrostatic 静电项),他书中的例子
-

------------- VARIATIONAL PCM RESULTS -------------
------------------------------------------------------
(a.u.) = -100.024068
(a.u.) = -100.028947
(a.u.) = -100.029186
(a.u.) = -100.023819
(a.u.) = -100.029187
Total free energy in sol.
(with non electrost.terms) (a.u.) = -100.028994
------------------------------------------------------
(Unpol.Solute)-Solvent (kcal/mol) = -3.06 静电项(包括极性
(Polar.Solute)-Solvent (kcal/mol) = -3.37 和非极性作用)
Solute Polarization (kcal/mol) = 0.16
Total Electrostatic (kcal/mol) = -3.21
------------------------------------------------------
Cavitation energy (kcal/mol) = 3.20 非静电项(包括
Dispersion energy (kcal/mol) = -4.12 空腔能量、分散
Repulsion energy (kcal/mol) = 1.04 能和排斥能)
Total non electr. (kcal/mol) = 0.12
------------------------------------------------------

但是我又查了中文手册发现没这一项 输出如下
SCF Done: E(RHF) = -100.029187240 A.U. after 5 cycles
Convg = 0.4249D-05 -V/T = 2.0033
S**2 = 0.0000
--------------------------------------------------------------------
Variational PCM results
=======================
(a.u.) = -98.568013
(a.u.) = -98.573228
Total free energy in solution:
with all non electrostatic terms (a.u.) = -98.569083
--------------------------------------------------------------------
(Polarized solute)-Solvent (kcal/mol) = -3.27
--------------------------------------------------------------------
Cavitation energy (kcal/mol) = 5.34
Dispersion energy (kcal/mol) = -3.08
Repulsion energy (kcal/mol) = 0.34
Total non electrostatic (kcal/mol) = 2.60
--------------------------------------------------------------------

这是怎么回事啊,计算溶剂效应时应该包含静电项能量么,是不是我的关键字输入有误,还是其他,敬请高手解答
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

S07111072

铁杆木虫 (著名写手)

★ ★
penghcp(金币+1,VIP+0):这个我知道,但是问题还是没解决,问什么我的输出文件没有这一项 9-24 12:32
yjcmwgk(金币+1,VIP+0): 9-24 13:38
Variational PCM results              
=======================              
            
            
            
            
            
Total free energy in solution:      
with all non electrostatic terms   
-------------------------------------
(Unpolarized solute)-Solvent         
(Polarized solute)-Solvent           
Solute polarization                  
Total electrostatic                  
-------------------------------------
Cavitation energy                    
Dispersion energy                    
Repulsion energy                     
Total non electrostatic              
DeltaG (solv)                        
包含静电项能量,一般PCM计算输出包含上面的能量项
2楼2009-09-24 09:28:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

S07111072

铁杆木虫 (著名写手)


yjcmwgk(金币+1,VIP+0): 9-25 09:16
引用回帖:
Originally posted by penghcp at 2009-9-23 21:27:
各位大虾,小弟在用pcm计算时发现输出文件少了一项能量,具体情况如下:
关键字:#P B3LYP/6-31G(d) scf=tight int=finegrid SCRF=(PCM,Read,Solvent=Water)

输出文件关于溶剂化能量如下


Variational PC ...

在你输入文件最后加上

PCMDOC      
RADII=UAHF  
SCFVAC
3楼2009-09-25 09:15:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 penghcp 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见