计算化学公社

 找回密码 Forget password
 注册 Register
Views: 4259|回复 Reply: 10
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 关于多金属氧酸盐氧化还原电位的计算问题

[复制链接 Copy URL]

70

帖子

0

威望

564

eV
积分
634

Level 4 (黑子)

选用高斯 使用连续溶剂化模型CPCM 对多金属氧酸盐 进行氧化还原电位计算时 实验测定值 -0.225V  我的理论计算值 -0.618V
计算细节:溶剂化能(g+s) : 泛函 M06-2X  (g+s)  基组:B O 6-31G(d)   + W SDD
                气态自由能 :泛函 B3LYP 基组:B O  6-31G(d,p)+ W  LANL2DZ
                气态单点能: 泛函 B3LYP  基组:B O   6-311++(d,p)  + W  SDD
                根据热力学循环 计算其氧化还原电位  考虑了参考电极电位
               
① 不知道偏差如此大 是计算细节上出现了问题 还是 选用溶剂化模型上 出现了原则上的问题
② 看到文献中关于 QM/MM 不连续溶剂化模型 计算氧化还原电位的方法 不知是否可行 如果可行 通常选用什么软件进行计算                                                                  

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

2#
发表于 Post on 2022-4-12 23:01:17 | 只看该作者 Only view this author
可能是泛函误差。根据https://pubs.acs.org/doi/10.1021/acs.jctc.6b00252这篇文章,B3LYP算CV的误差就是在0.4eV量级。解决方法一个是换其他泛函,比如刚才那篇文章测了M06-2X还不错,但是测的是有机分子,未必适用于金属体系;另一个是按https://pubs.acs.org/doi/10.1021/jp300485t的方法,不以参考电极的绝对电位为参比,而是用一个含W或者和W同周期过渡金属的常见配合物做参比,误差能抵消很多。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

3#
发表于 Post on 2022-4-13 01:31:57 | 只看该作者 Only view this author
误差来自于两方面:气相自由能计算误差,溶解自由能计算误差

从气相自由能计算角度来说,当前单点级别太low,B3LYP本身算过渡金属配合物就不是什么好的选择,基组也偏小。具体什么级别更合适,体系结构图贴出来才好说。

溶解自由能方面,通常误差大部分来自于带电状态的溶解自由能的计算,不管什么隐式溶剂模型,误差都不小。可以考虑显式溶剂模型。QM/MM或者Gaussian可以做的电子嵌入的ONIOM(QM:MM)也可以考虑。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

70

帖子

0

威望

564

eV
积分
634

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2022-4-13 14:55:21 | 只看该作者 Only view this author
wzkchem5 发表于 2022-4-12 23:01
可能是泛函误差。根据https://pubs.acs.org/doi/10.1021/acs.jctc.6b00252这篇文章,B3LYP算CV的误差就是在 ...

好的 感谢文章的分享  我一定回去仔细阅读  

70

帖子

0

威望

564

eV
积分
634

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-4-13 14:58:01 | 只看该作者 Only view this author
sobereva 发表于 2022-4-13 01:31
误差来自于两方面:气相自由能计算误差,溶解自由能计算误差

从气相自由能计算角度来说,当前单点级别太 ...

好的 感谢sob老师解答  我再试着用其他基组和泛函进行优化   并尝试使用显式溶剂化模型的计算

5

帖子

0

威望

19

eV
积分
24

Level 1 能力者

6#
发表于 Post on 2022-11-30 21:51:40 | 只看该作者 Only view this author
请问,您现在会算溶剂化自由能了吗

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

7#
发表于 Post on 2022-11-30 22:11:13 | 只看该作者 Only view this author
binghua 发表于 2022-11-30 14:51
请问,您现在会算溶剂化自由能了吗

http://sobereva.com/327的方法计算
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

5

帖子

0

威望

19

eV
积分
24

Level 1 能力者

8#
发表于 Post on 2022-12-3 18:47:27 | 只看该作者 Only view this author
谢谢,想知道用MS可以算出来氧化还原电位吗

57

帖子

0

威望

341

eV
积分
398

Level 3 能力者

9#
发表于 Post on 2024-9-6 21:52:29 | 只看该作者 Only view this author
wzkchem5 发表于 2022-4-12 23:01
可能是泛函误差。根据https://pubs.acs.org/doi/10.1021/acs.jctc.6b00252这篇文章,B3LYP算CV的误差就是在 ...

请问老师,关于第二个方法,我可以理解为:先通过DFT分别计算待测铁络合物和二茂铁的绝对氧化还原电位,并相减得到计算氧化还原电位差值,再把这个差值加到二茂铁的实验氧化还原电位上得到待测铁络合物的相对氧化还原电位吗?

451

帖子

9

威望

6037

eV
积分
6668

Level 6 (一方通行)

BSJ Institute

10#
发表于 Post on 2024-9-6 22:22:19 | 只看该作者 Only view this author
带金属的电位随着泛函能差很多,而且与HF成分并没有严格的平行关系。关于溶剂化,我的经验:高电荷阴离子多酸包含抗衡离子直到按照-1/-2电荷处理往往能得到比较好的效果

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

11#
发表于 Post on 2024-9-7 15:17:37 | 只看该作者 Only view this author
ljh123 发表于 2024-9-6 14:52
请问老师,关于第二个方法,我可以理解为:先通过DFT分别计算待测铁络合物和二茂铁的绝对氧化还原电位, ...

对,我之前有文章就是这么算的
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-23 23:29 , Processed in 0.219055 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list