版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

布丁3835

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 656.9
散金: 169
帖子: 168
在线: 79.9小时
虫号: 565873
注册: 2008-05-29
性别: MM
专业: 药物资源

[交流] 【求助】关于用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的溶剂化能？谢谢高手的指点

回复此楼

» 猜你喜欢

北京师范大学胡德平课题组以及袁萌伟课题组考研招生已经有0人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有296人回复
祈福论文小修已经有21人回复
香港理工大学计算材料与物理课题组招聘博后已经有15人回复
香港理工大学计算材料与物理课题组招聘博后已经有15人回复
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生已经有0人回复
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？已经有2人回复
英国利物浦大学屠昕和雍雪课题组招收CSC博士生已经有26人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）已经有0人回复
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人已经有44人回复

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

【求助】溶剂化效应自由能的输出问题已经有10人回复
用g03计算溶剂化效应时,出现了l502的错误,请高手指教已经有4人回复
如何在计算机群上安装G09 已经有9人回复
G09中计算AIM计算过渡金属原子和氢原子间作用的问题已经有7人回复
G09计算问题已经有9人回复
请问能否在溶剂中只算过渡态单点能的情况下把虚频消掉？已经有12人回复
G09集群计算结果的chk文件转换问题已经有7人回复
【求助】请教用高斯09算溶剂化单点时，怎样得到溶剂中的Gibbs自由能？已经有9人回复
【求助】按照G09说明书计算荧光出错。出在第二步上。已经有13人回复
【求助】G09 没有优化，直接算频率时出错已经有9人回复
【求助】同一个结构，G09和G03 opt计算的能量差别大吗？已经有7人回复
【求助】溶剂化模型中关于溶剂的问题？已经有5人回复
【求助】关于g09和g03的formchk和unfchk问题已经有11人回复

1楼2009-10-15 11:08:18

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

yytsnake

至尊木虫 (职业作家)

应助: 357 (硕士)
贵宾: 0.01
金币: 18204.3
散金: 5
红花: 36
帖子: 3189
在线: 444.6小时
虫号: 92066
注册: 2005-09-08
性别: GG

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

那算水在以水为溶剂的体系中结果如何？

赞一下(1人)

回复此楼

2楼2009-10-15 14:14:40

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

布丁3835

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 656.9
散金: 169
帖子: 168
在线: 79.9小时
虫号: 565873
注册: 2008-05-29
性别: MM
专业: 药物资源

算水的也差很多，我怀疑溶剂化能没找对，但是确实不知道该往哪里找了

赞一下

回复此楼

3楼2009-10-15 15:18:15

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

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

QC强帖: 38
应助: 165 (高中生)
贵宾: 16.999
金币: 553.8
散金: 39880
红花: 395
沙发: 2
帖子: 10406
在线: 1528.1小时
虫号: 448368
注册: 2007-11-01
专业: 理论和计算化学
管辖: 科研工具资源

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+1,VIP+0):谢谢 10-27 20:21

说实话我不喜欢scrf关键字
我喜欢在溶质分子周围摆上一圈儿溶剂分子

赞一下(11人)

回复此楼

4楼2009-10-15 16:19:09

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

yytsnake

至尊木虫 (职业作家)

应助: 357 (硕士)
贵宾: 0.01
金币: 18204.3
散金: 5
红花: 36
帖子: 3189
在线: 444.6小时
虫号: 92066
注册: 2005-09-08
性别: GG

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+1,VIP+0):谢谢回帖。又及：“化学对量化计算有不良影响”啥意思？呵呵 10-16 12:50

溶剂模型的算法肯定与不同量化算法有较大差别，以更加近似的方法计量溶剂与体系的作用。版主的设计才体现了溶剂分子或离子间的相对真实作用。虽然化学对量化计算有不良影响，但作为一个基础还是得先加加深。这样在使用量化讨论化学问题时，才能更加体会其奥秘之处。

赞一下(11人)

回复此楼

5楼2009-10-15 16:28:37

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

布丁3835

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 656.9
散金: 169
帖子: 168
在线: 79.9小时
虫号: 565873
注册: 2008-05-29
性别: MM
专业: 药物资源

自己顶一下，期待大家的帮助

赞一下

回复此楼

6楼2009-10-16 09:32:19

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

九层楼

木虫 (小有名气)

QC强帖: 2
应助: 0 (幼儿园)
金币: 3273.2
帖子: 260
在线: 204.5小时
虫号: 86618
注册: 2005-08-15
性别: GG
专业: 理论和计算化学

★ ★
小木虫(金币+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).

赞一下(11人)

回复此楼

7楼2009-11-24 05:37:08

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