碳酸酯氟化调控用于锂离子电池：一项理论研究
一、 研究团队与发表信息
本研究由来自北京化工大学有机无机复合材料国家重点实验室、化学工程学院的杨晓华、刘向阳、韩金辉、刘志平（通讯作者，zhiping liu@buct.edu.cn）和张现仁（通讯作者，zhangxr@mail.buct.edu.cn）共同完成。研究成果以题为“Tuning Fluorination of Carbonates for Lithium-Ion Batteries: A Theoretical Study”的论文形式，发表于*The Journal of Physical Chemistry B*期刊。文章于2023年1月20日接收，并于2023年3月5日修订后发表。
二、 学术背景与研究目标
本研究隶属于电化学储能材料与计算化学交叉领域，具体聚焦于锂离子电池（LIBs）电解质溶剂的设计与优化。随着电动汽车和可再生能源存储的快速发展，对高能量密度、高安全性和长循环寿命的锂离子电池需求日益迫切。锂金属负极因其极高的理论比容量被视为下一代高能量密度电池的理想选择，但其实际应用受到锂枝晶生长和循环过程中不稳定的固体电解质界面（SEI）膜严重制约。
电解质，作为电池的“血液”，对电池的稳定运行至关重要。商用锂离子电池普遍采用碳酸酯基电解质，但其与高活性锂金属负极会发生剧烈副反应，导致形成不均匀、力学性能差的SEI膜，加速电池失效。研究表明，氟代碳酸乙烯酯（FEC）作为一种有效的电解质添加剂，能够显著改善SEI膜的稳定性和弹性。近年来，多种氟代碳酸酯被证实对电池性能有积极影响，但其效果高度依赖于氟原子的取代数量和位置（即氟化模式）。氟化不仅能提升电解质的抗氧化稳定性以匹配高电压正极，还能通过形成富含无机成分（如LiF）的稳定SEI膜来稳定负极界面。此外，氟代碳酸酯较弱的锂离子配位能力也对改善电池性能具有独特作用。
尽管已有实验报道了不同氟化模式的影响，但其背后的分子机制尚不清晰。理论计算是揭示实验现象背后机理的强大工具。例如，已有研究通过理论计算发现，氟代碳酸酯在锂金属表面的反应性以及SEI膜的组成强烈受氟化程度和位置的影响。
基于此，本研究旨在通过系统的量子化学计算，深入探究氟化模式（氟原子的数量和位置）对一系列线性碳酸酯（如碳酸二甲酯DMC、碳酸二乙酯DEC、碳酸甲乙酯EMC）和环状碳酸酯（如碳酸乙烯酯EC、碳酸丙烯酯PC、2,3-丁二醇环状碳酸酯BC）的氧化还原稳定性、与锂离子的配位能力以及还原分解产物（特别是LiF的生成）的影响。研究目标是为设计新型高性能氟代碳酸酯溶剂或添加剂提供清晰的理论指导原则和分子设计策略。
三、 详细研究流程与方法
本研究采用了一套多层次、系统的理论计算工作流程，主要包含以下几个关键步骤：
1. 分子体系设计与建模： 研究者选取了三种环状碳酸酯（EC, PC, BC）和三种线性碳酸酯（DMC, DEC, EMC）作为基础分子骨架。在此基础上，系统性地设计了多种氟化模式，包括单氟取代（F）、双氟取代（df）、三氟甲基取代（fm）、二氟亚甲基取代（fme）、氟代乙基取代（fe）以及四氟乙基取代（fee）等（参见原文图1）。每种设计分子均根据其基础碳酸酯和取代基进行命名，共计研究了数十种氟代碳酸酯分子，构成了一个用于系统分析的数据集。
2. 氧化还原电位计算： 研究采用两种方法评估分子的电化学稳定性窗口（ESW）：传统的基于最高占据分子轨道（HOMO）与最低未占分子轨道（LUMO）能级差的方法，以及更精确的热力学循环法。 * 热力学循环法：这是本研究的核心计算方法。该方法通过构建溶液相中分子氧化/还原反应的热力学循环（参见原文Scheme 1），计算反应吉布斯自由能变（ΔG），再利用能斯特方程转换为相对于Li+/Li的氧化电位（Eox）和还原电位（Ered）。计算过程考虑了气相中的绝热电离势（AIP）、绝热电子亲和势（AEA）、零点振动能校正以及熵项，并通过隐式溶剂化模型（PCM，选用二甲基亚砜DMSO模拟EC/DEC混合电解质的介电环境）计算了溶剂化自由能。最终，根据IUPAC推荐的标准氢电极与Li+/Li电极的电位关系，将绝对电位转换为实验常用的vs Li+/Li标度。 * 计算细节：所有分子结构在B3LYP/6-311+G(d)理论水平下进行几何优化，并使用了Grimme的D3色散校正。单点能计算在更精确的双杂化泛函B2PLYP-D3/def2-TZVP水平下进行。热力学量通过准刚性转子谐振子模型计算。溶剂化自由能在M06-2X/6-31+G(d)水平下使用PCM模型计算。所有密度泛函理论计算均使用Gaussian 16程序完成。
3. 锂离子配位作用研究： 为了模拟实际电池中溶剂分子与锂离子的相互作用，研究者计算了溶剂分子与Li+的结合能。结合能定义为配位复合物能量与孤立溶剂分子和Li+能量之和的差值，并进行了基组重叠误差（BSSSE）校正。此外，为了深入理解结合能的本质，研究采用了对称性匹配微扰理论（SAPT） 对结合能进行了能量分解分析，将其拆分为静电、诱导、交换和色散等相互作用分量。SAPT计算使用PSI4软件包在SAPT0/aug-cc-pVDZ水平下完成。
4. 分子性质描述符计算与相关性分析： 为了建立氟化模式与宏观性质（氧化还原电位）之间的构效关系，研究者计算了一系列分子描述符。这些描述符包括： * 几何描述符：氟原子总数（Nf）、α位氟原子数（αf）、β位氟原子数（βf）、非酯基碳原子数（Nc）、分子对称性（sym）、是否含环（ring）。 * 量子化学描述符：HOMO/LUMO能级、分子范德华体积、静电势极值（ESPmin, ESPmax）、极性表面积（PSA）、分子极性指数（MPI）、平均局部电离能极值（ALIEmin, ALIEmax）、偶极矩（dm）。 随后，研究者进行了皮尔逊相关性分析，并进一步使用F检验和互信息（MI） 方法进行特征筛选，以识别对氧化还原电位预测最重要的描述符，为未来基于机器学习的溶剂筛选奠定基础。
5. 前沿分子轨道与静电势分析： 通过可视化HOMO和LUMO的等值面图，直观分析氟原子取代对分子得失电子能力的影响区域。利用Multiwfn程序计算并绘制分子表面静电势（ESP）分布图，从静电相互作用的角度解释氟化对锂离子配位能力影响的根源。
6. 电子转移动力学与还原分解路径探索： 应用马库斯电子转移理论，估算了溶剂分子（尤其是与Li+配位后）氧化和还原反应的自由能垒，将热力学稳定性与动力学过程联系起来。特别针对α位氟代的线性碳酸酯，计算了其还原反应路径，探索了LiF原位生成的可行性，通过寻找过渡态（TS）和计算反应能量变化来验证其反应趋势。
四、 主要研究结果
1. 氟化模式对氧化还原电位的显著影响： 计算得到的氧化还原电位数据汇总于原文表1和表2。结果显示，氟化能同时提高碳酸酯的氧化电位和还原电位，但还原电位对氟化模式的响应更为敏感。这意味着氟化能拓宽电解质的电化学稳定窗口，使其更耐高压氧化，但同时（在热力学上）也更易在负极被还原。然而，电位的提升并非简单地与氟原子数量成正比，而是强烈依赖于氟原子的取代位置。例如，具有相同氟原子数的EMC_fm_α和EMC_fm_β，其氧化和还原电位存在显著差异。相关性分析进一步证实，氧化还原电位与氟原子总数高度相关，与α位氟原子数中度相关，与β位氟原子数弱相关。研究还发现，基于HOMO-LUMO能隙估算的电化学窗口与热力学循环法结果存在较大偏差，后者与实验值更为吻合，凸显了热力学循环法在电化学研究中的重要性。
2. 关键量子描述符的识别： 通过特征工程分析，研究者发现HOMO能级、最小静电势（ESPmin）、最大/最小平均局部电离能（ALIEmax/ALIEmin）和偶极矩（dm） 是预测氧化还原电位最有效的描述符。这些描述符与氧化还原电位之间存在显著的线性或非线性依赖关系（参见原文图4）。这为未来使用计算成本较低的量子化学描述符构建可解释的机器学习模型，以快速筛选新型溶剂提供了重要依据。
3. 氟化影响氧化还原稳定性的电子结构根源： 前沿分子轨道分析揭示了氟化影响稳定性的内在机制。在HOMO中，氟原子的直接贡献很小，但氟化增强了酯基团的共轭效应，从而降低了HOMO能量，提升了抗氧化能力。这种效应在两侧对称氟化时尤为明显。在LUMO中，氟原子的贡献则非常显著。两侧氟化主要表现出共轭效应，而单侧氟化主要表现出诱导效应，这导致了两者在还原电位变化上的差异。
4. 锂离子配位后的行为变化： 研究发现，与Li+配位后，溶剂的氧化和还原电位均有所升高，表明配位使溶剂在正极更稳定，但在负极更易分解。结合能计算表明，氟化会削弱碳酸酯与Li+的配位能力，且氟原子越多，配位能力越弱。SAPT能量分解分析指出，这种削弱主要源于静电相互作用的减弱，其次是对诱导和交换相互作用的影响，而色散作用贡献很小。
5. 静电势分布揭示配位能力差异的直观图像： 静电势分布图（原文图8，9）直观地展示了氟化如何改变分子的电荷分布。在碳酸酯分子中，羰基氧原子提供最负的ESP区域（亲核位点），而烷基部分提供正的ESP区域。氟化会显著改变ESP分布：α位氟化使负的ESP区域更集中，而β位氟化则使ESP分布更均匀，后者导致了配位能力的普遍减弱。对于环状碳酸酯，环上氟化与环外烷基链氟化对ESP分布和配位能力的影响也截然不同。
6. α位氟代线性碳酸酯易于原位生成LiF： 动力学和反应路径计算表明，在还原反应中，α位氟代的线性碳酸酯（如DEC_fme, DEC_fmefm）与Li+配位后，能够以较低的能量势垒原位生成LiF（参见原文图11）。这是一个关键发现，因为富含LiF的SEI膜被认为具有高的界面能和机械强度，能有效抑制锂枝晶生长，提升电池循环性能。
五、 研究结论与价值
本研究通过系统的理论计算，深入阐明了氟原子的数量和位置（氟化模式）对不同类型碳酸酯溶剂的氧化还原稳定性、锂离子配位能力及SEI膜形成产物的影响机制，并提出了以下具体的设计指导原则：
本研究的科学价值在于从电子结构、热力学和动力学多个层面，揭示了氟代碳酸酯分子结构与电池关键性能（高压稳定性、界面稳定性、溶剂化结构）之间的内在联系，建立了清晰的“结构-性质”关系图谱。其应用价值在于为实验化学家设计合成新型高性能氟代碳酸酯电解质溶剂或添加剂提供了明确、可操作的理论指南，有望加速下一代高能量密度锂金属电池电解质的研发进程。
六、 研究亮点