| 查看: 2483 | 回复: 6 | ||
kaegi金虫 (小有名气)
|
[求助]
高斯09的溶剂化能到底包含哪几项？
|
|
以前用高斯03的时候溶剂化能都是用SCFVAC命令可以直接得到溶剂化能 ======================= 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 ] |
回复此楼» 猜你喜欢
苏州科技大学 物理学专业 接收调剂生
已经有0人回复
-
苏州科技大学 物理学专业 接收调剂生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有166人回复
-
苏州科技大学 物理学专业 接收调剂生
已经有0人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有0人回复
-
苏州科技大学 物理学专业 接收调剂
已经有0人回复
-
苏州科技大学 物理学专业 接收调剂生
已经有20人回复
-
河南农业大学化学学硕、材料与化工专硕二轮调剂还有少量名额
已经有0人回复
-
河南农业大学，化学学硕和材料化工专硕，二轮调剂
已经有0人回复
-
河南农业大学，化学学硕和材料化工专硕，二轮调剂
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【转载】窥探术 （美）萨姆·高斯林着
已经有34人回复
- 高斯计算微观动力学研究-计算结果：逆反应速率常数大于正反应速率常数 是否正常
已经有14人回复
- 请各位大侠看下高斯出错
已经有3人回复
- 高斯计算中zORA方法是否可以使用
已经有9人回复
- 【求助】请教用高斯09算溶剂化单点时，怎样得到溶剂中的Gibbs自由能？
已经有9人回复
- 高斯中键伸缩力常数与键角力常数的计算
已经有6人回复
- 高斯view中调节键长时怎么使四五个原子组成一个团（group）,和另一个原子之间进行调节
已经有4人回复
- 【求助】高斯计算出现Fatal Problem: The smallest alpha delta epsilon is -0.785373
已经有7人回复
- 【求助】运行高斯03or09的塔式服务器硬盘怎么配置合适呢！
已经有5人回复
- 【求助】高斯计算过渡态不收敛，请问如何解决？
已经有8人回复
- 【求助】高斯09 PCM指定pyridine溶剂的设定
已经有9人回复
kaegi
金虫 (小有名气)
- 应助: 9 (幼儿园)
- 金币: 1405.6
- 红花: 6
- 帖子: 276
- 在线: 87.4小时
- 虫号: 880820
- 注册: 2009-10-22
- 性别: GG
- 专业: 物理有机化学
2楼2011-06-03 10:11:26
invincibly
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 958.7
- 散金: 87
- 红花: 2
- 帖子: 145
- 在线: 82.8小时
- 虫号: 895734
- 注册: 2009-11-06
- 专业: 有机合成
3楼2011-06-03 10:56:12
kaegi(金币+10): 谢谢回答，在溶液中优化就没有涉及非静电相互作用吗？有没有相关参考资料可以学习一下？ 2011-06-03 14:26:47
kaegi(金币+20): 2011-06-03 14:27:12
kaegi(金币+20): 2011-06-03 14:27:12
|
本帖内容被屏蔽 |
4楼2011-06-03 13:34:40
|
本帖内容被屏蔽 |
5楼2011-06-07 13:42:21
kaegi
金虫 (小有名气)
- 应助: 9 (幼儿园)
- 金币: 1405.6
- 红花: 6
- 帖子: 276
- 在线: 87.4小时
- 虫号: 880820
- 注册: 2009-10-22
- 性别: GG
- 专业: 物理有机化学
6楼2011-06-08 06:02:37
7楼2012-10-26 10:38:54
