24小时热门版块排行榜    

查看: 2547  |  回复: 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的回帖

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的回帖

yytsnake

至尊木虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+1,VIP+0):谢谢回帖。又及:“化学对量化计算有不良影响”啥意思?呵呵 10-16 12:50
溶剂模型的算法肯定与不同量化算法有较大差别,以更加近似的方法计量溶剂与体系的作用。版主的设计才体现了溶剂分子或离子间的相对真实作用。虽然化学对量化计算有不良影响,但作为一个基础还是得先加加深。这样在使用量化讨论化学问题时,才能更加体会其奥秘之处。
5楼2009-10-15 16:28:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

布丁3835

金虫 (小有名气)

自己顶一下,期待大家的帮助
6楼2009-10-16 09:32:19
已阅   回复此楼   关注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的回帖

zhengpc3505

金虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by yjcmwgk at 2009-10-15 16:19:09:
说实话我不喜欢scrf关键字
我喜欢在溶质分子周围摆上一圈儿溶剂分子

这样弄,你也不知道要放多少溶剂分子啊,而且还不容易收敛嘛
8楼2011-03-15 10:16:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by zhengpc3505 at 2011-03-15 10:16:39:
这样弄,你也不知道要放多少溶剂分子啊,而且还不容易收敛嘛

哇 前年的帖子都被你挖出来了~~~~
9楼2011-03-15 16:28:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhengpc3505

金虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by yjcmwgk at 2011-03-15 16:28:24:
哇 前年的帖子都被你挖出来了~~~~

呵呵,版主能说说怎么做到的吗?很关心,望解答啊
10楼2011-03-15 23:19:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 布丁3835 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +6 Ncdx123456 2026-03-13 7/350 2026-03-16 07:44 by L135790
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 求调剂 +3 清风问长安 2026-03-09 3/150 2026-03-14 02:15 by JourneyLucky
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 云南财经大学信息学院计算机学硕专硕学位点 +3 zjptai 2026-03-10 5/250 2026-03-14 01:23 by 飞行琦
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 材料专硕274一志愿陕西师范大学求调剂 +4 薛云鹏 2026-03-13 4/200 2026-03-13 10:40 by 学员8dgXkO
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[考研] 298求调剂 +3 Vv呀! 2026-03-10 3/150 2026-03-10 22:40 by 剑诗杜康
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
信息提示
请填处理意见