24小时热门版块排行榜    

查看: 860  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-10-01 21:59 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 CfXuS1rDhLYN 2026-09-28 4/200 2026-10-01 09:23 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:45 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
请填处理意见