24小时热门版块排行榜    

查看: 3126  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

布丁3835

金虫 (小有名气)

[交流] 【求助】关于用G09算溶剂化能的问题 已有9人参与

我用G09计算DMSO在DMSO中的溶剂化能,输入文件的一部分如下,
opt freq hf/6-31+g**  scrf(read,pcm,solvent=dimethylsulfoxide),这样计算出来的结果找不到DeltaG,也没看到相关能量,如G03中类似的能量形式
Variational C-PCM results
=========================
                     (a.u.) =    -551.544999
                     (a.u.) =    -551.551263
                     (a.u.) =    -551.551943
                     (a.u.) =    -551.544233
                     (a.u.) =    -551.551944
Total free energy in solution:
  with all non electrostatic terms            (a.u.) =    -551.547465
--------------------------------------------------------------------
(Unpolarized solute)-Solvent             (kcal/mol) =      -3.93
(Polarized solute)-Solvent               (kcal/mol) =      -4.84
Solute polarization                      (kcal/mol) =       0.48
Total electrostatic                      (kcal/mol) =      -4.36
--------------------------------------------------------------------
Cavitation energy                        (kcal/mol) =      12.89
Dispersion energy                        (kcal/mol) =     -10.80
Repulsion energy                         (kcal/mol) =       0.72
Total non electrostatic                  (kcal/mol) =       2.81
DeltaG (solv)                            (kcal/mol) =      -1.55
看了下g09的帮助文献,将输入改成opt freq hf/6-31+g**  scrf(read,SMD,solvent=dimethylsulfoxide),结果在输出文件中出现了
Keep R1 ints in memory in canonical form, NReq=22711753.
Error on total polarization charges =  0.02051
SCF Done:  E(RHF) =  -551.571142309     A.U. after    9 cycles
             Convg  =    0.8499D-08             -V/T =  2.0010
SMD-CDS (non-electrostatic) energy             (kcal/mol) =      -0.02
(included in total energy above)
我想请问一下SMD.......=-0.02是不是g09算出来的溶剂化能,如果是的话,与G03的结果差得太大,如果不是的话,那该怎么写输入文件才能得到g09的溶剂化能?谢谢高手的指点
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

九层楼

木虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+1,VIP+0): 11-24 08:59
摘自G09 keyword scrf:

EXAMPLES
PCM Energy. In general, energy output from the default SCRF method appears in the normal way within the output file. For example, here are the sections of the output file containing the predicted energy from a Hartree-Fock and from an MP2 PCM calculation:

Hartree-Fock SCRF calculation:
SCF Done:  E(RHF) =  -99.4687828290     A.U. after    8 cycles
             Convg  =    0.2586D-08             -V/T =  2.0015

MP2 SCRF calculation:                       
E2 =    -0.1192799427D+00  EUMP2 =    -0.99584491345297D+02
The predicted energy in solution includes all computed corrections (unlike in Gaussian 03 output).
7楼2009-11-24 05:37:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

yytsnake

至尊木虫 (职业作家)

★
小木虫(金币+0.5):给个红包,谢谢回帖交流
那算水在以水为溶剂的体系中结果如何?
2楼2009-10-15 14:14:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

布丁3835

金虫 (小有名气)

算水的也差很多,我怀疑溶剂化能没找对,但是确实不知道该往哪里找了
3楼2009-10-15 15:18:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+1,VIP+0):谢谢 10-27 20:21
说实话我不喜欢scrf关键字
我喜欢在溶质分子周围摆上一圈儿溶剂分子
4楼2009-10-15 16:19:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[找工作] 售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一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 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文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 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文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 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 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
请填处理意见