求助：计算锂离子电池中酯溶剂还原电位的两条反应路径的选择

bluehour

各位老师好，我在做一个关于锂离子电池中酯类溶剂还原稳定性的机器学习课题，为更贴近负极界面真实反应环境，在 DFT 计算中引入 Li⁺参与酯分子的还原过程，目标计算相对于 Li⁺/Li 标度的还原电位。
目前我考虑了两个可能的反应路径：
(1) M + Li⁺ + e⁻ → [M–Li]•
(2) [M–Li]⁺ + e⁻ → [M–Li]•
我的Gaussian计算流程是：结构优化和频率热力学校正用 M062X/Def2SVP，气相高精度单点能用 M062X/Def2TZVPP，溶剂化能用 M052X/6-31G(d)，溶剂模型为 SMD(DMSO)。
此前已拜读过 sob 老师的 谈谈隐式溶剂模型下溶解自由能和体系自由能的计算 http://sobereva.com/327，了解到对于 Li⁺ 这样的离子，单独用隐式溶剂模型计算溶解自由能可能存在较大误差，往往需要引入显式溶剂才能得到合理结果，比较复杂。而对于路径 (2)，虽然避开了单独计算 Li⁺ 的溶解自由能，但实际电解液中溶剂分子并非完全处于与 Li⁺ 的配位状态，将 [M–Li]⁺ 作为标准反应物是否在热力学上合理，我对此也有些疑惑。

在此想请教各位老师：
1. 路径 (1) 和路径 (2) 哪种更适合用来评估酯类溶剂在负极表面的还原稳定性？
2. 如果我的整体计算方案存在不合理之处，恳请各位老师批评指正，也欢迎推荐更合适的计算方法或相关参考文献。
非常感谢！

sobereva

可以先用CP2K跑跑AIMD，看看Li+主要以什么状态存在

bluehour

sobereva 发表于 2026-7-31 00:43
可以先用CP2K跑跑AIMD，看看Li+主要以什么状态存在

感谢老师指点！在计算还原电位时我考虑用公式：E_red = −△G/(nF) − 1.37 ，这里 1.37 V 是标准氢电极绝对电势（4.44 V）与 Li⁺/Li 标准电极电势（−3.07 V）的差值，作为校正因子将绝对氧化还原电位标度转换为 Li⁺/Li 标度。请问老师这样在计算时还需要考虑电子的自由能吗？具体该如何计算和操作呢