计算化学公社

标题: 请教高斯计算redox potentials的问题 [打印本页]

作者
Author: yuda_li 时间: 2024-2-19 22:56
标题: 请教高斯计算redox potentials的问题
本帖最后由 yuda_li 于 2024-2-19 23:25 编辑

仔细看了sob老师的帖子谈谈隐式溶剂模型下溶解自由能和体系自由能的计算 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)，也大概懂了老师所说的那个算例。
对应SMD溶剂模型表现乙醇环境的情况下用B3LYP-D3(BJ)/6-311G*对体系进行优化和振动分析，得到几何结构，然后再用这个几何结构算出环氧乙烷在气相下、标况下（298.15 K、1 atm）的自由能以及真空下的电子能量，二者相加，就得到了此体系的气相标况自由能
接下来就是通过sob老师的3.gif和4.gif计算，环氧乙烷在乙醇溶剂中的溶解自由能，因为我们最终想得到的就是溶解自由能+气相标况自由能+1.89/627.51，所得到的就是1 M标准态的环氧乙烷在乙醇中的自由能。
对于上述算例我已经大概清楚老师的思路，但是我目前还是不太清楚应该如何计算一个分子或者一个溶剂化结构（例如3个DOL溶剂和一个Li+）这样的体系的氧化还原电势。因为老师给的这个溶解过程里，似乎并没有氧化还原反应发生。
我有看到sob老师在其他的帖子下回复如何计算redox potential的高斯是否可以计算oxidation and reduction potentials？ - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com)
但是老师描述的原话是：“需要计算的是得/失电子过程的气相自由能变（推荐用热力学组合方法算），以及得/失电子前后物种的溶解自由能（推荐用SMD模型），加在一起就是在溶液下得/失电子过程的自由能变，分别根据E_red=-ΔG_red/nF和E_ox=ΔG_ox/nF，就分别得到了绝对还原/氧化势，减去4.44V就是相对于标准氢电极的。”
请问这个得失电子前后，我应该在gaussian里如何刻画它？是只需要在输入文件那将电荷从0改成1或者-1吗？
请求各位老师指点！

作者
Author: fatpig 时间: 2024-2-20 09:46
得失电子前后的物种在气相中优化，用优化完的结构再算溶解自由能...

可以参考下本研究组的论文：

https://pubs.rsc.org/en/content/articlehtml/2023/cp/d3cp01904g

作者
Author: yuda_li 时间: 2024-2-20 10:58

fatpig 发表于 2024-2-20 09:46
得失电子前后的物种在气相中优化，用优化完的结构再算溶解自由能...

可以参考下本研究组的论文：

好的好的！马上去拜读！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)