24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2283  |  回复: 6

kaegi

金虫 (小有名气)

[求助] 高斯09的溶剂化能到底包含哪几项?

以前用高斯03的时候溶剂化能都是用SCFVAC命令可以直接得到溶剂化能
=======================
                     (a.u.) =     -40.516992
                     (a.u.) =     -40.517259
                     (a.u.) =     -40.517263
                     (a.u.) =     -40.516986
                     (a.u.) =     -40.517262
Total free energy in solution:
  with all non electrostatic terms            (a.u.) =     -40.515696
--------------------------------------------------------------------
(Unpolarized solute)-Solvent             (kcal/mol) =      -0.17
(Polarized solute)-Solvent               (kcal/mol) =      -0.17
Solute polarization                      (kcal/mol) =       0.00
Total electrostatic                      (kcal/mol) =      -0.17
--------------------------------------------------------------------
Cavitation energy                        (kcal/mol) =       5.62
Dispersion energy                        (kcal/mol) =      -6.15
Repulsion energy                         (kcal/mol) =       1.51
Total non electrostatic                  (kcal/mol) =       0.98
DeltaG (solv)                            (kcal/mol) =       0.81
例如上面这个用高斯03计算出来的溶剂化能取DeltaG一项。其中包括静电作用项(total electrostatic)和非静电作用项(Total non electrostatic)。
高斯09是推荐用smd溶剂化模型计算溶剂化能,但是输出显示的是
SCF Done:  E(RB3LYP) =  -40.5164428887     A.U. after    7 cycles
             Convg  =    0.7143D-09             -V/T =  2.0080
KE= 4.019441594569D+01 PE=-1.205021820324D+02 EE= 2.610394409625D+01
SMD-CDS (non-electrostatic) energy       (kcal/mol) =       0.54
似乎只包含了非静电相互作用能,这与高斯03相比岂不是少了一项能量?
看高斯手册上说L124可以进行state-specific and external iteration SCRF calculations。另一种办法是单独计算分子的气相能量和溶剂中的能量,然后相减,得到溶剂化能。

我的问题是:高斯09有没有什么办法可以通过只添加一个参数啊,命令啊就能计算得到溶剂化能?溶剂化能到底是否应该包含静电作用能,为什么高斯09推荐的命令里不包含这一项了?如果不能通过添加命令直接计算得到溶剂化能,那么如何修改IOP调用L124进行计算?

[ Last edited by kaegi on 2011-6-2 at 16:37 ]
回复此楼

» 猜你喜欢

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

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

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

kaegi

金虫 (小有名气)

再顶一个
2楼2011-06-03 10:11:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

invincibly

金虫 (小有名气)

不懂帮顶
3楼2011-06-03 10:56:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flowerge

禁虫 (小有名气)

kaegi(金币+10): 谢谢回答,在溶液中优化就没有涉及非静电相互作用吗?有没有相关参考资料可以学习一下? 2011-06-03 14:26:47
kaegi(金币+20): 2011-06-03 14:27:12
本帖内容被屏蔽

4楼2011-06-03 13:34:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flowerge

禁虫 (小有名气)

本帖内容被屏蔽

5楼2011-06-07 13:42:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaegi

金虫 (小有名气)

引用回帖:
Originally posted by flowerge at 2011-06-07 13:42:21:
在溶剂中优化,g03的收敛标准没有涉及非静电相互作用,只是把非静电相互作用列出来。但是g09的smd就把这些都整合在收敛标准里了,这是gaussian公司的官方回复。

谢谢回复!
6楼2011-06-08 06:02:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjliumanhua

新虫 (初入文坛)

我是分别计算气相和溶剂中的能量,可是输出文件有点看不懂:除了总能量以外还有Diel,不知道要不要计算
还有就是溶剂化能到底是怎么回事啊
解离能又是怎么计算的呢
希望高手帮忙解答一下啊
7楼2012-10-26 10:38:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kaegi 的主题更新
信息提示
请填处理意见