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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gamess/Q-Chem » GAMESS使用B3LYP+PCM尝试计算常温下水的离解常数遇到困难，求大神帮忙解惑，谢谢！
101/1返回列表
查看: 3610  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

张宇奇

金虫 (小有名气)

[求助] GAMESS使用B3LYP+PCM尝试计算常温下水的离解常数遇到困难，求大神帮忙解惑，谢谢！ 已有1人参与

最近想尝试计算溶液酸度对小分子酸分子离子浓度的影响，所以先尝试计算了水的离子积。计算过程是这样的。
1，选用密度泛函B3LYP在考虑容积效应下对H+、OH-和H2O进行构型优化，并通过热化学计算了吉布斯自由能。
数据分别为（单位都已经化为kj/mol）：H+，电子总能量-459.2206579，吉布斯自由能-26.254；OH，电子总能量-194791.2415，吉布斯自由能-20.172；H2O，电子总能量-196081.9112，吉布斯自由能7.319。
2，然后按照反应H2O=OH+H计算反应的吉布斯自由能变（其中其实有一个困惑，GAMESS中的吉布斯自由能是否包含了零点能？也希望得到解答。我默认它计算包含了零点能，就算在没有考虑零点能的情况下，反应能变仅差23kJ/mol）
3，反应的吉布斯自由能变为：deltaG=（E（H+）+G（H+））+（E（OH-）+G(OH-)）-（E（H2O）+G(H2O)）=777.704
    如果不考虑吉布斯自由能反应能量变化为：deltG=E（H+）+E（OH-）-E（H2O）=831.449
4，根据标准吉布斯自由能变deltaG=-RTln（【H+】*【OH-】/【H2O】）=-RTln（Kw/【H2O】）其中【】表示浓度，Kw表示水的离子积常数10^（-14）
5，代入数据可以得到Kw=exp【-deltaG/RT】*【H2O】=exp（-777.70*1000/（298.15*8.314））*55.6=3.093e-135.水的摩尔浓度以纯水为例55.6
由上边可以看到差别之大难以想象，如果给吉布斯机油能变式子里Kw取负对数，那么式子可以变成：
deltaG=RT*p【Kw】*ln（10/【H2O】），把ln项假定为C，也就是p【kw】=C*deltaG/RT，也就是p【Kw】和吉布斯自由能变成正比。计算所得pKw约为135，实际则为14，也就是吉布斯自由能变计算大了十倍！
由计算可以明显的发现电子能量的变化起到了反应能变的主要作用。反应能变的误差至少在6-700+KJ/mol，也就是大于150Kcal/mol，平均到每个粒子上误差在50Kcal/mol，我记得密度泛函误差不应该这么大的。这个误差有点大的离谱啊。就算气相和液相的差别在考虑了容积效应后也不应该差这么大吧？所以我认为我的计算中可能存在很多问题，但因为好长时间没有接触热力学和相关的计算，所以实在发现不了问题所在，还请大神指点一二。是我计算过程哪里疏忽了什么还是方法可能有点小问题？
另外我打算用CCSD（T）校正电子能量，是应该直接对溶解化效应计算中的除过溶剂化自由能的电子能量进行校正还是单独计算不考虑溶剂化效应的进行校正？因为GAMESS中CCSD（T）好像不支持PCM模块，所以这个也很困扰我！
总结一下，主要三个问题：1，GAMESS热力学计算中吉布斯自由能是否包含了零点能，计算反应吉布斯自由能应该怎样操作？2，我计算水的离子积常数时误差大的离谱，请问我的问题出在哪里？3，如果用CCSD（T）进行校正，应该怎样校正？是对考虑溶剂化效应后计算所得的电子能量还是未考虑溶剂化的能量进行校正。
以上这些问题困扰了我好些天，因为看了GAMESS手册和其他一些文献还是一头雾水，希望有朋友可以帮忙解答，万分感谢。我今天有事，收到解答后明天回复，再次感谢！祝学业有成、工作顺利！
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
在知识中充实；在奋发中图强
1楼 2017-09-25 20:42:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cxy3311


