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

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gamess/Q-Chem » GAMESS使用B3LYP+PCM尝试计算常温下水的离解常数遇到困难，求大神帮忙解惑，谢谢！
【调剂】北京石油化工学院2024年16个专业接受调剂
51/1返回列表
查看: 2571  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者张宇奇将赠送您 80 个金币
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

张宇奇

金虫 (小有名气)

[求助] 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的回帖

张宇奇

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

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

cxy3311


发自小木虫Android客户端
2楼2017-09-25 21:09:42
 已阅   回复此楼   关注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的回帖
查看全部 10 个回答
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 食品工程考数学283还能有什么调剂选择 +7 宇宙无敌美少女 2024-04-12 8/400 2024-04-16 10:11 by cyclodextrin726
[考研] 352求调剂 +6 歪比歪比888 2024-04-10 6/300 2024-04-15 23:05 by 可怜下孩子吧
[考研] 求调剂 +4 桃岸雪 2024-04-15 5/250 2024-04-15 18:49 by mthwyj
[论文投稿] 可以打电话问编辑部是否可以先发录用通知吗 +5 双倍好运锦鲤 2024-04-14 8/400 2024-04-15 15:14 by topedit
[论文投稿] Journal of The Science of Food and Agriculture投稿遇到格式问题 5+3 tongtongco 2024-04-12 3/150 2024-04-15 09:19 by shuoyi
[考研] 300求调剂 +5 @Xqy668800 2024-04-14 6/300 2024-04-14 23:03 by lincunhui
[考研] 338求调剂 +3 18280338551 2024-04-14 5/250 2024-04-14 10:03 by tcni
[考研] 334求调剂 +4 学药救人 2024-04-13 6/300 2024-04-13 20:27 by 献世的王
[考研] 314求调剂 +18 阿洛法 2024-04-12 20/1000 2024-04-13 13:40 by kxxmaqq
[教师之家] 表弟博士毕业，选30万安家费还是副教授？？？ +10 yingyu427 2024-04-11 12/600 2024-04-13 11:16 by karl_messi
[考研] 310求调剂 +17 a1845052693 2024-04-09 21/1050 2024-04-13 09:44 by 最爱懒羊羊。
[考研] 277求调剂 +10 尊的很想上学 2024-04-12 11/550 2024-04-13 07:11 by coco1981
[考研] 266求调剂 +3 ～酥糖 2024-04-12 15/750 2024-04-12 21:22 by 赵燕高仙兰
[考研] 275求调剂 +6 南幕星辰 2024-04-12 6/300 2024-04-12 18:56 by qikanlunwen
[考研] 化工学硕276求调剂 +35 星星陨落 2024-04-10 41/2050 2024-04-12 08:41 by lus902
[考研] 考研调剂 +5 夏末的晨曦12 2024-04-10 6/300 2024-04-11 18:51 by 17814699553
[考研] 生物与医药调剂！ +9 夏天快过去了吗 2024-04-10 10/500 2024-04-11 17:29 by 王有田有学上
[考研] 329化学调剂 +8 瞧瞧乔儿 2024-04-10 9/450 2024-04-11 09:46 by yqWu武
[考研] 求调剂 +3 LED绿光灯 2024-04-10 6/300 2024-04-10 17:22 by LED绿光灯
[考研] 085400调剂 +4 Dyqqqq 2024-04-10 4/200 2024-04-10 11:20 by 伤痕迹1990
信息提示
关闭
请填处理意见
关闭