24小时热门版块排行榜    

查看: 317  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 民办高校入职后稳定吗? +15 905452934 2024-05-16 23/1150 2024-05-21 10:20 by balabiubiu
[教师之家] 另一半是老师会加分哦 +14 acogw9004 2024-05-14 14/700 2024-05-21 08:54 by bxwisdom_886
[教师之家] 博士去高校就是为了有寒暑假吗? +19 wenwen0825 2024-05-16 25/1250 2024-05-21 08:49 by xiaoyucau
[考博] 双非博士还是985科研助理,然后再读博 +6 lxdatj123 2024-05-18 13/650 2024-05-21 08:19 by lxdatj123
[教师之家] 白天不懂夜的黑,90后青椒与60后老板 +5 zylfront 2024-05-18 6/300 2024-05-21 07:32 by 战狼?
[硕博家园] 博三一直没文章怎么办 +24 133456 2024-05-17 38/1900 2024-05-21 07:26 by 133456
[硕博家园] 海外博士,国内博后找工作求建议 +8 905452934 2024-05-16 22/1100 2024-05-20 21:42 by littlezl
[教师之家] 谁要做个人网站,可以交流啊! +12 hjc404 2024-05-14 36/1800 2024-05-20 21:11 by 枯井中抽水
[考博] 海南大学国家级人才田新龙教授团队氢燃料电池与器件方向招收科研助理 (重点支持方向) +6 Dreammzp 2024-05-17 9/450 2024-05-20 18:40 by Dreammzp
[考博] 考博跨专业可以跨多远 +5 竹外。。。 2024-05-18 6/300 2024-05-19 21:51 by 摩天思瑞
[教师之家] 来用亲身体会,一起说说年轻老师的辛苦 +19 zylfront 2024-05-17 25/1250 2024-05-19 16:19 by otani
[论文投稿] 推荐转投( transfer pending)是否有用? 50+3 lily5289 2024-05-17 7/350 2024-05-19 15:11 by wanghuawei
[无机/物化] 请问什么溶剂能溶解二氧化锰 (金币+5) +4 这是春天 2024-05-15 4/200 2024-05-19 11:24 by wanghuawei
[硕博家园] 又忍不住想吐槽文科的学术圈 +6 nvizqxuf 2024-05-18 8/400 2024-05-19 10:04 by nvizqxuf
[硕博家园] 耐高温垫片求购 +6 Sexyflea 2024-05-16 9/450 2024-05-18 08:31 by nanbeiY
[论文投稿] 一个审稿人已返回,也没邀请审稿人,也没人审稿中,编辑要闹啥? 3+3 林师妹 2024-05-14 12/600 2024-05-17 22:33 by 林师妹
[论文投稿] 投稿成功后又想撤回 +5 otani 2024-05-16 5/250 2024-05-17 16:02 by topedit
[电化学] 锂离子电池石墨负极用 1M LiPF6 in DEC:EC=1:1 Vol% 可以吗? 50+3 fffhhhhh 2024-05-15 8/400 2024-05-17 14:57 by 多点关心多点i
[基金申请] 青基 +3 变成超人 2024-05-15 4/200 2024-05-17 12:42 by ssxclkj
[考博] 深度学习申博 +4 酸菜鱼. 2024-05-15 4/200 2024-05-16 11:12 by 安塔瓦拉多
信息提示
请填处理意见