手机版

电解液模拟计算 — 锂离子电池电解液溶剂化结构与离子输运的MD模拟

发布时间：2026-07-21 来源：科研学术网

小中大

字号：

电解液模拟计算是锂离子电池研发中越来越离不开的工具。实验上测离子电导率用EIS（电化学阻抗谱），测溶剂化数用拉曼或NMR——但这些方法都只能给宏观平均值。想知道Li⁺的第一配位层到底是4个还是6个EC分子、溶剂化壳层的寿命是多长、游离离子和离子对的比例——这些微观信息只有MD模拟能直接给出。我做电解液模拟计算最早是为了解释一个实验现象：为什么1 M LiPF6在EC/DMC（体积比1:1）中的离子电导率在-20°C时急剧下降到常温的不到20%——实验只告诉你了结果，MD告诉你为什么。

一、电解液体系的分子动力学建模

电解液模拟计算的第一步是力场选择——这直接决定了模拟的物理可靠性。LiPF6/碳酸酯体系的MD模拟通常用OPLS-AA力场描述溶剂分子（EC、DMC、EMC等），用专门开发的离子力场描述Li⁺和PF6⁻。OPLS-AA的Lennard-Jones参数和部分电荷对碳酸酯分子的密度和介电常数预测较好（EC密度约1.32 g/cm³ vs 实验1.32，DMC密度约1.07 vs 实验1.07）。

Li⁺的力场参数对结果极其敏感。如果用标准OPLS-AA的Li⁺参数（σ约0.14 nm, ε约0.027 kcal/mol）——Li⁺的LJ半径太小——配位层过紧，第一峰g(r)位置偏低约0.2 Å，配位数偏高约1-2个。推荐使用专门针对碳酸酯体系优化的Li⁺参数（如Vatamanu-Borodin参数集，Li⁺的σ约0.18 nm）——这个参数是用更高精度的QM计算训练的，能更好地复现实验溶剂化结构。

体系构建：用Packmol按目标浓度（如1 M LiPF6 = 1 mol LiPF6 / 1 L EC/DMC混合溶剂）生成初始构型。1 M LiPF6在EC/DMC(1:1)中对应约40个Li⁺、40个PF6⁻、400个EC分子、400个DMC分子——总原子数约30000。盒子尺寸约5.5 nm——对于这样密集的离子体系，盒子至少需要5 nm以上才能避免周期性人为效应。

二、溶剂化结构的RDF解析

电解液模拟计算的核心输出之一是径向分布函数g(r)——它直接描述Li⁺周围的溶剂化结构。Li-O(EC)和Li-O(DMC)的RDF第一峰位置约2.0-2.1 Å（Li⁺与碳酸酯羰基氧的配位距离），第一峰积分到第一谷（约2.8 Å）得到Li⁺的第一配位层配位数。

100 ns MD对于1 M LiPF6/EC/DMC体系给出的典型结果：Li-O(EC)第一峰在2.05 Å，配位数约3.8，Li-O(DMC)第一峰在2.08 Å，配位数约1.2——总第一配位层配位数约5.0。这与拉曼光谱和NMR的实验估计（溶剂化数4-6）一致。EC的配位数显著高于DMC——原因是EC的羰基氧有更高的部分负电荷（OPLS力场中EC羰基O的电荷约-0.50 e，DMC的约-0.43 e）——EC对Li⁺的静电吸引力更强。

溶剂化壳层的寿命用Li-O（EC/DMC）的驻留时间自相关函数C(t) = <h(t)·h(0)>分析——h(t)在配体在配位层中时=1，离开时为0。第一配位层的平均驻留时间在常温（298 K）约1-2 ns——这意味着在100 ns的MD中，每个Li⁺会和周围溶剂分子交换配位约50-100次。驻留时间的温度依赖性可以提取溶剂化壳层解离的活化能——通常在15-25 kJ/mol范围内。

三、离子扩散系数与电导率

电解液模拟计算的第二个核心输出是离子扩散系数D——从均方位移（MSD）的线性拟合得到。Li⁺和PF6⁻的MSD需要分别计算：gmx msd -f traj.xtc -s topol.tpr -n index.ndx。MSD-t曲线在100-1000 ps区间做线性拟合，D = slope/6。

