已屏蔽 原因:{{ notice.reason }}已屏蔽
{{notice.noticeContent}}
~~空空如也
引用 last-order:
首先鼓励楼主发帖,但有一些存在的问题我必须讲一讲。(原文献没有看(其实我手头也有一篇计算pKa的综述,不过我也没有看))

我个人觉得使用DFT方法计算应该算“从头算”,半经验这种说法感觉不太合适。

本质上来讲,计算pKa就是要得到电离过程的准确的Gibbs自由能变。因此应当选择足够准确的理论水平进行计算。

楼主使用的泛函B3LYP应该没太大问题(不过计算有机物热力学量更推荐使用M06-2X),但是基组的使用很不合适。6-31+G这个基组没有极化函数,很不可取。通常而言,基组中的极化函数是不能缺少的,而弥散函数则在几何优化及频率计算时可以不加。为了得到准确的热力学量,还应当在高理论水平下计算单点能(你的文中提到没有算氢离子是因为该物种无电子,所以没有电子能,热力学校正部分只存在平动能\\(3/2 RT\\))。

所以我建议的计算的理论水平是,在M062X/defSVP水平下优化和频率分析(使用相应的校正因子),在B2PLYP-D3/def2TZVP水平下计算单点能(如果计算资源吃紧,那就用M062X/def2TZVP)。

关于溶剂模型,目前比较推荐使用SMD模型。而IEFPCM比Gaussian默认的PCM模型还差,尽量别用吧23333

其实我想说的是,很多文献使用的计算方法并不一定合理。除非是一定要重复文献,自己动手算的时候可以自己考虑一下有没有更合理的理论方法。

前面提到的那个综述是doi: 10.1002/wcms.1218,有兴趣可以看一看。
半经验法我也觉得不准确,文献中也没有这么说,但毕竟需要实验数据,所以我也拿不准,本来就像看看大家的意见
关于纯理论计算,文献中也说过“The computational determination of accurate pKa values is very demanding, as an error of only 1.36 kcal/mol in ΔGaq is quivalent to a unit difference in pKa value at room temperature.”

基组选择时加极化了,你可能没看见。关于弥散,因为有阴离子,所以还是加了。基组不大,分子也不大,所以在最后一两次优化时加上弥散什么的计算量也能承受。

溶剂模型同时选择IEFPCM和SMD是在Gaussian的手册中看到的,说是更好

很感谢如此认真的回答和建议
文号 / 838555

千古风流
名片发私信
学术分 4
总主题 39 帖总回复 571 楼拥有证书:学者 机友 笔友
注册于 2012-09-22 15:42最后登录 2023-08-25 22:23
主体类型:个人
所属领域:无
认证方式:手机号
IP归属地:未同步

个人简介

暂未填写
文件下载
加载中...
{{errorInfo}}
{{downloadWarning}}
你在 {{downloadTime}} 下载过当前文件。
文件名称:{{resource.defaultFile.name}}
下载次数:{{resource.hits}}
上传用户:{{uploader.username}}
所需积分:{{costScores}},{{holdScores}}下载当前附件免费{{description}}
积分不足,去充值
文件已丢失

当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
投诉或举报
加载中...
{{tip}}
请选择违规类型:
{{reason.type}}

空空如也

插入资源
全部
图片
视频
音频
附件
全部
未使用
已使用
正在上传
空空如也~
上传中..{{f.progress}}%
处理中..
上传失败,点击重试
等待中...
{{f.name}}
空空如也~
(视频){{r.oname}}
{{selectedResourcesId.indexOf(r.rid) + 1}}
处理中..
处理失败
插入表情
我的表情
共享表情
Emoji
上传
注意事项
最大尺寸100px,超过会被压缩。为保证效果,建议上传前自行处理。
建议上传自己DIY的表情,严禁上传侵权内容。
点击重试等待上传{{s.progress}}%处理中...已上传,正在处理中
空空如也~
处理中...
处理失败
加载中...
草稿箱
加载中...
此处只插入正文,如果要使用草稿中的其余内容,请点击继续创作。
{{fromNow(d.toc)}}
{{getDraftInfo(d)}}
标题:{{d.t}}
内容:{{d.c}}
继续创作
删除插入插入
插入公式
评论控制
加载中...
文号:{{pid}}
加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
收藏
取消收藏
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
管理提醒
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}
ID: {{user.uid}}