24小时热门版块排行榜    

查看: 331  |  回复: 2
【悬赏金币】回答本帖问题,作者万事顺意北将赠送您 20 个金币

万事顺意北

新虫 (正式写手)

[求助] 求助已有1人参与

想问一下如何用gaussian计算溶质的总静电能,可以用gaussian计算吗

发自小木虫IOS客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhou2009: 金币+10 2022-08-25 13:37:47
溶解自由能分解成不同的组成部分可以用Dovacuum关键词。比如写成下面的关键词

#p m052x/6-31g(d) SCRF=(smd,externaliteration,dovacuum)

这种写法算出来的溶解自由能就是跟自己算两次单点能(气相和溶剂中)相减得到的溶解自由能数值是一样的。

输出文件重会有类似下面的部分

Self-consistent PCM results
===========================
<psi(0)|   H    |psi(0)>                     (a.u.) =    -569.434121
<psi(0)|H+V(0)/2|psi(0)>                     (a.u.) =    -569.533542
<psi(0)|H+V(f)/2|psi(0)>                     (a.u.) =    -569.539399
<psi(f)|   H    |psi(f)>                     (a.u.) =    -569.431013
<psi(f)|H+V(f)/2|psi(f)>                     (a.u.) =    -569.536292
Total free energy in solution:
- with all non electrostatic terms           (a.u.) =    -569.526641
--------------------------------------------------------------------
(Unpolarized solute)-Solvent             (kcal/mol) =     -62.39
(Polarized solute)-Solvent               (kcal/mol) =     -66.06
Solute polarization                      (kcal/mol) =       1.95
Total electrostatic                      (kcal/mol) =     -64.11
--------------------------------------------------------------------
SMD-CDS (non-electrostatic) energy       (kcal/mol) =       6.06
Total non electrostatic                  (kcal/mol) =       6.06
DeltaG (solv)                            (kcal/mol) =     -58.06

其中 Total electrostatic  就是溶解自由能中总的静电部分的贡献。
2楼2022-08-25 12:44:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

万事顺意北

新虫 (正式写手)

引用回帖:
2楼: Originally posted by paramecium86 at 2022-08-25 12:44:16
溶解自由能分解成不同的组成部分可以用Dovacuum关键词。比如写成下面的关键词

#p m052x/6-31g(d) SCRF=(smd,externaliteration,dovacuum)

这种写法算出来的溶解自由能就是跟自己算两次单点能(气相和溶剂中) ...

非常感谢

发自小木虫IOS客户端
3楼2022-08-26 08:33:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 万事顺意北 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 讨厌酒桌文化 +10 hahamyid 2024-05-31 10/500 2024-06-03 17:50 by beefly
[硕博家园] 又想换工作 +16 brightmj 2024-05-27 27/1350 2024-06-03 16:11 by 小小冲天虎
[硕博家园] 每到中夜,情难自抑 +37 sioc-sunj 2024-05-28 61/3050 2024-06-03 15:16 by EllieWU
[基金申请] B口人才项目 +11 WOWO159357 2024-05-29 23/1150 2024-06-03 11:18 by XUAAAAA
[论文投稿] 《控制理论与应用》期刊收版面费吗? +5 ygj2015 2024-05-29 5/250 2024-06-03 09:19 by lmpmsm
[电化学] 电化学传感器 +3 哈哈嘿哟 2024-06-01 5/250 2024-06-03 08:46 by mecart2021
[教师之家] 职能部门工作人员态度不好是普遍的吗?怎么让他们态度好一些? +7 河西夜郎 2024-06-01 7/350 2024-06-03 08:37 by YY392619
[硕博家园] 论大家对6070后普通教授导师的看法 +8 SNaiL1995 2024-05-28 14/700 2024-06-03 05:07 by SNaiL1995
[论文投稿] 求助大神,Fe和Al离子对MOF都有淬灭,当两种离子共存时,怎么区分两种离子? 10+3 maoxiao 2024-06-02 3/150 2024-06-03 01:43 by nono2009
[教师之家] 选择 +8 459582015 2024-05-31 11/550 2024-06-02 18:30 by auvauv
[硕博家园] 各位同学能否分享一下实验室的学生劳务发放标准呀? +19 ma3252788 2024-05-30 19/950 2024-06-02 18:15 by Lcy6666
[论文投稿] 编辑对文章提的意见 40+3 qa8386136 2024-06-01 3/150 2024-06-02 13:28 by bobvan
[考博] 导师不让硕转博,让我去国外读博,能理解吗? +12 萧山幽谷 2024-05-29 20/1000 2024-06-02 12:01 by yuan0806
[考博] 求25博导,金属增材制造方向 +3 22机械 2024-06-01 3/150 2024-06-02 11:17 by Napoleonsky
[教师之家] 中年 (金币+3) +18 459582015 2024-05-28 19/950 2024-06-01 00:41 by 沈婉婷.Girl
[文学芳草园] 物是人非 +4 myrtle 2024-05-30 4/200 2024-05-30 15:05 by mapenggao
[论文投稿] 审稿专家比较坚定的让补充实验,但实在没法补充实验,修回还有希望吗? (EPI+1) 3+3 qweasd12345 2024-05-29 6/300 2024-05-30 08:11 by qweasd12345
[论文投稿] 真急着毕业,CPB主编终审17天了,邮件催稿了两次,就是一点动静没有 5+3 kkkk夏 2024-05-28 6/300 2024-05-29 11:18 by hitsdu
[论文投稿] 高手朋友推荐比较容易投稿和录用的SCI期刊,不投稿收费SC,对分区没有要求 5+3 xintangren 2024-05-28 4/200 2024-05-29 10:46 by xintangren
[基金申请] E05青基有几个评审 +4 KYXY123 2024-05-28 4/200 2024-05-28 19:25 by popt2t
信息提示
请填处理意见