1 M LiPF6/EC/DMC在298 K下的典型MD结果：D_Li⁺约1.5×10⁻¹⁰ m²/s，D_PF6⁻约3.2×10⁻¹⁰ m²/s。PF6⁻扩散更快——因为PF6⁻比Li⁺大但带电量小（-1 vs +1），且Li⁺被溶剂化壳层包裹后有效尺寸更大（含第一配位层的Li⁺有效半径约3.5 Å，远大于裸离子半径0.76 Å）。实验值（PFG-NMR）：D_Li⁺约2.3×10⁻¹⁰ m²/s，D_PF6⁻约4.5×10⁻¹⁰ m²/s——MD值偏低约35%。差值的主要来源是力场（OPLS-AA + 离子参数）对溶剂粘度的高估——OPLS力场模拟的EC/DMC粘度比实验值高约25-40%，导致所有离子的扩散系数系统性偏低。

从扩散系数估算离子电导率——用Nernst-Einstein近似：σ ≈ (e²/k_B T) × (N_Li×D_Li⁺ + N_PF6×D_PF6⁻)。但Nernst-Einstein假设所有离子独立运动——忽略了离子配对和交叉相关。更准确的计算是Green-Kubo公式——对体系的净电流自相关函数做时间积分。两者的比值（反比于Haven比H_R = D*/D_σ）反映了离子关联的程度——H_R<1说明有离子配对降低了有效导电离子数。

四、离子配对与聚集行为

电解液模拟计算中一个容易被忽略但很重要的分析是离子配对。Li⁺和PF6⁻在电解液中不只是独立溶剂化的自由离子——它们会形成接触离子对（CIP，Li⁺和PF6⁻直接接触，中间无溶剂分子）和溶剂分隔离子对（SSIP，Li⁺和PF6⁻之间隔着一层溶剂分子）。

从MD轨迹中识别CIP/SSIP的方法很简单：统计Li⁺-P(PF6)距离。如果d<3.5 Å——接触离子对；3.5<d<5.5 Å——溶剂分隔离子对；d>5.5 Å——自由离子。在1 M LiPF6/EC/DMC的MD中，接触离子对的比例约5-10%（常温），SSIP约15-20%，自由离子约70-80%。温度升高到333 K——CIP比例降到3-5%，自由离子比例升到约85%——这与电导率随温度升高的实验趋势一致（部分原因是离子解离增加了电荷载流子浓度）。

离子聚集的另一个形式是形成更大的离子团簇——在MD中可以观察到2-5个Li⁺和2-5个PF6⁻形成的瞬态团簇（寿命约100-500 ps）。这些团簇是中性的——对电导率没有贡献——它们的存在减少了有效载流子数目，是MD和实验Nernst-Einstein电导率差异的另一个来源。

五、温度依赖性——低温电解液的困境

回到开头的问题：为什么低温下电导率暴跌？MD模拟给出了清晰的微观解释。在238 K（-35°C）下重复模拟1 M LiPF6/EC/DMC——此时EC的熔点约309 K，DMC约275 K，混合溶剂在238 K处于过冷液态。MD结果显示：Li⁺的扩散系数从298 K的1.5×10⁻¹⁰骤降到238 K的约0.08×10⁻¹⁰ m²/s——下降了95%。

两个主要原因：第一，溶剂的粘度随温度呈Arrhenius型增长（EC的玻璃化转变温度约150 K），在238 K时粘度比298 K高约20倍——粘度直接反比于扩散系数（Stokes-Einstein关系）。第二，低温下Li⁺-EC配位键的寿命从298 K的~~1 ns增加到238 K的~~50 ns——配位层变得”僵硬”，Li⁺运动不再是独立扩散，而是拖着一整个溶剂化壳层在粘稠环境中缓慢移动。

这些MD结果指导了电解液配方优化方向：添加低粘度共溶剂（如乙酸乙酯EA，熔点和粘度都远低于EC）可以打乱Li⁺的规整溶剂化层，降低配位层寿命，从而在低温下维持一定的Li⁺扩散速率。