发自小木虫Android客户端
2楼2017-09-25 21:09:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mhhhwcf059

新虫 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
月只蓝: 金币-200, 屏蔽内容, 违规存档, 多次在技术性版块恶意灌水。 2017-10-20 08:48:16
本帖内容被屏蔽
3楼2017-09-25 23:51:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张宇奇

金虫 (小有名气)
引用回帖:
2楼: Originally posted by cxy3311 at 2017-09-25 21:09:42

请问有可以提示的地方么？谢谢哈！
赞一下 回复此楼
在知识中充实；在奋发中图强
4楼2017-09-26 08:45:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张宇奇

金虫 (小有名气)
引用回帖:
3楼: Originally posted by mhhhwcf059 at 2017-09-25 23:51:33
观看中

嘿嘿，观看怎么样了？有点提示么？谢谢！
赞一下 回复此楼
在知识中充实；在奋发中图强
5楼2017-09-26 08:54:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
张宇奇: 金币+20 2017-11-01 11:44:45
氢键你完全没考虑啊！
即使加pcm也是未包含氢键作用的，必须加显式溶剂。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

张宇奇
6楼2017-10-09 19:17:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张宇奇

金虫 (小有名气)
送红花一朵
引用回帖:
6楼: Originally posted by Jasminer at 2017-10-09 19:17:55
氢键你完全没考虑啊！
即使加pcm也是未包含氢键作用的，必须加显式溶剂。

您好，非常感谢您的回答，氢键确实我也是很困惑的，我只是想当然的认为因为电离前后氢键数量应该是相当的，所以就认为误差不至于像我计算的那么大。但我有点小白，对显式溶剂这个概念不了解，是要把溶剂直接加到计算粒子上么？可以麻烦您简单的介绍下或提供一篇文献介绍么？顺便可以告诉我在GAMESS中怎么实现么？ 不胜感激！
赞一下 回复此楼
在知识中充实；在奋发中图强
7楼2017-10-16 09:50:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)
引用回帖:
7楼: Originally posted by 张宇奇 at 2017-10-16 09:50:03
您好，非常感谢您的回答，氢键确实我也是很困惑的，我只是想当然的认为因为电离前后氢键数量应该是相当的，所以就认为误差不至于像我计算的那么大。但我有点小白，对显式溶剂这个概念不了解，是要把溶剂直接加到计 ...

显式溶剂就是优化前建模时直接在分子或离子周围摆一圈溶剂分子，后续计算过程无变化。
摆多少，怎么摆，一是靠化学直觉，二是靠收敛测试（三是靠信仰，四是靠人品 :-p）
赞一下 回复此楼
8楼2017-10-18 16:40:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张宇奇

金虫 (小有名气)
引用回帖:
8楼: Originally posted by Jasminer at 2017-10-18 16:40:51
显式溶剂就是优化前建模时直接在分子或离子周围摆一圈溶剂分子，后续计算过程无变化。
摆多少，怎么摆，一是靠化学直觉，二是靠收敛测试（三是靠信仰，四是靠人品 :-p）...

嗯嗯，谢谢啦，那样加水的话应该差不多，但因为构型太多，那样加的话工作量太大了，不过想想办法可以尝试下。总之非常感谢，祝工作顺利！
赞一下 回复此楼
在知识中充实；在奋发中图强
9楼2017-11-01 11:52:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蒲公英花
的

[ 发自手机版 http://muchong.com/3g ]
10楼2018-07-16 19:08:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 张宇奇 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 找导师 +4 yuanjiabao 2026-08-29 5/250 2026-08-30 09:52 by onesupeng
[基金申请] 国社科又开始会评了，不知道这次命运如何 +3 雨打竹帘 2026-08-30 3/150 2026-08-30 09:43 by KNIGHT2011
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 07:05 by ZPa0EcMwuECS
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 06:55 by ZPa0EcMwuECS
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:43 by ZPa0EcMwuECS
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 06:13 by ZPa0EcMwuECS
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 05:40 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:41 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗？ +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 基金不中，共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