24小时热门版块排行榜    

查看: 2764  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 01:04 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 00:28 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 8/400 2026-08-22 23:51 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[教师之家] 跳槽后在研项目怎么办? +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见